Coronal Mass Ejections (CMEs) often show different features in different band-passes. By combining data in white-light (WL) and ultraviolet (UV) bands, we have applied different techniques to derive plasma temperatures, electron density, internal radial speed, etc, within a fast CME. They serve as extensive tests of the diagnostic capabilities, developed for the observations provided by future multi-channel coronagraphs (such as Solar Orbiter/Metis, ASO-S/LST, PROBA-3/ASPIICS). The involved data include WL images acquired by SOHO/LASCO coronagraphs, and intensities measured by SOHO/UVCS at 2.45 R$_{\odot}$ in the UV (H I Ly$\alpha$ and O VI 1032 {AA} lines) and WL channels. Data from the UVCS WL channel have been employed for the first time to measure the CME position angle with polarization-ratio technique. Plasma electron and effective temperatures of the CME core and void are estimated by combining UV and WL data. Due to the CME expansion and the possible existence of prominence segments, the transit of the CME core results in decreases of the electron temperature down to $10^{5}$ K. The front is observed as a significant dimming in the Ly$\alpha$ intensity, associated with a line broadening due to plasma heating and flows along the line-of-sight. The 2D distribution of plasma speeds within the CME body is reconstructed from LASCO images and employed to constrain the Doppler dimming of Ly$\alpha$ line, and simulate future CME observations by Metis and LST.

Magnetic reconnection in the relativistic and trans-relativistic regimes is able to accelerate particles to hard power law energy spectra $f \propto \gamma^{-p}$ (approaching $p=1$). The underlying acceleration mechanism that determines the spectral shape is currently a topic of intense investigation. By means of fully kinetic plasma simulations, we carry out a study of particle acceleration during magnetic reconnection in the trans-relativistic regime of a proton-electron plasma. While earlier work in this parameter regime has focused on the effects of electric field parallel to the local magnetic field on the particle injection (from thermal energy to the lower energy bound of the power-law spectrum), here we examine the roles of both parallel and perpendicular electric fields to gain a more complete understanding on the injection process and further development of a power-law spectrum. We show that the parallel electric field does contribute significantly to particle injection, and is more important in the initial phase of magnetic reconnection. However, as the simulation proceeds, the acceleration by the perpendicular electric field becomes more important for particle injection and completely dominates the acceleration responsible for the high-energy power-law spectrum. This holds robustly, in particular for longer reconnection times and larger systems, i.e. in simulations that are more indicative of the processes in astrophysical sources.

We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

Magnetic reconnection in strongly magnetized astrophysical plasma environments is believed to be the primary process for fast energy release and particle energization. Currently there is strong interest in relativistic magnetic reconnection, in that it may provide a new explanation for high-energy particle acceleration and radiation in strongly magnetized astrophysical systems. We review recent advances in particle acceleration and reconnection physics in the magnetically-dominated regime. More discussion is focused on the physics of particle acceleration, power-law formation as well as the reconnection rate problem. In addition, we provide an outlook for studying reconnection acceleration mechanisms and kinetic physics in the next step.

We determine the inner density profiles of massive galaxy clusters (M$_{200}$ > $5 \times 10^{14}$ M$_{\odot}$) in the Cluster-EAGLE (C-EAGLE) hydrodynamic simulations, and investigate whether the dark matter density profiles can be correctly estimated from a combination of mock stellar kinematical and gravitational lensing data. From fitting mock stellar kinematics and lensing data generated from the simulations, we find that the inner density slopes of both the total and the dark matter mass distributions can be inferred reasonably well. We compare the density slopes of C-EAGLE clusters with those derived by Newman et al. for 7 massive galaxy clusters in the local Universe. We find that the asymptotic best-fit inner slopes of "generalized" NFW (gNFW) profiles, ${\gamma}_{\rm gNFW}$, of the dark matter haloes of the C-EAGLE clusters are significantly steeper than those inferred by Newman et al. However, the mean mass-weighted dark matter density slopes of the simulated clusters are in good agreement with the Newman et al. estimates. We also find that the estimate of ${\gamma}_{\rm gNFW}$ is very sensitive to the constraints from weak lensing measurements in the outer parts of the cluster and a bias can lead to an underestimate of ${\gamma}_{\rm gNFW}$.

We combine asteroseismology, optical high-resolution spectroscopy, and kinematic analysis for 26 halo red giant branch stars in the \textit{Kepler} field in the range of $-2.5<[\mathrm{{Fe}/{H}}]<-0.6$. After applying theoretically motivated corrections to the seismic scaling relations, we obtain an average mass of $0.96\pm 0.02\,\mathrm{M_{\odot}}$ for our sample of halo stars. This maps into an age of $\sim 8\,\mathrm{Gyr}$, significantly younger than independent age estimates of the Milky Way stellar halo. Although our results suggest that further improvements are needed to accurately determine ages of metal-poor stars, asteroseismology still provides robust relative ages. The lack of a mass dispersion among our sample constrains a relative age dispersion to $<18\%$, corresponding to $<2\,\mathrm{Gyr}$. The precise chemical abundances allow us to separate most of the stars with [{Fe}/{H}]$>-1.7$ into two [{Mg}/{Fe}] groups. While [$\alpha$/{Fe}] and [{Eu}/{Mg}] ratios are different between the two subsamples, [$s$/Eu], where $s$ stands for Ba, La, Ce, and Nd, does not show a significant difference. These abundance ratios suggest that the chemical evolution of the low-Mg population is contributed by type~Ia supernovae, but not by low-to-intermediate mass asymptotic giant branch stars, providing a constraint on its star formation timescale as $100\,\mathrm{Myr}<\tau<300\,\mathrm{Myr}$. We also do not detect any significant mass difference between the two [{Mg}/{Fe}] groups, thus suggesting that their formation epochs are not separated by more than 1.5 Gyr.

We investigate a simplified freeze-in dark-matter model in which the dark matter only interacts with the standard-model neutrinos via a light scalar. The extremely small coupling for the freeze-in mechanism is naturally realized in several neutrino-portal scenarios with the secret neutrino interactions. We study possible evolution history of the dark sector: the dark sector would undergo pure freeze-in production if the interactions between the dark-sector particles are negligible, while thermal equilibrium within the dark sector could occur if the reannihilation of the dark matter and the scalar mediator is rapid enough. We investigate the relic abundance in the freeze-in and dark freeze-out regimes, calculate evolution of the dark temperature, and study its phenomenological aspects on BBN and CMB constraints, the indirect-detection signature, as well as the potential to solve the small scale structure problem.

Sample selection is a necessary preparation for weak lensing measurement. It is well-known that selection itself may introduce bias in the measured shear signal. Using image simulation and the Fourier_Quad shear measurement pipeline, we quantify the selection bias in various commonly used selection function (signal-to-noise-ratio, magnitude, etc.). We proposed a new selection function defined in the power spectrum of the galaxy image. This new selection function has low selection bias, and it is particularly convenient for shear measurement pipelines based on Fourier transformation.

GWAC will have been built an integrated FOV of 5,000 $degree^2$ and have already built 1,800 square $degree^2$. The limit magnitude of a 10-second exposure image in the moonless night is 16R. In each observation night, GWAC produces about 0.7TB of raw data, and the data processing pipeline generates millions of single frame alerts. We describe the GWAC Data Processing and Management System (GPMS), including hardware architecture, database, detection-filtering-validation of transient candidates, data archiving, and user interfaces for the check of transient and the monitor of the system. GPMS combines general technology and software in astronomy and computer field, and use some advanced technologies such as deep learning. Practical results show that GPMS can fully meet the scientific data processing requirement of GWAC. It can online accomplish the detection, filtering and validation of millions of transient candidates, and feedback the final results to the astronomer in real-time. During the observation from October of 2018 to December of 2019, we have already found 102 transients.

The Chern-Simons term, through which the cosmic Axion-like field couples to the electromagnetic field, has the effect to rotate CMB polarization directions and to break the CPT symmetry. This rotation will change the CMB power spectra, no matter isotropic or anisotropic the rotation angle is. In this paper we revisit this issue by further considering the correlations between the (anisotropic) rotation angle $\alpha$ and the CMB temperature and (unrotated) $E$ polarization fields. These correlations could be generated in the Axion-like models with nonzero potential under the adiabatic initial condition. We first investigate how these correlations contribute further modifications to the CMB power spectra, then calculate the CMB bispectra for the temperature and rotated polarization fields. These bispectra would vanish if the $T\alpha$ and $E\alpha$ correlations are absent. So, they are useful in searching for CPT violation and the $T\alpha$ and $E\alpha$ correlations arisen in the Axion-like models.

We perform a suite of hydrodynamic simulations to investigate how initial density profiles of giant molecular clouds (GMCs) affect their subsequent evolution. We find that the star formation duration and integrated star formation efficiency of the whole clouds are not sensitive to the choice of different profiles but are mainly controlled by the interplay between gravitational collapse and stellar feedback. Despite this similarity, GMCs with different profiles show dramatically different modes of star formation. For shallower profiles, GMCs first fragment into many self-gravitation cores and form sub-clusters that distributed throughout the entire clouds. These sub-clusters are later assembled ``hierarchically'' to central clusters. In contrast, for steeper profiles, a massive cluster is quickly formed at the center of the cloud and then gradually grows its mass via gas accretion. Consequently, central clusters that emerged from clouds with shallower profiles are less massive and show less rotation than those with the steeper profiles. This is because 1) a significant fraction of mass and angular momentum in shallower profiles is stored in the orbital motion of the sub-clusters that are not able to merge into the central clusters 2) frequent hierarchical mergers in the shallower profiles lead to further losses of mass and angular momentum via violent relaxation and tidal disruption. Encouragingly, the degree of cluster rotations in steeper profiles is consistent with recent observations of young and intermediate-age clusters. We speculate that rotating globular clusters are likely formed via an ``accretion'' mode from centrally-concentrated clouds in the early Universe.

In this work, we investigate the impacts of the future gravitational-wave (GW) standard siren observation by the Einstein Telescope (ET) on constraining the interacting dark energy (IDE) models. We simulate 1000 GW events in the redshift range of $0\lesssim z \lesssim 5$ based on the 10-year observation of the ET. We combine the simulated GW data with the current mainstream cosmological electromagnetic observations including the cosmic microwave background anisotropies, the baryon acoustic oscillations, and the type Ia supernovae to constrain the IDE models. We consider typical IDE models in the context of a perturbed universe. To avoid the large-scale instability problem for IDE models, we apply the extended parameterized post-Friedmann approach to calculate the cosmological perturbations. We find that the addition of the GW standard siren data could significantly improve the constraint accuracies for most of the cosmological parameters (e.g., $H_{0}$, $w$, and $\Omega_{\rm m}$). For the coupling parameter $\beta$, the constraint errors could also be slightly improved when adding the GW data in the cosmological fit.

In this work, we investigate the constraints on the total neutrino mass in the scenario of vacuum energy interacting with cold dark matter (abbreviated as I$\Lambda$CDM) by using the latest cosmological observations. We consider four typical interaction forms, i.e., $Q=\beta H \rho_{\rm de}$, $Q=\beta H \rho_{\rm c}$, $Q=\beta H_{0} \rho_{\rm de}$, and $Q=\beta H_{0} \rho_{\rm c}$, in the I$\Lambda$CDM scenario. To avoid the large-scale instability problem in interacting dark energy models, we employ the extended parameterized post-Friedmann method for interacting dark energy to calculate the perturbation evolution of dark energy in these models. The observational data used in this work include the cosmic microwave background (CMB) measurements from the Planck 2018 data release, the baryon acoustic oscillation (BAO) data, the type Ia supernovae (SN) observation (Pantheon compilation), and the 2019 local distance ladder measurement of the Hubble constant $H_{0}$ from the Hubble Space Telescope. We find that, compared with those in the $\Lambda$CDM+$\sum m_{\nu}$ model, the constrains on $\sum m_{\nu}$ are looser in the four I$\Lambda$CDM+$\sum m_{\nu}$ models. When considering the three mass hierarchies of neutrinos, the constraints on $\sum m_{\nu}$ are tightest in the degenerate hierarchy case and loosest in the inverted hierarchy case. In addition, in the four I$\Lambda$CDM+$\sum m_{\nu}$ models, the values of coupling parameter $\beta$ are larger using the CMB+BAO+SN+$H_{0}$ data combination than that using the CMB+BAO+SN data combination, and $\beta>0$ is favored at more than 1$\sigma$ level when using CMB+BAO+SN+$H_{0}$ data combination. The issue of the $H_{0}$ tension is also discussed in this paper. We find that, compared with the $\Lambda$CDM+$\sum m_{\nu}$ model, the $H_{0}$ tension can be alleviated in the I$\Lambda$CDM+$\sum m_{\nu}$ model to some extent.

The absolute distance scale measurements in cosmology have always been an important mission. In particular, in recent years, the Hubble constant tension between the measurements from early and late universe has become a new crisis for cosmology, which calls for new, independent absolute cosmological distance measurements. Recently, a result of measuring the parallax distance to 3C 273 through spectroastrometry and reverberation mapping was reported. We comment on this novel method in this News & Views paper.

The omnipresence of super-Earths suggests that they are able to be retained in natal disks around low-mass stars, whereas exoplanets' mass distribution indicates that some cores have transformed into gas giants through runaway gas accretion at 1AU from solar-type stars. In this paper, we show that transition to runaway gas accretion by cores may be self-impeded by an increase of the grain opacity in their envelope after they have acquired sufficient mass (typically 10Mearth) to enter a pebble-isolation phase. The accumulation of mm-m size pebbles in their migration barriers enhances their local fragmentation rates. The freshly produced sub-mm grains pass through the barrier, elevate the effective dust opacity and reduce the radiative flux in the cores envelope. These effects alone are adequate to suppress the transition to runaway accretion and preserve super-Earths in the stellar proximity (0.1 AU), albeit entropy advection between the envelope and the disk can further reduce the accretion rate. At intermediate distance (1AU) from their host stars, the escalation in the dust opacity dominates over entropy advection in stalling the transition to runaway accretion for marginally pebble-isolated cores. Beyond a few AU, the transformation of more massive cores to gas giants is reachable before severe depletion of disk gas. This requirement can be satisfied either in extended disks with large scale height via orderly accretion of migrating pebbles or through the mergers of oligarchic protoplanetary embryos, and can account for the correlated occurrence of long-period gas giants and close-in super-Earths.

Since September 2018, LAMOST starts a new 5-year medium-resolution spectroscopic survey (MRS) using bright/gray nights. We present the scientific goals of LAMOST-MRS and propose a near optimistic strategy of the survey. A complete footprint is also provided. Not only the regular medium-resolution survey, but also a time-domain spectroscopic survey is being conducted since 2018 and will be end in 2023. According to the detailed survey plan, we expect that LAMOST-MRS can observe about 2 million stellar spectra with ~7500 and limiting magnitude of around G=15 mag. Moreover, it will also provide about 200 thousand stars with averagely 60-epoch observations and limiting magnitude of G~14 mag. These high quality spectra will give around 20 elemental abundances, rotational velocities, emission line profiles as well as precise radial velocity with uncertainty less than 1 km/s. With these data, we expect that LAMOST can effectively leverage sciences on stellar physics, e.g. exotic binary stars, detailed observation of many types of variable stars etc., planet host stars, emission nebulae, open clusters, young pre-main-sequence stars etc.

We investigate the impacts of dark energy on constraining massive (active/sterile) neutrinos in interacting dark energy (IDE) models by using the current observations. We employ two typical IDE models, the interacting $w$ cold dark matter (I$w$CDM) model and the interacting holographic dark energy (IHDE) model, to make an analysis. To avoid large-scale instability, we use the parameterized post-Friedmann approach to calculate the cosmological perturbations in the IDE models. The cosmological observational data used in this work include the Planck cosmic microwave background (CMB) anisotropies data, the baryon acoustic oscillation data, the type Ia supernovae data, the direct measurement of the Hubble constant, the weak lensing data, the redshift-space distortion data, and the CMB lensing data. We find that the dark energy properties could influence the constraint limits of active neutrino mass and sterile neutrino parameters in the IDE models. We also find that the dark energy properties could influence the constraints on the coupling strength parameter $\beta$, and a positive coupling constant, $\beta>0$, can be detected at the $2.5\sigma$ statistical significance for the IHDE+$\nu_s$ model by using the all-data combination. In addition, we also discuss the "Hubble tension" issue in these scenarios. We find that the $H_0$ tension can be effectively relieved by considering massive sterile neutrinos, and in particular in the IHDE+$\nu_s$ model the $H_0$ tension can be reduced to be at the $1.28\sigma$ level.

We discuss the merger process of binary black holes with Hawking radiation taken into account. Besides the redshifted radiation to infinity, binary black holes can exchange radiation between themselves, which is first redshifted and then blueshifted when it propagates from one hole to the other. The exchange rate should be large when the temperature-divergent horizons are penetrating each other to form a single horizon with unique temperature. This will cause non-negligible mass and angular momentum transfer between the black holes during the merging process of the horizons. We further argue in the large mass ratio limit that the light hole whose evaporation is enhanced by the competing redshift-blueshift effects will probably evaporate or decay completely before reaching the the horizon of the heavy one, generalising our previous Rindler/Tachyon speculations which state that a probe particle collapsing to a non-extreme black hole will completely decay into gravitons (or closed strings more generally) before reaching the horizon. We also discuss the possibility of testing Hawking radiation and even exploring the information loss puzzle in gravitational wave observations.

Before massive stars heavier than $(8 \cdots 10)$ solar masses evolve to the phase of a gravitational core collapse, they will emit a huge number of MeV-energy neutrinos that are mainly produced in the thermal processes and nuclear weak interactions. The detection of such pre-supernova (pre-SN) neutrinos could provide an important and independent early warning for the optical observations of core-collapse SNe. In this paper, we investigate the prospects of future large liquid-scintillator detectors for the observation of pre-SN neutrinos in both $\bar{\nu}^{}_e + p \to e^+ + n$ and $\nu (\bar{\nu}) + e^- \to \nu (\bar{\nu}) + e^-$ reaction channels, where $\nu$ ($\bar{\nu}$) denotes neutrinos (antineutrinos) of all three flavors. We propose a quantitative assessment of the capability in terms of three working criteria, namely, how far the SN distance can be covered, how long the early warning before the core collapse can be sent out, and how well the direction pointing to the SN can be determined. The dependence of the final results on the different models of progenitor stars, neutrino flavor conversions and the relevant backgrounds is also discussed.

We report the discovery of a highly dispersed fast radio burst, FRB~181123, from an analysis of $\sim$1500~hr of drift-scan survey data taken using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulse has three distinct emission components, which vary with frequency across our 1.0--1.5~GHz observing band. We measure the peak flux density to be $>0.065$~Jy and the corresponding fluence $>0.2$~Jy~ms. Based on the observed dispersion measure of 1812~cm$^{-3}$~pc, we infer a redshift of $\sim 1.9$. From this, we estimate the peak luminosity and isotropic energy to be $\lesssim 2\times10^{43}$~erg~s$^{-1}$ and $\lesssim 2\times10^{40}$~erg, respectively. With only one FRB from the survey detected so far, our constraints on the event rate are limited. We derive a 95\% confidence lower limit for the event rate of 900 FRBs per day for FRBs with fluences $>0.025$~Jy~ms. We performed follow-up observations of the source with FAST for four hours and have not found a repeated burst. We discuss the implications of this discovery for our understanding of the physical mechanisms of FRBs.

The WT scheme, a piecewise polynomial force interpolation scheme with time-step dependency, is proposed in this paper for relativistic particle-in-cell (PIC) simulations. The WT scheme removes the lowest order numerical Cherenkov instability (NCI) growth rate for arbitrary time steps allowed by the Courant condition. While NCI from higher order resonances is still present, the numerical tests show that for smaller time steps, the numerical instability grows much slower than using the optimal time step found in previous studies. The WT scheme is efficient for improving the quality and flexibility of relativistic PIC simulations.

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

We explore the bound of the trans-Planckian censorship conjecture on an inflation model with multiple stages. We show that if the first inflationary stage is responsible for the primordial perturbations in the cosmic microwave background window, the $e$-folding number of each subsequent stage will be bounded by the energy scale of the first stage. This seems to imply that the lifetime of the current era of accelerated expansion (regarded as one of the multiple inflationary stages) might be a probe for distinguishing inflation from its alternatives. We also present a multistage inflation model in a landscape consisting of anti-de Sitter vacua separated by potential barriers.

We have assembled the historical light curves of the BL Lac Object 3FGL J0449.4-4350 at optical and gamma-ray bands, the time spanning about 10 years, analyzed the periodic variability of the light curves by using four different methods (Lomb-Scargle periodogram, REDFIT38, Jurkevich and DACF). We detected a marginally possible quasi-periodic oscillation (QPO) of ~450 days. Assuming it originates from the helical motion jet in a supermassive binary black hole (SMBBH) system undergoing major merger, we estimate the primary black hole mass M~7.7*10^{9} M_sun. To explore the origin of the gamma-ray, we investigated the optical-gamma-ray correlations using discrete correlation function (DCF) method, and found that the correlation between the two bands is very significant. This strong correlation tends to imply lepton self-synchro-Compton (LSSC) model to produce the gamma-ray.

In this work, we study how the dust coagulation/fragmentation will influence the evolution and observational appearances of vortices induced by a massive planet embedded in a low viscosity disk by performing global 2D high-resolution hydrodynamical simulations. Within the vortex, due to its higher gas surface density and steeper pressure gradients, dust coagulation, fragmentation and drift (to the vortex center) are all quite efficient, producing dust particles ranging from micron to $\sim 1.0\ {\rm cm}$, as well as overall high dust-to-gas ratio (above unity). In addition, the dust size distribution is quite non-uniform inside the vortex, with the mass weighted average dust size at the vortex center ($\sim 4.0$ mm) being a factor of $\sim10$ larger than other vortex regions. Both large ($\sim$ mm) and small (tens of micron) particles contribute strongly to affect the gas motion within the vortex. As such, we find that the inclusion of dust coagulation has a significant impact on the vortex lifetime and the typical vortex lifetime is about 1000 orbits. After the initial gaseous vortex is destroyed, the dust spreads into a ring with a few remaining smaller gaseous vortices with a high dust concentration and a large maximum size ($\sim$ mm). At late time, the synthetic dust continuum images for the coagulation case show as a ring inlaid with several hot spots at 1.33 mm band, while only distinct hot spots remain at 7.0 mm.

One of the most unique phenomena in the Galactic center region is the existence of numerous long and narrow filamentary structures within a few hundred parsecs of Sgr A$^{\star}$. While more than one than one hundred radio filaments have been revealed by MeerKAT, about two dozens X-ray filaments have been discovered so far. In this article, we report our analysis on the deep Chandra and NuSTAR observations of a non-thermal X-ray filament, G0.13-0.11, which is located adjacent to the Radio arc. Chandra revealed a unique morphology of G0.13-0.11, which is an elongated (0.1 pc in width and 3.2 pc in length) structure slightly bended towards the Radio arc. A pulsar candidate ($\Gamma \sim 1.4$) is detected in the middle of the filament, with a tail of diffuse non-thermal X-ray emission on one side of the filament. The filament is detected by NuSTAR up to 79 keV, with the hard X-ray centroid consistent with the pulsar candidate. We found that the X-ray intensity decays along the filament farther away from the pulsar candidate, dropping to half of its peak value at 2.2 pc away. This system is mostly likely a Pulsar Wind Nebula interacting with ambient interstellar magnetic field, where the filaments are kinetic jets from PWN as recently proposed. The nature of this filament adds to complex origin of the X-ray filaments, which serve as powerful tools to probe local and global powerful particle accelerators in the Galactic center.

The nature of Milky Way dwarf spheroidals (MW dSphs) has been questioned, in particular whether they are dominated by dark matter (DM). Here we investigate an alternative scenario, for which tidal shocks are exerted by the MW to DM-free dSphs after a first infall of their gas-rich progenitors, and for which theoretical calculations have been verified by pure N-body simulations. Whether or not the dSphs are on their first infall cannot be resolved on the sole basis of their star formation history. In fact, gas removal may cause complex gravitational instabilities and near-pericenter passages can give rise to tidal disruptive processes. Advanced precision with the Gaia satellite in determining both their past orbital motions and the MW velocity curve is, however, providing crucial results. First, tidal shocks explain why DM-free dSphs are found preferentially near their pericenter, where they are in a destructive process, while their chance to be long-lived satellites is associated with a very low probability P~ 2 10^-7, which is at odds with the current DM-dominated dSph scenario. Second, most dSph binding energies are consistent with a first infall. Third, the MW tidal shocks that predict the observed dSph velocity dispersions are themselves predicted in amplitude by the most accurate MW velocity curve. Fourth, tidal shocks accurately predict the forces or accelerations exerted at half-light radius of dSphs, including the MW and the Magellanic System gravitational attractions. The above is suggestive of dSphs that are DM-free and tidally shocked near their pericenters, which may provoke a significant quake in our understanding of near-field cosmology.

We investigated the quasi-periodic pulsation (QPP) in Lyman-alpha, X-ray and extreme-ultraviolet (EUV) emissions during two solar flares, i.e., an X-class (SOL2012-01-27T) and a C-class (SOL2016-02-08T). The full-disk Lyman-alpha and X-Ray flux during these solar flares were recorded by the EUV Sensor and X-Ray Sensor on board the Geostationary Operational Environmental Satellite. The {\deg}are regions were located from the EUV images measured by the Atmospheric Imaging Assembly. The QPP could be identified as a series of regular and periodic peaks in the light curves, and its quasi-periodicity was determined from the global wavelet and Fourier power spectra. A quasi-periodicity at about 3 minutes is detected during the impulsive phase of the X-class flare, which could be explained as the acoustic wave in the chromosphere (e.g., Milligan et al. 2017). Interestingly, a quasi-periodicity at roughly 1 minute is discovered during the entire evolutionary phases of solar flares, including the precursor, impulsive, and gradual phases. This is the first report of 1-minute QPP in the Lyman-alpha emission during solar flares, in particular during the flare precursor. It may be interpreted as a self-oscillatory regime of the magnetic reconnection, such as magnetic dripping.

We investigate the axion like particle (ALP)-photon oscillation effect in the high energy $\gamma$-ray spectra of PG 1553+113 and PKS 2155-304 measured by Fermi-LAT and H.E.S.S.. The choice of extragalactic background light (EBL) model, which induces the attenuate effect in observed $\gamma$-ray spectra, would affect the ALP implication. For the ordinary EBL model that prefers a null hypothesis, we set constraint on the ALP-photon coupling constant at 95% C.L. as $g_{a\gamma}\lesssim 5\times 10^{-11} ~\rm{GeV}^{-1}$ for the ALP mass $\sim 10$ neV. We also consider the CIBER observation of the cosmic infrared radiation, which shows an excess at the wave wavelength of $\sim 1~\mu$m after the substraction of foregrounds. The high energy gamma-rays from extragalactic sources at high redshifts would suffer from a more significant attenuate effect caused by this excess. In this case, we find that the ALP-photon oscillation would improve the fit to the observed spectra of PKS 2155-304 and PG 1553+113 and find a favored parameter region at 95% C.L..

The Five-hundred-meter Aperture Spherical radio Telescope (FAST) passed national acceptance and is taking pilot cycle of 'Shared-Risk' observations. The 19-beam receiver covering 1.05-1.45 GHz was used for most of these observations. The electronics gain fluctuation of the system is better than 1\% over 3.5 hours, enabling enough stability for observations. Pointing accuracy, aperture efficiency and system temperature are three key parameters of FAST. The measured standard deviation of pointing accuracy is 7.9$''$, which satisfies the initial design of FAST. When zenith angle is less than 26.4$^\circ$, the aperture efficiency and system temperature around 1.4 GHz are $\sim$ 0.63 and less than 24 K for central beam, respectively. The measured value of these two parameters are better than designed value of 0.6 and 25 K, respectively. The sensitivity and stability of the 19-beam backend are confirmed to satisfy expectation by spectral HI observations toward N672 and polarization observations toward 3C286. The performance allows FAST to take sensitive observations in various scientific goals, from studies of pulsar to galaxy evolution.

Recent cosmological hydrodynamical simulations are able to reproduce numerous statistical properties of galaxies that are consistent with observational data. Yet, the adopted subgrid models strongly affect the simulation outcomes, limiting the predictive power of these simulations. In this work, we perform a suite of isolated galactic disk simulations under the {\it SMUGGLE} framework and investigate how different subgrid models affect the properties of giant molecular clouds (GMCs). We employ {\sc astrodendro}, a hierarchical clump-finding algorithm, to identify GMCs in the simulations. We find that different choices of subgrid star formation efficiency, $\epsilon_{\rm ff}$, and stellar feedback channels, yield dramatically different mass and spatial distributions for the GMC populations. Without feedback, the mass function of GMCs has a shallower power-law slope and extends to higher mass ranges compared to runs with feedback. Moreover, higher $\epsilon_{\rm ff}$ results in faster molecular gas consumption and steeper mass function slopes. Feedback also suppresses power in the two-point correlation function (TPCF) of the spatial distribution of GMCs. Specifically, radiative feedback strongly reduces the TPCF on scales below 0.2~kpc, while supernova feedback reduces power on scales above 0.2~kpc. Finally, runs with higher $\epsilon_{\rm ff}$ exhibit a higher TPCF than runs with lower $\epsilon_{\rm ff}$, because the dense gas is depleted more efficiently thereby facilitating the formation of well-structured supernova bubbles. We argue that comparing simulated and observed GMC populations can help better constrain subgrid models in the next-generation of galaxy formation simulations.

