Theoretically, CME kinematics are related to magnetic reconnection processes in the solar corona. However, the current quantitative understanding of this relationship is based on the analysis of only a handful of events. Here we report a statistical study of 60 CME-flare events from August 2010 to December 2013. We investigate kinematic properties of CMEs and magnetic reconnection in the low corona during the early phase of the eruptions, by combining limb observations from STEREO with simultaneous on-disk views from SDO. For a subset of 42 events with reconnection rate evaluated by the magnetic fluxes swept by the flare ribbons on the solar disk observed from SDO, we find a strong correlation between the peak CME acceleration and the peak reconnection rate. Also, the maximum velocities of relatively fast CMEs (> 600 km/s) are positively correlated with the reconnection flux, but no such correlation is found for slow CMEs. A time-lagged correlation analysis suggests that the distribution of the time lag of CME acceleration relative to reconnection rate exhibits three peaks, approximately 10 minutes apart, and on average, acceleration-lead events have smaller reconnection rates. We further compare the CME total mechanical energy with the estimated energy in the current sheet. The comparison suggests that, for small-flare events, reconnection in the current sheet alone is insufficient to fuel CMEs. Results from this study suggest that flare reconnection may dominate the acceleration of fast CMEs, but for events of slow CMEs and weak reconnection, other mechanisms may be more important.

Crater ellipticity determination is a complex and time consuming task that so far has evaded successful automation. We train a state of the art computer vision algorithm to identify craters in Lunar digital elevation maps and retrieve their sizes and 2D shapes. The computational backbone of the model is MaskRCNN, an "instance segmentation" general framework that detects craters in an image while simultaneously producing a mask for each crater that traces its outer rim. Our post-processing pipeline then finds the closest fitting ellipse to these masks, allowing us to retrieve the crater ellipticities. Our model is able to correctly identify 87\% of known craters in the longitude range we hid from the network during training and validation (test set), while predicting thousands of additional craters not present in our training data. Manual validation of a subset of these "new" craters indicates that a majority of them are real, which we take as an indicator of the strength of our model in learning to identify craters, despite incomplete training data. The crater size, ellipticity, and depth distributions predicted by our model are consistent with human-generated results. The model allows us to perform a large scale search for differences in crater diameter and shape distributions between the lunar highlands and maria, and we exclude any such differences with a high statistical significance. The predicted test set catalogue and trained model are available here: https://github.com/malidib/Craters_MaskRCNN/.

ZTF J153932.16+502738.8 (ZTFJ1539) is an eclipsing double-white-dwarf system with an orbital period of 6.91 minutes, and is a significant source of LISA detecting gravitational wave. However, the massive white dwarf (WD) with mass of about 0.61 M$_\odot$ has a high effective temperature (48900 K), and the lower mass WD with mass of about 0.21 M$_{\odot}$ has a low effective temperature($<$10000 K). It is challenging the popular theory of binary evolution. We investigate the formation of ZTFJ1539 via nova and Algol scenarios. Assuming that the massive WD in ZTFJ1539 just experiences a thermalnuclear runaway, nova scenario can explain the effective temperatures of two WDs in ZTFJ1539. However, in order to enlarging a semi-detached orbit of about 4---5 minutes to a detached orbit of about 7 minutes, nova scenario needs a much high kick velocity of about 200 km s$^{-1}$ during nova eruption. The high kick velocity can result in high eccentricity of about 0.2---0.6. Algol scenario can also produce ZTFJ1539 if we take a high efficient parameter for ejecting common envelope and enhance the mass-loss rate via stellar wind trigger by tidal effect.

The formation of bimodal dust species (namely the silicate and amorphous carbon dust grains coexistent) in a nova eruption is an open problem. According to the nova model simulated by Modules for Experiments in Stellar Astrophysics code, we calculate the formation and growth of carbon (C) and forsterite (Mg2SiO4) dust grains purely in nova ejecta for the fee-expansion model and the radiative shock model by assuming spherical geometry of the nova ejecta. In these models, the chemical properties of pre-existing circumstellar medium are not taken into account. In the free-expansion model, the nova ejecta is not an idea environment for dust nucleation. However, it can efficiently produce dust in the radiative shock model. We estimate that every nova can produce C grains with an average mass of about $10^{-9}$ and $10^{-8}$ ${\rm M_\odot}$, and Mg$_2$SiO$_4$ grains with an average mass of about $10^{-8}$ and $10^{-7}$ ${\rm M_\odot}$. Based on the mass of ejected gas, the ratio of dust to gas is about 1\%. The C grains form first after several or tens of days of nova eruption. After that, the Mg$_2$SiO$_4$ grains begin to grow in tens of days, which is consistent with observations.

We report measurements of the charged daughter fraction of $^{218}$Po as a result of the $^{222}$Rn alpha decay, and the mobility of $^{218}$Po$^+$ ions, using radon-polonium coincidences from the $^{238}$U chain identified in 532 live-days of DarkSide-50 WIMP-search data. The fraction of $^{218}$Po that is charged is found to be 0.37$\pm$0.03 and the mobility of $^{218}$Po$^+$ is (8.6$\pm$0.1)$\times$10$^{-4}$$\frac{\text{cm}^2}{\text{Vs}}$.

Subdwarf A-type stars (sdAs) are objects that have hydrogen-rich spectra with surface gravity similar to that of hot subdwarf stars but effective temperature below the zero-age horizontal branch (ZAHB). They are considered to be metal-poor main sequence (MS) stars or extremely low-mass white dwarfs (ELM WDs). In this work, using the stellar evolution code, Modules for Experiments in Stellar Astrophysics (MESA), we investigate the sdAs formed both by the evolution of (pre-)ELM WDs in double degenerate systems (DDs) and metal-poor main sequence with the single evolution models. We find that both of the evolutionary tracks of ELM WDs and metal-poor MSs can explain the observation properties of sdAs. However, the proportion between these two populations are uncertain. In this work, we adopt the method of binary population synthesis of both ELM WDs in the disk and metal-poor MS in the halo to obtain the populations of ELM WDs and metal-poor MSs at different stellar population ages and calculate their proportions. We find that the proportion of metal-poor MSs to sdA population for a stellar population of 10 Gyr is $\sim$98.5 percent, which is consistent with the conclusion that most sdAs ($>$95 percent) are metal-poor MSs. And the proportion of ELM WDs (metal-poor MSs) to sdA population increases (decreases) from 0.1\% (99.9\%) to 20\% (80\%) with the stellar population ages from 5 to 13.7 Gyr.

