Magnetic flux generated and intensified by the solar dynamo emerges into the solar atmosphere, forming active regions (ARs) including sunspots. Existing theories of flux emergence suggest that the magnetic flux can rise buoyantly through the convection zone but is trapped at the photosphere, while its further rising into the atmosphere resorts to the Parker buoyancy instability. To trigger such an instability, the Lorentz force in the photosphere needs to be as large as the gas pressure gradient to hold up an extra amount of mass against gravity. This naturally results in a strongly non-force-free photosphere, which is indeed shown in typical idealized numerical simulations of flux tube buoyancy from below the photosphere into the corona. Here we conduct a statistical study of the extents of normalized Lorentz forces and torques in the emerging photospheric magnetic field with a substantially large sample of SDO/HMI vector magnetograms. We found that the photospheric field has a rather small Lorentz force and torque on average, and thus is very close to a force-free state, which is not consistent with theories as well as idealized simulations of flux emergence. Furthermore, the small extents of forces and torques seem not to be influenced by the emerging AR's size, the emergence rate, or the non-potentiality of the field. This result puts an important constraint on future development of theories and simulations of flux emergence.

Thanks to the high quality data of space missions, the detection of mixed modes has become possible in numerous stars. In this work, we investigate how the mixing character of dipolar mixed modes changes with stellar evolution, as well as with frequency within each stellar model. This is achieved by monitoring the variations in the coupling strength and the period spacing of dipolar mixed modes in red-giant models. These parameters are measured by fitting the asymptotic expansion of mixed modes to the model frequencies of a grid of red-giant models with masses between 1.0 and 2.0 $M_\odot$ and three different chemical abundances. The coupling strength and the period spacing decrease with stellar evolution. We find that the slopes of their decreasing trends depend on the radial order of the pressure mode component. A non-negligible increase of the coupling strength with frequency by up to around 40\% is found in the observable frequency range for a set of red-giant models. On the contrary, no significant changes of the period spacing with frequency are found. The changes in the mixing character of the modes are in most cases affected by the model mass and metallicity. Buoyancy glitches also have an impact on the mixing character. Significant fluctuations in the estimated coupling strength and period spacing are found for models approaching the luminosity bump, if the glitch impact of the frequencies is not considered in the applied asymptotic expansion.

The Transiting Exoplanet Survey Satellite (TESS) is an all-sky survey mission aiming to search for exoplanets that transit bright stars. The high-quality photometric data of TESS are excellent for the asteroseismic study of solar-like stars. In this work, we present an asteroseismic analysis of the red-giant star HD~222076 hosting a long-period (2.4 yr) giant planet discovered through radial velocities. Solar-like oscillations of HD~222076 are detected around $203 \, \mu$Hz by TESS for the first time. Asteroseismic modeling, using global asteroseismic parameters as input, yields a determination of the stellar mass ($M_\star = 1.12 \pm 0.12\, M_\odot$), radius ($R_\star = 4.34 \pm 0.21\,R_\odot$), and age ($7.4 \pm 2.7\,$Gyr), with precisions greatly improved from previous studies. The period spacing of the dipolar mixed modes extracted from the observed power spectrum reveals that the star is on the red-giant branch burning hydrogen in a shell surrounding the core. We find that the planet will not escape the tidal pull of the star and be engulfed into it within about $800\,$Myr, before the tip of the red-giant branch is reached.

The orientation, chirality, and dynamics of solar eruptive filaments is a key to understanding the magnetic field of coronal mass ejections (CMEs) and therefore to predicting the geoeffectiveness of CMEs arriving at Earth. However, confusion and contention remain over the relationship between the filament chirality, magnetic helicity, and sense of rotation during eruption. To resolve the ambiguity in observations, in this paper, we used stereoscopic observations to determine the rotation direction of filament apex and the method proposed by Chen et al. (2014) to determine the filament chirality. Our sample of 12 eruptive active-region filaments establishes a strong one-to-one relationship, i.e., during the eruption, sinistral/dextral filaments (located in the southern/northern hemisphere) rotate clockwise/counterclockwise when viewed from above, and corroborates a weak hemispheric preference, i.e., a filament and related sigmoid both exhibit a forward (reverse) S shape in the southern (northern) hemisphere, which suggests that the sigmoidal filament is associated with a low-lying magnetic flux rope with its axis dipped in the middle. As a result of rotation, the projected S shape of a filament is anticipated to be reversed during eruption.

Solar eruptions are the most powerful drivers of space weather. To understand their cause and nature, it is crucial to know how the coronal magnetic field evolves before eruption. Here we study the formation process of a relatively large-scale magnetic flux rope (MFR) in active region NOAA~12371 that erupts with a major flare and coronal mass ejection on 2015 June 21. A data-driven numerical magnetohydrodynamic model is employed to simulate three-dimensional coronal magnetic field evolution of one-day duration before the eruption. Comparison between the observed features and our modeled magnetic field discloses how the pre-eruption MFR forms. Initially, the magnetic field lines were weakly twisted as being simple sheared arcades. Then a long MFR was formed along the polarity inversion line due to the complex photospheric motion, which is mainly shearing rather than twisting. The presence of the MFR is evidenced by a coherent set of magnetic field lines with twist number above unity. Below the MFR a current sheet is shown in the model, suggesting that tether-cutting reconnection plays a key role in the MFR formation. The MFR's flux grows as more and more field lines are twisted due to continuous injection of magnetic helicity by the photospheric motions. Meanwhile, the height of the MFR's axis increases monotonely from its formation. By an analysis of the decay index of its overlying field, we suggest that it is because the MFR runs into the torus instability regime and becomes unstable that finally triggers the eruption.

Context. The large quantity of high-quality asteroseismic data that obtained from space-based photometric missions and the accuracy of the resulting frequencies motivate a careful consideration of the accuracy of computed oscillation frequencies of stellar models, when applied as diagnostics of the model properties. Aims. Based on models of red-giant stars that have been independently calculated using different stellar evolution codes, we investigate the extent to which the differences in the model calculation affect the model oscillation frequencies. Methods. For each of the models, which cover four different masses and different evolution stages on the red-giant branch, we computed full sets of low-degree oscillation frequencies using a single pulsation code and, from these frequencies, typical asteroseismic diagnostics. In addition, we carried out preliminary analyses to relate differences in the oscillation properties to the corresponding model differences. Results. In general, the differences in asteroseismic properties between the different models greatly exceed the observational precision of these properties, in particular for the nonradial modes whose mixed acoustic and gravity-wave character makes them sensitive to the structure of the deep stellar interior. In some cases, identifying these differences led to improvements in the final models presented here and in Paper I; here we illustrate particular examples of this. Conclusions. Further improvements in stellar modelling are required in order fully to utilise the observational accuracy to probe intrinsic limitations in the modelling. However, our analysis of the frequency differences and their relation to stellar internal properties provides a striking illustration of the potential of the mixed modes of red-giant stars for the diagnostics of stellar interiors.