High spatial resolution observations of protoplanetary disks (PPDs) by ALMA have revealed many details that are providing interesting constraints on the disk physics as well as dust dynamics, both of which are essential for understanding planet formation. We carry out high-resolution, 2D global hydrodynamic simulations, including the effects of dust feedback, to study the stability of dusty rings. When the ring edges are relatively sharp and the dust surface density becomes comparable to the gas surface density, we find that dust feedback enhances the radial gradients of both the azimuthal velocity profile and the potential vorticity profile at the ring edges. This eventually leads to instabilities on meso-scales (spatial scales of several disk scale heights), causing dusty rings to be populated with many compact regions with highly concentrated dust densities on meso-scales. We also produce synthetic dust emission images using our simulation results and discuss the comparison between simulations and observations.

The observation of the transient sky through a multitude of astrophysical messengers hasled to several scientific breakthroughs these last two decades thanks to the fast evolution ofthe observational techniques and strategies employed by the astronomers. Now, it requiresto be able to coordinate multi-wavelength and multi-messenger follow-up campaign withinstruments both in space and on ground jointly capable of scanning a large fraction of thesky with a high imaging cadency and duty cycle. In the optical domain, the key challengeof the wide field of view telescopes covering tens to hundreds of square degrees is to dealwith the detection, the identification and the classification of hundreds to thousands of opticaltransient (OT) candidates every night in a reasonable amount of time. In the last decade, newautomated tools based on machine learning approaches have been developed to perform thosetasks with a low computing time and a high classification efficiency. In this paper, we presentan efficient classification method using Convolutional Neural Networks (CNN) to discard anybogus falsely detected in astrophysical images in the optical domain. We designed this toolto improve the performances of the OT detection pipeline of the Ground Wide field AngleCameras (GWAC) telescopes, a network of robotic telescopes aiming at monitoring the opticaltransient sky down to R=16 with a 15 seconds imaging cadency. We applied our trainedCNN classifier on a sample of 1472 GWAC OT candidates detected by the real-time detectionpipeline. It yields a good classification performance with 94% of well classified event and afalse positive rate of 4%.

We explore the isothermal total density profiles of early-type galaxies (ETGs) in the IllustrisTNG simulation. For the selected 559 ETGs at $z = 0$ with stellar mass $10^{10.7}\mathrm{M}_{\odot} \leqslant M_{\ast} \leqslant 10^{11.9}\mathrm{M}_{\odot}$, the total power-law slope has a mean of $\langle\gamma^{\prime}\rangle = 2.011 \pm 0.007$ and a scatter of $\sigma_{\gamma^{\prime}} = 0.171$ over the radial range 0.4 to 4 times the stellar half mass radius. Several correlations between $\gamma^{\prime}$ and galactic properties including stellar mass, effective radius, stellar surface density, central velocity dispersion, central dark matter fraction and in-situ-formed stellar mass ratio are compared to observations and other simulations, revealing that IllustrisTNG reproduces many correlation trends, and in particular, $\gamma^{\prime}$ is almost constant with redshift below $z = 2$. Through analyzing IllustrisTNG model variations we show that black hole kinetic winds are crucial to lowering $\gamma^{\prime}$ and matching observed galaxy correlations. The effects of stellar winds on $\gamma^{\prime}$ are subdominant compared to AGN feedback, and differ due to the presence of AGN feedback from previous works. The density profiles of the ETG dark matter halos are well-described by steeper-than-NFW profiles, and they are steeper in the full physics (FP) run than their counterparts in the dark matter only (DMO) run. Their inner density slopes anti-correlates (remain constant) with the halo mass in the FP (DMO) run, and anti-correlates with the halo concentration parameter $c_{200}$ in both types of runs. The dark matter halos of low-mass ETGs are contracted whereas high-mass ETGs are expanded, suggesting that variations in the total density profile occur through the different halo responses to baryons.

Magnetic reconnection, topological change in magnetic fields, is a fundamental process in magnetized plasmas. It is associated with energy release in regions of magnetic field annihilation, but this is only one facet of this process. Astrophysical flows normally have very large Reynolds numbers and are expected to be turbulent, in agreement with observations. In strong turbulence magnetic lines constantly reconnect everywhere at all scales, making magnetic reconnection an intrinsic part of turbulent cascade. We note that this is inconsistent with the usual practice of regarding magnetic lines as persistent dynamical elements. A number of theoretical, numerical, and observational studies, starting with Lazarian & Vishniac (1999), demonstrated that 3D turbulence makes magnetic reconnection fast and that these two processes are intrinsically connected. We discuss the dramatic violation of the textbook concept of magnetic flux-freezing in the presence of turbulence and demonstrate that in the presence of turbulence the plasma effects are subdominant to turbulence as far as the magnetic reconnection is concerned. This justifies an MHD-like treatment of magnetic reconnection at scales much larger than the relevant plasma scales. We discuss numerical and observational evidences supporting the turbulent reconnection model. In particular, we show that tearing reconnection is suppressed in 3D and, unlike the 2D case, the 3D reconnection induces turbulence that makes reconnection independent of resistivity. We show that turbulent reconnection dramatically affects the key astrophysical processes, e.g., star formation, turbulent dynamo, acceleration of cosmic rays. We provide criticism of the concept of "reconnection-mediated turbulence" and explain why turbulent reconnection is very different from enhanced turbulent resistivity and hyper-resistivity, and why the latter has fatal conceptual flaws.

Cosmological time crystal (TC) corresponds to a matter state where the periodic motion of field forms a limit cycle in its phase space. We explore what would happen if it existed in inflationary phase. It is found that the limit cycle responsible for TC will inevitably cause the periodic oscillation of the primordial perturbation spectrum. The oscillatory patterns of the spectrum depend on the TC parameters, and so encode the crystalline patterns of TC.

Tidal interactions between the embedded planets and their surrounding protoplanetary disks are often postulated to produce the observed complex dust substructures, including rings, gaps, and asymmetries. In this Letter, we explore the consequences of dust coagulation on the dust dynamics and ring morphology. Coagulation of dust grains leads to dust size growth which, under typical disk conditions, produces faster radial drifts, potentially threatening the dust ring formation. Utilizing 2D hydrodynamical simulations of protoplanetary disks which include a full treatment of dust coagulation, we find that if the planet does not open a gap quickly enough, the formation of an inner ring is impeded due to dust coagulation and subsequent radial drift. Furthermore, we find that a "buildup" of sub-mm sized grains often appears in the dust emission at the outer edge of the dust disk.

Many studies have considered the roles of clouds, outflows, and turbulence in producing the broad emission and absorption lines in the spectra of active galactic nuclei (AGNs). However, these are often treated as separate or even competing models. Here, we consider the possibility that AGN clouds are condensations formed within the thermally unstable zones of outflows and then compare the typical sizes of such condensations with the injection scale $k_0^{-1} \sim L_0$ of turbulence, where $L_0$ is assumed to be the scale height of a representative global outflow model. We find that for broad line region (BLR) parameters, clouds are many orders of magnitude smaller than $L_0$ and this has the following implication: BLR cloud dynamics can be modeled using a local approximation through the use of multiphase turbulence simulations of X-ray irradiated plasmas. We present the first such 3D local clumpy turbulent outflow simulations. We show that the condensations share the same type of selection effects characterizing the locally optimally emitting cloud (LOC) scenario, thereby offering a physical interpretation for the LOC model and accounting for its almost uncanny successes. The ubiquitous presence of emission line regions in AGNs can be simply explained as the natural outcome of there being a multiphase interval of $k$-space within the inertial range of a turbulent cascade.

We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code AREPO-RT, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically-symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in AREPO-RT, which treats UV and optical photons as single-scattered and IR photons as multi-scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

We present a pilot HI survey of 17 Planck Galactic Cold Clumps (PGCCs) with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). HI Narrow Self-Absorption (HINSA) is an effective method to detect cold HI being mixed with molecular hydrogen H$_2$ and improves our understanding of the atomic to molecular transition in the interstellar medium. HINSA was found in 58\% PGCCs that we observed. The column density of HINSA was found to have an intermediate correlation with that of $^{13}$CO, following $\rm log( N(HINSA)) = (0.52\pm 0.26) log(N_{^{13}CO}) + (10 \pm 4.1) $. HI abundance relative to total hydrogen [HI]/[H] has an average value of $4.4\times 10^{-3}$, which is about 2.8 times of the average value of previous HINSA surveys toward molecular clouds. For clouds with total column density N$\rm_H >5 \times 10^{20}$ cm$^{-2}$, an inverse correlation between HINSA abundance and total hydrogen column density is found, confirming the depletion of cold HI gas during molecular gas formation in more massive clouds. Nonthermal line width of $^{13}$CO is about 0-0.5 km s$^{-1}$ larger than that of HINSA. One possible explanation of narrower nonthermal width of HINSA is that HINSA region is smaller than that of $^{13}$CO. Based on an analytic model of H$_2$ formation and H$_2$ dissociation by cosmic ray, we found the cloud ages to be within 10$^{6.7}$-10$^{7.0}$ yr for five sources.

We study the evolutionary trend of the total density profile of early-type galaxies (ETGs) in IllustrisTNG. To this end, we trace ETGs from $z=0$ to $z=4$ and measure the power-law slope $\gamma^{\prime}$ of the total density profile for their main progenitors. We find that their $\gamma^{\prime}$ steepen on average during $z\sim4-2$, then becoming shallower until $z=1$, after which they remain almost constant, aside from a residual trend of becoming shallower towards $z=0$. We also compare to a statistical sample of ETGs at different redshifts, selected based on their luminosity profiles and stellar masses. Due to different selection effects, the average slopes of the statistical samples follow a modified evolutionary trend. They monotonically decrease since $z=3$, and after $z\approx 1$, they remain nearly invariant with a mild increase towards $z=0$. These evolutionary trends are mass-dependent for both samples, with low-mass galaxies having in general steeper slopes than their more massive counterparts. Galaxies that transitioned to ETGs more recently have steeper mean slopes as they tend to be smaller and more compact at any given redshift. By analyzing the impact of mergers and AGN feedback on the progenitors' evolution, we conjecture a multi-phase path leading to isothermality in ETGs: dissipation associated with rapid wet mergers tends to steepen $\gamma^{\prime}$ from $z=4$ to $z=2$, whereas subsequent AGN feedback (especially in the kinetic mode) makes $\gamma^{\prime}$ shallower again from $z=2$ to $z=1$. Afterwards, passive evolution from $z=1$ to $z=0$, mainly through gas-poor mergers, mildly decreases $\gamma^{\prime}$ and maintains the overall mass distribution close to isothermal.

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

We have studied the high-mass star-forming complex NGC 6334 with ALMA in the continuum emission at a frequency of 87.6 GHz, achieving a spatial resolution of 1300 au. Detecting 142 compact sources distributed over the whole observed area, we then used machine learning algorithms to group the compact cores in different clusters. A total of four main clusters were identified: NGC 6334-E, NGC 6334-I, NGC 6334-I(N) and NGC 6334-I(NW). The typical separations between cluster members (4000-12000 au) together with the core masses (0.2-100 $M_{\odot}$) are in agreement with turbulent fragmentation at scales of 0.1 pc. We find that the CMFs (core mass functions) show an excess of massive cores compared to the IMF. Typical uncertainties in temperature or unresolved multiplicity may mimic the observed excess of massive cores. Evidence of mass segregation in NGC 6334-I and NGC 6334-I(N), with the most massive cores located closer to the center, was found. Correlations between the physical properties of the clusters and their evolutionary stage were searched for, and found to show a larger separation between cores in the more evolved clusters, favouring the role of gas expulsion and stellar ejection over time. From these results, we suggest that NGC 6334-I(N) was primordially segregated, while NGC 6334-I may have become mass-segregated over time due to dynamical effects. Finally, the lack of massive cores in the most evolved cluster suggests that the gas reservoir may be already exhausted, while the less evolved clusters still have access to a larger gas mass reservoir. In general, the fragmentation in NGC 6334 at large scales (about 1 pc) is governed by turbulent pressure, while at smaller scales (a few hundred au), thermal pressure regulates the fragmentation process.

We report our photometric and spectroscopic follow-up observations of the superoutbursts of three dwarf novae (GWAC\,180415A, GWAC\,181017A and GWAC\,181211A) identified independently by the Ground Wide-angle Cameras system, one of the ground-based instruments of the China-France SVOM mission. Based on a combination of our photometry and that taken from the AAVSO, our period analysis of the superhumps enables us to determine the mass ratios to be 0.0967-0.1163, 0.1879-0.1883 and 0.0981-0.1173 for GWAC\,180415A, GWAC\,181017A and GWAC\,181211A, respectively. GWAC\,180415A can be firmly identified as a WZ sge-type dwarf novae due to its long duration ($\sim2$ weeks) multiple rebrightenings with amplitudes of 3-4 magnitudes, the early superhump associated with a double-wave modulation and the low mass ratio. The inferred low mass ratio and location in the $\varepsilon-P_{\mathrm{orb}}$ diagram suggest that GWAC\,181211A is a WZ sge-type dwarf novae candidate. The measured Balmer decrements suggest the Balmer line emission is produced from an optical thick region in GWAC\,180415A and GWAC\,181017A, and from an optical thin region in GWAC\,181211A.

Previous studies have analyzed the energy injection into the interstellar matter due to molecular bubbles. They found that the total kinetic energies of bubbles are comparable to, or even larger than, those of outflows but still less than the gravitational potential and turbulence energies of the hosting clouds. We examined the possibility that previous studies underestimated the energy injection due to being unable to detect dim or incomplete bubbles. We simulated typical molecular bubbles and inserted them into the $^{13}$CO Five College Radio Astronomical Observatory maps of the Taurus and Perseus Molecular Clouds. We determined bubble identification completeness by applying the same procedures to both simulated and real data sets. We proposed a detectability function for both the Taurus and Perseus molecular clouds based on a multivariate approach. In Taurus, bubbles with kinetic energy less than ~$1 \times 10^{44}$ erg are likely to be missed. We found that the total missing kinetic energy in Taurus is less than a couple of $10^{44}$ erg, which only accounts for around 0.2% of the total kinetic energy of identified bubbles. In Perseus, bubbles with kinetic energy less than ~$2 \times 10^{44}$ erg are likely to be missed. We found that the total missing kinetic energy in Perseus is less than $10^{45}$ erg, which only accounts for around 1% of the total kinetic energy of identified bubbles. We thus conclude that previous manual bubble identification routines used in Taurus and Perseus can be considered to be energetically complete. Therefore, we confirm that energy injection from dynamic structures, namely outflows and bubbles, produced by star formation feedback are sufficient to sustain turbulence at a spatial scale from ~0.1 to ~2.8 pc.

Using high spatial and temporal data from the New Vacuum Solar Telescope (NVST) and the Solar Dynamics Observatory (SDO), we present unambiguous observations of recurrent two-sided loop jets caused by magnetic reconnection between erupting minifilaments and nearby large filament. The observations demonstrate that three two-sided loop jets, which ejected along the large filament in opposite directions, had similar appearance and originated from the same region. We find that a minifilament erupted and drove the first jet. It reformed at the same neutral line later, and then underwent partial and total eruptions, drove the second and third jets, respectively. In the course of the jets, cool plasma was injected into the large filament. Furthermore, persistent magnetic flux cancelation occurred at the neutral line under the minifilament before its eruption and continued until the end of the observation. We infer that magnetic flux cancellation may account for building and then triggering the minifilament to erupt to produce the two-sided loop jets. This observation not only indicates that two-sided loop jets can be driven by minifilament eruptions, but also sheds new light on our understanding of the recurrent mechanism of two-sided loop jets.

The second observational campaign of gravitational waves organized by the LIGO/Virgo Collaborations has led to several breakthroughs such as the detection of gravitational wave signals from merger systems involving black holes or neutrons stars. During O2,14 gravitational wave alerts were sent to the astronomical community with sky regions covering mostly over hundreds of square degrees. Among them, 6 have been finally confirmed as real astrophysical events. Since 2013, a new set of ground-based robotic telescopes called GWAC and its pathfinder mini-GWAC have been developed to contribute to the various challenges of themulti-messenger and time domain astronomy. The GWAC system is built up in the framework of the ground-segment system of the SVOM mission that will be devoted to the study of the multi-wavelength transient sky in the next decade. During O2, only the mini-GWAC telescopenetwork was fully operational. Due to the wide field of view and fast automatic follow-up capabilities of the mini-GWAC telescopes, they were well adapted to efficiently cover the sky localization areas of the gravitational wave event candidates. In this paper, we present the mini-GWAC pipeline we have set up to respond to the GW alerts and we report our optical follow-up observations of 8 GW alerts detected during the O2 run. Our observations provided the largest coverage of the GW localization areas in a short latency made by any optical facility. We found tens of optical transient candidates in our images, but none of those could be securely associated with any confirmed black hole-black hole merger event. Based on this first experience and the near future technical improvements of our network system, we will be more competitive to detect the optical counterparts from some gravitational wave events that will be detected during the upcoming O3 run, especially those emerging from binary neutron star mergers.

The classic Bondi solution remains a common starting point both for studying black hole growth across cosmic time in cosmological simulations and for smaller scale simulations of AGN feedback. In nature, however, there will be inhomogenous distributions of rotational velocity and density along the outer radius ($R_o$) marking the sphere of influence of a black hole. While there have been many studies of how the Bondi solution changes with a prescribed angular momentum boundary condition, they have all assumed a constant density at $R_o$. In this Letter, we show that a non-uniform density at $R_o$ causes a meridional flow and due to conservation of angular momentum, the Bondi solution qualitatively changes into an inflow-outflow solution. Using physical arguments, we analytically identify the critical logarithmic density gradient $|\partial{\ln{\rho}}/\partial{\theta}|$ above which this change of the solution occurs. For realistic $R_o$, this critical gradient is less than 0.01 and tends to 0 as $R_o \rightarrow \infty$. We show using numerical simulations that, unlike for solutions with an imposed rotational velocity, the accretion rate for solutions under an inhomogenous density boundary condition remains constant at nearly the Bondi rate $\dot{M}_B$, while the outflow rate can greatly exceed $\dot{M}_B$.

POLAR is a Compton polarimeter sensitive in the 50 to 500 keV energy range. The Crab pulsar is a scientific target for POLAR on board the Chinese space laboratory Tiangong-2 (TG-2). With its large Field of View (FoV), POLAR detected significant pulsed signals from the Crab pulsar which is visible by POLAR in about half of observation time. In this work, we present the preliminary results including the pulse profile, timing and polarization measuring method. First, we show the highly significant pulse profile observed by POLAR which is compared to the results of other detectors including Fermi/LAT and INTEGRAL. And the pulse profile as a function of theta incident angle and as a function of channel number, which indicate that POLAR has a good detection performance, have been showed. Second, we find that the timing of the Crab pulses are accurately measured, which provides a unique verification and calibration to the POLAR timing system. Finally, the potential polarization measurement of the Crab pulsar is also discussed.

The POLAR detector is a space based Gamma-Ray Burst (GRB) polarimeter sensitive in the 15-500 keV energy range. Apart from its main scientific goal as a Gamma-Ray Burst polarimeter it is also able to detect photons from pulsars in orbit. By using the six-months in-orbit observation data, significant pulsation from the PSR B0531+21 (Crab pulsar) was obtained. In this work, we present the precise timing analysis of the Crab pulsar, together with a phase-resolved spectroscopic study using a joint-fitting method adapted for wide field of view instruments like POLAR. By using single power law fitting over the pulsed phase, we obtained spectral indices ranging from 1.718 to 2.315, and confirmed the spectral evolution in a reverse S shape which is homogenous with results from other missions over broadband. We will also show, based on the POLAR in-orbit performance and Geant4 Monte-Carlo simulation, the inferred capabilities of POLAR-2, the proposed follow-up mission of POLAR on board the China Space Station (CSS), for pulsars studies.

We search for dynamical substructures in the LAMOST DR3 very metal-poor (VMP) star catalog. After cross-matching with Gaia DR2, there are 3300 VMP stars with available high-quality astrometric information that have halo-like kinematics. We apply a method based on self-organizing maps to find groups clustered in the 4D space of orbital energy and angular momentum. We identify 57 dynamically tagged groups, which we label DTG-1 to DTG-57. Most of them belong to existing substructures in the nearby halo, such as the $Gaia$ Sausage or Sequoia. The stream identified by Helmi et al. is recovered, but the two disjoint portions of the substructure have distinct dynamical properties. The very retrograde substructure Rg5 found previously by Myeong et al. is also retrieved. We report 6 new DTGs with highly retrograde orbits, 2 with very prograde orbits, and 12 with polar orbits. By mapping other datasets (APOGEE halo stars, and catalogs of r-process-enhanced and CEMP stars) onto the trained neuron map, we can associate stars with detailed chemical abundances to the DTGs, and look for associations with chemically peculiar stars. The highly eccentric $Gaia$ Sausage groups contain representatives both of debris from the satellite itself (which is $\alpha$-poor) and the Splashed Disk, sent up into eccentric halo orbits from the encounter (and is $\alpha$-rich). The new prograde substructures also appear to be associated with the Splashed Disk. The DTGs belonging to the $Gaia$ Sausage host two relatively metal-rich $r$-II stars and six CEMP stars in different sub-classes, consistent with the idea that the $Gaia$ Sausage progenitor is a massive dwarf galaxy. Rg5 is dynamically associated with two highly $r$-process-enhanced stars with [Fe/H] $\sim -$3. This finding indicates that its progenitor might be an ultra-faint dwarf galaxy that has experienced $r$-process enrichment from neutron star mergers.

We investigate the impact of a highly eccentric 10 $M_{\rm \oplus}$ (where $M_{\rm \oplus}$ is the Earth mass) planet embedded in a dusty protoplanetary disk on the dust dynamics and its observational implications. By carrying out high-resolution 2D gas and dust two-fluid hydrodynamical simulations, we find that the planet's orbit can be circularized at large radii. After the planet's orbit is circularized, partial gap opening and dust ring formation happen close to the planet's circularization radius, which can explain the observed gaps/rings at the outer region of disks. When the disk mass and viscosity become low, we find that an eccentric planet can even open gaps and produce dust rings close to the pericenter and apocenter radii before its circularization. This offers alternative scenarios for explaining the observed dust rings and gaps in protoplanetary disks. A lower disk viscosity is favored to produce brighter rings in observations. An eccentric planet can also potentially slow down the dust radial drift in the outer region of the disk when the disk viscosity is low ($\alpha \lesssim2\times10^{-4}$) and the circularization is faster than the dust radial drift.

Future observations with next generation radio telescopes will help us to understand the presence and the evolution of magnetic fields in galaxy clusters through the determination of the so-called Rotation Measure (RM). In this work, we applied the RM-synthesis technique on synthetic SKA1-MID radio images of a pair of merging galaxy clusters, measured between 950 and 1750 MHz with a resolution of 10 arcsec and a thermal noise of 0.1$\mu$Jy/beam. The results of our RM-synthesis analysis are compared to the simulations' input parameters. We study two cases: one with radio haloes at the cluster centres, and another without. We found that the information obtained with the RM-synthesis is in general agreement with the input information. Some discrepancies are however present. We characterise them in this work, with the final goal of determining the potential impact of SKA1-MID on the study of cluster magnetic fields.

Enhancement of minor ions such as $^3$He and heavy ions in flare-associated solar energetic particle (SEP) events remains one of the major puzzles in heliophysics. In this work, we use 3D hybrid simulations (kinetic protons and fluid electrons) to investigate particle energization in a turbulent low-beta environment similar to solar flares. It is shown that in this regime the injected large-amplitude Alfv\'en waves develop into compressible and anisotropic turbulence, which efficiently heats thermal ions of different species. We find that temperature increase of heavy ions is inversely proportional to the charge to mass ratio, which is consistent with observations of impulsive SEP events. Further analysis reveals that ions are energized by interacting with nearly perpendicular magnetosonic waves near proton inertial scale.

Dust growth is often neglected when building models of protoplanetary disks due to its complexity and computational expense. However, it does play a major role in shaping the evolution of protoplanetary dust and planet formation. In this paper, we present a numerical model coupling 2-D hydrodynamic evolution of a protoplanetary disk, including a Jupiter-mass planet, and dust coagulation. This is obtained by including multiple dust fluids in a single grid-based hydrodynamic simulation and solving the Smoluchowski equation for dust coagulation on top of solving for the hydrodynamic evolution. We find that fragmentation of dust aggregates trapped in a pressure bump outside of the planetary gap leads to an enhancement in density of small grains. We compare the results obtained from the full coagulation treatment to the commonly used, fixed dust size approach and to previously applied, less computationally intensive methods for including dust coagulation. We find that the full coagulation results cannot be reproduced using the fixed-size treatment, but some can be mimicked using a relatively simple method for estimating the characteristic dust size in every grid cell.

This is the fourth in a series of companion papers showing that, when an efficient dynamo can be maintained by accretion disks around supermassive black holes in Active Galactic Nuclei (AGNs), it will lead to the formation of a powerful, magnetically-collimated helix that could explain both the observed jet/radiolobe structures on very large scales and ultimately the enormous power inferred from the observed ultra high energy cosmic rays (UHECRs) with energies > 10^19 eV. Many timescales are involved in this process. Our hyper-resistive magnetohydrodynamic (MHD) model provides a bridge between General Relativistic MHD simulations of dynamo formation, on the short accretion timescale, and observational evidence of magnetic collimation of large-scale jets on astrophysical timescales. Given the final magnetic structure, we apply hyper-resistive kinetic theory to show how instability causes slowly-evolving magnetically-collimated jets to become the most powerful relativistic accelerators in the Universe. The model yields nine observables in reasonable agreement with observations: the jet length, radiolobe radius and apparent opening angle as observed by synchrotron radiation; the synchrotron total power, synchrotron wavelengths and maximum electron energy (TeVs); and the maximum UHECR energy, the cosmic ray energy spectrum and the cosmic ray intensity on Earth.

While observations have suggested that power-law electron energy spectra are a common outcome of strong energy release during magnetic reconnection, e.g., in solar flares, kinetic simulations have not been able to provide definite evidence of power-laws in energy spectra of non-relativistic reconnection. By means of 3D large-scale fully kinetic simulations, we study the formation of power-law electron energy spectra in non-relativistic low-$\beta$ reconnection. We find that both the global spectrum integrated over the entire domain and local spectra within individual regions of the reconnection layer have power-law tails with a spectral index $p \sim 4$ in the 3D simulation, which persist throughout the non-linear reconnection phase until saturation. In contrast, the spectrum in the 2D simulation rapidly evolves and quickly becomes soft. We show that 3D effects such as self-generated turbulence and chaotic magnetic field lines enable the transport of high-energy electrons across the reconnection layer and allow them to access several main acceleration regions. This leads to a sustained and nearly constant acceleration rate for electrons at different energies. We construct a model that explains the observed power-law spectral index in terms of the dynamical balance between particle acceleration and escape from main acceleration regions, which are defined based upon a threshold for the curvature drift acceleration term. This result could be important for explaining the formation of power-law energy spectrum in solar flares.

Some circumstellar disks are observed to show prominent spiral arms in infrared scattered light or (sub-)millimeter dust continuum. The spirals might be formed from self-gravity, shadows, or planet-disk interactions. Recently, it was hypothesized that massive vortices can drive spiral arms in protoplanetary disks in a way analogous to planets. In this paper, we study the basic properties of vortex-driven spirals by the Rossby Wave Instability in 2D hydrodynamics simulations. We study how the surface density contrast, the number, and the shape of vortex-driven spirals depend on the properties of the vortex. We also compare vortex-driven spirals with those induced by planets. The surface density contrast of vortex-driven spirals in our simulations are comparable to those driven by a sub-thermal mass planet, typically a few to a few tens of Earth masses. In addition, different from the latter, the former is not sensitive to the mass of the vortex. Vortex-driven spiral arms are not expected to be detectable in current scattered light observations, and the prominent spirals observed in scattered light in a few protoplanetary disks, such as SAO 206462 (HD 135344B), MWC 758, and LkHa 330, are unlikely to be induced by the candidate vortices in them.

In a series of three papers, we introduced a novel cluster formation model that describes the formation, growth, and disruption of star clusters in high-resolution cosmological simulations. We tested this model on a Milky Way-sized galaxy and found that various properties of young massive clusters, such as the mass function and formation efficiency, are consistent with observations in the local universe. Interestingly, most massive clusters -- globular cluster candidates -- are preferentially formed during major merger events. We follow the dynamical evolution of clusters in the galactic tidal field. Due to tidal disruption, the cluster mass function evolves from initial power law to a peaked shape. The surviving clusters at $z=0$ show a broad range of metallicity [Fe/H] from -3 to -0.5. A robust prediction of the model is the age--metallicity relation, in which metal-rich clusters are systematically younger than metal-poor clusters by up to 3 Gyr.

