Papers for Thursday, Jan 21 2021

The dominant site of production of $r$-process elements remains unclear despite recent observations of a neutron star merger. Observational constraints on the properties of the sites can be obtained by comparing $r$-process abundances in different environments. The recent Gaia data releases and large samples from high-resolution optical spectroscopic surveys are enabling us to compare $r$-process element abundances between stars formed in an accreted dwarf galaxy, Gaia-Enceladus, and those formed in the Milky Way. We aim to understand the origin of $r$-process elements in Gaia-Enceladus. We first construct a sample of stars to study Eu abundances without being affected by the detection limit. We then kinematically select 71 Gaia-Enceladus stars and 93 in-situ stars from the Galactic Archaeology with HERMES (GALAH) DR3, of which 50 and 75 stars can be used to study Eu reliably. Gaia-Enceladus stars clearly show higher ratios of [{Eu}/{Mg}] than in-situ stars. High [{Eu}/{Mg}] along with low [{Mg}/{Fe}] are also seen in relatively massive satellite galaxies such as the LMC, Fornax, and Sagittarius dwarfs. On the other hand, unlike these galaxies, Gaia-Enceladus does not show enhanced [{Ba}/{Eu}] or [{La}/{Eu}] ratios suggesting a lack of significant $s$-process contribution. From comparisons with simple chemical evolution models, we show that the high [{Eu}/{Mg}] of Gaia-Enceladus can naturally be explained by considering $r$-process enrichment by neutron-star mergers with delay time distribution that follows a similar power-law as type~Ia supernovae but with a shorter minimum delay time.

Since the first signal in 2015, the gravitational-wave (GW) detections of merging binary black holes (BBHs) by the LIGO and Virgo collaborations (LVC) have completely transformed our understanding of the lives and deaths of compact object binaries, and have motivated an enormous amount of theoretical and phenomenological work on the astrophysical origin of these objects. In this work, we show that the phenomenological fit to the redshift-dependent merger rate of BBHs from Abbott et al. (2020) is consistent with a purely dynamical origin for these objects, and that the current merger rate of BBHs from the LVC could be explained entirely with globular clusters alone. While this does not prove that globular clusters (GCs) are the only (or even dominant) formation channel, we emphasize that many potential formation scenarios could contribute a significant fraction of the current LVC rate, and that any analysis that assumes a single (or dominant) mechanism for producing BBH mergers is implicitly using a specious astrophysical prior.

Cloud-cloud collision (CCC) has been suggested as a mechanism to induce massive star formation. Recent simulations suggest that a CCC speed is different among galactic-scale environments, which is responsible for observed differences in star formation activity. In particular, a high-speed CCC is proposed as a cause of star formation suppression in the bar regions in barred spiral galaxies. Focusing on the strongly barred galaxy NGC1300, we investigate the CCC speed. We find the CCC speed in the bar and bar-end tend to be higher than that in the arm. The estimated CCC speed is $\sim20~\rm km~s^{-1}$, $\sim16~\rm km~s^{-1}$, and $\sim11~\rm km~s^{-1}$ in the bar, bar-end, and arm, respectively. Although the star formation activity is different in the bar and bar-end, the CCC speed and the number density of high-speed CCC with $> 20~\rm km~s^{-1}$ are high in both regions, implying the existence of other parameters that control the star formation. The difference in molecular gas mass (average density) of the giant molecular clouds (GMCs) between the bar (lower mass and lower density) and bar-end (higher mass and higher density) may be cause for the different star formation activity. Combining with our previous study (Maeda et al.), the leading candidates of causes for the star formation suppression in the bar in NGC1300 are the presence of a large amount of diffuse molecular gases and high-speed CCCs between low mass GMCs.

Ultra Diffuse Galaxies (UDGs), a type of large Low Surface Brightness (LSB) galaxies with particularly large effective radii (r_eff > 1.5 kpc), are now routinely studied in the local (z<0.1) universe. While they are found to be abundant in clusters, groups, and in the field, their formation mechanisms remain elusive and an active topic of debate. New insights may be found by studying their counterparts at higher redshifts (z>1.0), even though cosmological surface brightness dimming makes them particularly diffcult to detect and study there. This work uses the deepest Hubble Space Telescope (HST) imaging stacks of z > 1 clusters, namely: SPT-CL J2106-5844 and MOO J1014+0038. These two clusters, at z=1.13 and z=1.23, were monitored as part of the HST See-Change program. Compared to the Hubble Extreme Deep Field (XDF) as reference field, we find statistical over-densities of large LSB galaxies in both clusters. Based on stellar population modelling and assuming no size evolution, we find that the faintest sources we can detect are about as bright as expected for the progenitors of the brightest local UDGs. We find that the LSBs we detect in SPT-CL J2106-5844 and MOO J1014-5844 already have old stellar populations that place them on the red sequence. Correcting for incompleteness, and based on an extrapolation of local scaling relations, we estimate that distant UDGs are relatively under-abundant compared to local UDGs by a factor ~ 3. A plausible explanation for the implied increase with time would be a significant size growth of these galaxies in the last ~ 8 Gyr, as also suggested by hydrodynamical simulations.

We quantify two main pathways through which baryonic physics biases cluster count cosmology. We create mock cluster samples that reproduce the baryon content inferred from X-ray observations. The clusters are linked to their counterparts in a dark matter-only universe, whose abundances can be predicted robustly, by assuming the shape of the dark matter density profile does not change significantly due to baryons. We derive halo masses from different weak lensing fitting methods and infer the best-fitting cosmological parameters $\Omega_\mathrm{m}$, $S_8=\sigma_8(\Omega_\mathrm{m}/0.3)^{0.2}$, and $w_0$ from the mock cluster sample. We find that because of the need to accommodate the change in the density profile due to the ejection of baryons, weak lensing mass calibrations are only unbiased if the concentration is left free when fitting the reduced shear with NFW profiles. However, even unbiased total mass estimates give rise to biased cosmological parameters if the measured mass functions are compared with predictions from dark matter-only simulations. This is the dominant bias for haloes with $m_\mathrm{500c} < 10^{14.5} \, h^{-1} \, \mathrm{M}_\odot$. For a stage IV-like cluster survey with area $\approx 15000 \, \mathrm{deg^2}$ and a constant mass cut of $m_\mathrm{200m,min} = 10^{14} \, h^{-1} \, \mathrm{M}_\odot$, the biases are $-11 \pm 1 \, \%$ in $\Omega_\mathrm{m}$, $-3.29 \pm 0.04 \, \%$ in $S_8$, and $9 \pm 1.5 \, \%$ in $w_0$. These systematic biases exceed the statistical uncertainties by factors of 11, 82, and 6, respectively. We suggest that rather than the total halo mass, the (re-scaled) dark matter mass inferred from the combination of weak lensing and observations of the hot gas, should be used for cluster count cosmology.

