We report the detection of four new hot corino sources, G211.47-19.27S, G208.68-19.20N1, G210.49-19.79W and G192.12-11.10 from a survey study of Planck Galactic Cold Clumps in the Orion Molecular Cloud Complex with the Atacama Compact Array (ACA). Three sources had been identified as low mass Class 0 protostars in the Herschel Orion Protostar Survey (HOPS). One source in the lambda Orionis region is firstly reported as a protostellar core. We have observed abundant complex organic molecules (COMs), primarily methanol but also other oxygen-bearing COMs (in G211.47-19.27S and G208.68-19.20N1) and the molecule of prebiotic interest NH2CHO (in G211.47-19.27S), signifying the presence of hot corinos. While our spatial resolution is not sufficient for resolving most of the molecular emission structure, the large linewidth and high rotational temperature of COMs suggest that they likely reside in the hotter and innermost region immediately surrounding the protostar. In G211.47-19.27S, the D/H ratio of methanol ([CH2DOH]/[CH3OH]) and the 12C/13C ratio of methanol ([CH3OH]/[13CH3OH]) are comparable to those of other hot corinos. Hydrocarbons and long carbon-chain molecules such as c-C3H2 and HCCCN are also detected in the four sources, likely tracing the outer and cooler molecular envelopes.

The VRI passbands light curves of V344 Lac were presented and analyzed by using the latest version of the W-D code. The observed spectrum reveals that V344 Lac is not an A3 type but would be a later F type star according to the yielded temperature. The results of solution show that V344 Lac is an A-subtype contact binary, with a mediate photometric mass ratio of 0.387(0.003) and a mediate contact factor of 44.6(3.0)%. Based on the parallax given by Gaia, the parameters of the components are estimated as: M1 = 1.16 Ms, M2 = 0.45 Ms, R1 = 1.31 Rs, R2 = 0.88 Rs, L1 = 2.512 Ls, L2 = 1.057 Ls. The period investigation indicates that V344 Lac may have an eccentric orbital oscillation, with P3 = 12.4(0.5) yr, A3 = 0.0020(0.0002) d, and e = 0.38(0.16). Analysis shows such oscillation would be caused by a magnetic activity which can be explained by the Applegate mechanism. Meanwhile, according to the value of l3 and the estimated physical parameters of V344 Lac, the mass of the third companion may be 0.79 Ms. This third body could be a wide company.

A statistical survey of spectral anisotropy of space plasma turbulence is performed using five years measurements from MMS in the magnetosheath. By measuring the five-point second-order structure functions of the magnetic field, we have for the first time quantified the three-dimensional anisotropies and scalings at sub-ion-scales ($<$ 100 km). In the local reference frame $(\hat L_{\perp}, \hat l_{\perp}, \hat l_{\parallel})$ defined with respect to local mean magnetic field $\boldsymbol {B}_0$ (Chen et al. 2012), the "statistical eddies" are found to be mostly elongated along $\boldsymbol {B}_0$ and shortened in the direction perpendicular to both $\boldsymbol {B}_0$ and local field fluctuations. From several $d_i$ (ion inertial length) toward $\sim$ 0.05 $d_i$, the ratio between eddies' parallel and perpendicular lengths features a trend of rise then fall, whereas the anisotropy in the perpendicular plane appears scale-invariant. Specifically, the anisotropy relations for the total magnetic field at 0.1-1.0 $d_i$ are obtained as $l_{\parallel} \simeq 2.44 \cdot l_{\perp}^{0.71}$, and $L_{\perp} \simeq 1.58 \cdot l_{\perp}^{1.08}$, respectively. Our results provide new observational evidence to compare with phenomenological models and numerical simulations, which may help to better understand the nature of kinetic scale turbulence.

Using the galaxy catalog built from ELUCID N-body simulation and the semi-analytical galaxy formation model, we have built a mock HI intensity mapping map. We have implemented the Finger-of-God (FoG) effect in the map by considering the galaxy HI gas velocity dispersion. By comparing the HI power spectrum in the redshift space with the measurement from IllustrisTNG simulation, we have found that such FoG effect can explain the discrepancy between current mock map built from N-body simulation and Illustris TNG simulation. Then we built a parameter-free FoG model and a shot-noise model to calculate the HI power spectrum. We found that our model can accurately fit both the monopole and quadrupole moments of the HI matter power spectrum. Our method of building the mock HI intensity map and the parameter-free FoG model will be widely useful for the up-coming 21cm intensity mapping experiments, such as CHIME, Tianlai, BINGO, FAST and SKA. It is also crucial for us to study the non-linear effects in 21cm intensity mapping.

We present the first numerical simulations of gravitational waves (GWs) passing through a potential well generated by a compact object in 3-D space, with a realistic source waveform derived from numerical relativity for the merger of two black holes. Unlike the previous work, our analyses focus on the time-domain, in which the propagation of GWs is a well-posed "initial-value" problem for the hyperbolic equations with rigorous rooting in mathematics and physics. Based on these simulations, we investigate for the first time in realistic 3-D space how the wave nature of GWs affects the speed and waveform of GWs in a potential well. We find that GWs travel faster than the prediction of the Shapiro time-delay in the geometric limit due to the effects of diffraction and wavefront geometry. As the wave speed of GWs is closely related to the locality and wavefront geometry of GWs, which are inherently difficult to be addressed in the frequency-domain, our analyses in the time-domain, therefore, provide the first robust analyses to date on this issue based on solid physics. Moreover, we also investigate, for the first time, the interference between the incident and the scattered waves (the "echoes" of the incident waves). We find that such interference makes the total lensed waveforms dramatically different from those of the original incident ones not only in the amplitude but also in the phase and pattern, especially for signals near the merger of the two back holes.

Methane is one of the simplest stable molecules that is both abundant and widely distributed across space. It is thought to have partial origin from interstellar molecular clouds, which are near the beginning of the star formation cycle. Observational surveys of CH$_4$ ice towards low- and high-mass young stellar objects showed that much of the CH$_4$ is expected to be formed by the hydrogenation of C on dust grains, and that CH$_4$ ice is strongly correlated with solid H$_2$O. Yet, this has not been investigated under controlled laboratory conditions, as carbon-atom chemistry of interstellar ice analogues has not been experimentally realized. In this study, we successfully demonstrate with a C-atom beam implemented in an ultrahigh vacuum apparatus the formation of CH$_4$ ice in two separate co-deposition experiments: C + H on a 10 K surface to mimic CH$_4$ formation right before H$_2$O ice is formed on the dust grain, and C + H + H$_2$O on a 10 K surface to mimic CH$_4$ formed simultaneously with H$_2$O ice. We confirm that CH$_4$ can be formed by the reaction of atomic C and H, and that the CH$_4$ formation rate is 2 times greater when CH$_4$ is formed within a H$_2$O-rich ice. This is in agreement with the observational finding that interstellar CH$_4$ and H$_2$O form together in the polar ice phase, i.e., when C- and H-atoms simultaneously accrete with O-atoms on dust grains. For the first time, the conditions that lead to interstellar CH$_4$ (and CD$_4$) ice formation are reported, and can be incorporated into astrochemical models to further constrain CH$_4$ chemistry in the interstellar medium and in other regions where CH$_4$ is inherited.

To study the high-transition dense-gas tracers and their relationships to the star formation of the inner $\sim$ 2 kpc circumnuclear region of NGC253, we present HCN $J=4-3$ and HCO$^+ J=4-3$ maps obtained with the James Clerk Maxwell Telescope (JCMT). With the spatially resolved data, we compute the concentration indices $r_{90}/r_{50}$ for the different tracers. HCN and HCO$^+$ 4-3 emission features tend to be centrally concentrated, which is in contrast to the shallower distribution of CO 1-0 and the stellar component. The dense-gas fraction ($f_\text{dense}$, traced by the velocity-integrated-intensity ratios of HCN/CO and HCO$^+$/CO) and the ratio $R_\text{31}$ (CO 3-2/1-0) decline towards larger galactocentric distances, but increase with higher SFR surface density. The radial variation and the large scatter of $f_\text{dense}$ and $R_\text{31}$ imply distinct physical conditions in different regions of the galactic disc. The relationships of $f_\text{dense}$ versus $\Sigma_\text{stellar}$, and SFE$_\text{dense}$ versus $\Sigma_\text{stellar}$ are explored. SFE$_\text{dense}$ increases with higher $\Sigma_\text{stellar}$ in this galaxy, which is inconsistent with previous work that used HCN 1-0 data. This implies that existing stellar components might have different effects on the high-$J$ HCN and HCO$^+$ than their low-$J$ emission. We also find that SFE$_\text{dense}$ seems to be decreasing with higher $f_\text{dense}$, which is consistent with previous works, and it suggests that the ability of the dense gas to form stars diminishes when the average density of the gas increases. This is expected in a scenario where only the regions with high-density contrast collapse and form stars.

We investigate how a property of a galaxy correlates most tightly with a property of its host dark matter halo, using state-of-the-art hydrodynamical simulations of galaxy formation EAGLE, Illustris, and IllustrisTNG. Unlike most of the previous work, our analyses focus on all types of galaxies, including both central and satellite galaxies. We find that the stellar mass of a galaxy at the epoch of the peak circular velocity with an evolution correction gives the tightest such correlation to the peak circular velocity $V_{\rm peak}$ of the galaxy's underling dark matter halo. The evolution of galaxy stellar mass reduces rather than increases scatter in such a relation. We also find that one major source of scatter comes from star stripping due to the strong interactions between galaxies. Even though, we show that the size of scatter predicted by hydrodynamical simulations has a negligible impact on the clustering of dense $V_{\rm peak}$-selected subhalo from simulations, which suggests that even the simplest subhalo abundance matching (SHAM), without scatter and any additional free parameter, can provide a robust prediction of galaxy clustering that can agree impressively well with the observations from the SDSS main galaxy survey.

Offsets of molecular line emission peaks from continuum peaks are very common but frequently difficult to explain with a single spherical cloud chemical model. We propose that the spatial projection effects of an irregular three dimensional (3D) cloud structure can be a solution. This work shows that the idea can be successfully applied to the Planck cold clump G224.4-0.6 by approximating it with four individual spherically symmetric cloud cores whose chemical patterns overlap with each other to produce observable line maps. With the empirical physical structures inferred from the observation data of this clump and a gas-grain chemical model, the four cores can satisfactorily reproduce its 850 $\mu$m continuum map and the diverse peak offsets of CCS, HC$_3$N and N$_2$H$^+$ simultaneously at chemical ages of about $8\times 10^5\sim 3\times 10^6$ yrs. The 3D projection effects on chemistry has the potential to explain such asymmetrical distributions of chemicals in many other molecular clouds.

To deepen our understanding of the chemical properties of the Planck Galactic Cold Clump (PGCC) G168.72-15.48, we performed observations of nine molecular species, namely, \ce{c-C3H}, \ce{H2CO}, \ce{HC5N}, \ce{HC7N}, \ce{SO}, \ce{CCH}, \ce{N2H+}, \ce{CH3OH}, and \ce{CH3CCH}, toward two dense cores in PGCC G168.72-15.48 using the Tianma Radio Telescope and Purple Mountain Observatory Telescope. We detected \ce{c-C3H}, \ce{H2CO}, \ce{HC5N}, \ce{N2H+}, \ce{CCH}, and \ce{CH3OH} in both G168-H1 and G168-H2 cores, whereas \ce{HC7N} and \ce{CH3CCH} were detected only in G168-H1 and SO was detected only in G168-H2. Mapping observations reveal that the \ce{CCH}, \ce{N2H+}, \ce{CH3OH}, and \ce{CH3CCH} emissions are well coupled with the dust emission in G168-H1. Additionally, \ce{N2H+} exhibits an exceptionally weak emission in the denser and more evolved G168-H2 core, which may be attributed to the \ce{N2H+} depletion. We suggest that the \ce{N2H+} depletion in G168-H2 is dominated by \ce{N2} depletion, rather than the destruction by CO. The local thermodynamic equilibrium calculations indicate that the carbon-chain molecules of \ce{CCH}, \ce{HC5N}, \ce{HC7N}, and \ce{CH3CCH} are more abundant in the younger G168-H1 core. We found that starless core G168-H1 may have the properties of cold dark clouds based on its abundances of carbon-chain molecules. While, the prestellar core G168-H2 exhibits lower carbon-chain molecular abundances than the general cold dark clouds. With our gas-grain astrochemical model calculations, we attribute the observed chemical differences between G168-H1 and G168-H2 to their different gas densities and different evolutionary stages.

We present a CO(3-2) survey of selected regions in the M31 disc as part of the JCMT large programme, HARP and SCUBA-2 High-Resolution Terahertz Andromeda Galaxy Survey (HASHTAG). The 12 CO(3-2) fields in this survey cover a total area of 60 square arcminutes, spanning a deprojected radial range of 2 - 14 kpc across the M31 disc. Combining these observations with existing IRAM 30m CO(1-0) observations and JCMT CO(3-2) maps of the nuclear region of M31, as well as dust temperature and star formation rate surface density maps, we are able to explore the radial distribution of the CO(3-2)/CO(1-0) integrated intensity ratio (R31) and its relationship with dust temperature and star formation. We find that the value of R31 between 2 - 9 kpc galactocentric radius is 0.14, significantly lower than what is seen in the nuclear ring at ~1 kpc (R31 ~ 0.8), only to rise again to 0.27 for the fields centred on the 10 kpc star forming ring. We also found that R31 is positively correlated with dust temperature, with Spearman's rank correlation coefficient $\rho$ = 0.55. The correlation between star formation rate surface density and CO(3--2) intensity is much stronger than with CO(1-0), with $\rho$ = 0.54 compared to -0.05, suggesting that the CO(3-2) line traces warmer and denser star forming gas better. We also find that R31 correlates well with star formation rate surface density, with $\rho$ = 0.69.

We investigate the wave effects of gravitational waves (GWs) using numerical simulations with the finite element method (FEM) based on the publicly available code {\it deal.ii}. We robustly test our code using a point source monochromatic spherical wave. We examine not only the waveform observed by a local observer but also the global energy conservation of the waves. We find that our numerical results agree very well with the analytical predictions. Based on our code, we study the scattering of GWs by compact objects. Using monochromatic waves as the input source, we find that if the wavelength of GWs is much larger than the Schwarzschild radius of the compact object, the amplitude of the total scattered GWs does not change appreciably due to the strong diffraction effect, for an observer far away from the scatterer. This finding is consistent with the results reported in the literature. However, we also find that, near the scatterer, not only the amplitude of the scattered waves is very large, comparable to that of the incident waves, but also the phase of the GWs changes significantly due to the interference between the scattered and incident waves. As the evolution of the phase of GWs plays a crucial role in the matched filtering technique in extracting GW signals from the noisy background, our findings suggest that wave effects should be taken into account in the data analysis in the future low-frequency GW experiments, if GWs are scattered by nearby compact objects in our local environment.

We started a multi-scale analysis of G202.3+2.5, an intertwined filamentary region of Monoceros OB1. In Paper I, we examined the distributions of dense cores and protostars and found enhanced star formation (SF) activity in the junction region of the filaments. In this second paper, we aim to unveil the connections between the core and filament evolutions, and between the filament dynamics and the global evolution of the cloud. We characterise the gas dynamics and energy balance using Herschel and WISE observations and molecular tracers observed with the IRAM 30m and TRAO 14m telescopes. The velocity field of the cloud is examined and velocity-coherent structures are put in perspective with the cloud environment. Two main velocity components (VCs) are revealed, well separated in the north and merged around the location of intense N2H+ emission where Paper I found the peak of SF activity. The relative position of the two VCs along the sightline, and the velocity gradient in N2H+ emission imply that the VCs have been undergoing collision for ~10^5 yrs. The dense gas where N2H+ is detected is interpreted as the compressed region between the two filaments, which corresponds to a high mass inflow rate of ~1e-3 Msun/yr and possibly leads to an increase in its SF efficiency. We identify a protostar in the junction region that possibly powers two crossed intermittent outflows. We show that the HII region around the nearby cluster NCG 2264 is still expanding and its role in the collision is examined. However, we cannot rule out the idea that the collision arises mostly from the global collapse of the cloud. The (sub-)filament-scale observables examined in this paper reveal a collision between G202.3+2.5 sub-structures and its probable role in feeding the cores in the junction region. One must now characterise the cloud morphology, its fragmentation, and magnetic field, all at high resolution.

Current theories and models attempt to explain star formation globally, from core scales to giant molecular cloud scales. A multi-scale observational characterisation of an entire molecular complex is necessary to constrain them. We investigate star formation in G202.3+2.5, a ~10x3 pc sub-region of the Monoceros OB1 cloud with a complex morphology harbouring interconnected filamentary structures. We aim to connect the evolution of cores and filaments in G202.3+2.5 with the global evolution of the cloud and to identify the engines of the cloud dynamics. In this first paper, the star formation activity is evaluated by surveying the distributions of dense cores and protostars, and their evolutionary state, as characterised using both infrared observations from the Herschel and WISE telescopes and molecular line observations with the IRAM 30-m telescope. We find ongoing star formation in the whole cloud, with a local peak in star formation activity around the centre of G202.3+2.5 where a chain of massive cores (10-50 Msun) forms a massive ridge (>150 Msun). All evolutionary stages from starless cores to Class II protostars are found in G202.3+2.5, including a possibly starless, large column density (8x10^{22} cm-2), and massive (52 Msun) core. All the core-scale observables examined in this paper point to an enhanced star formation activity centred on the junction between the three main branches of the ramified structure of G202.3+2.5. This suggests that the increased star-formation activity results from the convergence of these branches. To further investigate the origin of this enhancement, it is now necessary to extend the analysis to larger scales, in order to examine the relationship between cores, filaments and their environment. We address these points through the analysis of the dynamics of G202.3+2.5 in a joint paper.

Temporal variation of millimeter lines is a new direction of research for evolved stars. It has the potential to probe the dynamical wind launching processes from time dimension. We report here the first ALMA (Atacama Large Millimeter Array) results that cover 817 days of an on-going monitoring of 1.1 mm lines in the archetypal carbon star IRC +10216. The monitoring is done with the compact 7-m array (ACA) and in infrared with a 1.25 m telescope in Crimea. A high sensitivity of the cumulative spectra covering a total of ~7.2 GHz between 250 - 270 GHz range has allowed us to detect about 148 known transitions of 20 molecules, together with more of their isotopologues, and 81 unidentified lines. An overview of the variabilities of all detected line features are presented in spectral plots. Although a handful of lines are found to be very possibly stable in time, most other lines are varying either roughly in phase or in anti-correlation with the near-infrared light. Several lines have their variations in the ALMA data coincident with existing single dish monitoring results, while several others do not, which requires an yet unknown mechanism in the circumstellar envelop to explain.

We present the highest resolution single-dish submillimetre observations of the detached shell source U Antliae to date. The observations were obtained at $450~\micron$ and $850~\micron$ with SCUBA-2 instrument on the James Clerk Maxwell Telescope as part of the Nearby Evolved Stars Survey. The emission at $850~\micron$ peaks at $40\arcsec$ with hints of a second peak seen at $\sim 20\arcsec$. The emission can be traced out to a radius of $56\arcsec$ at a $3\sigma$ level. The outer peak observed at $850~\micron$ aligns well with the peak observed at Herschel/PACS wavelengths. With the help of spectral energy distribution fitting and radiative transfer calculations of multiple-shell models for the circumstellar envelope, we explore the various shell structures and the variation of grain sizes along the in the circumstellar envelope. We determine a total shell dust mass of $(2.0 \pm 0.3) \times 10^{-5}$ M$_{\odot}$ and established that the thermal pulse which gave rise to the detached shell occurred 3500 $\pm$ 500 years ago.

A grand-challenge problem at the forefront of physics is to understand how energy is transported and transformed in plasmas. This fundamental research priority encapsulates the conversion of plasma-flow and electromagnetic energies into particle energy, either as heat or some other form of energisation. The smallest characteristic scales, at which electron dynamics determines the plasma behaviour, are the next frontier in space and astrophysical plasma research. The analysis of astrophysical processes at these scales lies at the heart of the field of electron-astrophysics. Electron scales are the ultimate bottleneck for dissipation of plasma turbulence, which is a fundamental process not understood in the electron-kinetic regime. Since electrons are the most numerous and most mobile plasma species in fully ionised plasmas and are strongly guided by the magnetic field, their thermal properties couple very efficiently to global plasma dynamics and thermodynamics.

This white paper gives a brief summary of Galactic continuum surveys with the next generation 850 micron camera at the James Clerk Maxwell Telescope (JCMT) in the next decade. This new camera will have mapping speeds at least 10-20 times faster than the present SCUBA-2 camera, and will enable deep (<10 mJy/beam) and extensive continuum surveys such as a wider (-5 deg.<l<240 deg., |b|<2-5 deg.) Galactic Plane survey, a larger Gould Belt survey, and a follow-up survey of all Planck compact objects visible from the northern hemisphere. These surveys will provide a complete census of molecular clouds, filaments, and dense cores across the Galaxy, vital for studying star formation in various environments.

We propose a novel approach to accurately pin down the systematics due to the peculiar velocities of galaxies in measuring the Hubble constant from nearby galaxies in current and future gravitational-wave (GW) standard-siren experiments. Given the precision that future GW standard-siren experiments aim to achieve, the peculiar velocities of nearby galaxies will be a major source of uncertainty. Unlike the conventional backward reconstruction that requires additional redshift-independent distance indicators to recover the peculiar velocity field, we forwardly model the peculiar velocity field by using a high-fidelity mock galaxy catalog built from high-resolution dark matter only (DMO) N-body simulations with a physically motivated subhalo abundance matching technique without introducing any free parameters. Our mock galaxy catalog can impressively well reproduce the observed spectroscopic redshift space distortions (RSDs) in highly non-linear regimes down to very small scales, which is a robust test of the velocity field of our mock galaxy catalog. Based on this mock galaxy catalog, we accurately, for the first time, measure the peculiar velocity probability distributions for the SDSS main galaxy samples. We find that the systematics induced by the peculiar velocities of SDSS like galaxies on the measured Hubble constant can be reduced to below $1\%$($1\sigma$) for GW host galaxies with a Hubble flow redshift just above $0.13$, a distance that can be well probed by future GW experiments and galaxy surveys.

Complex organic molecules (COMs) have been detected in the gas-phase in cold and lightless molecular cores. Recent solid-state laboratory experiments have provided strong evidence that COMs can be formed on icy grains through 'non-energetic' processes. In this contribution, we show that propanal and 1-propanol can be formed in this way at the low temperature of 10 K. Propanal has already been detected in space. 1-propanol is an astrobiologically relevant molecule, as it is a primary alcohol, and has not been astronomically detected. Propanal is the major product formed in the C2H2 + CO + H experiment, and 1-propanol is detected in the subsequent propanal + H experiment. The results are published in Qasim et al. (2019c). ALMA observations towards IRAS 16293-2422B are discussed and provide a 1-propanol:propanal upper limit of < 0.35 - 0.55, which are complemented by computationally-derived activation barriers in addition to the performed laboratory experiments.