We use the Illustris simulation to study the distributions of ellipticities and kinematic misalignments of galactic projections, as well as the intrinsic shapes and rotation of the simulated galaxies. Our results for the projections of galaxies display clear trends of an overall increase of kinematic misalignment and a slight decrease of ellipticity for fast rotators with increasing masses, while revealing no clear distinction between slow rotators of different mass. It is also found that the number of very slow rotators with large ellipticities is much larger than found in observations. The intrinsic properties of the galaxies are then analysed. The results for the intrinsic shapes of the galaxies are mostly consistent with previous results inferred from observational data. The distributions of intrinsic misalignment of the galaxies suggest that some of the galaxies produced by Illustris have significant rotation around their medium axes. Further analysis reveals that most of these galaxies display signs of non-equilibrium. We then study the evolution of the intrinsic misalignments and shapes of three specific Illustris galaxies, which we consider as typical ones, along the main progenitor line of their merger trees, revealing how mergers influence the intrinsic shapes and kinematics: the spin axis in general stays close to the shortest axis, and tends to quickly relax to such an equilibrium state within a few dynamical times of the galaxy after major perturbations; triaxiality and intrinsic flatness in general decrease with time, however, sometimes increases occur that are clearly seen to correlate with major merger events.

With the accumulation of big data of CME observations by coronagraphs, automatic detection and tracking of CMEs has proven to be crucial. The excellent performance of convolutional neural network in image classification, object detection and other computer vision tasks motivates us to apply it to CME detection and tracking as well. We have developed a new tool for CME Automatic detection and tracking with MachinE Learning (CAMEL) techniques. The system is a three-module pipeline. It is first a supervised image classification problem. We solve it by training a neural network LeNet with training labels obtained from an existing CME catalog. Those images containing CME structures are flagged as CME images. Next, to identify the CME region in each CME-flagged image, we use deep descriptor transforming to localize the common object in an image set. A following step is to apply the graph cut technique to finely tune the detected CME region. To track the CME in an image sequence, the binary images with detected CME pixels are converted from cartesian to polar coordinate. A CME event is labeled if it can move in at least two frames and reach the edge of coronagraph field of view. For each event, a few fundamental parameters are derived. The results of four representative CMEs with various characteristics are presented and compared with those from four existing automatic and manual catalogs. We find that CAMEL can detect more complete and weaker structures, and has better performance to catch a CME as early as possible.

Base on about 4,500 large tangential velocity ($V_\mathrm{tan}>0.75V_\mathrm{esc}$) with high-precision proper motions and $5\sigma$ parallaxes in Gaia DR2 5D information derived from parallax and proper motion, we identify more than 600 high velocity stars with $50\%$ unbound probability. Of these, 28 nearby (less than 6 kpc) late-type Hypervelocity stars (HVSs) with over $99\%$ possibility of unbound are discovered. In order to search for the unbound stars from the full Gaia DR2 6D phase space information derived from parallax, proper motion and radial velocity, we also identify 28 stars from the total velocity ($V_\mathrm{gc}>0.75V_\mathrm{esc}$) that have probabilities greater than $50\%$ of being unbound from the Galaxy. Of these, only three have a nearly $99\%$ probabilities of being unbound. On the whole HVSs subsample, there is 12 sources reported by other surveys. We study the spatial distribution of angular positions and angular separation of HVSs. We find the unbound HVSs are spatially anisotropic that is most significant in the Galactic longitude at more than $3\sigma$ level, and lower unbound probability HVSs are systematically more isotropic. The spatial distribution can reflect the origin of HVSs and we discuss the possible origin link with the anisotropy.

We search a large parameter space of the LkCa 15's disk density profile to fit its observed radial intensity profile of $^{12}$CO (J = 3-2) obtained from ALMA. The best-fit model within the parameter space has a disk mass of 0.1 $M_{\odot}$ (using an abundance ratio of $^{12}$CO/H$_2$ $=$ 1.4 $\times 10^{-4}$ in mass), an inner cavity of 45 AU in radius, an outer edge at $\sim$ 600 AU, and a disk surface density profile follows a power-law of the form $\rho_r \propto r^{-4}$. For the disk density profiles that can lead to a small reduced $\chi^2$ of goodness-of-fit, we find that there is a clear linear correlation between the disk mass and the power-law index $\gamma$ in the equation of disk density profile. This suggests that the $^{12}$CO disk of LkCa 15 is optically thick and we can fit its $^{12}$CO radial intensity profile using either a lower disk mass with a smaller $\gamma$ or a higher disk mass with a bigger $\gamma$. By comparing the $^{12}$CO channel maps of the best-fit model with disk models with higher or lower masses, we find that a disk mass of $\sim$ 0.1 $M_{\odot}$ can best reproduce the observed morphology of the $^{12}$CO channel maps. The dust continuum map at 0.87 mm of the LkCa 15 disk shows an inner cavity of the similar size of the best-fit gas model, but its out edge is at $\sim$ 200 AU, much smaller than the fitted gas disk. Such a discrepancy between the outer edges of the gas and dust disks is consistent with dust drifting and trapping models.

Quantum electrodynamics (QED) effects may be included in physical processes of magnetar and pulsar magnetospheres with strong magnetic fields. Involving the quantum corrections, the Maxwell electrodynamics is modified to non-linear electrodynamics. In this work, we study the force-free magnetosphere in non-linear electrodynamics in a general framework. The pulsar equation describing a steady and axisymmetric magnetosphere is derived, which now admits solutions with corrections. We derive the first-order non-linear corrections to the near-zone dipole magnetosphere in some popular non-linear effective theories. The field lines of the corrected dipole tend to converge on the rotational axis so that the fields in the polar region are stronger compared to the pure dipole case.

In this paper, we investigate the observational constraints on the scenario of vacuum energy interacting with cold dark matter. We consider eight typical interaction forms in such an interacting vacuum energy scenario. The observational data used in this work to constrain these models include the JLA sample of type Ia supernovae observation, the Planck 2015 distance priors data of cosmic microwave background anisotropies observation, the baryon acoustic oscillations data, and the Hubble constant direct measurement. We find that the current observational data almost equally favor these interacting vacuum energy models. We also find that for all these models of vacuum energy interacting with cold dark matter the case of no interaction is actually well consistent with the current observational data within 1$\sigma$ range.

By comparing cumulative linear mass profiles of 12 Gould Belt molecular clouds within 500 pc, we study how the linear mass distributions of molecular clouds vary with the angles between the molecular cloud long axes and the directions of the local magnetic fields (cloud-field direction offsets). We find that molecular clouds with the long axes perpendicular to the magnetic field directions show more even distributions of the linear mass. The result supports that magnetic field orientations can affect the fragmentation of molecular clouds (Li et al. 2017).

A recent proposal to associate 60~TeV to 2~PeV IceCube neutrino events with gamma-ray bursts~(GRBs) indicates the Lorentz violation of cosmic neutrinos and leads further to the $CPT$ symmetry violation between neutrinos and antineutrinos. Here we find that another 12 northern hemisphere track events possibly correlated with GRBs from three-year IceCube data satisfy the same regularity at a lower energy scale around 1~TeV. The combined fitting indicates a Lorentz violation scale ${E}_{\rm LV}=(6.4\pm 1.5)\times10^{17}~{ \rm GeV}$ and an intrinsic time difference ${\Delta {t}_{\rm in}=(-2.8\pm 0.7)\times10^2~{\rm s}}$, from which we find an earlier emission of neutrinos than photons at the GRB source. We also suggest analyzing neutrino events detected a few minutes before the GRB trigger time to test the $CPT$ violation of ultrahigh-energy neutrinos.

While a growing body of research indicates that relativistic magnetic reconnection is a prodigious source of particle acceleration in high-energy astrophysical systems, the dominant acceleration mechanism remains controversial. Using a combination of fully kinetic simulations and theoretical analysis, we demonstrate that Fermi-type acceleration within the large-scale motional electric fields dominates over direct acceleration from non-ideal electric fields within small-scale diffusion regions. This result has profound implications for modeling particle acceleration in large-scale astrophysical problems, since it opens up the possiblity of modeling the energetic spectra without resolving microscopic diffusion regions.

In this work, we study the chemical compositions and kinematic properties of six metal-poor stars with [Fe/H] $< -2.5$ in the Galactic halo. From high-resolution (R $\sim$~110,000) spectroscopic observations obtained with the Lick/APF, we determined individual abundances for up to 23 elements, to quantitatively evaluate our sample. We identify two carbon-enhanced metal-poor stars (J1630+0953 and J2216+0246) without enhancement in neutron-capture elements (CEMP-no stars), while the rest of our sample stars are carbon-intermediate. By comparing the light-element abundances of the CEMP stars with predicted yields from non-rotating zero-metallicity massive-star models, we find that possible the progenitors of J1630+0953 and J2216+0246 could be in the 13-25 M$_{\odot}$ mass range, with explosion energies 0.3-1.8$ \times 10^{51}$ erg. In addition, the detectable abundance ratios of light and heavy elements suggest that our sample stars are likely formed from a well-mixed gas cloud, which is consistent with previous studies. We also present a kinematic analysis, which suggests that most of our program stars likely belong to the inner-halo population, with orbits passing as close as $\sim$ 2.9 kpc from the Galactic center. We discuss the implications of these results on the critical constraints on the origin and evolution of CEMP stars, as well as the nature of the Population III progenitors of the lowest metallicity stars in our Galaxy.

We report on the highest energy photons from the Crab Nebula observed by the Tibet air shower array with the underground water-Cherenkov-type muon detector array. Based on the criterion of muon number measured in an air shower, we successfully suppress 99.92% of the cosmic-ray background events with energies $E>100$ TeV. As a result, we observed 24 photon-like events with $E>100$ TeV against 5.5 background events, which corresponds to 5.6$\sigma$ statistical significance. This is the first detection of photons with $E>100$ TeV from an astrophysical source.

We determined the chemical and kinematic properties of the Galactic thin and thick disk using a sample of 307,246 A/F/G/K-type giant stars from the LAMOST spectroscopic survey and Gaia DR2 survey. Our study found that the thick disk globally exhibits no metallicity radial gradient, but the inner disk ($R \le 8$ kpc) and the outer disk ($R>8$ kpc) have different gradients when they are studied separately. The thin disk also shows two different metallicity radial gradients for the inner disk and the outer disk, and has steep metallicity vertical gradient of d[Fe/H]/d$|z|$ $=-0.12\pm0.0007$ dex kpc$^{-1}$, but it becomes flat when it is measured at increasing radial distance, while the metallicity radial gradient becomes weaker with increasing vertical distance. Adopting a galaxy potential model, we derived the orbital eccentricity of sample stars and found a downtrend of average eccentricity with increasing metallicity for the thick disk. The variation of the rotation velocity with the metallicity shows a positive gradient for the thick disk stars and a negative one for the thin disk stars. Comparisons of our observed results with models of disk formation suggest that radial migration could have influenced the chemical evolution of the thin disk. The formation of the thick disk could be affected by more than one processes: the accretion model could play an indispensable role, while other formation mechanisms, such as the radial migration or heating scenario model could also have a contribution.

The determination of the speed of Coronal Mass Ejections (CMEs) is usually done by tracking brighter features (such as the CME front and core) in visible light coronagraphic images and by deriving unidimensional profiles of the CME speed as a function of altitude or time. Nevertheless, CMEs are usually characterized by the presence of significant density inhomogeneities propagating outward with different radial and latitudinal projected speeds, resulting in a complex evolution eventually forming the Interplanetary CME. In this work, we demonstrate for the first time how coronagraphic image sequences can be analyzed with cross-correlation technique to derive 2D maps of the almost instantaneous plasma speed distribution within the body of CMEs. The technique is first tested with the analysis of synthetic data, and then applied to real observations. Results from this work allow to characterize the distribution and time evolution of kinetic energy inside CMEs, as well as the mechanical energy (combined with the kinetic and potential energy) partition between the core and front of the CME. In the future, CMEs will be observed by two channels (VL and UV Ly-$\alpha$) coronagraphs, such as Metis on-board ESA Solar Orbiter mission as well as Ly-$\alpha$ Solar Telescope (LST) on-board Chinese Advanced Space-based Solar Observatory (ASO-S) mission. Our results will help the analysis of these future observations, helping in particular to take into account the 2D distribution of Ly-$\alpha$ Doppler dimming effect.

In this paper,we report a detailed analysis of recurrent jets originated from a location with emerging, canceling and converging negative magnetic field at the east edge of NOAA active region AR11166 from 2011 March 09 to 10. The event presented several interesting features. First, a satellite sunspot appeared and collided with a pre-existing opposite polarity magnetic field and caused a recurrent solar jet event. Second, the evolution of the jets showed blowout-like nature and standard characteristics. Third, the satellite sunspot exhibited a motion toward southeast of AR11166 and merged with the emerging flux near the opposite polarity sunspot penumbra, which afterward, due to flux convergence and cancellation episodes, caused recurrent jets. Fourth, three of the blowout jets associated with coronal mass ejections (CMEs), were observed from field of view of the Solar Terrestrial Relations Observatory. Fifth, almost all the blowout jet eruptions were accompanied with flares or with more intense brightening in the jet base region, while almost standard jets did not manifest such obvious feature during eruptions. The most important, the blowout jets were inclined to faster and larger scale than the standard jets. The standard jets instead were inclined to relative longer-lasting. The obvious shearing and twisting motions of the magnetic field may be interpreted as due to the shearing and twisting motions for a blowout jet eruption. From the statistical results, about 30% blowout jets directly developed into CMEs. It suggests that the blowout jets and CMEs should have a tight relationship.

We present different signatures of chromospheric evaporation in two solar flares observed by the Interface Region Imaging Spectrograph (IRIS). In the B1.6 flare on 2016 December 6 (SOL2016-12-06T10:40), the transition region Si IV line and the chromospheric C II and Mg II lines show blueshifts with low velocities up to 20 km/s at the flare loop footpoints in the rise phase, indicative of a gentle chromospheric evaporation. While in the C1.6 flare on 2015 December 19 (SOL2015-12-19T10:51), the Si IV, C II, and Mg II lines exhibit redshifts with velocities from several to tens of km/s at the footpoints, which might suggest an explosive chromospheric evaporation. Explosive evaporation has been observed in many flares that were captured by IRIS; however, gentle evaporation, especially manifested as blueshifts in the cool Si IV, C II, and Mg II lines, has scarcely been reported. Our results bring some new insights into chromospheric evaporation in the IRIS era.

By means of fully kinetic particle-in-cell simulations, we study whether the proton-to-electron mass ratio $m_i/m_e$ influences the energy spectrum and underlying acceleration mechanism during magnetic reconnection. While kinetic simulations are essential for studying particle acceleration during magnetic reconnection, a reduced $m_i/m_e$ is often used to alleviate the demanding computing resources, which leads to artificial scale separation between electron and proton scales. Recent kinetic simulations with high-mass-ratio have suggested new regimes of reconnection, as electron pressure anisotropy develops in the exhaust region and supports extended current layers. In this work, we study whether different $m_i/m_e$ changes the particle acceleration processes by performing a series of simulations with different mass ratio ($m_i/m_e=25-400$) and guide-field strength in a low-$\beta$ plasma. We find that mass ratio does not strongly influence reconnection rate, magnetic energy conversion, ion internal energy gain, plasma energization processes, ion energy spectra, and the acceleration mechanisms for high-energy ions. Simulations with different mass ratios are different in electron acceleration processes, including electron internal energy gain, electron energy spectrum and the acceleration efficiencies for high-energy electrons. We find that high-energy electron acceleration becomes less efficient when the mass ratio gets larger because the \textit{Fermi}-like mechanism associated with particle curvature drift becomes less efficient. These results indicate that when particle curvature drift dominates high-energy particle acceleration, the further the particle kinetic scales are from the magnetic field curvature scales ($\sim d_i$), the weaker the acceleration will be.

Energetic feedback from star clusters plays a pivotal role in shaping the dynamical evolution of giant molecular clouds (GMCs). To study the effects of stellar feedback on the star formation efficiency of the clouds and the dynamical response of embedded star clusters, we perform a suite of isolated GMC simulations with star formation and momentum feedback subgrid models using the moving-mesh hydrodynamics code \textsc{Arepo}. The properties of our simulated GMCs span a wide range of initial mass, radius, and velocity configurations. We find that the ratio of the final stellar mass to the total cloud mass, $\epsilon_{\rm int}$, scales strongly with the initial cloud surface density and momentum feedback strength. This correlation is explained by an analytic model that considers force balancing between gravity and momentum feedback. For all simulated GMCs, the stellar density profiles are systematically steeper than that of the gas at the epochs of the peaks of star formation, suggesting a centrally concentrated stellar distribution. We also find that star clusters are always in a sub-virial state with a virial parameter $\sim0.6$ prior to gas expulsion. Both the sub-virial dynamical state and steeper stellar density profiles prevent clusters from dispersal during the gas removal phase of their evolution. The final cluster bound fraction is a continuously increasing function of $\epsilon_{\rm int}$. GMCs with star formation efficiency smaller than 0.5 are still able to form clusters with large bound fractions.

In previous papers of this series, we developed a new algorithm for modeling the formation of star clusters in galaxy formation simulations. Here we investigate how dissolution of bound star clusters affects the shape of the cluster mass function and the metallicity distribution of surviving clusters. Cluster evolution includes the loss of stars that become unbound due to tidal disruption as well as mass-loss due to stellar evolution. We calculate the tidal tensor along cluster trajectories and use it to estimate the instantaneous mass-loss rate. The typical tidal tensor exhibits large variations on a time-scale of $\sim100$~Myr, with maximum eigenvalue of $10^7$~Gyr$^{-2}$, and median value of $10^4$~Gyr$^{-2}$ for the first Gyr after cluster formation. As a result of dynamical disruption, at the final available output of our simulations at redshift $z\approx1.5$, the cluster mass function has an approximately log-normal shape peaked at $\sim10^{4.3}M_\odot$. Extrapolation of the disruption to $z=0$ results in too many low-mass clusters compared to the observed Galactic globular clusters (GCs). Over 70\% of GC candidates are completely disrupted before the present; only 10\% of the total GC candidate mass remains in surviving clusters. The total mass of surviving clusters at $z=0$ varies from run to run in the range $(2-6)\times10^7M_\odot$, consistent with the observed mass of GC systems in Milky Way-sized galaxies. The metallicity distributions of all massive star clusters and of the surviving GCs have similar shapes but different normalization because of cluster disruption. The model produces a larger fraction of very metal-poor clusters than observed. A robust prediction of the model is the age-metallicity relation, in which metal-rich clusters are systematically younger than metal-poor clusters by up to 3~Gyr.

Pulsar Wind Nebula (PWN) DA 495 (G65.7+1.2) was detected in TeV gamma-rays by the High Altitude Water Cherenkov Observatory (HAWC) in 2017 (2HWC J1953+294). Follow-up observations by the Very Energetic Radiation Imaging Telescope Array System (VERITAS) confirmed the association between 2HWC J1953+294 and DA 495 and found the TeV emission to be spatially coincident with the radio emission first reported in 1968. The detection of TeV gamma-rays from DA 495, along with past X-ray detection up to 10 keV, prompted high energy X-ray observations as part of the NuSTAR Galactic Legacy Survey. We present the results of these NuSTAR observations, combined with archival Chandra and XMM-Newton observations, and confirm the previous X-ray photon index of $\Gamma_{2-20 \rm\ keV} = 2.0 \pm 0.1$. We find no spectral cutoff up to 20 keV. With the spectral information for DA 495 extended to TeV gamma-rays, we were able to perform analytical modeling to test leptonic and hadronic emission scenarios. The leptonic models can explain the broadband emission, but also imply a diffuse X-ray nebula of similar extent to the radio and TeV nebulae, which cannot be confirmed by our observations. The hadronic models can simultaneously explain the spectrum and the spatial extent in all wavelengths; however, we need a very high magnetic field strength pervading the radio and TeV nebulae and a surprisingly high particle kinetic energy. These requirements deepen the mystery of the physical nature of DA 495. Future observations in radio to infrared bands and spatially resolved $\gamma$-rays can further constrain the physical conditions and radiation mechanisms in DA 495.

We investigate how the dark energy properties impact the constraints on the total neutrino mass in interacting dark energy (IDE) models. In this study, we focus on two typical interacting dynamical dark energy models, i.e., the interacting $w$ cold dark matter (I$w$CDM) model and the interacting holographic dark energy (IHDE) model. To avoid the large-scale instability problem in IDE models, we apply the parameterized post-Friedmann approach to calculate the perturbation of dark energy. We employ the Planck 2015 cosmic microwave background temperature and polarization data, combined with low-redshift measurements on baryon acoustic oscillation distance scales, type Ia supernovae, and the Hubble constant, to constrain the cosmological parameters. We find that the dark energy properties could influence the constraint limits on the total neutrino mass. Once dynamical dark energy is considered in the IDE models, the upper bounds of $\sum m_\nu$ will be changed. By considering the values of $\chi^2_{\rm min}$, we find that in these IDE models the normal hierarchy case is slightly preferred over the inverted hierarchy case; for example, $\Delta\chi^2=2.720$ is given in the IHDE+$\sum m_\nu$ model. In addition, we also find that in the I$w$CDM+$\sum m_\nu$ model $\beta=0$ is consistent with current observational data inside the 1$\sigma$ range, and in the IHDE+$\sum m_\nu$ model $\beta>0$ is favored at more than 2$\sigma$ level.

Small stellar systems like dwarf galaxies are suggested to be the main building blocks of our Galaxy by numerical simulations in Lambda CDM models. The existence of star streams like Sagittarius tidal stream indicates that dwarf galaxies play a role in the formation of the Milky Way. However, it is unclear how many and what kind of stars in our Galaxy are originated from satellite dwarf galaxies, which could be constrained by chemical abundances of metal-poor stars. Here we report on the discovery of a metal-poor star with an extreme r-process enhancement and alpha-element deficiency. In this star, the abundance ratio of the r-process element Eu with respect to Fe is more than one order of magnitude higher than the Sun and the metallicity is 1/20 of the solar one. Such kind of stars have been found in present-day dwarf galaxies, providing the clearest chemical signature of past accretion events. The long timescale of chemical evolution of the host dwarf galaxy expected from the abundance of alpha element with respect to Fe suggests that the accretion occurred in a relatively late phase compared to most of the accretions that formed the bulk of the Milky Way halo.

The growth of solids from sub-micron to millimeter and centimeter sizes is the early step toward the formation of planets inside protoplanetary disks (PPDs). However, such processes and their potential impact on the later stages of solid growth are still poorly understood. In this work, we test the hypothesis that most disks contain at least one ringed structure with a relative small radius. We have carried out a large family of 1D two-fluid (gas+dust) hydrodynamical simulations by evolving the gas and dust motion self-consistently while allowing dust size to evolve via coagulation and fragmentation. We investigate the joint effects of ringed structures and dust size growth on the overall sub-millimeter and millimeter (mm) flux and spectral index of PPDs. Ringed structures slow down the dust radial drift and speed up the dust growth. In particular, we find that those unresolved disks with a high fragmentation velocity ($\sim10\ {\rm m~s^{-1}}$) and a high dust surface density ($\sim10\ {\rm g\ cm^{-2}}$ in the ring) can have mm spectral indices as low as $\sim2.0$, consistent with mm observations of faint disks in nearby star forming regions. Furthermore, disks with more than one ringed structure can potentially reproduce brighter disks with spectral indices lower than $\sim2.5$. Future multi-wavelength high-resolution observations of these low spectral index sources can be used to test the existence of the ringed structures in the unresolved disks and differentiate the effects of dust size growth from optical depth.

We evaluate the capabilities of Schwarzschild's orbit-superposition method by applying it to galaxies from the large scale, high resolution Illustris simulation. Nine early-type galaxies with a range of triaxiality are selected, and we create mock integral field unit data for five line-of-sight projections of each galaxy. Each of the 45 mock data sets is taken as an independent observed galaxy. Using van den Bosch's 2008 triaxial Schwarzschild implementation, we assess model estimates of various galaxy properties, covering mass profiles, intrinsic shapes, stellar orbit distributions and velocity anisotropies. Total mass within $\overline{R_{\rm e}}$ is recovered well with average deviations within $\pm15$ percent. Stellar mass is underestimated by $\sim24$ percent and dark matter overestimated by $\sim38$ percent (assuming an NFW dark matter profile and allowing for degeneracy between stellar mass and dark matter mass). Using a gNFW profile, these values improve to $\sim13$ percent for stellar mass and $\sim18$ percent for dark matter. Axis ratio estimates show a moderate bias of $\Delta (b/a)=0.07$ and $\Delta (c/a)=0.14$ ($a\ge b\ge c$). Distributions of the orbit circularities $\lambda_z$ and $\lambda_x$, representing rotation about the minor and major axes, are well reconstructed. Separating orbits into thermal categories, our models match the average fractions of these categories to within $10$ percent. Velocity anisotropy is well estimated with values matching in the inner regions but becoming slightly radially biased in the outer regions. Overall, the galaxy property estimates we obtained using Schwarzschild modelling are not implausible and are representative of the simulated galaxies we modelled.

Line lists for magnesium oxide are computed and extensive comparisons are made with existing experimental spectra. The LiTY line lists cover all ro-vibration transitions within the five lowest-lying electronic states ($X\,{}^1\Sigma^+$, $a\,{}^3\Pi$, $A\,{}^1\Pi$, $B\,{}^1\Sigma^+$ and $b\,{}^3\Sigma^{+}$) and five isotopologues: $^{24}$Mg$^{16}$O, $^{25}$Mg$^{16}$O, $^{26}$Mg$^{16}$O, $^{24}$Mg$^{17}$O, $^{24}$Mg$^{18}$O, $^{24}$Mg$^{17}$O and $^{24}$Mg$^{18}$O. The calculation use potential energy cures, spin-orbit and electronic angular momentum couplings curves determined by fitting to empirical energy levels; these levels are reproduced to within 0.01 \cm\ in most cases. Computed nuclear-motion wavefunctions are combined with {\it ab initio} dipole moment curves to give transition intensities and excited state radiative lifetimes which are compared with laboraroty measurements. The $^{24}$Mg$^{16}$O line list comprises 186 842 ($J\le 320$) ro-vibronic states and 72 833 173 transitions with angular momenta, $J$, up to 300 and covering wavenumbers up to 33 000 cm$^{-1}$ ($\lambda > 0.3$ $\mu$m). The line lists are suitable for temperatures up to about 5000 K. They are relevant to astrophysical studies of exoplanet atmospheres, cool stars and brown dwarfs, and are made available in electronic form at the CDS and ExoMol databases.

There is a longstanding discussion about whether low mass stars can form from pristine gas in the early Universe. A particular point of interest is whether we can find surviving pristine stars from the first generation in our local neighbourhood. We present here a simple analytical estimate that puts tighter constraints on the existence of such stars. In the conventional picture, should these stars have formed in significant numbers and have preserved their pristine chemical composition until today, we should have found them already. With the presented method most current predictions for survivor counts larger than zero can be ruled out.

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

We report on the discovery of five carbon-enhanced metal-poor (CEMP) stars in the metallicity range of $-3.3<$ [Fe/H] $<-2.4$. These stars were selected from the LAMOST DR3 low-resolution (R$ \sim$ 2,000) spectroscopic database as metal-poor candidates and followed-up with high-resolution spectroscopy (R$ \sim$110,000) with the LICK/APF. Stellar parameters and individual abundances for 25 chemical elements (from Li to Eu) are presented for the first time. These stars exhibit chemical abundance patterns that are similar to those reported in other literature studies of very and extremely metal-poor stars. One of our targets, J2114$-$0616, shows high enhancement in carbon ([C/Fe]=1.37), nitrogen ([N/Fe]= 1.88), barium ([Ba/Fe]=1.00), and europium ([Eu/Fe]=0.84). Such chemical abundance pattern suggests that J2114$-$0616 can be classified as CEMP-r/s star. In addition, the star J1054+0528 can be classified as a CEMP-rI star, with [Eu/Fe]=0.44 and [Ba/Fe]=$-$0.52. The other stars in our sample show no enhancements in neutron-capture elements and can be classified as CEMP-no stars. We also performed a kinematic and dynamical analysis of the sample stars based on Gaia DR2 data. The kinematic parameters, orbits, and binding energy of these stars, show that J2114$-$0616 is member of the outer halo population, while the remaining stars belong to the inner halo population but with an accreted origin. Collectively, these results add important constraints on the origin and evolution of CEMP stars as well as on their possible formation scenarios.