Massive star-forming regions exhibit an extremely rich and diverse chemistry, which in principle provides a wealth of molecular probes, as well as laboratories for interstellar prebiotic chemistry. Since the chemical structure of these sources displays substantial spatial variation among species on small scales (${\lesssim}10^4$ au), high angular resolution observations are needed to connect chemical structures to local environments and inform astrochemical models of massive star formation. To address this, we present ALMA 1.3 mm observations toward OB cluster-forming region G10.6-0.4 (hereafter "G10.6") at a resolution of 0.14$^{\prime\prime}$ (700 au). We find highly-structured emission from complex organic molecules (COMs) throughout the central 20,000 au, including two hot molecular cores and several shells or filaments. We present spatially-resolved rotational temperature and column density maps for a large sample of COMs and warm gas tracers. These maps reveal a range of gas substructure in both O- and N-bearing species. We identify several spatial correlations that can be explained by existing models of COM formation, including NH$_2$CHO/HNCO and CH$_3$OCHO/CH$_3$OCH$_3$, but also observe unexpected distributions and correlations which suggest that our current understanding of COM formation is far from complete. Importantly, complex chemistry is observed throughout G10.6, rather than being confined to hot cores. The COM composition appears to be different in the cores compared to the more extended structures, which illustrates the importance of high spatial resolution observations of molecular gas in elucidating the physical and chemical processes associated with massive star formation.

Anomalies in the flux-ratios of the images of quadruply-lensed quasars have been used to constrain the nature of dark matter. Assuming these lensing perturbations are caused by dark matter haloes, it is possible to constrain the mass of a hypothetical Warm Dark Matter (WDM) particle to be $m_\chi > 5.2$ keV. However, the assumption that perturbations are only caused by DM haloes might not be correct as other structures, such as filaments and pancakes, exist and make up a significant fraction of the mass in the universe, ranging between 5$\%$ -- 50$\%$ depending on the dark matter model. Using novel fragmentation-free simulations of 1 and 3keV WDM cosmologies we study these "non-halo" structures and estimate their impact on flux-ratio observations. We find that these structures display sharp density gradients with short correlation lengths, and can contribute more to the lensing signal than all haloes up to the half-mode mass combined, thus reducing the differences expected among WDM models. We estimate that this becomes especially important for any flux-ratio based constraint sensitive to haloes of mass $M \sim 10^8 M_\odot$. We conclude that accounting for all types structures in strong-lensing observations is required to improve the accuracy of current and future constraints.

The gas content of the complete compilation of Local Group dwarf galaxies (119 within 2 Mpc) is presented using HI survey data. Within the virial radius of the Milky Way (224 kpc here), 53 of 55 dwarf galaxies are devoid of gas to limits of M$_{\rm HI}<10^4$ M$_\odot$. Within the virial radius of M31 (266 kpc), 27 of 30 dwarf galaxies are devoid of gas (with limits typically $<10^5$ M$_\odot$). Beyond the virial radii of the Milky Way and M31, the majority of the dwarf galaxies have detected HI gas and have HI masses higher than the limits. When the relationship between gas content and distance is investigated using a Local Group virial radius, more of the non-detected dwarf galaxies are within this radius (85$\pm1$ of the 93 non-detected dwarf galaxies) than within the virial radii of the Milky Way and M31. Using the Gaia proper motion measurements available for 38 dwarf galaxies, the minimum gas density required to completely strip them of gas is calculated. Halo densities between $10^{-5}$ and $5 \times 10^{-4}$ cm$^{-3}$ are typically required for instantaneous stripping at perigalacticon. When compared to halo density with radius expectations from simulations and observations, 80% of the dwarf galaxies with proper motions are consistent with being stripped by ram pressure at Milky Way pericenter. The results suggest a diffuse gaseous galactic halo medium is important in quenching dwarf galaxies, and that a Local Group medium also potentially plays a role.

Joint analyses of small-scale cosmological structure probes are relatively unexplored and promise to advance measurements of microphysical dark matter properties using heterogeneous data. Here, we present a multidimensional analysis of dark matter substructure using strong gravitational lenses and the Milky Way (MW) satellite galaxy population, accounting for degeneracies in model predictions and using covariances in the constraining power of these individual probes for the first time. We simultaneously infer the projected subhalo number density and the half-mode mass describing the suppression of the subhalo mass function in thermal relic warm dark matter (WDM), $M_{\mathrm{hm}}$, using the semi-analytic model $\texttt{Galacticus}$ to connect the subhalo population inferred from MW satellite observations to the strong lensing host halo mass and redshift regime. Combining MW satellite and strong lensing posteriors in this parameter space yields $M_{\mathrm{hm}}<10^{7.0}\ \mathrm{M}_{\mathrm{\odot}}$ (WDM particle mass $m_{\mathrm{WDM}}>9.7\ \mathrm{keV}$) at $95\%$ confidence and disfavors $M_{\mathrm{hm}}=10^{7.4}\ \mathrm{M}_{\mathrm{\odot}}$ ($m_{\mathrm{WDM}}=7.4\ \mathrm{keV}$) with a 20:1 marginal likelihood ratio, improving limits on $m_{\mathrm{WDM}}$ set by the two methods independently by $\sim 30\%$. These results are marginalized over the line-of-sight contribution to the strong lensing signal, the mass of the MW host halo, and the efficiency of subhalo disruption due to baryons, and are robust to differences in the disruption efficiency between the MW and strong lensing regimes at the $\sim 10\%$ level. This work paves the way for unified analyses of next-generation small-scale structure measurements covering a wide range of scales and redshifts.

We study primordial black holes (PBHs) formation in the excursion set theory (EST) in a vast range of PBHs masses with and without confirmed constraints on their abundance. In this work, a new concept of the first touch in the EST is introduced for PBHs formation which takes into account the earlier horizon reentry of smaller masses. Our study shows that in the EST, it is possible to produce PBHs in different mass ranges which could make up all dark matter. We also show that in a broad blue-tilted power spectrum, the production of PBHs is dominated towards a smaller mass. Our analysis put upper limit $\sim\,$0.1 on the amplitude of the curvature power spectrum at length scales relevant for PBHs formation.

We present the measured gas-phase metal column densities in 155 sub-damped Lyman alpha systems (subDLAs) with the aim to investigate the contribution of subDLAs to the chemical evolution of the Universe. The sample was identified within the absorber-blind XQ-100 quasar spectroscopic survey over the redshift range 2.4<=z<=4.3. Using all available column densities of the ionic species investigated (mainly CIV, SiII, MgII, SiIV, AlII, FeII, CII, and OI; in order of decreasing detection frequency), we estimate the ionization-corrected gas-phase metallicity of each system using Markov Chain Monte Carlo techniques to explore a large grid of Cloudy ionization models. Without accounting for ionization and dust depletion effects, we find that the HI-weighted gas-phase metallicity evolution of subDLAs are consistent with damped Lyman alpha systems (DLAs). When ionization corrections are included, subDLAs are systematically more metal-poor than DLAs (between ~0.5 sigma and ~3 sigma significance) by up to ~1.0 dex over the redshift range 3<=z<=4.3. The correlation of gas-phase [Si/Fe] with metallicity in subDLAs appears to be consistent with that of DLAs, suggesting that the two classes of absorbers have a similar relative dust depletion pattern. As previously seen for Lyman limit systems, the gas-phase [C/O] in subDLAs remains constantly solar for all metallicities indicating that both subDLAs and Lyman limit systems could trace carbon-rich ejecta, potentially in circumgalactic environments.