1-propanol (CH3CH2CH2OH) is a three carbon-bearing representative of primary linear alcohols that may have its origin in the cold dark cores in interstellar space. To test this, we investigated in the laboratory whether 1-propanol ice can be formed along pathways possibly relevant to the prestellar core phase. We aim to show in a two-step approach that 1-propanol can be formed through reaction steps that are expected to take place during the heavy CO freeze-out stage by adding C2H2 into the CO + H hydrogenation network via the formation of propanal (CH3CH2CHO) as an intermediate and its subsequent hydrogenation. Temperature programmed desorption-quadrupole mass spectrometry (TPD-QMS) is used to identify the newly formed propanal and 1-propanol. Reflection absorption infrared spectroscopy (RAIRS) is used as a complementary diagnostic tool. The mechanisms that can contribute to the formation of solid-state propanal and 1-propanol, as well as other organic compounds, during the heavy CO freeze-out stage are constrained by both laboratory experiments and theoretical calculations. Here it is shown that recombination of HCO radicals, formed upon CO hydrogenation, with radicals formed upon C2H2 processing - H2CCH and H3CCH2 - offers possible reaction pathways to solid-state propanal and 1-propanol formation. This extends the already important role of the CO hydrogenation chain in the formation of larger COMs (complex organic molecules). The results are used to compare with ALMA observations. The resulting 1-propanol:propanal ratio concludes an upper limit of < 0:35-0:55, which is complemented by computationally-derived activation barriers in addition to the experimental results.

Formic acid (HCOOH) and carbon dioxide (CO2) are simple species that have been detected in the interstellar medium. The solid-state formation pathways of these species under experimental conditions relevant to prestellar cores are primarily based off of weak infrared transitions of the HOCO complex and usually pertain to the H2O-rich ice phase, and therefore more experimental data are desired. In this article, we present a new and additional solid-state reaction pathway that can form HCOOH and CO2 ice at 10 K 'non-energetically' in the laboratory under conditions related to the "heavy" CO freeze-out stage in dense interstellar clouds, i.e., by the hydrogenation of an H2CO:O2 ice mixture. This pathway is used to piece together the HCOOH and CO2 formation routes when H2CO or CO reacts with H and OH radicals. Temperature programmed desorption - quadrupole mass spectrometry (TPD-QMS) is used to confirm the formation and pathways of newly synthesized ice species as well as to provide information on relative molecular abundances. Reflection absorption infrared spectroscopy (RAIRS) is additionally employed to characterize reaction products and determine relative molecular abundances. We find that for the conditions investigated in conjunction with theoretical results from the literature, H+HOCO and HCO+OH lead to the formation of HCOOH ice in our experiments. Which reaction is more dominant can be determined if the H+HOCO branching ratio is more constrained by computational simulations, as the HCOOH:CO2 abundance ratio is experimentally measured to be around 1.8:1. H+HOCO is more likely than OH+CO (without HOCO formation) to form CO2. Isotope experiments presented here further validate that H+HOCO is the dominant route for HCOOH ice formation in a CO-rich CO:O2 ice mixture that is hydrogenated. These data will help in the search and positive identification of HCOOH ice in prestellar cores.

We present a novel fitting formula for the halo concentration enhancement in chameleon $f(R)$ gravity relative to General Relativity (GR). The formula is derived by employing a large set of $N$-body simulations of the Hu-Sawicki $f(R)$ model which cover a wide range of model and cosmological parameters, resolutions and simulation box sizes. The complicated dependence of the concentration on halo mass $M$, redshift $z$, and the $f(R)$ and cosmological parameters can be combined into a simpler form that depends only on a rescaled mass $M/10^{p_2}$, with $p_2\equiv1.5\log_{10}\left[|{\bar{f}_R(z)}|/(1+z)\right]+21.64$ and $\bar{f}_R(z)$ the background scalar field at $z$, irrespective of the $f(R)$ model parameter. Our fitting formula can describe the concentration enhancement well for redshifts $z\leq3$, nearly 7 orders of magnitude in $M/10^{p_2}$ and five decades in halo mass. This is part of a series of works which aims to provide a general framework for self-consistent and unbiased tests of gravity using present and upcoming galaxy cluster surveys. The fitting formula, which is the first quantitative model for the concentration enhancement due to chameleon type modified gravity, is an important part in this framework and will allow continuous exploration of the parameter space. It can also be used to model other statistics such as the matter power spectrum.

A number of recent experimental studies have shown that solid-state complex organic molecules (COMs) can form under conditions that are relevant to the CO freeze-out stage in dense clouds. In this work, we show that alcohols can be formed well before the CO freeze-out stage (i.e., in the H2O-rich ice phase). This joint experimental and computational investigation shows that the isomers, n- and i-propanol (H3CCH2CH2OH and H3CCHOHCH3) and n- and i-propenol (H3CCHCHOH and H3CCOHCH2), can be formed in radical-addition reactions starting from propyne (H3CCCH) + OH at the low temperature of 10 K, where H3CCCH is one of the simplest representatives of stable carbon chains already identified in the interstellar medium (ISM). The resulting average abundance ratio of 1:1 for n-propanol:i-propanol is aligned with the conclusions from the computational work that the geometric orientation of strongly interacting species is influential to the extent of which 'mechanism' is participating, and that an assortment of geometries leads to an averaged-out effect. Three isomers of propanediol are also tentatively identified in the experiments. It is also shown that propene and propane (H3CCHCH2 and H3CCH2CH3) are formed from the hydrogenation of H3CCCH. Computationally-derived activation barriers give additional insight into what types of reactions and mechanisms are more likely to occur in the laboratory and in the ISM. Our findings not only suggest that the alcohols studied here share common chemical pathways and therefore can show up simultaneously in astronomical surveys, but also that their extended counterparts that derive from polyynes containing H3C(CC)nH structures may exist in the ISM. Such larger species, like fatty alcohols, are the possible constituents of simple lipids that primitive cell membranes on the early Earth are thought to be partially composed of.

The need to characterize ices coating dust grains in dense interstellar clouds arises from the importance of ice morphology in facilitating the diffusion and storage of radicals and reaction products in ices, a well-known place for the formation of complex molecules. Yet, there is considerable uncertainty about the structure of ISM ices, their ability to store volatiles and under what conditions. We measured the infrared absorption spectra of CO on the pore surface of porous amorphous solid water (ASW), and quantified the effective pore surface area of ASW. Additionally, we present results obtained from a Monte Carlo model of ASW in which the morphology of the ice is directly visualized and quantified. We found that 200 ML of ASW annealed to 20 K has a total pore surface area that is equivalent to 46 ML. This surface area decreases linearly with temperature to about 120 K. We also found that (1) dangling OH bonds only exist on the surface of pores; (2) almost all of the pores in the ASW are connected to the vacuum--ice interface, and are accessible for adsorption of volatiles from the gas phase; there are few closed cavities inside ASW at least up to a thickness of 200 ML; (3) the total pore surface area is proportional to the total 3-coordinated water molecules in the ASW in the temperature range 60--120 K. We also discuss the implications on the structure of ASW and surface reactions in the ice mantle in dense clouds.

Ellerman bombs (EBs) are small-scale intense brightenings in H$\alpha$ wing images, which are generally believed to be signatures of magnetic reconnection events around the temperature minimum region of the solar atmosphere. They have a flame-like morphology when observed near the solar limb. Recent observations from the Interface Region Imaging Spectrograph (IRIS) reveal another type of small-scale reconnection events, termed UV bursts, in the lower atmosphere. Though previous observations have shown a clear coincidence of some UV bursts and EBs, the exact relationship between these two phenomena is still under debate. We investigate the spatial and temporal relationship between flame-like EBs and UV bursts using joint near-limb observations between the 1.6--meter Goode Solar Telescope (GST) and IRIS. In total 161 EBs have been identified from the GST observations, and 20 of them reveal signatures of UV bursts in the IRIS images. Interestingly, we find that these UV bursts have a tendency to appear at the upper parts of their associated flame-like EBs. The intensity variations of most EB-related UV bursts and their corresponding EBs match well. Our results suggest that these UV bursts and EBs are likely formed at different heights during a common reconnection process.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

Solar ultraviolet (UV) bursts are small-scale compact brightenings in transition region images. The spectral profiles of transition region lines in these bursts are significantly enhanced and broadened, often with chromospheric absorption lines such as Ni~{\sc{ii}} 1335.203 and 1393.330 {\AA} superimposed. We investigate the properties of several UV bursts using a coordinated observation of the Interface Region Imaging Spectrograph (IRIS), Solar Dynamics Observatory (SDO), and \textit{Hinode} on 2015 February 7. We have identified 12 UV bursts, and 11 of them reveal small blueshifts of the Ni~{\sc{ii}} absorption lines. However, the Ni~{\sc{ii}} lines in one UV burst exhibit obvious redshifts of $\sim$20 km s$^{-1}$, which appear to be related to the cold plasma downflows observed in the IRIS slit-jaw images. We also examine the three-dimensional magnetic field topology using a magnetohydrostatic model, and find that some UV bursts are associated with magnetic null points or bald patches. In addition, we find that these UV bursts reveal no obvious coronal signatures from the observations of the Atmospheric Imaging Assembly (AIA) on board SDO and the EUV Imaging Spectrometer (EIS) on board \textit{Hinode}.

Kinetic Alfv\'{e}n waves (KAWs) are the short-wavelength extension of the MHD Alfv\'{e}n-wave branch in the case of highly-oblique propagation with respect to the background magnetic field. Observations of space plasma show that small-scale turbulence is mainly KAW-like. We apply two theoretical approaches, collisional two-fluid theory and collisionless kinetic theory, to obtain predictions for the KAW polarizations depending on $\beta_\mathrm{p}$ (the ratio of the proton thermal pressure to the magnetic pressure) at the ion gyroscale in terms of fluctuations in density, bulk velocity, and pressure. We perform a wavelet analysis of MMS magnetosheath measurements and compare the observations with both theories. We find that the two-fluid theory predicts the observations better than kinetic theory, suggesting that the small-scale KAW-like fluctuations exhibit a fluid-like behavior in the magnetosheath although the plasma is weakly collisional. We also present predictions for the KAW polarizations in the inner heliosphere that are testable with Parker Solar Probe and Solar Orbiter.

The Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Survey (BGS) will be a survey of bright, low redshift galaxies, which is planned to cover an area of ~14,000 sq deg in 3 passes. Each pass will cover the survey area with ~2000 pointings, each of area ~8 sq deg. The BGS is currently proposed to consist of a bright high priority sample to an r-band magnitude limit r ~ 19.5, with a fainter low priority sample to r ~ 20. The geometry of the DESI fibre positioners in the focal plane of the telescope affects the completeness of the survey, and has a non-trivial impact on clustering measurements. Using a BGS mock catalogue, we show that completeness due to fibre assignment primarily depends on the surface density of galaxies. Completeness is high (>95%) in low density regions, but very low (<10%) in the centre of massive clusters. We apply the pair inverse probability (PIP) weighting correction to clustering measurements from a BGS mock which has been through the fibre assignment algorithm. This method is only unbiased if it is possible to observe every galaxy pair. To facilitate this, we randomly promote a small fraction of the fainter sample to be high priority, and dither the set of tile positions by a small angle. We show that inverse pair weighting combined with angular upweighting provides an unbiased correction to galaxy clustering measurements for the complete 3 pass survey, and also after 1 pass, which is highly incomplete.

Four-color charge-coupled device (CCD) light curves in $B$, $V$, $Rc$ and $Ic$ bands of the total-eclipsing binary system, V1853 Ori, are presented. By comparing our light curves with those published by previous investigators, it is detected that the O'Connell effect on the light curves is disappeared. By analyzing those multi-color light curves with the Wilson-Devinney code (W-D code), it is discovered that V1853 Ori is an A-type intermediate-contact binary with a degree of contact factor of $f=33.3\%(3.7\%)$ and a mass ratio of $q=0.1896(0.0013)$. Combining our 10 new determined times of light minima together with the others published in the literature, the period changes of the system is investigated. We found that the general trend of the observed-calculated $(O-C)$ curve shows a downward parabolic variation that corresponds to a long-term decrease in the orbital period with a rate of $dP/dt=-1.96(0.46)\times{10^{-7}}$ d yr$^{-1}$. The long-term period decrease could be explained by mass transfer from the more-massive component to the less-massive one. By combining our photometric solutions with the Gaia DR 2 data, absolute parameters were derived as $M_{1}$ = 1.20 M$_{\odot}$, $M_{2}$ = 0.23 M$_{\odot}$, $R_{1}$ = 1.36 R$_{\odot}$, and $R_{2}$ = 0.66 R$_{\odot}$. The long-term period decrease and the intermediate-contact configuration suggest that V1853 Ori will evolve into a high fill-out overcontact binary.