Low density regions are less affected by the nonlinear structure formation and baryonic physics. They are ideal places for probing the nature of dark energy, a possible explanation for the cosmic acceleration. Unlike void lensing, which requires identifications of individual voids, we study the stacked lensing signals around the low-density-positions (LDP), defined as places that are devoid of foreground bright galaxies in projection. The method allows a direct comparison with numerical results by drawing correspondence between the bright galaxies with halos. It leads to lensing signals that are significant enough for differentiating several dark energy models. In this work, we use the CFHTLenS catalogue to define LDPs, as well as measuring their background lensing signals. We consider several different definitions of the foreground bright galaxies (redshift range \& magnitude cut). Regarding the cosmological model, we run six simulations: the first set of simulations have the same initial conditions, with $\rm{w_{de}=-1,-0.5,-0.8,-1.2}$; the second set of simulations include a slightly different $\Lambda$CDM model and a w(z) model from \cite{2017NatAs...1..627Z}. The lensing results indicate that the models with $\rm{w_{de}=-0.5,-0.8}$ are not favored, and the other four models all achieve comparable agreement with the data.

We describe PSR J1926-0652, a pulsar recently discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Using sensitive single-pulse detections from FAST and long-term timing observations from the Parkes 64-m radio telescope, we probed phenomena on both long and short time scales. The FAST observations covered a wide frequency range from 270 to 800 MHz, enabling individual pulses to be studied in detail. The pulsar exhibits at least four profile components, short-term nulling lasting from 4 to 450 pulses, complex subpulse drifting behaviours and intermittency on scales of tens of minutes. While the average band spacing P3 is relatively constant across different bursts and components, significant variations in the separation of adjacent bands are seen, especially near the beginning and end of a burst. Band shapes and slopes are quite variable, especially for the trailing components and for the shorter bursts. We show that for each burst the last detectable pulse prior to emission ceasing has different properties compared to other pulses. These complexities pose challenges for the classic carousel-type models.

The tilt, rotation, or offset of each CCD with respect to the focal plane, as well as the distortion of the focal plane itself, cause shape distortions to the observed objects, an effect typically known as field distortion (FD). We point out that FD provides a unique way of quantifying the accuracy of cosmic shear measurement. The idea is to stack the shear estimators from galaxies that share similar FD-induced shape distortions. Given that the latter can be calculated with parameters from astrometric calibrations, the accuracy of the shear estimator can be directly tested on real images. It provides a way to calibrate the multiplicative and additive shear recovery biases within the scientific data itself, without requiring simulations or any external data sets. We use the CFHTLenS images to demonstrate the accuracy of the Fourier_Quad shear recovery method. We highlight some details in our image processing pipeline, including background removal, source identification and deblending, astrometric calibration, star selection for PSF reconstruction, noise reduction, etc.. We show that in the shear ranges of -0.005 < g_1 < 0.005 and -0.008 < g_2 < 0.008, the multiplicative biases are at the level of < 0.04. Slight additive biases on the order of 5E-4 (6 sigma) are identified for sources provided by the official CFHTLenS catalog (not using its shear catalog), but are minor (4 sigma) for source catalog generated by our Fourier_Quad pipeline.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

The Five-hundred-meter Aperture Spherical radio Telescope(FAST) is expected to complete its commissioning in 2019. FAST will soon begin the Commensal Radio Astronomy FasT Survey(CRAFTS), a novel and unprecedented commensal drift scan survey of the entire sky visible from FAST. The goal of CRAFTS is to cover more than 20000 $deg^{2}$ and reach redshift up to about 0.35. We provide empirical measurements of the beam size and sensitivity of FAST across the 1.05 to 1.45 GHz frequency range of the FAST L-band Array of 19-beams(FLAN). Using a simulated HI-galaxy catalogue based on the HI Mass Function(HIMF), we estimate the number of galaxies that CRAFTS may detect. At redshifts below 0.35, over $6\, \times \, 10^{5}$ HI galaxies may be detected. Below the redshift of 0.07, the CRAFTS HIMF will be complete above a mass threshold of $10^{9.5}\,M_{\odot}$. FAST will be able to investigate the environmental and redshift dependence of the HIMF to an unprecedented depth, shedding light onto the missing baryon and missing satellite problems.

We derive the mass weighted total density slopes within the effective (half-light) radius, $\gamma'$, for more than 2000 nearby galaxies from the SDSS-IV MaNGA survey using Jeans-anisotropic-models applied to IFU observations. Our galaxies span a wide range of the stellar mass ($10^9$ $M_{\rm \odot}< M_* < 10^{12}$ M$_{\odot}$) and the velocity dispersion (30 km/s $< \sigma_v <$ 300 km/s). We find that for galaxies with velocity dispersion $\sigma_v>100$ km/s, the density slope has a mean value $\langle \gamma^{\prime} \rangle = 2.24$ and a dispersion $\sigma_{\gamma}=0.22$, almost independent of velocity dispersion. A clear turn over in the $\gamma'-\sigma_v$ relation is present at $\sigma\sim 100$ km/s, below which the density slope decreases rapidly with $\sigma_v$. Our analysis shows that a large fraction of dwarf galaxies (below $M_* = 10^{10}$ M$_{\odot}$) have total density slopes shallower than 1, which implies that they may reside in cold dark matter halos with shallow density slopes. We compare our results with that of galaxies in hydrodynamical simulations of EAGLE, Illustris and IllustrisTNG projects, and find all simulations predict shallower density slopes for massive galaxies with high $\sigma_v$. Finally, we explore the dependence of $\gamma'$ on the positions of galaxies in halos, namely centrals vs. satellites, and find that for the same velocity dispersion, the amplitude of $\gamma'$ is higher for satellite galaxies by about 0.1.

TeV Halos, extended regions of TeV gamma-ray emission around middle-aged pulsars, have recently been established as a new source class in gamma-ray astronomy. These halos have been attributed to relativistic electrons and positrons that have left the acceleration region close to the pulsar and are diffusing in the surrounding medium. Measuring the morphology of TeV Halos enables for the first time a direct measurement of the electron diffusion on scales of tens of parsecs. There are hints that the presence of relativistic particles affects the diffusion rate in the pulsars' surroundings. Understanding electron diffusion is necessary to constrain the origins of the apparent `excess' of cosmic-ray positrons at tens of GeV. TeV Halos can also be used to find mis-aligned pulsars, as well as study certain properties of the Galaxy's pulsar population. Future VHE gamma-ray instruments will detect more of those TeV Halos and determine how much pulsars contribute to the observed cosmic-ray electron and positron fluxes, and how they affect diffusion in their environments.

This paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for such an investigation, and it now appears to be technologically feasible to carry out such an investigation, following the pioneering Cluster and MMS missions. Quantities that need to be measured using multipoint measurements include the two-point, two-time second correlation function of velocity, magnetic field and density, and higher order statistical objects such as third and fourth order structure functions. Some details of these requirements are given here, with a eye towards achieving closure on fundamental questions regarding the cascade rate, spectral anisotropy, characteristic coherent structures, intermittency, and dissipation mechanisms that describe plasma turbuelence, as well as its variability with plasma parameters in the solar wind. The motivation for this discussion is the current planning for a proposed Helioswarm mission that would be designed to make these measurements,leading to breakthrough understanding of the physics of space and astrophysical turbulence.

To assist with the commissioning (Jiang et al. 2019) of the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we performed a pulsar search, with the primary goal of developing and testing the pulsar data acquisition and processing pipelines. We tested and used three pipelines, two (P1 and P2 hereafter) searched for the periodic signature of pulsars whereas the other one was used to search for bright single pulses (P3 hereafter). A pulsar candidate was discovered in the observation on the 22nd August, 2017, and later confirmed by the Parkes radio telescope on the 10th September, 2017. This pulsar, named PSR J1900-0134, was the first pulsar discovered by FAST. The pulsar has a pulse period of 1.8 s and a dispersion measure (DM) of 188\,pc\,cm$^{-3}$.

Much of the baryonic matter in the Universe, including the most active and luminous sources, are best studied in the X-ray band. Key advances in X-ray optics and detectors have paved the way for the Advanced X-ray Imaging Satellite (AXIS), a Probe-class mission that is a major improvement over Chandra, which has generated a steady stream of important discoveries for the past 2 decades. AXIS can be launched in the late 2020s and will transform our understanding in several major areas of astrophysics, including the growth and fueling of supermassive black holes, galaxy formation and evolution, the microphysics of cosmic plasmas, the time-variable universe, and a wide variety of cutting-edge studies. Relative to Chandra, the AXIS PSF is nearly twice as sharp on-axis; its field of view for subarcsecond imaging 70 times larger by area; its effective area at 1 keV is 10 times larger. The low-Earth orbit ensures a low and stable detector background, resulting in 50 times greater sensitivity than Chandra for extended sources. AXIS has a rapid repointing response with operations similar to Swift, but is 100 times more sensitive for time-domain science. These capabilities open up a vast discovery space and complement the next generation of astronomical observatories. A high-spectral-resolution mission (Athena) operating at the same time as a high-angular-resolution mission (AXIS) greatly increases the range of scientific discovery. AXIS will use lightweight X-ray optics made of thin single-crystal silicon mirrors developed at NASA Goddard. The detector array builds on a long legacy of X-ray CCD and provides improved photon localization, much faster readout time, and broader energy band. The estimated mission costs are consistent with the $1B Probe mission cost guideline.

MeV blazars are the most luminous persistent sources in the Universe and emit most of their energy in the MeV band. These objects display very large jet powers and accretion luminosities and are known to host black holes with a mass often exceeding $10^9 M_{\odot}$. An MeV survey, performed by a new generation MeV telescope which will bridge the entire energy and sensitivity gap between the current generation of hard X-ray and gamma-ray instruments, will detect $>$1000 MeV blazars up to a redshift of $z=5-6$. Here we show that this would allow us: 1) to probe the formation and growth mechanisms of supermassive black holes at high redshifts, 2) to pinpoint the location of the emission region in powerful blazars, 3) to determine how accretion and black hole spin interplay to power the jet.

The Five-hundred-meter Aperture Spherical radio Telescope (FAST) was completed with its main structure installed on September 25, 2016, after which it entered the commissioning phase. This paper aims to introduce the commissioning progress of the FAST over the past two years. To improve its operational reliability and ensure effective observation time, FAST has been equipped with a real-time information system for the active reflector system and hierarchical commissioning scheme for the feed support system, which ultimately achieves safe operation of the two systems. For meeting the high-performance indices, a high-precision measurement system was set up based on the effective control methods that were implemented for the active reflector system and feed support system. Since the commissioning of the FAST, a low-frequency ultra-wideband receiver and 19-beam 1.05-1.45 GHz receiver have been mainly used. Telescope efficiency, pointing accuracy, and system noise temperature were completely tested and ultimately achieved the acceptance indices of the telescope. The FAST has been in the process of national acceptance preparations and has begun to search for pulsars. In the future, it will still strive to improve its capabilities and expand its application prospects.

The possible association of the blazar TXS 0506+056 with a high-energy neutrino detected by IceCube holds the tantalizing potential to answer three astrophysical questions: 1. Where do high-energy neutrinos originate? 2. Where are cosmic rays produced and accelerated? 3. What radiation mechanisms produce the high-energy {\gamma}-rays in blazars? The MeV gamma-ray band holds the key to these questions, because it is an excellent proxy for photo-hadronic processes in blazar jets, which also produce neutrino counterparts. Variability in MeV gamma-rays sheds light on the physical conditions and mechanisms that take place in the particle acceleration sites in blazar jets. In addition, hadronic blazar models also predict a high level of polarization fraction in the MeV band, which can unambiguously distinguish the radiation mechanism. Future MeV missions with a large field of view, high sensitivity, and polarization capabilities will play a central role in multi-messenger astronomy, since pointed, high-resolution telescopes will follow neutrino alerts only when triggered by an all-sky instrument.

The model-independent reconstruction of the energy spectra of $\overline{\nu}^{}_e$, $\nu^{}_e$ and $\nu^{}_x$ (i.e., $\nu^{}_\mu$, $\nu^{}_\tau$ and their antiparticles) from the future observation of a galactic core-collapse supernova (SN) is of crucial importance to understand the microscopic physics of SN explosions. To this end, we propose a practically useful method to combine the multi-channel detection of SN neutrinos in a large liquid-scintillator detector (e.g., JUNO), namely, the inverse beta decay $\overline{\nu}^{}_e + p \to e^+ + n$, the elastic neutrino-proton scattering $\nu + p \to \nu + p$ and the elastic neutrino-electron scattering $\nu + e^- \to \nu + e^-$, and reconstruct the energy spectra of $\overline{\nu}^{}_e$, $\nu^{}_e$ and $\nu^{}_x$ by making the best use of the observational data in those three channels. In addition, the neutrino energy spectra from the numerical simulations of the delayed neutrino-driven SN explosions are implemented to demonstrate the robustness of our method. Taking the ordinary matter effects into account, we also show how to extract the initial neutrino energy spectra in the presence of neutrino flavor conversions.

The constantly improving sensitivity of ground-based and space-borne observatories has made possible the detection of high-energy emission (X-rays and gamma-rays) from several thousands of extragalactic sources. Enormous progress has been made in measuring the continuum flux enabling us to perform imaging, spectral and timing studies. An important remaining challenge for high-energy astronomy is measuring polarization. The capability to measure polarization is being realized currently at X-ray energies (e.g. with IXPE), and sensitive gamma-ray telescopes capable of measuring polarization, such as AMEGO, AdEPT, e-ASTROGAM, etc., are being developed. These future gamma-ray telescopes will probe the radiation mechanisms and magnetic fields of relativistic jets from active galactic nuclei at spatial scales much smaller than the angular resolution achieved with continuum observations of the instrument. In this white paper, we discuss the scientific potentials of high-energy polarimetry, especially gamma-ray polarimetry, including the theoretical implications, and observational technology advances being made. In particular, we will explore the primary scientific opportunities and wealth of information expected from synergy of multi-wavelength polarimetry that will be brought to multi-messenger astronomy.

The high-energy universe has revealed that energetic particles are ubiquitous in the cosmos and play a vital role in the cultivation of cosmic environments on all scales. Energetic particles in our own galaxy, galactic cosmic rays (GCRs), engage in a complex interplay with the interstellar medium and magnetic fields in the galaxy, giving rise to many of its key characteristics. This White Paper is the first of a two-part series highlighting the most well-known high-energy cosmic accelerators and contributions that MeV gamma-ray astronomy will bring to understanding their energetic particle phenomena. The focus of this white paper is galactic cosmic rays, supernova remnants, protostellar jets and superbubbles, and colliding wind binaries.

The high-energy universe has revealed that energetic particles are ubiquitous in the cosmos and play a vital role in the cultivation of cosmic environments on all scales. Though they play a key role in cultivating the cosmological environment and/or enabling our studies of it, there is still much we do not know about AGNs and GRBs, particularly the avenue in which and through which they supply radiation and energetic particles, namely their jets. This White Paper is the second of a two-part series highlighting the most well-known high-energy cosmic accelerators and contributions that MeV gamma-ray astronomy will bring to understanding their energetic particle phenomena. The focus of this white paper is active galactic nuclei and gamma-ray bursts.

The characteristic two-component blazar spectral energy distribution (SED) can be of either leptonic and/or hadronic origins. The potential association of the high-energy neutrino event IceCube-170922A with the flaring blazar TXS~0506+056 indicates that hadronic processes may operate in a blazar jet. Despite multi-wavelength follow-ups of the event and extensive theoretical modelings, the radiation mechanisms and the underlying magnetic field strength and configuration remain poorly understood. In this paper, we consider generic leptonic and hadronic blazar spectral models with distinct magnetic field strengths and radiation mechanisms. We analytically reproduce the SEDs and the neutrino flux of hadronic models, and predict their X-ray to $\gamma$-ray polarization degrees. Furthermore, by performing relativistic magnetohydrodynamic (RMHD) simulations taking into account the polarization-dependent radiation transfer, we study the time-dependent multi-wavelength polarization variability of the proton synchrotron model under a shock scenario. Our results suggest that the high-energy polarization degree and the neutrino flux can be jointly used to pinpoint the leptonic and/or hadronic blazar radiation mechanisms in the X-ray and $\gamma$-ray bands, and to infer the magnetic field strength in the emission region. Additionally, the temporal multi-wavelength polarization signatures in the proton synchrotron model shed light on the jet energy composition and the dynamical importance of magnetic fields in the blazar emission region. Future multi-wavelength polarimetry facilities such as {\it IXPE} and {\it AMEGO} together with neutrino telescopes such as {\it IceCube} can provide unprecedented observational constraints to probe the blazar radiation mechanisms and jet dynamics.

We explore the properties of dense star clusters that are likely to be the nurseries of stellar black holes pairing in close binaries. We combine a cosmological model of globular cluster formation with analytic prescriptions for the dynamical assembly and evolution of black hole binaries to constrain which types of clusters are most likely to form binaries tight enough to coalesce within a Hubble time. We find that black hole binaries which are ejected and later merge ex-situ form in clusters of a characteristic mass $M_{\rm cl} \sim 10^{5.3}M_{\odot}$, whereas binaries which merge in-situ form in more massive clusters, $M_{\rm cl} \sim 10^{5.7}M_{\odot}$. The clusters which dominate the production of black hole binaries are similar in age and metallicity to the entire population. Finally, we estimate an approximate cosmic black hole merger rate of dynamically assembled binaries using the mean black hole mass for each cluster given its metallicity. We find an intrinsic rate of $\sim 6\,\mathrm{Gpc\,yr^{-1}}$ at $z=0$, a weakly increasing merger rate out to $z=1.5$, and then a decrease out to $z=4$. Our results can be used to provide a cosmological context and choose initial conditions in numerical studies of black hole binaries assembled in star clusters.

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.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

Based on the second Gaia data release (Gaia DR2) and spectroscopy from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) Data, we identified 20,089 halo stars kinematically and chemically. The halo streams in the solar neighborhood could be detected in the space of energy and angular momentum. We reshuffle the velocities of these stars to determine the significance of substructure. Finally, we identify 4 statistically significant substructures that are labeled GL-1 through 4. Among these substructures, GL-1 is previously known stream ("N2" stream) and the rest 3 substructures are new. These substructures may be the debris of dwarf galaxies accretion event, their dynamical and chemical information can help to understand the history of the Milky Way.

In this work, we consider the possibility of energy release in pulsar magnetospheres deformed by gravitational waves from nearby sources. The strong electromagnetic fields in the magnetospheres may release non-negligible energy despite the weakness of the gravitational wave. When the background spacetime is perturbed due to the passage of a gravitational wave, the original force-free state of the inner magnetosphere will be slightly violated. The plasma-filled magnetosphere tends to evolve into new force-free states as the spacetime varies with time. During this process, a small portion of the electromagnetic energy stored in the magnetosphere will be released to the acceleration of charged particles along the magnetic field lines. When the pulsar is close enough to the gravitational wave source (e.g., $\sim10^{-2}$ pc to the gravitational wave sources observed recently), the resulting energy loss rate is comparable with the radio luminosity of the pulsar. It is also noticed that, under very stringent conditions (for magnetars with much shorter distance to the sources), the released energy can reach the typical energy observed from fast radio bursts (FRBs).

The linewidth of ions has been observed to be systematically narrower than that of the coexisting neutrals in molecular clouds (Houde et al. 2000) and been interpreted as the signature of the decoupling of the neutral turbulence from magnetic fields in partially ionized medium (Li & Houde 2008; Paper I, hereafter). As a sequel of Paper I, here we present further observational evidence that lend support to these earlier proposals with the velocity coordinate spectrum analysis (Lazarian & Pogosyan 2006). We recover the turbulent energy spectra of HCN and HCO+(4-3) in a starless molecular filament in NGC 6334 where magnetic fields play a dynamically important role (Li et al. 2015). Our analysis showed that the neutral spectrum is consistent with Kolmogorov-type (k^-5/3, where k is the wave number), while that of the ions is the same on the large scale but steeper (k^-2) for scales smaller than 0.404pc. We carefully ruled out the possibilities that the spectrum difference can stem from the differences of ion and neutral optical depth and hyper-fine structures.

The hadronic interaction of cosmic rays with solar atmosphere can produce high energy gamma rays. The gamma-ray luminosity is correlated both with the flux of primary cosmic rays and the intensity of the solar magnetic field. The gamma rays below 200 GeV have been observed by $Fermi$ without any evident energy cutoff. The bright gamma-ray flux above 100 GeV has been detected only during solar minimum. The only available data in TeV range come from the HAWC observations, however outside the solar minimum. The ARGO-YBJ dataset has been used to search for sub-TeV/TeV gamma rays from the Sun during the solar minimum from 2008 to 2010, the same time period covered by the Fermi data. A suitable model containing the Sun shadow, solar disk emission and inverse-Compton emission has been developed, and the chi-square minimization method was used to quantitatively estimate the disk gamma-ray signal. The result shows that no significant gamma-ray signal is detected and upper limits to the gamma-ray flux at 0.3$-$7 TeV are set at 95\% confidence level. In the low energy range these limits are consistent with the extrapolation of the Fermi-LAT measurements taken during solar minimum and are compatible with a softening of the gamma-ray spectrum below 1 TeV. They provide also an experimental upper bound to any solar disk emission at TeV energies. Models of dark matter annihilation via long-lived mediators predicting gamma-ray fluxes > $10^{-7}$ GeV $cm^{-2}$ $s^{-1}$ below 1 TeV are ruled out by the ARGO-YBJ limits.

Gamma-ray bursts are the strongest explosions in the Universe since the Big Bang, believed to be produced either in forming black holes at the end of massive star evolution or merging of compact objects. Spectral and timing properties of gamma-ray bursts suggest that the observed bright gamma-rays are produced in the most relativistic jets in the Universe; however, the physical properties, especially the structure and magnetic topologies in the jets are still not well known, despite several decades of studies. It is widely believed that precise measurements of the polarization properties of gamma-ray bursts should provide crucial information on the highly relativistic jets. As a result there have been many reports of gamma-ray burst polarization measurements with diverse results, see, however many such measurements suffered from substantial uncertainties, mostly systematic. After the first successful measurements by the GAP and COSI instruments, here we report a statistically meaningful sample of precise polarization measurements, obtained with the dedicated gamma-ray burst polarimeter, POLAR onboard China's Tiangong-2 spacelab. Our results suggest that the gamma-ray emission is at most polarized at a level lower than some popular models have predicted; although our results also show intrapulse evolution of the polarization angle. This indicates that the low polarization degrees could be due to an evolving polarization angle during a gamma-ray burst.

We (Li et al. 2009; Paper-I) compared the magnetic field directions inferred from polarimetry data obtained from 100-pc scale inter-cloud media (ICM) and from sub-pc scale molecular cloud cores. The highly correlated result led us to conclude that cloud turbulence must be sub-Alfvenic. Here we extend the study with 0.01-pc cores observed by interferometers. The inferred field directions at this scale significantly deviate from that of the surrounding ICM. An obvious question to ask is whether this high-resolution result contradicts the sub-Alfvenic picture concluded earlier. We performed MHD simulations of a slightly super-critical (magnetic criticality = 2) clouds with Alfvenic Mach number $M_A = 0.63$, which can reproduce the Paper-I results, and observed the development towards smaller scales. Interestingly, all subregions hosting cores with $n_H$$_2$ > $10^{5}$/cc (the typical density observed by interferometers) possess $M_A = 2-3$. Not too surprisingly, these slightly super-Alfvenic cores result in B-field orientation offsets comparable to the interferometer observations. The result suggests that gravity can concentrate (and maybe also contribute to, which takes more study to confirm) turbulent energy and create slightly super-Alfvenic cores out from sub-Alfvenic clouds. The results of our simulations also agree with the observed velocity-scale (Kauffmann et al. 2013), mass-scale (Lombardi et al. 2010) and field strength-density (Li et al. 2015; Crutcher et al. 2010) relations.

We investigate the resolved relation between local extinction and star formation surface density within nearby star-forming galaxies selected from the MaNGA survey. Balmer decrement measurements imply an extinction of the H{\alpha} line emission which scales approximately linearly with the logarithm of the star formation surface density: $ A_{H{\alpha}} = 0.46 \log(\Sigma_{SFR}) + 1.53$. Secondary dependencies are observed such that, at a given $\Sigma_{SFR}$, regions of lower metallicity and/or enhanced H{\alpha} equivalent width (EW) suffer less obscuration than regions of higher metallicity and/or lower H{\alpha} EW. Spaxels lying above the mean relation also tend to belong to galaxies that are more massive, larger and viewed under higher inclination than average. We present a simple model in which the observed trends can be accounted for by a metallicity-dependent scaling between $\Sigma_{SFR}$ and $\Sigma_{dust}$ via a super-linear Kennicutt-Schmidt relation ($n_{KS} \sim 1.47$) and a dust-to-gas ratio which scales linearly with metallicity (DGR($Z_{\odot}$) = 0.01). The relation between the resulting total dust column and observed effective extinction towards nebular regions requires a geometry for the relative distribution of H{\alpha} emitting regions and dust that deviates from a uniform foreground screen and also from an entirely homogeneous mixture of dust and emitting sources. The best-fit model features an H{\alpha} EW and galactocentric distance dependent fraction of the dust mass in a clumpy foreground screen in front of a homogeneous mixture.

Quasi-periodic fast-propagating (QFP) magnetosonic waves and extreme ultraviolet (EUV) waves were proposed to be driven by solar flares and coronal mass ejections (CMEs), respectively. In this Letter, we present a detailed analysis of an interesting event in which we find that both QFP magnetosonic waves and EUV waves are excited simultaneously in one solar eruption event. The co-existence of the two wave phenomena offers an excellent opportunity to explore their driving mechanisms. The QFP waves propagate in a funnel-like loop system with a speed of 682--837 \speed{} and a lifetime of 2 minutes. On the contrary, the EUV waves, which present a faster component and a slower component, propagate in a wide angular extent, experiencing reflection and refraction across a magnetic quasi-separatrix layer. The faster component of the EUV waves travels with a speed of 412--1287 \speed{}, whereas the slower component travels with a speed of 246--390 \speed{}. The lifetime of the EUV waves is $\sim$15 minutes. It is revealed that the faster component of the EUV waves is cospatial with the first wavefront of the QFP wave train. Besides, The QFP waves have a period of about 45$\pm$5 seconds, which is absent in the associated flares. All these results imply that QFP waves can also be excited by mass ejections, including CMEs or jets.

In this paper, we present a detailed analysis of a coronal blowout jet eruption which was associated with an obvious extreme-ultraviolet (EUV) wave and one complicated coronal mass ejection (CME) event based on the multi-wavelength and multi-view-angle observations from {\sl Solar Dynamics Observatory} and {\sl Solar Terrestrial Relations Observatory}. It is found that the triggering of the blowout jet was due to the emergence and cancellation of magnetic fluxes on the photosphere. During the rising stage of the jet, the EUV wave appeared just ahead of the jet top, lasting about 4 minutes and at a speed of 458 - \speed{762}. In addition, obvious dark material is observed along the EUV jet body, which confirms the observation of a mini-filament eruption at the jet base in the chromosphere. Interestingly, two distinct but overlapped CME structures can be observed in corona together with the eruption of the blowout jet. One is in narrow jet-shape, while the other one is in bubble-shape. The jet-shaped component was unambiguously related with the outwardly running jet itself, while the bubble-like one might either be produced due to the reconstruction of the high coronal fields or by the internal reconnection during the mini-filament ejection according to the double-CME blowout jet model firstly proposed by Shen et al. (2012b), suggesting more observational evidence should be supplied to clear the current ambiguity based on large samples of blowout jets in future studies.

Based on the second Gaia data (Gaia DR2) and spectroscopy from the LAMOST Data Release 5, we defined the high-velocity (HiVel) stars sample as those stars with $v_{\mathrm{gc}} > 0.85 v_{\mathrm{esc}}$, and derived the final sample of 24 HiVel stars with stellar astrometric parameters and radial velocities. Most of the HiVel stars are metal-poor and $\alpha$-enhanced. In order to further explore the origin of these HiVel stars, we traced the backwards orbits of each HiVel star in the Galactic potential to derive probability parameters which are used to classify these HiVel stars. Of these, 5 stars are from the tidal debris of disrupted dwarf galaxy and 19 stars are runaway-star candidates which originate from the stellar disk.

We explore orbit properties of 35 prolate-triaxial galaxies selected from the Illustris cosmological hydrodynamic simulation. We present a detailed study of their orbit families, and also analyse relations between the relative abundance of the orbit families and the spin parameter, triaxiality, the ratio of the angular momentum and the baryon fraction. We find that box orbits dominate the orbit structure for most prolate-triaxial galaxies, especially in the central region. The fraction of irregular orbits in the prolate-triaxial galaxies is small, and it increases with galaxy radius. Both the x-tube and z-tube orbits are important in prolate-triaxial galaxies, especially in the outer regions. The fraction of box orbits for prolate-triaxial galaxies ($0.7<T<1$) decreases with the triaxiality of the stars, while the fraction of x-tube orbits increases with the triaxiality for a given axis ratio. The fraction of box orbits increases and the fractions of x-tube and z-tube orbits weakly decrease with increasing baryon fraction. These results help to understand the structure of prolate-triaxial galaxies and provide cross-checks for constructing dynamical models of prolate-triaxial galaxies by using the Schwarzschild or the Made-to-Measure methods.