We present the analysis of the diffuse, low column density HI environment of 18 MHONGOOSE galaxies. We obtained deep observations with the Robert C. Byrd Green Bank Telescope, and reached down to a 3sigma column density detection limit of NHI=6.3x10^{17} cm^{-2} over a 20 km/s linewidth. We analyze the environment around these galaxies, with a focus on HI gas that reaches column densities below NHI=10^{19} cm^{-2}. We calculate the total amount of HI gas in and around the galaxies revealing that nearly all of these galaxies contained excess HI outside of their disks. We quantify the amount of diffuse gas in the maps of each galaxy, defined by HI gas with column densities below 10^{19} cm^{-2}, and find a large spread in percentages of diffuse gas. However, by binning the percentage of diffuse HI into quarters, we find that the bin with the largest number of galaxies is the lowest quartile (0-25\% diffuse HI). We identified several galaxies which may be undergoing gas accretion onto the galaxy disk using multiple methods of analysis, including azimuthally averaging column densities beyond the disk, and identifying structure within our integrated intensity (Moment 0) maps. We measured HI mass outside the disks of most of our galaxies, with rising cumulative flux even at large radii. We also find a strong correlation between the fraction of diffuse gas in a galaxy and its baryonic mass, and test this correlation using both Spearman and Pearson correlation coefficients. We see evidence of a dark matter halo mass threshold of M_{halo}~10^{11.1} \msun{} in which galaxies with high fractions of diffuse HI all reside below. It is in this regime in which cold-mode accretion should dominate. Finally, we suggest a rotation velocity of v_{rot}~80 km\s as an upper threshold to find diffuse gas-dominated galaxies.

The evolution of circumstellar discs is highly influenced by their surroundings, in particular by external photoevaporation due to nearby stars and dynamical truncations. The impact of these processes on disc populations depends on the dynamical evolution of the star-forming region. Here we implement a simple model of molecular cloud collapse and star formation to obtain primordial positions and velocities of young stars and follow their evolution in time, including that of their circumstellar discs. Our disc model takes into account viscous evolution, internal and external photoevaporation, dust evolution, and dynamical truncations. The disc evolution is resolved simultaneously with the star cluster dynamics and stellar evolution. Our results show that an extended period of star formation allows for massive discs formed later in the simulations to survive for several million years. This could explain massive discs surviving in regions of high UV radiation.

Brightness-weighted differential source counts $S^2 n(S)$ spanning the eight decades of flux density between $0.25\,\mu\mathrm{Jy}$ and 25 Jy at 1.4 GHz were measured from (1) the confusion brightness distribution in the MeerKAT DEEP2 image below $10\,\mu\mathrm{Jy}$, (2) counts of DEEP2 sources between $10\,\mu\mathrm{Jy}$ and $2.5\,\mathrm{mJy}$, and (3) counts of NVSS sources stronger than $2.5\,\mathrm{mJy}$. We present our DEEP2 catalog of $1.7 \times 10^4$ discrete sources complete above $S = 10\,\mu\mathrm{Jy}$ over $\Omega = 1.04\,\mathrm{deg}^2$. The brightness-weighted counts converge as $S^2 n(S) \propto S^{1/2}$ below $S = 10\,\mu\mathrm{Jy}$, so $>99\%$ of the $\Delta T_\mathrm{b} \sim 0.06\,\mathrm{K}$ sky brightness produced by active galactic nuclei and $\approx96\%$ of the $\Delta T_\mathrm{b} \sim 0.04\,\mathrm{K}$ added by star-forming galaxies has been resolved into sources with $S \geq 0.25\,\mu\mathrm{Jy}$. The $\Delta T_\mathrm{b} \approx 0.4\,\mathrm{K}$ excess brightness measured by ARCADE 2 cannot be produced by faint sources smaller than $\approx 50\,\mathrm{kpc}$ if they cluster like galaxies.

High-speed downflows have been observed in the solar transition region (TR) and lower corona for many decades. Despite their abundance, it has been hard to find signatures of such downflows in the solar chromosphere. In this work, we target an enhanced network region that shows ample occurrences of rapid spicular downflows in the \halpha\ spectral line that could potentially be linked to high-speed TR downflowing counterparts. We used the $k$-means algorithm to classify the spectral profiles of on-disk spicules in \halpha{} and \cak{} data observed from the Swedish 1-m Solar Telescope (SST) and employed an automated detection method based on advanced morphological image processing operations to detect such downflowing features, in conjunction with rapid blue-shifted and red-shifted excursions (RBEs and RREs). We report the existence of a new category of RREs (termed as downflowing RRE) for the first time that, contrary to earlier interpretation, are associated with chromospheric field-aligned downflows moving towards the strong magnetic field regions. Statistical analysis performed on nearly 20,000 RBEs and 15,000 RREs (including the downflowing counterparts), detected in our 97~min long dataset, shows that the downflowing RREs are very similar to RBEs and RREs except for their oppositely directed plane-of-sky motion. Furthermore, we also find that RBEs, RREs and downflowing RREs can be represented by a wide range of spectral profiles with varying Doppler offsets, and \halpha{} line core widths, both along and perpendicular to the spicule axis, that causes them to be associated with multiple substructures that evolve together. We speculate that these rapid plasma downflows could well be the chromospheric counterparts of the commonly observed TR downflows.

LISA-Pathfinder is an ESA space mission flown between 2015 and 2017 to demonstrate a technological maturity sufficient for building a gravitational waves telescope in space, such as the Laser Interferometer Space Antenna (LISA). A pair of cubic test masses is hosted inside the LISA-Pathfinder spacecraft and shielded from any force other than the interplanetary gravitational field. The purity of the shielding gives the performance of the mission. There are a number of aspects that had to be confirmed in-flight. One of them is the transition phase from the launch configuration, when the test masses are locked, to the science free-falling configuration. Each test mass is initially released from the mechanical constraints via a dedicated mechanism and then captured by an electrostatic control system. In fact, each test mass is surrounded by a set of electrodes for actuation and sensing purposes. The performance criterion of the release is the final velocity of the test mass relative to the spacecraft, with an upper threshold set to 5 $\mu m/s$. The LISA-Pathfinder first in-flight release velocities highlighted an unexpected dynamics with large linear and angular velocities. The electrostatic control was successful, but only relying on a manual procedure that cannot be considered as baseline for LISA. This paper helps investigating the in-flight non-compliance by dealing with the modeling of the electrostatic environment around each test mass and its contribution to the release and capture dynamics. The electrostatic model is based on the method of moments, a boundary element numerical technique suitable for estimating forces and capacitances between conductors. We also provide a short overview of the method, which can be used for the analysis of other phenomena within LISA and for the design of future gravitational waves telescopes and space projects.