Analyses of infrared signatures of CO$_2$ in water dominated ices in the ISM can give information on the physical state of CO$_2$ in icy grains and on the thermal history of the ices themselves. In many sources, CO$_2$ was found in the `pure' crystalline form, as signatured by the splitting in the bending mode absorption profile. To a large extent, pure CO$_2$ is likely to have formed from segregation of CO$_2$ from a CO$_2$:H$_2$O mixture during thermal processing. Previous laboratory studies quantified the temperature dependence of segregation, but no systematic measurement of the concentration dependence of segregation is available. In this study, we measured both the temperature dependence and concentration dependence of CO$_2$ segregation in CO$_2$:H$_2$O mixtures, and found that no pure crystalline CO$_2$ forms if the CO$_2$:H$_2$O ratio is less than 23%. Therefore the segregation of CO$_2$ is not always a good thermal tracer of the ice mantle. We found that the position and width of the broad component of the asymmetric stretching vibrational mode of $^{13}$CO$_2$ change linearly with the temperature of CO$_2$:H$_2$O mixtures, but are insensitive to the concentration of CO$_2$. We recommend using this mode, which will be observable towards low mass protostellar envelopes and dense clouds with the James Webb Space Telescope, to trace the thermal history of the ice mantle, especially when segregated CO$_2$ is unavailable. We used the laboratory measured $^{13}$CO$_2$ profile to analyze the ISO-SWS observations of ice mantles towards Young Stellar Objects, and the astrophysical implications are discussed.

The recent discovery of gravitational waves marks the culmination of a sequence of successful tests of the general theory of relativity (GR) since its formulation in 1915. Yet these tests remain confined to the scale of stellar systems or the strong gravity regime. A departure from GR on larger, cosmological scales has been advocated by the proponents of modified gravity theories as an alternative to the Cosmological Constant to account for the observed cosmic expansion history. While indistinguishable in these terms by construction, such models on the other hand yield distinct values for the linear growth rate of density perturbations and, as a consequence, for the associated galaxy peculiar velocity field. Measurements of the resulting anisotropy of galaxy clustering, when spectroscopic redshifts are used to derive distances, have thus been proposed as a powerful probe of the validity of GR on cosmological scales. However, despite significant effort in modelling such redshift space distortions, systematic errors remain comparable to current statistical uncertainties. Here, we present the results of a different forward-modelling approach, which fully exploits the sensitivity of the galaxy velocity field to modifications of GR. We use state-of-the-art, high-resolution N-body simulations of a standard GR and a compelling f(R) model, one of GR's simplest variants, to build simulated catalogues of stellar-mass-selected galaxies through a robust match to the Sloan Digital Sky Survey observations. We find that, well within the uncertainty of this technique, f(R) fails to reproduce the observed redshift-space clustering on scales 1-10 Mpc/h. Instead, the standard LCDM GR model agrees impressively well with the data. This result provides a strong confirmation, on cosmological scales, of the robustness of Einstein's general theory of relativity.

JINGLE is a new JCMT legacy survey designed to systematically study the cold interstellar medium of galaxies in the local Universe. As part of the survey we perform 850um continuum measurements with SCUBA-2 for a representative sample of 193 Herschel-selected galaxies with M*>10^9Msun, as well as integrated CO(2-1) line fluxes with RxA3m for a subset of 90 of these galaxies. The sample is selected from fields covered by the Herschel-ATLAS survey that are also targeted by the MaNGA optical integral-field spectroscopic survey. The new JCMT observations combined with the multi-wavelength ancillary data will allow for the robust characterization of the properties of dust in the nearby Universe, and the benchmarking of scaling relations between dust, gas, and global galaxy properties. In this paper we give an overview of the survey objectives and details about the sample selection and JCMT observations, present a consistent 30 band UV-to-FIR photometric catalog with derived properties, and introduce the JINGLE Main Data Release (MDR). Science highlights include the non-linearity of the relation between 850um luminosity and CO line luminosity, and the serendipitous discovery of candidate z>6 galaxies.

We present unprecedented high-resolution TiO images and Fe I 1565 nm spectropolarimetric data of two light bridges taken by the 1.6-m Goode Solar Telescope at Big Bear Solar Observatory. In the first light bridge (LB1), we find striking knot-like dark structures within the central dark lane. Many dark knots show migration away from the penumbra along the light bridge. The sizes, intensity depressions and apparent speeds of their proper motion along the light bridges of 33 dark knots identified from the TiO images are mainly in the ranges of 80$\sim$200~km, 30\%$\sim$50\%, and 0.3$\sim$1.2~km~s$^{-1}$, respectively. In the second light bridge (LB2), a faint central dark lane and striking transverse intergranular lanes were observed. These intergranular lanes have sizes and intensity depressions comparable to those of the dark knots in LB1, and also migrate away from the penumbra at similar speeds. Our observations reveal that LB2 is made up of a chain of evolving convection cells, as indicated by patches of blue shift surrounded by narrow lanes of red shift. The central dark lane generally corresponds to blueshifts, supporting the previous suggestion of central dark lanes being the top parts of convection upflows. In contrast, the intergranular lanes are associated with redshifts and located at two sides of each convection cell. The magnetic fields are stronger in intergranular lanes than in the central dark lane. These results suggest that these intergranular lanes are manifestations of convergent convective downflows in the light bridge. We also provide evidence that the dark knots observed in LB1 may have a similar origin.

The diffusion of atoms and molecules in ices covering dust grains in dense clouds in interstellar space is an important but poorly characterized step in the formation of complex molecules in space. Here we report the measurement of diffusion of simple molecules in amorphous solid water (ASW), an analog of interstellar ices, which are amorphous and made mostly of water molecules. The new approach that we used relies on measuring in situ the change in band strength and position of mid-infrared features of OH dangling bonds as molecules move through pores and channels of ASW. We obtained the Arrhenius pre-exponents and activation energies for diffusion of CO, O$_2$, N$_2$, CH$_4$, and Ar in ASW. The diffusion energy barrier of H$_2$ and D$_2$ were also measured, but only upper limits were obtained. These values constitute the first comprehensive set of diffusion parameters of simple molecules on the pore surface of ASW, and can be used in simulations of the chemical evolution of ISM environments, thus replacing unsupported estimates. We also present a set of argon temperature programmed desorption experiments to determine the desorption energy distribution of argon on non-porous ASW.

We present HCN J=4-3 and HCO^+ J=4-3 maps of six nearby star-forming galaxies, NGC 253, NGC 1068, IC 342, M82, M83, and NGC 6946, obtained with the James Clerk Maxwell Telescope as part of the MALATANG survey. All galaxies were mapped in the central 2 arcmin $\times$ 2 arcmin region at 14 arcsec (FWHM) resolution (corresponding to linear scales of ~ 0.2-1.0 kpc). The L_IR-L'_dense relation, where the dense gas is traced by the HCN J=4-3 and the HCO^+ J=4-3 emission, measured in our sample of spatially-resolved galaxies is found to follow the linear correlation established globally in galaxies within the scatter. We find that the luminosity ratio, L_IR/L'_dense, shows systematic variations with L_IR within individual spatially resolved galaxies, whereas the galaxy-integrated ratios vary little. A rising trend is also found between L_IR/L'_dense ratio and the warm-dust temperature gauged by the 70 \mu m/100 \mu m flux ratio. We find the luminosity ratios of IR/HCN(4-3) and IR/HCO^+(4-3), which can be taken as a proxy for the efficiency of star formation in the dense molecular gas (SFE_dense), appears to be nearly independent of the dense-gas fraction (f_dense) for our sample of galaxies. The SFE of the total molecular gas (SFE_mol) is found to increase substantially with f_dense when combining our data with that on local (ultra)luminous infrared galaxies and high-z quasars. The mean L'_HCN(4-3)/L'_HCO^+(4-3) line ratio measured for the six targeted galaxies is 0.9+/-0.6. No significant correlation is found for the L'_HCN(4-3)/L'_HCO^+(4-3) ratio with the SFR as traced by L_IR, nor with the warm-dust temperature, for the different populations of galaxies.

Using isochronous mass spectrometry at the experimental storage ring CSRe in Lanzhou, the masses of $^{82}$Zr and $^{84}$Nb were measured for the first time with an uncertainty of $\sim 10$ keV, and the masses of $^{79}$Y, $^{81}$Zr, and $^{83}$Nb were re-determined with a higher precision. %The latter differ significantly from their literature values. The latter are significantly less bound than their literature values. Our new and accurate masses remove the irregularities of the mass surface in this region of the nuclear chart. Our results do not support the predicted island of pronounced low $\alpha$ separation energies for neutron-deficient Mo and Tc isotopes, making the formation of Zr-Nb cycle in the $rp$-process unlikely. The new proton separation energy of $^{83}$Nb was determined to be 490(400)~keV smaller than that in the Atomic Mass Evaluation 2012. This partly removes the overproduction of the $p$-nucleus $^{84}$Sr relative to the neutron-deficient molybdenum isotopes in the previous $\nu p$-process simulations.

We present an investigation of a coronal cavity observed above the western limb in the coronal red line Fe X 6374 {\AA} using a telescope of Peking University and in the green line Fe XIV 5303 {\AA} using a telescope of Yunnan Observatories, Chinese Academy of Sciences during the total solar eclipse on 2017 August 21. A series of magnetic field models are constructed based on the magnetograms taken by the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO) one week before the eclipse. The model field lines are then compared with coronal structures seen in images taken by the Atmospheric Imaging Assembly on board SDO and in our coronal red line images. The best-fit model consists of a flux rope with a twist angle of 3.1$\pi$, which is consistent with the most probable value of the total twist angle of interplanetary flux ropes observed at 1 AU. Linear polarization of the Fe XIII 10747 {\AA} line calculated from this model shows a "lagomorphic" signature that is also observed by the Coronal Multichannel Polarimeter of the High Altitude Observatory. We also find a ring-shaped structure in the line-of-sight velocity of Fe XIII 10747 {\AA}, which implies hot plasma flows along a helical magnetic field structure, in the cavity. These results suggest that the magnetic structure of the cavity is a highly twisted flux rope, which may erupt eventually. The temperature structure of the cavity has also been investigated using the intensity ratio of Fe XIII 10747 {\AA} and Fe X 6374 {\AA}.

We propose a new framework for testing gravity using cluster observations, which aims to provide an unbiased constraint on modified gravity models from Sunyaev Zel'dovich (SZ) and X-ray cluster counts and the cluster gas fraction, among other possible observables. Focusing on a popular $ f(R) $ model of gravity, we propose a novel procedure to recalibrate mass scaling relations from $ \Lambda $CDM to $ f(R) $ gravity for SZ and X-ray cluster observables. We find that the complicated modified gravity effects can be simply modelled as a dependence on a combination of the background scalar field and redshift, $ f_R(z)/(1+z) $, regardless of the $ f(R) $ model parameter. By employing a large suite of N-body simulations, we demonstrate that a theoretically derived tanh fitting formula is in excellent agreement with the dynamical mass enhancement of dark matter haloes for a large range of background field parameters and redshifts. Our framework is sufficiently flexible to allow for tests of other models and inclusion of further observables. The one-parameter description of the dynamical mass enhancement can have important implications on the theoretical modelling of observables and on practical tests of gravity.

Solid evidence of magnetic reconnection is rarely reported within sunspots, the darkest regions with the strongest magnetic fields and lowest temperatures in the solar atmosphere. Using the world's largest solar telescope, the 1.6-meter Goode Solar Telescope, we detect prevalent reconnection through frequently occurring fine-scale jets in the H${\alpha}$ line wings at light bridges, the bright lanes that may divide the dark sunspot core into multiple parts. Many jets have an inverted Y-shape, shown by models to be typical of reconnection in a unipolar field environment. Simultaneous spectral imaging data from the Interface Region Imaging Spectrograph show that the reconnection drives bidirectional flows up to 200~km~s$^{-1}$, and that the weakly ionized plasma is heated by at least an order of magnitude up to $\sim$80,000 K. Such highly dynamic reconnection jets and efficient heating should be properly accounted for in future modeling efforts of sunspots. Our observations also reveal that the surge-like activity previously reported above light bridges in some chromospheric passbands such as the H${\alpha}$ core has two components: the ever-present short surges likely to be related to the upward leakage of magnetoacoustic waves from the photosphere, and the occasionally occurring long and fast surges that are obviously caused by the intermittent reconnection jets.

About 766 $\delta$ Scuti stars were observed by LAMOST by June 16, 2017. Stellar atmospheric parameters of 525 variables were determined, while spectral types were obtained for all samples. In the paper, those spectroscopic data are catalogued. We detect a group of 131 unusual and cool variable stars (UCVs) that are distinguished from the normal $\delta$ Scuti stars (NDSTs). On the H-R diagram and the $\log g-T$ diagram, the UCVs are far beyond the red edge of pulsational instability trip. Their effective temperatures are lower than 6700\,K with periods in the range from 0.09 to 0.22\,days. NDSTs have metallicity close to that of the Sun as expected, while UCVs are slightly metal poor than NDSTs. The two peaks on the distributions of the period and stellar atmospheric parameters are all caused by the existence of UCVs. When those UCVs are excluded, it is discovered that the effective temperature, the gravitational acceleration and the metallicity all are correlated with the pulsating period for NDSTs and their physical properties and evolutionary states are discussed. The detecting of those UCVs indicates that about 25\% of the known $\delta$ Scuti stars may be misclassified. However, if some of them are confirmed to be pulsating stars, they will be a new-type pulsator and their investigations will shed light on theoretical instability domains and on the theories of interacting between the pulsation and the convection of solar-type stars. Meanwhile, 88 $\delta$ Scuti stars are detected to be the candidates of binary or multiple systems.

Analysis of all-sky Planck submillimetre observations and the IRAS 100um data has led to the detection of a population of Galactic cold clumps. The clumps can be used to study star formation and dust properties in a wide range of Galactic environments. Our aim is to measure dust spectral energy distribution (SED) variations as a function of the spatial scale and the wavelength. We examine the SEDs at large scales using IRAS, Planck, and Herschel data. At smaller scales, we compare with JCMT/SCUBA-2 850um maps with Herschel data that are filtered using the SCUBA-2 pipeline. Clumps are extracted using the Fellwalker method and their spectra are modelled as modified blackbody functions. According to IRAS and Planck data, most fields have dust colour temperatures T_C ~ 14-18K and opacity spectral index values of beta=1.5-1.9. The clumps/cores identified in SCUBA-2 maps have T~ 13K and similar beta values. There are some indications of the dust emission spectrum becoming flatter at wavelengths longer than 500um. In fits involving Planck data, the significance is limited by the uncertainty of the corrections for CO line contamination. The fits to the SPIRE data give a median beta value slightly above 1.8. In the joint SPIRE and SCUBA-2 850um fits the value decreases to beta ~1.6. Most of the observed T-beta anticorrelation can be explained by noise. The typical submillimetre opacity spectral index beta of cold clumps is found to be ~1.7. This is above the values of diffuse clouds but lower than in some previous studies of dense clumps. There is only tentative evidence of T-beta anticorrelation and beta decreasing at millimetre wavelengths.

About 3196 EA-type binaries (EAs) were observed by LAMOST by June 16, 2017 and their spectral types were derived. Meanwhile stellar atmospheric parameters of 2020 EAs were determined. In the paper, those EAs are catalogued and their physical properties and evolutionary states are investigated. The period distribution of EAs suggests that the period limit of tidal locking for the close binaries is about 6 days. It is found that the metallicity of EAs is higher than that of EWs indicating that EAs are generally younger than EWs and they are the progenitors of EWs. The metallicities of long-period EWs (0.4 < P < 1 days) are the same as those of EAs with the same periods, while their values of Log (g) are usually smaller than those of EAs. These support the evolutionary process that EAs evolve into long-period EWs through the combination of angular momentum loss (AML) via magnetic braking and case A mass transfer. For short-period EWs, their metallicities are lower than those of EAs, while their gravitational accelerations are higher. These reveal that they may be formed from cool short-period EAs through AML via magnetic braking with little mass transfer. For some EWs with high metallicities, they may be contaminated by material from the evolution of unseen neutron stars and black holes or they have third bodies that may help them to form rapidly through a short timescale of pre-contact evolution. The present investigation suggests that the modern EW populations may be formed through the combination of aforementioned mechanisms.

Hydrogen, being the most abundant element, is the driver of many if not most reactions occurring on interstellar dust grains. In hydrogen atom addition reactions, the rate is usually determined by the surface kinetics of the hydrogen atom instead of the other reaction partner. Three mechanisms exist to explain hydrogen addition reactions on surfaces: Langmuir-Hinshelwood, Eley-Rideal, and hot-atom. In gas-grain models, which mechanism is assumed greatly affects the simulation results. In this work, we quantify the temperature dependence of the rates of atomic hydrogen addition reactions by studying the reaction of H+O$_3$$\rightarrow$O$_2$+OH on the surface of a film of non-porous amorphous solid water (np-ASW) in the temperature range from 10 K to 50 K. The reaction rate is found to be temperature independent. This disagrees with the results of simulations with a network of rate equations that assume Langmuir-Hinshelwood mechanism through either thermal diffusion or tunneling diffusion; the reaction rates assuming such mechanism possesses a strong temperature dependence, either explicitly or implicitly, that is not seen experimentally. We suggest that the Eley-Rideal and/or hot-atom mechanism play a key role in hydrogen atom addition reactions, and should be included in gas-grain models. We also suggest that our newly developed time-resolved reactive scattering can be utilized to measure the chemical desorption efficiency in grain surface reactions.

The thermonuclear $^{19}$F($p$,$\alpha_0$)$^{16}$O reaction rate in a temperature region of 0.007--10 GK has been derived by re-evaluating the available experimental data, together with the low-energy theoretical $R$-matrix extrapolations. Our new rate deviates up to about 30\% compared to the previous ones, although all rates are consistent within the uncertainties. At very low temperature (e.g. 0.01 GK) our reaction rate is about 20\% smaller than the most recently published rate, because of a difference in the low energy extrapolated $S$-factor and a more accurate estimate of the reduced mass entering in the calculation of the reaction rate. At temperatures above $\sim$1 GK, our rate is smaller, for instance, by about 20\% around 1.75 GK, because we have re-evaluated in a meticulous way the previous data (Isoya et al., Nucl. Phys. 7, 116 (1958)). The present interpretation is supported by the direct experimental data. The uncertainties of the present evaluated rate are estimated to be about 20\% in the temperature region below 0.2 GK, which are mainly caused by the lack of low-energy experimental data and the large uncertainties of the existing data. The asymptotic giant branch (AGB) star evolves at temperatures below 0.2 GK, where the $^{19}$F($p$,$\alpha$)$^{16}$O reaction may play a very important role. However, the current accuracy of the reaction rate is insufficient to help to describe, in a careful way, for the fluorine overabundances phenomenon observed in AGB stars. Precise cross section (or $S$ factor) data in the low energy region are therefore mandatory for astrophysical nucleosynthesis studies.

Carbon dioxide is abundant in ice mantles of dust grains; some is found in the pure crystalline form as inferred from the double peak splitting of the bending profile at about 650 cm$^{-1}$. To study how CO$_2$ segregates into the pure form from water-rich mixtures of ice mantles and how it then crystallizes, we used Reflection Absorption InfraRed Spectroscopy (RAIRS) to study the structural change of pure CO$_2$ ice as a function of both ice thickness and temperature. We found that the $\nu_1+\nu_3$ combination mode absorption profile at 3708~cm$^{-1}$ provides an excellent probe to quantify the degree of crystallinity in CO$_2$ ice. We also found that between 20 and 30~K, there is an ordering transition that we attribute to reorientation of CO$_2$ molecules, while the diffusion of CO$_2$ becomes significant at much higher temperatures. In the formation of pure crystalline CO$_2$ in ISM ices, the rate limiting process is the diffusion/segregation of CO$_2$ molecules in the ice instead of the phase transition from amorphous to crystalline after clusters/islands of CO$_2$ are formed.

A single dish monitoring of millimeter maser lines SiS J=14-13 and HCN nu_2 = 1^f J=3-2 and several other rotational lines is reported for the archetypal carbon star IRC+10216. Relative line strength variations of 5%~30% are found for eight molecular line features with respect to selected reference lines. Definite line-shape variation is found in limited velocity intervals of the SiS and HCN line profiles. The asymmetrical line profiles of the two lines are mainly due to the varying components. Their dominant varying components of the line profiles have similar periods and phases as the IR light variation, although both quantities show some degree of velocity dependence; there is also variability asymmetry between the blue and red line wings of both lines. Combining the velocities and amplitudes with a wind velocity model, we suggest that the line profile variations are due to SiS and HCN masing lines emanating from the wind acceleration zone. The possible link of the variabilities to thermal, dynamical and/or chemical processes within or under this region is also discussed.

Numerous EWs were discovered by several deep photometric survey and there are about 40785 EW-type binary systems listed in the international variable star index (VSX) by March 13, 2017. 7938 of them were observed by LAMOST by November 30, 2016 and their spectral types were given. Stellar atmospheric parameters of 5363 EW-type binary stars were determined based on good spectroscopic observations. In the paper, those EWs were catalogued and their properties are analyzed. The distributions of the orbital period (P), the effect temperature (T), the gravitational acceleration (Log(g)), the metallicity ([Fe/H]) and the radial velocity (RV) are presented for those observed EW-type systems. It is shown that about 80.6\% sample stars have metallicity below zero indicating that EW-type systems are old stellar population. This is in agreement with the the conclusion that the EW binaries are formed from moderately close binaries through angular momentum loss via magnetic braking that takes a few hundred million to a few billion years. The unusual high metallicities of a few percent of EWs may be caused by contaminating of material from the evolution of unseen neutron stars and black holes in the systems. The correlations between the orbital period and the effect temperature, the gravitational acceleration and the metallicity are presented and their scatters are mainly caused by (i) the presence of the third bodies and (ii) the wrong determined periods sometimes. It is shown that some EW contain evolved component stars and the physical properties of EWs are mainly depending on their orbital periods. It is found that the extremely short-period EWs may be older than their long-period cousins because they have lower metallicities. This reveals that they have a longer timescale of pre-contact evolution and their formation and evolution are mainly driven by angular momentum loss via magnetic braking.

We extend previous studies of big bang nucleosynthesis, with the assumption that ordinary matter and dark matter sectors are entangled through the number of degrees of freedom entering the Friedmann equations. This conjecture allows us to find a relation between the temperatures of the weakly interacting matter and dark-matter sectors. The constraints imposed by observations are studied by comparison with calculations of big bang nucleosynthesis for the abundance of light elements. It is found that these constraints are compatible with cold dark matter and a wide range number of dark sectors.

Within one billion years of the Big Bang, intergalactic hydrogen was ionized by sources emitting ultraviolet and higher energy photons. This was the final phenomenon to globally affect all the baryons (visible matter) in the Universe. It is referred to as cosmic reionization and is an integral component of cosmology. It is broadly expected that intrinsically faint galaxies were the primary ionizing sources due to their abundance in this epoch. However, at the highest redshifts ($z>7.5$; lookback time 13.1 Gyr), all galaxies with spectroscopic confirmations to date are intrinsically bright and, therefore, not necessarily representative of the general population. Here, we report the unequivocal spectroscopic detection of a low luminosity galaxy at $z>7.5$. We detected the Lyman-$\alpha$ emission line at $\sim 10504$ {\AA} in two separate observations with MOSFIRE on the Keck I Telescope and independently with the Hubble Space Telescope's slit-less grism spectrograph, implying a source redshift of $z = 7.640 \pm 0.001$. The galaxy is gravitationally magnified by the massive galaxy cluster MACS J1423.8+2404 ($z = 0.545$), with an estimated intrinsic luminosity of $M_{AB} = -19.6 \pm 0.2$ mag and a stellar mass of $M_{\star} = 3.0^{+1.5}_{-0.8} \times 10^8$ solar masses. Both are an order of magnitude lower than the four other Lyman-$\alpha$ emitters currently known at $z > 7.5$, making it probably the most distant representative source of reionization found to date.

As a space-borne detector POLAR is designed to conduct hard X-ray polarization measurements of gamma-ray bursts on the statistically significant sample of events and with an unprecedented accuracy. During its development phase a number of tests, calibrations runs and verification measurements were carried out in order to validate instrument functionality and optimize operational parameters. In this article we present results on gain optimization togeter with verification data obtained in the course of broad laboratory and environmental tests. In particular we focus on exposures to the $^{137}$Cs radioactive source and determination of the gain dependence on the high voltage for all 1600 detection channels of the polarimeter. Performance of the instrument is described in detail with respect to the dynamic range, energy resolution and temperature dependence. Gain optimization algorithms and response non-uniformity studies are also broadly discussed. Results presented below constitute important parts for development of the POLAR calibration and operation database.

We describe and demonstrate the potential of a new and very efficient method for simulating certain classes of modified gravity theories, such as the widely studied $f(R)$ gravity models. High resolution simulations for such models are currently very slow due to the highly nonlinear partial differential equation that needs to be solved exactly to predict the modified gravitational force. This nonlinearity is partly inherent, but is also exacerbated by the specific numerical algorithm used, which employs a variable redefinition to prevent numerical instabilities. The standard Newton-Gauss-Seidel iterative method used to tackle this problem has a poor convergence rate. Our new method not only avoids this, but also allows the discretised equation to be written in a form that is analytically solvable. We show that this new method greatly improves the performance and efficiency of $f(R)$ simulations. For example, a test simulation with $512^3$ particles in a box of size $512 \, \mathrm{Mpc}/h$ is now 5 times faster than before, while a Millennium-resolution simulation for $f(R)$ gravity is estimated to be more than 20 times faster than with the old method. Our new implementation will be particularly useful for running very high resolution, large-sized simulations which, to date, are only possible for the standard model, and also makes it feasible to run large numbers of lower resolution simulations for covariance analyses. We hope that the method will bring us to a new era for precision cosmological tests of gravity.

High-resolution observations of the solar chromosphere and transition region often reveal surge-like oscillatory activities above sunspot light bridges. These oscillations are often interpreted as intermittent plasma jets produced by quasi-periodic magnetic reconnection. We have analyzed the oscillations above a light bridge in a sunspot using data taken by the Interface Region Imaging Spectrograph (IRIS). The chromospheric 2796\AA{}~images show surge-like activities above the entire light bridge at any time, forming an oscillating wall. Within the wall we often see that the Mg~{\sc{ii}}~k 2796.35\AA{}~line core first experiences a large blueshift, and then gradually decreases to zero shift before increasing to a red shift of comparable magnitude. Such a behavior suggests that the oscillations are highly nonlinear and likely related to shocks. In the 1400\AA{}~passband which samples emission mainly from the Si~{\sc{iv}}~ion, the most prominent feature is a bright oscillatory front ahead of the surges. We find a positive correlation between the acceleration and maximum velocity of the moving front, which is consistent with numerical simulations of upward propagating slow-mode shock waves. The Si~{\sc{iv}} 1402.77\AA{}~line profile is generally enhanced and broadened in the bright front, which might be caused by turbulence generated through compression or by the shocks. These results, together with the fact that the oscillation period stays almost unchanged over a long duration, lead us to propose that the surge-like oscillations above light bridges are caused by shocked p-mode waves leaked from the underlying photosphere.

Solar filaments/prominences are one of the most common features in the corona, which may lead to energetic coronal mass ejections (CMEs) and flares when they erupt. Filaments are about one hundred times cooler and denser than the coronal material, and physical understanding of their material origin remains controversial. Two types of scenarios have been proposed: one argues that the filament plasma is brought into the corona from photosphere or chromosphere through a siphon or evaporation/injection process, while the other suggests that the material condenses from the surrounding coronal plasma due to thermal instability. The elemental abundance analysis is a reasonable clue to constrain the models, as the siphon or evaporation/injection model would predict that the filament material abundances are close to the photospheric or chromospheric ones, while the condensation model should have coronal abundances. In this letter, we analyze the elemental abundances of a magnetic cloud that contains the ejected filament material. The corresponding filament eruption occurred on 1998 April 29, accompanying an M6.8 class soft X-ray flare located at the heliographic coordinates S18E20 (NOAA 08210) and a fast halo CME with the linear velocity of 1374 km s$^{-1}$ near the Sun. We find that the abundance ratios of elements with low and high First Ionization Potential such as Fe/O, Mg/O, and Si/O are 0.150, 0.050, and 0.070, respectively, approaching their corresponding photospheric values 0.065, 0.081, and 0.066, which does not support the coronal origin of the filament plasma.

Big Bang nucleosynthesis (BBN) theory predicts the abundances of the light elements D, $^3$He, $^4$He and $^7$Li produced in the early universe. The primordial abundances of D and $^4$He inferred from observational data are in good agreement with predictions, however, the BBN theory overestimates the primordial $^7$Li abundance by about a factor of three. This is the so-called "cosmological lithium problem". Solutions to this problem using conventional astrophysics and nuclear physics have not been successful over the past few decades, probably indicating the presence of new physics during the era of BBN. We have investigated the impact on BBN predictions of adopting a generalized distribution to describe the velocities of nucleons in the framework of Tsallis non-extensive statistics. This generalized velocity distribution is characterized by a parameter $q$, and reduces to the usually assumed Maxwell-Boltzmann distribution for $q$ = 1. We find excellent agreement between predicted and observed primordial abundances of D, $^4$He and $^7$Li for $1.069\leq q \leq 1.082$, suggesting a possible new solution to the cosmological lithium problem.

We report the observations of an electron vortex magnetic hole corresponding to a new type of coherent structures in the magnetosheath turbulent plasma using the Magnetospheric Multiscale (MMS) mission data. The magnetic hole is characterized by a magnetic depression, a density peak, a total electron temperature increase (with a parallel temperature decrease but a perpendicular temperature increase), and strong currents carried by the electrons. The current has a dip in the center of the magnetic hole and a peak in the outer region of the magnetic hole. The estimated size of the magnetic hole is about 0.23 \r{ho}i (~ 30 \r{ho}e) in the circular cross-section perpendicular to its axis, where \r{ho}i and \r{ho}e are respectively the proton and electron gyroradius. There are no clear enhancement seen in high energy electron fluxes, but an enhancement in the perpendicular electron fluxes at ~ 90{\deg} pitch angles inside the magnetic hole is seen, implying that the electron are trapped within it. The variations of the electron velocity components Vem and Ven suggest that an electron vortex is formed by trapping electrons inside the magnetic hole in the circular cross-section (in the M-N plane). These observations demonstrate the existence of a new type of coherent structures behaving as an electron vortex magnetic hole in turbulent space plasmas as predicted by recent kinetic simulations.

Gamma-ray polarimetry is a new powerful tool to study the processes responsible for the emission from astrophysical sources and the environments in which this emission takes place. Few successful polarimetric measurements have however been performed thus far in the gamma-ray energy band due to the difficulties involved. POLAR is a dedicated polarimeter designed to perform high precision measurements of the polarization of the emission from gamma-ray burst in the 50-500 keV energy range. This new polarimeter is expected to detect approximately 50 gamma-ray bursts per year while performing high precision polarization measurements on approximately 10 bursts per year. The instrument was launched into lower earth orbit as part of the second Chinese space lab, the Tiangong-2, on September 15th 2016 and has been taking data successfully since being switched on one week after. The instrument uses a segmented scintillator array consisting of 1600 plastic scintillator bars, read out by 25 flat-panel multi-anode photomultipliers, to measure the Compton scattering angles of incoming photons. The small segmentation and relatively large uniform effective area allow the instrument to measure the polarization of a large number of transient events, such as gamma-ray bursts, with an unprecedented precision during its two year life-time. The final flight model underwent detailed calibration prior to launch as well as intensive space qualification tests, a summary of which will be presented in this paper. The instrument design will be discussed first followed by an overview of the on-ground tests, finally the in-orbit behavior as measured during the first weeks of the mission will be presented.

We observed thirteen Planck cold clumps with the James Clerk Maxwell Telescope/SCUBA-2 and with the Nobeyama 45 m radio telescope. The N$_2$H$^+$ distribution obtained with the Nobeyama telescope is quite similar to SCUBA-2 dust distribution. The 82 GHz HC$_3$N, 82 GHz CCS, and 94 GHz CCS emission are often distributed differently with respect to the N$_2$H$^+$ emission. The CCS emission, which is known to be abundant in starless molecular cloud cores, is often very clumpy in the observed targets. We made deep single-pointing observations in DNC, HN$^{13}$C, N$_2$D$^+$, cyclic-C$_3$H$_2$ toward nine clumps. The detection rate of N$_2$D$^+$ is 50\%. Furthermore, we observed the NH$_3$ emission toward 15 Planck cold clumps to estimate the kinetic temperature, and confirmed that most of targets are cold ($\lesssim$ 20 K). In two of the starless clumps observe, the CCS emission is distributed as it surrounds the N$_2$H$^+$ core (chemically evolved gas), which resembles the case of L1544, a prestellar core showing collapse. In addition, we detected both DNC and N$_2$D$^+$. These two clumps are most likely on the verge of star formation. We introduce the Chemical Evolution Factor (CEF) for starless cores to describe the chemical evolutionary stage, and analyze the observed Planck cold clumps.

Observations of solar wind turbulence indicate the existence of multi-scale pressure-balanced structures (PBSs) in the solar wind. In this work, we conduct a numerical simulation to investigate multi-scale PBSs and in particular their formation in compressive MHD turbulence. By the use of a higher order Godunov code Athena,a driven compressible turbulence with an imposed uniform guide field is simulated. The simulation results show that both the magnetic pressure and the thermal pressure exhibit a turbulent spectrum with a Kolmogorov-like power law, and that in many regions of the simulation domain they are anti-correlated. The computed wavelet cross-coherence spectrum of the magnetic pressure and the thermal pressure, as well as their space series, indicate the existence of multi-scale PBSs, with the small PBSs being embedded in the large ones. These multi-scale PBSs are likely to be related with the highly oblique-propagating slow-mode waves, as the traced multi-scale PBS is found to be traveling in a certain direction at a speed consistent with that predicted theoretically for a slow-mode wave propagating in the same direction.

Nuclear fusion reactions are the most important processes in nature to power stars and produce new elements, and lie at the center of the understanding of nucleosynthesis in the universe. It is critically important to study the reactions in full plasma environments that are close to true astrophysical conditions. By using laser-driven counter-streaming collisionless plasmas, we studied the fusion D$+$D$\rightarrow n +^3$He in a Gamow-like window around 27 keV. The results show that astrophysical nuclear reaction yield can be modulated significantly by the self-generated electromagnetic fields and the collective motion of the plasma. This plasma-version mini-collider may provide a novel tool for studies of astrophysics-interested nuclear reactions in plasma with tunable energies in earth-based laboratories.

Observations by ISO and Spitzer towards young stellar objects (YSOs) showed that CO$_2$ segregates in the icy mantles covering dust grains. Thermal processing of ice mixture was proposed as responsible for the segregation. Although several laboratory studied thermally induced segregation, a satisfying quantification is still missing. We propose that the diffusion of CO$_2$ along pores inside water ice is the key to quantify segregation. We combined Temperature Programmed Desorption (TPD) and Reflection Absorption InfraRed Spectroscopy (RAIRS) to study how CO$_2$ molecules interact on a non-porous amorphous solid water (np-ASW) surface. We found that CO$_2$ diffuses significantly on a np-ASW surface above 65~K and clusters are formed at well below one monolayer. A simple rate equation simulation finds that the diffusion energy barrier of CO$_2$ on np-ASW is 2150$\pm$50 K, assuming a diffusion pre-exponential factor of 10$^{12}$ s$^{-1}$. This energy should also apply to the diffusion of CO$_2$ on wall of pores. The binding energy of CO$_2$ from CO$_2$ clusters and CO$_2$ from H$_2$O ice have been found to be $2415\pm20$ and $2250\pm20$~K, respectively, assuming the same prefactor for desorption. CO$_2$-CO$_2$ interaction is stronger than CO$_2$-H$_2$O interaction, in agreement with the experimental finding that CO$_2$ does not wet np-ASW surface. For comparison, we carried out similar experiments with CO on np-ASW, and found that the CO-CO interaction is always weaker than CO-H$_2$O. As a result, CO wets np-ASW surface. This study should be of help to uncover the thermal history of CO$_2$ on the icy mantles of dust grains.

Using the liminality N-body simulations of Shi et. al., we present the first predictions for galaxy clustering in $f(R)$ gravity using subhalo abundance matching. We find that, for a given galaxy density, even for an $f(R)$ model with $f_{R0}=-10^{-6}$, for which the cold dark matter clustering is very similar to $\Lambda$CDM, the predicted clustering of galaxies in the $f(R)$ model is very different from $\Lambda$CDM. The deviation can be as large as $40\%$ for samples with mean densities close to that of $L_*$ galaxies. This large deviation is testable given the accuracy that future large-scale galaxy surveys aim to achieve. Our result demonstrates that galaxy surveys can provide a stringent test of General Relativity on cosmological scales, which is comparable to the tests from local astrophysical observations.

Four candidates of eclipsing multiples, based on new extraneous eclipses found on Kepler binary light curves, are presented and studied. KIC 7622486 is a double eclipsing binary candidate with orbital period of 2.2799960 days and 40.246503 days. The two binary systems do not eclipse each other in the line of sight, but there is mutual gravitational influence between them which leads to the small but definite eccentricity 0.0035(0.0022) on the short 2.2799960 days period orbit. KIC 7668648 is hierarchical quadruple system candidate, with two sets of solid 203(+-5) days period extraneous eclipses and another independent set of extraneous eclipses. A clear and credible extraneous eclipse is found on the binary light curve of KIC 7670485 which made it a triple system candidates. Two sets of extraneous eclipse of about 390 days and 220 days period are found on KIC 8938628 binary curves, which not only confirms the previous conclusion of $388.5(+-0.3) triple system, but also proposed a new additional objects that make KIC 8938628 a hierarchical quadruple system candidates. It is wished that the four candidates here will make some contributions to the field of eclipsing multiples.

The thermonuclear rate of the 50Fe(p,gamma)51Co reaction in the Type I X-ray bursts (XRBs) temperature range has been reevaluated based on a recent precise mass measurement at CSRe lanzhou, where the proton separation energy Sp=142+/-77 keV has been determined firstly for the 51Co nucleus. Comparing to the previous theoretical predictions, the experimental Sp value has much smaller uncertainty. Based on the nuclear shell model and mirror nuclear structure information, we have calculated two sets of thermonuclear rates for the 50Fe(p,gamma)51Co reaction by utilizing the experimental Sp value. It shows that the statistical-model calculations are not ideally applicable for this reaction primarily because of the low density of low-lying excited states in 51Co. In this work, we recommend that a set of new reaction rate based on the mirror structure of 51Cr should be incorporated in the future astrophysical network calculations.

Modified gravity theories are a popular alternative to dark energy as a possible explanation for the observed accelerating cosmic expansion, and their cosmological tests are currently an active research field. Studies in recent years have been increasingly focused on testing these theories in the nonlinear regime, which is computationally demanding. Here we show that, under certain circumstances, a whole class of theories can be ruled out by using background cosmology alone. This is possible because certain classes of models (i) are fundamentally incapable of producing specific background expansion histories, and (ii) said histories are incompatible with local gravity tests. As an example, we demonstrate that a popular class of models, $f(R)$ gravity, would not be viable if observations suggest even a slight deviation of the background expansion history from that of the $\Lambda$CDM paradigm.

In the solar atmosphere, jets are prevalent and they are significant for the mass and energy transport. Here we conduct numerical simulations to investigate the mass and energy contributions of the recently observed high-speed jets to the solar wind. With a one-dimensional hydrodynamic solar wind model, the time-dependent pulses are imposed at the bottom to simulate the jets. The simulation results show that without other energy source, the injected plasmas are accelerated effectively to be a transonic wind with a substantial mass flux. The rapid acceleration occurs close to the Sun, and the resulting asymptotic speed, number density at 0.3 AU, as well as mass flux normalized to 1 AU are compatible with in situ observations. As a result of the high speed, the imposed pulses generate a train of shocks traveling upward. By tracing the motions of the injected plasma, it is found that these shocks heat and accelerate the injected plasmas successively step by step to push them upward and eventually allow them to escape. The parametric studies show that increasing the speed of the imposed pulses or their temperature gives a considerably faster, and hotter solar wind, while increasing their number density or decreasing their recurring period only bring a denser solar wind. These studies provide a possibility that the ubiquitous high-speed jets are a substantial mass and energy contributions to the solar wind.

In the solar atmosphere, the jets are ubiquitous and found to be at various spatia-temporal scales. They are significant to understand energy and mass transport in the solar atmosphere. Recently, the high-speed transition region jets are reported from the observation. Here we conduct a numerical simulation to investigate the mechanism in their formation. Driven by the supergranular convection motion, the magnetic reconnection between the magnetic loop and the background open flux occurring in the transition region is simulated with a two-dimensional magnetohydrodynamics model. The simulation results show that not only a fast hot jet, much resemble the found transition region jets, but also a adjacent slow cool jet, mostly like classical spicules, is launched. The force analysis shows that the fast hot jet is continually driven by the Lorentz force around the reconnection region, while the slow cool jet is induced by an initial kick through the Lorentz force associated with the emerging magnetic flux. Also, the features of the driven jets change with the amount of the emerging magnetic flux, giving the varieties of both jets. These results will inspire our understanding of the formation of the prevalence of both the fast hot jet and slow cool jet from the solar transition region and chromosphere.

The $B$ $V$ $R_c$ $I_c$ bands light curves of the newly discovered binary system \astrobj{GSC 03122-02426} are obtained and analyzed using the Wilson-Devinney (W-D) code. The solutions suggest that the mass ratio of the binary system is $q = 2.70$ and the less massive component is $422K$ hotter than the more massive one. We conclude that \astrobj{GSC 03122-02426} is a W-subtype shallow contact (with a contact degree of $f = 15.3\,\%$) binary system. It may be a newly formed contact binary system which is just under geometrical contact and will evolve to be a thermal contact binary system. The high orbital inclination ($i = 81.6^{\circ}$) implies that \astrobj{GSC 03122-02426} is a total eclipsing binary system and the photometric parameters obtained by us are quite reliable. We also estimate the absolute physical parameters of the two components in \astrobj{GSC 03122-02426}, which will provide fundamental information for the research of contact binary systems. The formation and evolutionary scenario of \astrobj{GSC 03122-02426} is discussed.

Two sets of light curves in $V$ $R_c$ $I_c$ bands for a newly discovered binary system UCAC4 436-062932 are obtained and analyzed using the Wilson-Devinney (W-D) code. The two sets of light curves get almost consistent results. The determined mass ratio is about $q = 2.7$ and the less massive component is nearly $250K$ hotter than the more massive one. The solutions conclude that UCAC4 436-062932 is a W-subtype shallow contact (with a contact degree of $f = 20\,\%$) binary system. Since the O'Connell effect appears on one set of the light curves, theories proposed to explain the effect are discussed. We assume that spot model may be the more plausible one to the O'Connell effect appeared on the asymmetric light curves of the binary system UCAC4 436-062932. Therefore, we add a cool spot on the surface of the more massive star (component with lower effective temperature) and get a quite approving results for the light curve fitting. It will provide evidence to support the spot model in the explanatory mechanism of O'Connell effect.

In this work we investigate the effects of ion accretion and size-dependent dust temperatures on the abundances of both gas-phase and grain-surface species. While past work has assumed a constant areal density for icy species, we show that this assumption is invalid and the chemical differentiation over grain sizes are significant. We use a gas-grain chemical code to numerically demonstrate this in two typical interstellar conditions: dark cloud (DC) and cold neutral medium (CNM). It is shown that, although the grain size distribution variation (but with the total grain surface area unchanged) has little effect on the gas-phase abundances, it can alter the abundances of some surface species by factors up to $\sim2-4$ orders of magnitude. The areal densities of ice species are larger on smaller grains in the DC model as the consequence of ion accretion. However, the surface areal density evolution tracks are more complex in the CNM model due to the combined effects of ion accretion and dust temperature variation. The surface areal density differences between the smallest ($\sim 0.01\mu$m) and the biggest ($\sim 0.2\mu$m) grains can reach $\sim$1 and $\sim$5 orders of magnitude in the DC and CNM models, respectively.

New high-resolution far-infrared (FIR) observations of both ortho- and para-NH3 transitions toward IRC+10216 were obtained with Herschel, with the goal of determining the ammonia abundance and constraining the distribution of NH3 in the envelope of IRC+10216. We used the Heterodyne Instrument for the Far Infrared (HIFI) on board Herschel to observe all rotational transitions up to the J=3 level (three ortho- and six para-NH3 lines). We conducted non-LTE multilevel radiative transfer modelling, including the effects of near-infrared (NIR) radiative pumping through vibrational transitions. We found that NIR pumping is of key importance for understanding the excitation of rotational levels of NH3. The derived NH3 abundances relative to molecular hydrogen were (2.8+-0.5)x10^{-8} for ortho-NH3 and (3.2^{+0.7}_{-0.6})x10^{-8} for para-NH3, consistent with an ortho/para ratio of 1. These values are in a rough agreement with abundances derived from the inversion transitions, as well as with the total abundance of NH3 inferred from the MIR absorption lines. To explain the observed rotational transitions, ammonia must be formed near to the central star at a radius close to the end of the wind acceleration region, but no larger than about 20 stellar radii (1 sigma confidence level).

We propose a new method to model cluster scaling relations in modified gravity. Using a suite of non-radiative hydrodynamical simulations, we show that the scaling relations of cumulative gas quantities, such as the Sunyaev Zel'dovich effect (Compton-y parameter) and the x-ray Compton-Y parameter, can be accurately predicted using the known results in the $\Lambda$CDM model with a precision of $\sim3\%$. This method provides a reliable way to analyze the gas physics in modified gravity using the less-demanding and much more efficient pure cold dark matter simulations. Our results therefore have important theoretical and practical implications in constraining gravity using cluster surveys.

An observation from the Interface Region Imaging Spectrograph reveals coherent oscillations in the loops of an M1.6 flare on 2015 March 12. Both the intensity and Doppler shift of Fe~{\sc{xxi}}~1354.08\AA{}~show clear oscillations with a period of $\sim$25 seconds. Remarkably similar oscillations were also detected in the soft X-ray flux recorded by the Geostationary Operational Environmental Satellites (GOES). With an estimated phase speed of $\sim$2420~km~s$^{-1}$ and a derived electron density of at least 5.4$\times$10$^{10}$ cm$^{-3}$, the observed short-period oscillation is most likely the global fast sausage mode of a hot flare loop. We find a phase shift of $\sim$$\pi$/2 (1/4 period) between the Doppler shift oscillation and the intensity/GOES oscillations, which is consistent with a recent forward modeling study of the sausage mode. The observed oscillation requires a density contrast between the flare loop and coronal background of a factor $\geqslant$42. The estimated phase speed of the global mode provides an lower limit of the Alfv\'en speed outside the flare loop. We also find an increase of the oscillation period, which might be caused by the separation of the loop footpoints with time.

Accurate modeling of physical and chemical processes in the interstellar medium requires detailed knowledge of how atoms and molecule adsorb on dust grains. However, the sticking coefficient, a number between 0 and 1 that measures the first step in the interaction of a particle with a surface, is usually assumed in simulations of ISM environments to be either 0.5 or 1. Here we report on the determination of the sticking coefficient of H$_2$, D$_2$, N$_2$, O$_2$, CO, CH$_4$, and CO$_2$ on non-porous amorphous solid water (np-ASW). The sticking coefficient was measured over a wide range of surface temperatures using a highly collimated molecular beam. We showed that the standard way of measuring the sticking coefficient --- the King-Wells method --- leads to the underestimation of trapping events in which there is incomplete energy accommodation of the molecule on the surface. Surface scattering experiments with the use of a pulsed molecular beam are used instead to measure the sticking coefficient. Based on the values of the measured sticking coefficient we suggest a useful general formula of the sticking coefficient as a function of grain temperature and molecule-surface binding energy. We use this formula in a simulation of ISM gas-grain chemistry to find the effect of sticking on the abundance of key molecules both on grains and in the gas-phase.

Recent observations by the Interface Region Imaging Spectrograph (IRIS) have revealed pockets of hot gas ($\sim$2--8$\times$10$^{4}$ K) potentially resulting from magnetic reconnection in the partially ionized lower solar atmosphere (IRIS bombs; IBs). Using joint observations between IRIS and the Chinese New Vacuum Solar Telescope, we have identified ten IBs. We find that three are unambiguously and three others are possibly connected to Ellerman bombs (EBs), which show intense brightening of the extended H$_{\alpha}$ wings without leaving an obvious signature in the H$_{\alpha}$ core. These bombs generally reveal the following distinct properties: (1) The O~{\sc{iv}}~1401.156\AA{} and 1399.774\AA{} lines are absent or very weak; (2) The Mn~{\sc{i}}~2795.640\AA{} line manifests as an absorption feature superimposed on the greatly enhanced Mg~{\sc{ii}}~k line wing; (3) The Mg~{\sc{ii}}~k and h lines show intense brightening in the wings and no dramatic enhancement in the cores; (4) Chromospheric absorption lines such as Ni~{\sc{ii}}~1393.330\AA{} and 1335.203\AA{} are very strong; (5) The 1700\AA{} images obtained with the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory reveal intense and compact brightenings. These properties support the formation of these bombs in the photosphere, demonstrating that EBs can be heated much more efficiently than previously thought. We also demonstrate that the Mg~{\sc{ii}}~k and h lines can be used to investigate EBs similarly to H$_{\alpha}$, which opens a promising new window for EB studies. The remaining four IBs obviously have no connection to EBs and they do not have the properties mentioned above, suggesting a higher formation layer possibly in the chromosphere.

We measured the binding energy of N$_2$, CO, O$_2$, CH$_4$, and CO$_2$ on non-porous (compact) amorphous solid water (np-ASW), of N$_2$ and CO on porous amorphous solid water (p-ASW), and of NH$_3$ on crystalline water ice. We were able to measure binding energies down to a fraction of 1\% of a layer, thus making these measurements more appropriate for astrochemistry than the existing values. We found that CO$_2$ forms clusters on np-ASW surface even at very low coverages. The binding energies of N$_2$, CO, O$_2$, and CH$_4$ decrease with coverage in the submonolayer regime. Their values at the low coverage limit are much higher than what is commonly used in gas-grain models. An empirical formula was used to describe the coverage dependence of the binding energies. We used the newly determined binding energy distributions in a simulation of gas-grain chemistry for cold cloud and hot core models. We found that owing to the higher value of desorption energy in the sub-monlayer regime a fraction of all these ices stays much longer and up to higher temperature on the grain surface compared to the single value energies currently used in the astrochemical models.

The extent of nucleosynthesis in models of type I X-ray bursts and the associated impact on the energy released in these explosive events are sensitive to nuclear masses and reaction rates around the $^{64}$Ge waiting point. Using the well known mass of $^{64}$Ge, the recently measured $^{65}$As mass, and large-scale shell model calculations, we have determined new thermonuclear rates of the $^{64}$Ge($p$,$\gamma$)$^{65}$As and $^{65}$As($p$,$\gamma$)$^{66}$Se reactions with reliable uncertainties. The new reaction rates differ significantly from previously published rates. Using the new data we analyze the impact of the new rates and the remaining nuclear physics uncertainties on the $^{64}$Ge waiting point in a number of representative one-zone X-ray burst models. We find that in contrast to previous work, when all relevant uncertainties are considered, a strong $^{64}$Ge $rp$-process waiting point cannot be ruled out. The nuclear physics uncertainties strongly affect X-ray burst model predictions of the synthesis of $^{64}$Zn, the synthesis of nuclei beyond $A=64$, energy generation, and burst light curve. We also identify key nuclear uncertainties that need to be addressed to determine the role of the $^{64}$Ge waiting point in X-ray bursts. These include the remaining uncertainty in the $^{65}$As mass, the uncertainty of the $^{66}$Se mass, and the remaining uncertainty in the $^{65}$As($p$,$\gamma$)$^{66}$Se reaction rate, which mainly originates from uncertain resonance energies.

Quasi-periodic disturbances of emission-line parameters are frequently observed in the corona. These disturbances propagate upward along the magnetic field with speeds $\sim100~\rm{km~s}^{-1}$. This phenomenon has been interpreted as evidence of the propagation of slow magnetosonic waves or argued to be signature of the intermittent outflows superposed on the background plasmas. Here we aim to present a new "wave + flow" model to interpret these observations. In our scenario, the oscillatory motion is a slow mode wave, and the flow is associated with a beam created by the wave-particle interaction owing to Landau resonance. With the help of a Vlasov model, we simulate the propagation of the slow mode wave and the generation of the beam flow. We find that weak periodic beam flows can be generated owing to Landau resonance in the solar corona, and the phase with strongest blueward asymmetry is ahead of that with strongest blueshift by about 1/4 period. We also find that the slow wave damps to the level of 1/e after the transit time of two wave periods, owing to Landau damping and Coulomb collisions in our simulation. This damping time scale is similar to that resulting from thermal-conduction in the magnetohydrodynamics regime. The beam flow is weakened/attenuated with increasing wave period and decreasing wave amplitude since Coulomb collision becomes more and more dominant over the wave action. We suggest that this "wave + flow" kinetic model provides an alternative explanation for the observed quasi-periodic propagating perturbations in various parameters in the solar corona.

We are performing a series of observations with ground-based telescopes toward Planck Galactic cold clumps (PGCCs) in the $\lambda$ Orionis complex in order to systematically investigate the effects of stellar feedback. In the particular case of PGCC G192.32-11.88, we discovered an extremely young Class 0 protostellar object (G192N) and a proto-brown dwarf candidate (G192S). G192N and G192S are located in a gravitationally bound bright-rimmed clump. The velocity and temperature gradients seen in line emission of CO isotopologues indicate that PGCC G192.32-11.88 is externally heated and compressed. G192N probably has the lowest bolometric luminosity ($\sim0.8$ L$_{\sun}$) and accretion rate (6.3$\times10^{-7}$ M$_{\sun}$~yr$^{-1}$) when compared with other young Class 0 sources (e.g. PACS Bright Red sources (PBRs)) in the Orion complex. It has slightly larger internal luminosity ($0.21\pm0.01$ L$_{\sun}$) and outflow velocity ($\sim$14 km~s$^{-1}$) than the predictions of first hydrostatic cores (FHSCs). G192N might be among the youngest Class 0 sources, which are slightly more evolved than a FHSC. Considering its low internal luminosity ($0.08\pm0.01$ L$_{\odot}$) and accretion rate (2.8$\times10^{-8}$ M$_{\sun}$~yr$^{-1}$), G192S is an ideal proto-brown dwarf candidate. The star formation efficiency ($\sim$0.3\%-0.4\%) and core formation efficiency ($\sim$1\%) in PGCC G192.32-11.88 are significantly smaller than in other giant molecular clouds or filaments, indicating that the star formation therein is greatly suppressed due to stellar feedback.

Theoretical studies have revealed that dust grains are usually moving fast through the turbulent interstellar gas, which could have significant effects upon interstellar chemistry by modifying grain accretion. This effect is investigated in this work on the basis of numerical gas-grain chemical modeling. Major features of the grain motion effect in the typical environment of dark clouds (DC) can be summarised as follows: 1) decrease of gas-phase (both neutral and ionic) abundances and increase of surface abundances by up to 2-3 orders of magnitude; 2) shifts of the existing chemical jumps to earlier evolution ages for gas-phase species and to later ages for surface species by factors of about ten; 3) a few exceptional cases in which some species turn out to be insensitive to this effect and some other species can show opposite behaviors too. These effects usually begin to emerge from a typical DC model age of about 10^5 yr. The grain motion in a typical cold neutral medium (CNM) can help overcome the Coulomb repulsive barrier to enable effective accretion of cations onto positively charged grains. As a result, the grain motion greatly enhances the abundances of some gas-phase and surface species by factors up to 2-6 or more orders of magnitude in the CNM model. The grain motion effect in a typical molecular cloud (MC) is intermediate between that of the DC and CNM models, but with weaker strength. The grain motion is found to be important to consider in chemical simulations of typical interstellar medium.

We study the dynamic changes of electron energy distribution (EED) through systematically analysing the quasi-simultaneous spectral energy distributions (SEDs) of the flat spectrum radio quasar 3C 279 in different states. With Markov chain Monte Carlo (MCMC) technique we model fourteen SEDs of 3C 279 using a leptonic model with a three-parameter log-parabola electron energy distribution (EED). The 14 SEDs can be satisfactorily fitted with the one-zone leptonic model. The observed $\gamma$ rays in 13 states are attributed to Compton scattering of external infrared photons from a surrounding dusty torus. The curved $\gamma$-ray spectrum observed during 2-8 April 2014 is well explained by the external Compton of dust radiation. It is found that there is a clear positive correlation between the curvature parameter $b$ of the EED and the electron peak energy $\gamma'_{\rm pk}$. No clear correlation between $b$ and the synchrotron peak frequency $\nu_{\rm s}$ is found, due to the varied product of Doppler factor and fluid magnetic field from state to state. We interpret the correlation of $b-\gamma'_{\rm pk}$ in a stochastic acceleration scenario. This positive correlation is in agreement with the prediction in the stage when the balance between acceleration and radiative cooling of the electrons is nearly established in the case of the turbulence spectral index $q=2$.

In this paper, we try to explain the observed correlation between the covering factor (CF) of hot dust and the properties of active galactic nuclei (AGNs), e.g., the bolometric luminosity ($L_{\rm{bol}}$) and black hole mass ($M_{\rm{BH}}$). Combining the possible dust distribution in the torus, the angular dependence of the radiation of the accretion disc, and the relation between the critical angle of torus and the Eddington ratio, there are eight possible models investigated in our work. We fit the observed CF with these models to determine the parameters of them. As a result, clumpy torus models can generally explain the observed correlations of tori, while the smooth models fail to produce the required CFs. However, there is still significant scatter even for the best-fitting model, which is the combination of a clumpy torus illuminated by the anisotropic radiation of accretion disc in an AGN. Although some of the observed scatter is due to the uncertainties in measuring $L_{\rm{bol}}$ and $M_{\rm{BH}}$, other factors are required in more realistic model. The models examined in this paper are not necessary to be the physical model of tori. However, the reasonable assumptions selected during this process should be helpful in constructing physical models of tori.

We present the first study of the spatial distribution of star formation in z~0.5 cluster galaxies. The analysis is based on data taken with the Wide Field Camera 3 as part of the Grism Lens-Amplified Survey from Space (GLASS). We illustrate the methodology by focusing on two clusters (MACS0717.5+3745 and MACS1423.8+2404) with different morphologies (one relaxed and one merging) and use foreground and background galaxies as field control sample. The cluster+field sample consists of 42 galaxies with stellar masses in the range 10^8-10^11 M_sun, and star formation rates in the range 1-20 M_sun/yr. Both in clusters and in the field, H{\alpha} is more extended than the rest-frame UV continuum in 60% of the cases, consistent with diffuse star formation and inside out growth. In ~20% of the cases, the H{\alpha} emission appears more extended in cluster galaxies than in the field, pointing perhaps to ionized gas being stripped and/or star formation being enhanced at large radii. The peak of the H{\alpha} emission and that of the continuum are offset by less than 1 kpc. We investigate trends with the hot gas density as traced by the X-ray emission, and with the surface mass density as inferred from gravitational lens models and find no conclusive results. The diversity of morphologies and sizes observed in H_alpha illustrates the complexity of the environmental process that regulate star formation. Upcoming analysis of the full GLASS dataset will increase our sample size by almost an order of magnitude, verifying and strengthening the inference from this initial dataset.

Using N-body simulations, we measure the power spectrum of the effective dark matter density field, which is defined through the modified Poisson equation in $f(R)$ cosmologies. We find that when compared to the conventional dark matter power spectrum, the effective power spectrum deviates more significantly from the $\Lambda$CDM model. For models with $f_{R0}=-10^{-4}$, the deviation can exceed 150\% while the deviation of the conventional matter power spectrum is less than 50\%. Even for models with $f_{R0}=-10^{-6}$, for which the conventional matter power spectrum is very close to the $\Lambda$CDM prediction, the effective power spectrum shows sizeable deviations. Our results indicate that traditional analyses based on the dark matter density field may seriously underestimate the impact of $f(R)$ gravity on galaxy clustering. We therefore suggest the use of the effective density field in such studies. In addition, based on our findings, we also discuss several possible methods of making use of the differences between the conventional and effective dark matter power spectra in $f(R)$ gravity to discriminate the theory from the $\Lambda$CDM model.

Intensive studies have been conducted to understand the anisotropy of solar wind turbulence. However, the anisotropy of Els\"asser variables ($\textbf{Z}^\pm$) in 2D wave-vector space has yet to be investigated. Here we first verify the transformation based on the projection-slice theorem between the power spectral density PSD$_{2D}(k_\parallel,k_\perp )$ and the spatial correlation function CF$_{2D} (r_\parallel,r_\perp )$. Based on the application of the transformation to the magnetic field and the particle measurements from the WIND spacecraft, we investigate the spectral anisotropy of Els\"asser variables ($\textbf{Z}^\pm$), and the distribution of residual energy E$_{R}$, Alfv\'en ratio R$_{A}$ and Els\"asser ratio R$_{E}$ in the $(k_\parallel,k_\perp)$ space. The spectra PSD$_{2D}(k_\parallel,k_\perp )$ of $\textbf{B}$, $\textbf{V}$, and $\textbf{Z}_{major}$ (the larger of $\textbf{Z}^\pm$) show a similar pattern that PSD$_{2D}(k_\parallel,k_\perp )$ is mainly distributed along a ridge inclined toward the $k_\perp$ axis. This is probably the signature of the oblique Alfv\'enic fluctuations propagating outwardly. Unlike those of $\textbf{B}$, $\textbf{V}$, and $\textbf{Z}_{major}$, the spectrum PSD$_{2D}(k_\parallel,k_\perp )$ of $\textbf{Z}_{minor}$ is distributed mainly along the $k_\perp$ axis. Close to the $k_\perp$ axis, $\left| {E}_{R}\right|$ becomes larger while R$_{A}$ becomes smaller, suggesting that the dominance of magnetic energy over kinetic energy becomes more significant at small $k_\parallel$. R$_{E}$ is larger at small $k_\parallel$, implying that PSD$_{2D}(k_\parallel,k_\perp )$ of $\textbf{Z}_{minor}$ is more concentrated along the $k_\perp$ direction as compared to that of $\textbf{Z}_{major}$. The residual energy condensate at small $k_\parallel$ is consistent with simulation results in which E$_{R}$ is spontaneously generated by Alfv\'en wave interaction.

Magnetohydronamic turbulence is believed to play a crucial role in heating the laboratorial, space, and astrophysical plasmas. However, the precise connection between the turbulent fluctuations and the particle kinetics has not yet been established. Here we present clear evidence of plasma turbulence heating based on diagnosed wave features and proton velocity distributions from solar wind measurements by the Wind spacecraft. For the first time, we can report the simultaneous observation of counter-propagating magnetohydrodynamic waves in the solar wind turbulence. Different from the traditional paradigm with counter-propagating Alfv\'en waves, anti-sunward Alfv\'en waves (AWs) are encountered by sunward slow magnetosonic waves (SMWs) in this new type of solar wind compressible turbulence. The counter-propagating AWs and SWs correspond respectively to the dominant and sub-dominant populations of the imbalanced Els\"asser variables. Nonlinear interactions between the AWs and SMWs are inferred from the non-orthogonality between the possible oscillation direction of one wave and the possible propagation direction of the other. The associated protons are revealed to exhibit bi-directional asymmetric beams in their velocity distributions: sunward beams appearing in short and narrow patterns and anti-sunward broad extended tails. It is suggested that multiple types of wave-particle interactions, i.e., cyclotron and Landau resonances with AWs and SMWs at kinetic scales, are taking place to jointly heat the protons perpendicularly and parallel.

We propose a new cosmological test of gravity, by using the observed mass fraction of X-ray emitting gas in massive galaxy clusters. The cluster gas fraction, believed to be a fair sample of the average baryon fraction in the Universe, is a well-understood observable, which has previously mainly been used to constrain background cosmology. In some modified gravity models, such as $f(R)$ gravity, gas temperature in a massive cluster is determined by the effective mass of that cluster, which can be larger than its true mass. On the other hand, X-ray luminosity is determined by the true gas density, which in both modified gravity and $\Lambda$CDM models depends mainly on $\Omega_{\rm b}/\Omega_{\rm m}$ and hence the true total cluster mass. As a result, the standard practice of combining gas temperatures and X-ray surface brightnesses of clusters to infer their gas fractions can, in modified gravity models, lead to a larger - in $f(R)$ gravity this can be $1/3$ larger - value of $\Omega_{\rm b}/\Omega_{\rm m}$ than that inferred from other observations such as the CMB. A quick calculation shows that the Hu-Sawicki $n=1$ $f(R)$ model with $|\bar{f}_{R0}|=3\sim5\times10^{-5}$ is in tension with the gas fraction data of the 42 clusters analysed by Allen et al. (2008). We also discuss the implications for other modified gravity models.

Transient brightenings in the transition region of the Sun have been studied for decades and are usually related to magnetic reconnection. Recently, absorption features due to chromospheric lines have been identified in transition region emission lines raising the question of the thermal stratification during such reconnection events. We analyse data from the Interface Region Imaging Spectrograph (IRIS) in an emerging active region. Here the spectral profiles show clear self-absorption features in the transition region lines of Si\,{\sc{iv}}. While some indications existed that opacity effects might play some role in strong transition region lines, self-absorption has not been observed before. We show why previous instruments could not observe such self-absorption features, and discuss some implications of this observation for the corresponding structure of reconnection events in the atmosphere. Based on this we speculate that a range of phenomena, such as explosive events, blinkers or Ellerman bombs, are just different aspects of the same reconnection event occurring at different heights in the atmosphere.

We introduce the idea of {\it effective} dark matter halo catalog in $f(R)$ gravity, which is built using the {\it effective} density field. Using a suite of high resolution N-body simulations, we find that the dynamical properties of halos, such as the distribution of density, velocity dispersion, specific angular momentum and spin, in the effective catalog of $f(R)$ gravity closely mimic those in the $\Lambda$CDM model. Thus, when using effective halos, an $f(R)$ model can be viewed as a $\Lambda$CDM model. This effective catalog therefore provides a convenient way for studying the baryonic physics, the galaxy halo occupation distribution and even semi-analytical galaxy formation in $f(R)$ cosmologies.

The SEDs of 18 GeV FSRQs are collected and compiled from literature, in which both the jet emission and the accretion disk radiation can be observed, in order to investigate the correlations among their jet power (P_jet), accretion disk luminosity (L_disk), and luminosity of broad line region (BLR, L_BLR). On the basis of the SED fits with the jet radiation and accretion disk radiation models, we calculate P_jet and L_disk. No correlation between P_jet with either L_disk or L_BLR is found. With a sub-sample of L_BLR for 13 GeV FSRQs, it is observed that L_BLR is strongly correlated with their L_disk. We also study the BLR covering factors of the GeV FSRQs in our sample, averagely which are smaller than that of the large samples of radio-loud and radio-quiet quasars. P_jet of some GeV FSRQs is higher than L_disk, but P_jet of all the GeV FSRQs is lower than the accretion power of black hole (BH), which is estimated by \dot{M}c^2=L_disk/0.1, indicating that the total accretion power of BH is sufficient to drive the jets in these sources; however the uncorrelation between L_disk and P_jet of the GeV FSRQs may suggest that their jets are launched by the BZ process via extracting the rotational energy of BH. Using the L_BLR-L_disk relation of the GeV FSRQs, we estimate L_disk of a BL Lac sample with their L_BLR. A comparison of L_BLR and Eddington ratio (L_disk/L_Edd) among BL Lacs, very radio-loud NLS1 galaxies, and FSRQs is also presented. It is found that along with the BL Lac-NLS1-FSRQ sequence L_BLR and L_disk/L_Edd increase, which may correspond to the change of the accretion disk structure and the transformation of the dominant mechanism for jet launching. This is also consistent with the division of their parent populations, i.e., low-excitation radio galaxies and high-excitation radio galaxies.

Measurements of solar wind turbulence reveal the ubiquity of discontinuities. In this study, we investigate how the discontinuities, especially rotational discontinuities (RDs), are formed in magnetohydrodynamic (MHD) turbulence. In a simulation of the decaying compressive three-dimensional (3-D) MHD turbulence with an imposed uniform background magnetic field, we detect RDs with sharp field rotations and little variations of magnetic field intensity as well as mass density. At the same time, in the de Hoffman-Teller (HT) frame, the plasma velocity is nearly in agreement with the Alfv\'{e}n speed, and is field-aligned on both sides of the discontinuity. We take one of the identified RDs to analyze in details its 3-D structure and temporal evolution. By checking the magnetic field and plasma parameters, we find that the identified RD evolves from the steepening of the Alfv\'{e}n wave with moderate amplitude, and that steepening is caused by the nonuniformity of the Alfv\'{e}n speed in the ambient turbulence.

In solar wind, magnetohydrodynamic (MHD) discontinuities are ubiquitous and often found to be at the origin of turbulence intermittency. They may also play a key role in the turbulence dissipation and heating of the solar wind. The tangential (TD) and rotational (RD) discontinuities are the two most important types of discontinuities. Recently, the connection between turbulence intermittency and proton thermodynamics has been being investigated observationally. Here we present numerical results from three-dimensional MHD simulation with pressure anisotropy and define new methods to identify and to distinguish TDs and RDs. Three statistical results obtained about the relative occurrence rates and heating effects are highlighted: (1) RDs tend to take up the majority of the discontinuities along with time; (2) the thermal states embedding TDs tend to be associated with extreme plasma parameters or instabilities, while RDs do not; (3) TDs have a higher average T as well as perpendicular temperature $T_\perp$. The simulation shows that TDs and RDs evolve and contribute to solar wind heating differently. These results will inspire our understanding of the mechanisms that generate discontinuities and cause plasma heating.

In the standard Big-Bang nucleosynthesis (BBN) model, the primordial $^7$Li abundance is overestimated by about a factor of 2--3 comparing to the astronomical observations, so called the pending cosmological lithium problem. The $^7$Be($n$,$\alpha$)$^4$He reaction, which may affect the $^7$Li abundance, was regarded as the secondary important reaction in destructing the $^7$Be nucleus in BBN. However, the thermonuclear rate of $^7$Be($n$,$\alpha$)$^4$He has not been well studied so far. This reaction rate was firstly estimated by Wagoner in 1969, which has been generally adopted in the current BBN simulations and the reaction rate library. This simple estimation involved only a direct-capture reaction mechanism, but the resonant contribution should be also considered according to the later experimental results. In this work, we have revised this rate based on the indirect cross-section data available for the $^4$He($\alpha$,$n$)$^7$Be and $^4$He($\alpha$,$p$)$^7$Li reactions, with the charge symmetry and detailed-balance principle. Our new result shows that the previous rate (acting as an upper limit) is overestimated by about a factor of ten. The BBN simulation shows that the present rate leads to a 1.2\% increase in the final $^7$Li abundance compared to the result using the Wagoner rate, and hence the present rate even worsens the $^7$Li problem. By the present estimation, the role of $^7$Be($n$,$\alpha$)$^4$He in destroying $^7$Be is weakened from the secondary importance to the third, and the $^7$Be($d$,$p$)2$^4$He reaction becomes of secondary importance in destructing $^7$Be.

The energy to desorb atomic oxygen from an interstellar dust grain surface, $E_{\rm des}$, is an important controlling parameter in gas-grain models; its value impacts the temperature range over which oxygen resides on a dust grain. However, no prior measurement has been done of the desorption energy. We report the first direct measurement of $E_{\rm des}$ for atomic oxygen from dust grain analogs. The values of $E_{\rm des}$ are $1660\pm 60$~K and $1850\pm 90$~K for porous amorphous water ice and for a bare amorphous silicate film, respectively, or about twice the value previously adopted in simulations of the chemical evolution of a cloud. We use the new values to study oxygen chemistry as a function of depth in a molecular cloud. For $n=10^4$ cm$^{-3}$ and $G_0$=10$^2$ ($G_0$=1 is the average local interstellar radiation field), the main result of the adoption of the higher oxygen binding energy is that H$_2$O can form on grains at lower visual extinction $A_{\rm V}$, closer to the cloud surface. A higher binding energy of O results in more formation of OH and H$_2$O on grains, which are subsequently desorbed by FUV radiation, with consequences for gas-phase chemistry. For higher values of $n$ and $G_0$, the higher binding energy can lead to a large increase in the column of H$_2$O but a decrease in the column of O$_2$.

The quest to detect prebiotic molecules in space, notably amino acids, requires an understanding of the chemistry involving nitrogen atoms. Hydroxylamine (NH$_2$OH) is considered a precursor to the amino acid glycine. Although not yet detected, NH$_2$OH is considered a likely target of detection with ALMA. We report on an experimental investigation of the formation of hydroxylamine on an amorphous silicate surface via the oxidation of ammonia. The experimental data are then fed into a simulation of the formation of NH$_2$OH in dense cloud conditions. On ices at 14 K and with a modest activation energy barrier, NH$_2$OH is found to be formed with an abundance that never falls below a factor 10 with respect to NH$_3$. Suggestions of conditions for future observations are provided.

In the primordial Big Bang nucleosynthesis (BBN), only the lightest nuclides (D, $^3$He, $^4$He, and $^7$Li) were synthesized in appreciable quantities, and these relics provide us a unique window on the early universe. Currently, BBN simulations give acceptable agreement between theoretical and observed abundances of D and $^4$He, but it is still difficult to reconcile the predicted $^7$Li abundance with the observation for the Galactic halo stars. The BBN model overestimates the primordial $^7$Li abundance by about a factor of three, so called the cosmological lithium problem, a long-lasting pending issue in BBN. Great efforts have been paid in the past decades, however, the conventional nuclear physics seems unable to resolve such problem. It is well-known that the classical Maxwell-Boltzmann (MB) velocity distribution has been usually assumed for nuclei in the Big-Bang plasma. In this work, we have thoroughly investigated the impact of non-extensive Tsallis statistics (deviating from the MB) on thermonuclear reaction rates involved in standard models of BBN. It shows that the predicted primordial abundances of D, $^4$He, and $^7$Li agree very well with those observed ones by introducing a non-extensive parameter $q$. It is discovered that the velocities of nuclei in a hot Big-Bang plasma indeed violate the classical Maxwell-Boltzmann (MB) distribution in a very small deviation of about 6.3--8.2%. Thus, we have for the first time found a new solution to the cosmological lithium problem without introducing any mysterious theories. Furthermore, the implications of non-extensive statistics in other exotic high-temperature and density astrophysical environments should be explored, which might offer new insight into the nucleosynthesis of heavy elements.

In order to investigate the physical and chemical properties of massive star forming cores in early stages, we analyse the excitation and abundance of four organic species, CH3OH, CH3OCH3, HCOOCH3 and CH3CH2CN, toward 29 Extended Green Object (EGO) cloud cores that were observed by our previous single dish spectral line survey. The EGO cloud cores are found to have similar methanol J_3-J_2 rotation temperatures of ~44 K, a typical linear size of ~0.036 pc, and a typical beam averaged methanol abundance of several 10^(-9) (the beam corrected value could reach several 10^(-7)). The abundances of the latter three species, normalized by that of methanol, are found to be correlated also across a large variety of clouds such as EGO cloud cores, hot corinos, massive hot cores and Galactic Center clouds. The chemical properties of the EGO cloud cores lie between that of hot cores and hot corinos. However, the abundances and abundance ratios of the four species can not be satisfactorily explained by recent chemical models either among the EGO cloud cores or among the various types of cloud cores from literature.

We reexamine the screening mechanism in $f(R)$ gravity using N-body simulations. By explicitly examining the relation between the extra scalar field $\delta f_R$ and the gravitational potential $\phi$ in the perturbed Universe, we find that the relation between these two fields plays an important role in understanding the screening mechanism. We show that the screening mechanism in $f(R)$ gravity depends mainly on the depth of the potential well, and find a useful condition for identifying unscreened halos in simulations. We also discuss the potential application of our results to real galaxy surveys.

We provide the most stringent constraint to date on possible deviations from the usually-assumed Maxwell-Boltzmann (MB) velocity distribution for nuclei in the Big-Bang plasma. The impact of non-extensive Tsallis statistics on thermonuclear reaction rates involved in standard models of Big-Bang Nucleosynthesis (BBN) has been investigated. We find that the non-extensive parameter $q$ may deviate by, at most, $|\delta q|$=6$\times$10$^{-4}$ from unity for BBN predictions to be consistent with observed primordial abundances; $q$=1 represents the classical Boltzmann-Gibbs statistics. This constraint arises primarily from the {\em super}sensitivity of endothermic rates on the value of $q$, which is found for the first time. As such, the implications of non-extensive statistics in other astrophysical environments should be explored. This may offer new insight into the nucleosynthesis of heavy elements.

The $^{14}$O($\alpha$,$p$)$^{17}$F reaction is one of the key reactions involved in the breakout from the hot-CNO cycle to the rp-process in type I x-ray bursts (XRBs). The resonant properties in the compound nucleus $^{18}$Ne have been investigated through resonant elastic scattering of $^{17}$F+$p$. The radioactive $^{17}$F beam was separated by the CNS Radioactive Ion Beam separator (CRIB) and bombarded a thick H$_2$ gas target at 3.6 MeV/nucleon. The recoiling light particles were measured by three ${\Delta}$E-E silicon telescopes at laboratory angles of $\theta$$_{lab}$$\approx$3$^\circ$, 10$^\circ$ and 18$^\circ$, respectively. Five resonances at $E_{x}$=6.15, 6.28, 6.35, 6.85, and 7.05 MeV were observed in the excitation functions, and their spin-parities have been determined based on an $R$-matrix analysis. In particular, $J^{\pi}$=1$^-$ was firmly assigned to the 6.15-MeV state which dominates the thermonuclear $^{14}$O($\alpha$,$p$)$^{17}$F rate below 2 GK. As well, a possible new excited state in $^{18}$Ne was observed at $E_{x}$=6.85$\pm$0.11 MeV with tentative $J$=0 assignment. This state could be the analog state of the 6.880 MeV (0$^{-}$) level in the mirror nucleus $^{18}$O, or a bandhead state (0$^+$) of the six-particle four-hole (6$p$-4$h$) band. A new thermonuclear $^{14}$O($\alpha$,$p$)$^{17}$F rate has been determined, and the astrophysical impact of multiple recent rates has been examined using an XRB model. Contrary to previous expectations, we find only modest impact on predicted nuclear energy generation rates from using reaction rates differing by up to several orders of magnitude.

We constrain the neutrino mass in $f(R)$ gravity using the latest observations from the Planck, BAO and BICEP2 data. We find that the measurement on the B-modes can break the degeneracy between the massive neutrinos and the $f(R)$ gravity. We find a non-zero value of the Compton wavelengths $B_{0}$ at a $68\%$ confidence level for the $f(R)$ model in the presence of massive neutrinos when the BICEP2 data is used. Furthermore, the tension on the tensor-to-scalar ratios between the measured values from Plank and BICEP2 is significantly reconciled in our model.

Superhalo electrons appear to be continuously present in the interplanetary medium, even at very quiet times, with a power-law spectrum at energies above $\sim$2 keV. Here we numerically investigate the generation of superhalo electrons by magnetic reconnection in the solar wind source region, using the MHD and test particle simulations for both single X-line reconnection and multiple X-line reconnection. We find that the direct current electric field, produced in the magnetic reconnection region, can accelerate electrons from an initial thermal energy of T $\sim10^5$ K up to hundreds of keV. After acceleration, some of the accelerated electrons, together with the nascent solar wind flow driven by the reconnection, propagate upwards along the newly-opened magnetic field lines into the interplanetary space, while the rest move downwards into the lower atmosphere. Similar to the observed superhalo electrons at 1 AU, the flux of the upward-traveling accelerated electrons versus energy displays a power-law distribution at $\sim$ 2 $-$ 100 keV, $f(E) \sim E^{-\delta}$, with a $\delta$ of $\sim$ 1.5 $-$ 2.4. For single (multiple) X-line reconnection, the spectrum becomes harder (softer) as the anomalous resistivity parameter $\alpha$ (uniform resistivity $\eta$) increases. These modeling results suggest that the acceleration in the solar wind source region may contribute to superhalo electrons.

The current Big-Bang Nucleosynthesis (BBN) model has been constructed based on a nuclear reaction network operating with thermal reactivities of Maxwell-Boltzmann (MB) distribution plasma. However, does the classical MB distribution still hold for the extremely high-temperature (in order of 10$^9$ K) plasma involved in the Big-Bang environment? In this work, we have investigated the impact of non-extensive Tsallis statistics (in $q$-Guassian distribution) on the thermonuclear reaction rates. We show for the first time that the reverse rates are extremely sensitive to the non-extensive $q$ parameter. Such sensitivity does not allow a large deviation of non-extensive distribution from the usual MB distribution. With a newly developed BBN code, the impact of primordial light-element abundances on $q$ values has been studied by utilizing the most recent BBN cosmological parameters and the available nuclear cross-section data. For the first time, we have accurately verified the microscopic MB distribution with the macroscopic BBN theory and bservation. By comparing the recent observed primordial abundances with our predictions, only a tiny deviation of $\pm$6$\times$10$^{-4}$ at most can be allowed for the MB distribution. However, validity of the classical statistics needs to be studied further for the self-gravitating stars and binaries of high-density environment, with the extreme sensitivity of reverse rate on $q$ found here.

New multi-color light curves of the very short period K-type eclipsing binary V1799 Ori were obtained and analyzed with the W-D code. The photometric solutions reveal that the system is a W-type shallow-contact binary with a mass ratio of $q=1.335(\pm0.005)$ and a degree of contact about $f = 3.5(\pm1.1)\%$. In general, the results are in good agreement with which is reported by Samec. The remarkable O'Connell effects in the light curves are well explained by employing star spots on the binary surface, which confirms that the system is active at present. Several new times of light minimum were obtained. All the available times of light minimum were collected, along with the recalculated and new obtained. Applying a least-squares method to the constructed O-C diagram, a new ephemeris was derived for V1799 Ori. The orbital period is found to show a continuous weak increase at a rate of $1.8(\pm0.6)\times10^{-8}$ days$\cdot$yr$^{-1}$. The extreme-shallow contact, together with the period increase, suggests that the binary may be at a critical stage predicted by the TRO theory. \keywords{binaries : close -- binaries : eclipsing -- stars: individual (V1799 Ori)

The $^{18}$Ne($\alpha$,$p$)$^{21}$Na reaction is thought to be one of the key breakout reactions from the hot CNO cycles to the rp-process in type I x-ray bursts. In this work, the resonant properties of the compound nucleus $^{22}$Mg have been investigated by measuring the resonant elastic scattering of $^{21}$Na+$p$. An 89 MeV $^{21}$Na radioactive beam delivered from the CNS Radioactive Ion Beam Separator bombarded an 8.8 mg/cm$^2$ thick polyethylene (CH$_{2}$)$_{n}$ target. The $^{21}$Na beam intensity was about 2$\times$10$^{5}$ pps, with a purity of about 70% on target. The recoiled protons were measured at the center-of-mass scattering angles of $\theta_{c.m.}$$\approx$175.2${^\circ}$, 152.2${^\circ}$, and 150.5${^\circ}$ by three sets of $\Delta E$-$E$ telescopes, respectively. The excitation function was obtained with the thick-target method over energies $E_x$($^{22}$Mg)=5.5--9.2 MeV. In total, 23 states above the proton-threshold in $^{22}$Mg were observed, and their resonant parameters were determined via an $R$-matrix analysis of the excitation functions. We have made several new $J^{\pi}$ assignments and confirmed some tentative assignments made in previous work. The thermonuclear $^{18}$Ne($\alpha$,$p$)$^{21}$Na rate has been recalculated based on our recommended spin-parity assignments. The astrophysical impact of our new rate has been investigated through one-zone postprocessing x-ray burst calculations. We find that the $^{18}$Ne($\alpha$,$p$)$^{21}$Na rate significantly affects the peak nuclear energy generation rate, reaction fluxes, as well as the onset temperature of this breakout reaction in these astrophysical phenomena.

We constrain the neutrino properties in $f(R)$ gravity using the latest observations from cosmic microwave background(CMB) and baryon acoustic oscillation(BAO) measurements. We first constrain separately the total mass of neutrinos $\sum m_\nu$ and the effective number of neutrino species $N_{\rm eff}$. Then we constrain $N_{\rm eff}$ and $\sum m_\nu$ simultaneously. We find $\sum m_\nu<0.462 {\rm eV}$ at a $95\%$ confidence level for the combination of Planck CMB data, WMAP CMB polarization data, BAO data and high-$l$ data from the Atacama Cosmology Telescope and the South Pole Telescope. We also find $N_{\rm eff}=3.32^{+0.54}_{-0.51}$ at a $95\%$ confidence level for the same data set. When constraining $N_{\rm eff}$ and $\sum m_\nu$ simultaneously, we find $N_{\rm eff}=3.58^{+0.72}_{-0.69}$ and $\sum m_\nu<0.860{\rm eV}$ at a $95\%$ confidence level, respectively.

A CO\,2-1 line survey is performed toward a sample of 58 high Galactic latitude post-AGB (pAGB) stars. To complement the observations, a compilation of literature CO\,2-1 line data of known pAGB stars is done. After combining the datasets, CO\,2-1 line data are available for 133 pAGB stars (about 34 per cent of known pAGB stars) among which 44 are detections. The CO line strengths are compared with infrared dust emission for these pAGB stars by defining a ratio between the integrated CO\,2-1 line flux and {\it IRAS} 25\,mu flux density (CO-IR ratio). The relationship between the CO-IR ratio and the {\it IRAS} color C23 (defined with the 25 and 60\,mu flux densities) is called here the CO-IR diagram. The pAGB objects are found to be located between AGB stars and planetary nebulae (PNe), and segregate into three distinctive groups (I, II and III) on the CO-IR diagram. By analyzing their various properties such as chemical types, spectral types, binarity, circumstellar envelope expansion velocities, and pAGB sub-types on the CO-IR diagram, it is argued that the group-I objects are mainly intermediate mass C-rich pAGB stars in early pAGB stage (almost all of the considered carbon rich `21\,mu' stars belong to this group); the group-II objects are massive or intermediate mass pAGB stars which already follow the profound trend of PNe; and the group-III objects are mainly low mass binary pAGB stars with very weak CO\,2-1 line emission (almost all of the considered RV\,Tau variables belong to this group). The CO-IR diagram is proven to be a powerful tool to investigate the co-evolution of circumstellar gas and dust during the short pAGB stage of stellar evolution.