Since the onset of the `space revolution' of high-precision high-cadence photometry, asteroseismology has been demonstrated as a powerful tool for informing Galactic archaeology investigations. The launch of the NASA TESS mission has enabled seismic-based inferences to go full sky -- providing a clear advantage for large ensemble studies of the different Milky Way components. Here we demonstrate its potential for investigating the Galaxy by carrying out the first asteroseismic ensemble study of red giant stars observed by TESS. We use a sample of 25 stars for which we measure their global asteroseimic observables and estimate their fundamental stellar properties, such as radius, mass, and age. Significant improvements are seen in the uncertainties of our estimates when combining seismic observables from TESS with astrometric measurements from the Gaia mission compared to when the seismology and astrometry are applied separately. Specifically, when combined we show that stellar radii can be determined to a precision of a few percent, masses to 5-10% and ages to the 20% level. This is comparable to the precision typically obtained using end-of-mission Kepler data

We present our measurements of the H$\alpha$, [OIII], and [OII] luminosity functions as part of the Lyman Alpha Galaxies at Epoch of Reionization (LAGER) survey using our samples of 1577 $z = 0.47$ H$\alpha$-, 3933 $z = 0.93$ [OIII]-, and 5367 $z = 1.59$ [OII]-selected emission line galaxies in a single 3 deg$^2$ CTIO/Blanco DECam pointing of the COSMOS field. Our observations reach 5$\sigma$ depths of $8.2\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ and comoving volumes of $(1-7)\times10^{5}$ Mpc$^3$ making our survey one of the deepest narrowband surveys. We measure the observed luminosity functions and find best-fits of $\phi^\star = 10^{-3.16\pm0.09}$ Mpc$^{-3}$ and $L^\star = 10^{41.72\pm0.09}$ erg s$^{-1}$ for H$\alpha$, $\phi^\star = 10^{-2.16^{+0.10}_{-0.12}}$ Mpc$^{-3}$ and $L^\star = 10^{41.38^{+0.07}_{-0.06}}$ erg s$^{-1}$ for [OIII], and $\phi^\star = 10^{-1.97^{+0.07}_{-0.07}}$ Mpc$^{-3}$ and $L^\star = 10^{41.66\pm0.03}$ erg s$^{-1}$ for [OII], with $\alpha$ fixed to $-1.75$, $-1.6$, and $-1.3$, respectively. An excess of bright $> 10^{42}$ erg s$^{-1}$ [OIII] emitters is observed and may be due to AGN contamination. Dust corrections are applied assuming $A_{\rm{H}\alpha} = 1$ mag. We also design our own empirical rest-frame $g - r$ calibration using SDSS DR12 data, test it against our $z = 0.47$ H$\alpha$ emitters with $z$COSMOS $1$D spectra, and calibrate it for $(g - r)$ between $-0.8$ and $1.3$ mag. Dust and AGN-corrected star formation rate densities (SFRDs) are measured as $\log_{10}\rho_{\rm{SFR}}/(\rm{M}_\odot\ \rm{yr}^{-1}\ \rm{Mpc}^{-3}) = -1.63\pm0.04$, $-1.07\pm0.06$, and $-0.90\pm0.10$ for H$\alpha$, [OIII], and [OII], respectively. We find our [OIII] and [OII] samples fully trace cosmic star formation activity at their respective redshifts in comparison to multi-wavelength SFRDs, while the H$\alpha$ sample traces $\sim 70$ percent of the total $z = 0.47$ SFRD.

Two X-class solar flares occurred on 2017 September 6 from active region NOAA 12673: the first one is a confined X2.2 flare, and it is followed only $\sim 3$ hours later by the second one, which is the strongest flare in solar cycle 24, reaching X9.3 class and accompanied with a coronal mass ejection. Why these two X-class flares occurred in the same position with similar magnetic configurations, but one is eruptive while the other is not? Here we track the coronal magnetic field evolution via nonlinear force-free field extrapolations from a time sequence of vector magnetograms with high cadence. A detailed analysis of the magnetic field shows that a magnetic flux rope (MFR) forms and grows gradually before the first flare, and shortly afterwards, the MFR's growth is significantly enhanced with a much faster rise in height, from far below the threshold of torus instability to above it, while the magnetic twist only increases mildly. Combining EUV observations and the magnetic field extrapolation, we found that overlying the MFR is a null-point magnetic topology, where recurrent brightening is seen after the first flare. We thus suggest a scenario to interpret the occurrence of the two flares. The first flare occurred since the MFR reached a high enough height to activate the null point, and its continuous expansion forces the null-point reconnection recurrently. Such reconnection weakens the overlying field, allowing the MFR to rise faster, which eventually crosses the threshold of torus instability and triggers the second, eruptive flare.

For a better understanding of magnetic field in the solar corona and dynamic activities such as flares and coronal mass ejections, it is crucial to measure the time-evolving coronal field and accurately estimate the magnetic energy. Recently, a new modeling technique called the data-driven coronal field model, in which the time evolution of magnetic field is driven by a sequence of photospheric magnetic and velocity field maps, has been developed and revealed the dynamics of flare-productive active regions. Here we report on the first qualitative and quantitative assessment of different data-driven models using a magnetic flux emergence simulation as a ground-truth (GT) data set. We compare the GT field with those reconstructed from the GT photospheric field by four data-driven algorithms. It is found that, at least, the flux rope structure is reproduced in all coronal field models. Quantitatively, however, the results show a certain degree of model dependence. In most cases, the magnetic energies and relative magnetic helicity are comparable to or at most twice of the GT values. The reproduced flux ropes have a sigmoidal shape (consistent with GT) of various sizes, a vertically-standing magnetic torus, or a packed structure. The observed discrepancies can be attributed to the highly non-force-free input photospheric field, from which the coronal field is reconstructed, and to the modeling constraints such as the treatment of background atmosphere, the bottom boundary setting, and the spatial resolution.

We report a detailed observational study of two quasi-periodic fast-propagating (QFP) magnetosonic wave events occurred on 2011 March 09 and 10, respectively. Interestingly, both the two events have two wave trains (WTs): one main and strong (WT-1) whereas the second appears small and weak (WT-2). Peculiar and common characteristics of the two events are observed, namely: 1) the two QFP waves are accompanied with brightenings during the whole stage of the eruptions; 2) both the two main wave trains are nearly propagating along the same direction; 3) EUV waves are found to be associated with the two events. Investigating various aspects of the target events, we argue that: 1) the second event is accompanied with a flux rope eruption during the whole stage; 2) the second event eruption produces a new filament-like (FL) dark feature; 3) the ripples of the two WT-2 QFP waves seem to result from different triggering mechanisms. Based on the obtained observational results, we propose that the funnel-like coronal loop system is indeed playing an important role in the two WT-1 QFP waves. The development of the second WT-2 QFP wave can be explained as due to the dispersion of the main EUV front. The co-existence of the two events offer thereby a significant opportunity to reveal what driving mechanisms and structures are tightly related to the waves.

(Abridged). We introduce the Aarhus Red Giants Challenge, a series of detailed comparisons between widely used stellar evolution and oscillation codes aiming at establishing the minimum level of uncertainties in properties of red giants arising solely from numerical implementations. Using 9 state-of-the-art stellar evolution codes, we defined a set of input physics and physical constants for our calculations and calibrated the convective efficiency to a specific point on the main sequence. We produced evolutionary tracks and stellar structure models at fixed radius along the red-giant branch for masses of 1.0 M$_\odot$, 1.5 M$_\odot$, 2.0 M$_\odot$, and 2.5 M$_\odot$, and compared the predicted stellar properties. Once models have been calibrated on the main sequence we find a residual spread in the predicted effective temperatures across all codes of ~20 K at solar radius and ~30-40 K in the RGB regardless of the considered stellar mass. The predicted ages show variations of 2-5% (increasing with stellar mass) which we track down to differences in the numerical implementation of energy generation. The luminosity of the RGB-bump shows a spread of about 10% for the considered codes, which translates into magnitude differences of ~0.1 mag in the optical V-band. We also compare the predicted [C/N] abundance ratio and found a spread of 0.1 dex or more for all considered masses. Our comparisons show that differences at the level of a few percent still remain in evolutionary calculations of red giants branch stars despite the use of the same input physics. These are mostly due to differences in the energy generation routines and interpolation across opacities, and call for further investigations on these matters in the context of using properties of red giants as benchmarks for astrophysical studies.

We present the largest-ever sample of 79 Ly$\alpha$ emitters (LAEs) at $z\sim$ 7.0 selected in the COSMOS and CDFS fields of the LAGER project (the Lyman Alpha Galaxies in the Epoch of Reionization). Our newly amassed ultradeep narrowband exposure and deeper/wider broadband images have more than doubled the number of LAEs in COSMOS, and we have selected 30 LAEs in the second field CDFS. We detect two large-scale LAE-overdense regions in the COSMOS that are likely protoclusters at the highest redshift to date. We perform injection and recovery simulations to derive the sample incompleteness. We show significant incompleteness comes from blending with foreground sources, which however has not been corrected in LAE luminosity functions in {the} literature. The bright end bump in the Ly$\alpha$ luminosity function in COSMOS is confirmed with 6 (2 newly selected) luminous LAEs (L$_{Ly\alpha}$ $>$ 10$^{43.3}$ erg s$^{-1}$). Interestingly, the bump is absent in CDFS, in which only one luminous LAE is detected. Meanwhile, the faint end luminosity functions from the two fields well agree with each other. The 6 luminous LAEs in COSMOS coincide with 2 LAE-overdense regions, while such regions are not seen in CDFS. The bright-end luminosity function bump could be attributed to ionized bubbles in a patchy reionization. It appears associated with cosmic overdensities, thus supports an inside-out reionization topology at $z$ $\sim$ 7.0, i.e., the high density peaks were ionized earlier compared to the voids. An average neutral hydrogen fraction of $x_{HI}$ $\sim$ 0.2 -- 0.4 is derived at $z\sim$ 7.0 based on the cosmic evolution of the Ly$\alpha$ luminosity function.

The characterisation of stellar cores may be accomplished through the modelling of asteroseismic data from stars exhibiting either gravity-mode or mixed-mode pulsations, potentially shedding light on the physical processes responsible for the production, mixing, and segregation of chemical elements. In this work we validate against model data an analytical expression for the period spacing that will facilitate the inference of the properties of stellar cores, including the detection and characterisation of buoyancy glitches (strong chemical gradients). This asymptotically-based analytical expression is tested both in models with and without buoyancy glitches. It does not assume that glitches are small and, consequently, predicts non-sinusoidal glitch-induced period-spacing variations, as often seen in model and real data. We show that the glitch position and width inferred from the fitting of the analytical expression to model data consisting of pure gravity modes are in close agreement (typically better than 7$\%$ relative difference) with the properties measured directly from the stellar models. In the case of fitting mixed-mode model data, the same expression is shown to reproduce well the numerical results, when the glitch properties are known a priori. In addition, the fits performed to mixed-mode model data reveal a frequency dependence of the coupling coefficient, $q$, for a moderate-luminosity red-giant-branch model star. Finally, we find that fitting the analytical expression to the mixed-mode period spacings may provide a way to infer the frequencies of the pure acoustic dipole modes that would exist if no coupling took place between acoustic and gravity waves.

Magnetic flux ropes (MFRs) are thought to be the central structure of solar eruptions, and their ideal MHD instabilities can trigger the eruption. Here we performed a study of all the MFR configurations that lead to major solar flares, either eruptive or confined, from 2011 to 2017 near the solar disk center. The coronal magnetic field is reconstructed from observed magnetograms, and based on magnetic twist distribution, we identified the MFR, which is defined as a coherent group of magnetic field lines winding an axis with more than one turn. It is found that 90% of the events possess pre-flare MFRs, and their three-dimensional structures are much more complex in details than theoretical MFR models. We further constructed a diagram based on two parameters, the magnetic twist number which controls the kink instability (KI), and the decay index which controls the torus instability (TI). It clearly shows lower limits for TI and KI thresholds, which are $n_{\rm crit} = 1.3$ and $|T_w|_{\rm crit} = 2$, respectively, as all the events above $n_{\rm crit}$ and nearly 90% of the events above $|T_w|_{\rm crit}$ erupted. Furthermore, by such criterion, over 70% of the events can be discriminated between eruptive and confined flares, and KI seems to play a nearly equally important role as TI in discriminating between the two types of flare. There are more than half of events with both parameters below the lower limits, and 29% are eruptive. These events might be triggered by magnetic reconnection rather than MHD instabilities.

Coronal mass ejections (CMEs) play a decisive role in driving space weather, especially, the fast ones (e.g., with speeds above $800$~km~s$^{-1}$). Understanding the trigger mechanisms of fast CMEs can help us gaining important information in forecasting them. The filament eruptions accompanied with CMEs provide a good tracer in studying the early evolution of CMEs. Here we surveyed 66 filament-accompanied fast CMEs to analyse the correlation between the trigger mechanisms, namely either magnetic reconnection or ideal MHD process, associated flares, and CME speeds. Based on the data gathering from SDO, GONG and STEREO, we find that: (1) Active region (AR) filament and intermediate filaments (IFs) eruptions show a higher probability for producing fast CMEs than quiet Sun (QS) filaments, while the probability of polar crown (PC) filament eruptions is zero in our statistic; (2) AR filament eruptions that produce fast CMEs are more likely triggered by magnetic reconnection, while QS and IFs are more likely triggered by ideal MHD process; (3) For AR filaments and IFs, it seems that the specific trigger mechanism does not have a significant influence on the resulted CME speeds, while for the QS filaments, the ideal MHD mechanism can more likely generate a faster CME; (4) Comparing with previous statistic study, the onset heights of filament eruptions and the decay indexes of the overlying field show some differences: for AR filaments and IFs, the decay indexes are larger and much closer to the theoretical threshold, while for QS filaments, the onset heights are higher than those obtained in previous results.

We present the discovery of TOI-197.01, the first transiting planet identified by the Transiting Exoplanet Survey Satellite (TESS) for which asteroseismology of the host star is possible. TOI-197 (HIP116158) is a bright (V=8.2 mag), spectroscopically classified subgiant which oscillates with an average frequency of about 430 muHz and displays a clear signature of mixed modes. The oscillation amplitude confirms that the redder TESS bandpass compared to Kepler has a small effect on the oscillations, supporting the expected yield of thousands of solar-like oscillators with TESS 2-minute cadence observations. Asteroseismic modeling yields a robust determination of the host star radius (2.943+/-0.064 Rsun), mass (1.212 +/- 0.074 Msun) and age (4.9+/-1.1 Gyr), and demonstrates that it has just started ascending the red-giant branch. Combining asteroseismology with transit modeling and radial-velocity observations, we show that the planet is a "hot Saturn" (9.17+/-0.33 Rearth) with an orbital period of ~14.3 days, irradiance of 343+/-24 Fearth, moderate mass (60.5 +/- 5.7 Mearth) and density (0.431+/-0.062 gcc). The properties of TOI-197.01 show that the host-star metallicity - planet mass correlation found in sub-Saturns (4-8 Rearth) does not extend to larger radii, indicating that planets in the transition between sub-Saturns and Jupiters follow a relatively narrow range of densities. With a density measured to ~15%, TOI-197.01 is one of the best characterized Saturn-sized planets to date, augmenting the small number of known transiting planets around evolved stars and demonstrating the power of TESS to characterize exoplanets and their host stars using asteroseismology.

We present the design for the first narrowband filter NB964 for the Dark Energy Camera (DECam), which is operated on the 4m Blanco Telescope at the Cerro Tololo Inter-American Observatory. The NB964 filter profile is essentially defined by maximizing the power of searching for Lyman alpha emitting galaxies (LAEs) in the epoch of reionization, with the consideration of the night sky background in the near-infrared and the DECam quantum efficiency. The NB964 filter was manufactured by Materion in 2015. It has a central wavelength of 964.2 nm and a full width at half maximum (FWHM) of 9.2 nm. An NB964 survey named LAGER (Lyman Alpha Galaxies in the Epoch of Reionization) has been ongoing since December 2015. Here we report results of lab tests, on-site tests and observations with the NB964 filter. The excellent performances of this filter ensure that the LAGER project is able to detect LAEs at z~7 with a high efficiency.

We spectroscopically confirmed two narrow-band selected redshift 7.0 Ly$\alpha$ galaxies and studied their restframe UV spectra. The Ly$\alpha$ and other UV nebular lines are very useful to confirm the galactic redshifts and diagnose the different mechanisms driving the ionizing emission. We observed two narrowband-selected $z$=7.0 Ly$\alpha$ candidates in the LAGER Chandra Deep Field South (CDFS) field with IMACS at Magellan telescope and confirmed they are Ly$\alpha$ emitters at $z$=6.924 and 6.931. In one galaxy, we also obtained deep NIR spectroscopy, which yields non-detections of the high-ionization UV nebular lines. We measured upper-limits of the ratios of CIV$\lambda$1548/Ly$\alpha$, HeII$\lambda$1640/Ly$\alpha$, OIII]$\lambda$1660/Ly$\alpha$, and CIII]$\lambda$1909/Ly$\alpha$ from the NIR spectra. These upper-limits imply that the ionizing emission in this galaxy is dominated by normal star formation instead of AGN.

We observed an H$\alpha$ surge occurred in the active region NOAA 12401 on 2015 August 17, and discuss its trigger mechanism, kinematic and thermal properties. It is suggested that this surge is caused by a chromospheric reconnection which ejects cool and dense material with the transverse velocity about 21-28 km s$^{-1}$ and the initial doppler velocity of 12 km s$^{-1}$. This surge is similar to the injection of newly formed filament materials from their footpoints, except that the surge here occurred in a relatively weak magnetic environment of ~100 G. Thus we discuss the possibility of filament material replenishment via the erupting mass in such a weak magnetic field, which is often associated with quiescent filaments. It is found that the local plasma can be heated up to about 1.3 times of original temperature, which results in an acceleration about -0.017 km s$^{-2}$. It can lift the dense material up to 10 Mm and higher with a inclination angle smaller than 50$^{\circ}$, namely typical height of active region filaments. But it can hardly inject the material up to those filaments higher than 25 Mm, namely some quiescent filaments. Thus we think injection model does not work well in the formation of quiescent filaments.

Solar flares are often associated with coronal eruptions, but there are confined ones without eruption, even for some X-class flares. How such large flares occurred and why they are confined are still not well understood. Here we studied a confined X2.2 flare in NOAA 12673 on 2017 September 6. It exhibits two episodes of flare brightening with rather complex, atypical ribbons. Based on topology analysis of extrapolated coronal magnetic field, we revealed that there is a two-step magnetic reconnection process during the flare. Prior to the flare, there is a magnetic flux rope (MFR) with one leg rooted in a rotating sunspot. Neighboring to the leg is a magnetic null-point structure. The sunspot drives the MFR to expand, pushing magnetic flux to the null point, and reconnection is first triggered there. The disturbance from the null-point reconnection triggers the second reconnection, i.e., a tether-cutting reconnection below the rope. However, these two reconnections failed to produce an eruption, because the rope is firmly held by its strapping flux. Furthermore, we compared this flare with an eruptive X9.3 flare in the same region with 2 hours later, which has a similar MFR configuration. The key difference between them is that, for the confined flare, the MFR is fully below the threshold of torus instability, while for the eruptive one, the MFR reaches entirely above the threshold. This study provides a good evidence supporting that reconnection alone may not be able to trigger eruption, rather, MHD instability plays a more important role.

Three-dimensional magnetic topology is crucial to understanding the explosive release of magnetic energy in the corona during solar flares. Much attention has been given to the pre-flare magnetic topology to identify candidate sites of magnetic reconnection, yet it is unclear how the magnetic reconnection and its attendant topological changes shape the eruptive structure and how the topology evolves during the eruption. Here we employed a realistic, data-constrained magnetohydrodynamic simulation to study the evolving magnetic topology for an X9.3 eruptive flare that occurred on 2017 September 6. The simulation successfully reproduces the eruptive features and processes in unprecedented detail. The numerical results reveal that the pre-flare corona contains multiple twisted flux systems with different connections, and during the eruption, these twisted fluxes form a coherent flux rope through tether-cutting-like magnetic reconnection below the rope. Topological analysis shows that the rising flux rope is wrapped by a quasi-separatrix layer, which intersects itself below the rope, forming a topological structure known as hyperbolic flux tube, where a current sheet develops, triggering the reconnection. By mapping footpoints of the newly-reconnected field lines, we are able to reproduce both the spatial location and, for the first time, the temporal separation of the observed flare ribbons, as well as the dynamic boundary of the flux rope's feet. Futhermore, the temporal profile of the total reconnection flux is comparable to the soft X-ray light curve. Such a sophisticated characterization of the evolving magnetic topology provides important insight into the eventual understanding and forecast of solar eruptions.

We develop a novel approach in exploring the joint dependence of halo bias on multiple halo properties using Gaussian process regression. Using a $\Lambda$CDM $N$-body simulation, we carry out a comprehensive study of the joint bias dependence on halo structure, formation history and environment. We show that the bias is a multivariate function of halo properties that falls into three regimes. For massive haloes, halo mass explains the majority of bias variation. For early-forming haloes, bias depends sensitively on the recent mass accretion history. For low-mass and late-forming haloes, bias depends more on the structure of a halo such as its shape and spin. Our framework enables us to convincingly prove that $V_\mathrm{max}/V_\mathrm{vir}$ is a lossy proxy of formation time for bias modelling, whereas the mass, spin, shape and formation time variables are non-redundant with respect to each other. Combining mass and formation time largely accounts for the mass accretion history dependence of bias. Combining all the internal halo properties fully accounts for the density profile dependence inside haloes, and predicts the clustering variation of individual haloes to a $20\%$ level at $\sim 10\mathrm{Mpc}h^{-1}$. When an environmental density is measured outside $1\mathrm{Mpc}h^{-1}$ from the halo centre, it outperforms and largely accounts for the bias dependence on the internal halo structure, explaining the bias variation above a level of $30\%$.

Ly$\alpha$ Emitters (LAEs) may represent an important galaxy population in the low mass regime. We present our deep narrowband imaging surveys in the COSMOS and ECDF-S fields and study the properties of LAEs at $z=2.23\pm0.03$. The narrowband surveys conducted at Magellan II telescope allow us to obtain a sample of 452 LAEs reaching a $5\sigma$ limiting magnitude of $\sim26$ mag. Our Ly$\alpha$ luminosity functions extend to $10^{41.8}$ erg s$^{-1}$ with steep faint-end slope. Using multi-wavelength ancillary data, especially the deep Spitzer/IRAC 3.6$\mu$m and 4.5$\mu$m photometric data, we obtained reliable stellar mass estimates for 130 IRAC-detected LAEs, spanning a range of $8 < {\rm log}(M_\star/M_\odot)< 11.5$. For the remaining IRAC-undetected LAEs, the median-stacked spectral energy distribution yields a stellar mass of ${\rm log}(M_\star/M_\odot)=7.97^{+0.05}_{-0.07}$ and the rest-frame ultraviolet emission indicates a median star formation rate of ${\rm log} (SFR/M_\odot$ yr$^{-1})=-0.14\pm0.35$. There are six LAEs detected by the Spitzer/MIPS 24$\mu$m or even Herschel far-infrared observations. Taking into account the six MIR/FIR detected LAEs, our LAEs cover a wide range in the star formation rate (${\rm 1<SFR<2000}$ M$_\odot$ yr$^{-1}$). Although LAEs as a population are diverse in their stellar properties, they are mostly low-mass star-forming galaxies and follow the star formation main sequence relations or their extrapolations to the low-mass end, implying a normal star-forming nature of LAEs. The clustering analysis indicates that our LAEs reside in dark matter halos with ${\rm <\log(M_{h}/M_{\odot})> =10.8^{+0.56}_{-1.1}}$, suggesting that they are progenitors of local Large Magellanic Cloud-like galaxies.

The Stellar Observation Network Group (SONG) is an initiative to build a worldwide network of 1m telescopes with highprecision radial-velocity spectrographs. Here we analyse the first radial-velocity time series of a red-giant star measured by the SONG telescope at Tenerife. The asteroseismic results demonstrate a major increase in the achievable precision of the parameters for redgiant stars obtainable from ground-based observations. Reliable tests of the validity of these results are needed, however, before the accuracy of the parameters can be trusted. We analyse the first SONG time series for the star 46 LMi, which has a precise parallax and an angular diameter measured from interferometry, and therefore a good determination of the stellar radius. We use asteroseismic scaling relations to obtain an accurate mass, and modelling to determine the age. A 55-day time series of high-resolution, high S/N spectra were obtained with the first SONG telescope. We derive the asteroseismic parameters by analysing the power spectrum. To give a best guess on the large separation of modes in the power spectrum, we have applied a new method which uses the scaling of Kepler red-giant stars to 46 LMi. Several methods have been applied: classical estimates, seismic methods using the observed time series, and model calculations to derive the fundamental parameters of 46 LMi. Parameters determined using the different methods are consistent within the uncertainties. We find the following values for the mass M (scaling), radius R (classical), age (modelling), and surface gravity (combining mass and radius): M = 1.09 +- 0.04Msun, R = 7.95 +- 0.11 Rsun, age t = 8.2 +- 1.9Gy, and log g = 2.674 +- 0.013.

Three-dimensional magnetic topology of solar flare plays a crucial role in understanding its explosive release of magnetic energy in the corona. However, such three-dimensional coronal magnetic field is still elusive in direct observation. Here we realistically simulate the magnetic evolution during the eruptive process of a great flare, using a numerical magnetohydrodynamic model constrained by observed solar vector magnetogram. The numerical results reveal that the pre-flare corona contains multi-set twisted magnetic flux, which forms a coherent rope during the eruption. The rising flux rope is wrapped by a quasi-separatrix layer, which intersects itself below the rope, forming a hyperbolic flux tube and magnetic reconnection is triggered there. By tracing the footprint of the newly-reconnected field lines, we reproduce both the spatial location and its temporal evolution of flare ribbons with an expected accuracy in comparison of observed images. This scenario strongly confirms the three-dimensional version of standard flare model.

We investigate the star formation history and the dust attenuation in the galaxy merger Mrk848. Thanks to the multiwavelength photometry from the ultraviolet (UV) to the infrared (IR), and MaNGA's integral field spectroscopy, we are able to study this merger in a detailed way. We divide the whole merger into the core and tail regions, and fit both the optical spectrum and the multi-band spectral energy distribution (SED) to models to obtain the star formation properties for each region respectively. We find that the color excess of stars in the galaxy $E(B-V)_s^{SED}$ measured with the multi-band SED fitting is consistent with that estimated both from the infrared excess (the ratio of IR to UV flux) and from the slope of the UV continuum. Furthermore, the reliability of the $E(B-V)_s^{SED}$ is examined with a set of mock SEDs, showing that the dust attenuation of the stars can be well constrained by the UV-to-IR broadband SED fitting. The dust attenuation obtained from optical continuum $E(B-V)_s^{spec}$ is only about half of $E(B-V)_s^{SED}$. The ratio of the $E(B-V)_s^{spec}$ to the $E(B-V)_g$ obtained from the Balmer decrement is consistent with the local value (around 0.5). The difference between the results from the UV-to-IR data and the optical data is consistent with the picture that younger stellar populations are attenuated by an extra dust component from the birth clouds compared to older stellar populations which are only attenuated by the diffuse dust. Both with the UV-to-IR SED fitting and the spectral fitting, we find that there is a starburst younger than 100~Myr in one of the two core regions, consistent with the scenario that the interaction-induced gas inflow can enhance the star formation in the center of galaxies.

Mixed modes have been extensively observed in post-main-sequence stars by the Kepler and CoRoT space missions. The mixture of the p and g modes can be measured by the dimensionless coefficient $q$, the so-called coupling strength factor. In this paper we discuss the utility of the phase shifts $\theta$ from the eigenvalue condition for mixed modes as a tool to characterize dipolar mixed modes from the theoretical as well as the practical point of view. Unlike the coupling strength, whose variation in a given star is very small over the relevant frequency range, the phase shifts vary significantly for different modes. The analysis in terms of $\theta$ can also provide a better understanding of the pressure and gravity radial order for a given mixed mode. Observed frequencies of the Kepler red-giant star KIC 3744043 are used to test the method. The results are very promising.

Solar flare ribbons provide an important clue to the magnetic reconnection process and associated magnetic field topology in the solar corona. We detected a large-scale secondary flare ribbon of a circular shape that developed in association with two successive M-class flares and one CME. The ribbon revealed interesting properties such as 1) a quasi-circular shape and enclosing the central active region; 2) the size as large as 500\arcsec\, by 650\arcsec\,, 3) successive brightenings in the clockwise direction at a speed of \kms{160} starting from the nearest position to the flaring sunspots, 4) radial contraction and expansion in the northern and the southern part, respectively at speeds of $\leq$ \kms{10}. Using multi-wavelength data from \textit{SDO}, \textit{RHESSI}, XRT, and Nobeyama, along with magnetic field extrapolations, we found that: 1) the secondary ribbon location is consistent with the field line footpoints of a fan-shaped magnetic structure that connects the flaring region and the ambient decaying field; 2) the second M2.6 flare occurred when the expanding coronal loops driven by the first M2.0 flare encountered the background decayed field. 3) Immediately after the second flare, the secondary ribbon developed along with dimming regions. Based on our findings, we suggest that interaction between the expanding sigmoid field and the overlying fan-shaped field triggered the secondary reconnection that resulted in the field opening and formation of the quasi-circular secondary ribbon. We thus conclude that interaction between the active region and the ambient large-scale fields should be taken into account to fully understand the entire eruption process.

The high planetary multiplicity revealed by Kepler implies that Transit Time Variations (TTVs) are intrinsically common. The usual procedure for detecting these TTVs is biased to long-period, deep transit planets whereas most transiting planets have short periods and shallow transits. Here we introduce the Spectral Approach to TTVs technique that allows expanding the TTVs catalog towards lower TTV amplitude, shorter orbital period, and shallower transit depth. In the Spectral Approach we assume that a sinusoidal TTV exists in the data and then calculate the improvement to $\chi^2$ this model allows over that of linear ephemeris model. This enables detection of TTVs even in cases where the transits are too shallow so individual transits cannot be timed. The Spectral Approach is more sensitive due to the reduced number of free parameters in its model. Using the Spectral Approach, we: (a) detect 131 new periodic TTVs in Kepler data (an increase of ~2/3 over a previous TTV catalog); (b) Constrain the TTV periods of 34 long-period TTVs and reduce amplitude errors of known TTVs; (c) Identify cases of multi-periodic TTVs, for which absolute planetary mass determination may be possible. We further extend our analysis by using perturbation theory assuming small TTV amplitude at the detection stage, which greatly speeds up our detection (to a level of few seconds per star). Our extended TTVs sample shows no deficit of short period or low amplitude transits, in contrast to previous surveys in which the detection schemes were significantly biased against such systems.

As a fundamental magnetic structure in the solar corona, electric current sheets (CSs) can form either prior to or during solar flare, and they are essential for magnetic energy dissipation in the solar corona by enabling magnetic reconnection. However static reconstruction of CS is rare, possibly due to limitation inherent in available coronal field extrapolation codes. Here we present the reconstruction of a large-scale pre-flare CS in solar active region 11967 using an MHD-relaxation model constrained by SDO/HMI vector magnetogram. The CS is found to be associated with a set of peculiar homologous flares that exhibit unique X-shaped ribbons and loops occurring in a quadrupolar magnetic configuration. This is evidenced by that the field lines traced from the CS to the photosphere form an X shape which nearly precisely reproduces the shape of the observed flare ribbons, suggesting that the flare is a product of the dissipation of the CS through reconnection. The CS forms in a hyperbolic flux tube, which is an intersection of two quasi-separatrix layers. The recurrence of the X-shaped flares might be attributed to the repetitive formation and dissipation of the CS, as driven by the photospheric footpoint motions. These results demonstrate the power of data-constrained MHD model in reproducing CS in the corona as well as providing insight into the magnetic mechanism of solar flares.

We carried out targeted ALMA observations of 129 fields in the COSMOS region at 1.25 mm, detecting 152 galaxies at S/N$\geq$5 with an average continuum RMS of 150 $\mu$Jy. These fields represent a S/N-limited sample of AzTEC / ASTE sources with 1.1 mm S/N$\geq$4 over an area of 0.72 square degrees. Given ALMA's fine resolution and the exceptional spectroscopic and multiwavelength photometric data available in COSMOS, this survey allows us unprecedented power in identifying submillimeter galaxy counterparts and determining their redshifts through spectroscopic or photometric means. In addition to 30 sources with prior spectroscopic redshifts, we identified redshifts for 113 galaxies through photometric methods and an additional nine sources with lower limits, which allowed a statistically robust determination of the redshift distribution. We have resolved 33 AzTEC sources into multi-component systems and our redshifts suggest that nine are likely to be physically associated. Our overall redshift distribution peaks at $z\sim$2.0 with a high redshift tail skewing the median redshift to $\tilde{z}$=2.48$\pm$0.05. We find that brighter millimeter sources are preferentially found at higher redshifts. Our faintest sources, with S$_{1.25 \rm mm}$<1.25 mJy, have a median redshift of $\tilde{z}$=2.18$\pm$0.09, while the brightest sources, S$_{1.25 \rm mm}$>1.8 mJy, have a median redshift of $\tilde{z}$=3.08$\pm$0.17. After accounting for spectral energy distribution shape and selection effects these results are consistent with several previous submillimeter galaxy surveys, and moreover, support the conclusion that the submillimeter galaxy redshift distribution is sensitive to survey depth.

Narrowband imaging is a highly successful approach for finding large numbers of high redshift Lya emitting galaxies (LAEs) up to z~6.6. However, at z>~7 there are as yet only 3 narrowband selected LAEs with spectroscopic confirmations (two at z~6.9-7.0, one at z~7.3), which hinders extensive studies on cosmic reionization and galaxy evolution at this key epoch. We have selected 23 candidate z~6.9 LAEs in COSMOS field with the large area narrowband survey LAGER (Lyman-Alpha Galaxies at the End of Reionization). In this work we present spectroscopic followup observations of 12 candidates using IMACS on Magellan. For 9 of these, the observations are sufficiently deep to detect the expected lines. Lya emission lines are identified in six sources (yielding a success rate of 2/3), including 3 luminous LAEs with Lya luminosities of L(Lya) ~ 10^{43.5} erg/s, the highest among known spectroscopically confirmed galaxies at >~7.0. This triples the sample size of spectroscopically confirmed narrowband selected LAEs at z>~7, and confirms the bright end bump in the Lya luminosity function we previously derived based on the photometric sample, supporting a patchy reionization scenario. Two luminous LAEs appear physically linked with projected distance of 1.1 pMpc and velocity difference of ~ 170 km/s. They likely sit in a common ionized bubble produced by themselves or with close neighbors, which reduces the IGM attenuation of Lya. A tentative narrow NV${\lambda}$1240 line is seen in one source, hinting at activity of a central massive black hole with metal rich line emitting gas.

The effect of metallicity on the granulation activity in stars is still poorly understood. Available spectroscopic parameters from the updated APOGEE-\textit{Kepler} catalog, coupled with high-precision photometric observations from NASA's \textit{Kepler} mission spanning more than four years of observation, make oscillating red giant stars in open clusters crucial testbeds. We determine the role of metallicity on the stellar granulation activity by discriminating its effect from that of different stellar properties such as surface gravity, mass, and temperature. We analyze 60 known red giant stars belonging to the open clusters NGC 6791, NGC 6819, and NGC 6811, spanning a metallicity range from [Fe/H] $\simeq -0.09$ to $0.32$. The parameters describing the granulation activity of these stars and their $\nu_\mathrm{max}$, are studied by considering the different masses, metallicities, and stellar evolutionary stages. We derive new scaling relations for the granulation activity, re-calibrate existing ones, and identify the best scaling relations from the available set of observations. We adopted the Bayesian code DIAMONDS for the analysis of the background signal in the Fourier spectra of the stars. We performed a Bayesian parameter estimation and model comparison to test the different model hypotheses proposed in this work and in the literature. Metallicity causes a statistically significant change in the amplitude of the granulation activity, with a dependency stronger than that induced by both stellar mass and surface gravity. We also find that the metallicity has a significant impact on the corresponding time scales of the phenomenon. The effect of metallicity on the time scale is stronger than that of mass. A higher metallicity increases the amplitude of granulation and meso-granulation signals and slows down their characteristic time scales toward longer periods.

We determine the physical properties of a sample of SMGs in the COSMOS field that were pre-selected at the observed wavelength of $\lambda_{\rm obs}=1.1$ mm, and followed up at $\lambda_{\rm obs}=1.3$ mm with ALMA. We used MAGPHYS to fit the panchromatic (ultraviolet to radio) SEDs of 124 of the target SMGs, 19.4% of which are spectroscopically confirmed. The SED analysis was complemented by estimating the gas masses of the SMGs by using the $\lambda_{\rm obs}=1.3$ mm emission as a tracer of the molecular gas. The sample median and 16th-84th percentile ranges of the stellar masses, SFRs, dust temperatures, and dust and gas masses were derived to be $\log(M_{\star}/{\rm M}_{\odot})=11.09^{+0.41}_{-0.53}$, ${\rm SFR}=402^{+661}_{-233}$ ${\rm M}_{\odot}~{\rm yr}^{-1}$, $T_{\rm dust}=39.7^{+9.7}_{-7.4}$ K, $\log(M_{\rm dust}/{\rm M}_{\odot})=9.01^{+0.20}_{-0.31}$, and $\log(M_{\rm gas}/{\rm M}_{\odot})=11.34^{+0.20}_{-0.23}$, respectively. The median gas-to-dust ratio and gas fraction were found to be $120^{+73}_{-30}$ and $0.62^{+0.27}_{-0.23}$, respectively. We found that 57.3% of our SMGs populate the main sequence (MS) of star-forming galaxies, while 41.9% of the sources lie above the MS by a factor of >3 (one source lies below the MS). The largest 3 GHz radio sizes are found among the MS sources. Those SMGs that appear irregular in the rest-frame UV are predominantly starbursts, while the MS SMGs are mostly disk-like. The larger radio-emitting sizes of the MS SMGs compared to starbursts is a likely indication of their more widespread, less intense star formation. The irregular UV morphologies of the starburst SMGs are likely to echo their merger nature. Our results suggest that the transition from high-$z$ SMGs to local ellipticals via compact, quiescent galaxies (cQGs) at $z \sim 2$ might not be universal, and the latter population might also descend from the so-called blue nuggets.

Solar flares and coronal mass ejections (CMEs) are the most powerful explosions in the Sun. They are major sources of potentially destructive space weather conditions. However, the possible causes of their initiation remain controversial. By using high resolution data observed by NST of BBSO, supplemented by Solar Dynamics Observatory (SDO) observations, we present unusual observations of a small-scale emerging flux rope near a large sunspot, whose eruption produced an M-class flare and a coronal mass ejection. The presence of the small-scale flux rope was indicated by static nonlinear force free field (NLFFF) extrapolation as well as data-driven MHD modeling of the dynamic evolution of the coronal 3D magnetic field. During the emergence of the flux rope, rotation of satellite sunspots at the footpoints of the flux rope was observed. Meanwhile, the Lorentz force, magnetic energy, vertical current, and transverse fields were increasing during this phase. The free energy from the magnetic flux emergence and twisting magnetic fields is sufficient to power the M-class flare. These observations for the first time present the complete process from the emergence of the small-scale flux rope to the production of solar eruptions.

We present the first results from the ongoing LAGER project (Lyman Alpha Galaxies in the Epoch of Reionization), which is the largest narrowband survey for $z \sim$ 7 galaxies to date. Using a specially built narrowband filter NB964 for the superb large-area Dark-Energy Camera (DECam) on the NOAO/CTIO 4m Blanco telescope, LAGER has collected 34 hours NB964 narrowband imaging data in the 3 deg$^2$ COSMOS field. We have identified 23 Lyman Alpha Emitter (LAE) candidates at $z$ = 6.9 in the central 2-deg$^2$ region, where DECam and public COSMOS multi-band images exist. The resulting luminosity function can be described as a Schechter function modified by a significant excess at the bright end (4 galaxies with $L_{Ly\alpha} \sim $ 10$^{43.4\pm0.2}$ erg s$^{-1}$). The number density at $L_{Ly\alpha}\sim$ 10$^{43.4\pm0.2}$ erg s$^{-1}$ is little changed from z= 6.6, while at fainter $L_{Ly\alpha}$ it is substantially reduced. Overall, we see a fourfold reduction in Ly$\alpha$ luminosity density from $z$ = 5.7 to 6.9. Combined with a more modest evolution of the continuum UV luminosity density, this suggests a factor of $\sim 3$ suppression of Ly$\alpha$ by radiative transfer through the $z \sim$ 7 intergalactic medium (IGM). It indicates an IGM neutral fraction $x_{HI}$ $\sim$ 0.4--0.6 (assuming Ly$\alpha$ velocity offsets of 100-200 km s$^{-1}$). The changing shape of the Ly$\alpha$ luminosity function between $z\lesssim 6.6$ and $z=6.9$ supports the hypothesis of ionized bubbles in a patchy reionization at $z\sim$ 7.

Magnetic field extrapolation is an important tool to study the three-dimensional (3D) solar coronal magnetic field which is difficult to directly measure. Various analytic models and numerical codes exist but their results often drastically differ. Thus a critical comparison of the modeled magnetic field lines with the observed coronal loops is strongly required to establish the credibility of the model. Here we compare two different non-potential extrapolation codes, a non-linear force-free field code (CESE-MHD-NLFFF) and a non-force-free field (NFFF) code in modeling a solar active region (AR) that has a sigmoidal configuration just before a major flare erupted from the region. A 2D coronal-loop tracing and fitting method is employed to study the 3D misalignment angles between the extrapolated magnetic field lines and the EUV loops as imaged by SDO/AIA. It is found that the CESE-MHD-NLFFF code with preprocessed magnetogram performs the best, outputting a field which matches the coronal loops in the AR core imaged in AIA 94 {\AA} with a misalignment angle of ~10 degree. This suggests that the CESE-MHD-NLFFF code, even without using the information of coronal loops in constraining the magnetic field, performs as good as some coronal-loop forward-fitting models. For the loops as imaged by AIA 171 {\AA} in the outskirts of the AR, all the codes including the potential-field give comparable results of mean misalignment angle (~30 degree). Thus further improvement of the codes is needed for a better reconstruction of the long loops enveloping the core region.

With SDO observations and a data-constrained MHD model, we identify a confined multi-ribbon flare occurred on 2010 October 25 in solar active region 11117 as a magnetic bald patch (BP) flare with strong evidences. From the photospheric magnetic field observed by SDO/HMI, we find there is indeed magnetic BPs on the PILs which match parts of the flare ribbons. From the 3D coronal magnetic field derived from a MHD relaxation model constrained by the vector magnetograms, we find strikingly good agreement of the BP separatrix surface (BPSS) footpoints with the flare ribbons, and the BPSS itself with the hot flaring loop system. Moreover, the triggering of the BP flare can be attributed to a small flux emergence under the lobe of the BPSS, and the relevant change of the coronal magnetic field through the flare is well reproduced by the pre-flare and post-flare MHD solutions, which match the corresponding pre and post-flare AIA observations, respectively. Our work contributes to the study of non-typical flares that constitute the majority of solar flares but cannot be explained by the standard flare model.

We determine the radio size distribution of a large sample of 152 SMGs in COSMOS that were detected with ALMA at 1.3 mm. For this purpose, we used the observations taken by the VLA-COSMOS 3 GHz Large Project. One hundred and fifteen of the 152 target SMGs were found to have a 3 GHz counterpart. The median value of the major axis FWHM at 3 GHz is derived to be $4.6\pm0.4$ kpc. The radio sizes show no evolutionary trend with redshift, or difference between different galaxy morphologies. We also derived the spectral indices between 1.4 and 3 GHz, and 3 GHz brightness temperatures for the sources, and the median values were found to be $\alpha=-0.67$ and $T_{\rm B}=12.6\pm2$ K. Three of the target SMGs, which are also detected with the VLBA, show clearly higher brightness temperatures than the typical values. Although the observed radio emission appears to be predominantly powered by star formation and supernova activity, our results provide a strong indication of the presence of an AGN in the VLBA and X-ray-detected SMG AzTEC/C61. The median radio-emitting size we have derived is 1.5-3 times larger than the typical FIR dust-emitting sizes of SMGs, but similar to that of the SMGs' molecular gas component traced through mid-$J$ line emission of CO. The physical conditions of SMGs probably render the diffusion of cosmic-ray electrons inefficient, and hence an unlikely process to lead to the observed extended radio sizes. Instead, our results point towards a scenario where SMGs are driven by galaxy interactions and mergers. Besides triggering vigorous starbursts, galaxy collisions can also pull out the magnetised fluids from the interacting disks, and give rise to a taffy-like synchrotron-emitting bridge. This provides an explanation for the spatially extended radio emission of SMGs, and can also cause a deviation from the well-known IR-radio correlation.

The detection and analysis of oscillations in binary star systems is critical in understanding stellar structure and evolution. This is partly because such systems have the same initial chemical composition and age. Solar-like oscillations have been detected by Kepler in both components of the asteroseismic binary HD 176465. We present an independent modelling of each star in this binary system. Stellar models generated using MESA (Modules for Experiments in Stellar Astrophysics) were fitted to both the observed individual frequencies and complementary spectroscopic parameters. The individual theoretical oscillation frequencies for the corresponding stellar models were obtained using GYRE as the pulsation code. A Bayesian approach was applied to find the probability distribution functions of the stellar parameters using AIMS (Asteroseismic Inference on a Massive Scale) as the optimisation code. The ages of HD 176465 A and HD 176465 B were found to be 2.81 $\pm$ 0.48 Gyr and 2.52 $\pm$ 0.80 Gyr, respectively. These results are in agreement when compared to previous studies carried out using other asteroseismic modelling techniques and gyrochronology.

We study the physical mechanism of a major X-class solar flare that occurred in the super NOAA active region (AR) 12192 using a data-driven numerical magnetohydrodynamic (MHD) modeling complemented with observations. With the evolving magnetic fields observed at the solar surface as bottom boundary input, we drive an MHD system to evolve self-consistently in correspondence with the realistic coronal evolution. During a two-day time interval, the modeled coronal field has been slowly stressed by the photospheric field evolution,which gradually created a large-scale coronal current sheet, i.e., a narrow layer with intense current, in the core of the AR. The current layer was successively enhanced until it became so thin that a tether-cutting reconnection between the sheared magnetic arcades was set in, which led to a flare. The modeled reconnecting field lines and their footpoints match well the observed hot flaring loops and the flare ribbons, respectively, suggesting that the model has successfully "reproduced" the macroscopic magnetic process of the flare. In particular, with simulation, we explained why this event is a confined eruption-the consequent of the reconnection is the shared arcade instead of a newly formed flux rope. We also found much weaker magnetic implosion effect comparing to many other X-class flares

We present a narrowband survey with three adjacent filters for z=2.8--2.9 Lyman Alpha Emitter (LAE) galaxies in the Extended Chandra Deep Field South (ECDFS), along with spectroscopic followup. With a complete sample of 96 LAEs in the narrowband NB466, we confirm a large-scale structure at z~ 2.8. Compared to the blank field in NB470 and NB475, the LAE density excess in the NB466 field is ~6.0+/-0.8 times the standard deviation expected at z~2.8, assuming a linear bias of 2. The overdense large scale structure in NB466 can be decomposed into 4 protoclusters, whose overdensities are 4.6 - 6.6. These 4 protoclusters are expected to evolve into a Coma-like cluster at z~ 0. In the meanwhile, we investigate the average star-formation rates derived from Ly{\alpha}, rest-frame UV and X-ray, the Ly{\alpha} luminosity functions, the Ly{\alpha} photon densities and their dependence on the environment. We find that the Ly{\alpha} photon density in the overdense field (NB466) is ~50\% higher than that in the blank field (NB470+NB475). The 3 brightest LAEs, including a quasar at z=2.81, are all detected in X-ray and in NB466. These three LAE-AGNs contribute an extra 20--30\% Ly{\alpha} photon density. Furthermore, we find that LAEs in overdense regions are younger and less dusty. We conclude that the structure we found is a significant and rare density peak, and narrowband imaging is an efficient method to detect and study such structures in the high-z universe.

In Zombie Vortex Instability (ZVI), perturbations excite critical layers in stratified, rotating shear flow (as in protoplanetary disks), causing them to generate vortex layers, which roll-up into anticyclonic zombie vortices and cyclonic vortex sheets. The process is self-sustaining as zombie vortices perturb new critical layers, spawning a next generation of zombie vortices. Here, we focus on two issues: the minimum threshold of perturbations that trigger self-sustaining vortex generation, and the properties of the late-time zombie turbulence on large and small scales. The critical parameter that determines whether ZVI is triggered is the magnitude of the vorticity on the small scales (and not velocity), the minimum Rossby number needed for instability is $Ro_{crit}\sim0.2$ for $\beta\equiv N/\Omega = 2$, where $N$ is the Brunt-V\"ais\"al\"a frequency. While the threshold is set by vorticity, it is useful to infer a criterion on the Mach number, for Kolmogorov noise, the critical Mach number scales with Reynolds number: $Ma_{crit}\sim Ro_{crit}Re^{-1/2}$. In protoplanetary disks, this is $Ma_{crit}\sim10^{-6}$. On large scales, zombie turbulence is characterized by anticyclones and cyclonic sheets with typical Rossby number $\sim$0.3. The spacing of the cyclonic sheets and anticyclones appears to have a "memory" of the spacing of the critical layers. On the small scales, zombie turbulence has no memory of the initial conditions and has a Kolmogorov-like energy spectrum. While our earlier work was in the limit of uniform stratification, we have demonstrated that ZVI works for non-uniform Brunt-V\"ais\"al\"a frequency profiles that may be found in protoplanetary disks.

Solar active region (AR) 11283 is a very magnetically complex region and it has produced many eruptions. However, there exists a non-eruptive filament in the plage region just next to an eruptive one in the AR, which gives us an opportunity to perform a comparison analysis of these two filaments. The coronal magnetic field extrapolated using a CESE-MHD-NLFFF code (Jiang & Feng 2013) reveals that two magnetic flux ropes (MFRs) exist in the same extrapolation box supporting these two filaments, respectively. Analysis of the magnetic field shows that the eruptive MFR contains a bald-patch separatrix surface (BPSS) co-spatial very well with a pre-eruptive EUV sigmoid, which is consistent with the BPSS model for coronal sigmoids. The magnetic dips of the non-eruptive MFRs match H{\alpha} observation of the non-eruptive filament strikingly well, which strongly supports the MFR-dip model for filaments. Compared with the non-eruptive MFR/filament (with a length of about 200 Mm), the eruptive MFR/filament is much smaller (with a length of about 20 Mm), but it contains most of the magnetic free energy in the extrapolation box and holds a much higher free energy density than the non-eruptive one. Both the MFRs are weakly twisted and cannot trigger kink instability. The AR eruptive MFR is unstable because its axis reaches above a critical height for torus instability, at which the overlying closed arcades can no longer confine the MFR stably. On the contrary, the quiescent MFR is very firmly held by its overlying field, as its axis apex is far below the torus-instability threshold height. Overall, this comparison investigation supports that MFR can exist prior to eruption and the ideal MHD instability can trigger MFR eruption.

In the solar corona, magnetic flux rope is believed to be a fundamental structure accounts for magnetic free energy storage and solar eruptions. Up to the present, the extrapolation of magnetic field from boundary data is the primary way to obtain fully three-dimensional magnetic information of the corona. As a result, the ability of reliable recovering coronal magnetic flux rope is important for coronal field extrapolation. In this paper, our coronal field extrapolation code (CESE-MHD-NLFFF, Jiang & Feng 2012) is examined with an analytical magnetic flux rope model proposed by Titov & Demoulin (1999), which consists of a bipolar magnetic configuration holding an semi-circular line-tied flux rope in force-free equilibrium. By using only the vector field in the bottom boundary as input, we test our code with the model in a representative range of parameter space and find that the model field is reconstructed with high accuracy. Especially, the magnetic topological interfaces formed between the flux rope and the surrounding arcade, i.e., the "hyperbolic flux tube" and "bald patch separatrix surface", are also reliably reproduced. By this test, we demonstrate that our CESE-MHD-NLFFF code can be applied to recovering magnetic flux rope in the solar corona as long as the vector magnetogram satisfies the force-free constraints.

There is considerable interest in hydrodynamic instabilities in dead zones of protoplanetary disks as a mechanism for driving angular momentum transport and as a source of particle-trapping vortices to mix chondrules and incubate planetesimal formation. We present simulations with a pseudo-spectral anelastic code and with the compressible code Athena, showing that stably stratified flows in a shearing, rotating box are violently unstable and produce space-filling, sustained turbulence dominated by large vortices with Rossby numbers of order 0.2-0.3. This Zombie Vortex Instability (ZVI) is observed in both codes and is triggered by Kolmogorov turbulence with Mach numbers less than 0.01. It is a common view that if a given constant density flow is stable, then stable vertical stratification should make the flow even more stable. Yet, we show that sufficient vertical stratification can be unstable to ZVI. ZVI is robust and requires no special tuning of boundary conditions, or initial radial entropy or vortensity gradients (though we have studied ZVI only in the limit of infinite cooling time). The resolution of this paradox is that stable stratification allows for a new avenue to instability: baroclinic critical layers. ZVI has not been seen in previous studies of flows in rotating, shearing boxes because those calculations frequently lacked vertical density stratification and/or sufficient numerical resolution. Although we do not expect appreciable angular momentum transport from ZVI in the small domains in this study, we hypothesize that ZVI in larger domains with compressible equations may lead to angular transport via spiral density waves.

We report the discovery of Kepler-432b, a giant planet ($M_b = 5.41^{+0.32}_{-0.18} M_{\rm Jup}, R_b = 1.145^{+0.036}_{-0.039} R_{\rm Jup}$) transiting an evolved star $(M_\star = 1.32^{+0.10}_{-0.07} M_\odot, R_\star = 4.06^{+0.12}_{-0.08} R_\odot)$ with an orbital period of $P_b = 52.501129^{+0.000067}_{-0.000053}$ days. Radial velocities (RVs) reveal that Kepler-432b orbits its parent star with an eccentricity of $e = 0.5134^{+0.0098}_{-0.0089}$, which we also measure independently with asterodensity profiling (AP; $e=0.507^{+0.039}_{-0.114}$), thereby confirming the validity of AP on this particular evolved star. The well-determined planetary properties and unusually large mass also make this planet an important benchmark for theoretical models of super-Jupiter formation. Long-term RV monitoring detected the presence of a non-transiting outer planet (Kepler-432c; $M_c \sin{i_c} = 2.43^{+0.22}_{-0.24} M_{\rm Jup}, P_c = 406.2^{+3.9}_{-2.5}$ days), and adaptive optics imaging revealed a nearby (0\farcs87), faint companion (Kepler-432B) that is a physically bound M dwarf. The host star exhibits high signal-to-noise asteroseismic oscillations, which enable precise measurements of the stellar mass, radius and age. Analysis of the rotational splitting of the oscillation modes additionally reveals the stellar spin axis to be nearly edge-on, which suggests that the stellar spin is likely well-aligned with the orbit of the transiting planet. Despite its long period, the obliquity of the 52.5-day orbit may have been shaped by star-planet interaction in a manner similar to hot Jupiter systems, and we present observational and theoretical evidence to support this scenario. Finally, as a short-period outlier among giant planets orbiting giant stars, study of Kepler-432b may help explain the distribution of massive planets orbiting giant stars interior to 1 AU.

We carry out a comprehensive investigation comparing the three-dimensional magnetic field restructuring, flare energy release, and the helioseismic response, of two homologous flares, the 2011 September 6 X2.1 (FL1) and September 7 X1.8 (FL2) flares in NOAA AR 11283. In our analysis, (1) a twisted flux rope (FR) collapses onto the surface at a speed of 1.5 km/s after a partial eruption in FL1. The FR then gradually grows to reach a higher altitude and collapses again at 3 km/s after a fuller eruption in FL2. Also, FL2 shows a larger decrease of the flux-weighted centroid separation of opposite magnetic polarities and a greater change of the horizontal field on the surface. These imply a more violent coronal implosion with corresponding more intense surface signatures in FL2. (2) The FR is inclined northward, and together with the ambient fields, it undergoes a southward turning after both events. This agrees with the asymmetric decay of the penumbra observed in the peripheral regions. (3) The amounts of free magnetic energy and nonthermal electron energy released during FL1 are comparable to those of FL2 within the uncertainties of the measurements. (4) No sunquake was detected in FL1; in contrast, FL2 produced two seismic emission sources S1 and S2 both lying in the penumbral regions. Interestingly, S1 and S2 are connected by magnetic loops, and the stronger source S2 has weaker vertical magnetic field. We discuss these results in relation to the implosion process in the low corona and the sunquake generation.

High-precision space observations, such as made by the \kepler\ and \corot\ missions, allow us to detect mixed modes for $l = 1$ modes in their high signal-to-noise photometry data. By means of asteroseismology, the inner structure of red giant (RG) stars is revealed the first time with the help of mixed modes. We analyse these mixed modes of a 1.3 M$_{\sun}$ RG model theoretically from the approximate asymptotic descriptions of oscillations. While fitting observed frequencies with the eigenvalue condition for mixed modes, a good estimate of period spacing and coupling strength is also acquired for more evolved models. We show that the behaviour of the mode inertia in a given mode varies dramatically when the coupling is strong. An approximation of period spacings is also obtained from the asymptotic dispersion relation, which provides a good estimate of the coupling strength as well as period spacing when g-mode-like mixed modes are sufficiently dense. By comparing the theoretical coupling strength from the integral expression with the ones from fitting methods, we confirmed that the theoretical asymptotic equation is problematic in the evanescent region due to the potential singularities as well as the use of the Cowling approximation.

We study how to measure the galaxy merger rate from the observed close pair count. Using a high-resolution N-body/SPH cosmological simulation, we find an accurate scaling relation between galaxy pair counts and merger rates down to a stellar mass ratio of about 1:30. The relation explicitly accounts for the dependence on redshift (or time), on pair separation, and on mass of the two galaxies in a pair. With this relation, one can easily obtain the mean merger timescale for a close pair of galaxies. The use of virial masses, instead of stellar masses, is motivated by the fact that the dynamical friction time scale is mainly determined by the dark matter surrounding central and satellite galaxies. This fact can also minimize the error induced by uncertainties in modeling star formation in the simulation. Since the virial mass can be read from the well-established relation between the virial masses and the stellar masses in observation, our scaling relation can be easily applied to observations to obtain the merger rate and merger time scale. For major merger pairs (1:1-1:4) of galaxies above a stellar mass of 4*10^10 M_sun/h at z=0.1, it takes about 0.31 Gyr to merge for pairs within a projected distance of 20 kpc/h with stellar mass ratio of 1:1, while the time taken goes up to 1.6 Gyr for mergers with stellar mass ratio of 1:4. Our results indicate that a single timescale usually used in literature is not accurate to describe mergers with the stellar mass ratio spanning even a narrow range from 1:1 to 1:4.

Solar filament are commonly thought to be supported in magnetic dips, in particular, of magnetic flux ropes (FRs). In this Letter, from the observed photospheric vector magnetogram, we implement a nonlinear force-free field (NLFFF) extrapolation of a coronal magnetic FR that supports a large-scale intermediate filament between an active region and a weak polarity region. This result is the first in that current NLFFF extrapolations with presence of FRs are limited to relatively small-scale filaments that are close to sunspots and along main polarity inversion line (PIL) with strong transverse field and magnetic shear, and the existence of a FR is usually predictable. In contrast, the present filament lies along the weak-field region (photospheric field strength $\lesssim 100$ G), where the PIL is very fragmented due to small parasitic polarities on both side of the PIL and the transverse field has a low value of signal-to-noise ratio. Thus it represents a far more difficult challenge to extrapolate a large-scale FR in such case. We demonstrate that our CESE--MHD--NLFFF code is competent for the challenge. The numerically reproduced magnetic dips of the extrapolated FR match observations of the filament and its barbs very well, which supports strongly the FR-dip model for filaments. The filament is stably sustained because the FR is weakly twisted and strongly confined by the overlying closed arcades.

We present a magnetic analysis of the formation and eruption of an active region (AR) sigmoid in AR 11283 from 2011 September 4 to 6. To follow the quasi-static evolution of the coronal magnetic field, we reconstruct a time sequence of static fields using a recently developed nonlinear force-free field model constrained by the SDO/HMI vector magnetograms. A detailed analysis of the fields compared with the SDO/AIA observations suggests the following scenario for the evolution of the region. Initially, a new bipole emerges into the negative polarity of a pre-existing bipolar AR, forming a null point topology between the two flux systems. A weakly twisted flux rope (FR) is then built up slowly in the embedded core region, largely through flux-cancellation photospheric reconnections, forming a bald patch separatrix surface (BPSS) separating the FR from its ambient field. The FR grows gradually until its axis runs into a torus instability (TI) domain near the end of the third day, and the BPSS also develops a fully S-shape. Unlike in the case of standard TI, the FR does not erupt instantly since it is still attached at the photosphere along the bald patch (BP) portion of the polarity inversion line. The combined effects of the TI-driven expansion of the FR and the line-tying at the BP tear the FR into two parts with the upper portion freely expelled and the lower portion remaining behind the post-flare arcades. This process dynamically perturbs the BPSS and results in the transient enhanced brightening of the sigmoid. The accelerated expansion of the upper portion of the FR strongly pushes its envelope flux near the null point and triggers breakout reconnection at the null, as evidenced by a remarkable circular flare ribbon, which further facilitates the eruption. We discuss the important implications of these results for the formation and disruption of sigmoid region with FR.

A previously unknown instability creates space-filling lattices of 3D vortices in linearly-stable, rotating, stratified shear flows. The instability starts from an easily-excited critical layer. The layer intensifies by drawing energy from the background shear and rolls-up into vortices that excite new critical layers and vortices. The vortices self-similarly replicate to create lattices of turbulent vortices. The vortices persist for all time. This self-replication occurs in stratified Couette flows and in the dead zones of protoplanetary disks where it can de-stabilize Keplerian flows.

Numerical reconstruction/extrapolation of coronal nonlinear force-free magnetic field (NLFFF) usually takes the photospheric vector magnetogram as input at the bottom boundary. Magnetic field observed at the photosphere, however, contains force which is in conflict with the fundamental assumption of the force-free model and measurement noise which is unfavorable for practical computation. Preprocessing of the raw magnetogram has been proposed by Wiegelmann, Inhester, and Sakurai (2006) to remove the force and noise for providing better input for NLFFF modeling. In this paper we develop a new code of magnetogram preprocessing which is consistent with our extrapolation method CESE-MHD-NLFFF (Jiang, Feng, and Xiang, 2012; Jiang and Feng, 2012). Basing on a magnetic-splitting rule that a magnetic field can be split into a potential field part and a non-potential part, we split the magnetogram and deal with the two parts separately. Preprocessing of the magnetogram's potential part is based on a numerical potential field model, and the non-potential part is preprocessed using the similar optimization method of Wiegelmann et al (2006). The code is applied to the SDO/HMI data and results show that the method can remove efficiently the force and noise and improve the quality of extrapolation.

Current magnetohydrodynamic (MHD) simulations of the initiation of solar eruptions are still commonly carried out with idealized magnetic field models, whereas the realistic coronal field prior to eruptions can possibly be reconstructed from the observable photospheric field. Using a nonlinear force-free field extrapolation prior to a sigmoid eruption in AR 11283 as the initial condition in a MHD model, we successfully simulate the realistic initiation process of the eruption event, as is confirmed by a remarkable resemblance to the SDO/AIA observations. Analysis of the pre-eruption field reveals that the envelope flux of the sigmoidal core contains a coronal null and furthermore the flux rope is prone to a torus instability. Observations suggest that reconnection at the null cuts overlying tethers and likely triggers the torus instability of the flux rope, which results in the eruption. This kind of simulation demonstrates the capability of modeling the realistic solar eruptions to provide the initiation process.

Due to the absence of direct measurement, the magnetic field in the solar corona is usually extrapolated from the photosphere in numerical way. At the moment, the nonlinear force-free field (NLFFF) model dominates the physical models for field extrapolation in the low corona. Recently we have developed a new NLFFF model with MHD relaxation to reconstruct the coronal magnetic field. This method is based on CESE--MHD model with the conservation-element/solution-element (CESE) spacetime scheme. In this paper, we report the application of the CESE--MHD--NLFFF code to \SDO/HMI data with magnetograms sampled for two active regions (ARs), NOAA AR 11158 and 11283, both of which were very non-potential, producing X-class flares and eruptions. The raw magnetograms are preprocessed to remove the force and then inputted into the extrapolation code. Qualitative comparison of the results with the \SDO/AIA images shows that our code can reconstruct magnetic field lines resembling the EUV-observed coronal loops. Most important structures of the active regions are reproduced excellently, like the highly-sheared field lines that suspend filaments in AR 11158 and twisted flux rope which corresponds to a sigmoid in AR 11283. Quantitative assess of the results shows that the force-free constraint is fulfilled very well in the strong-field regions but apparently not that well in the weak-field regions because of data noise and numerical errors in the small currents.

A proposal submitted to the FNAL PAC is described to search for light sub-GeV WIMP dark matter at MiniBooNE. The possibility to steer the beam past the target and into an absorber leads to a significant reduction in neutrino background, allowing for a sensitive search for elastic scattering of WIMPs off nucleons or electrons in the detector. Dark matter models involving a vector mediator can be probed in a parameter region consistent with the required thermal relic density, and which overlaps the region in which these models can resolve the muon g-2 discrepancy. Estimates of signal significance are presented for various operational modes and parameter points. The experimental approach outlined for applying MiniBooNE to a light WIMP search may also be applicable to other neutrino facilities.

We present preliminary results on modelling KIC 7693833, the so far most metal-poor red-giant star observed by {\it Kepler}. From time series spanning several months, global oscillation parameters and individual frequencies were obtained and compared to theoretical calculations. Evolution models are calculated taking into account spectroscopic and asteroseismic constraints. The oscillation frequencies of the models were computed and compared to the {\it Kepler} data. In the range of mass computed, there is no preferred model, giving an uncertainty of about 30 K in $T_eff$, 0.02 dex in $\log g$, $0.7 \,R_\odot$ in radius and of about 2.5 Gyr in age.

We investigate the radial number density profile and the abundance distribution of faint satellites around central galaxies in the low redshift universe using the CFHT Legacy Survey. We consider three samples of central galaxies with magnitudes of M_r=-21, -22, and -23 selected from the Sloan Digital Sky Survey (SDSS) group catalog of Yang et al.. The satellite distribution around these central galaxies is obtained by cross-correlating these galaxies with the photometric catalogue of the CFHT Legacy Survey. The projected radial number density of the satellites obeys a power law form with the best-fit logarithmic slope of -1.05, independent of both the central galaxy luminosity and the satellite luminosity. The projected cross correlation function between central and satellite galaxies exhibits a non-monotonic trend with satellite luminosity. It is most pronounced for central galaxies with M_r=-21, where the decreasing trend of clustering amplitude with satellite luminosity is reversed when satellites are fainter than central galaxies by more than 2 magnitudes. A comparison with the satellite luminosity functions in the Milky Way and M31 shows that the Milky Way/M31 system has about twice as many satellites as around a typical central galaxy of similar luminosity. The implications for theoretical models are briefly discussed.

We apply a data-driven MHD model to investigate the three-dimensional (3D) magnetic field of NOAA active region (AR) 11117 around the time of a C-class confined flare occurred on 2010 October 25. The MHD model, based on the spacetime conservation-element and solution-element (CESE) scheme, is designed to focus on the magnetic-field evolution and to consider a simplified solar atomsphere with finite plasma $\beta$. Magnetic vector-field data derived from the observations at the photoshpere is inputted directly to constrain the model. Assuming that the dynamic evolution of the coronal magnetic field can be approximated by successive equilibria, we solve a time sequence of MHD equilibria basing on a set of vector magnetograms for AR 11117 taken by the Helioseismic and Magnetic Imager (HMI) on board the {\it Solar Dynamic Observatory (SDO)} around the time of flare. The model qualitatively reproduces the basic structures of the 3D magnetic field, as supported by the visual similarity between the field lines and the coronal loops observed by the Atmospheric Imaging Assembly (AIA), which shows that the coronal field can indeed be well characterized by the MHD equilibrium in most time. The magnetic configuration changes very limited during the studied time interval of two hours. A topological analysis reveals that the small flare is correlated with a bald patch (BP, where the magnetic field is tangent to the photoshpere), suggesting that the energy release of the flare can be understood by magnetic reconnection associated with the BP separatrices. The total magnetic flux and energy keep increasing slightly in spite of the flare, while the computed magnetic free energy drops during the flare with an amount of $\sim 10^{30}$ erg, which seems to be adequate to provide the energy budget of the minor C-class confined flare.

Reliable measurements of the solar magnetic field are still restricted to the photosphere, and our present knowledge of the three-dimensional coronal magnetic field is largely based on extrapolation from photospheric magnetogram using physical models, e.g., the nonlinear force-free field (NLFFF) model as usually adopted. Most of the currently available NLFFF codes have been developed with computational volume like Cartesian box or spherical wedge while a global full-sphere extrapolation is still under developing. A high-performance global extrapolation code is in particular urgently needed considering that Solar Dynamics Observatory (SDO) can provide full-disk magnetogram with resolution up to $4096\times 4096$. In this work, we present a new parallelized code for global NLFFF extrapolation with the photosphere magnetogram as input. The method is based on magnetohydrodynamics relaxation approach, the CESE-MHD numerical scheme and a Yin-Yang spherical grid that is used to overcome the polar problems of the standard spherical grid. The code is validated by two full-sphere force-free solutions from Low & Lou's semi-analytic force-free field model. The code shows high accuracy and fast convergence, and can be ready for future practical application if combined with an adaptive mesh refinement technique.

The extinct p-process nuclide 146Sm serves as an astrophysical and geochemical chronometer through measurements of isotopic anomalies of its alpha-decay daughter 142Nd. Based on analyses of 146Sm/147Sm alpha-activity and atom ratios, we determined the half-life of 146Sm to be 68 \pm 7 (1sigma) million years (Ma), which is shorter than the currently used value of 103 \pm 5 Ma. This half-life value implies a higher initial 146Sm abundance in the early solar system, (146Sm/144Sm_0 = 0.0094\pm0.0005 (2sigma), than previously estimated. Terrestrial, Lunar and Martian planetary silicate mantle differentiation events dated with 146Sm-142Nd converge to a shorter time span and in general to earlier times, due to the combined effect of the new 146Sm half-life and (146Sm/144Sm)_0 values.

Magnetic field in the solar corona is usually extrapolated from photospheric vector magnetogram using a nonlinear force-free field (NLFFF) model. NLFFF extrapolation needs a considerable effort to be devoted for its numerical realization. In this paper we present a new implementation of the magnetohydrodynamics (MHD)-relaxation method for NLFFF extrapolation. The magneto-frictional approach which is introduced for speeding the relaxation of the MHD system is novelly realized by the spacetime conservation-element and solution-element (CESE) scheme. A magnetic field splitting method is used to further improve the computational accuracy. The bottom boundary condition is prescribed by changing the transverse field incrementally to match the magnetogram, and all other artificial boundaries of the computational box are simply fixed. We examine the code by two types of NLFFF benchmark tests, the Low & Lou (1990) semi-analytic force-free solutions and a more realistic solar-like case constructed by van Ballegooijen et al. (2007). The results show that our implementation are successful and versatile for extrapolations of either the relatively simple cases or the rather complex cases which need significant rebuilding of the magnetic topology, e.g., a flux rope. We also compute a suite of metrics to quantitatively analyze the results and demonstrate that the performance of our code in extrapolation accuracy basically reaches the same level of the present best-performing code, e.g., that developed by Wiegelmann (2004).

Fusion data for $^{13}$C+$^{13}$C, $^{12}$C+$^{13}$C and $^{12}$C+$^{12}$C are analyzed by coupled-channels calculations that are based on the M3Y+repulsion, double-folding potential. The fusion is determined by ingoing-wave-boundary conditions (IWBC) that are imposed at the minimum of the pocket in the entrance channel potential. Quadrupole and octupole transitions to low-lying states in projectile and target are included in the calculations, as well as mutual excitations of these states. The effect of one-neutron transfer is also considered but the effect is small in the measured energy regime. It is shown that mutual excitations to high-lying states play a very important role in developing a comprehensive and consistent description of the measurements. Thus the shapes of the calculated cross sections for $^{12}$C+$^{13}$C and $^{13}$C+$^{13}$C are in good agreement with the data. The fusion cross sections for $^{12}$C+$^{12}$C determined by the IWBC are generally larger than the measured cross sections but they are consistent with the maxima of some of the observed peak cross sections. They are therefore expected to provide an upper limit for the extrapolation into the low-energy regime of interest to astrophysics.

We have analysed oscillations of the red giant star HD 186355 observed by the NASA Kepler satellite. The data consist of the first five quarters of science operations of Kepler, which cover about 13 months. The high-precision time-series data allow us to accurately extract the oscillation frequencies from the power spectrum. We find the frequency of the maximum oscillation power, {\nu}_max, and the mean large frequency separation, {\Delta}{\nu}, are around 106 and 9.4 {\mu}Hz respectively. A regular pattern of radial and non-radial oscillation modes is identified by stacking the power spectra in an echelle diagram. We use the scaling relations of {\Delta}{\nu} and {\nu}_max to estimate the preliminary asteroseismic mass, which is confirmed with the modelling result (M = 1.45 \pm 0.05 M_sun) using the Yale Rotating stellar Evolution Code (YREC7). In addition, we constrain the effective temperature, luminosity and radius from comparisons between observational constraints and models. A number of mixed l = 1 modes are also detected and taken into account in our model comparisons. We find a mean observational period spacing for these mixed modes of about 58 s, suggesting that this red giant branch star is in the shell hydrogen-burning phase.

Carbon burning is a crucial process for a number of important astrophysical scenarios. The lowest measured energy is around E$_{\rm c.m.}$=2.1 MeV, only partially overlapping with the energy range of astrophysical interest. The currently adopted reaction rates are based on an extrapolation which is highly uncertain because of potential resonances existing in the unmeasured energy range and the complication of the effective nuclear potential. By comparing the cross sections of the three carbon isotope fusion reactions, ${}^{12}$C+${}^{12}$C, ${}^{12}$C+${}^{13}$C and ${}^{13}$C+${}^{13}$C, we have established an upper limit on the molecular resonance strengths in ${}^{12}$C+${}^{12}$C fusion reaction. The preliminary results are presented and the impact on nuclear astrophysics is discussed.

We have measured solar-like oscillations in red giants using time-series photometry from the first 34 days of science operations of the Kepler Mission. The light curves, obtained with 30-minute sampling, reveal clear oscillations in a large sample of G and K giants, extending in luminosity from the red clump down to the bottom of the giant branch. We confirm a strong correlation between the large separation of the oscillations (Delta nu) and the frequency of maximum power (nu_max). We focus on a sample of 50 low-luminosity stars (nu_max > 100 muHz, L <~ 30 L_sun) having high signal-to-noise ratios and showing the unambiguous signature of solar-like oscillations. These are H-shell-burning stars, whose oscillations should be valuable for testing models of stellar evolution and for constraining the star-formation rate in the local disk. We use a new technique to compare stars on a single echelle diagram by scaling their frequencies and find well-defined ridges corresponding to radial and non-radial oscillations, including clear evidence for modes with angular degree l=3. Measuring the small separation between l=0 and l=2 allows us to plot the so-called C-D diagram of delta nu_02 versus Delta nu. The small separation delta nu_01 of l=1 from the midpoint of adjacent l=0 modes is negative, contrary to the Sun and solar-type stars. The ridge for l=1 is notably broadened, which we attribute to mixed modes, confirming theoretical predictions for low-luminosity giants. Overall, the results demonstrate the tremendous potential of Kepler data for asteroseismology of red giants.

Following our previous work(Jiang et al.(2008)), in which we studied the merger time-scale of galaxies in a high-resolution cosmological hydro/N-body simulation, we investigate the potential influence of uncertainties in the numerical implementation of baryonic physics on the merger time-scale. The simulation used in the previous work suffers from the overcooling problem which causes the central galaxies of large halos too massive. This may result in a shortened merger time-scale compared to that in the real universe. We run a similar simulation, but the stellar mass is significantly reduced to model another extreme case of low stellar mass. Our result shows that the merger time-scale is little affected by the star formation recipes, except for the satellites in nearly radial orbits which show a 22 percent higher time-scale in the lower stellar mass case. Since the radial orbits only account for a small part of the satellites' orbits, the fitting formula in Jiang et al.(2008) is still applicable to a reasonable accuracy, nearly immune to the uncertainty in the baryonic physics.

We study galaxy mergers using a high-resolution cosmological hydro/N-body simulation with star formation, and compare the measured merger timescales with theoretical predictions based on the Chandrasekhar formula. In contrast to Navarro et al., our numerical results indicate, that the commonly used equation for the merger timescale given by Lacey and Cole, systematically underestimates the merger timescales for minor mergers and overestimates those for major mergers. This behavior is partly explained by the poor performance of their expression for the Coulomb logarithm, \ln (m_pri/m_sat). The two alternative forms \ln (1+m_pri/m_sat) and 1/2\ln [1+(m_pri/m_sat)^2] for the Coulomb logarithm can account for the mass dependence of merger timescale successfully, but both of them underestimate the merger time scale by a factor 2. Since \ln (1+m_pri/m_sat) represents the mass dependence slightly better we adopt this expression for the Coulomb logarithm. Furthermore, we find that the dependence of the merger timescale on the circularity parameter \epsilon is much weaker than the widely adopted power-law \epsilon^{0.78}, whereas 0.94*{\epsilon}^{0.60}+0.60 provides a good match to the data. Based on these findings, we present an accurate and convenient fitting formula for the merger timescale of galaxies in cold dark matter models.

The 44Ti(t1/2 = 59 y) nuclide, an important signature of supernova nucleosynthesis, has recently been observed as live radioactivity by gamma-ray astronomy from the Cas A remnant. We investigate in the laboratory the major 44Ti production reaction, 40Ca(alpha,gamma)44Ti (E_cm = 0.6-1.2 MeV/u), by direct off- line counting of 44Ti nuclei. The yield, significantly higher than inferred from previous experiments, is analyzed in terms of a statistical model using microscopic nuclear inputs. The associated stellar rate has important astrophysical consequences, increasing the calculated supernova 44Ti yield by a factor ~2 over previous estimates and bringing it closer to Cas A observations.