We introduce a new binary detection technique, Binary INformation from Open Clusters using SEDs (binocs), which we show is able to determine reliable stellar multiplicity and masses over a much larger mass range than current approaches. This new technique determines accurate component masses of binary and single systems of the open clusters main sequence by comparing observed magnitudes from multiple photometric filters to synthetic star spectral energy distributions (SEDs) allowing systematically probing the binary population for low mass stars in clusters for 8 well-studied open clusters. We provide new deep, infrared photometric catalogs (1.2 - 8.0 microns) for the key open clusters NGC 1960 (M36), NGC 2099 (M37), NGC 2420, and NGC2682 (M67), using observation from NOAO/NEWFIRM and Spitzer}/IRAC. Using these deep multi-wavelength catalogs, the binocs method is applied to these clusters to determine accurate component masses for unresolved cluster binaries. We explore binary fractions as a function of cluster age, Galactic location and metallicity.

We present and analyse 120 spectroscopic binary and triple cluster members of the old (4 Gyr) open cluster M67 (NGC 2682). As a cornerstone of stellar astrophysics, M67 is a key cluster in the WIYN Open Cluster Study (WOCS); radial-velocity (RV) observations of M67 are ongoing and extend back over 45 years, incorporating data from seven different telescopes, and allowing us to detect binaries with orbital periods <~10^4 days. Our sample contains 1296 stars (604 cluster members) with magnitudes of 10 <= V <= 16.5 (about 1.3 to 0.7 Msolar), from the giants down to ~4 mag below the main-sequence turnoff, and extends in radius to 30 arcminutes (7.4 pc at a distance of 850 pc, or ~7 core radii). This paper focuses primarily on the main-sequence binaries, but orbital solutions are also presented for red giants, yellow giants and sub-subgiants. Out to our period detection limit and within our magnitude and spatial domain, we find a global main-sequence incompleteness-corrected binary fraction of 34% +/- 3%, which rises to 70% +/- 17% in the cluster center. We derive a tidal circularization period of P_circ = 11.0 +1.1 -1.0 days. We also analyze the incompleteness-corrected distributions of binary orbital elements and masses. The period distribution rises toward longer periods. The eccentricity distribution, beyond P_circ, is consistent with a uniform distribution. The mass-ratio distribution is also consistent with a uniform distribution. Overall, these M67 binaries are closely consistent with similar binaries in the galactic field, as well as the old (7 Gyr) open cluster NGC 188. WIYN Open Cluster Study. 83.

In the search for life in the cosmos, NASA's Transiting Exoplanet Survey Satellite (TESS) mission has already monitored about 74% of the sky for transiting extrasolar planets, including potentially habitable worlds. However, TESS only observed a fraction of the stars long enough to be able to find planets like Earth. We use the primary mission data - the first two years of observations - and identify 4,239 stars within 210pc that TESS observed long enough to see 3 transits of an exoplanet that receives similar irradiation to Earth: 738 of these stars are located within 30pc. We provide reliable stellar parameters from the TESS Input Catalog that incorporates Gaia DR2 and also calculate the transit depth and radial velocity semi-amplitude for an Earth-analog planet. Of the 4,239 stars in the Revised TESS HZ Catalog, 9 are known exoplanet hosts - GJ 1061, GJ 1132, GJ 3512, GJ 685, Kepler-42, LHS 1815, L98-59, RR Cae, TOI 700 - around which TESS could identify additional Earth-like planetary companions. 37 additional stars host yet unconfirmed TESS Objects of Interest: three of these orbit in the habitable zone - TOI 203, TOI 715, and TOI 2298. For a subset of 614 of the 4,239 stars, TESS has observed the star long enough to be able to observe planets throughout the full temperate, habitable zone out to the equivalent of Mars' orbit. Thus, the Revised TESS Habitable Zone Catalog provides a tool for observers to prioritize stars for follow-up observation to discover life in the cosmos. These stars are the best path towards the discovery of habitable planets using the TESS mission data.

In the manuscript, we report evidence on broad \oiii components apparently obscured in Type-2 AGN under the framework of the Unified model, after checking properties of broad \oiii emissions in large samples of Type-1 and Type-2 AGN in SDSS DR12. We can well confirm the statistically lower flux ratios of the broad to the core \oiii components in Type-2 AGN than in Type-1 AGN, which can be naturally explained by stronger obscured broad \oiii components by central dust torus in Type-2 AGN, unless the Unified model for AGN was not appropriate to the narrow emission lines. The results provide further evidence to support broad \oiii components coming from emission regions nearer to central BHs, and also indicate the core \oiii component as the better indicator for central activities in Type-2 AGN, due to few effects of obscuration on the core \oiii component. Considering the broad \oiii components as signs of central outflows, the results provide evidence for strong central outflows being preferentially obscured in Type-2 AGN. Furthermore, the obscured broad \oiii component can be applied to explain the different flux ratios of \oiiihb between Type-1 and Type-2 AGN in the BPT diagram.

Magnetars are highly magnetized neutron stars that are characterized by recurrent emission of short-duration bursts in soft gamma-rays/hard X-rays. Recently, FRB 200428 were found to be associated with an X-ray burst from a Galactic magnetar. Two fast radio bursts (FRBs) show mysterious periodic activity. However, whether magnetar X-ray bursts are periodic phenomena is unclear. In this paper, we investigate the period of SGR 1806-20 activity. More than 3000 short bursts observed by different telescopes are collected, including the observation of RXTE, HETE-2, ICE and Konus. We consider the observation windows and divide the data into two sub-samples to alleviate the effect of unevenly sample. The epoch folding and Lomb-Scargle methods are used to derive the period of short bursts. We find a possible period about $ 398.20 \pm 25.45 $ days. While other peaks exist in the periodograms. If the period is real, the connection between short bursts of magnetars and FRBs should be extensively investigated.

We observe 1.3~mm spectral lines at 2000~AU resolution toward four massive molecular clouds in the Central Molecular Zone of the Galaxy to investigate their star formation activities. We focus on several potential shock tracers that are usually abundant in protostellar outflows, including SiO, SO, CH$_3$OH, H$_2$CO, HC$_3$N, and HNCO. We identify 43 protostellar outflows, including 37 highly likely ones and 6 candidates. The outflows are found toward both known high-mass star forming cores and less massive, seemingly quiescent cores, while 791 out of the 834 cores identified based on the continuum do not have detected outflows. The outflow masses range from less than 1~$M_\odot$ to a few tens of $M_\odot$, with typical uncertainties of a factor of 70. We do not find evidence of disagreement between relative molecular abundances in these outflows and in nearby analogs such as the well-studied L1157 and NGC7538S outflows. The results suggest that i) protostellar accretion disks driving outflows ubiquitously exist in the CMZ environment, ii) the large fraction of candidate starless cores is expected if these clouds are at very early evolutionary phases, with a caveat on the potential incompleteness of the outflows, iii) high-mass and low-mass star formation is ongoing simultaneously in these clouds, and iv) current data do not show evidence of difference between the shock chemistry in the outflows that determines the molecular abundances in the CMZ environment and in nearby clouds.