We present a first statistical study of subproton and electron scales turbulence in the terrestrial magnetosheath using the Cluster Search Coil Magnetometer (SCM) waveforms of the STAFF instrument measured in the frequency range [1,180] Hz. It is found that clear spectral breaks exist near the electron scale, which separate two power-law like frequency bands referred to as the dispersive and the electron dissipation ranges. The frequencies of the breaks f_b are shown to be well correlated with the electron gyroscale \rho_e rather than with the electron inertial length de. The distribution of the slopes below fb was found to be narrow and peaks near -2.9, while that of the slopes above fb was found broader, peaks near -5.2 and has values as low as -7.5. This is the first time that such steep power-law spectra are reported in space plasma turbulence. These observations provide strong constraints on theoretical modeling of kinetic turbulence and dissipation in collisionless magnetized plasmas.

In this paper, we study the non-linear matter power spectrum in a specific family of $f(R)$ models that can reproduce the $\Lambda$CDM background expansion history, using high resolution $N$-body simulations based on the {\sc ecosmog} code. We measure the matter power spectrum in the range of $0.05h{\rm Mpc}^{-1}<k<10h{\rm Mpc}^{-1}$ from simulations for our $f(R)$ models and give theoretical explanations to their behaviour and evolution patterns. We also examine the chameleon mechanism for our models and find that it works throughout the cosmic history in dense regions, for our $f(R)$ models with $|f_{R0}|<10^{-4}$. On the other hand, for models with $|f_{R0}|>10^{-3}$, we find no chameleon screening in dense regions at late times ($z<3$), which means that those models could be ruled out due to the factor-of-$1/3$ enhancement to the strength of Newtonian gravity. We also give the best-fit parameters for a generalised PPF fitting formula which works well for the models studied here.