Various studies have implied the existence of a gaseous halo around the Galaxy extending out to 100 kpc. Galactic cosmic rays (CRs) that propagate to the halo, either by diffusion or by convection with the possibly existing large-scale Galactic wind, can interact with the gas therein and produce gamma-rays via proton-proton collision. We calculate the cosmic ray distribution in the halo and the gamma-ray flux, and explore the dependence of the result on model parameters such as diffusion coefficient, CR luminosity, CR spectral index. We find that the current measurement of isotropic gamma-ray background at $\lesssim$TeV with Fermi Large Area Telescope already approaches a level that can provide interesting constraints on the properties of Galactic cosmic ray (e.g., with CR luminosity $L_{CR}\leq 10^{41}$erg/s). We also discuss the possibilities of the Fermi bubble and IceCube neutrinos originating from the proton-proton collision between cosmic rays and gas in the halo, as well as the implication of our results for the baryon budget of the hot circumgalactic medium of our Galaxy. Given that the isotropic gamma-ray background is likely to be dominated by unresolved extragalactic sources, future telescopes may extract more individual sources from the IGRB, and hence put even more stringent restriction on the relevant quantities (such as Galactic cosmic ray luminosity and baryon budget in the halo) in the presence of a turbulent halo that we consider.

The Advanced Space-based Solar Observatory (ASO-S) mission aims to explore two most spectacular eruptions in the Sun: solar flares and coronal mass ejections (CMEs), and their magnetism. For the studies of CMEs, the payload Lyman-alpha Solar Telescope (LST) has been proposed. It includes a traditional white-light coronagraph and a Lyman-alpha coronagraph which opens a new window to CME observations. Polarization measurements taken by white-light coronagraphs are crucial to derive fundamental physical parameters of CMEs. To make such measurements, there are two options of Stokes polarimeter which have been used by existing white-light coronagraphs for space missions. One uses a single or triple linear polarizers, the other involves both a half-wave plate and a linear polarizer. We find that the former option subjects to less uncertainty in the derived Stokes vector propagated from detector noise. The latter option involves two plates which are prone to internal reflections and may have a reduced transmission factor. Therefore, the former option is adopted as our Stokes polarimeter scheme for LST. Based on the parameters of the intended linear polarizer(s) colorPol provided by CODIXX and the half-wave plate 2-APW-L2-012C by Altechna, it is further shown that the imperfect maximum transmittance of the polarizer significantly increases the variance amplification of Stokes vector by at least about 50% when compared with the ideal case. The relative errors of Stokes vector caused by the imperfection of colorPol polarizer and the uncertainty due to the polarizer assembling in the telescope are estimated to be about 5%. Among the considered parameters, we find that the dominant error comes from the uncertainty in the maximum transmittance of the polarizer.

We investigate the structure of galaxies formed in a suite of high-resolution cosmological simulations. Consistent with observations of high-redshift galaxies, our simulated galaxies show irregular, prolate shapes with thick stellar disks, which are dominated by turbulent motions instead of rotation. Yet molecular gas and young stars are restricted to relatively thin disks. We examine the accuracy of applying the Toomre linear stability analysis to predict the location and amount of gas available for star formation. We find that the Toomre criterion still works for these irregular galaxies, after correcting for multiple gas and stellar components: the $Q$ parameter in $\rm{H_2}$ rich regions is in the range $0.5-1$, remarkably close to unity. Due to the violent stellar feedback from supernovae and strong turbulent motions, young stars and molecular gas are not always spatially associated. Neither the $Q$ map nor the $\rm{H_2}$ surface density map coincide with recent star formation exactly. We argue that the Toomre criterion is a better indicator of future star formation than a single $\rm{H_2}$ surface density threshold because of the smaller dynamic range of $Q$. The depletion time of molecular gas is below 1~Gyr on kpc scale, but with large scatter. Centering the aperture on density peaks of gas/young stars systematically biases the depletion time to larger/smaller values and increases the scatter.

The discovery of thousands of exoplanets has shown that the birth of planets is a very efficient process in nature. Several physical mechanisms have been proposed to describe the assembly of planets in disks surrounding young stars. However, observational constraints have been sparse on account of insufficient sensitivity and resolution, especially for imaging the inner $10 - 20$ au from the star, where most of planets are observed. Thanks to its unprecedented angular resolution and sensitivity at wavelengths where the emission from the circumstellar material is optically thin, the ngVLA has the potential to transform our understanding of planet formation. In this chapter, we present state-of-the-art theoretical models of planetary systems in the act of forming that demonstrate ngVLA capabilities of imaging and follow the temporal evolution of young solar system analogues up to the distance of the Orion Nebula. These images will unveil how planets form and interact with their parental disks, and will shed light on the diverse properties of exoplanetary systems.

The MaNGA project has obtained IFU data for several thousand nearby galaxies, including barred galaxies. With the two dimensional spectral and kinematic information provided by IFUs, we can measure the pattern speed of a barred galaxy, which determines the bar dynamics. We apply the non-parametric method proposed by Tremaine \& Weinberg to estimate the bar pattern speed for 53 barred galaxies, making this the largest sample studied so far in this way. Our sample is selected from the MaNGA first public data release as part of SDSS Data Release 13 according mainly to the axis ratio and position angle difference between the bar and disc, while kinematic data is from the later SDSS Data Release 14. We have used both the photometric position angle from the photometric image and the kinematic position angle from the stellar velocity map to derive the pattern speed. Combining three independent bar length measurements and the circular velocity from Jeans Anisotropic modelling (JAM), we also determine the dimensionless ratio $\cal{R}$ of the corotation radius to the bar length. We find that the galaxy's position angle is the main uncertainty in determining the bar pattern speed. The kinematic position angle leads to fewer ultrafast bars than the photometric position angle, and this could be due to the method of measuring the kinematic position angle. We study the dependence of $\cal{R}$ values on galaxy properties such as the dark matter fraction from JAM modelling and the stellar age and metallicity from stellar population synthesis (SPS). A positive correlation between the bar length and bar strength is found: the longer the bar, the stronger the bar. However, no other significant correlations are found. This may result from errors in deriving the $\cal{R}$ values or from the complex formation and slowdown processes of galactic bars.

In solar flares and other astrophysical systems, a major challenge for solving particle acceleration problem associated with magnetic reconnection is the enormous scale separation between kinetic scales and observed reconnection scale. Because of this, it has been difficult to draw any definite conclusions by just using kinetic simulations. Particle acceleration model that solves energetic particle transport equation can capture the main acceleration physics found in kinetic simulations, and thus provide a practical way to make observable predictions and directly compare model results with observations. Here we study compression particle acceleration in magnetic reconnection by solving Parker (diffusion-advection) transport equation using velocity and magnetic fields from two-dimensional high-Lundquist-number magnetohydrodynamics (MHD) simulations of a low-$\beta$ reconnection layer. We show that the compressible reconnection layer can give significant particle acceleration, leading to the formation of power-law particle energy distributions. We analyze the acceleration rate and find that the acceleration in the reconnection layer is a mixture of first order and second order Fermi processes. When including a guide field, we find the spectrum becomes steeper and both the power-law cutoff energy and maximum particle energy decrease as plasma becomes less compressible. This model produces 2D particle distribution that one can use to generate radiation map and directly compare with solar flare observations. This provides a framework to explain particle acceleration at large-scale astrophysical reconnection sites, such as solar flares.

SS 433 is a binary system containing a supergiant star that is overflowing its Roche lobe with matter accreting onto a compact object (either a black hole or neutron star). Two jets of ionized matter with a bulk velocity of $\sim0.26c$ extend from the binary, perpendicular to the line of sight, and terminate inside W50, a supernova remnant that is being distorted by the jets. SS 433 differs from other microquasars in that the accretion is believed to be super-Eddington, and the luminosity of the system is $\sim10^{40}$ erg s$^{-1}$. The lobes of W50 in which the jets terminate, about 40 pc from the central source, are expected to accelerate charged particles, and indeed radio and X-ray emission consistent with electron synchrotron emission in a magnetic field have been observed. At higher energies (>100 GeV), the particle fluxes of $\gamma$ rays from X-ray hotspots around SS 433 have been reported as flux upper limits. In this energy regime, it has been unclear whether the emission is dominated by electrons that are interacting with photons from the cosmic microwave background through inverse-Compton scattering or by protons interacting with the ambient gas. Here we report TeV $\gamma$-ray observations of the SS 433/W50 system where the lobes are spatially resolved. The TeV emission is localized to structures in the lobes, far from the center of the system where the jets are formed. We have measured photon energies of at least 25 TeV, and these are certainly not Doppler boosted, because of the viewing geometry. We conclude that the emission from radio to TeV energies is consistent with a single population of electrons with energies extending to at least hundreds of TeV in a magnetic field of $\sim16$~micro-Gauss.

In the present paper, we investigate three scalar fields, \qu field, phantom field and tachyon field, to explore the source of dark energy, using the Gaussian processes method from the background data and perturbation growth rate data. The corresponding reconstructions all suggest that the dark energy should be dynamical. Moreover, the quintom field, a combination between \qu field and phantom field, is powerfully favored by the data within 68\% confidence level. Using the mean values of scalar field $\phi$ and potential $V$, we fit the function $V(\phi)$ in different fields. The fitted results imply that potential $V(\phi)$ in each scalar field may be a double exponential function or Gaussian function. The Gaussian processes reconstructions also indicate that the tachyon scalar field cannot be convincingly favored by the data and is at a disadvantage to describe the dark energy.

Supermassive black holes (SMBHs) are the central engines of luminous quasars and are found in most massive galaxies today. But the recent discoveries of ULAS J1120+0641, a $2 \times 10^9$ M$_{\odot}$ BH at $z =$ 7.1, and ULAS J1342+0928, a $8.0 \times 10^{8}$ M$_{\odot}$ BH at $z =$ 7.5, now push the era of quasar formation up to just 690 Myr after the Big Bang. Here we report new cosmological simulations of SMBHs with X-rays fully coupled to primordial chemistry and hydrodynamics that show that J1120 and J1342 can form from direct collapse black holes (DCBHs) if their growth is fed by cold, dense accretion streams, like those thought to fuel rapid star formation in some galaxies at later epochs. Our models reproduce all of the observed properties of J1120: its mass, luminosity, and H II region as well as star formation rates and metallicities in its host galaxy. They also reproduce the dynamical mass of the innermost 1.5 kpc of its emission region recently measured by ALMA and J-band magnitudes that are in good agreement with those found by the VISTA Hemisphere Survey.

Several nearby protoplanetary disks have been observed to display large scale crescents in the (sub)millimeter dust continuum emission. One interpretation is that these structures correspond to anticyclonic vortices generated by the Rossby wave instability within the gaseous disk. Such vortices have local gas over-densities and are expected to concentrate dust particles with Stokes number around unity. This process might catalyze the formation of planetesimals. Whereas recent observations showed that dust crescent are indeed regions where millimeter-size particles have abnormally high concentration relative to the gas and smaller grains, no observations have yet shown that the gas within the crescent region counter-rotates with respect to the protoplanetary disk. Here we investigate the detectability of anticyclonic features through measurement of the line-of-sight component of the gas velocity obtained with ALMA. We carry out 2D hydrodynamic simulations and 3D radiative transfer calculation of a protoplanetary disk characterized by a vortex created by the tidal interaction with a massive planet. As a case study, the disk parameters are chosen to mimic the IRS 48 system, which has the most prominent crescent observed to date. We generate synthetic ALMA observations of both the dust continuum and 12CO emission around the frequency of 345 GHz. We find that the anticyclonic features of vortex are weak but can be detected if both the source and the observational setup are properly chosen. We provide a recipe for maximizing the probability to detect such vortex features and present an analysis procedure to infer their kinematic properties.

ALMA has found multiple dust gaps and rings in a number of protoplanetary disks in continuum emission at millimeter wavelengths. The origin of such structures is in debate. Recently, we documented how one super-Earth planet can open multiple (up to five) dust gaps in a disk with low viscosity ($\alpha\lesssim10^{-4}$). In this paper, we examine how the positions, depths, and total number of gaps opened by one planet depend on input parameters, and apply our results to real systems. Gap locations (equivalently, spacings) are the easiest metric to use when making comparisons between theory and observations, as positions can be robustly measured. We fit the locations of gaps empirically as functions of planet mass and disk aspect ratio. We find that the locations of the double gaps in HL Tau and TW Hya, and of all three gaps in HD 163296, are consistent with being opened by a sub-Saturn mass planet. This scenario predicts the locations of other gaps in HL Tau and TW Hya, some of which appear consistent with current observations. We also show how the Rossby wave instability may develop at the edges of several gaps and result in multiple dusty vortices, all caused by one planet. A planet as low in mass as Mars may produce multiple dust gaps in the terrestrial planet forming region.

About one percent of giants\upcite{Brown1989} are detected to have anomalously high lithium (Li) abundances in their atmospheres, conflicting directly with the prediction of the standard stellar evolution models\upcite{Iben1967}, and making the production and evolution of Li more intriguing, not only in the sense of the Big Bang nucleosynthesis\upcite{Cyburt2016,Spite1982} or the Galactic medium\upcite{Tajitsu2015}, but also the evolution of stars. Decades of efforts have been put into explaining why such outliers exist\upcite{Sackmann1999, Denissenkov2004, Charbonnel2010}, yet the origins of Li-rich giants are still being debated. Here we report the discovery of the most Li-rich giant known to date, with a super-high Li abundance of 4.51. This rare phenomenon was snapshotted together with another short-term event that the star is experiencing its luminosity bump on the red giant branch. Such high Li abundance indicates that the star might be at the very beginning of its Li-rich phase, which provides a great opportunity to investigate the origin and evolution of Li in the Galaxy. A detailed nuclear simulation is presented with up-to-date reaction rates to recreate the Li enriching process in this star. Our results provide tight constraints on both observational and theoretical points of view, suggesting that low-mass giants can produce Li inside themselves to a super high level via $^{7}$Be transportation during the red giant phase.

Globular cluster (GC) systems around galaxies of a vast mass range show remarkably simple scaling relations. The combined mass of all GCs is a constant fraction of the total galaxy mass and the mean metallicity and metallicity dispersion of the GC system scale up weakly with galaxy mass. The metallicity of massive, metal-poor ("blue") clusters increases with cluster mass, while that of metal-rich ("red") clusters does not. A significant age-metallicity relation emerges from analysis of resolved stellar populations in Galactic GCs and unresolved populations in nearby galaxies. Remarkably, all these trends can be explained by a simple merger-based model developed in previous work and updated here using recent observations of galaxy scaling relations at high redshift. We show that the increasing dispersion of GC metallicity distributions with galaxy mass is a robust prediction of the model. It arises from more massive galaxies having more mergers that combine satellite GC systems. The average metallicity also increases by 0.6~dex over 3~dex in halo mass. The models show a non-linear trend between the GC system mass and host galaxy mass which is consistent with the data. The model does not consider GC self-enrichment, yet predicts a correlation between cluster mass and metallicity for massive blue clusters. The age-metallicity relation is another robust prediction of the model. Half of all clusters are predicted to form within the redshift range $5<z<2.3$, corresponding to ages of $10.8-12.5$~Gyr, in halos of masses $10^{11}-10^{12.5} M_{\odot}$.

We present extensive spectroscopic observations for one of the closest type Ia supernovae (SNe Ia), SN 2014J discovered in M82, ranging from 10.4 days before to 473.2 days after B-band maximum light. The diffuse interstellar band (DIB) features detected in a high-resolution spectrum allow an estimate of line-of-sight extinction as Av=1.9+/-0.6 mag. Spectroscopically, SN 2014J can be put into the high-velocity (HV) subgroup in Wang's classification with a velocity of Si~II 6355 at maximum light as about 12200 km/s, but has a low velocity gradient (LVG, following Benetti's classification) as 41+/-2 km/s/day, which is inconsistent with the trend that HV SNe Ia generally have larger velocity gradients. We find that the HV SNe Ia with LVGs tend to have relatively stronger Si III (at ~4400 Angstrom) absorptions in early spectra, larger ratios of S II 5468 to S II 5640, and weaker Si II 5972 absorptions compared to their counterparts with similar velocities but high velocity gradients. This shows that the HV+LVG subgroup of SNe Ia may have intrinsically higher photospheric temperature, which indicates that their progenitors may experience more complete burning in the explosions relative to the typical HV SNe Ia.

The association of the high-energy neutrino event IceCube-170922A with the flaring blazar TXS~0506+056 indicates that hadronic processes may operate in a blazar jet. We perform semi-analytical spectral fitting of the multi-wavelength emission to obtain estimates of the jet physical parameters, and find that the multi-wavelength emission can be explained by either a proton synchrotron scenario or an electron inverse Compton scattering scenario. In the proton synchrotron scenario, a strong magnetic field of $10-100$~G is required, implying that the particle acceleration is likely driven by magnetic energy dissipation such as magnetic reconnection events. The inverse Compton scenario implies a magnetic field of $0.1-1$~G. Thus the particle acceleration is likely driven by the kinetic energy dissipation such as shocks. We also discuss the neutrino production in the context of single-zone and multi-zone models based on the above two scenarios. We demonstrate that the variability and optical polarization signatures can be used to distinguish the two scenarios due to their drastically different magnetic field. Specifically, the proton synchrotron scenario may show orphan fast variability in the low-energy spectral component on top of the active state, with an optical polarization degree $\lesssim 10\%$ throughout the active state. The inverse Compton scattering scenario instead predicts co-variability of the low- and high-energy components on both short and long time scales, as well as a strongly variable optical polarization degree that can reach $\gtrsim 20\%$. Our results suggest that optical polarization measurements and well-resolved multi-wavelength light curves can be used to understand the electromagnetic and high-energy neutrino emissions by TXS~0506+056 and similar events in the future.

We report on the elemental abundances of the carbon-enhanced metal-poor (CEMP) star J2217+2104 discovered by our metal-poor star survey with LAMOST and Subaru. This object is a red giant having extremely low Fe abundance ([Fe/H]=-4.0) and very large enhancement of C, N, and O with excesses of Na, Mg, Al, and Si. This star is a new example of a small group of such CEMP stars identified by previous studies. We find a very similar abundance pattern for O-Zn in this class of objects that shows enhancement of elements up to Si and normal abundance of Ca and Fe-group elements. Whereas the C/N ratio is different among these stars, the (C+N)/O ratio is similar. This suggests that C was also yielded with similar abundance ratios relative to O-Zn in progenitors, and was later affected by the CN-cycle. By contrast, the heavy neutron-capture elements Sr and Ba are deficient in J2217+2104, compared to the four objects in this class previously studied. This indicates that the neutron-capture process in the early Galaxy, presumably the r-process, has no direct connection to the phenomenon that has formed such CEMP stars. Comparisons of the abundance pattern well determined for such CEMP stars with those of supernova nucleosynthesis models constrain the progenitor mass to be about 25Msun, which is not particularly different from typical mass of progenitors expected for extremely metal-poor stars in general.

We present the observations of an extreme ultraviolet (EUV) wave, a quasi-periodic fast-propagating (QFP) magnetosonic wave, and a kink wave that were simultaneously associated with the impingement of a coronal jet upon a group of coronal loops. After the interaction, the coronal loop showed obvious kink oscillation that had a period of about 428 seconds. In the meantime, a large-scale EUV wave and a QFP wave are observed on the west of the interaction position. It is interesting that the QFP wave showed refraction effect during its passing through two strong magnetic regions. The angular extent, speed, and lifetime of the EUV (QFP) wave were about 140 (40) degree, 423 (322) km/s, and 6 (26) minutes, respectively. It is measured that the period of the QFP wave was about 390 +/- 100 second. Based on the observational analysis results, we propose that the kink wave was probably excited by the interaction of the jet; the EUV was probably launched by the sudden expansion of the loop system due to the impingement of the coronal jet; and the QFP wave was possibly formed through the dispersive evolution of the disturbance caused by the jet-loop interaction.

By inverting the distributions of galaxies' apparent ellipticities and misalignment angles (measured around the projected half-light radius $R_{\rm e}$) between their photometric and kinematic axes, we study the intrinsic shape distribution of 189 slow rotator early-type galaxies with stellar masses $2\times 10^{11} M_{\odot}<M_\ast<2\times 10^{12} M_{\odot}$, extracted from a sample of about 2200 galaxies with integral-field stellar kinematics from the DR14 of the SDSS-IV MaNGA IFU survey. Thanks to the large sample of slow rotators, Graham+18 showed that there is clear structure in the misalignment angle distribution, with two peaks at both $0^{\circ}$ and $90^{\circ}$ misalignment (characteristic of oblate and prolate rotation respectively). Here we invert the observed distribution from Graham+18. The large sample allows us to go beyond the known fact that slow rotators are weakly triaxial and to place useful constraints on their intrinsic triaxiality distribution (around $1R_{\rm e}$) for the first time. The shape inversion is generally non-unique. However, we find that, for a wide set of model assumptions, the observed distribution clearly requires a dominant triaxial-oblate population. For some of our models, the data suggest a hint for a minor triaxial-prolate population, but a dominant prolate population is ruled out.

We investigate the possible connection between the most massive globular clusters, such as $\omega$ Cen and M54, and nuclear star clusters of dwarf galaxies that exhibit similar spreads in age and metallicity. We examine galactic nuclei in cosmological galaxy formation simulations at $z \approx 1.5$ to explore whether their age and metallicity spreads could explain these massive globular clusters. We derive structural properties of these nuclear regions, including mass, size, rotation, and shape. By using theoretical supernova yields to model the supernova enrichment in the simulations, we obtain individual elemental abundances for Fe, O, Na, Mg, and Al. Our nuclei are systematically more metal-rich than their host galaxies, which lie on the expected mass-metallicity relation. Some nuclei have a spread in Fe and age comparable to the massive globular clusters of the Milky Way, lending support to the hypothesis that nuclear star clusters of dwarf galaxies could be the progenitors of these objects. None of our nuclear regions contain the light element abundance spreads that characterize globular clusters, even when a large age spread is present. Our results demonstrate that extended star formation history within clusters, with metal pollution provided solely by supernova ejecta, is capable of replicating the metallicity spreads of massive globular clusters, but still requires another polluter to produce the light element variations.

Based on the first Gaia data release and spectroscopy from the LAMOST Data Release 4, we study the kinematics and chemistry of the local halo stars. The halo stars are identified kinematically with a relative speed of at least 220 km s$^{-1}$ with respect to the local standard of rest. In total, 436 halo stars are identified. From this halo sample, 16 high velocity (HiVel) stars are identified. We studied the metallicity and [$\alpha$/Fe] distribution of these HiVel stars. Though most of HiVel stars are metal-poor, there are several stars that have metallicity above $-0.5$ dex. To understand the origin of high velocity stars, we evolve the trajectory of the star backwards along the orbit in our adopted Galaxy potential model to determine the orbital parameters and assess whether the star could have originated in the Galactic center. We found that some high velocity stars could have originated from the Galactic center, but other stars were probably kicked up from the Galactic disk.

This study presents a study of equilibrium points, periodic orbits, stabilities, and manifolds in a rotating plane symmetric potential field. It has been found that the dynamical behaviour near equilibrium points is completely determined by the structure of the submanifolds and subspaces. The non-degenerate equilibrium points are classified into twelve cases. The necessary and sufficient conditions for linearly stable, non resonant unstable and resonant equilibrium points are established. Furthermore, the results show that a resonant equilibrium point is a Hopf bifurcation point. In addition, if the rotating speed changes, two non degenerate equilibria may collide and annihilate each other. The theory developed here is lastly applied to two particular cases, motions around a rotating, homogeneous cube and the asteroid 1620 Geographos. We found that the mutual annihilation of equilibrium points occurs as the rotating speed increases, and then the first surface shedding begins near the intersection point of the x axis and the surface. The results can be applied to planetary science, including the birth and evolution of the minor bodies in the Solar system, the rotational breakup and surface mass shedding of asteroids, etc.

The implementation of star formation and stellar feedback in cosmological simulations plays a critical role in shaping galaxy properties. In the first paper of the series, we presented a new method to model star formation as a collection of star clusters. In this paper, we improve the algorithm by eliminating accretion gaps, boosting momentum feedback, and introducing a subgrid initial bound fraction, $f_i$, that distinguishes cluster mass from stellar particle mass. We perform a suite of simulations with different star formation efficiency per freefall time $\epsilon_{\rm ff}$ and supernova momentum feedback intensity $f_{\rm boost}$. We find that the star formation history of a Milky Way-sized galaxy is sensitive to $f_{\rm boost}$, which allows us to constrain its value, $f_{\rm boost}\approx5$, in the current simulation setup. Changing $\epsilon_{\rm ff}$ from a few percent to 200\% has little effect on global galaxy properties. However, on smaller scales, the properties of star clusters are very sensitive to $\epsilon_{\rm ff}$. We find that $f_i$ increases with $\epsilon_{\rm ff}$ and cluster mass. Through the dependence on $f_i$, the shape of the cluster initial mass function varies strongly with $\epsilon_{\rm ff}$. The fraction of clustered star formation and maximum cluster mass increase with the star formation rate surface density, with the normalization of both relations dependent on $\epsilon_{\rm ff}$. The cluster formation timescale systematically decreases with increasing $\epsilon_{\rm ff}$. Local variations in the gas accretion history lead to a 0.25~dex scatter for the integral cluster formation efficiency. Joint constraints from all the observables prefer the runs that produce a median integral efficiency of 16\%.

Magnetic field is one of the key agents that play a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained due to the limitations in observations. We study the magnetic field in the massive infrared dark cloud G035.39-00.33 from dust continuum polarization observations at 850 $\micron$ with SCUBA-2/POL-2 at JCMT. The magnetic field tends to be perpendicular to the densest part of the main filament (F$_{M}$), whereas it has a less defined relative orientation in the rest of the structure, where it tends to be parallel to some diffuse regions. A mean plane-of-the-sky magnetic field strength of $\sim$50 $\mu$G for F$_{M}$ is obtained using Davis-Chandrasekhar-Fermi method. Based on $^{13}$CO (1-0) line observations, we suggest a formation scenario of F$_{M}$ due to large-scale ($\sim$10 pc) cloud-cloud collision. Using additional NH$_3$ line data, we estimate that F$_{M}$ will be gravitationally unstable if it is only supported by thermal pressure and turbulence. The northern part of F$_{M}$, however, can be stabilized by a modest additional support from the local magnetic field. The middle and southern parts of F$_{M}$ are likely unstable even if the magnetic field support is taken into account. We claim that the clumps in F$_{M}$ may be supported by turbulence and magnetic fields against gravitational collapse. Finally, we identified for the first time a massive ($\sim$200 M$_{\sun}$), collapsing starless clump candidate, "c8", in G035.39-00.33. The magnetic field surrounding "c8" is likely pinched, hinting at an accretion flow along the filament.

Using high spatial and temporal data from the \emph{Solar Dynamics Observatory} (\emph{SDO}) and the \emph{Interface Region Imaging Spectrograph} (\emph{IRIS}), several observational signatures of magnetic reconnection in the course of magnetic flux cancellation are presented, including two loop-loop interaction processes, multiple plasma blob ejections, and a sheet-like structure that appeared above the flux cancellation sites with a Y-shaped and an inverted Y-shaped ends. The \emph{IRIS} 1400 \AA\ observations show that the plasma blobs were ejected from the tip of the Y-shaped ends of the sheet-like structure. Obvious photospheric magnetic flux cancellation occurred after the first loop-loop interaction and continued until the end of the observation. Complemented by the nonlinear force-free field extrapolation, we found that two sets of magnetic field lines, which reveal an X-shaped configuration, align well with the interacted coronal loops. Moreover, a magnetic null point is found to be situated at about $0.9$ Mm height right above the flux cancellation sites and located between the two sets of magnetic field lines. These results suggest that the flux cancellation might be a result of submergence of magnetic field lines following magnetic reconnection that occurs in the lower atmosphere of the Sun, and the ejected plasma blobs should be plasmoids created in the sheet-like structure due to the tearing-mode instability. This observation reveals detailed magnetic field structure and dynamic process above the flux cancellation sites and will help us to understand magnetic reconnection in the lower atmosphere of the Sun.

The aim of this paper is to probe the features of the bouncing cosmology with the current observational data. Basing on bounce inflation model, with high derivative term, we propose a general parametrization of primordial power spectrum which includes the typical bouncing parameters, such as bouncing time-scale, and energy scale. By applying Markov Chain Monto Carlo analysis with current data combination of Planck 2015, BAO and JLA, we report the posterior probability distributions of the parameters. We find that, bouncing models can well explain CMB observations, especially the deficit and oscillation on large scale in TT power spectrum.

In the present paper, we investigate the dark energy equation of state using the Gaussian processes analysis method, without confining a particular parametrization. The reconstruction is carried out by adopting the background data including supernova and Hubble parameter, and perturbation data from the growth rate. It suggests that the background and perturbation data both present a hint of dynamical dark energy. However, the perturbation data have a more promising potential to distinguish non-evolution dark energy including the cosmological constant model. We also test the influence of some parameters on the reconstruction. We find that the matter density parameter $\Omega_{m0}$ has a slight effect on the background data reconstruction, but has a notable influence on the perturbation data reconstruction. While the Hubble constant presents a significant influence on the reconstruction from background data.