NASA's Transiting Exoplanet Survey Satellite (TESS) mission is expected to discover hundreds of planets via single transits first identified in their light curves. Determining the orbital period of these single transit candidates typically requires a significant amount of follow-up work to observe a second transit or measure a radial velocity orbit. In Yao et al. (2019), we developed simulations that demonstrated the ability to use archival photometric data in combination with TESS to "precover" the orbital period for these candidates with a precision of several minutes, assuming circular orbits. In this work, we incorporate updated models for TESS single transits, allowing for eccentric orbits, along with an updated methodology to improve the reliability of the results. Additionally, we explore how radial velocity (RV) observations can be used to follow up single transit events, using strategies distinct from those employed when the orbital period is known. We find that the use of an estimated period based on a circular orbit to schedule reconnaissance RV observations can efficiently distinguish eclipsing binaries from planets. For candidates that pass reconnaissance RV observations, we simulate RV monitoring campaigns that enable one to obtain an approximate orbital solution. We find this method can regularly determine the orbital periods for planets more massive than 0.5 M_J with orbital periods as long as 100 days.

We constrain the cosmological parameters from cosmic microwave background (CMB) and gravitational wave observations. We not only combine LIGO observations with CMB to measure cosmological parameters, but also forecast the potential abilities of the LISA detector and PTA projects. In the $\Lambda$CDM+$r$ model, the constraint from BK15 and FAST projects have a significant impact on tensor-to-scalar ratio.

We analyse 55 ks of Chandra X-ray observations of the Galactic globular cluster M13. Using the latest radio timing positions of six known millisecond pulsars (MSPs) in M13 from Wang et al. (2020), we detect confident X-ray counterparts to five of the six MSPs at X-ray luminosities of $L_X$(0.3-8 keV)$\sim 3 \times 10^{30} - 10^{31}~{\rm erg~s^{-1}}$, including the newly discovered PSR J1641+3627F. There are limited X-ray counts at the position of PSR J1641+3627A, for which we obtain an upper limit $L_X<1.3 \times 10^{30}~{\rm erg~s^{-1}}$. We analyse X-ray spectra of all six MSPs, which are well-described by either a single blackbody or a single power-law model. We also incorporate optical/UV imaging observations from the Hubble Space Telescope (HST) and find optical counterparts to PSR J1641+3627D and J1641+3627F. Our colour-magnitude diagrams indicate the latter contains a white dwarf, consistent with the properties suggested by radio timing observations. The counterpart to J1641+3627D is only visible in the V band; however, we argue that the companion to J1641+3627D is also a white dwarf, since we see a blackbody-like X-ray spectrum, while MSPs with nondegenerate companions generally show non-thermal X-rays from shocks between the pulsar and companion winds. Our work increases the sample of known X-ray and optical counterparts of MSPs in globular clusters.

Fast radio bursts (FRBs) are extremely powerful sources of radio waves observed at cosmological distances. We use a sophisticated model of FRB observations -- presented in detail in a companion paper -- to fit FRB population parameters using large samples of FRBs detected by ASKAP and Parkes, including seven sources with confirmed host galaxies. Our fitted parameters demonstrate that the FRB population evolves with redshift in a manner consistent with the star-formation rate, ruling out a non-evolving population at 98\% C.L., and 95\% even when considering alternative models. Our estimated maximum FRB energy is $\log_{10} E_{\rm max} [{\rm erg}] = 41.85_{-0.20}^{+0.44}$ (68\% C.L.) assuming a 1\,GHz emission bandwidth, with slope of the cumulative luminosity distribution $\gamma=-1.16_{-0.12}^{+0.11}$. We find a log-mean host DM contribution of $138_{-51}^{+71}$\,pc\,cm$^{-3}$ on top of a typical local (ISM and halo) contribution of $\sim80$\,pc\,cm$^{-3}$, which is higher than most literature values. These results are consistent with the model of FRBs arising as the high-energy limit of magnetar bursts.

We develop a sophisticated model of FRB observations, accounting for the intrinsic cosmological gas distribution and host galaxy contributions, and give a full account of observational biases due to burst width, dispersion measure, and the exact telescope beamshape. Our results offer a significant increase in both accuracy and precision beyond those previously obtained. Using results from ASKAP and Parkes, we present our best-fit FRB population parameters in a companion paper. Here, we consider in detail the expected and fitted distributions in redshift, dispersion measure, and signal-to-noise. We estimate that the unlocalised ASKAP FRBs arise from $z<0.5$, with between a third and a half within $z<0.1$. Our predicted source-counts ("logN--logS") distribution confirms previous indications of a steepening index near the Parkes detection threshold of $1$\,Jy\,ms. We find no evidence for a minimum FRB energy, and rule out $E_{\rm min} > 10^{39}$\,erg at 90\% C.L. Importantly, we find that above a certain DM, observational biases cause the Macquart (DM--z) relation to become inverted, implying that the highest-DM events detected in the unlocalised Parkes and ASKAP samples are unlikely to be the most distant. We do not expect our quantitative estimates in this region to be accurate until it is directly probed with localised FRBs. Since the cause of this effect is a well-understood observational bias however, it is guaranteed to be present to some degree. Works assuming a 1--1 DM--z relation may therefore derive erroneous results.

We solve the two-dimensional hydrodynamic equations of hot accretion flow in the presence of the thermal conduction. The flow is assumed to be in steady-state and axisymmetric, and self-similar approximation is adopted in the radial direction. In this hydrodynamic study, we consider the viscous stress tensor to mimic the effects of the magnetorotational instability for driving angular momentum. We impose the physical boundary conditions at both the rotation axis and the equatorial plane and obtain the solutions in the full $ r-\theta $ space. We have found that thermal conduction is indispensable term for investigating the inflow-wind structure of the hot accretion flows with very low mass accretion rates. One of the most interesting results here is that the disc is convectively stable in hot accretion mode and in the presence of the thermal conduction. Furthermore, the properties of wind and also its driving mechanisms are studied. Our analytical results are consistent with previous numerical simulations of hot accretion flow.

We present our photometric and spectroscopic observations on the peculiar transient AT2018cow. The multi-band photometry covers from peak to $\sim$70 days and the spectroscopy ranges from 5 to $\sim$50 days. The rapid rise ($t_{\mathrm{r}}$$\lesssim$2.9 days), high luminosity ($M_{V,\mathrm{peak}}\sim-$20.8 mag) and fast decline after peak make AT2018cow stand out of any other optical transients. While we find that its light curves show high resemblance to those of type Ibn supernovae. Moreover, the spectral energy distribution remains high temperature of $\sim$14,000 K after $\sim$15 days since discovery. The spectra are featureless in the first 10 days, while some broad emission lines due to H, He, C and O emerge later, with velocity declining from $\sim$14,000 km s$^{-1}$ to $\sim$3000 km s$^{-1}$ at the end of our observations. Narrow and weak He I emission lines emerge in the spectra at $t>$20 days since discovery. These emission lines are reminiscent of the features seen in interacting supernovae like type Ibn and IIn subclasses. We fit the bolometric light curves with a model of circumstellar interaction (CSI) and radioactive decay (RD) of \Ni and find a good fit with ejecta mass $M_{\mathrm{ej}}\sim$3.16 M$_{\odot}$, circumstellar material mass $M_{\mathrm{CSM}}\sim$0.04 M$_{\odot}$, and ejected \Ni mass $M_{^{56}\mathrm{Ni}}\sim$0.23 M$_{\odot}$. The CSM shell might be formed in an eruptive mass ejection of the progenitor star. Furthermore, host environment of AT2018cow implies connection of AT2018cow with massive stars. Combining observational properties and the light curve fitting results, we conclude that AT2018cow might be a peculiar interacting supernova originated from a massive star.