The $^{18}$Ne($\alpha$,$p$)$^{21}$Na reaction provides a pathway for breakout from the hot CNO cycles to the $rp$-process in type I x-ray bursts. To better determine this astrophysical reaction rate, the resonance parameters of the compound nucleus $^{22}$Mg have been investigated by measuring the resonant elastic scattering of $^{21}$Na+$p$. An 89 MeV $^{21}$Na radioactive ion beam was produced at the CNS Radioactive Ion Beam Separator and bombarded an 8.8 mg/cm$^2$ thick polyethylene target. The recoiled protons were measured at scattering angles of $\theta_{c.m.}$$\approx 175 {^\circ}$ and 152${^\circ}$ by three $\Delta E$-$E$ silicon telescopes. The excitation function was obtained with a thick-target method over energies $E_x$($^{22}$Mg)=5.5--9.2 MeV. The resonance parameters have been determined through an $R$-matrix analysis. For the first time, the $J^{\pi}$ values for ten states above the alpha threshold in $^{22}$Mg have been experimentally determined in a single consistent measurement. We have made three new $J^{\pi}$ assignments and confirmed seven of the ten tentative assignments in the previous work. The $^{18}$Ne($\alpha$,$p$)$^{21}$Na reaction rate has been recalculated, and the astrophysical impact of our new rate has been investigated through one-zone postprocessing x-ray burst calculations. We find that the $^{18}$Ne($\alpha$,$p$)$^{21}$Na rate significantly affects the peak nuclear energy generation rate and the onset temperature of this breakout reaction in these phenomena.