We examine the possible influence of Earth-directed coronal mass ejections (ECMEs) on the Sun's shadow in the 3~TeV cosmic-ray intensity observed by the Tibet-III air shower (AS) array. We confirm a clear solar-cycle variation of the intensity deficit in the Sun's shadow during ten years between 2000 and 2009. This solar-cycle variation is overall reproduced by our Monte Carlo (MC) simulations of the Sun's shadow based on the potential field model of the solar magnetic field averaged over each solar rotation period. We find, however, that the magnitude of the observed intensity deficit in the Sun's shadow is significantly less than that predicted by MC simulations, particularly during the period around solar maximum when a significant number of ECMEs is recorded. The $\chi^2$ tests of the agreement between the observations and the MC simulations show that the difference is larger during the periods when the ECMEs occur, and the difference is reduced if the periods of ECMEs are excluded from the analysis. This suggests the first experimental evidence of the ECMEs affecting the Sun's shadow observed in the 3~TeV cosmic-ray intensity.

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

POLAR is a compact wide-field space-borne detector for precise measurements of the linear polarisation of hard X-rays emitted by transient sources in the energy range from 50 keV to 500 keV. It consists of a 40$\times$40 array of plastic scintillator bars used as a detection material. The bars are grouped in 25 detector modules. The energy range sensitivity of POLAR is optimized to match with the prompt emission photons from the gamma-ray bursts (GRBs). Polarization measurements of the prompt emission would probe source geometries, emission mechanisms and magnetic structures in GRB jets. The instrument can also detect hard X-rays from solar flares and be used for precise measurement of their polarisation. POLAR was launched into a low Earth orbit on-board the Chinese space-lab TG-2 on September 15th, 2016. To achieve high accuracies in polarisation measurements it is essential to assure both before and during the flight a precise energy calibration. Such calibrations are performed with four low activity $^{22}$Na radioactive sources placed inside the instrument. Energy conversion factors are related to Compton edge positions from the collinear annihilation photons from the sources. This paper presents main principles of the in-flight calibration, describes studies of the method based on Monte Carlo simulations and its laboratory verification and finally provides some observation results based on the in-flight data analysis.

We demonstrate the dragging of the magnetic field by the super-Alfvenic shear flows out of the reconnection plane can strongly localize the reconnection x-line in collisionless plasmas, reversing the current direction at the x-line. Reconnection with this new morphology, which is impossible in resistive-magnetohydrodynamic (MHD), is enabled by electron inertia. Surprisingly, the quasi-steady reconnection rate remains of order 0.1 even though the aspect ratio of the local x-line geometry is larger than unity. We explain this by examining the transport of the reconnected magnetic flux and the opening angle made by the upstream magnetic field, concluding that the reconnection rate is still limited by the constraint imposed at the inflow region. This study further suggests the nearly universal fast rate value of order 0.1 cannot be explained by the physics of tearing modes, nor can it be explained by a universal localization mechanism.

We are interested in stable periodic orbits for spacecrafts in the gravitational field of minor celestial bodies. The stable periodic orbits around minor celestial bodies are useful not only for the mission design of the deep space exploration, but also for studying the long-time stability of small satellites in the large-size-ratio binary asteroids. The irregular shapes and gravitational fields of the minor celestial bodies are modeled by the polyhedral model. Using the topological classifications of periodic orbits and the grid search method, the stable periodic orbits can be calculated and the topological cases can be determined. Furthermore, we find five different types of stable periodic orbits around minor celestial bodies: A) Stable periodic orbits generated from the stable equilibrium points outside the minor celestial body, B) Stable periodic orbits continued from the unstable periodic orbits around the unstable equilibrium points, C) Retrograde and nearly circular periodic orbits with zero-inclination around minor celestial bodies, D) Resonance periodic orbits, E) Near-surface inclined periodic orbits. We take asteroid 243 Ida, 433 Eros, 6489 Golevka, 101955 Bennu, and the comet 1P/Halley for examples.

POLAR is a compact space-borne detector designed to perform reliable measurements of the polarization for transient sources like Gamma-Ray Bursts in the energy range 50-500keV. The instrument works based on the Compton Scattering principle with the plastic scintillators as the main detection material along with the multi-anode photomultiplier tube. POLAR has been launched successfully onboard the Chinese space laboratory TG-2 on 15th September, 2016. In order to reliably reconstruct the polarization information a highly detailed understanding of the instrument is required for both data analysis and Monte Carlo studies. For this purpose a full study of the in-orbit performance was performed in order to obtain the instrument calibration parameters such as noise, pedestal, gain nonlinearity of the electronics, threshold, crosstalk and gain, as well as the effect of temperature on the above parameters. Furthermore the relationship between gain and high voltage of the multi-anode photomultiplier tube has been studied and the errors on all measurement values are presented. Finally the typical systematic error on polarization measurements of Gamma-Ray Bursts due to the measurement error of the calibration parameters are estimated using Monte Carlo simulations.

This paper reports on the observation of the sidereal large-scale anisotropy of cosmic rays using data collected by the ARGO-YBJ experiment over 5 years (2008$-$2012). This analysis extends previous work limited to the period from 2008 January to 2009 December,near the minimum of solar activity between cycles 23 and 24.With the new data sample the period of solar cycle 24 from near minimum to maximum is investigated. A new method is used to improve the energy reconstruction, allowing us to cover a much wider energy range, from 4 to 520 TeV. Below 100 TeV, the anisotropy is dominated by two wide regions, the so-called "tail-in" and "loss-cone" features. At higher energies, a dramatic change of the morphology is confirmed. The yearly time dependence of the anisotropy is investigated. Finally, no noticeable variation of cosmic-ray anisotropy with solar activity is observed for a median energy of 7 TeV.

We report on a broadband study of a complex X-ray source (1SAX J0618.0+2227) associated with the interaction site of the supernova remnant (SNR) IC 443 and ambient molecular cloud (MC) using NuSTAR, XMM_Newton, and Chandra observations. Its X-ray spectrum is composed of both thermal and non-thermal components. The thermal component can be equally well represented by either a thin plasma model with kT=0.19 keV or a blackbody model with kT=0.11 keV. The non-thermal component can be fit with either a power-law with Gamma~1.7 or a cutoff power-law with Gamma~1.5 and a cutoff energy at E_cut~18 keV. Using the newly obtained NuSTAR dataset, we test three possible scenarios for isolated X-ray sources in the SNR-MC interaction site: 1) pulsar wind nebula (PWN); 2) SNR ejecta fragment; 3) shocked molecular clump. We conclude that this source is most likely composed of a SNR ejecta (or a PWN) and surrounding shocked molecular clumps. The nature of this hard X-ray source in the SNR-MC interaction site of IC 443 may shed light on unidentified X-ray sources with hard X-ray spectra in rich environments for star forming regions, such as the Galactic center.

Using mock spectra based on Vazdekis/MILES library fitted within the wavelength region 3600-7350\AA, we analyze the bias and scatter on the resulting physical parameters induced by the choice of fitting algorithms and observational uncertainties, but avoid effects of those model uncertainties. We consider two full-spectrum fitting codes: pPXF and STARLIGHT, in fitting for stellar population age, metallicity, mass-to-light ratio, and dust extinction. With pPXF we find that both the bias in the population parameters and the scatter in the recovered logarithmic values follows the expected trend. The bias increases for younger ages and systematically makes recovered ages older, $M_*/L_r$ larger and metallicities lower than the true values. For reference, at S/N=30, and for the worst case ($t=10^8$yr), the bias is 0.06 dex in $M_*/L_r$, 0.03 dex in both age and [M/H]. There is no significant dependence on either E(B-V) or the shape of the error spectrum. Moreover, the results are consistent for both our 1-SSP and 2-SSP tests. With the STARLIGHT algorithm, we find trends similar to pPXF, when the input E(B-V)<0.2 mag. However, with larger input E(B-V), the biases of the output parameter do not converge to zero even at the highest S/N and are strongly affected by the shape of the error spectra. This effect is particularly dramatic for youngest age, for which all population parameters can be strongly different from the input values, with significantly underestimated dust extinction and [M/H], and larger ages and $M_*/L_r$. Results degrade when moving from our 1-SSP to the 2-SSP tests. The STARLIGHT convergence to the true values can be improved by increasing Markov Chains and annealing loops to the "slow mode". For the same input spectrum, pPXF is about two order of magnitudes faster than STARLIGHT's "default mode" and about three order of magnitude faster than STARLIGHT's "slow mode".

The abundance scatter of heavier r-process elements (Z > 56) relative to Fe ([r/Fe]) in metal-poor stars preserves excellent information of the star formation history and provides important insights into the various situations of the Galactic chemical enrichment. In this respect, the upper and lower boundaries of [r/Fe] could present useful clues for investigating the extreme situations of the star formation history and the early Galactic chemical evolution. In this paper, we investigate the formation of the upper and lower boundaries of [r/Fe] for the gas clouds. We find that, for a cloud from which metal-poor stars formed, the formation of the upper limits of [r/Fe] is mainly due to the pollution from a single main r-process event. For a cloud from which metal-poor stars formed, the formation of the lower limits of [r/Fe] is mainly due to the pollution from a single SN II event that ejects primary Fe.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

The Beijing Faint Object Spectrograph and Camera (BFOSC) is one of the most important instruments of the 2.16-m telescope of the Xinglong Observatory. Every year there are ~ 20 SCI-papers published based on the observational data of this telescope. In this work, we have systemically measured the total efficiency of the BFOSC of the 2.16-m reflector, based on the observations of two ESO flux standard stars. We have obtained the total efficiencies of the BFOSC instrument of different grisms with various slit widths in almost all ranges, and analysed the factors which effect the efficiency of telescope and spectrograph. For the astronomical observers, the result will be useful for them to select a suitable slit width, depending on their scientific goals and weather conditions during the observation; For the technicians, the result will help them systemically find out the real efficiency of telescope and spectrograph, and further to improve the total efficiency and observing capacity of the telescope technically.

We present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

We present our optical photometric and spectroscopical observations of GRB 140629A. A redshift of $z=2.275\pm0.043$ is measured through the metal absorption lines in our spectroscopic data. Using our photometric data and multiple observational data observed with other telescopes, we show that its optical light curve is well interpreted with the standard forward shock models in the thin shell case. Its optical-X-ray afterglow spectrum is jointly fitted with a single power-law function, yielding a photon index of $-1.90\pm 0.05$. The optical extinction and neutral hydrogen absorption of the GRB host galaxy are negligible. The fit to the light curve with the standard models shows that the ambient density is $60\pm 9$ cm$^{-3}$ and the GRB radiating efficiency is as low as $\sim 0.24\%$, likely indicating a baryonic-dominated ejecta of this GRB. This burst agrees well with the $L_{\rm p, iso}-E_p^{'}-\Gamma_0$ relation, but confidently violates those empirical relations involving geometric corrections (or jet break time). This gives rise to an issue of possible selection effect on these relations since the jet opening angle of this GRB is extremely narrow (0.04 rad).

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.

In this Letter, we investigate the long-duration quasi-static evolution of 12 pre-eruptive filaments (4 active region and 8 quiescent filaments), mainly focusing on the evolution of the filament height in three dimension (3D) and the decay index of the background magnetic field. The filament height in 3D is derived through two-perspective observations of \textit{Solar Dynamics Observatory} and \textit{Solar TErrestrial RElations Observatory}. The coronal magnetic field is reconstructed using the potential field source surface model. A new finding is that the filaments we studied show two types of long-duration evolution: one type is comprised of a long-duration static phase and a short slow rise phase with a duration of less than 12 hours and a speed of 0.1--0.7 km s$^{-1}$, while the other one only presents a slow rise phase but with an extremely long duration of more than 60 hours and a smaller speed of 0.01--0.2 km s$^{-1}$. At the moment approaching the eruption, the decay index of the background magnetic field at the filament height is similar for both active region and quiescent filaments. The average value and upper limit are $\sim$0.9 and $\sim$1.4, close to the critical index of torus instability. Moreover, the filament height and background magnetic field strength are also found to be linearly and exponentially related with the filament length, respectively.

In this paper, we show how to carry out a relatively more realistic and complete reconstruction of supernova neutrino spectra in the future large liquid-scintillator detectors, by implementing the method of singular value decomposition with a proper regularization. For a core-collapse supernova at a distance of $10~{\rm kpc}$ in the Milky Way, its $\overline{\nu}^{}_e$ spectrum can be precisely determined from the inverse beta-decay process $\overline{\nu}^{}_e + p \to e^+ + n$, for which a $20~{\rm kiloton}$ liquid-scintillator detector with the resolution similar to the Jiangmen Underground Neutrino Observatory (JUNO) may register more than 5000 events. We have to rely predominantly on the elastic neutrino-electron scattering $\nu + e^- \to \nu + e^-$ and the elastic neutrino-proton scattering $\nu + p \to \nu + p$ for the spectra of $\nu^{}_e$ and $\nu^{}_x$, where $\nu$ denotes collectively neutrinos and antineutrinos of all three flavors and $\nu^{}_x$ for $\nu^{}_\mu$ and $\nu^{}_\tau$ as well as their antiparticles. To demonstrate the validity of our approach, we also attempt to reconstruct the neutrino spectra by using the time-integrated neutrino data from the latest numerical simulations of delayed neutrino-driven supernova explosions.

We investigate the orbital stability close to the unique L4-point Jupiter binary Trojan asteroid 624 Hektor. The gravitational potential of 624 Hektor is calculated using the polyhedron model with observational data of 2038 faces and 1021 vertexes. Previous studies have presented three different density values for 624 Hektor. The equilibrium points in the gravitational potential of 624 Hektor with different density values have been studied in detail. There are five equilibrium points in the gravitational potential of 624 Hektor no matter the density value. The positions, Jacobian, eigenvalues, topological cases, stability, as well as the Hessian matrix of the equilibrium points are investigated. For the three different density values the number, topological cases, and the stability of the equilibrium points with different density values are the same. However, the positions of the equilibrium points vary with the density value of the asteroid 624 Hektor. The outer equilibrium points move away from the asteroid s mass center when the density increases, and the inner equilibrium point moves close to the asteroid s mass center when the density increases. There exist unstable periodic orbits near the surface of 624 Hektor. We calculated an orbit near the primary s equatorial plane of this binary Trojan asteroid, the results indicate that the orbit remains stable after 28.8375 d.

Recent Atacama Large Millimeter and Submillimeter Array (ALMA) observations of the protoplanetary disk around the Herbig Ae star HD 163296 revealed three depleted dust gaps at 60, 100 and 160 au in the 1.3 mm continuum as well as CO depletion in the middle and outer dust gaps. However, no CO depletion was found in the inner dust gap. To examine the planet--disk interaction model, we present results of two-dimensional two fluid (gas + dust) hydrodynamic simulations coupled with three-dimensional radiative transfer simulations. In order to fit the high gas-to-dust ratio of the first gap, we find the Shakura--Sunyaev viscosity parameter $\alpha$ must be very small ($\lesssim 10^{-4}$) in the inner disk. On the other hand, a relatively large $\alpha$ ($\sim 7.5\times 10^{-3}$) is required to reproduce the dust surface density in the outer disk. We interpret the variation of $\alpha$ as an indicator of the transition from an inner dead zone to the outer magnetorotational instability (MRI) active zone. Within $\sim 100$ au, the HD 163296 disk's ionization level is low, and non-ideal magnetohydrodynamic (MHD) effects could suppress the MRI, so the disk can be largely laminar. The disk's ionization level gradually increases toward larger radii, and the outermost disk ($r > 300$ au) becomes turbulent due to MRI. Under this condition, we find that the observed dust continuum and CO gas line emissions can be reasonably fit by three half-Jovian-mass planets (0.46, 0.46 and 0.58 $M_\textrm{J}$) at 59, 105 and 160 au, respectively.

Particle acceleration in space and astrophysical reconnection sites is an important unsolved problem in studies of magnetic reconnection. Earlier kinetic simulations have identified several acceleration mechanisms that are associated with particle drift motions. Here, we show that, for sufficient large systems, the energization processes due to particle drift motions can be described as fluid compression and shear, and that the shear energization is proportional to the pressure anisotropy of energetic particles. By analyzing results from fully kinetic simulations, we show that the compression energization dominates the acceleration of high-energy particles in reconnection with a weak guide field, and the compression and shear effects are comparable when the guide field is 50\% of the reconnecting component. Spatial distributions of those energization effects reveal that reconnection exhausts, contracting islands, and island merging regions are the three most important regions for compression and shear acceleration. This study connects particle energization by particle guiding-center drift motions with that due to background fluid motions, as in the energetic particle transport theory. It provides foundations for building particle transport models for large-scale reconnection acceleration such as those in solar flares.

The discovery of thousands of exoplanets over the last couple of decades has shown that the birth of planets is a very efficient process in nature. Theories invoke a multitude of mechanisms to describe the assembly of planets in the disks around pre-main-sequence stars, but observational constraints have been sparse on account of insufficient sensitivity and resolution. Understanding how planets form and interact with their parental disk is crucial also to illuminate the main characteristics of a large portion of the full population of planets that is inaccessible to current and near-future observations. This White Paper describes some of the main issues for our current understanding of the formation and evolution of planets, and the critical contribution expected in this field by the Next Generation Very Large Array.

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

We present the observed H$\alpha$ flux and derived star formation rates (SFRs) for a fall sample of low$-$surface$-$brightness galaxies (LSBGs). The sample is selected from the fall sky region of the 40$\%$ ALFALFA {\ion{H}{1}} survey $-$ SDSS DR7 photometric data, and all the $H\alpha$ images were obtained using the 2.16 m telescope, operated by the National Astronomy Observatories, Chinese Academy of Sciences. A total of 111 LSBGs were observed and $H\alpha$ flux was measured in 92 of them. Though almost all the LSBGs in our sample are {\ion{H}{1}}$-$rich, their SFRs derived from the extinction and filter$-$transmission$-$corrected $H\alpha$ flux, are less than 1$M_{\sun}$$yr^{-1}$. LSBGs and star forming galaxies have similar {\ion{H}{1}} surface densities, but LSBGs have much lower SFRs and SFR surface densities than star$-$forming galaxies. Our results show that LSBGs deviate from the Kennicutt-Schmidt law significantly, which indicate that they have low star formation efficiency. The SFRs of LSBGs are close to average SFRs in Hubble time and support the previous arguments that most of the LSBGs are stable systems and they tend to seldom contain strong interactions or major mergers during their star formation histories.

We study the radial acceleration relation (RAR) for early-type galaxies (ETGs) in the SDSS MaNGA MPL5 dataset. The complete ETG sample show a slightly offset RAR from the relation reported by McGaugh et al. (2016) at the low-acceleration end; we find that the deviation is due to the fact that the slow rotators show a systematically higher acceleration relation than the McGaugh's RAR, while the fast rotators show a consistent acceleration relation to McGaugh's RAR. There is a 1\sigma significant difference between the acceleration relations of the fast and slow rotators, suggesting that the acceleration relation correlates with the galactic spins, and that the slow rotators may have a different mass distribution compared with fast rotators and late-type galaxies. We suspect that the acceleration relation deviation of slow rotators may be attributed to more galaxy merger events, which would disrupt the original spins and correlated distributions of baryons and dark matter orbits in galaxies.

Rapid response and short time latency are very important for Time Domain Astronomy, such as the observations of Gamma-ray Bursts (GRBs) and electromagnetic (EM) counterparts of gravitational waves (GWs). Based on the near real-time Fermi/GBM data, we developed a low-latency pipeline to automatically calculate the temporal and spectral properties of GRBs. With this pipeline, some important parameters can be obtained, such as T90 and fluence, within ~20 minutes after the GRB trigger. For ~90% GRBs, T90 and fluence are consistent with the GBM catalog results within 2 sigma errors. This pipeline has been used by the Gamma-ray Bursts Polarimeter (POLAR) and the Insight Hard X-ray Modulation Telescope (Insight-HXMT) to follow up the bursts of interest. For GRB 170817A, the first EM counterpart of GW events detected by Fermi/GBM and INTEGRAL/SPI-ACS, the pipeline gave T90 and spectral information in 21 minutes after the GBM trigger, providing important information for POLAR and Insight-HXMT observations.

As an essential part of China’s Gravitational Waves Program, the Ali CMB Polarization Telescope (AliCPT) is a ground-based experiment aiming at the Primordial Gravitational Waves (PGWs) by measuring B-mode polarization of Cosmic Microwave Background (CMB). First proposed in 2014 and currently in fast construction phase, AliCPT is China’s first CMB project that plans for commissioning in 2019. Led by the Institute of High Energy Physics (IHEP) under the Chinese Academy of Sciences (CAS), the project is a worldwide collaboration of more than fifteen universities and research institutes. Ali CMB Project is briefly introduced.

We measure $\lambda_{R_e}$, a proxy for galaxy specific stellar angular momentum within one effective radius, and the ellipticity, $\epsilon$, for about 2300 galaxies of all morphological types observed with integral field spectroscopy as part of the MaNGA survey, the largest such sample to date. We use the $(\lambda_{R_e}, \epsilon)$ diagram to separate early-type galaxies into fast and slow rotators. We also visually classify each galaxy according to its optical morphology and two-dimensional stellar velocity field. Comparing these classifications to quantitative $\lambda_{R_e}$ measurements reveals tight relationships between angular momentum and galaxy structure. In order to account for atmospheric seeing, we use realistic models of galaxy kinematics to derive a general approximate analytic correction for $\lambda_{R_e}$. Thanks to the size of the sample and the large number of massive galaxies, we unambiguously detect a clear bimodality in the $(\lambda_{R_e}, \epsilon)$ diagram which may result from fundamental differences in galaxy assembly history. There is a sharp secondary density peak inside the region of the diagram with low $\lambda_{R_e}$ and $\epsilon < 0.4$, previously suggested as the definition for slow rotators. Most of these galaxies are visually classified as non-regular rotators and have high velocity dispersion. The intrinsic bimodality must be stronger, as it tends to be smoothed by noise and inclination. The large sample of slow rotators allows us for the first time to unveil a secondary peak at +/-90 degrees in their distribution of the misalignments between the photometric and kinematic position angles. We confirm that genuine slow rotators start appearing above a stellar mass of 2\times10^{11} M_{\odot}$ where a significant number of high-mass fast rotators also exist.

A correlation between the secondary cosmic ray flux and the near-earth electric field intensity, measured during thunderstorms, has been found by analyzing the data of the ARGO-YBJ experiment, a full coverage air shower array located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China). The counting rates of showers with different particle multiplicities, have been found to be strongly dependent upon the intensity and polarity of the electric field measured during the course of 15 thunderstorms. In negative electric fields (i.e. accelerating negative charges downwards), the counting rates increase with increasing electric field strength. In positive fields, the rates decrease with field intensity until a certain value of the field EFmin (whose value depends on the event multiplicity), above which the rates begin increasing. By using Monte Carlo simulations, we found that this peculiar behavior can be well described by the presence of an electric field in a layer of thickness of a few hundred meters in the atmosphere above the detector, which accelerates/decelerates the secondary shower particles of opposite charge, modifying the number of particles with energy exceeding the detector threshold. These results, for the first time, give a consistent explanation for the origin of the variation of the electron/positron flux observed for decades by high altitude cosmic ray detectors during thunderstorms.

In this paper, we will give a general introduction to the project of Ali CMB Polarization Telescope (AliCPT), which is a Sino-US joint project led by the Institute of High Energy Physics (IHEP) and has involved many different institutes in China. It is the first ground-based Cosmic Microwave Background (CMB) polarization experiment in China and an integral part of China's Gravitational Waves Program. The main scientific goal of AliCPT project is to probe the primordial gravitational waves (PGWs) originated from the very early Universe. The AliCPT project includes two stages. The first stage referred to as AliCPT-1, is to build a telescope in the Ali region of Tibet with an altitude of 5,250 meters. Once completed, it will be the worldwide highest ground-based CMB observatory and open a new window for probing PGWs in northern hemisphere. AliCPT-1 telescope is designed to have about 7,000 TES detectors at 90GHz and 150GHz. The second stage is to have a more sensitive telescope (AliCPT-2) with the number of detectors more than 20,000. Our simulations show that AliCPT will improve the current constraint on the tensor-to-scalar ratio $r$ by one order of magnitude with 3 years' observation. Besides the PGWs, the AliCPT will also enable a precise measurement on the CMB rotation angle and provide a precise test on the CPT symmetry. We show 3 years' observation will improve the current limit by two order of magnitude.