Recently, an indicative evidence of a stochastic process, reported by the NANOGrav Collaboration based on the analysis of 12.5-year pulsar timing array data which might be interpreted as a potential stochastic gravitational wave signal, has aroused keen interest of theorists. The first-order cosmological confinement/deconfinement phase transition at the QCD scale could be one of the cosmological sources for the NANOGrav signal. If the phase transition is flavor dependent and happens sequentially, it is important to find the dominant phase transition producing the gravitational waves that match the NANOGrav signal during the evolution of the universe. In this paper, we would like to illustrate that the NANOGrav signal could be generated from confinement/deconfinement transition in either heavy static quarks with a zero baryon chemical potential, or quarks with a finite baryon chemical potential, or a pure gluon system. Moreover, the quark confinement as compared with the gluon confinement is more likely to be the source for the NANOGrav signal regardless of whether the chemical potential is finite or not. Future observation will help to distinguish between different scenarios.

Among exoplanets, the small-size population constitutes the dominant one, with a diversity of properties and compositions ranging from rocky to gas dominated envelope. While a large fraction of them have masses and radii similar to or smaller than Neptune, yet none share common properties in term of orbital period and insulation with our ice giants. These exoplanets belong to multi-planet systems where planets are closely packed within the first tenth of AU and often exposed to strong irradiation from their host star. Their formation process, subsequent evolution, and fate are still debated and trigger new developments of planet formation models. This paper reviews the characteristics and properties of this extended sample of planets with radii between $\sim$ 1.6 and 4.0$_\oplus$. Even though we still lack real Neptune/Uranus analogues, these exoplanets provide us with key observational constraints that allow the formation of our ice giants to be placed in a more general framework than the sole example of our solar system.

Since the discovery of exoplanets around pulsars, there has been a debate on their origin. Popular scenarios include in situ formation or the dynamical capture of a planet in a dense stellar system. The possibility of a planet surviving its host star's supernova is often neglected, because a planet in orbit around a single exploding star is not expected to survive the supernova. A circum-binary planet, however, may stand a chance in staying bound when one of the binary components explodes. We investigate the latter and constrain the distribution of post-supernova orbital parameters of circum-binary planets. This is done by performing population synthesis calculations of binary stars until the first supernova. Just before the supernova, we add a planet in orbit around the binary to study its survivability. In our supernova model, the exploding star's mass is assumed to change instantaneously, and we apply a velocity kick to the newly formed remnant. The mass loss and velocity kick affect the orbits of the two stars and the planet. Only $2 \cdot 10^{-3}$ of systems survive the supernova while keeping the circum-binary planet bound. The surviving planetary orbits are wide ($a \apgt 10$ au) and eccentric ($e \apgt 0.3$). It turns out much more likely ($3\cdot 10^{-2}$ system fraction) that the newly formed compact object is ejected from the system, leaving the planet bound to its companion star in a highly eccentric orbit (typically $\apgt 0.9$). We expect that the Milky way Galaxy hosts at most $10$ x-ray binaries that are still orbited by a planet, and $\aplt 150$ planets that survived in orbit around the compact object's companion. These numbers should be convolved with the fraction of massive binaries that is orbited by a planet.

We present the long-term X-ray spectral and temporal analysis of a 'bare-type AGN' Ark 120. We consider the observations from XMM-Newton, Suzaku, Swift, and NuSTAR from 2003 to 2018. The spectral properties of this source are studied using various phenomenological and physical models present in the literature. We report (a) the variations of several physical parameters, such as the temperature and optical depth of the electron cloud, the size of the Compton cloud, and accretion rate for the last fifteen years. The spectral variations are explained from the change in the accretion dynamics; (b) the X-ray time delay between 0.2-2 keV and 3-10 keV light-curves exhibited zero-delay in 2003, positive delay of 4.71 \pm 2.1 ks in 2013, and negative delay of 4.15 \pm 1.5 ks in 2014. The delays are explained considering Comptonization, reflection, and light-crossing time; (c) the long term intrinsic luminosities obtained using nthcomp, of the soft-excess and the primary continuum show a correlation with a Pearson Correlation Coefficient of 0.922. This indicates that the soft-excess and the primary continuum are originated from the same physical process. From a physical model fitting, we infer that the soft excess for Ark 120 could be due to a small number of scatterings in the Compton cloud. Using Monte-Carlo simulations, we show that indeed the spectra corresponding to fewer scatterings could provide a steeper soft-excess power-law in the 0.2-3 keV range. Simulated luminosities are found to be in agreement with the observed values.

Galaxies in dense environments are subject to interactions and mechanisms which directly affect their evolution by lowering their gas fractions and reducing their star-forming capacity earlier than their isolated counterparts. The aim of our project is to get new insights about the role of environment on the stellar and baryonic content of galaxies using a kinematic approach, through the study of the Tully-Fisher relation (TFR). We study a sample of galaxies in 8 groups spanning a redshift range of $0.5<z<0.8$ and located in 10 pointings of the MAGIC MUSE Guaranteed Time Observations program. We perform a morpho-kinematics analysis of this sample and set up a selection based on galaxy size, [OII] emission line doublet signal-to-noise ratio, bulge-to-disk ratio and nuclear activity to construct a robust kinematic sample of 67 star-forming galaxies. This selection considerably reduces the number of outliers in the TFR, which are predominantly dispersion-dominated galaxies. Our results suggest a significant offset of the TFR zero-point between galaxies in low- and high-density environments, whatever kinematics estimator is used. This can be interpreted as a decrease of either stellar mass by $\sim 0.05 - 0.3$ dex or an increase of rotation velocity by $\sim 0.02 - 0.06$ dex for galaxies in groups, depending on the samples used for comparison. We also studied the stellar and baryon mass fractions within stellar disks and found they both increase with stellar mass, the trend being more pronounced for the stellar component alone. These fractions do not exceed 50%. We show that this evolution of the TFR is consistent either with a decrease of star formation or with a contraction of the mass distribution due to the environment. These two effects probably act together with their relative contribution depending on the mass regime.

The celebrated CLEAN algorithm has been the cornerstone of deconvolution algorithms in radio interferometry almost since its conception in the 1970s. For all its faults, CLEAN is remarkably fast, robust to calibration artefacts and in its ability to model point sources. We demonstrate how the same assumptions that afford CLEAN its speed can be used to accelerate more sophisticated deconvolution algorithms.