Nuclear reactions in stars occur between nuclei in the high-energy tail of the energy distribution and are sensitive to possible deviations from the standard equilibrium thermal-energy distribution, the well-known Maxwell-Boltzmann Distribution (\textsf{MBD}). In a previous paper published in Physics Letters 441B(1998)291, Degl'Innocenti {\it et al}. made strong constrains on such deviations with the detailed helioseismic information of the solar structure. With a small deviation parameterized with a factor exp$[{-\delta (E/kT)^2}]$, it was shown $\delta$ restricted between -0.005 and +0.002. These constrains have been carefully reexamined in the present work. We find that a normalization factor was missed in the previous modified \textsf{MBD}. In this work, the normalization factor $c$ is calculated as a function of $\delta$. It shows the factor $c$ is almost unity within the range 0$< \delta \leq$0.002, which supports the previous conclusion. However, it demonstrates that $\delta$ cannot take a negative value from the normalization point of view. As a result, a stronger constraint on $\delta$ is defined as 0$\leq \delta \leq$0.002. The astrophysical implication on the solar neutrino fluxes is simply discussed based on a positive $\delta$ value of 0.003. The reduction of the $^7$Be and $^8$B neutrino fluxes expected from the modified \textsf{MBD} can possibly shed alternative light on the solar neutrino problem. In addition, the resonant reaction rates for the $^{14}$N($p$,$\gamma$)$^{15}$O reaction are calculated with a standard \textsf{MBD} and a modified \textsf{MBD}, respectively. It shows that the rates are quite sensitive even to a very small $\delta$. This work demonstrates the importance and necessity of experimental verification or test of the well-known \textsf{MBD} at high temperatures.

The solar transition region (TR), in which above the photosphere the tempera- ture increases rapidly and the density drops dramatically, is believed to play an important role in coronal heating and solar wind acceleration. Long-lasting up-flows are present in the upper TR and interpreted as signatures of mass supply to large coronal loops in the quiet Sun. Coronal bright points (BPs) are local heating phenomena and we found a different Doppler-shift pattern at TR and coronal temperatures in one BP, which might be related to the twisted loop system. The dominant energy loss in the lower TR is the Ly-alpha emission. It has been found that most Ly-alpha radiance profiles are stronger in the blue peak, an asymmetry opposite to higher order Lyman lines. This asymmetry is stronger when the downflow in the middle TR is stronger, indicating that the TR flows play an important role in the line formation process. The peak separation of Ly-alpha is found to be larger in coronal holes than in the quiet Sun, reflecting the different magnetic structures and radiation fields between the two regions. The Lyman line profiles are found to be not reversed in sunspot plume and umbra regions, while they are obviously reversed in the surrounding plage region. At TR temperatures, the densities of the sunspot plume and umbra are a factor of 10 lower than of the plage, indicating that the sunspot plasma emitting at TR temperatures is higher and possibly more extended above sunspots than above the plage region.

The Sun is the only star that we can spatially resolve and it can be regarded as a fundamental plasma laboratory of astrophysics. The solar transition region (TR), the layer between the solar chromosphere and corona, plays an important role in solar wind origin and coronal heating. Recent high-resolution observations made by SOHO, TRACE, and Hinode indicate that the TR is highly nonuniform and magnetically structured. Through a combination of spectroscopic observations and magnetic field extrapolations, the TR magnetic structures and plasma properties have been found to be different in coronal holes and in the quiet Sun. In active regions, the TR density and temperature structures also differ in sunspots and the surrounding plage regions. Although the TR is believed to be a dynamic layer, quasi-steady flows lasting from several hours to several days are often present in the quiet Sun, coronal holes, and active regions, indicating some kind of plasma circulation/convection in the TR and corona. The emission of hydrogen Lyman lines, which originates from the lower TR, has also been intensively investigated in the recent past. Observations show clearly that the flows and dynamics in the middle and upper TR can greatly modify the Lyman line profiles.

We present observations of the power spectral anisotropy in wave-vector space of solar wind turbulence, and study how it evolves in interplanetary space with increasing heliocentric distance. For this purpose we use magnetic field measurements made by the Helios-2 spacecraft at three positions between 0.29 and 0.9 AU. To derive the power spectral density (PSD) in $(k_\parallel, k_\bot)$-space based on single-satellite measurements is a challenging task not yet accomplished previously. Here we derive the spectrum $\rm{PSD}_{\rm{2D}}$($\rm{k}_\parallel$, $\rm{k}_\bot$) from the spatial correlation function $\rm{CF}_{\rm{2D}}(r_\parallel, r_\bot)$ by a transformation according to the projection-slice theorem. We find the so constructed PSDs to be distributed in k-space mainly along a ridge that is more inclined toward the $\rm{k}_\bot$ than $\rm{k}_\parallel$ axis, a new result which probably indicates preferential cascading of turbulent energy along the $\rm{k}_\bot$ direction. Furthermore, this ridge of the distribution is found to gradually get closer to the $\rm{k}_\bot$ axis, as the outer scale length of the turbulence becomes larger while the solar wind flows further away from the Sun. In the vicinity of the $\rm{k}_\parallel$ axis, there appears a minor spectral component that probably corresponds to quasi-parallel Alfv\'enic fluctuations. Their relative contribution to the total spectral density tends to decrease with radial distance. These findings suggest that solar wind turbulence undergoes an anisotropic cascade transporting most of its magnetic energy towards larger $\rm{k}_\bot$, and that the anisotropy in the inertial range is radially developing further at scales that are relatively far from the ever increasing outer scale.

We reconstruct an $f(R)$ gravity model that gives rise to the particular $\Lambda$CDM background evolution of the universe. We find well-defined, real-valued analytical forms for the $f(R)$ model to describe the universe both in the early epoch from the radiation to matter dominated eras and the late time acceleration period. We further examine the viability of the derived $f(R)$ model and find that it is viable to describe the evolution of the universe in the past and there does not exist the future singularity in the Lagrangian.

In this work, we further investigate the family of $f(R)$ dark energy models that can exactly mimic the same background expansion history as that of the $\Lambda$CDM model. We study the large scale structure in the $f(R)$ gravity using the full set cosmological perturbation equations. We investigate the structure formation in both the spatially flat and curved Universe. We also confront our model with the latest observations and conduct a Markov Chain Monte Carlo analysis on the parameter space.

A northern subsample of 89 Spitzer GLIMPSE extended green objects (EGOs), the candidate massive young stellar objects, are surveyed for molecular lines in two 1-GHz ranges: 251.5- 252.5 and 260.188-261.188 GHz. A comprehensive catalog of observed molecular line data and spectral plots are presented. Eight molecular species are undoubtedly detected: H13CO+, SiO, SO, CH3OH, CH3OCH3, CH3CH2CN, HCOOCH3, and HN13C. H13CO+ 3-2 line is detected in 70 EGOs among which 37 ones also show SiO 6-5 line, demonstrating their association to dense gas and supporting the outflow interpretation of the extended 4.5 um excess emission. Our major dense gas and outflow tracers (H13CO+, SiO, SO and CH3OH) are combined with our previous survey of 13CO, 12CO and C18O 1-0 toward the same sample of EGOs for a multi-line multi- cloud analysis of line width and luminosity correlations. Good log-linear correlations are found among all considered line luminosities, which requires a universal similarity of density and thermal structures and probably of shock properties among all EGO clouds to explain. It also requires that the shocks should be produced within the natal clouds of the EGOs. Diverse degrees of correlation are found among the line widths. However, both the line width and luminosity correlations tend to progressively worsen across larger cloud subcomponent size-scales, depicting the increase of randomness across cloud subcomponent sizes. Moreover, the line width correlations among the three isotopic CO 1-0 lines show data scatter as linear functions of the line width itself, indicating that the velocity randomness also increases with whole-cloud sizes and has some regularity behind.

We report here the discovery of an optical flare observed in R band from the red-dwarf eclipsing binary CU Cnc whose component stars are at the upper boundary of full convection (M1=0.43 and M2=0.4M0, M0 is the solar mass). The amplitude of the flare is the largest among those detected in R band (~0.52mag) and the duration time is about 73 minutes. As those observed on the Sun, quasi-periodic oscillations were seen during and after the flare. Three more R-band flares were found by follow up monitoring. In total, this binary was monitored photometrically by using R filter for 79.9 hours, which reveals a R-band flare rate about 0.05 flares per hour. These detections together with other strong chromospheric and coronal activities, i.e., very strong H_alpha and H_beta emission features and an EUV and X-ray source, indicate that it has very strong magnetic activity. Therefore, the apparent faintness (~1.4 magnitude in V) of CU Cnc compared with other single red dwarfs of the same mass can be plausibly explained by the high coverage of the dark spots.

We observed four "water fountain" sources in the CO J=3-2 line emission with the Atacama Submillimeter Telescope Experiment (ASTE) 10 m telescope in 2010-2011. The water fountain sources are evolved stars that form high-velocity collimated jets traced by water maser emission. The CO line was detected only from IRAS 16342-3814. The present work confirmed that the ^{12}CO to ^{13}CO line intensity ratio is ~1.5 at the systemic velocity. We discuss the origins of the very low ^{12}CO to ^{13}CO intensity ratio, as possible evidence for the "hot-bottom burning" in an oxygen-rich star, and the CO intensity variation in IRAS 16342-3814.

We review the conformal equivalence in describing the background expansion of the universe by $f(R)$ gravity both in the Jordan frame and the Einstein frame. In the Jordan frame, we present the general analytic expression for $f(R)$ models that have the same expansion history as the $\Lambda$CDM model. This analytic form can provide further insights on how cosmology can be used to test the $f(R)$ gravity at the largest scales. Moreover we present a systematic and self-consistent way to construct the viable $f(R)$ model in Jordan frame using the mass dilation rate function from the Einstein frame through the conformal transformation. In addition, we extend our study to the linear perturbation theories and we further exhibit the equivalence of the $f(R)$ gravity presented in the Jordan frame and Einstein frame in the perturbed space-time. We argue that this equivalence has solid physics root.