We perform full spectrum fitting stellar population analysis and Jeans Anisotropic modelling (JAM) of the stellar kinematics for about 2000 early-type galaxies (ETGs) and spiral galaxies from the MaNGA DR14 sample. Galaxies with different morphologies are found to be located on a remarkably tight mass plane which is close to the prediction of the virial theorem, extending previous results for ETGs. By examining an inclined projection (`the mass-size' plane), we find that spiral and early-type galaxies occupy different regions on the plane, and their stellar population properties (i.e. age, metallicity and stellar mass-to-light ratio) vary systematically along roughly the direction of velocity dispersion, which is a proxy for the bulge fraction. Galaxies with higher velocity dispersions have typically older ages, larger stellar mass-to-light ratios and are more metal rich, which indicates that galaxies increase their bulge fractions as their stellar populations age and become enriched chemically. The age and stellar mass-to-light ratio gradients for low-mass galaxies in our sample tend to be positive ($\rm centre<outer$), while the gradients for most massive galaxies are negative. The metallicity gradients show a clear peak around velocity dispersion $\log_{10} \sigma_{\rm e}\approx 2.0$, which corresponds to the critical mass $\sim 3\times 10^{10}M_{\odot}$ of the break in the mass-size relation. Spiral galaxies with large mass and size have the steepest gradients, while the most massive ETGs, especially above the critical mass $M_{\rm crit}\ge 2\times 10^{11} M_{\odot}$, where slow rotator ETGs start dominating, have much flatter gradients. This may be due to differences in their evolution histories, e.g. mergers.

We analyze the Sun's shadow observed with the Tibet-III air shower array and find that the shadow's center deviates northward (southward) from the optical solar disc center in the "Away" ("Toward") IMF sector. By comparing with numerical simulations based on the solar magnetic field model, we find that the average IMF strength in the "Away" ("Toward") sector is $1.54 \pm 0.21_{\rm stat} \pm 0.20_{\rm syst}$ ($1.62 \pm 0.15_{\rm stat} \pm 0.22_{\rm syst}$) times larger than the model prediction. These demonstrate that the observed Sun's shadow is a useful tool for the quantitative evaluation of the average solar magnetic field.

We have studied the variability of \object{S5 0716+714} at radio 15 GHz and $\gamma$-ray band using three different methods. A possible periodicity of $P_{15 GHz}=266.0\pm11.5$ and $P_{\gamma}=344.0\pm16.4$ days are obtained for radio 15 GHz and $\gamma$-ray light curves, respectively. The variability may be related to the intrinsically emission mechanism. The difference between the variability timescales of radio 15 GHz and $\gamma$-ray may be due to that the emission of radio 15 GHz is produced via the synchrotron process, while the $\gamma$-ray is produced by both the SSC and EC processes.

Current sheet is believed to be the region of energy dissipation via magnetic reconnection in solar flares. However, its properties, for example, the dynamic process, have not been fully understood. Here we report a current sheet in a solar flare (SOL2017-09-10T16:06) that was clearly observed by the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory as well as the EUV Imaging Spectrometer on Hinode. The high-resolution imaging and spectroscopic observations show that the current sheet is mainly visible in high temperature (>10 MK) passbands, particularly in the Fe XXIV 192.03 line with a formation temperature of ~18 MK. The hot Fe XXIV 192.03 line exhibits very large nonthermal velocities up to 200 km/s in the current sheet, suggesting that turbulent motions exist there. The largest turbulent velocity occurs at the edge of the current sheet, with some offset with the strongest line intensity. At the central part of the current sheet, the turbulent velocity is negatively correlated with the line intensity. From the line emission and turbulent features we obtain a thickness in the range of 7--11 Mm for the current sheet. These results suggest that the current sheet has internal fine and dynamic structures that may help the magnetic reconnection within it proceeds efficiently.

We investigate the potential of the Atacama Large Millimeter/submillimeter Array (ALMA) and the Next Generation Very Large Array (ngVLA) to observe substructures in nearby young disks which are due to the gravitational interaction between disk material and planets close to the central star. We simulate the gas and dust dynamics in the disk using the LA-COMPASS hydrodynamical code. We generate synthetic images for the dust continuum emission at sub-millimeter to centimeter wavelengths and simulate ALMA and ngVLA observations. We explore the parameter space of some of the main disk and planet properties that would produce substructures that can be visible with ALMA and the ngVLA. We find that ngVLA observations with an angular resolution of 5 milliarcsec at 3 mm can reveal and characterize gaps and azimuthal asymmetries in disks hosting planets with masses down to $\approx~5 M_{\oplus}$ $\approx 1 - 5$ au from a Solar-like star in the closest star forming regions, whereas ALMA can detect gaps down to planetary masses of $\approx 20~M_{\oplus}$ at 5 au. Gaps opened by super-Earth planets with masses $\approx 5 - 10 M_{\oplus}$ are detectable by the ngVLA in the case of disks with low viscosity ($\alpha \sim 10^{-5}$) and low pressure scale height ($h \approx 0.025$ au at 5 au). The ngVLA can measure the proper motion of azimuthal asymmetric structures associated with the disk-planet interaction, as well as possible circumplanetary disks on timescales as short as one to a few weeks for planets at $1 - 5$ au from the star.

Li abundances in the bulk of low-mass metal-poor stars are well reproduced by stellar evolution models adopting a constant initial abundance. However, a small number of stars have exceptionally high Li abundances, for which no convincing models have been established. We report on the discovery of 12 very metal-poor stars that have large excesses of Li, including an object having more than 100 times higher Li abundance than the values found in usual objects, which is the the largest excess in metal-poor stars known to date. The sample is distributed over a wide range of evolutionary stages, including five unevolved stars, showing no abundance anomaly in other elements. The results indicate the existence of an efficient process to enrich Li in a small fraction of low-mass stars at the main-sequence or subgiant phase. The wide distribution of Li-rich stars along the red giant branch could be explained by dilution of surface Li by mixing that occurs when the stars evolve into red giants. Our study narrows down the problem to be solved to understand the origins of Li-excess found in low-mass stars, suggesting the presence of unknown process that affects the surface abundances preceding red giant phases.

In this work, we present a novel centroiding method based on Fourier space Phase Fitting(FPF) for Point Spread Function(PSF) reconstruction. We generate two sets of simulations to test our method. The first set is generated by GalSim with elliptical Moffat profile and strong anisotropy which shifts the center of the PSF. The second set of simulation is drawn from CFHT i band stellar imaging data. We find non-negligible anisotropy from CFHT stellar images, which leads to $\sim$0.08 scatter in unit of pixels using polynomial fitting method Vakili and Hogg (2016). And we apply FPF method to estimate the centroid in real space, this scatter reduces to $\sim$0.04 in SNR=200 CFHT like sample. In low SNR (50 and 100) CFHT like samples, the background noise dominates the shifting of the centroid, therefore the scatter estimated from different methods are similar. We compare polynomial fitting and FPF using GalSim simulation with optical anisotropy. We find that in all SNR$\sim$50, 100 and 200) samples, FPF performs better than polynomial fitting by a factor of $\sim$3. In general, we suggest that in real observations there are anisotropy which shift the centroid, and FPF method is a better way to accurately locate it.

The hybrid Tibet AS array was successfully constructed in 2014. It has 4500 m$^{2}$ underground water Cherenkov pools used as the muon detector (MD) and 789 scintillator detectors covering 36900 m$^{2}$ as the surface array. At 100 TeV, cosmic-ray background events can be rejected by approximately 99.99\%, according to the full Monte Carlo (MC) simulation for $\gamma$-ray observations. In order to use the muon detector efficiently, we propose to extend the surface array area to 72900 m$^{2}$ by adding 120 scintillator detectors around the current array to increase the effective detection area. A new prototype scintillator detector is developed via optimizing the detector geometry and its optical surface, by selecting the reflective material and adopting dynode readout. This detector can meet our physics requirements with a positional non-uniformity of the output charge within 10\% (with reference to the center of the scintillator), time resolution FWHM of $\sim$2.2 ns, and dynamic range from 1 to 500 minimum ionization particles.

Further studies on the corner singularity of GS reconstruction are are compiled in this paper. It's focused on solution of the Data Completion (DC) problem with the Extended Hilbert Transform (EHT) over plane rectangular region. Optimal selections of the regular parameters in{\it Tikhonov} solution of corresponding DC problem are developed in this study. The 6-parameter regular solutions and the {\it Jacobian} matrix and one {\it Hessian} tensor to the regular parameters are derived in this work. A concise formula for EHT in the near field of corners, which shows property for EHT near the corner, are also provided. It serves as the additional constraints for our parameter optimization problem (OP). Third, a nonlinear convex function defined by the regular solution and the corner constraints is introduced and is taken as the object function for the OP of the 6 regular parameters on half-space $(\mathbf{p}>0)$. Given an initial guess of $\mathbf{p}_0$, the optimal parameters are solved from the OP through a well known constrained nonlinear optimization method. Last, the benchmark tests to the proposed solution approach are carried out, and detailed results from totally 9 different bench-cases are tabulated. In contrast to solutions with given regular parameters, our bench results demonstrate that an objective way for selection of the optimal $\mathbf{p}$ is successfully laid out here. Robustness and efficiency of the suggested new approach are also highlighted in this study.

The Cosmic Microwave Background (CMB) Polarization plays an important role in current cosmological studies. CMB B-mode polarization is the most effective probe to primordial gravitational waves (PGWs) and a test of the inflation as well as other theories of the early universe such as bouncing and cyclic universe. So far, major ground-based CMB polarization experiments are located in the southern hemisphere.Recently, China has launched the Ali CMB Polarization Telescope (AliCPT) in Tibetan Plateau to measure CMB B mode polarization and detect the PGWs in northern hemisphere. AliCPT include two stages, the first one is to build a telescope at the 5250m site (AliCPT-1) and the second one is to have a more sensitive telescope at a higher altitude of about 6000m (AliCPT-2). In this paper, we report the atmospherical conditions, sky coverage and the current infrastructure associated with AliCPT. We analyzed the reanalysis data from MERRA-2 together with radiosonde data from the Ali Meteorological Service and found that the amount of water vapor has a heavy seasonal variation and October to March is the suitable observation time. We also found 95/150 GHz to be feasible for AliCPT-1 and higher frequencies to be possible for AliCPT-2. Then we analyzed the observable sky and the target fields, and showed that Ali provides us a unique opportunity to observe CMB with less foreground contamination in the northern hemisphere and is complementary to the existed southern CMB experiments. Together with the developed infrastructure, we point out that Ali opens a new window for CMB observation and will be one of the major sites in the world along with Antarctic and Atacama.

The unexpectedly high flux of cosmic ray positrons detected at Earth may originate from nearby astrophysical sources, dark matter, or unknown processes of cosmic-ray secondary production. We report the detection, using the HighAltitude Water Cherenkov Observatory (HAWC), of extended tera-electron volt gamma-ray emission coincident with the locations of two nearby middle-aged pulsars (Geminga and PSR B0656+14). The HAWC observations demonstrate that these pulsars are indeed local sources of accelerated leptons, but the measured tera-electron volt emission profile constrains the diffusion of particles away from these sources to be much slower than previously assumed. We demonstrate that the leptons emitted by these objects are therefore unlikely to be the origin of the excess positrons, which may have a more exotic origin.

In this paper, we make a deep analysis for the five typical interacting holographic dark energy models with the interaction terms $Q=3\beta H_{0}\rho_{\rm{de}}$, $Q=3\beta H_{0}\rho_{\rm{c}}$, $Q=3\beta H_{0}(\rho_{\rm{de}}+\rho_{\rm c})$, $Q=3\beta H_{0}\sqrt{\rho_{\rm{de}}\rho_{\rm c}}$, and $Q=3\beta H_{0}\frac{\rho_{\rm{de}}\rho_{c}}{\rho_{\rm{de}}+\rho_{\rm c}}$, respectively. We obtain observational constraints on these models by using the type Ia supernova data (the Joint Light-curve Analysis sample), the cosmic microwave background data (Planck 2015 distance priors), the baryon acoustic oscillations data, and the direct measurement of the Hubble constant. We find that the values of $\chi_{\rm min}^2$ for all the five models are almost equal (around~699), indicating that the current observational data equally favor these IHDE models. In addition, a comparison with the cases of interaction term involving the Hubble parameter $H$ is also made.

Recent observations have shown that coronal shocks driven by coronal mass ejections can develop and accelerate particles within several solar radii in large solar energetic particle (SEP) events. Motivated by this, we present an SEP acceleration study including the process that a fast shock propagates through a streamer-like magnetic field with both closed and open field lines in the low corona region. The acceleration of protons is modeled by numerically solving the Parker transport equation with spatial diffusion both along and across the magnetic field. We show that particles can be sufficiently accelerated to up to several hundred MeV within 2-3 solar radii. When the shock propagates through a streamer-like magnetic field, particles are more efficiently accelerated compared to the case with a simple radial magnetic field, mainly due to perpendicular shock geometry and the natural trapping effect of closed magnetic fields. Our results suggest that the coronal magnetic field configuration is an important factor for producing large SEP events. We further show that the coronal magnetic field configuration strongly influences the distribution of energetic particles, leading to different locations of source regions along the shock front where most of high energy particles are concentrated. This work may have strong implications to SEP observations. The upcoming Parker Solar Probe will provide in situ observations for the distribution of energetic particles in the coronal shock region, and test the results of the study.

In this paper we analyze the dynamical behavior of large dust grains in the vicinity of a cometary nucleus. To this end we consider the gravitational field of the irregularly shaped body, as well as its electric and magnetic fields. Without considering the effect of gas friction and solar radiation, we find that there exist grains which are static relative to the cometary nucleus; the positions of these grains are the stable equibria. There also exist grains in the stable periodic orbits close to the cometary nucleus. The grains in the stable equibria or the stable periodic orbits will not escape or impact on the surface of the cometary nucleus. The results are applicable for large charge dusts with small area-mass ratio which are near the cometary nucleus and far from the Solar. It is found that the resonant periodic orbit can be stable, and there exist stable non-resonant periodic orbits, stable resonant periodic orbits and unstable resonant periodic orbits in the potential field of cometary nuclei. The comet gravity force, solar gravity force, electric froce, magnetic force, solar radiation pressure, as well as the gas drag force are all considered to analyze the order of magnitude of these forces acting on the grains with different parameters. Let the distance of the dust grain relative to the mass centre of the cometary nucleus, the charge and the mass of the dust grain vary, respectively, fix other parameters, we clacluated the strengths of different forces.

We study galaxy shapes in the Illustris cosmological hydrodynamic simulation. We find that massive galaxies have a higher probability of being prolate. For galaxies with stellar mass larger than $10^{11}\rm M_{\odot}$, 35 out of total 839 galaxies are prolate. For 21 galaxies with stellar mass larger than $10^{12}\rm M_{\odot}$, 9 are prolate, 4 are triaxial while the others are close to being oblate. There are almost no prolate galaxies with stellar mass smaller than $3\times10^{11}\rm M_{\odot}$. We check the merger history of the prolate galaxies, and find that they are formed by major dry mergers. All the prolate galaxies have at least one such merger, with most having mass ratios between $1:1$ and $1:3$. The gas fraction (gas mass to total baryon mass) of the progenitors is 0-3 percent for nearly all these mergers, except for one whose second progenitor contains $\sim 15\%$ gas mass, while its main progenitor still contains less than $5\%$. For the 35 massive prolate galaxies that we find, 18 of them have minor axis rotation, and their angular momenta mostly come from the spin angular momenta of the progenitors (usually that of the main progenitor). We analyse the merger orbits of these prolate galaxies and find that most of them experienced a nearly radial merger orbit. Oblate galaxies with major dry mergers can have either radial or circular merger orbits. We further discuss various properties of these prolate galaxies, such as spin parameter $\lambda_{\rm R}$, spherical anisotropy parameter $\beta$, dark matter fraction, as well as inner density slopes for the stellar, dark matter and total mass distributions.

The post-Newtonian formulation of a general class of f(R) theories is set up to 3rd order approximation. It turns out that the information of a specific form of f(R) gravity is encoded in the Yukawa potential, which is contained in the perturbative expansion of the metric components. Although the Yukawa potential is canceled in the 2nd order expression of the effective refraction index of light, detailed analysis shows that the difference of the lensing effect between the f(R) gravity and general relativity does appear at the 3rd order when $\sqrt{f''(0)/f'(0)}$ is larger than the distance $d_0$ to the gravitational source. However, the difference between these two kinds of theories will disappear in the axially symmetric spacetime region. Therefore only in very rare case the f(R) theories are distinguishable from general relativity by gravitational lensing effect at the 3rd order post-Newtonian approximation.

POLAR is a new satellite-born detector aiming to measure the polarization of an unprecedented number of Gamma-Ray Bursts in the 50-500 keV energy range. The instrument, launched on-board the Tiangong-2 Chinese Space lab on the 15th of September 2016, is designed to measure the polarization of the hard X-ray flux by measuring the distribution of the azimuthal scattering angles of the incoming photons. A detailed understanding of the polarimeter and specifically of the systematic effects induced by the instrument's non-uniformity are required for this purpose. In order to study the instrument's response to polarization, POLAR underwent a beam test at the European Synchrotron Radiation Facility in France. In this paper both the beam test and the instrument performance will be described. This is followed by an overview of the Monte Carlo simulation tools developed for the instrument. Finally a comparison of the measured and simulated instrument performance will be provided and the instrument response to polarization will be presented.

The early optical afterglow emission of several gamma-ray bursts (GRBs) shows a high linear polarization degree (PD) of tens of percent, suggesting an ordered magnetic field in the emission region. The light curves are consistent with being of a reverse shock (RS) origin. However, the magnetization parameter, $\sigma$, of the outflow is unknown. If $\sigma$ is too small, an ordered field in the RS may be quickly randomized due to turbulence driven by various perturbations so that the PD may not be as high as observed. Here we use the "Athena++" relativistic MHD code to simulate a relativistic jet with an ordered magnetic field propagating into a clumpy ambient medium, with a focus on how density fluctuations may distort the ordered magnetic field and reduce PD in the RS emission for different $\sigma$ values. For a given density fluctuation, we discover a clear power-law relationship between the relative PD reduction and the $\sigma$ value of the outflow. Such a relation may be applied to estimate $\sigma$ of the GRB outflows using the polarization data of early afterglows.

In the present paper, we investigate the cosmographic problem using the bias-variance trade-off. We find that both the z-redshift and the $y=z/(1+z)$-redshift can present a small bias estimation. It means that the cosmography can describe the supernova data more accurately. Minimizing risk, it suggests that cosmography up to the second order is the best approximation. Forecasting the constraint from future measurements, we find that future supernova and redshift drift can significantly improve the constraint, thus having the potential to solve the cosmographic problem. We also exploit the values of cosmography on the deceleration parameter and equation of state of dark energy $w(z)$. We find that supernova cosmography cannot give stable estimations on them. However, much useful information was obtained, such as that the cosmography favors a complicated dark energy with varying $w(z)$, and the derivative $dw/dz<0$ for low redshift. The cosmography is helpful to model the dark energy.

In this work, expanded solutions of force-free magnetospheres on general Kerr black holes are derived through a radial distance expansion method. From the regular conditions both at the horizon and at spatial infinity, two previously known asymptotical solutions (one of them is actually an exact solution) are identified as the only solutions that satisfy the same conditions at the two boundaries. Taking them as initial conditions at the boundaries, expanded solutions up to the first few orders are derived by solving the stream equation order by order. It is shown that our extension of the exact solution can (partially) cure the problems of the solution: it leads to magnetic domination and a mostly timelike current for restricted parameters.

The ARGO-YBJ experiment, a full coverage extensive air shower (EAS) detector located at high altitude (4300 m a.s.l.) in Tibet, China, has smoothly taken data, with very high stability, since November 2007 to the beginning of 2013. The array consisted of a carpet of about 7000 m$^2$ Resistive Plate Chambers (RPCs) operated in streamer mode and equipped with both digital and analog readout, providing the measurement of particle densities up to few particles per cm$^2$. The unique detector features (full coverage, readout granularity, wide dynamic range, etc) and location (very high altitude) allowed a detailed study of the lateral density profile of charged particles at ground very close to the shower axis and its description by a proper lateral distribution function (LDF). In particular, the information collected in the first 10 m from the shower axis have been shown to provide a very effective tool for the determination of the shower development stage ("age") in the energy range 50 TeV - 10 PeV. The sensitivity of the age parameter to the mass composition of primary Cosmic Rays is also discussed.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

Understanding star formation rates (SFR) is a central goal of modern star-formation models, which mainly involve gravity, turbulence and, in some cases, magnetic fields (B-fields). However, a connection between B-fields and SFR has never been observed. Here, a comparison between the surveys of SFR and a study of cloud-field alignment - which revealed a bimodal (parallel or perpendicular) alignment - shows consistently lower SFR per solar mass for clouds almost perpendicular to the B-fields. This is evidence of B-fields being a primary regulator of SFR. The perpendicular alignment possesses a significantly higher magnetic flux than the parallel alignment and thus a stronger support of the gas against self-gravity. This results in overall lower masses of the fragmented components, which are in agreement with the lower SFR.

Based on the Effective Field Theory (EFT) of cosmological perturbations, we explicitly clarify the pathology in nonsingular cubic Galileon models and show how to cure it in EFT with new insights into this issue. With the least set of EFT operators that are capable to avoid instabilities in nonsingular cosmologies, we construct a nonsingular model dubbed the Genesis-inflation model, in which a slowly expanding phase (namely, Genesis) with increasing energy density is followed by slow-roll inflation. The spectrum of the primordial perturbation may be simulated numerically, which shows itself a large-scale cutoff, as the large-scale anomalies in CMB might be a hint for.

We investigate the observational signatures of super-Earths (i.e., Earth-to-Neptune mass planets) in their natal disks of gas and dust. Combining two-fluid global hydrodynamics simulations with a radiative transfer code, we calculate the distributions of gas and of sub-mm-sized dust in a disk perturbed by a super-Earth, synthesizing images in near-infrared scattered light and the mm-wave thermal continuum for direct comparison with observations. In low viscosity gas ($\alpha\lesssim10^{-4}$), a super-Earth opens two annular gaps to either side of its orbit by the action of Lindblad torques. This double gap and its associated gas pressure gradients cause dust particles to be dragged by gas into three rings: one ring sandwiched between the two gaps, and two rings located at the gap edges farthest from the planet. Depending on system parameters, additional rings may manifest for a single planet. A double gap located at tens of AUs from a host star in Taurus can be detected in the dust continuum by the Atacama Large Millimeter Array (ALMA) at an angular resolution of ~0".03 after two hours of integration. Ring and gap features persist in a variety of background disk profiles, last for thousands of orbits, and change their relative positions and dimensions depending on the speed and direction of planet migration. Candidate double gaps have been observed by ALMA in systems like HL Tau (D5 and D6) and TW Hya (at 37 and 43 AU); we submit that each double gap is carved by one super-Earth in nearly inviscid gas.

In this paper, we consider the singular isothermal sphere lensing model that has a spherically symmetric power-law mass distribution $\rho_{tot}(r)\sim r^{-\gamma}$. We investigate whether the mass density power-law index $\gamma$ is cosmologically evolutionary by using the strong gravitational lensing (SGL) observation, in combination with other cosmological observations. We also check whether the constraint result of $\gamma$ is affected by the cosmological model, by considering several simple dynamical dark energy models. We find that the constraint on $\gamma$ is mainly decided by the SGL observation and independent of the cosmological model, and we find no evidence for the evolution of $\gamma$ from the SGL observation.

Multi-color ($B$ $V$ $R_c$ $I_c$) CCD photometric light curves of the contact binary V502 Oph are analyzed using the Wilson-Devinney (W-D) program. The solutions reveal that V502 Oph is a W-subtype contact ($f = 35.3\,\%$) binary system. The temperature difference between its two components is $240K$ and the more massive star has a lower surface temperature. A cool spot is added in our model to account for the light curves' asymmetry (O'Connell effect) and third light is detected for the first time in the light curves' modeling. Combining the orbital inclination ($i = 76.4^{\circ}$) with the published mass function of V502 Oph, the absolute physical parameters of the two components are determined, which are $M_1= 0.46(\pm0.02)M_\odot, M_2=1.37(\pm0.02)M_\odot, R_1=0.94(\pm0.01)R_\odot, R_2=1.51(\pm0.01)R_\odot, L_1=1.13(\pm0.02)L_\odot$ and $L_2=2.49(\pm0.03)L_\odot$. The formation and evolutionary status of V502 Oph are discussed. All photoelectric and CCD times of light minimum about V502 Oph are gathered and its orbital period variations are analyzed. The results show that the orbital period of V502 Oph is decreasing continuously at a rate of $dP/dt=-1.69\times{10^{-7}}day\cdot year^{-1}$, which corresponds to a conservative mass transfer rate of $\frac{dM_{2}}{dt}=-3.01\times{10^{-8}}M_\odot/year$. The light-travel time effect (LTTE) is due to the presence of a close-in tertiary component with a period of $P_3=18.7$ years and an amplitute of 0.00402days. V502 Oph is an ideal target to test the formation and evolution theories of binary and multiple systems in which the light curves, $O - C$ curve and spectroscopic observations are comprehensively researched.

We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

The multi-TeV $\gamma$-rays from the Galactic Center (GC) have a cutoff at tens of TeV, whereas the diffuse emission has no such cutoff, which is regarded as an indication of PeV proton acceleration by the HESS experiment. It is important to understand the inconsistency and study the possibility that PeV cosmic-ray acceleration could account for the apparently contradictory point and diffuse $\gamma$-ray spectra. In this work, we propose that the cosmic rays are accelerated up to $>$PeV in GC. The interaction between cosmic rays and molecular clouds is responsible for the multi-TeV $\gamma$-ray emissions from both the point source and diffuse sources today. Enhanced by the small volume filling factor (VFF) of the clumpy structure, the absorption of the $\gamma$-rays leads to a sharp cutoff spectrum at tens of TeV produced in the GC. Away from galactic center, the VFF grows and the absorption enhancement becomes negligible. As a result, the spectra of $\gamma$-ray emissions for both point source and diffuse sources can be successfully reproduced under such self-consistent picture. In addition, a "surviving-tail" at $\sim$100 TeV is expected from the point source, which can be observed by future projects CTA and LHAASO. Neutrinos are simultaneously produced during proton-proton (PP) collision. With 5-10 years observations, the KM3Net experiment will be able to detect the PeV source according to our calculation.

The distinctive morphology of head-tail radio galaxies reveals strong interactions between the radio jets and their intra-cluster environment, the general consensus on the morphology origin of head-tail sources is that radio jets are bent by violent intra-cluster weather. We demonstrate in this paper that such strong interactions provide a great opportunity to study the jet properties and also the dynamics of intra-cluster medium (ICM). By three-dimensional magnetohydrodynamical simulations, we analyse the detailed bending process of a magnetically dominated jet, based on the magnetic tower jet model. We use stratified atmospheres modulated by wind/shock to mimic the violent intra-cluster weather. Core sloshing is found to be inevitable during the wind-cluster core interaction, which induces significant shear motion and could finally drive ICM turbulence around the jet, making it difficult for jet to survive. We perform detailed comparison between the behaviour of pure hydrodynamical jets and magnetic tower jet, and find that the jet-lobe morphology could not survive against the violent disruption in all of our pure hydrodynamical jet models. On the other hand, the head-tail morphology is well reproduced by using a magnetic tower jet model bent by wind, in which hydrodynamical instabilities are naturally suppressed and the jet could always keep its integrity under the protection of its internal magnetic fields. Finally, we also check the possibility for jet bending by shock only. We find that shock could not bend jet significantly, so could not be expected to explain the observed long tails in head-tail radio galaxies.

The key to photoelectric X-ray polarimetry is the determination of the emission direction of photoelectrons. Because of the low mass of an electron, the ionisation trajectory is not straight and the useful information needed for polarimetry is stored mostly in the initial part of the track where less energy is deposited. We present a new algorithm, based on the shortest path problem in graph theory, to reconstruct the 2D electron track from the measured image that is blurred due to transversal diffusion along drift and multiplication in the gas chamber. Compared with previous methods based on moment analysis, this algorithm allows us to identify the photoelectric interaction point more accurately and precisely for complicated tracks resulting from high energy photons or low pressure chambers. This leads to a better position resolution and a higher degree of modulation toward high energy X-rays. The new algorithm is justified using simulations and measurements with the gas pixel detector (GPD), and it should also work for other polarimetric techniques such as a time projection chamber (TPC). As the improvement is restricted in the high energy band, this new algorithm shows limited improvement for the sensitivity of GPD polarimeters, but it may have a larger potential for low-pressure TPC polarimeters.

The ARGO-YBJ detector, located at the Yangbajing Cosmic Ray Laboratory (4300 m a. s. l., Tibet, China), was a full coverage air shower array dedicated to gamma ray astronomy and cosmic ray studies. The wide field of view (~ 2 sr) and high duty cycle (> 86%), made ARGO-YBJ suitable to search for short and unexpected gamma ray emissions like gamma ray bursts (GRBs). Between 2007 November 6 and 2013 February 7, 156 satellite-triggered GRBs (24 of them with known redshift) occurred within the ARGO-YBJ field of view. A search for possible emission associated to these GRBs has been made in the two energy ranges 10-100 GeV and 10-1000 GeV. No significant excess has been found in time coincidence with the satellite detections nor in a time window of one hour after the bursts. Taking into account the EBL absorption, upper limits to the energy fluence at 99% of confidence level have been evaluated,with values ranging from ~ 10-5 erg cm-2 to ~10-1 erg cm-2.

We perform Jeans anisotropic modeling (JAM) on elliptical and spiral galaxies from the MaNGA DR13 sample. By comparing the stellar mass-to-light ratios estimated from stellar population synthesis (SPS) and from JAM, we find a similar systematic variation of the initial mass function (IMF) as in the earlier $\rm ATLAS^{3D}$ results. Early type galaxies (elliptical and lenticular) with lower velocity dispersions within one effective radius are consistent with a Chabrier-like IMF while galaxies with higher velocity dispersions are consistent with a more bottom heavy IMF such as the Salpeter IMF. Spiral galaxies have similar systematic IMF variations, but with slightly different slopes and larger scatters, due to the uncertainties caused by higher gas fractions and extinctions for these galaxies. Furthermore, we examine the effects of stellar mass-to-light ratio gradients on our JAM modeling, and find that the trends from our results becomes stronger after considering the gradients.

By using a new method of forced-field extrapolation, we study the emerging flux region AR 11850 observed by the Interface Region Imaging Spectrograph (IRIS) and Solar Dynamical Observatory (SDO). Our results suggest that the bright points (BPs) in this emerging region have responses in lines formed from the upper photosphere to the transition region, with a relatively similar morphology. They have an oscillation of several minutes according to the Atmospheric Imaging Assembly (AIA) data at 1600 and 1700 A . The ratio between the BP intensities measured in 1600 A and 1700 A filtergrams reveals that these BPs are heated differently. Our analysis of the Helioseismic and Magnetic Imager (HMI) vector magnetic field and the corresponding topology in AR11850 indicates that the BPs are located at the polarity inversion line (PIL) and most of them related with magnetic reconnection or cancelation. The heating of the BPs might be different due to different magnetic topology. We find that the heating due to the magnetic cancelation would be stronger than the case of bald patch reconnection. The plasma density rather than the magnetic field strength could play a dominant role in this process. Based on physical conditions in the lower atmosphere, our forced-field extrapolation shows consistent results between the bright arcades visible in slit-jaw image (SJI) 1400 A and the extrapolated field lines that pass through the bald patches. It provides a reliable observational evidence for testing the mechanism of magnetic reconnection for the BPs and arcades in emerging flux region, as proposed in simulation works.

In this paper, we explore the nonsingular cosmology within the framework of the Effective Field Theory(EFT) of cosmological perturbations. Due to the recently proved no-go theorem, any nonsingular cosmological models based on the cubic Galileon suffer from pathologies. We show how the EFT could help us clarify the origin of the no-go theorem, and offer us solutions to break the no-go. Particularly, we point out that the gradient instability can be removed by using some spatial derivative operators in EFT. Based on the EFT description, we obtain a realistic healthy nonsingular cosmological model, and show the perturbation spectrum can be consistent with the observations.

In the present work, we consider the cosmological constant model $\Lambda\propto\alpha^{-6}$, which is well motivated from three independent approaches. As is well known, the hint of varying fine structure constant $\alpha$ was found in 1998. If $\Lambda\propto\alpha^{-6}$ is right, it means that the cosmological constant $\Lambda$ should also be varying. Here, we try to develop a suitable framework to model this varying cosmological constant $\Lambda\propto\alpha^{-6}$, in which we view it from an interacting vacuum energy perspective. Then, we consider the observational constraints on these models by using the 293 $\Delta\alpha/\alpha$ data from the absorption systems in the spectra of distant quasars. We find that the model parameters can be tightly constrained to the very narrow ranges of ${\cal O}(10^{-5})$ typically. On the other hand, we can also view the varying cosmological constant model $\Lambda\propto\alpha^{-6}$ from another perspective, namely it can be equivalent to a model containing "dark energy" and "warm dark matter", but there is no interaction between them. We find that this is also fully consistent with the observational constraints on warm dark matter.

The complex interplay between turbulence, magnetic fields, and self-gravity leads to the formation of molecular clouds out of the diffuse interstellar medium (ISM). One avenue of studying this interplay is by analyzing statistical features derived from observations, where the interpretation of these features is greatly facilitated by comparisons with numerical simulations. Here we focus on the statistical anisotropy present in synthetic maps of velocity centroid data, which we derive from three-dimensional magnetohydrodynamic simulations of a turbulent, magnetized, self-gravitating patch of ISM. We study how the orientation and magnitude of the velocity anisotropy correlate with the magnetic field and with the structures generated by gravitational collapse. Motivated by recent observational constraints, our simulations focus on the supersonic (sonic Mach number $\mathcal{M} \approx 2 - 17$) but sub- to trans-alfvenic (alfvenic Mach number $\mathcal{M}_A \approx 0.2 - 1.2$) turbulence regime, and we consider clouds which are barely to mildly magnetically supercritical (mass-to-flux ratio equal to once or twice the critical value). Additionally we explore the impact of the turbulence driving mechanism (solenoidal or compressive) on the velocity anisotropy. While we confirm previous findings that the velocity anisotropy generally aligns well with the plane-of-sky magnetic field, our inclusion of the effects of self-gravity reveals that in regions of higher column density, the velocity anisotropy may be destroyed or even reoriented to align with the gravitationally formed structures. We provide evidence that this effect is not necessarily due to the increase of $\mathcal{M}_A$ inside the high-density regions.

Most of relatively warm, unevolved metal-poor stars ($T_{\rm eff}\gtrsim 5800\,\mathrm{K}$ and $[{\rm Fe/H}]\lesssim -1.5$) exhibit almost constant lithium abundances, irrespective of metallicity or effective temperature, and thus form the so-called Spite plateau. This was originally interpreted as arising from lithium created by the Big Bang nucleosynthesis. Recent observations, however, have revealed that ultra metal-poor stars (UMP stars; $[\mathrm{Fe/H}]<-4.0$) have significantly lower lithium abundances than that of the plateau. Since most of the UMP stars are CEMP-no stars, carbon-enhanced metal-poor stars with no excess of neutron-capture elements, a connection between the carbon enhancement and lithium depletion is suspected. A straightforward approach to this question is to investigate carbon-normal UMP stars. However only one object is known in this class. As an alternative, we have determined lithium abundances for two CEMP-no main-sequence turn-off stars with metallicities $[{\rm Fe/H}]\sim -3.0$, where there are numerous carbon-normal stars with available lithium abundances that can be considered. Our 1D LTE analysis indicates that the two CEMP-no stars have lithium abundances that are consistent with values near the plateau, which suggests that carbon enhancement and lithium depletion are not directly related. Instead, our results suggest that extremely low iron abundance is a fundamental cause to depleted lithium in UMP stars.

A number of transition disks exhibit significant azimuthal asymmetries in thermal dust emission. One possible origin for these asymmetries is dust trapping in vortices formed at the edges of dead zones. We carry out high-resolution, two-dimensional hydrodynamic simulations of this scenario, including the effects of dust feedback. We find that, although feedback weakens the vortices and slows down the process of dust accumulation, the dust distribution in the disk can nonetheless remain asymmetric for many thousands of orbits. We show that even after $10^4$ orbits, or $2.5$ Myr when scaled to the parameters of Oph IRS 48 (a significant fraction of its age), the dust is not dispersed into an axisymmetric ring, in contrast to the case of a vortex formed by a planet. This is because accumulation of mass at the dead zone edge constantly replenishes the vortex, preventing it from being fully destroyed. We produce synthetic dust emission images using our simulation results. We find that multiple small clumps of dust may be distributed azimuthally. These clumps, if not resolved from one another, appear as a single large feature. A defining characteristic of a disk with a dead zone edge is that an asymmetric feature is accompanied by a ring of dust located about twice as far from the central star.

High-resolution optical spectroscopy was conducted for the metal-poor post-AGB star CC Lyr to determine its chemical abundances and spectral line profiles. Our standard abundance analysis confirms its extremely low metallicity ([Fe/H]<-3.5) and a clear correlation between abundance ratios and the condensation temperature for 11 elements, indicating that dust depletion is the cause of the abundance anomaly of this object. The very low abundances of Sr and Ba, which are detected for the first time for this object, suggest that heavy neutron-capture elements are not significantly enhanced in this object by the s-process during its evolution through AGB phase. Radial velocity of this object and profiles of some atomic absorption lines show variations depending on pulsation phases, which could be formed by dynamics of the atmosphere rather than by binarity or contributions of circumstellar absorption. On the other hand, the H-alpha emission with double peaks shows no evident velocity shift, suggesting that the emission is originating from the circumstellar matter, presumably the rotating disk around the object.

Kink instabilities are likely to occur in the current-carrying magnetized plasma jets. Recent observations of the blazar radiation and polarization signatures suggest that the blazar emission region may be considerably magnetized. While the kink instability has been studied with first-principle magnetohydrodynamic (MHD) simulations, the corresponding time-dependent radiation and polarization signatures have not been investigated. In this paper, we perform comprehensive polarization-dependent radiation modeling of the kink instability in the blazar emission region based on relativistic MHD (RMHD) simulations. We find that the kink instability may give rise to strong flares with polarization angle (PA) swings or weak flares with polarization fluctuations, depending on the initial magnetic topology and magnetization. These findings are consistent with observations. Compared with the shock model, the kink model generates polarization signatures that are in better agreement with the general polarization observations. Therefore, we suggest that kink instabilities may widely exist in the jet environment, and provide an efficient way to convert the magnetic energy and produce multiwavelength flares and polarization variations.

We present a study on the stellar age and metallicity distributions for 1105 galaxies using the STARLIGHT software on MaNGA integral field spectra. We derive age and metallicity gradients by fitting straight lines to the radial profiles, and explore their correlations with total stellar mass M*, NUV-r colour and environments, as identified by both the large scale structure (LSS) type and the local density. We find that the mean age and metallicity gradients are close to zero but slightly negative, which is consistent with the inside-out formation scenario. Within our sample, we find that both the age and metallicity gradients show weak or no correlation with either the LSS type or local density environment. In addition, we also study the environmental dependence of age and metallicity values at the effective radii. The age and metallicity values are highly correlated with M* and NUV-r and are also dependent on LSS type as well as local density. Low-mass galaxies tend to be younger and have lower metallicity in low-density environments while high-mass galaxies are less affected by environment.

Simulations suggest collisionless steady-state magnetic reconnection of Harris-type current sheets proceeds with a rate of order 0.1, independent of dissipation mechanism. We argue this long-standing puzzle is a result of constraints at the magnetohydrodynamic (MHD) scale. We perform a scaling analysis of the reconnection rate as a function of the opening angle made by the upstream magnetic fields, finding a maximum reconnection rate close to 0.2. The predictions compare favorably to particle-in-cell simulations of relativistic electron-positron and non-relativistic electron-proton reconnection. The fact that simulated reconnection rates are close to the predicted maximum suggests reconnection proceeds near the most efficient state allowed at the MHD-scale. The rate near the maximum is relatively insensitive to the opening angle, potentially explaining why reconnection has a similar fast rate in differing models.

We present a new implementation of star formation in cosmological simulations, by considering star clusters as a unit of star formation. Cluster particles grow in mass over several million years at the rate determined by local gas properties, with high time resolution. The particle growth is terminated by its own energy and momentum feedback on the interstellar medium. We test this implementation for Milky Way-sized galaxies at high redshift, by comparing the properties of model clusters with observations of young star clusters. We find that the cluster initial mass function is best described by a Schechter function rather than a single power law. In agreement with observations, at low masses the logarithmic slope is $\alpha\approx 1.8-2$, while the cutoff at high mass scales with the star formation rate. A related trend is a positive correlation between the surface density of star formation rate and fraction of stars contained in massive clusters. Both trends indicate that the formation of massive star clusters is preferred during bursts of star formation. These bursts are often associated with major merger events. We also find that the median timescale for cluster formation ranges from 0.5 to 4 Myr and decreases systematically with increasing star formation efficiency. Local variations in the gas density and cluster accretion rate naturally lead to the scatter of the overall formation efficiency by an order of magnitude, even when the instantaneous efficiency is kept constant. Comparison of the formation timescale with the observed age spread of young star clusters provides an additional important constraint on the modeling of star formation and feedback schemes.

The coupled dark energy model provides a possible approach to mitigate the coincidence problem of cosmological standard model. Here, the coupling term is assumed as $\bar{Q}=3H\xi_x\bar{\rho}_x$, which is related to the interaction rate and energy density of dark energy. We derive the background and perturbation evolution equations for several coupled models. Then, we test these models by currently available cosmic observations which include cosmic microwave background radiation from Planck 2015, baryon acoustic oscillation, type Ia supernovae, $f\sigma_8(z)$ data points from redshift-space distortions, and weak gravitational lensing. The constraint results tell us there is no evidence of interaction at 2$\sigma$ level, it is very hard to distinguish different coupled models from other ones.

We present spectroscopic measurements on the night sky of Xinglong Observatory for a period of 12 years from 2004 to 2015. The spectra were obtained on moonless clear nights using the OMR spectrograph mounted on a 2.16-m reflector with a wavelength coverage of 4000-7000A. The night sky spectrum shows the presence of emission lines from Hg I and Na I due to local artificial sources, along with the atmospheric emission lines, i.e., O I and OH molecules, indicating the existence of light pollution. We have monitored the night sky brightness during the whole night and found some decrement in the sky brightness with time, but the change is not significant. Also, we monitored the light pollution level in different azimuthal directions and found that the influence of light pollution from the direction of Beijing is stronger compared with that from the direction of Tangshan and other areas. An analysis of night sky spectra for the entire data set suggested that the zenith sky brightness of Xinglong Observatory has brightened by about 0.5 mag arcsec -2 in the V and B bands from 2004 to 2015. We recommend consecutive spectroscopic measurements of the night sky brightness at Xinglong Observatory in the future, not only for monitoring but also for scientific reference.

Historically, the weak s-process contribution to metal-poor stars is thought to be extremely small, due to the effect of the secondary-like nature of the neutron source 22Ne(a;n)25Mg in massive stars, which means that metal-poor weak s-process stars could not be found. ET0097 is the first observed carbon-enhanced metal-poor (CEMP) star in the Sculptor dwarf spheroidal galaxy. Because C is enriched and the elements heavier than Ba are not overabundant, ET0097 can be classified as a CEMP-no star. However, this star shows overabundances of lighter n-capture elements (i.e., Sr, Y and Zr). In this work, having adopted the abundance decomposition approach, we investigate the astrophysical origins of the elements in ET0097. We find that the light elements and iron-peak elements (from O to Zn) of the star mainly originate from the primary process of massive stars and the heavier n-capture elements (heavier than Ba) mainly come from the main r-process. However, the lighter n-capture elements such as Sr, Y and Zr should mainly come from the primary weak s-process. The contributed fractions of the primary weak s-process to the Sr, Y and Zr abundances of ET0097 are about 82%, 84% and 58% respectively, suggesting that the CEMP star ET0097 should have the footprints of the weak s-process. The derived result should be a significant evidence that the weak s-process elements can be produced in metal-poor massive stars.

Additional physics beyond standard hydrodynamics is needed to fully model the intracluster medium (ICM); however, as we move to more sophisticated models, it is important to consider the role of magnetic fields and the way the fluid approximation breaks down. This paper represents a first step towards developing a self-consistent model of the ICM by characterizing the statistical properties of magnetic fields in cosmological simulations of galaxy clusters. We find that plasma conditions are largely homogeneous across a range of cluster masses and relaxation states. We also find that the magnetic field length scales are resolution dependent and not based on any particular physical process. Energy transfer mechanisms and scales are also identified, and imply the existence of small scale dynamo action. The scales of the small scale dynamo are resolution limited and driven by numerical resistivity and viscosity.

The MaNGA Survey (Mapping Nearby Galaxies at Apache Point Observatory) is one of three core programs in the Sloan Digital Sky Survey IV. It is obtaining integral field spectroscopy (IFS) for 10K nearby galaxies at a spectral resolution of R~2000 from 3,622-10,354A. The design of the survey is driven by a set of science requirements on the precision of estimates of the following properties: star formation rate surface density, gas metallicity, stellar population age, metallicity, and abundance ratio, and their gradients; stellar and gas kinematics; and enclosed gravitational mass as a function of radius. We describe how these science requirements set the depth of the observations and dictate sample selection. The majority of targeted galaxies are selected to ensure uniform spatial coverage in units of effective radius (Re) while maximizing spatial resolution. About 2/3 of the sample is covered out to 1.5Re (Primary sample), and 1/3 of the sample is covered to 2.5Re (Secondary sample). We describe the survey execution with details that would be useful in the design of similar future surveys. We also present statistics on the achieved data quality, specifically, the point spread function, sampling uniformity, spectral resolution, sky subtraction, and flux calibration. For our Primary sample, the median r-band signal-to-noise ratio is ~73 per 1.4A pixel for spectra stacked between 1-1.5 Re. Measurements of various galaxy properties from the first year data show that we are meeting or exceeding the defined requirements for the majority of our science goals.

In bounce inflation scenario, the inflation is singularity-free, while the advantages of inflation are reserved. We analytically calculate the power spectrum of its primordial gravitational waves (GWs), and show a universal result including the physics of the bounce phase. The spectrum acquires a cutoff at large scale, while the oscillation around the cutoff scale is quite drastic, which is determined by the details of bounce. Our work highlights that the primordial GWs at large scale may encode the physics of the bounce ever happened at about $\sim 60$ efolds before inflation.

The Xinglong 2.16-m reflector is the first 2-meter class astronomical telescope in China. It was jointly designed and built by the Nanjing Astronomical Instruments Factory (NAIF), Beijing Astronomical Observatory (now National Astronomical Observatories, Chinese Academy of Sciences, NAOC) and Institute of Automation, Chinese Academy of Sciences in 1989. It is Ritchey-Chr\'{e}tien (R-C) reflector on an English equatorial mount and the effective aperture is 2.16 meters. It had been the largest optical telescope in China for $\sim18$ years until the Guoshoujing Telescope (also called Large Sky Area Multi-Object Fiber Spectroscopic Telescope, LAMOST) and the Lijiang 2.4-m telescope were built. At present, there are three main instruments on the Cassegrain focus available: the Beijing Faint Object Spectrograph and Camera (BFOSC) for direct imaging and low resolution ($R\sim500-2000$) spectroscopy, the spectrograph made by Optomechanics Research Inc. (OMR) for low resolution spectroscopy (the spectral resolutions are similar to those of BFOSC) and the fiber-fed High Resolution Spectrograph (HRS, $R\sim30000-65000$). The telescope is widely open to astronomers all over China as well as international astronomical observers. Each year there are more than 40 ongoing observing projects, including 6-8 key projects. Recently, some new techniques and instruments (e.g., astro-frequency comb calibration system, polarimeter and adaptive optics) have been or will be tested on the telescope to extend its observing abilities.

A model of global magnetic reconnection rate in relativistic collisionless plasmas is developed and validated by the fully kinetic simulation. Through considering the force balance at the upstream and downstream of the diffusion region, we show that the global rate is bounded by a value $\sim 0.3$ even when the local rate goes up to $\sim O(1)$ and the local inflow speed approaches the speed of light in strongly magnetized plasmas. The derived model is general and can be applied to magnetic reconnection under widely different circumstances.

We report on a measurement of thermal neutrons, generated by the hadronic component of extensive air showers (EAS), by means of a small array of EN-detectors developed for the PRISMA project (PRImary Spectrum Measurement Array), novel devices based on a compound alloy of ZnS(Ag) and $^{6}$LiF. This array has been operated within the ARGO-YBJ experiment at the high altitude Cosmic Ray Observatory in Yangbajing (Tibet, 4300 m a.s.l.). Due to the tight correlation between the air shower hadrons and thermal neutrons, this technique can be envisaged as a simple way to estimate the number of high energy hadrons in EAS. Coincident events generated by primary cosmic rays of energies greater than 100 TeV have been selected and analyzed. The EN-detectors have been used to record simultaneously thermal neutrons and the air shower electromagnetic component. The density distributions of both components and the total number of thermal neutrons have been measured. The correlation of these data with the measurements carried out by ARGO-YBJ confirms the excellent performance of the EN-detector.

Plasma energization through magnetic reconnection in the magnetically-dominated regime featured by low plasma beta ($\beta = 8 \pi nkT_0/B^2 \ll 1$) and/or high magnetization ($\sigma = B^2/(4 \pi nmc^2) \gg 1$) is important in a series of astrophysical systems such as solar flares, pulsar wind nebula, and relativistic jets from black holes, etc. In this paper, we review the recent progress on kinetic simulations of this process and further discuss plasma dynamics and particle acceleration in a low-$\beta$ reconnection layer that consists of electron-positron pairs. We also examine the effect of different initial thermal temperatures on the resulting particle energy spectra. While earlier papers have concluded that the spectral index is smaller for higher $\sigma$, our simulations show that the spectral index approaches $p=1$ for sufficiently low plasma $\beta$, even if $\sigma \sim 1$. Since this predicted spectral index in the idealized limit is harder than most observations, it is important to consider effects that can lead to a softer spectrum such as open boundary simulations. We also remark that the effects of 3D reconnection physics and turbulence on reconnection need to be addressed in the future.

We studied the magnetic topology of active region 12158 on 2014 September 10 and compared it with the observations before and early in the flare which begins at 17:21 UT (SOL2014-09-10T17:45:00). Our results show that the sigmoidal structure and flare ribbons of this active region observed by SDO/AIA can be well reproduced from a Grad-Rubin non linear force free field extrapolation method. Various inverse-S and -J shaped magnetic field lines, that surround a coronal flux rope, coincide with the sigmoid as observed in different extreme ultraviolet wavelengths, including its multi-threaded curved ends. Also, the observed distribution of surface currents in the magnetic polarity where it was not prescribed is well reproduced. This validates our numerical implementation and set-up of the Grad-Rubin method. The modeled double inverse-J shaped Quasi-Separatrix Layer (QSL) footprints match the observed flare ribbons during the rising phase of the flare, including their hooked parts. The spiral-like shape of the latter may be related to a complex pre-eruptive flux rope with more than one turn of twist, as obtained in the model. These ribbon-associated flux-rope QSL-footprints are consistent with the new standard flare model in 3D, with the presence of a hyperbolic flux tube located below an inverse tear drop shaped coronal QSL. This is a new step forward forecasting the locations of reconnection and ribbons in solar flares, and the geometrical properties of eruptive flux ropes.

The jet composition and energy dissipation mechanism of Gamma-ray bursts (GRBs) and Blazars are fundamental questions which remain not fully understood. One plausible model is to interpret the $\gamma$-ray emission of GRBs and optical emission of blazars as synchrotron radiation of electrons accelerated from the collision-induced magnetic dissipation regions in Poynting-flux-dominated jets. The polarization observation is an important and independent information to test this model. Based on our recent 3D relativistic MHD simulations of collision-induced magnetic dissipation of magnetically dominated blobs, here we perform calculations of the polarization properties of the emission in the dissipation region and apply the results to model the polarization observational data of GRB prompt emission and blazar optical emission. We show that the same numerical model with different input parameters can reproduce well the observational data of both GRBs and blazars, especially the $90^{\circ}$ polarization angle (PA) change in GRB 100826A and the $180^{\circ}$ PA swing in Blazar 3C279. This supports a unified model for GRB and blazar jets, suggesting that collision-induced magnetic reconnection is a common physical mechanism to power the relativistic jet emission from events with very different black hole masses.

We have assembled the historical variability data of Mrk 421 at radio 15 GHz, X-ray and gamma-ray bands, spanning about 6.3, 10.3 and 7.5 yr, respective- ly. We analyzed the variability by using three methods. The results indicated that there is a period of 287.6 +- 4.4 days for 15 GHz, 309.5 +- 5.8 days for X- ray and 283.4+- 4.7 days for gamma-ray, respectively. This period can be reasonably explained by the nonballistic helical motion of the emitting material. The cor- relation analysis suggested that the variabilities of radio 15 GHz, X-ray and 0-ray are remarkable correlated, and the emission of radio 15 GHz lags behind that of X-ray, and the X-ray flux lags behind the gamma-ray. This suggests that the gamma-ray derives from inverse Compton (IC) scattering of the synchrotron photons, supporting the synchrotron self-Compton (SSC) model. Moreover, the time delay between different wavebands could be explained by the shock-in-jet models, in which a moving emission region produces the radio to gamma-ray activity, implying that the emission region of gamma-ray is closer to the center than ones of X-ray and radio emission.

Magnetic reconnection is a leading mechanism for dissipating magnetic energy and accelerating nonthermal particles in Poynting-flux dominated flows. In this letter, we investigate nonthermal particle acceleration during magnetic reconnection in a magnetically-dominated ion-electron plasma using fully kinetic simulations. For an ion-electron plasma with the total magnetization $\sigma_0=B^2/(4\pi n(m_i+m_e)c^2)$, the magnetization for each species is $\sigma_i \sim \sigma_0$ and $\sigma_e \sim (m_i/m_e) \sigma_0$, respectively. We have studied the magnetically dominated regime by varying $\sigma_{e} = 10^3 - 10^5$ with initial ion and electron temperatures $T_i = T_e = 5 - 20 m_ec^2$ and mass ratio $m_i/m_e = 1 - 1836$. The results demonstrate that reconnection quickly establishes power-law energy distributions for both electrons and ions within several ($2-3$) light-crossing times. For the cases with periodic boundary conditions, the power-law index is $1<s<2$ for both electrons and ions. The hard spectra limit the power-law energies for electrons and ions to be $\gamma_{be} \sim \sigma_e$ and $\gamma_{bi} \sim \sigma_i$, respectively. The main acceleration mechanism is a Fermi-like acceleration through the drift motions of charged particles. When comparing the spectra for electrons and ions in momentum space, the spectral indices $s_p$ are identical as predicted in Fermi acceleration. We also find that the bulk flow can carry a significant amount of energy during the simulations. We discuss the implication of this study in the context of Poynting-flux dominated jets and pulsar winds especially the applications for explaining the nonthermal high-energy emissions.

We use extensive global two-dimensional hydrodynamic disk gas+dust simulations with embedded planets, coupled with three dimensional radiative transfer calculations, to model the dust ring and gap structures in the HL Tau protoplanetary disk observed with the Atacama Large Millimeter/Submillimeter Array (ALMA). We include the self-gravity of disk gas and dust components and make reasonable choices of disk parameters, assuming an already settled dust distribution and no planet migration. We can obtain quite adequate fits to the observed dust emission using three planets with masses 0.35, 0.17, and 0.26 $M_{Jup}$ at 13.1, 33.0, and 68.6 AU, respectively. Implications for the planet formation as well as the limitations of this scenario are discussed.

Nowadays, $f(R)$ theory has been one of the leading modified gravity theories to explain the current accelerated expansion of the universe, without invoking dark energy. It is of interest to find the exact cosmological solutions of $f(R)$ theories. Besides other methods, symmetry has been proved as a powerful tool to find exact solutions. On the other hand, symmetry might hint the deep physical structure of a theory, and hence considering symmetry is also well motivated. As is well known, Noether symmetry has been extensively used in physics. Recently, the so-called Hojman symmetry was also considered in the literature. Hojman symmetry directly deals with the equations of motion, rather than Lagrangian or Hamiltonian, unlike Noether symmetry. In this work, we consider Hojman symmetry in $f(R)$ theories in both the metric and Palatini formalisms, and find the corresponding exact cosmological solutions of $f(R)$ theories via Hojman symmetry. There exist some new solutions significantly different from the ones obtained by using Noether symmetry in $f(R)$ theories. To our knowledge, they also have not been found previously in the literature. This work confirms that Hojman symmetry can bring new features to cosmology and gravity theories.

The study of reionization history plays an important role in understanding the evolution of our universe. It is commonly believed that the intergalactic medium (IGM) in our universe are fully ionized today, however the reionizing process remains to be mysterious. A simple instantaneous reionization process is usually adopted in modern cosmology without direct observational evidence. However, the history of ionization fraction, $x_e(z)$ will influence cosmic microwave background (CMB) observables and constraints on optical depth $\tau$. With the mocked future data sets based on featured reionization model, we find the bias on $\tau$ introduced by instantaneous model can not be neglected. In this paper, we study the cosmic reionization history in a model independent way, the so called principle component analysis (PCA) method, and reconstruct $x_e (z)$ at different redshift $z$ with the data sets of Planck, WMAP 9 years temperature and polarization power spectra, combining with the baryon acoustic oscillation (BAO) from galaxy survey and type Ia supernovae (SN) Union 2.1 sample respectively. The results show that reconstructed $x_e(z)$ is consistent with instantaneous behavior, however, there exists slight deviation from this behavior at some epoch. With PCA method, after abandoning the noisy modes, we get stronger constraints, and the hints for featured $x_e(z)$ evolution could become a little more obvious.

It is of great interest to connect cosmology in the early universe to the Standard Model of particle physics. In this paper, we try to construct a bounce inflation model with the standard model Higgs boson, where the one loop correction is taken into account in the effective potential of Higgs field. In this model, a Galileon term has been introduced to eliminate the ghost mode when bounce happens. Moreover, due to the fact that the Fermion loop correction can make part of the Higgs potential negative, one naturally obtains a large equation of state(EoS) parameter in the contracting phase, which can eliminate the anisotropy problem. After the bounce, the model can drive the universe into the standard higgs inflation phase, which can generate nearly scale-invariant power spectrum.

We demonstrate that particles are regularly accelerated while experiencing curvature drift in flows driven by magnetic tension. Some examples of such flows include spontaneous turbulent reconnection and decaying magnetohydrodynamic (MHD) turbulence, where magnetic field relaxes to a lower-energy configuration and transfers part of its energy to kinetic motions. The opposite process, such as dynamo, will actually result in the net cooling of particles by the curvature drift. Being very generic, this acceleration mechanism is likely to be responsible in production of non-thermal particle distribution in magnetically-dominant environments such as solar chromosphere, pulsar magnetosphere, jets from supermassive black holes, $\gamma$-ray bursts, etc.

The optical radiation and polarization signatures in blazars are known to be highly variable during flaring activities. It is frequently argued that shocks are the main driver of the flaring events. However, the spectral variability modelings generally lack detailed considerations of the self-consistent magnetic field evolution modeling, thus so far the associated optical polarization signatures are poorly understood. We present the first simultaneous modeling of the optical radiation and polarization signatures based on 3D magnetohydrodynamic simulations of relativistic shocks in the blazar emission environment, with the simplest physical assumptions. By comparing the results with observations, we find that shocks in a weakly magnetized environment will largely lead to significant changes in the optical polarization signatures, which are seldom seen in observations. Hence an emission region with relatively strong magnetization is preferred. In such an environment, slow shocks may produce minor flares with either erratic polarization fluctuations or considerable polarization variations, depending on the parameters; fast shocks can produce major flares with smooth PA rotations. In addition, the magnetic fields in both cases are observed to actively revert to the original topology after the shocks. All these features are consistent with observations. Future observations of the radiation and polarization signatures will further constrain the flaring mechanism and the blazar emission environment.

The Ngari (Ali) prefecture of Tibet, one of the highest areas in the world, has recently emerged as a promising site for future astronomical observation. Here we use 31 years of reanalysis data from the Climate Forecast System Reanalysis (CFSR) to examine the astroclimatology of Ngari, using the recently-erected Ali Observatory at Shiquanhe (5~047~m above mean sea level) as the representative site. We find the percentage of photometric night, median atmospheric seeing and median precipitable water vapor (PWV) of the Shiquanhe site to be $57\%$, $0.8"$ and 2.5~mm, comparable some of the world's best astronomical observatories. Additional calculation supports the Shiquanhe region as one of the better sites for astronomical observations over the Tibetan Plateau. Based on the studies taken at comparable environment at Atacama, extraordinary observing condition may be possible at the few vehicle-accessible 6~000~m heights in the Shiquanhe region. Such possibility should be thoroughly investigated in future.

We are interested in the periodic motion and bifurcations near the surface of an asteroid. The gravity field of an irregular asteroid and the equation of motion of a particle near the surface of an asteroid are studied. The periodic motions around the major body of triple asteroid 216 Kleopatra and the OSIRIS REx mission target asteroid 101955 Bennu are discussed. We find that motion near the surface of an irregular asteroid is quite different from the motion near the surface of a homoplastically spheroidal celestial body. The periodic motions around the asteroid 101955 Bennu and 216 Kleopatra indicate that the geometrical shapes of the orbits are probably very sophisticated. There exist both stable periodic motions and unstable periodic motions near the surface of the same irregular asteroid. This periodic motion which is unstable can be resonant or non resonant. The period doubling bifurcation and pseudo period doubling bifurcation of periodic orbits coexist in the same gravity field of the primary of the triple asteroid 216 Kleopatra. It is found that both of the period doubling bifurcations of periodic orbits and pseudo period-doubling bifurcation of periodic orbits have four different paths. The pseudo period doubling bifurcation found in the potential field of primary of triple asteroid 216 Kleopatra shows that there exist stable periodic orbits near the primary s equatorial plane, which gives an explanation for the motion stability of the triple asteroid 216 Kleopatra s two moonlets, Alexhelios and Cleoselene.

We highlight distinct and systematic observational features of magnetic field morphologies in polarized submm dust continuum. We illustrate this with specific examples and show statistical trends from a sample of 50 star-forming regions.