We study the structure of galactic halos within a scalar dark matter model, endowed with a repulsive quartic self-interaction, capable of undergoing the superfluid phase transition in high-density regions. We demonstrate that the thermalized cores are prone to fragmentation into superfluid droplets due to the Jeans instability. Furthermore, since cores of astrophysical size may be generated only when most of the particles comprising the halo reside in a highly degenerate phase-space, the well-known bound on the dark matter self-interaction cross section inferred from the collision of clusters needs to be revised, accounting for the enhancement of the interaction rate due to degeneracy. As a result, generation of kpc-size superfluid solitons, within the parameter subspace consistent with the Bullet Cluster bound, requires dark matter particles to be ultra-light.

Evidence is reported for an 8 MeV neutrino line associated with SN 1987A. This discovery is based on an analysis of 997 events recorded in the Kamiokande-II detector on the day of the supernova, as well as other data, and it has very high statistical significance.

Using five year monitoring observations, we did a blind search for pulses for rotating radio transient (RRAT) J0139+33 and PSR B0320+39. At the interval \pm 1.5m of the time corresponding to the source passing through the meridian, we detected 39377 individual pulses for the pulsar B0320+39 and 1013 pulses for RRAT J0139+33. The share of registered pulses from the total number of observed periods for the pulsar B0320+39 is 74%, and for the transient J0139+33 it is 0.42%. Signal-to-noise ratio (S/N) for the strongest registered pulses is approximately equal to: S/N = 262 (for B0320+39) and S/N = 154 (for J0139+33). Distributions of the number of detected pulses in S/N units for the pulsar and for the rotating transient are obtained. The distributions could be approximated with a lognormal and power dependencies. For B0320+39 pulsar, the dependence is lognormal, it turns into a power dependence at high values of S/N, and for RRAT J0139+33, the distribution of pulses by energy is described by a broken (bimodal) power dependence with an exponent of about 0.4 and 1.8 (S/N < 19 and S/N > 19). We have not detected regular (pulsar) emission of J0139+33. Analysis of the obtained data suggests that RRAT J0139+33 is a pulsar with giant pulses.

The design and implementation of astronomical cameras based on the large-format CCD and CMOS detectors is described in this paper. The Dinacon-5 controller is used for work with the CCDs and to achieve high performance and low noise. A new controller is designed for CMOS sensors. The main characteristics of the provided systems are estimated on the basis of experimental data. The spatial autocorrelation analysis is applied for PSF estimation. The obtained test results are presented.

We present a comprehensive analysis of the formation and evolution of a fan-spine-like configuration that developed over a complex photospheric configuration where dispersed negative polarity regions were surrounded by positive polarity regions. This unique photospheric configuration, analogous to the geological "atoll" shape, hosted four homologous flares within its boundary. Computation of the degree of squashing factor (Q) maps clearly revealed an elongated region of high Q-values between the inner and outer spine-like lines, implying the presence of an hyperbolic flux tube (HFT). The coronal region associated with the photospheric atoll configuration was distinctly identified in the form of a diffused dome-shaped bright structure directly observed in EUV images. A filament channel resided near the boundary of the atoll region. The activation and eruption of flux ropes from the filament channel led to the onset of four eruptive homologous quasi-circular ribbon flares within an interval of $\approx$11 hours. During the interval of the four flares, we observed continuous decay and cancellation of negative polarity flux within the atoll region. Accordingly, the apparent length of the HFT gradually reduced to a null-point-like configuration before the fourth flare. Prior to each flare, we observed localised brightening beneath the filaments which, together with flux cancellation, provided support for the tether-cutting model of solar eruption. The analysis of magnetic decay index revealed favourable conditions for the eruption, once the pre-activated flux ropes attained the critical heights for torus instability.

Context. Planetary mass and radius data are showing a wide variety in densities of low-mass exoplanets. This includes sub-Neptunes, whose low densities can be explained with the presence of a volatile-rich layer. Water is one of the most abundant volatiles, which can be in the form of different phases depending on the planetary surface conditions. To constrain their composition and interior structure, it is required to develop models that calculate accurately the properties of water at its different phases. Aims. We present an interior structure model that includes a multiphase water layer with steam, supercritical and condensed phases. We derive the constraints for planetary compositional parameters and their uncertainties, focusing on the multiplanetary system TRAPPIST-1, which presents both warm and temperate planets. Methods. We use a 1D steam atmosphere in radiative-convective equilibrium with an interior whose water layer is in supercritical phase self-consistently. For temperate surface conditions, we implement liquid and ice Ih to ice VII phases in the hydrosphere. We adopt a MCMC inversion scheme to derive the probability distributions of core and water compositional parameters Results. We refine the composition of all planets and derive atmospheric parameters for planets b and c. The latter would be in a post-runaway greenhouse state and could be extended enough to be probed by space mission such as JWST. Planets d to h present condensed ice phases, with maximum water mass fractions below 20%. Conclusions. The derived amounts of water for TRAPPIST-1 planets show a general increase with semi-major axis, with the exception of planet d. This deviation from the trend could be due to formation mechanisms, such as migration and an enrichment of water in the region where planet d formed, or an extended CO$_{2}$-rich atmosphere.

Relic gravitational waves (GWs) can be produced by primordial magnetic fields. However, not much is known about the resulting GW amplitudes and their dependence on the details of the generation mechanism. Here we treat magnetic field generation through the chiral magnetic effect (CME) as a generic mechanism and explore its dependence on the speed of generation (the product of magnetic diffusivity and characteristic wavenumber) and the speed characterizing the limiting magnetic field strength expected from the CME. When the latter exceeds the former (regime I), the regime applicable to the early universe, we obtain an inverse cascade with moderate GW energy that scales with the third power of the magnetic energy. When the generation speed exceeds the CME limit (regime II), the GW energy continues to increase without a corresponding increase of magnetic energy. In the early kinematic phase, the GW energy spectrum (per linear wavenumber interval) has opposite slopes in both regimes and is characterized by an inertial range spectrum in regime I and a white noise spectrum in regime II. The occurrence of these two slopes is shown to be a generic consequence of a nearly monochromatic exponential growth of the magnetic field. The resulting GW energy is found to be proportional to the fifth power of the limiting CME speed and the first (third) power of the growth speed regime I (II).

We use simple models of the spatial structure of the quasar broad line region (BLR) to investigate the properties of so-called ghostly damped Lyman-{\alpha} (DLA) systems detected in SDSS data. These absorbers are characterized by the presence of strong metal lines but no Hi Lyman-{\alpha} trough is seen in the quasar spectrum indicating that, although the region emitting the quasar continuum is covered by an absorbing cloud, the BLR is only partially covered. One of the models has a spherical geometry, another one is the combination of two wind flows whereas the third model is a Keplerian disk. The models can reproduce the typical shape of the quasar Lyman-{\alpha} emission and different ghostly configurations. We show that the DLA Hi column density can be recovered precisely independently of the BLR model used. The size of the absorbing cloud and its distance to the centre of the AGN are correlated. However it may be possible to disentangle the two using an independent estimate of the radius from the determination of the particle density. Comparison of the model outputs with SDSS data shows that the wind and disk models are more versatile than the spherical one and can be more easily adapted to the observations. For all the systems we derive logN(Hi)(cm^{-2})>20.5. With higher quality data it may be possible to distinguish between the models.