We report a survey with the Purple Mountain Observatory (PMO) 13.7-m radio telescope for class I methanol masers from the 95 GHz (8_0 - 7_1 A^+) transition. The 214 target sources were selected by combining information from both the Spitzer GLIMPSE and 1.1 mm BGPS survey catalogs. The observed sources satisfy both the GLIMPSE mid-IR criteria of [3.6]-[4.5]>1.3, [3.6]-[5.8]>2.5, [3.6]-[8.0]>2.5 and 8.0 um magnitude less than 10, and also have an associated 1.1 mm BGPS source. Class I methanol maser emission was detected in 63 sources, corresponding to a detection rate of 29% for this survey. For the majority of detections (43), this is the first identification of a class I methanol maser associated with these sources. We show that the intensity of the class I methanol maser emission is not closely related to mid-IR intensity or the colors of the GLIMPSE point sources, however, it is closely correlated with properties (mass and beam-averaged column density) of the BGPS sources. Comparison of measures of star formation activity for the BGPS sources with and without class I methanol masers indicate that the sources with class I methanol masers usually have higher column density and larger flux density than those without them. Our results predict that the criteria log(S_{int})<-38.0+1.72log(N_{H_{2}}^{beam}) and log(N_{H_{2}}^{beam})>22.1, which utilizes both the integrated flux density (S_{int}) and beam-averaged column density (N_{H_{2}}^{beam}) of the BGPS sources, are very efficient for selecting sources likely to have an associated class I methanol maser. Our expectation is that searches using these criteria will detect 90% of the predicted number of class I methanol masers from the full BGPS catalog (~ 1000), and do so with a high detection efficiency (~75%).

We analyze several data sets obtained by Hinode/EIS and find various types of flows during CMEs and EUV jet eruptions. CME-induced dimming regions are found to be characterized by significant blueshift and enhanced line width by using a single Gaussian fit. While a red-blue (RB) asymmetry analysis and a RB-guided double Gaussian fit of the coronal line profiles indicate that these are likely caused by the superposition of a strong background emission component and a relatively weak (~10%) high-speed (~100 km s-1) upflow component. This finding suggests that the outflow velocity in the dimming region is probably of the order of 100 km s-1, not ~20 km s-1 as reported previously. Density and temperature diagnostics suggest that dimming is primarily an effect of density decrease rather than temperature change. The mass losses in dimming regions as estimated from different methods are roughly consistent with each other and they are 20%-60% of the masses of the associated CMEs. With the guide of RB asymmetry analysis, we also find several temperature-dependent outflows (speed increases with temperature) immediately outside the (deepest) dimming region. In an erupted CME loop and an EUV jet, profiles of emission lines formed at coronal and transition region temperatures are found to exhibit two well-separated components, an almost stationary component accounting for the background emission and a highly blueshifted (~200 km s-1) component representing emission from the erupting material. The two components can easily be decomposed through a double Gaussian fit and we can diagnose the electron density, temperature and mass of the ejecta. Combining the speed of the blueshifted component and the projected speed of the erupting material derived from simultaneous imaging observations, we can calculate the real speed of the ejecta.

We report here the tentative discovery of a Jovian planet in orbit around the rapidly pulsating subdwarf B-type (sdB-type) eclipsing binary NY Vir. By using new determined eclipse times together with those collected from the literature, we detect that the observed-calculated (O-C) curve of NY Vir shows a small-amplitude cyclic variation with a period of 7.9\,years and a semiamplitude of 6.1\,s, while it undergoes a downward parabolic change (revealing a period decrease at a rate of $\dot{P}=-9.2\times{10^{-12}}$). The periodic variation was analyzed for the light-travel time effect via the presence of a third body. The mass of the tertiary companion was determined to be $M_3\sin{i^{\prime}}=2.3(\pm0.3)$\,$M_{Jupiter}$ when a total mass of 0.60\,$M_{\odot}$ for NY Vir is adopted. This suggests that it is most probably a giant circumbinary planet orbiting NY Vir at a distance of about 3.3 astronomical units (AU). Since the rate of period decrease can not be explained by true angular momentum loss caused by gravitational radiation or/and magnetic braking, the observed downward parabolic change in the O-C diagram may be only a part of a long-period (longer than 15 years) cyclic variation, which may reveal the presence of another Jovian planet ($\sim2.5$$M_{Jupiter}$) in the system.

We examine the conformal equivalence between the $f(R)$ gravity and the interacting dark sector model. We review the well-known result that the conformal transformation physically corresponds to the mass dilation which marks the strength of interaction between dark sectors. Instead of modeling f(R) gravity in the Jordan frame, we construct the $f(R)$ gravity in terms of mass dilation function in the Einstein frame. We find that the condition to keep $f(R)$ gravity consistent with CMB observations ensures the energy flow from dark energy to dark matter in the corresponding interacting model, which meets the requirement to alleviate the coincidence problem in the Einstein framework.

We study the interacting dark energy model with time varying dark energy equation of state. We examine the stability in the perturbation formalism and the degeneracy among the coupling between dark sectors, the time-dependent dark energy equation of state and dark matter abundance in the cosmic microwave background radiation. Further we discuss the possible ways to break such degeneracy by doing global fitting using the latest observational data and we get a tight constraint on the interaction between dark sectors.

The reaction rates of the key stellar reaction of $^{14}$O($\alpha$,p)$^{17}$F have been reevaluated. It is thought that the previous 1$^-$ assignment for the 6.15-MeV state is incorrect by a careful reanalysis of the previous experimental data [J. G\'{o}mez del Campo {\it et al.}, Phys. Rev. Lett. {\bf 86}, 43 (2001)]. Most probably, the 6.286-MeV state is the key 1$^-$ state and the 6.15-MeV state is a 2$^-$ one, and hence the resonance at $E_x$=6.286 MeV ($J^{\pi}$=1$^-$) actually dominates the reaction rates in the temperature region of astrophysical interests. The newly calculated reaction rates for the $^{14}$O($\alpha$,p)$^{17}$F reaction are quite different from the previous ones, for instance, it's only about 1/6 of the previous value around 0.4 GK, while it's about 2.4 times larger than the previous value around 2 GK. The astrophysical implications have been briefly discussed based on the present conclusions.

Observations from the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) reveal ubiquitous episodic outflows (jets) with an average speed around 120 km s-1 at temperatures often exceeding a million degree in plume-like structures, rooted in magnetized regions of the quiet solar atmosphere. These outflows are not restricted to the well-known plumes visible in polar coronal holes, but are also present in plume-like structures originating from equatorial coronal holes and quiet-Sun regions. Outflows are also visible in the "interplume" regions throughout the atmosphere. Furthermore, the structures traced out by these flows in both plume and inter-plume regions continually exhibit transverse (Alfvenic) motion. Our finding suggests that high-speed outflows originate mainly from the magnetic network of the quiet Sun and coronal holes, and that the plume flows observed are highlighted by the denser plasma contained therein. These outflows might be an efficient means to provide heated mass into the corona and serve as an important source of mass supply to the solar wind. We demonstrate that the quiet-Sun plume flows can sometimes significantly contaminate the spectroscopic observations of the adjacent coronal holes - greatly affecting the Doppler shifts observed, thus potentially impacting significant investigations of such regions.

The new discovered short-period close binary star, XY LMi, was monitored photometrically since 2006. It is shown that the light curves are typical EW-type and show complete eclipses with an eclipse duration of about 80 minutes. By analyzing the complete B, V, R, and I light curves with the 2003 version of the W-D code, photometric solutions were determined. It is discovered that XY LMi is a high fill-out, extreme mass ratio overcontact binary system with a mass ratio of q=0.148 and a fill-out factor of f=74.1%, suggesting that it is on the late evolutionary stage of late-type tidal-locked binary stars. As observed in other overcontact binary stars, evidence for the presence of two dark spots on both components are given. Based on our 19 epoches of eclipse times, it is found that the orbital period of the overcontact binary is decreasing continuously at a rate of dP/dt=-1.67\times10^{-7}\,days/year, which may be caused by the mass transfer from the primary to the secondary or/and angular momentum loss via magnetic stellar wind. The decrease of the orbital period may result in the increase of the fill-out, and finally, it will evolve into a single rapid-rotation star when the fluid surface reaching the outer critical Roche Lobe.

Using the precise times of mid-egress of the eclipsing polar HU Aqr, we discovered that this polar is orbited by two or more giant planets. The two planets detected so far have masses of at least 5.9 and 4.5\,M_{Jup}. Their respective distances from the polar are 3.6 AU and 5.4 AU with periods of 6.54 and 11.96 years, respectively. The observed rate of period decrease derived from the downward parabolic change in O-C curve is a factor 15 larger than the value expected for gravitational radiation. This indicates that it may be only a part of a long-period cyclic variation, revealing the presence of one more planet. It is interesting to note that the two detected circumbinary planets follow the Titus-Bode law of solar planets with n=5 and 6. We estimate that another 10 years of observations will reveal the presence of the predicted third planet.

Cosmological analysis based on currently available observations are unable to rule out a sizeable coupling between dark energy and dark matter. However, the signature of the coupling is not easy to grasp, since the coupling is degenerate with other cosmological parameters, such as the dark energy equation of state and the dark matter abundance. We discuss the possible ways to break such degeneracy. Based on the perturbation formalism, we carry out the global fitting by using latest observational data and get a tight constraint on the interaction between dark sectors. We find that the appropriate interaction can alleviate the coincidence problem.

We report a spectral-line survey of the extreme carbon star IRC+10216 carried out between 293.9 and 354.8 GHz with the Submillimeter Array. A total of 442 lines were detected, more than 200 for the first time; 149 are unassigned. Maps at an angular resolution of ~3" were obtained for each line. A substantial new population of narrow lines with an expansion velocity of ~4 km/s (i.e. ~ 30% of the terminal velocity) was detected. Most of these are attributed to rotational transitions within vibrationally excited states, emitted from energy levels above the v=0, J=0 ground state with excitation energy of 1000-3000 K. Emission from these lines appears to be centered on the star with an angular extent of <1". We use multiple transitions detected in several molecules to derive physical conditions in this inner envelope of IRC+10216.

Based on perturbation theory, we study the dynamics of how dark matter and dark energy in the collapsing system approach dynamical equilibrium while interacting. We find that the interaction between dark sectors cannot ensure the dark energy to fully cluster along with dark, leading to the energy non-conservation problem in the collapsing system We examine the cluster number counts dependence on the interaction between dark sectors. Furthermore, we analyze how dark energy inhomogeneities affect cluster abundances. It is shown that cluster number counts can provide specific signature of dark sectors interaction and dark energy inhomogeneities.

Properties of proton resonances in $^{18}$Ne have been investigated efficiently by utilizing a technique of proton resonant elastic scattering with a $^{17}$F radioactive ion (RI) beam and a thick proton target. A 4.22~MeV/nucleon $^{17}$F RI beam was produced via a projectile-fragmentation reaction, and subsequently separated by a Radioactive Ion Beam Line in Lanzhou ({\tt RIBLL}). Energy spectra of the recoiled protons were measured by two sets of $\Delta$E-E silicon telescope at center-of-mass scattering angles of $\theta_{c.m.}$$\approx$175${^\circ}$$\pm$5${^\circ}$, $\theta_{c.m.}$$\approx$152${^\circ}$$\pm$8${^\circ}$, respectively. Several proton resonances in $^{18}$Ne were observed, and their resonant parameters have been determined by an $R$-matrix analysis of the differential cross sections in combination with the previous results. The resonant parameters are related to the reaction-rate calculation of the stellar $^{14}$O($\alpha$,$p$)$^{17}$F reaction, which was thought to be the breakout reaction from the hot CNO cycles into the $rp$-process in x-ray bursters. Here, $J^\pi$=(3$^-$, 2$^-$) are tentatively assigned to the 6.15-MeV state which was thought the key 1$^-$ state previously. In addition, a doublet structure at 7.05 MeV are tentatively identified, and its contribution to the resonant reaction rate of $^{14}$O($\alpha$,$p$)$^{17}$F could be enhanced by at least factors of about 4$\sim$6 in comparison with the previous estimation involving only a singlet. The present calculated resonant rates are much larger than those previous values, and it may imply that this breakout reaction could play a crucial role under x-ray bursters conditions.

Why is there little dust in elliptical galaxies? Firstly, the universe is rational and the rationality is the proportional density distribution in any independent galaxy. Proportion means that the density distribution of stars is orderly. For example, there are four giants standing in array. Their heights are A, B, C, D, respectively, and A, B stand in the first row from left to right, C, D in the second row from left to right. Proportion requires that A divided by C be equal to B divided by D. Rational structures are smooth 2- (or 3-) dimensional structures, which are either circularly (spherically) symmetric with respect to the center point, or bilaterally symmetric. Secondly, spiral galaxy arms are linearly shaped, which are neither circularly symmetric nor bilaterally symmetric. Therefore, arms are not rational structures. Arm patterns exist only in spiral galaxies and they are weak compared with the main disk structure. Therefore, the presence of arm structure is the disturbance to the rational structure. Because arm patterns exist only in spiral galaxies and only spiral galaxies present dust, I come to the critical answer to the above question: any disturbance to rational structure produces cosmic dust. From now on we call the disturbance to rational structure the cosmic impulse. Therefore, the universe is originated not only rationally but also impulsively.

The interaction between interplanetary small-scale magnetic flux ropes and the magnetic field in the ambient solar wind is an important topic to understand- ing the evolution of magnetic structures in the heliosphere. Through a survey of 125 previously reported small flux ropes from 1995 to 2005, we find that 44 of them reveal clear signatures of Alfvenic fluctuations, and thus classify them into Alfven wave trains rather than flux ropes. Signatures of magnetic reconnection, generally including a plasma jet of ~30 km/s within a magnetic field rotational region, are clearly present at boundaries of about 42% of the flux ropes and 14% of the wave trains. The reconnection exhausts are often observed to show a local increase in the proton temperature, density and plasma beta. About 66% of the reconnection events at flux rope boundaries are associated with a magnetic field shear angle larger than 90 deg and 73% of them reveal a decrease by 20% or more in the magnetic field magnitude, suggesting a dominance of anti-parallel reconnec- tion at flux rope boundaries. The occurrence rate of magnetic reconnection at flux rope boundaries through the year of 1995 to 2005 is also investigated and we find that it is relatively low around solar maximum and much higher when ap- proaching solar minima. The average magnetic field depression and shear angle for reconnection events at flux rope boundaries also reveal a similar trend from 1995 to 2005. Our results demonstrate for the first time that boundaries of a substantial fraction of small-scale flux ropes have properties similar to those of magnetic clouds, in the sense that both of them exhibit signatures of magnetic reconnection. The observed reconnection signatures could be related either to the formation of small flux ropes, or to the interaction between flux ropes and the interplanetary magnetic fields.

Persistent outflows have recently been detected at boundaries of some active regions. Although these outflows are suggested to be possible sources of the slow solar wind, the nature of these outflows is poorly understood. Through an analysis of an image sequence obtained by the X-Ray Telescope onboard the Hinode spacecraft, we found that quasi-period outflows are present in the boundary of an active region. The flows are observed to occur intermittently, often with a period of 5-10 minutes. The projected flow speed can reach more than 200 km/s, while its distribution peaks around 50 km/s. This sporadic high-speed outflow may play an important role in the mass loading process of the slow solar wind. Our results may imply that the outflow of the slow solar wind in the boundary of the active region is intermittent and quasi-periodic in nature.

We study horizontal supergranule-scale motions revealed by TRACE observation of the chromospheric emission, and investigate the coupling between the chromosphere and the underlying photosphere. A highly efficient feature-tracking technique called balltracking has been applied for the first time to the image sequences obtained by TRACE (Transition Region and Coronal Explorer) in the passband of white light and the three ultraviolet passbands centered at 1700 {\AA}, 1600 {\AA}, and 1550 {\AA}. The resulting velocity fields have been spatially smoothed and temporally averaged in order to reveal horizontal supergranule-scale motions that may exist at the emission heights of these passbands. We find indeed a high correlation between the horizontal velocities derived in the white-light and ultraviolet passbands. The horizontal velocities derived from the chromospheric and photospheric emission are comparable in magnitude. The horizontal motions derived in the UV passbands might indicate the existence of a supergranule-scale magnetoconvection in the chromosphere, which may shed new light on the study of mass and energy supply to the corona and solar wind at the height of the chromosphere. However, it is also possible that the apparent motions reflect the chromospheric brightness evolution as produced by acoustic shocks which might be modulated by the photospheric granular motions in their excitation process, or advected partly by the supergranule-scale flow towards the network while propagating upward from the photosphere. To reach a firm conclusion, it is necessary to investigate the role of granular motions in the excitation of shocks through numerical modeling, and future high-cadence chromospheric magnetograms must be scrutinized.

Mach Principle presents the absolute universe. For example, when Einstein stood on the ground and relaxed, his arms fell down naturally. However, if he rotated his body then his arms were lifted up as the rotation was faster and faster. Mach principle is that the matter of the whole universe can affect local dynamic systems. That is, the matter of the whole universe sets up the local absolute reference frames. However, Einstein's general theory of relativity is against the absolute reference frames of Mach Principle.

The starting point of quantum mechanics is the relationship between energy and momentum: energy is proportional to the squared momentum! As a result, energy and momentum have not been treated equally. The wave equation required by quantization is a differential equation which accordingly, does not treat time and space equally. As a result, the Planck constant is not canceled out from the two sides of the equation. Dr. He's gravity is the local bending of background space-time which as suggested by Einstein, can be described by a differential form of second order, and treats time and space equally. Therefore, the Planck constant is completely canceled out in the resulting wave equation. In other words, the quantization of gravity simply does not need the Planck constant!