We selected 90 massive star-forming clumps with strong N2H+, HCO+, HCN, and HNC emission from the Millimetre Astronomy Legacy Team 90 GHz survey. We obtained Herschel data for all 90 sources and NRAO VLA Sky Survey data for 51 of them. We convolved and regridded all images to the same resolution and pixel size and derived the temperature, H2 column density, molecules' abundances and abundance, and ratios of each pixel. Our analysis yields three main conclusions. First, the abundances of N2H+, HCO+, HCN, and HNC increase when the column density decreases and the temperature increases, with spatial variations in their abundances dominated by changes in the H2 column density. Second, the abundance ratios between N2H+, HCO+, HCN, and HNC also display systemic variations as a function of the column density due to the chemical properties of these molecules. Third, the sources associated with the 20 cm continuum emission can be classified into four types based on the behavior of the abundances of the four molecules considered here as a function of this emission. The variations of the first three types could also be attributed to the variation of the H2 column density.

Isolated millisecond pulsars (IMSPs) are a topic of academic contention. There are various models to explain their formation. We explore the formation of IMSP via quark novae (QN). During this formation process, low-mass X-ray binaries (LMXBs) are disrupted when the mass of the neutron star (NS) reaches 1.8 $\rm M_\odot$. Using population synthesis, this work estimates that the Galactic birthrate of QN-produced IMSPs lies between $\sim 9.5\times10^{-6}$ and $\sim1.7\times10^{-4}$ $\rm yr^{-1}$. The uncertainties shown in our experiment model is due to the QN's kick velocity. Furthermore, our findings not only show that QN-produced IMSPs are statistically more significant than those produced by mergers, but also that millisecond pulsar binaries with a high eccentricity may originate from LMXBs that have been involved in, yet not disrupted by, a QN.

The formation of bimodal dust species (that is, both of silicate and amorphous carbon dust grains are observed in a nova eruption) in nova ejecta is still debated. Using the Modules for Experiments in Stellar Astrophysics code and considering the effects of WD's mass, mass-accretion rate and the chemical profiles of WD which are described by new parameter --- mixing depth on the chemical abundances of nova ejecta, we investigate the possibility that bimodal dust species are produced in a nova eruption. We find that $C/O$ (the ratio of carbon number density to oxygen number density) of nova ejecta is affected by the mixing depth. For the model with a small mixing depth, the $C/O$ of nova ejecta can evolute from lager than 1.0 to smaller than 1.0 in a whole eruption, which provides the chemical condition for the formation of bimodal dust species.

R-mode oscillations have been identified as viable and promising targets for continuous gravitational wave searches, meanwhile, it would allow us to probe the interior of compact stars directly. As well as emitting gravitational wave, r-modes would strongly affect the thermal and spin evolution of compact stars. In this paper, we reviewed the theory behind the gravitational wave driven r-mode instability in a rapidly rotating compact star. In particular, we will focus on r-mode instability window, r-mode evolution and detectability of r-mode.

This paper uses population synthesis to investigate the possible origin of isolated millisecond pulsars (IMSPs) as born from the coalescence of a neutron star (NS) and a white dwarf (WD). Results show that the galactic birth - rate of IMSPs is likely to lie between $5.8\times10^{-5}$ $\rm yr^{-1}$ and $2.0\times10^{-4}$ $\rm yr^{-1}$, depending on critical variables, such as the stability of mass transfer via the Roche lobe and the value of kick velocity. In addition to this, this paper estimates that the solar mass of IMSPs can range from $1.5$ and 2.0 $M_{\odot}$, making them more massive than other 'normal' pulsars. Finally, the majority of IMSPs in our simulations have spin periods ranging from several to 20 milliseconds, which is consistent with previous observations.

It is well accepted that a magnetic flux rope (MFR) is a critical component of many coronal mass ejections (CMEs), yet how it evolves toward eruption remains unclear. Here we investigate the continuous evolution of a pre-existing MFR, which is rooted in strong photospheric magnetic fields and electric currents. The evolution of the MFR is observed by the Solar Terrestrial Relations Observatory (STEREO) and the Solar Dynamics Observatory (SDO) from multiple viewpoints. From STEREO's perspective, the MFR starts to rise slowly above the limb five hours before it erupts as a halo CME on 2012 June 14. In SDO observations, conjugate dimmings develop on the disk, simultaneously with the gradual expansion of the MFR, suggesting that the dimmings map the MFR's feet. The evolution comprises a two-stage gradual expansion followed by another stage of rapid acceleration/eruption. Quantitative measurements indicate that magnetic twist of the MFR increases from 1.0 +/- 0.5 to 2.0 +/- 0.5 turns during the five-hour expansion, and further increases to about 4.0 turns per AU when detected as a magnetic cloud at 1 AU two day later. In addition, each stage is preceded by flare(s), implying reconnection is actively involved in the evolution and eruption of the MFR. The implications of these measurements on the CME initiation mechanisms are discussed.

The DarkSide-50 direct-detection dark matter experiment is a dual-phase argon time projection chamber operating at Laboratori Nazionali del Gran Sasso. This paper reports on the blind analysis of a (16,660+-270) kg d exposure using a target of low-radioactivity argon extracted from underground sources. We find no events in the dark matter selection box and set a 90% C.L. upper limit on the dark matter-nucleon spin-independent cross section of 1.14E-44 cm^2 (3.78E-44 cm^2, 3.43E-43 cm^2) for a WIMP mass of 100 GeV/c^2 (1 TeV/c^2, 10 TeV/c^2).

Highly precise pulsar timing is very important for understanding the nature of a neutron star, and it can even be used to detect gravitational waves. Unfortunately, the accuracy of the pulsar timing is seriously affected by the spin-down irregularities of pulsars.