An algorithm has been developed for finding the global minimum of a multidimensional error function by fitting model spectral maps into observed ones. Principal component analysis is applied to reduce the dimensionality of the model and the coupling degree between the parameters, and to determine the region of the minimum. The k-nearest neighbors method is used to calculate the optimal parameter values. The algorithm is used to estimate the physical parameters of the contracting dense star-forming core of L1287. Maps in the HCO+(1-0), H13CO+(1-0), HCN(1-0), and H13CN(1-0) lines, calculated within a 1D microturbulent model, are fitted into the observed ones. Estimates are obtained for the physical parameters of the core, including the radial profiles of density ($\propto r^{-1.7}$), turbulent velocity ($\propto r^{-0.4}$), and contraction velocity ($\propto r^{-0.1}$). Confidence intervals are calculated for the parameter values. The power-law index of the contraction-velocity radial profile, considering the determination error, is lower in absolute terms than the expected one in the case of gas collapse onto the protostar in free fall. This result can serve as an argument in favor of a global contraction model for the L1287 core.

OCTO-TIGER is an astrophysics code to simulate the evolution of self-gravitating and rotat-ing systems of arbitrary geometry based on the fast multipole method, using adaptive mesh refinement. OCTO-TIGER is currently optimised to simulate the merger of well-resolved stars that can be approximated by barotropic structures, such as white dwarfs or main sequence stars. The gravity solver conserves angular momentum to machine precision, thanks to a correction algorithm. This code uses HPX parallelization, allowing the overlap of work and communication and leading to excellent scaling properties, allowing for the computation of large problems in reasonable wall-clock times. In this paper, we investigate the code performance and precision by running benchmarking tests. These include simple problems, such as the Sod shock tube, as well as sophisticated, full, white-dwarf binary simulations. Results are compared to analytic solutions, when known, and to other grid based codes such as FLASH. We also compute the interaction between two white dwarfs from the early mass transfer through to the merger and compare with past simulations of similar systems. We measure OCTO-TIGERs scaling properties up to a core count of 80,000, showing excellent performance for large problems. Finally, we outline the current and planned areas of development aimed at tackling a number of physical phenomena connected to observations of transients.

Operating experience from fusion research shows how Spitzer resistivity may render ohmic heating in the chromosphere self limiting and thus serve to define the lower margin of the transition region. Its upper margin is at about 6000 K, where radiative cooling of He:H plasma decelerates sharply. The third and last stage in the proposed scheme is expansion into the tenuous plasma of space, which leads to the acceleration of ions to high energies, long recorded by spacecraft instruments. There is thus dynamic continuity all the way from the solar interior, the energy source for spinning columns in the Rayleigh Benard setting of the convection zone, to the coronal exhalation of the solar wind, a finding which should benefit the analysis of space weather, witness the association between helium in the solar wind and the incidence of coronal mass ejections.

We present the first study on the amplification of magnetic fields by the turbulent dynamo in the highly subsonic regime, with Mach numbers ranging from $10^{-3}$ to $0.4$. We find that for the lower Mach numbers the saturation efficiency of the dynamo, $(E_{\mathrm{mag}}/E_{\mathrm{kin}})_{\mathrm{sat}}$, increases as the Mach number decreases. Even in the case when injection of energy is purely through longitudinal forcing modes, $(E_{\mathrm{mag}}/E_{\mathrm{kin}})_{\mathrm{sat}}$ $\gtrsim 10^{-2}$ at a Mach number of $10^{-3}$. We apply our results to magnetic field amplification in the early Universe and predict that a turbulent dynamo can amplify primordial magnetic fields to $\gtrsim$ $10^{-16}$ Gauss on scales up to 0.1 pc and $\gtrsim$ $10^{-13}$ Gauss on scales up to 100 pc. This produces fields compatible with lower limits of the intergalactic magnetic field inferred from blazar $\gamma$-ray observations.

An one-parameter regularization freedom of the Hamiltonian constraint for loop quantum gravity is analyzed. The corresponding spatially flat, homogenous and isotropic model includes the two well-known models of loop quantum cosmology as special cases. The quantum bounce nature is tenable in the generalized cases. For positive value of the regularization parameter, the effective Hamiltonian leads to an asymptotic de-Sitter branch of the Universe connecting to the standard Friedmann branch by the quantum bounce. Remarkably, by suitably choosing the value of the regularization parameter, the observational cosmological constant can emerge at large volume limit from the effect of quantum gravity, and the effective Newtonian constant satisfies the experimental restrictions in the meantime.

Nuclear clusters or voids in the inner crust of neutron stars were predicted to have various shapes collectively nicknamed nuclear pasta. The recent review in Ref. \cite{Lopez1} by L\'opez, Dorso and Frank summarized their systematic investigations into properties especially the morphological and thermodynamical phase transitions of the nuclear pasta within a Classical Molecular Dynamics model, providing further stimuli to find more observational evidences of the predicted nuclear pasta in neutron stars.

In light of the growing interest in the Hayward black hole solution, a detailed study on the corresponding lapse function and its roots is presented. The lapse function is expressed in terms of the classical Schwarzschild radius $r_s$ and the Hayward's parameter $l$. Both of these quantities are used as thermodynamic variables to find related thermodynamic quantities. In this context, the variable $l$ is associated with a canonical conjugate variable $\mathcal{F}_H$, and a free energy $\Xi$. Moreover, a second order phase transition is found to appears at $l\approx0.333\,r_s$.

A general covariant local field theory of the holographic dark energy model is presented. It turns out the low energy effective theory of the holographic dark energy is the massive gravity theory whose graviton has 3 polarisations, including one scalar mode and two tensor modes. The Compton wavelength is the size of the future event horizon of the universe. The physical interpretation for the UV-IR correspondence $ \Lambda\sim\sqrt{M_p/L}$ in the holographic dark energy model is provided in the framework of our effective field theory, where $L$ is interpreted as the graviton's Compton wavelength, and $\Lambda$ is interpreted as the energy scale where the scalar graviton strongly couples to itself.

We discuss deep learning inference for the neutron star equation of state (EoS) using the real observational data of the mass and the radius. We make a quantitative comparison between the conventional polynomial regression and the neural network approach for the EoS parametrization. For our deep learning method to incorporate uncertainties in observation, we augment the training data with noise fluctuations corresponding to observational uncertainties. Deduced EoSs can accommodate a weak first-order phase transition, and we make a histogram for likely first-order regions. We also find that our observational data augmentation has a byproduct to tame the overfitting behavior. To check the performance improved by the data augmentation, we set up a toy model as the simplest inference problem to recover a double-peaked function and monitor the validation loss. We conclude that the data augmentation could be a useful technique to evade the overfitting without tuning the neural network architecture such as inserting the dropout.