The reaction rates of the key stellar reaction of $^{14}$O($\alpha$,p)$^{17}$F have been reexamined. The previous conclusion, the 6.15-MeV state ($J^{\pi}$=1$^-$) dominating this reaction rate, has been overthrown by a careful reanalysis of the previous experimental data [J. G\'{o}mez del Campo {\it et al.}, Phys. Rev. Lett. {\bf 86}, 43 (2001)]. According to the present $R$-matrix analysis, the previous 1$^-$ assignment for the 6.15-MeV state is definitely wrong. Most probably, the 6.286-MeV state is the 1$^-$ state and the 6.15-MeV state is a 3$^-$ one, and hence the resonance at $E_x$=6.286 MeV ($J^{\pi}$=1$^-$) dominates the reaction rates in the temperature region of astrophysical interests. The newly calculated reaction rates for the $^{14}$O($\alpha$,p)$^{17}$F reaction are quite different from the previous ones, for instance, it's only about 1/6 of the previous value around 0.4 GK, while it's about 2.4 times larger than the previous value around 2 GK. The astrophysical implications have been discussed based on the present conclusions.

The origin of the solar wind is one of the most important unresolved prob- lems in space and solar physics. We report here the first spectroscopic signatures of the nascent fast solar wind on the basis of observations made by the EUV Imaging Spectrometer (EIS) on Hinode in a polar coronal hole, in which patches of blueshift are clearly present on dopplergrams of coronal emission lines with a formation temperature of lg(T/K)=6.0. The corresponding upflow is associated with open field lines in the coronal hole, and seems to start in the solar transition region and becomes more prominent with increasing temperature. This temperature-dependent plasma outflow is interpreted as evidence of the nascent fast solar wind in the polar coronal hole. The patches with significant upflows are still isolated in the upper transition region but merge in the corona, in agreement with the scenario of solar wind outflow being guided by expanding magnetic funnels.

We report the first systematic survey of molecular lines (including HCO+ (1-0) and 12CO, 13CO, C18O (1-0) lines at 3 mm band) towards a new sample of 88 massive young stellar object (MYSO) candidates associated with ongoing outflows (known as extended green objects or EGOs) identified from the Spitzer GLIMPSE survey in the northern hemisphere with the PMO-13.7 m radio telescope. By analyzing the asymmetries of the optically thick line HCO+ for 69 of 72 EGOs with HCO+ detection, we found 29 sources with blue asymmetric profiles and 19 sources with red asymmetric profiles. This results in a blue excess of 0.14, seen as a signature of collapsing cores in the observed EGO sample. The relatively small blue excess measured in our full sample due to that the observed EGOs are mostly dominated by outflows and at an earlier evolutionary phase associated with IRDCs and 6.7 GHz methanol masers. The physical properties of clouds surrounding EGOs derived from CO lines are similar to those of massive clumps wherein the massive star forming cores associated with EGOs possibly embedded. The infall velocities and mass infall rates derived for 20 infall candidates are also consistent with the typical values found in MYSOs. Thus our observations further support the speculation of Cyganowski et al. (2008) that EGOs trace a population with ongoing outflow activity and active rapid accretion stage of massive protostellar evolution from a statistical view, although there maybe have limitations due to single-pointing survey with a large beam.

The well-known astrophysical resonant reaction rate equations for an isolated narrow resonance have been reexamined. The validity of those `look reliable' assumptions used in deriving the analytic reaction rate equations has been checked, and the reality is they only hold for certain circumstances. Importantly an integration-range issue hidden in the classical integration equations has been revealed and it hints us not to use those analytic equations in the future. This work may influence all those work in which the classical narrow resonant equations were used for calculating the stellar reaction rates, especially at low stellar temperatures, and may affect the results of stellar evolution and nucleosynthesis network calculations considerably.

We investigate the physical meaning of the prominent blue shifts of Ne VIII, which is observed to be associated with quiet-Sun network junctions (boundary intersections), through data analyses combining force-free-field extrapolations with EUV spectroscopic observations. For a middle-latitude region, we reconstruct the magnetic funnel structure in a sub-region showing faint emission in EIT-Fe 195. This funnel appears to consist of several smaller funnels that originate from network lanes, expand with height and finally merge into a single wide open-field region. However, the large blue shifts of Ne VIII are generally not associated with open fields, but seem to be associated with the legs of closed magnetic loops. Moreover, in most cases significant upflows are found in both of the funnel-shaped loop legs. These quasi-steady upflows are regarded as signatures of mass supply to the coronal loops rather than the solar wind. Our observational result also reveals that in many cases the upflows in the upper transition region (TR) and the downflows in the middle TR are not fully cospatial. Based on these new observational results, we suggest different TR structures in coronal holes and in the quiet Sun.

We observed CO J=3-2 emission from the "water fountain" sources, which exhibit high-velocity collimated stellar jets traced by water maser emission, with the Atacama Submillimeter Telescope Experiment (ASTE) 10 m telescope. We detected the CO emission from two sources, IRAS 16342-3814 and IRAS 18286-0959. The IRAS 16342-3814 CO emission exhibits a spectrum that is well fit to a Gaussian profile, rather than to a parabolic profile, with a velocity width (FWHM) of 158+/-6 km/s and an intensity peak at VLSR = 50+/-2 km/s. The mass loss rate of the star is estimated to be ~2.9x10^-5 M_sun/yr. Our morpho-kinematic models suggest that the CO emission is optically thin and associated with a bipolar outflow rather than with a (cold and relatively small) torus. The IRAS 18286-0959 CO emission has a velocity width (FWHM) of 3.0+/-0.2 km/s, smaller than typically seen in AGB envelopes. The narrow velocity width of the CO emission suggests that it originates from either an interstellar molecular cloud or a slowly-rotating circumstellar envelope that harbors the water maser source.

The prominent blue shifts of Ne viii associated with the junctions of the magnetic network in the quiet Sun are still not well understood. By comparing the coronal magnetic-field structures as obtained by a potential-field reconstruction with the conspicuous blue-shift patches on the dopplergram of Ne viii as observed in an equatorial quiet-Sun region, we find that most of the regions with significant upflow are associated with the funnel-like legs of magnetic loops and co-spatial with increments of the line width. These quasi-steady upflows can be regarded as the signatures of mass supply to coronal loops. By using the square-root of the line intensity as a proxy for the plasma density, the mass flux of the upflow in each funnel can be estimated. We find that the mass flux is anti-correlated with the funnel's expansion factor as determined from the extrapolated magnetic field. One of the loop systems is associated with a coronal bright point, which was observed by several instruments and exhibited various morphologies in different wavelengths and viewing directions. A remarkable agreement between its magnetic structure and the associated EUV emission pattern was found, suggesting an almost potential-field nature of the coronal magnetic field. We also report the direct detection of a small-scale siphon flow by both STEREO satellites. However, this transient siphon flow occurred in a weak mixed-polarity-field region, which was outside the adjacent magnetic funnel, and thus it is perhaps not related to plasma upflow in the funnel. Based on these observations, we suggest that at upper-TR temperatures the dominant flows in quiet-Sun coronal loops are long-lasting upflows rather than siphon flows. We also discuss the implications for coronal heating and unresolved magnetic structures.

Dark energy interacting with dark matter is a promising model to solve the cosmic coincidence problem. We study the signature of such interaction on large scale cosmic microwave background (CMB) temperature anisotropies. Based on the detail analysis in perturbation equations of dark energy and dark matter when they are in interaction, we find that the large scale CMB, especially the late Integrated Sachs Wolfe effect, is a useful tool to measure the coupling between dark sectors. We also discuss the possibility to detect the coupling by cross-correlating CMB maps with tracers of the large scale structure. We finally perform the global fitting to constrain the coupling by using the CMB power spectrum data together with other observational data. We find that in the $1\sigma$ range, the constrained coupling between dark sectors can solve the coincidence problem.

We study the emission heights of the coronal bright points (BPs) above the photosphere in the bipolar magnetic loops that are apparently associated with them. As BPs are seen in projection against the disk their true emission heights are unknown. The correlation of the BP locations on the Fe XII radiance map from EIT with the magnetic field features (in particular neutral lines) was investigated in detail. The coronal magnetic field was determined by an extrapolation of the photospheric field to different altitudes above the disk. It was found that most BPs sit on or near a photospheric neutral line, but that the emission occurs at a height of about 5 Mm. Some BPs, while being seen in projection, still seem to coincide with neutral lines, although their emission takes place at heights of more than 10 Mm. Such coincidences almost disappear for emissions above 20 Mm. We also projected the upper segments of the 3-D magnetic field lines above different heights, respectively, on to the x-y plane. The shape of each BP was compared with the respective field-line segment nearby. This comparison suggests that most coronal BPs are actually located on the top of their associated magnetic loops. Finally, we calculated for each selected BP region the correlation coefficient between the Fe XII intensity enhancement and the horizontal component of the extrapolated magnetic field vector at the same x-y position in planes of different heights, respectively. We found that for almost all the BP regions we studied the correlation coefficient, with increasing height, increases to a maximal value and then decreases again. The height corresponding to this maximum was defined as the correlation height, which for most bright points was found to range below 20 Mm.

We present a general formalism to study the growth of dark matter perturbations when dark energy perturbations and interactions between dark sectors are present. We show that dynamical stability of the growth of structure depends on the type of coupling between dark sectors. By taking the appropriate coupling to ensure the stable growth of structure, we observe that the effect of the dark sectors' interaction overwhelms that of dark energy perturbation on the growth function of dark matter perturbation. Due to the influence of the interaction, the growth index can differ from the value without interaction by an amount within the observational sensibility, which provides a possibility to disclose the interaction between dark sectors through future observations on the growth of large structure.

The radiance and Doppler-shift distributions across the solar network provide observational constraints of two-dimensional modeling of transition-region emission and flows in coronal funnels. Two different methods, dispersion plots and average-profile studies, were applied to investigate these distributions. In the dispersion plots, we divided the entire scanned region into a bright and a dark part according to an image of Fe xii; we plotted intensities and Doppler shifts in each bin as determined according to a filtered intensity of Si ii. We also studied the difference in height variations of the magnetic field as extrapolated from the MDI magnetogram, in and outside network. For the average-profile study, we selected 74 individual cases and derived the average profiles of intensities and Doppler shifts across the network. The dispersion plots reveal that the intensities of Si ii and C iv increase from network boundary to network center in both parts. However, the intensity of Ne viii shows different trends, namely increasing in the bright part and decreasing in the dark part. In both parts, the Doppler shift of C iv increases steadily from internetwork to network center. The average-profile study reveals that the intensities of the three lines all decline from the network center to internetwork region. The binned intensities of Si ii and Ne viii have a good correlation. We also find that the large blue shift of Ne viii does not coincide with large red shift of C iv. Our results suggest that the network structure is still prominent at the layer where Ne viii is formed in the quiet Sun, and that the magnetic structures expand more strongly in the dark part than in the bright part of this quiet Sun region.

We performed a systematic study of the Doppler shifts and electron densities measured in an EUV bright point (hereafter BP) observed in more than 10 EUV lines with formation temperatures from log (T/K) p 4.5 to 6.3. Those parts of a BP seen in transition region and coronal lines are defined as its cool and hot components, respectively. We find that the transition from cool to hot occurs at a temperature around log (T/K) p 5.7. The two components of the BP reveal a totally different orientation and Doppler-shift pattern, which might result from a twist of the associated magnetic loop system. The analysis of magnetic field evolution and topology seems to favor a two-stage heating process, in which magnetic cancellation and separator reconnection are powering, respectively, the cool and hot components of the BP. We also found that the electron densities of both components of the BP are higher than those of the surrounding quiet Sun, and comparable to or smaller than active region densities.

We study the height variations of the sizes of chromospheric and transition-region features in a small coronal hole and the adjacent quiet Sun, considering images of the intensity, Doppler shift, and non-thermal motion of ultraviolet emission lines as measured by SUMER, together with the magnetic field as obtained by extrapolation from photospheric magnetograms. In order to estimate the characteristic sizes of the different features present in the chromosphere and transition region, we have calculated the autocorrelation function for the images as well as the corresponding extrapolated magnetic field at different heights. The HWHM of the autocorrelation function is considered to be the characteristic size of the feature shown in the corresponding image. Our results indicate that, in both the coronal hole and quiet Sun, the HWHM of the intensity image is larger than that of the images of Doppler-shift and non-thermal width at any given altitude. The HWHM of the intensity image is smaller in the chromosphere than in the TR, where the sizes of intensity features of lines at different temperatures are almost the same. But in the upper part of the transition region, the intensity size increases more strongly with temperature in the coronal hole than in the quiet Sun. We also studied the height variations of the HWHM of the magnetic field magnitude B and its component |Bz|, and found they are equal to each other at a certain height below 40 Mm in the coronal hole. The height variations of the HWHM of |Bz/B| seem to be consistent with the temperature variations of the intensity size. Our results suggest that coronal loops are much lower, and magnetic structures expand through the upper TR and lower corona much more strongly with height in the coronal hole than in the quiet Sun.

Aims. The physical implication of large blue shift of Ne viii in the quiet Sun region is investigated in this paper. Methods. We compare the significant Ne viii blue shifts, which are visible as large blue patches on the Doppler-shift map of a middlelatitude quiet-Sun region observed by SUMER, with the coronal magnetic-field structures as reconstructed from a simultaneous photospheric magnetogram by means of a force-free-field extrapolation. Results. We show for the first time that coronal funnels also exist in the quiet Sun. The region studied contains several small funnels that originate from network lanes, expand with height and finally merge into a single wide open-field region. However, the large blue shifts of the Ne viii line are not generally associated with funnels. A comparison between the projections of coronal loops onto the solar x-y-plane and the Ne viii dopplergram indicates that there are some loops that reveal large Ne viii blue shifts in both legs, and some loops with upflow in one and downflow in the other leg. Conclusions. Our results suggest that strong plasma outflow, which can be traced by large Ne viii blue shift, is not necessarily associated with the solar wind originating in coronal funnels but appears to be a signature of mass supply to coronal loops. Under the assumption that the measured Doppler shift of the Ne viii line represents the real outflow velocity of the neon ions being markers of the proton flow, we estimate the mass supply rate to coronal loops to be about 10\^{34} s\^{-1}.

From the Main galaxy sample of the Sloan Digital Sky Survey Data Release 6 (SDSS DR6), we construct two volume-limited samples above and below the value of M, to explore the difference of the environmental dependence of galaxy properties between galaxies above and below the value of M . We measure the local three-dimensional galaxy density in a comoving sphere with a radius of the distance to the 5th nearest galaxy for each galaxy, and compare basic properties of galaxies in the lowest density regime with those of galaxies in the densest regime. It is found that the galaxy luminosity strongly depend on local environments only for galaxies above M, but this dependence is very weak for galaxies below M . It is worth noting that g-r color, concentration index ci and galaxy morphologies strongly depend on local environments for all galaxies with different luminosities. This shows that M is an characteristic parameter only for the environmental dependence of galaxy luminosity.

We consider perturbations in a cosmological model with a small coupling between dark energy and dark matter. We prove that the stability of the curvature perturbation depends on the type of coupling between dark sectors. When the dark energy is of quintessence type, if the coupling is proportional to the dark matter energy density, it will drive the instability in the curvature perturbations; however if the coupling is proportional to the energy density of dark energy, there is room for the stability in the curvature perturbations. When the dark energy is of phantom type, the perturbations are always stable, no matter whether the coupling is proportional to the one or the other energy density.

[Context] Water fountain stars are very young post-AGB stars with high velocity water maser jets. They are the best objects to study the onset of bipolar jets from evolved stars due to their young dynamical ages. [Methods] We use the Arizona Radio Observatory 10m telescope to observe the CO J=2-1 line and compare the line parameters with that of masers. [Results] We report the detection of 12CO and 13COJ=2-1 lines from IRAS 16342-3814. The inferred 12CO mass loss rate is an order of magnitude lower than the infrared and OH mass loss rates, indicating a very cold and thick O-rich circumstellar envelope around the star. We also find a 12CO expansion velocity of Vexp = 46 +- 1 km/s that is too high for an AGB wind and confirm the systemic velocity of 44 +- 1 km/s. In addition we measure a very low 12CO/13CO line ratio of 1.7. [Conclusions] The first detection of CO lines has provided a new way to investigate the water fountain stars. Given the high expansion velocity of the CO gas and its relation to maser velocities, we infer that the CO emission region is co-located with the OH mainline masers in the warm base of the optical bipolar lobes, while the high velocity OH1612MHz and H2O masers are located in the side walls and at the farthest ends of the bipolar lobes, respectively. Further observations are highly desired to understand the very low 12CO/13CO line ratio.

We examine the effects of possible phenomenological interactions between dark energy and dark matter on cosmological parameters and their efficiency in solving the coincidence problem. We work with two simple parameterizations of the dynamical dark energy equation of state and the constant dark energy equation of state. Using observational data coming from the new 182 Gold type Ia supernova samples, the shift parameter of the Cosmic Microwave Background given by the three-year Wilkinson Microwave Anisotropy Probe observations, and the baryon acoustic oscillation measurement from the Sloan Digital Sky Survey, we perform a statistical joint analysis of different forms of phenomenological interactions between dark energy and dark matter.

We present the results of our spectral line surveys in the 2 mm and 1.3 mm windows toward the carbon rich envelope of IRC +10216. Totally 377 lines are detected, among which 360 lines are assigned to 57 known molecules (including 29 rare isotopomers and 2 cyclic isomers). Only 17 weak lines remain unidentified. Rotational lines of isotopomers 13CCH and HN13C are detected for the first time in IRC +10216. The detection of the formaldehyde lines in this star is also confirmed. Possible abundance difference among the three 13C substituted isotopic isomers of HC3N is reported. Isotopic ratios of C and O are confirmed to be non-solar while those of S and Si to be nearly solar. Column densities have been estimated for 15 molecular species. Modified spectroscopic parameters have been calculated for NaCN, Na13CN, KCN and SiC2. Transition frequencies from the present observations were used to improve the spectroscopic parameters of Si13CC, 29SiC2 and 30SiC2.

Gravity whose nature is fundamental to the understanding of solar system, galaxies and the structure and evolution of the Universe, is theorized by the assumption of curved spacetime, according to Einstein`s general theory of relativity (EGR). Particles move on curved spacetime along straight lines (geodesics). In the last year, I proposed the mirrored version of EGR, the flat-spacetime general relativity (FGR), in which particles move along curved lines on flat spacetime. This puts gravitational study back to the traditional Lagrangian formulation. In fact, all claimed accurate verification of general relativity is the verification of FGR, because people when confronting GR to observational data, calculate time, distance, or angle by directly using the coordinates in Schwarzschild solution or in post Newtonian formulation. For example, people calculate the angles by directly using the coordinate $\phi$. However, only when spacetime is flat does there exists one coordinate system which has direct meaning of time, distance, angle, and vice verse. This is the famous Riemann theorem. Therefore, the more claims are made that classical tests of general relativity fit data with great accuracy, the more falsified is the curved-spacetime assumption. People made three such specious claims to EGR as collected in the present paper. However, FGR predicts observationally verified results consistently for solar system, galaxies, and the universe on the whole. I show that the possibility of curved spacetime, black holes, and big bang is less than one billionth. An experiment is proposed whose results will completely decide the fate of curved spacetime assumption. with the original article `Einstein`s Geometrization vs. Holonomic Cancellation of Gravity via Spatial Coordinate-rescale` attached.

Using the data coming from the new 182 Gold type Ia supernova samples, the baryon acoustic oscillation measurement from the Sloan Digital Sky Survey and the H(z) data, we have performed a statistical joint analysis of the DGP brane-world model with a high curvature Gauss-Bonnet term in the bulk. Consistent parameters estimations show that the Gauss-Bonnet-Induced Gravity model is a viable candidate to explain the observed acceleration of our universe.

The goal is to build up a simple dynamical model for the out-flowing circumstellar envelope around AGB stars in a wide binary system to explore the parameter dependence of the geometrical characteristics of column density patterns. An AGB star in a wide binary system is considered as a 3-D piston model that can induce a 3-D quasi-spherical density structure in the circumstellar envelope by orbital motion of the AGB star. The column density pattern only depends on two parameters: eccentricity of the orbit e and the terminal outflow velocity to mean orbital velocity ratio gamma. When viewed perpendicular to the orbital plane, spiral, broken spiral, and incomplete concentric shell patterns can be seen, while when viewed along the orbital plane, alternative concentric half-shell, egg-shell, and half-shell half-gap patterns will develop. Non-zero eccentricity causes asymmetry, while larger gamma makes a weaker pattern and helps bring out asymmetry. A spiral pattern may becomes broken when e > 0.4. The spiral center is always less than 12% of spiral pitch away from the orbit center. One should have more chances (~ 80%) seeing spiral-like patterns than seeing concentric shells (~ 20%) in the circumstellar envelope of wide binary AGB stars.

In order to explore the statistical properties of galaxy morphological types in compact groups (CGs), we construct a random group sample which has the same distributions of redshift and number of member galaxies as those of the CG sample. It turns out that the proportion of early-type galaxies in different redshift bins for the CG sample is statistically higher than that for random group sample, and with growing redshift z this kind of difference becomes more significant. This may be due to the existence of interactions and mergers within a significant fraction of SDSS CGs. We also compare statistical results of CGs with those of more compact groups and pairs, but do not observe as large statistical difference as Hickson (1982)'results.

Using a volume-limited sample of Main Galaxies from the SDSS Data Release 5, we investigate the dependence of galaxy properties on local environment. For each galaxy, the local three-dimensional density is calculated. We find that galaxy morphologies strongly depend on local environment: galaxies in dense environments have predominantly early type morphologies, but other galaxy properties do not present significant dependence on local environment. Clearly, this puts a important constraint on proposed physical mechanisms.

S\'ersic (1968) generalized the de Vaucouleurs law which follows the projected (observed) one dimensional radial profile of elliptical galaxies closely and Dehnen (1993) proposed an analytical formula of the 3-dimensional light distributions whose projected line profile resembles the de Vaucouleurs law. This paper is involved to recover the Dehnen model and generalize the model to account for galaxy elliptical shapes by means of curvilinear coordinate systems and employing a symmetry principle. The symmetry principle maps an orthogonal coordinate system to a light distribution pattern. The coordinate system for elliptical galaxy patterns turns out to be the one which is formed by the complex-plane reciprocal transformation $Z=1/W$. The resulting spatial (3-dimensional) light distribution is spherically symmetric and has infinite gradient at its centre, which is called spherical-nucleus solution and is used to model galaxy central area. We can make changes of the coordinate system by cutting out some column areas of its definition domain, the areas containing the galaxy centre. The resulting spatial (3-dimensional) light distributions are axisymmetric or triaxial and have zero gradient at the centre, which are called elliptical-shape solutions and are used to model global elliptical patterns. The two types of logarithmic light distributions are added together to model full elliptical galaxy patterns. The model is a generalization of the Dehnen model. One of the elliptical-shape solutions permits realistic numerical calculation and is fitted to all R-band elliptical images from the Frei {\it et al.}(1996)'s galaxy sample. The fitting is satisfactory. This suggests that elliptical galaxy patterns can be represented in terms of a few basic parameters.

The observed SEDs of a sample of 60 OH/IR stars are fitted using a radiative transfer model of a dusty envelope. Among the whole sample, 21 stars have reliable phase-lag distances while the others have less accurate distances. L*-P,Mlr-P and Mlr-L* relations have been plotted for these stars. It is found that type E (with emission feature at 10um and type A (with absorption feature at 10um) OH/IR stars have different L*-P and Mlr-L* relations while both of them follow a single Mlr-P relation. The type E stars are proven to be located in the area without large scale dense interstellar medium while the type A stars are located probably in dense interstellar medium. It is argued here that this may indicate the two types of OH/IR stars have different chemical composition or zero age main sequence mass and so evolve in different ways. This conclusion has reinforced the argument by Chen et al.(2001) who reached a similar conclusion from the galactic distribution of about 1000 OH/IR stars with the IRAS low-resolution spectra (LRS).

(Pseudo) radiative pumprate of OH 1612 MHz masers is defined for a sample of 44 OH/IR sources (infrared sources with OH 1612 MHz maser), irrespective of the real maser pumping mechanisms. The correlation between the (pseudo) maser pumprates and the evolutionary status of the maser sources reveals that the radiative pumprates of stellar OH masers are nearly fixed, which agrees with the theoretical prediction for radiatively pumped OH maser. The (pseudo) radiative pumprates of interstellar OH masers are not only very small but also varying broadly over two orders of magnitude, which is argued to be the manifestation of varying number of quiet absorbing OH cloudlets and/or various OH maser pumping mechanisms and/or competitive gain between mainline and 1612 MHz OH masers and/or anisotropy of the maser emission. The radiative pumprates of post-AGB OH masers very possibly decrease with increasing IRAS C32 color indices and distribute in an interim region between the stellar and interstellar OH masers in the pumprate-color diagram.

The 1612 MHz hydroxyl maser in circumstellar envelopes has long been thought to be pumped by 34.6um photons. Only recently, the Infrared Space Observatory has made possible spectroscopic observations which enable the direct confirmation of this pumping mechanism in a few cases. To look for the presence of this pumping line, we have searched the Infrared Space Observatory Data Archive and found 178 spectra with data around 34.6um for 87 galactic 1612MHz masers. The analysis performed showed that the noise level and the spectral resolution of the spectra are the most important factors affecting the detection of the 34.6um absorption line. Only 5 objects from the sample (3 red supergiants and 2 galactic center sources) are found to show clear 34.6um absorption (all of them already known) while two additional objects only tentatively show this line. The 3 supergiants show similar pump rates and their masers might be purely radiatively pumped. The pump rates of OH masers in late type stars are found to be about 0.05, only 1/5 of the theoretical value of 0.25 derived by Elitzur (1992). We have also found 16 maser sources which, according to the analysis assuming Elitzur's pump rate, should show the 34.6 $\mu$m absorption line but do not. These non-detections can be tentatively explained by far-infrared photon pumping, clumpy nature of the OH masing region or a limb-filling emission effect in the OH shell.

This paper analyzes the 53.3um line in the ISO LWS spectra towards a similar sample of OH/IR sources. We find 137 LWS spectra covering 53.3um and associated with 47 galactic OH/IR sources. Ten of these galactic OH/IR sources are found to show and another 5 ones tentatively show the 53.3um absorption while another 7 sources highly probably do not show this line. The source class is found to be correlated with the type of spectral profile: red supergiants (RSGs) and AGB stars tend to show strong blue-shifted filling emission in their 53.3um absorption line profiles while HII regions tend to show a weak red-shifted filling emission in the line profile. GC sources and megamasers do not show filling emission feature. It is argued that the filling emission might be the manifestation of an unresolved half emission half absorption profile of the 53.3um doublet. The 53.3 to 34.6um equivalent width (EW) ratio is close to unity for RSGs but much larger than unity for GC sources and megamasers while H II regions only show the 53.3um line. The pump rate defined as maser to IR photon flux ratio is approximately 5% for RSGs. The pump rates of GC sources are three order of magnitude smaller. Both the large 53.3 to 34.6um EW ratio and the small pump rate of the GC OH masers reflect that the two detected `pumping lines' in these sources are actually of interstellar origin. The pump rate of Arp 220 is 32%--much larger than that of RSGs, which indicates that the contribution of other pumping mechanisms to this megamaser is important.

The population inversion which leads to the 1612 MHz OH maser emission has long been thought to be radiatively pumped. Since OH rotational lines involved in this pumping scheme lie in the far-infrared they became observable only after the launch of the ISO satellite. With the aim to investigate the pumping conditions of the 1612 MHz OH maser emission in more details we have searched the ISO Archive for SWS observations around 34.6 um of 1024 OH/IR sources with IRAS LRS spectra from compilation of Chen et al. (2001). Surprisingly, among 81 OH/IR sources which have appropriate SWS data only already reported objects: VY CMa, IRC+10420 and the Galactic center, show clear 34.6 um absorption line. We discuss possible reasons for non-detection of this pumping line.