Crater counting on the Moon and other bodies is crucial to constrain the dynamical history of the Solar System. This has traditionally been done by visual inspection of images, thus limiting the scope, efficiency, and/or accuracy of retrieval. In this paper we demonstrate the viability of using convolutional neural networks (CNNs) to determine the positions and sizes of craters from Lunar digital elevation maps (DEMs). We recover 92% of craters from the human-generated test set and almost double the total number of crater detections. Of these new craters, 15% are smaller in diameter than the minimum crater size in the ground-truth dataset. Our median fractional longitude, latitude and radius errors are 11% or less, representing good agreement with the human-generated datasets. From a manual inspection of 361 new craters we estimate the false positive rate of new craters to be 11%. Moreover, our Moon-trained CNN performs well when tested on DEM images of Mercury, detecting a large fraction of craters in each map. Our results suggest that deep learning will be a useful tool for rapidly and automatically extracting craters on various Solar System bodies. We make our code and data publicly available at https://github.com/silburt/DeepMoon.git and https://doi.org/10.5281/zenodo.1133969 .

Highly precise pulsar timing is very important for understanding the nature of neutron stars and the implied physics, even being used to detect directly gravitational waves. Unfortunately, the accuracy of the pulsar timing is seriously affected by the spin-down irregularities, such as spin fluctuations with a manifestation of low frequency structures (the so-called red noise processes), and various activities of magnetospheres. Except from the random timing noise, significant periodic or quasi-periodic components could also be found in some timing residuals of pulsars, indicating the presence of unmodelled deterministic effects. The physical origins of these effects still remain unexplained. In this paper, we suggest a new mechanism involving the de Haas-van Alphen magnetic oscillation, which could trigger the observed low frequency structures. We find that the proposed magnetic oscillation period is about 1-$10^{2}$ yr, which is about $10^{-4}$ times as long as the classical characteristic time scale of interior magnetic field evolution for a normal neutron star. Due to the magnetic oscillation, we estimate both the braking index range between $10^{-5}$ and $10^{5}$ and the residuals range between 4 to 906 ms for some specific samples of quasi-periodic pulsars. Those are consistent with the pulsar timing observations. Nonlinear phenomena of de Haas-van Alphen effect and the Condon domain structure are considered. Avalanche-like magnetized process of domain wall motion (rapidly magnetic energy release) could associated with some special emitting manifestations of neuron star radiation.

The precise observations of Galactic cosmic ray fluxes of the secondary family, such as Li, Be, B, are expected to have significant implications on our understanding of the cosmic ray origin and propagation. Here we employ the recent very precise measurements of those species by the Alpha Magnetic Spectrometer on the International Space Station, together with their parent species (C and O), as well as the data collected by the Voyager-1 spacecraft outside the heliosphere and the Advanced Composition Explorer, we investigate the propagation of cosmic rays in the Milky Way. We find that these new data strongly favor the existence of reacceleration of cosmic rays during the propagation rather than the convective transport. We further find that for the reacceleration model, the slope of rigidity dependence of the diffusion coefficient is close to the Kolmogorov-type of turbulence, and the spectral hardenings of both the primary and secondary particles can be well described by the injection hardening rather than the propagation hardening.

We present new constraints on sub-GeV dark matter particles scattering off electrons in argon based on an analysis of ionization signal data from the DarkSide-50 detector.

We present the results of a search for dark matter WIMPs in the mass range below 20 GeV/c^2 using a target of low-radioactivity argon. The data were obtained using the DarkSide-50 apparatus at Laboratori Nazionali del Gran Sasso (LNGS). The analysis is based on the ionization signal, for which the DarkSide-50 time projection chamber is fully efficient at 0.1 keVee. The observed rate in the detector at 0.5 keVee is about 1.5 events/keVee/kg/day and is almost entirely accounted for by known background sources. We obtain a 90% C.L. exclusion limit above 1.8 GeV/c^2 for the spin-independent cross section of dark matter WIMPs on nucleons, extending the exclusion region for dark matter below previous limits in the range 1.8-6 GeV/c^2.

The acceleration of cosmic ray particles and their propagation in the Milky Way and the heliosphere tangle with each other, leading to complexity and degeneracy of the modeling of Galactic cosmic rays (GCRs). The recent measurements of the GCR spectra by Voyager-1 from outside of the heliosphere gave a first direct observation of GCRs in the local interstellar (LIS) environment. Together with the high-precision data near the Earth taken by ACE and AMS-02, we derive the LIS spectra of Helium, Lithium, Beryllium, Boron, Carbon, and Oxygen nuclei from a few MeV/n to TeV/n, using a non-parameterization method. These LIS spectra are helpful in further studying the injection and propagation parameters of GCRs. The nearly 20 years of data recorded by ACE are used to determine the solar modulation parameters over the solar cycles 23 and 24, based on the force-field approximation. We find general agreements of the modulation potential with the results inferred from neutron monitors and other cosmic ray data.

By a method of population synthesis we construct a model of dark matter (DM) accretion onto binary black holes (BHs), and investigate the merger rate of the binary BHs. We find that the merger rate can weakly increase (less than 10\%). However, the DM accretion can efficiently enhance the masses of binary BHs. In our model, the result for $Z = 0.01$ without the DM accretion can not explain the GW170104, GW170814, and GW150914, while with the DM accretion, it can cover all observations well. For the higher metallicity ($Z = 0.02$), our model can not explain the mergers of high-mass binary BHs like GW170104, GW170814, and GW150914. We estimate that the merger rate of binary BHs lies between $\rm55Gpc^{-3}yr^{-1}$ to $\rm197Gpc^{-3}yr^{-1}$.

Using the stellar evolution code---Modules for Experiments in Stellar Astrophysics (MESA), we investigate the evolution of massive stars with different rotational velocities and metallicities towards Wolf-Rayet stars. In our simulations, the initial rotating velocities are taken as 0, 250, 500 and 650 km s$^{-1}$, and the metallicities equal to 0.02, 0.014, 0.008, 0.006, 0.004 and 0.002. We show our rapid rotation models in the HR diagram compared with the observations. We find that the rotational mixing is less efficient at high metallicity, and these stars become Wolf-Rayet (WR) stars when the helium in their center is ignited. However, rapid rotating massive stars at low metallicity can easily evolve into WR stars due to the rotation resulted in chemically homogeneous evolution. This can explain the origin of single WR stars in galaxy at low metallicity. In our models, the observed SMC WR stars are consistent with the single-star evolution models. However at higher metallicities our single-star evolution models can only explain the luminous, hydrogen-rich WN stars and O stars (which are classified as WR stars previously).

Coronal loops interconnecting two active regions, called as interconnecting loops (ILs), are prominent large-scale structures in the solar atmosphere. They carry a significant amount of magnetic flux, therefore are considered to be an important element of the solar dynamo process. Earlier observations show that eruptions of ILs are an important source of CMEs. It is generally believed that ILs are formed through magnetic reconnection in the high corona (>150-200"), and several scenarios have been proposed to explain their brightening in soft X-rays (SXRs). Yet, the detailed IL formation process has not been fully explored and the associated energy release in the corona still remains unresolved. Here we report the complete formation process of a set of ILs connecting two nearby active regions, with successive observations by STEREO-A on the far side of the Sun and SDO and Hinode on the Earth side. We conclude that ILs are formed by gradual reconnection high in the corona, in line with earlier postulations. In addition, we show evidence supporting that ILs become brightened in SXRs and EUVs through heating at or close to the reconnection site in the corona (i.e., through direct heating process of reconnection), a process that has been largely overlooked in earlier studies on ILs

We analyze and model a C5.7 two-ribbon solar flare observed by SDO, Hinode and GOES on 2011 December 26. The flare is made of many loops formed and heated successively over one and half hours, and their footpoints are brightened in the UV 1600 A before enhanced soft X-ray and EUV missions are observed in flare loops. Assuming that anchored at each brightened UV pixel is a half flaring loop, we identify more than 6,700 half flaring loops, and infer the heating rate of each loop from the UV light curve at the foot-point. In each half loop, the heating rate consists of two phases, an intense impulsive heating followed by a low-rate heating persistent for more than 20 minutes. Using these heating rates, we simulate the evolution of their coronal temperatures and densities with the model of "enthalpy-based thermal evolution of loops" (EBTEL). In the model, suppression of thermal conduction is also considered. This model successfully reproduces total soft X-ray and EUV light curves observed in fifteen pass-bands by four instruments GOES, AIA, XRT, and EVE. In this flare, a total energy of 4.9x10^30 ergs is required to heat the corona, around 40% of this energy is in the slow heating phase. About two fifth of the total energy used to heat the corona is radiated by the coronal plasmas, and the other three fifth transported to the lower atmosphere by thermal conduction.

Using the Modules for Experiments in Stellar Astrophysics code, we investigate the influences of irradiation on ultra-compact X-ray binary (UCXB) evolution. Although the persistent UCXBs have short orbital periods which result in high irradiation flux, the irradiation hardly affects the evolution of persistent sources because the WDs in these binaries have large masses which lead to very low irradiation depth. The irradiation has a significant effect on the transient sources during outburst phase. At the beginning of the outburst, high X-ray luminosity produces high radiation flux, which results in the significant expansion of WD. Then, the irradiation triggers high mass-transfer rates, which can last several days for the transient sources with WDs whose masses are larger than $\sim0.015 M_\odot$ or several hundred years for these sources with WDs whose masses are less than $\sim0.012 M_\odot$. The observed three persistent UCXBs, XTE J0929-314, 4U 1916-05 and SWIFT J1756.9-2508, may belong to the latter.

Pulsar death line can be defined in $P-\dot{P}$ diagram. Traditionally, radio-loud pulsars are supposed to locate above the death line where is the radio-loud region. With the development of observational equipment, the observational properties of the neutron star are remarkably diverse. In $P-\dot{P}$ diagram, some of the special sources are radio-quiet but lie above the death line. From the definition of the pulsar death line, different equation of states for neutron star or strange star results in different death lines. We discuss the influence of the equation of state on the pulsar death line and the possible link between different neutron star groups. The results show that central compact objects would be the small mass of self-bound strange stars, and rotating radio transients might be old pulsars on the verge of death. We suggest that PSR J2144-3933 is likely to be a large mass pulsar, which would be larger than $2.0\rm{M_{\odot}}$. Multiple observational facts would help us to reveal the nature of pulsars.

We report the intriguing large-scale dynamic phenomena associated with the M6.5 flare~(SOL2015-06-22T18:23) in NOAA active region 12371, observed by RHESSI, Fermi, and the Atmospheric Image Assembly (AIA) and Magnetic Imager (HMI) on the Solar Dynamic Observatory (SDO). The most interesting feature of this event is a third ribbon (R3) arising in the decay phase, propagating along a dimming channel (seen in EUV passbands) towards a neighboring sunspot. The propagation of R3 occurs in the presence of hard X-ray footpoint emission, and is broadly visible at temperatures from 0.6~MK to over 10~MK through the Differential Emission Measure (DEM) analysis. The coronal loops then undergo an apparent slipping motion following the same path of R3, after a $\sim$ 80~min delay. To understand the underlying physics, we investigate the magnetic configuration and the thermal structure of the flaring region. Our results are in favor of a slipping-type reconnection followed by the thermodynamic evolution of coronal loops. In comparison with those previously reported slipping reconnection events, this one proceeds across a particularly long distance ($\sim$60 Mm) over a long period of time ($\sim$50 min), and shows two clearly distinguished phases: the propagation of the footpoint brightening driven by nonthermal particle injection and the apparent slippage of loops governed by plasma heating and subsequent cooling.

We investigate the effects of the core-collapse supernova ejecta on a rapidly rotating and massive companion star. We show that the stripped mass raises by twice when compare with a massive but non-rotating companion star. In close binaries with orbital periods of about 1 day, the stripped masses reach up to $\sim 1 M_\odot$. By simulating the evolutions of the rotational velocities of the massive companion stars based on different stripped masses, we find that the rotational velocity decreases greatly for stripped mass that is higher than about $1 M_\odot$. Of all the known high mass X-ray binaries (HMXBs), Cygnus X-3 and 1WGA J0648.024418 have the shortest orbital periods of 0.2 and 1.55 days, respectively. The optical counterpart of the former is a Wolf-Rayet star, whereas it is a hot subdwarf for the latter. Applying our model to the two HMXBs, we suggest that the hydrogen-rich envelopes of their optical counterparts may have been stripped by CCSN ejecta.

Neutron stars are ideal astrophysical laboratories for testing theories of the de Haas-van Alphen (dHvA) effect and diamagnetic phase transition which is associated with magnetic domain formation. The "magnetic interaction" between delocalized magnetic moments of electrons (the Shoenberg effect), can result in an effect of the diamagnetic phase transition into domains of alternating magnetization (Condon's domains). Associated with the domain formation are prominent magnetic field oscillation and anisotropic magnetic stress which may be large enough to fracture the crust of magnetar with a super-strong field. Even if the fracture is impossible as in "low-field" magnetar, the depinning phase transition of domain wall motion driven by low field rate (mainly due to the Hall effect) in the randomly perturbed crust can result in a catastrophically variation of magnetic field. This intermittent motion, similar to the avalanche process, makes the Hall effect be dissipative. These qualitative consequences about magnetized electron gas are consistent with observations of magnetar emission, and especially the threshold critical dynamics of driven domain wall can partially overcome the difficulties of "low-field" magnetar bursts and the heating mechanism of transient, or "outbursting" magnetar.

According to the nova model from \citet{Yaron2005} and \citet{Jose1998} and using Monte Carlo simulation method, we investigate the contribution of chemical abundances in nova ejecta to the interstellar medium (ISM) of the Galaxy. We find that the ejected mass by classical novae (CNe) is about $2.7\times10^{-3}$ $ \rm M_\odot\ {\rm yr^{-1}}$. In the nova ejecta, the isotopic ratios of C, N and O, that is, $^{13}$C/$^{12}$C, $^{15}$N/$^{14}$N and $^{17}$O/$^{16}$O, are higher about one order of magnitude than those in red giants. We estimate that about 10$\%$, 5$\%$ and 20$\%$ of $^{13}$C, $^{15}$N and $^{17}$O in the ISM of the Galaxy come from nova ejecta, respectively. However, the chemical abundances of C, N and O calculated by our model can not cover all of observational values. This means that there is still a long way to go for understanding novae.

In this paper, we investigate the electron Landau-level stability and its influence on the electron Fermi energy, $E_{\rm F}(e)$, in the circumstance of magnetars, which are powered by magnetic field energy. In a magnetar, the Landau levels of degenerate and relativistic electrons are strongly quantized. A new quantity $g_{n}$, the electron Landau-level stability coefficient is introduced. According to the requirement that $g_{n}$ decreases with increasing the magnetic field intensity $B$, the magnetic-field index $\beta$ in the expression of $E_{\rm F}(e)$ must be positive. By introducing the Dirac$-\delta$ function, we deduce a general formulae for the Fermi energy of degenerate and relativistic electrons, and obtain a particular solution to $E_{\rm F}(e)$ in a superhigh magnetic field (SMF). This solution has a low magnetic-field index of $\beta=1/6$, compared with the previous one, and works when $\rho\geq 10^{7}$~g cm$^{-3}$ and $B_{\rm cr}\ll B\leq 10^{17}$~Gauss. By modifying the phase space of relativistic electrons, a SMF can enhance the electron number density $n_e$, and decrease the maximum of electron Landau level number, which results in a redistribution of electrons. According to Pauli exclusion principle, the degenerate electrons will fill quantum states from the lowest Landau level to the highest Landau level. As $B$ increases, more and more electrons will occupy higher Landau levels, though $g_{n}$ decreases with the Landau level number $n$. The enhanced $n_{e}$ in a SMF means an increase in the electron Fermi energy and an increase in the electron degeneracy pressure. The results are expected to facilitate the study of the weak-interaction processes inside neutron stars and the magnetic-thermal evolution mechanism for megnetars.

Liquid argon is a bright scintillator with potent particle identification properties, making it an attractive target for direct-detection dark matter searches. The DarkSide-50 dark matter search here reports the first WIMP search results obtained using a target of low-radioactivity argon. DarkSide-50 is a dark matter detector, using two-phase liquid argon time projection chamber, located at the Laboratori Nazionali del Gran Sasso. The underground argon is shown to contain Ar-39 at a level reduced by a factor (1.4 +- 0.2) x 10^3 relative to atmospheric argon. We report a background-free null result from (2616 +- 43) kg d of data, accumulated over 70.9 live-days. When combined with our previous search using an atmospheric argon, the 90 % C.L. upper limit on the WIMP-nucleon spin-independent cross section based on zero events found in the WIMP search regions, is 2.0 x 10^-44 cm^2 (8.6 x 10^-44 cm^2, 8.0 x 10^-43 cm^2) for a WIMP mass of 100 GeV/c^2 (1 TeV/c^2 , 10 TeV/c^2).

We report multi-wavelength and multi-viewpoint observations of a solar eruptive event which involves loop-loop interactions. During a C2.0 flare, motions associated with inflowing and outflowing plasma provide evidence for ongoing magnetic reconnection. The flare loop top and a rising "concave-up" feature are connected by a current-sheet-like structure (CSLS). The physical properties (thickness, length, temperature, and density) of the CSLS are evaluated. In regions adjacent to the CSLS, the EUV emission (characteristic temperature at 1.6 MK) begins to increase more than ten minutes prior to the onset of the flare, and steeply decreases during the decay phase. The reduction of the emission resembles that expected from coronal dimming. The dynamics of this event imply a magnetic reconnection rate in the range 0.01 -- 0.05.

Accurate numerical solutions of the equations of hydrodynamics play an ever more important role in many fields of astrophysics. In this work, we reinvestigate the accuracy of the moving-mesh code \textsc{Arepo} and show how its convergence order can be improved for general problems. In particular, we clarify that for certain problems \textsc{Arepo} only reaches first-order convergence for its original formulation. This can be rectified by simple modifications we propose to the time integration scheme and the spatial gradient estimates of the code, both improving the accuracy of the code. We demonstrate that the new implementation is indeed second-order accurate under the $L^1$ norm, and in particular substantially improves conservation of angular momentum. Interestingly, whereas these improvements can significantly change the results of smooth test problems, we also find that cosmological simulations of galaxy formation are unaffected, demonstrating that the numerical errors eliminated by the new formulation do not impact these simulations. In contrast, simulations of binary stars followed over a large number of orbital times are strongly affected, as here it is particularly crucial to avoid a long-term build up of errors in angular momentum conservation.

We report on an eruption involving a relatively rare filament-filament interaction on 2013 June 21, observed by SDO and STEREO-B. The two filaments were separated in height with a 'double-decker' configuration. The eruption of the lower filament began simultaneously with a descent of the upper filament resulting in a convergence and direct interaction of the two filaments. The interaction was accompanied by the heating of surrounding plasma and an apparent crossing of a loop-like structure through the upper filament. The subsequent coalescence of the filaments drove a bright front ahead of the erupting structures. The whole process was associated with a C3.0 flare followed immediately by an M2.9 flare. Shrinking loops and descending dark voids were observed during the M2.9 flare at different locations above a C-shaped flare arcade as part of the energy release, giving us a unique insight into the flare dynamics.

As the evolutionary link between the radio millisecond pulsars (MSPs) and the low mass X-ray binaries or intermediate mass X-ray binaries, the millisecond X-ray pulsars (MSXPs) are important objects in testing theories of pulsar formation and evolution. In general, neutron stars in MSXPs can form via core collapse supernova (CC channel) of massive stars or accretion induced collapse (AIC channel) of an accreting ONeMg WD whose mass reaches the Chandrasekhar limit. Here, in addition to CC and AIC channels we also consider another channel, i.e., evolution induced collapse (EIC channel) of a helium star with mass between $1.4$ and $2.5 M_{\odot}$. Using a population synthesis code, we have studied MSXPs arising from three different evolutionary channels. We find that the Galactic birthrates of transient MSXPs and persistent MSXPs are about 0.7---$1.4\times 10^{-4}$ yr$^{-1}$. Our population synthesis calculations have shown that about 50\%---90\% of the MSXPs have undergone CC channel, about 10\%---40\% of them have undergone EIC channel, and the MSXPs via AIC channel are the least.

White dwarf (WD) binary mergers are possible progenitors to a number of unusual stars and transient phenomena, including type Ia supernovae. To date, simulations of mergers have not included magnetic fields, even though they are believed to play a significant role in the evolution of the merger remnant. We simulated a 0.625 - 0.65 $M_{\odot}$ carbon-oxygen WD binary merger in the magnetohydrodynamic moving mesh code Arepo. Each WD was given an initial dipole field with a surface value of $\sim10^3$ G. As in simulations of merging double neutron star binaries, we find exponential field growth within Kelvin-Helmholtz instability-generated vortices during the coalescence of the two stars. The final field has complex geometry, and a strength $>10^{10}$ G at the center of the merger remnant. Its energy is $\sim2\times10^{47}$ ergs, $\sim0.2$% of the remnant's total energy. The strong field likely influences further evolution of the merger remnant by providing a mechanism for angular momentum transfer and additional heating, potentially helping to ignite carbon fusion.

In this paper, four sets of evolutionary models are computed with different values of the mixing length parameter $\alpha_{\rm p}$ and the overshooting parameter $\delta_{\rm ov}$. The properties of the convective cores and the convective envelopes are studied in the massive stars. We get three conclusions: First, the larger $\alpha_{\rm p}$ leads to enhancing the convective mixing, removing the chemical gradient, and increasing the convective heat transfer efficiency. Second, core potential $\phi_{\rm c} = M_{\rm c} / R_{\rm c}$ describes sufficiently the evolution of a star, whether it is a red or blue supergiant at central helium ignition. Third, the discontinuity of hydrogen profile above the hydrogen burning shell seriously affect the occurrence of blue loops in the Hertzsprung--Russell diagram.

The origin of dust in a galaxy is poorly understood. Recently, the surveys of the Large Magellanic Cloud (LMC) provide astrophysical laboratories for the dust studies. By a method of population synthesis, we investigate the contributions of dust produced by asymptotic giant branch (AGB) stars, common envelope (CE) ejecta and type II supernovae (SNe II) to the total dust budget in the LMC. Based on our models, the dust production rates (DPRs) of AGB stars in the LMC are between about $2.5\times10^{-5}$ and $4.0\times10^{-6}M_\odot{\rm yr^{-1}}$. The uncertainty mainly results from different models for the dust yields of AGB stars. The DPRs of CE ejecta are about $6.3\times10^{-6}$(The initial binary fraction is 50\%). These results are within the large scatter of several observational estimates. AGB stars mainly produce carbon grains, which is consistent with the observations. Most of dust grains manufactured by CE ejecta are silicate and iron grains. The contributions of SNe II are very uncertain. Compared with SNe II without reverse shock, the DPRs of AGB stars and CE ejecta are negligible. However, if only 2 \% of dust grains produced by SNe II can survive after reverse shock, the contributions of SNe II are very small. The total dust masses produced by AGB stars in the LMC are between $2.8\times10^4$ and $3.2\times10^5M_\odot$, and those produced by CE ejecta are about $6.3\times10^4$. They are much lower than the values estimated by observations. Therefore, there should be other dust sources in the LMC.

We report the first results of DarkSide-50, a direct search for dark matter operating in the underground Laboratori Nazionali del Gran Sasso (LNGS) and searching for the rare nuclear recoils possibly induced by weakly interacting massive particles (WIMPs). The dark matter detector is a Liquid Argon Time Projection Chamber with a (46.4+-0.7) kg active mass, operated inside a 30 t organic liquid scintillator neutron veto, which is in turn installed at the center of a 1 kt water Cherenkov veto for the residual flux of cosmic rays. We report here the null results of a dark matter search for a (1422+-67) kg d exposure with an atmospheric argon fill. This is the most sensitive dark matter search performed with an argon target, corresponding to a 90% CL upper limit on the WIMP-nucleon spin-independent cross section of 6.1x10^-44 cm^2 for a WIMP mass of 100 GeV/c^2.

It is generally inferred from astronomical measurements that Dark Matter (DM) comprises approximately 27\% of the energy-density of the universe. If DM is a subatomic particle, a possible candidate is a Weakly Interacting Massive Particle (WIMP), and the DarkSide-50 (DS) experiment is a direct search for evidence of WIMP-nuclear collisions. DS is located underground at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy, and consists of three active, embedded components; an outer water veto (CTF), a liquid scintillator veto (LSV), and a liquid argon (LAr) time projection chamber (TPC). This paper describes the data acquisition and electronic systems of the DS detectors, designed to detect the residual ionization from such collisions.

The merger of two carbon-oxygen white dwarfs has long been theorized to lead to a massive carbon-oxygen or oxygen-neon white dwarf, accretion-induced collapse to a neutron star, or a type Ia supernova. Determining which mergers lead to a particular outcome requires hydrodynamic simulations of the merging process. I give a brief overview of the current understanding of mergers and their end-products derived from simulations, and show how temperature, rather than density or mass, most strongly determines a merging binary's subsequent evolution. I then describe recent simulations that show mergers generate strong magnetic fields that could help drive a merger remnant to ignition.

The origin of dust grains in the interstellar medium is still open problem. \cite{Nicholls2013} found the presence of a significant amount of dust around V1309 Sco which maybe originate from the merger of a contact binary. We investigate the origin of dust around V1309 Sco, and suggest that these dust grains are efficiently produced in the binary-merger ejecta. By means of \emph{AGBDUST} code, we estimate that $\sim 5.2\times10^{-4} M_\odot$ of dust grains are produced, and their radii are $\sim 10^{-5}$ cm. These dust grains mainly are composed of silicate and iron grains. Because the mass of the binary-merger ejecta is very small, the contribution of dust produced by binary-merger ejecta to the overall dust production in the interstellar medium is negligible. However, it is the most important that the discovery of a significant amount of dust around V1309 Sco offers a direct support for the idea---common-envelope ejecta provides an ideal environment for dust formation and growth. Therefore, we confirm that common-envelope ejecta can be important source of cosmic dust.

The electron gas inside a neutron star is highly degenerate and relativistic. Due to the electron-electron magnetic interaction, the differential susceptibility can equal or exceed 1, which causes the magnetic system of the neutron star to become metastable or unstable. The Fermi liquid of nucleons under the crust can be in a metastable state, while the crust is unstable to the formation of layers of alternating magnetization. The change of the magnetic stress acting on adjacent domains can result in a series of shifts or fractures in the crust. The releasing of magnetic free energy and elastic energy in the crust can cause the bursts observed in magnetars. Simultaneously, a series of shifts or fractures in the deep crust which is closed to the Fermi liquid of nucleons can trigger the phase transition of the Fermi liquid of nucleons from a metastable state to a stable state. The released magnetic free energy in the Fermi liquid of nucleons corresponds to the giant flares observed in some magnetars.

The merger of two carbon-oxygen white dwarfs can lead either to a spectacular transient, stable nuclear burning or a massive, rapidly rotating white dwarf. Simulations of mergers have shown that the outcome strongly depends on whether the white dwarfs are similar or dissimilar in mass. In the similar-mass case, both white dwarfs merge fully and the remnant is hot throughout, while in the dissimilar case, the more massive, denser white dwarf remains cold and essentially intact, with the disrupted lower mass one wrapped around it in a hot envelope and disk. In order to determine what constitutes "similar in mass" and more generally how the properties of the merger remnant depend on the input masses, we simulated unsynchronized carbon-oxygen white dwarf mergers for a large range of masses using smoothed-particle hydrodynamics. We find that the structure of the merger remnant varies smoothly as a function of the ratio of the central densities of the two white dwarfs. A density ratio of 0.6 approximately separates similar and dissimilar mass mergers. Confirming previous work, we find that the temperatures of most merger remnants are not high enough to immediately ignite carbon fusion. During subsequent viscous evolution, however, the interior will likely be compressed and heated as the disk accretes and the remnant spins down. We find from simple estimates that this evolution can lead to ignition for many remnants. For similar-mass mergers, this would likely occur under sufficiently degenerate conditions that a thermonuclear runaway would ensue.

The material that is ejected in a common-envelope (CE) phase in a close binary system provides an ideal environment for dust formation. By constructing a simple toy model to describe the evolution of the density and the temperature of CE ejecta and using the \emph{AGBDUST} code to model dust formation, we show that dust can form efficiently in this environment. The actual dust masses produced in the CE ejecta depend strongly on their temperature and density evolution. We estimate the total dust masses produced by CE evolution by means of a population synthesis code and show that, compared to dust production in AGB stars, the dust produced in CE ejecta may be quite significant and could even dominate under certain circumstances.

Strange quark stars (SSs) may originate from accreting neutron stars (NSs) in low-mass X-ray binaries (LMXBs). Assuming that conversion of NS matter to SSs occurs when the core density of accreting NS reaches to the density of quark deconfinement, $\sim 5 \rho_0$, where $\rho_0\sim 2.7\times 10^{14}$g cm$^{-3}$ is nuclear saturation density, we investigate LMXBs with SSs (qLMXBs). In our simulations, about 1\permil --- 10% of LMXBs can produce SSs, which greatly depends on the masses of nascent NSs and the fraction of transferred matter accreted by the NSs. If the conversion does not affect binaries systems, LMXBs evolve into qLMXBs. We find that some observational properties (spin periods, X-ray luminosities and orbital periods) of qLMXBs are similar with those of LMXBs, and it is difficult to differ them. If the conversion disturbs the binaries systems, LMXBs can produce isolated SSs. These isolated SSs could be submillisecond pulsars, and their birthrate in the Galaxy is $\sim$5--70 per Myr.

Properties of X-ray luminosities in low-mass X-ray binaries (LMXBs) mainly depend on donors. We have carried out a detailed study of donors in persistent neutron-star LMXBs (PLMXBs) by means of a population synthesis code. PLMXBs with different donors have different formation channels. Our numerical simulations show that more than 90% of PLMXBs have main sequence (MS) donors, and PLMXBs with red giant (RG) donors via stellar wind (Wind) are negligible. In our model, most of neutron stars (NSs) in PLMXBs with hydrogen-rich donors form via core-collapse supernovae, while more than 90% of NSs in PLMXBs with naked helium star (He) donors or white dwarf (WD) donors form via an evolution-induced collapse via helium star ($1.4 \leq M_{\rm He}/M_\odot \leq 2.5$) or an accretion-induced collapses for an accreting ONeMg WD.

Symbiotic X-ray binaries (SyXBs) comprise a rare class of low-mass X-ray binaries. We study the Galactic SyXBs, which we consider as detached binaries composed of low-mass giants and wind-fed neutron star companions, by simulation of the interaction of a magnetized neutron star (NS) with its environment and utilizing a population synthesis code. We focus mainly on the parameters that influence observational appearance of the SyXB: the donor wind velocity (vw) and the angular momentum distribution in the shell of matter settling onto NS. We estimate the birthrate of SyXB as $\sim 4.1\times 10^{-5}$ yr$^{-1}$ to $ \sim 6.6\times 10^{-6}$ yr$^{-1}$ and their number in the Galaxy as $\sim$(100 -- 1000). Assumed stellar wind velocity from cool giants is the input parameter that influences the model SyXBs population most. Among known SyXBs or candidate systems, 4U 1954+31 and IGR J16358-4724 in which NS have very long spin periods may host quasi-spherically accreting NSs. GX 1+4 has a peculiar long-term spin behaviour and it may also be a quasi-spherical wind-accreting source. We cannot identify whether there are wind-fed accretion disks in 4U 1700+24, Sct X-1, IRXS J180431.1-273932 and 2XMM J174016.0-290337.

Recent studies have shown that for suitable initial conditions both super- and sub-Chandrasekhar mass carbon-oxygen white dwarf mergers produce explosions similar to observed SNe Ia. The question remains, however, how much fine tuning is necessary to produce these conditions. We performed a large set of SPH merger simulations, sweeping the possible parameter space. We find trends for merger remnant properties, and discuss how our results affect the viability of our recently proposed sub-Chandrasekhar merger channel for SNe Ia.

Using a simple accelerating model and an assumption that $\gamma$-rays originate from $p-p$ collisions for a $\pi^0$ model, we investigate $\gamma$-ray sources like V407 Cygni in symbiotic stars. The upper limit of their occurrence rate in the Galaxy is between 0.5 and 5 yr$^{-1}$, indicating that they may be an important source of the high-energy $\gamma$-rays. The maximum energies of the accelerated protons mainly distribute around $10^{11}$ eV, and barely reach $10^{15}$ eV. The novae occurring in D-type SSs with ONe WDs and long orbital periods are good candidates for $\gamma$-ray sources. Due to a short orbital period which results in a short acceleration duration, the nova occurring in symbiotic star RS Oph can not produce the $\gamma$-ray emission like that in V407 Cygni.

We employ the DAVE (differential affine velocity estimator, Schuck 2005; 2006) tracking technique on a time series of MDI/1m high spatial resolution line- of-sight magnetograms to measure the photospheric flow velocity for three newly emerging bipolar active regions. We separately calculate the magnetic helicity injection rate of the leading and following polarities to confirm or refute the magnetic helicity asymmetry, found by Tian & Alexander (2009) using MDI/96m low spatial resolution magnetograms. Our results demonstrate that the magnetic helicity asymmetry is robust being present in the three active regions studied, two of which have an observed balance of the magnetic flux. The magnetic helicity injection rate measured is found to depend little on the window size selected, but does depend on the time interval used between the two successive magnetograms tracked. It is found that the measurement of the magnetic helicity injection rate performs well for a window size between 12x10 and 18x15 pixels, and at a time interval {\Delta}t=10 minutes. Moreover, the short-lived magnetic structures, 10-60 minutes, are found to contribute 30-50% of the magnetic helicity injection rate. Comparing with the results calculated by MDI/96m data, we find that the MDI/96m data, in general, can outline the main trend of the magnetic properties, but they significantly underestimate the magnetic flux in strong field region and are not appropriate for quantitative tracking studies, so provide a poor estimate of the amount of magnetic helicity injected into the corona.

Assuming that soft X-ray sources in symbiotic stars result from strong thermonuclear runaways, and supersoft X-ray sources from weak thermonuclear runaways or steady hydrogen burning symbiotic stars, we investigate the Galactic soft and supersoft X-ray sources in symbiotic stars by means of population synthesis. The Galactic occurrence rates of soft X-ray sources and supersoft X-ray sources are from $\sim$ 2 to 20 $\rm yr^{-1}$, and $\sim$ 2 to 17 $\rm yr^{-1}$, respectively. The numbers of X-ray sources in symbiotic stars range from 2390 to 6120. We simulate the distribution of X-ray sources over orbital periods, masses and mass-accretion rates of white dwarfs. The agreement with observations is reasonable.

By assuming an aspherical stellar wind with an equatorial disk from a red giant, we investigate the production of Type Ia supernovae (SNe Ia) via symbiotic channel. We estimate that the Galactic birthrate of SNe Ia via symbiotic channel is between $1.03\times 10^{-3}$ and $2.27\times 10^{-5}$ yr$^{-1}$, the delay time of SNe Ia has wide range from $\sim$ 0.07 to 5 Gyr. The results are greatly affected by the outflow velocity and mass-loss rate of the equatorial disk. Using our model, we discuss the progenitors of SN 2002ic and SN 2006X.

We have carried out a study of the chemical abundances of $^1$H, $^4$He, $^{12}$C, $^{13}$C, $^{14}$N, $^{15}$N, $^{16}$O, $^{17}$O, $^{20}$Ne and $^{22}$Ne in symbiotic stars (SSs) by means of a population synthesis code. We find that the ratios of the number of O-rich SSs to that of C-rich SSs in our simulations are between 3.4 and 24.1, depending on the third dredge-up efficiency $\lambda$ and the terminal velocity of the stellar wind $v(\infty)$. The fraction of SSs with $extrinsic$ C-rich cool giants in C-rich cool giants ranges from 2.1% to 22.7%, depending on $\lambda$, the common envelope algorithm and the mass-loss rate. Compared with the observations, the distributions of the relative abundances of $^{12}$C/$^{13}$C vs. [C/H] of the cool giants in SSs suggest that the thermohaline mixing in low-mass stars may exist. The distributions of the relative abundances of C/N vs. O/N, Ne/O vs. N/O and He/H vs. N/O in the symbiotic nebulae indicate that it is quite common that the nebular chemical abundances in SSs are modified by the ejected materials from the hot components. Helium overabundance in some symbiotic nebulae may be relevant to a helium layer on the surfaces of white dwarf accretors.