Abstract visibility
Abstract text size

Papers for Friday, Jun 25 2021

Papers with local authors

Peter Harrington, Mustafa Mustafa, Max Dornfest, Benjamin Horowitz, Zarija Lukić

10 pages, 6 figures

3 votes
Paper 10 — arXiv:2106.12662
3 votes
Paper 10 — arXiv:2106.12662

Full-physics cosmological simulations are powerful tools for studying the formation and evolution of structure in the universe but require extreme computational resources. Here, we train a convolutional neural network to use a cheaper N-body-only simulation to reconstruct the baryon hydrodynamic variables (density, temperature, and velocity) on scales relevant to the Lyman-$\alpha$ (Ly$\alpha$) forest, using data from Nyx simulations. We show that our method enables rapid estimation of these fields at a resolution of $\sim$20kpc, and captures the statistics of the Ly$\alpha$ forest with much greater accuracy than existing approximations. Because our model is fully-convolutional, we can train on smaller simulation boxes and deploy on much larger ones, enabling substantial computational savings. Furthermore, as our method produces an approximation for the hydrodynamic fields instead of Ly$\alpha$ flux directly, it is not limited to a particular choice of ionizing background or mean transmitted flux.

Benjamin Horowitz, Max Dornfest, Zarija Lukić, Peter Harrington

13 pages, 11 figures

0 votes
Paper 12 — arXiv:2106.12675
0 votes
Paper 12 — arXiv:2106.12675

Generating large volume hydrodynamical simulations for cosmological observables is a computationally demanding task necessary for next generation observations. In this work, we construct a novel fully convolutional variational auto-encoder (VAE) to synthesize hydrodynamic fields conditioned on dark matter fields from N-body simulations. After training the model on a single hydrodynamical simulation, we are able to probabilistically map new dark matter only simulations to corresponding full hydrodynamical outputs. By sampling over the latent space of our VAE, we can generate posterior samples and study the variance of the mapping. We find that our reconstructed field provides an accurate representation of the target hydrodynamical fields as well as a reasonable variance estimates. This approach has promise for the rapid generation of mocks as well as for implementation in a full Bayesian inverse model of observed data.

Haopeng Zhang, Yuqin Chen, Gang Zhao

19 pages, 11 figures, accepted for publication in the Astrophysical Journal

0 votes
Paper 25 — arXiv:2106.12841
0 votes
Paper 25 — arXiv:2106.12841

Radial migration is an important process in the evolution of the Galactic disk. The metallicity gradient of open clusters and its outliers provide an effective way to probe for this process. In this work, we compile metallicity, age, and kinematic parameters for 225 open clusters and carry out a quantitative analysis of radial migration via the calculated migration distances. Based on clusters with age $< 0.5$ Gyr, we obtain the present-day metallicity gradient of $-0.074 \pm 0.007$ dex/kpc. Along this gradient distributes three sequences, and clusters in the upper, the middle, and the lower groups are found to be old outward-migrators, in-situ clusters, and inward-migrators, respectively. The migration distance increases with age, but its most effective time is probably less than 3 Gyr. The metallicity gradient breaks out at $R_g$ (guiding center radius) $\sim11.5$ kpc, which is caused by the lack of young open clusters in the outer disk and the presence of old outward-migrators in the upper sequence. It shows that this boundary is related to the different effects of radial migration between the inner and outer disks. We also found many special open clusters in and near the outer disk of $R > 11$ kpc and a steeper metallicity gradient from the inner disk of $R_g < 7$ kpc, which tells a complicated evolution history of the Galactic disk by different effects of stellar radial migration.

All other papers

We present cosmological parameter measurements from the effective field theory-based full-shape analysis of the power spectrum of emission line galaxies (ELGs). First, we perform extensive tests on simulations and determine appropriate scale cuts for the perturbative description of the ELG power spectrum. We study in detail non-linear redshift-space distortions ('fingers-of-God') for this sample and show that they are somewhat weaker than those of luminous red galaxies. This difference is not significant for current data, but may become important for future surveys like Euclid/DESI. Then we analyze recent measurements of the ELG power spectrum from the extended Baryon acoustic Oscillation Spectroscopic Survey (eBOSS) within the $\nu\Lambda$CDM model. Combined with the BBN baryon density prior, the ELG power spectrum alone constrains the matter density $\Omega_m=0.287_{-0.061}^{+0.04}$ and the mass fluctuation amplitude $\sigma_8=0.73_{-0.13}^{+0.09}$. Combining with other full-shape and BAO data we measure $\Omega_m=0.320_{-0.016}^{+0.014}$, $\sigma_8=0.700_{-0.044}^{+0.04}$, and the Hubble constant $H_0=69.14_{-1.1}^{+1.0}$ km/s/Mpc. The total neutrino mass is constrained to be $M_{\rm tot}<0.61$ eV (95% CL) from the BBN, full-shape and BAO data only.

A. Veropalumbo, I. Saez Casares, E. Branchini, B. Granett, L. Guzzo, F. Marulli, M. Moresco, L. Moscardini, A. Pezzotta, S. de la Torre

18 pages, 12 figures, 4 tables. Submitted to MNRAS. Comments welcome

We measure the galaxy 2- and 3-point correlation functions at $z=[0.5,0.7]$ and $z=[0.7, 0.9]$, from the final data release of the VIPERS survey (PDR2). We model the two statistics including a nonlinear 1-loop model for the 2-point function and a tree-level model for the 3-point function, and perform a joint likelihood analysis. The entire process and nonlinear corrections are tested and validated through the use of the 153 highly realistic VIPERS mock catalogues, showing that they are robust down to scales as small as 10 $h^{-1} \, \mathrm{Mpc}$. The mocks are also adopted to compute the covariance matrix that we use for the joint 2- and 3-point analysis. Despite the limited statistics of the two (volume-limited) sub-samples analysed, we demonstrate that such a combination successfully breaks the degeneracy existing at 2-point level between clustering amplitude $\sigma_8$, linear bias $b_1$ and the linear growth rate of fluctuations $f$. For the latter, in particular, we measure $f(z=0.61)=0.64^{+0.55}_{-0.37}$ and $f(z=0.8)=1.0\pm1.0$, while the amplitude of clustering is found to be $\sigma_8(z=0.61)=0.50\pm 0.12$ and $\sigma_8(z=0.8)=0.39^{+0.11}_{-0.13}$. These values are in excellent agreement with the extrapolation of a Planck cosmology.

Seungsoo Hong, Dongwook Lim, Chul Chung, Jaeyeon Kim, Sang-Il Han, Young-Wook Lee

Accepted for publication in AJ

We performed low-resolution spectroscopy for the red giant branch stars in an intriguing globular cluster (GC) NGC 2808, which hosts subpopulations with extreme helium and light-element abundances. In order to trace N, C, and Ca abundance differences among subpopulations, we measured CN, CH, and Ca II H&K spectral indices, respectively. We identified four subpopulations (G1, G2, G3, and G4) from CN and CH strength, with CN-weak/CH-strong G1, CN-intermediate/CH-strong G2, CN-strong/CH-intermediate G3, and CN-strong/CH-weak G4. Compared to [Na/O] from high-resolution spectroscopy, we show that CN index can more clearly separate G1 and G2. Since CN traces N abundance in a GC, it implies that G1 and G2 would show a larger difference in [N/Fe] compared to [Na/Fe], as predicted by chemical evolution models. Later generation stars G3 and G4, however, are better separated with high-resolution spectroscopy. We also found that G4 shows a stronger Ca II H&K line strength compared to that of G1, but we suspect this to be a result of unusually strong He enhancement and/or Mg depletion in G4 of this GC. This work illustrates that combining low- and high-resolution spectroscopic studies can improve the separation of subpopulations in GCs.

Alec J. M. Thomson, T. L. Landecker, N. M. McClure-Griffiths, John M. Dickey, J. L. Campbell, Ettore Carretti, S. E. Clark, Christoph Federrath, B. M. Gaensler, J. L. Han, Marijke Haverkorn, Alex. S. Hill, S.A. Mao, Anna Ordog, Luke Pratley, Wolfgang Reich, Cameron L. Van Eck, J. L. West, M. Wolleben

25 pages, 14 figures. Accepted for publication in MNRAS

Using the Global Magneto-Ionic Medium Survey (GMIMS) Low-Band South (LBS) southern sky polarization survey, covering 300 to 480 MHz at 81 arcmin resolution, we reveal the brightest region in the Southern polarized sky at these frequencies. The region, G150-50, covers nearly 20deg$^2$, near (l,b)~(150 deg,-50 deg). Using GMIMS-LBS and complementary data at higher frequencies (~0.6--30 GHz), we apply Faraday tomography and Stokes QU-fitting techniques. We find that the magnetic field associated with G150-50 is both coherent and primarily in the plane of the sky, and indications that the region is associated with Radio Loop II. The Faraday depth spectra across G150-50 are broad and contain a large-scale spatial gradient. We model the magnetic field in the region as an expanding shell, and we can reproduce both the observed Faraday rotation and the synchrotron emission in the GMIMS-LBS band. Using QU-fitting, we find that the Faraday spectra are produced by several Faraday dispersive sources along the line-of-sight. Alternatively, polarization horizon effects that we cannot model are adding complexity to the high-frequency polarized spectra. The magnetic field structure of Loop II dominates a large fraction of the sky, and studies of the large-scale polarized sky will need to account for this object. Studies of G150-50 with high angular resolution could mitigate polarization horizon effects, and clarify the nature of G150-50.

Philipp Girichidis

20 pages, 16 figures, submitted to MNRAS, comments welcome

Magnetic fields are an important component of the interstellar medium (ISM) and exhibit strongly varying field strengths and a non-trivial correlation with the gas density. Its dynamical impact varies between individual regions of the ISM and correlates with the orientation of the field with respect to the gas structures. Using high-resolution magneto-hydrodynamical simulations of the ISM we explore the connection between the orientation of the field and the dynamical state of the gas. We find that the onset of gravitational instability in molecular gas above a density of $\rho\sim10^{-21}\,\mathrm{g\,cm^{-3}}$ $(n\sim400\,\mathrm{cm^{-3}})$ coincides with an alignment of the magnetic field lines and the gas flow. At this transition the gradient of the density changes from mainly perpendicular to preferentially parallel to the field lines. A connection between the three-dimensional alignment and projected two-dimensional observables is non-trivial, because of a large dispersion of the magnetic field orientation along the line of sight. The turbulent correlation lengths can be small compared to the typical integration lengths. As a consequence the small scale signal of the orientation can sensitively depend on the line of sight or the dynamical state of the cloud, can fluctuate stochastically or be completely averaged out. With higher spatial resolution more small scale structures are resolved, which aggravates the link between magneto-hydrodynamical quantities and projected observables.

Siyao Xu, Alex Lazarian

28 pages, 12 figures, published in Reviews of Modern Plasma Physics

Small-scale turbulent dynamo is responsible for the amplification of magnetic fields on scales smaller than the driving scale of turbulence in diverse astrophysical media. Most earlier dynamo theories concern the kinematic regime and small-scale magnetic field amplification. Here we review our recent progress in developing the theories for the nonlinear dynamo and the dynamo regime in a partially ionized plasma. The importance of reconnection diffusion of magnetic fields is identified for both the nonlinear dynamo and magnetic field amplification during gravitational contraction. For the dynamo in a partially ionized plasma, the coupling state between neutrals and ions and the ion-neutral collisional damping can significantly affect the dynamo behavior and the resulting magnetic field structure. We present both our analytical predictions and numerical tests with a two-fluid dynamo simulation on the dynamo features in this regime. In addition, to illustrate the astrophysical implications, we discuss several examples for the applications of the dynamo theory to studying magnetic field evolution in both preshock and postshock regions of supernova remnants, in weakly magnetized molecular clouds, during the (primordial) star formation, and during the first galaxy formation.

J. A. Rencoret, C. Aguilera-Gómez, A. Reisenegger

17 pages, 11 figures, XVI SOCHIAS Annual Meeting, and APS conference 2021: Neutron starquakes caused by spindown

Pulsars show a steady decrease in their rotational frequency, occasionally interrupted by sudden spin-ups called glitches, whose physical origin is still a mystery. One suggested explanation for at least the small glitches are starquakes, i.e., failures of the solid neutron star crust, in which the progressive reduction of the centrifugal force deforms the star, stressing the solid until it breaks. This produces a spin-up, dissipating energy inside the star. We explore this suggestion by analyzing a mostly analytical model in order to understand the possible consequences of starquakes, particularly whether they can explain at least the small glitches. We analyze the deformations and strains produced by the decreasing centrifugal force, modeling the neutron star with a fluid core and a solid crust, each with uniform density and with the core possibly denser than the crust. The deformation of a star with very different densities in the core and crust is qualitatively different from the previously studied case of equal densities. The former more closely resembles the behavior of a fluid star, in which the core-crust interface is a surface of constant gravitational plus centrifugal potential. Regardless of the uncertain breaking strain, the glitch activity in this model is several orders of magnitude smaller than observed, even if only small glitches are considered. For a large breaking strain, suggested by simulations, glitches due to starquakes could be roughly of the correct size, but much less frequent than observed glitches. The energy released in each such glitch is much larger than in the standard model of angular momentum transfer from a faster rotating superfluid in the inner crust. We also confirm that stresses in the crust can in principle support an ellipticity much larger than some observational upper limits from pulsar timing and continuous gravitational wave searches.

Clare Higgs, Alan McConnachie

15 pages, 6 figures, 5 tables. Accepted to MNRAS

We compare and contrast the stellar structures of isolated Local Group dwarf galaxies, as traced by their oldest stellar populations, with the satellite dwarf galaxies of the Milky Way and M31. All Local Group dwarfs with Mv < -6 and surface brightness < 26.5 mags. per square arcsec. are considered, taking advantage of measurements from surveys that use similar observations and analysis techniques. For the isolated dwarfs, we use the results from Solitary Local (Solo) Dwarf Galaxy Survey. We begin by confirming that the structural and dynamical properties of the two satellite populations are not obviously statistically different from each other, but we note that there many more satellites around M31 than around the Milky Way down to equivalent magnitude and surface brightness limits. We find that dwarfs in close proximity to a massive galaxy generally show more scatter in their Kormendy relations than those in isolation. Specifically, isolated Local Group dwarf galaxies show a tighter trend of half-light radius versus magnitude than the satellite populations, and similar effects are also seen for related parameters. There appears to be a transition in the structural and dynamical properties of the dwarf galaxy population around ~400 kpc from the Milky Way and M31, such that the smallest, faintest, most circular dwarf galaxies are found closer than this separation. We discuss the impact of selection effects on our analysis, and we argue that our results point to the significance of tidal interactions on the population of systems within approximately 400 kpc from the MW and M31.

Alexander P. Ji, Sergey E. Koposov, Ting S. Li, Denis Erkal, Andrew B. Pace, Joshua D. Simon, Vasily Belokurov, Lara R. Cullinane, Gary S. Da Costa, Kyler Kuehn, Geraint F. Lewis, Dougal Mackey, Nora Shipp, Jeffrey D. Simpson, Daniel B. Zucker, Terese T. Hansen, Joss Bland-Hawthorn, S5 Collaboration

23 pages (+appendix), 10 figures, submitted to ApJ. Video of Antlia 2 tidal disruption here: this https URL

We present new spectroscopic observations of the diffuse Milky Way satellite galaxies Antlia 2 and Crater 2, taken as part of the Southern Stellar Stream Spectroscopic Survey (S5). The new observations approximately double the number of confirmed member stars in each galaxy and more than double the spatial extent of spectroscopic observations in Antlia 2. A full kinematic analysis, including Gaia EDR3 proper motions, detects a clear velocity gradient in Antlia 2 and a tentative velocity gradient in Crater 2. The velocity gradient magnitudes and directions are consistent with particle stream simulations of tidal disruption. Furthermore, the orbit and kinematics of Antlia 2 require a model that includes the reflex motion of the Milky Way induced by the Large Magellanic Cloud. We also find that Antlia 2's metallicity was previously overestimated, so it lies on the empirical luminosity-metallicity relation and is likely only now experiencing substantial stellar mass loss. This low stellar mass loss contrasts with current dynamical models of Antlia 2's size and velocity dispersion, which require it to have lost more than 90% of its stars to tides. Overall, the new kinematic measurements support a tidal disruption scenario for the origin of these large and extended dwarf spheroidal galaxies.

Jonathan S. Gómez, Nelson D. Padilla, John C. Helly, Cedric G. Lacey, Carlton M. Baugh, Claudia del P. Lagos

20 pages, 16 figures

We examine the effect of using different halo finders and merger tree building algorithms on galaxy properties predicted using the GALFORM semi-analytical model run on a high resolution, large volume dark matter simulation. The halo finders/tree builders HBT, ROCKSTAR, SUBFIND and VELOCIRAPTOR differ in their definitions of halo mass, on whether only spatial or phase-space information is used, and in how they distinguish satellite and main haloes; all of these features have some impact on the model galaxies, even after the trees are post-processed and homogenised by GALFORM. The stellar mass function is insensitive to the halo and merger tree finder adopted. However, we find that the number of central and satellite galaxies in GALFORM does depend slightly on the halo finder/tree builder. The number of galaxies without resolved subhaloes depends strongly on the tree builder, with VELOCIRAPTOR, a phase-space finder, showing the largest population of such galaxies. The distributions of stellar masses, cold and hot gas masses, and star formation rates agree well between different halo finders/tree builders. However, because VELOCIRAPTOR has more early progenitor haloes, with these trees GALFORM produces slightly higher star formation rate densities at high redshift, smaller galaxy sizes, and larger stellar masses for the spheroid component. Since in all cases these differences are small we conclude that, when all of the trees are processed so that the main progenitor mass increases monotonically, the predicted GALFORM galaxy populations are stable and consistent for these four halo finders/tree builders.

Antoine Marchal, Peter G. Martin, Munan Gong

32 pages, 31 figures, 7 tables. Accepted for publication in ApJ

The physical properties of galactic halo gas have a profound impact on the life cycle of galaxies. As gas travels through a galactic halo, it undergoes dynamical interactions, influencing its impact on star formation and the chemical evolution of the galactic disk. In the Milky-Way halo, considerable effort has been made to understand the spatial distribution of neutral gas, which are mostly in the form of large complexes. However, the internal variations of their physical properties remains unclear. In this study, we investigate the thermal and dynamical state of the neutral gas in HVCs. High-resolution observations (1.'1) of the 21 cm line emission in the EN field of the DHIGLS HI survey are used to analyze the physical properties of the bright concentration C I B located at an edge of complex C. We use the Gaussian decomposition code ROHSA to model its multiphase content, and perform a power spectrum analysis to analyze its multi-scale structure. Physical properties of some 200 structures extracted using dendrograms are examined. We identify two distinct regions, one of which has a prominent protrusion extending from the edge of complex C that exhibits an ongoing phase transition from warm diffuse gas to cold dense gas and filaments. The scale at which the warm gas becomes unstable and undergoes a thermal condensation is about 15 pc, corresponding to a cooling time about 1.5 Myr. We find that a transition from subsonic to trans-sonic turbulence is associated with the thermal condensation. A large scale perspective of complex C suggests that hydrodynamic instabilities are involved in creating the structured concentration C I B and the phase transition therein. However, the details of the dynamical and thermal processes remain unclear and will require further investigation, through both observations and numerical simulations. (Shortened for arxiv)

Zhong-Rui Bai, Hao-Tong Zhang, Hai-Long Yuan, Dong-Wei Fan, Bo-Liang He, Ya-Juan Lei, Yi-Qiao Dong, Si-Cheng Yu, Yong-Heng Zhao, Yong Zhang, Yong-Hui Hou, Yao-Quan Chu

21 pages, 15 figures

LAMOST Data Release 5, covering $\sim$17,000 $deg^2$ from $-10^{\circ}$ to $80^{\circ}$ in declination, contains 9 millions co-added low resolution spectra of celestial objects, each spectrum combined from repeat exposure of two to tens of times during Oct 2011 to Jun 2017. In this paper, We present the spectra of individual exposures for all the objects in LAMOST Data Release 5. For each spectrum, equivalent width of 60 lines from 11 different elements are calculated with a new method combining the actual line core and fitted line wings. For stars earlier than F type, the Balmer lines are fitted with both emission and absorption profiles once two components are detected. Radial velocity of each individual exposure is measured by minimizing ${\chi}^2$ between the spectrum and its best template. Database for equivalent widths of spectral lines and radial velocities of individual spectra are available online. Radial velocity uncertainties with different stellar type and signal-to-noise ratio are quantified by comparing different exposure of the same objects. We notice that the radial velocity uncertainty depends on the time lag between observations. For stars observed in the same day and with signal-to-noise ratio higher than 20, the radial velocity uncertainty is below 5km/s, and increase to 10km/s for stars observed in different nights.

S.C. Schuler, J.J. Andrews, V.R. Clanzy II, M. Mourabit, J. Chanamé, M.A. Agüeros

21 pages, 12 figures, seven tables; accepted to AJ

While most moving groups are young and nearby, a small number have been identified in the Galactic halo. Understanding the origin and evolution of these groups is an important piece of reconstructing the formation history of the halo. Here we report on our analysis of three putative halo moving groups: G03-37, G18-39, and G21-22. Based on Gaia EDR3 data, the stars associated with each group show some scatter in velocity (e.g., Toomre diagram) and integrals of motion (energy, angular momentum) spaces, counter to expectations of moving-group stars. We choose the best candidate of the three groups, G21-22, for follow-up chemical analysis based on high-resolution spectroscopy of six presumptive members. Using a new Python code that uses a Bayesian method to self-consistently propagate uncertainties from stellar atmosphere solutions in calculating individual abundances and spectral synthesis, we derive the abundances of $\alpha$- (Mg, Si, Ca, Ti), Fe-peak (Cr, Sc, Mn, Fe, Ni), odd-$Z$ (Na, Al, V), and neutron-capture (Ba, Eu) elements for each star. We find that the G21-22 stars are not chemically homogeneous. Based on the kinematic analysis for all three groups and the chemical analysis for G21-22, we conclude the three are not genuine moving groups. The case for G21-22 demonstrates the benefit of combining kinematic and chemical information in identifying conatal populations when either alone may be insufficient. Comparing the integrals of motion and velocities of the six G21-22 stars with those of known structures in the halo, we tentatively associate them with the Gaia-Enceladus accretion event.

Haruka Washinoue, Takeru K. Suzuki

11 pages, 5 figures. Accepted for publication in MNRAS

We systematically investigated the heating of coronal loops on metal-free stars with various stellar masses and magnetic fields by magnetohydrodynamic simulations. It is found that the coronal property is dependent on the coronal magnetic field strength $B_{\rm c}$ because it affects the difference of the nonlinearity of the Alfv\'{e}nic waves. Weaker $B_{\rm c}$ leads to cooler and less dense coronae because most of the input waves dissipate in the lower atmosphere on account of the larger nonlinearity. Accordingly EUV and X-ray luminosities also correlate with $B_{\rm c}$, while they are emitted in a wide range of the field strength. Finally we extend our results to evaluating the contribution from low-mass Population III coronae to the cosmic reionization. Within the limited range of our parameters on magnetic fields and loop lengths, the EUV and X-ray radiations give a weak impact on the ionization and heating of the gas at high redshifts. However, there still remains a possibility of the contribution to the reionization from energetic flares involving long magnetic loops.

S. R. Federman (1), Johnathan S. Rice (1), A. M. Ritchey (2), Hwihyun Kim (3), John H. Lacy (4), Paul F. Goldsmith (5), Nicolas Flagey (6 and 7), Gregory N. Mace (4), David L. Lambert (4) ((1) University of Toledo, (2) Eureka Scientific, (3) Gemini Observatory, (4) University of Texas at Austin, (5) Jet Propulsion Laboratory, (6) Canada France Hawaii Telescope, (7) Space Telescope Science Institute)

11 figures, 14 tables, 75 pages; higher quality maps appear in the published version

We study four lines of sight that probe the transition from diffuse molecular gas to molecular cloud material in Taurus. Measurements of atomic and molecular absorption are used to infer the distribution of species and the physical conditions toward stars behind the Taurus Molecular Cloud (TMC). New high-resolution spectra at visible and near infrared wavelengths of interstellar Ca II, Ca I, K I, CH, CH^+, C2, CN, and CO toward HD28975 and HD29647 are combined with data at visible wavelengths and published CO results from ultraviolet measurements for HD27778 and HD30122. Gas densities and temperatures are inferred from C2, CN, and CO excitation and CN chemistry. Our results for HD29647 are noteworthy because the CO column density is 10^{18} cm^{-2} while C2 and CO excitation reveals a temperature of 10 K and density about 1000 cm^{-3}, more like conditions found in dark molecular clouds. Similar results arise from our chemical analysis for CN through reactions involving observations of CH, C2, and NH. Enhanced potassium depletion and a reduced CH/H2 column density ratio also suggest the presence of a dark cloud. The directions toward HD27778 and HD30122 probe molecule-rich diffuse clouds, which can be considered CO-dark gas, while the sight line toward HD28975 represents an intermediate case. Maps of dust temperature help refine the description of the material along the four sight lines and provide an estimate of the distance between HD29647 and a clump in the TMC. An Appendix provides results for the direction toward HD26571; this star also probes diffuse molecular gas.

Takahito Sakaue, Kazunari Shibata

31 pages, 24 figures. accepted for publication in The Astrophysical Journal. for associated python code, see this https URL

M dwarf's atmosphere is expected to be highly magnetized. The magnetic energy can be responsible for heating the stellar chromosphere and corona, and driving the stellar wind. The nonlinear propagation of Alfv\'en wave is the promising mechanism for both heating stellar atmosphere and driving stellar wind. Based on this Alfv\'en wave scenario, we carried out the one-dimensional compressive magnetohydrodynamic (MHD) simulation to reproduce the stellar atmospheres and winds of TRAPPIST-1, Proxima Centauri, YZ CMi, AD Leo, AX Mic, as well as the Sun. The nonlinear propagation of Alfv\'en wave from the stellar photosphere to chromosphere, corona, and interplanetary space is directly resolved in our study. The simulation result particularly shows that the slow shock generated through the nonlinear mode coupling of Alfv\'en wave is crucially involved in both dynamics of stellar chromosphere (stellar spicule) and stellar wind acceleration. Our parameter survey further revealed the following general trends of physical quantities of stellar atmosphere and wind. (1) The M dwarfs' coronae tend to be cooler and denser than solar corona. (2) M dwarfs' stellar winds can be characterized with relatively faster velocity and much smaller mass-loss rate compared to those of solar wind. The physical mechanisms behind these tendencies are clarified in this paper, where the stronger stratification of M dwarf's atmosphere and relatively smaller Alfv\'en wave energy input from the M dwarf's photosphere are remarkable.

Cong-Ming Shi, Hui Deng, Feng Wang, Ying Mei, Shao-Guang Guo, Chen Yang, Chen Wu, Shou-Lin Wei, Andreas Wicenec

11 pages, 10 figures

Data archiving is one of the most critical issues for modern astronomical observations. With the development of a new generation of radio telescopes, the transfer and archiving of massive remote data have become urgent problems to be solved. Herein, we present a practical and robust file-level flow-control approach, called the Unlimited Sliding-Window (USW), by referring to the classic flow-control method in TCP protocol. Basing on the USW and the Next Generation Archive System (NGAS) developed for the Murchison Widefield Array telescope, we further implemented an enhanced archive system (ENGAS) using ZeroMQ middleware. The ENGAS substantially improves the transfer performance and ensures the integrity of transferred files. In the tests, the ENGAS is approximately three to twelve times faster than the NGAS and can fully utilize the bandwidth of network links. Thus, for archiving radio observation data, the ENGAS reduces the communication time, improves the bandwidth utilization, and solves the remote synchronous archiving of data from observatories such as Mingantu spectral radioheliograph. It also provides a better reference for the future construction of the Square Kilometer Array (SKA) Science Regional Center.

We present a high fidelity snapshot spectroscopic radio imaging study of a weak type I solar noise storm which took place during an otherwise exceptionally quiet time. Using high fidelity images from the Murchison Widefield Array, we track the observed morphology of the burst source for 70 minutes and identify multiple instances where its integrated flux density and area are strongly anti-correlated with each other. The type I radio emission is believed to arise due to electron beams energized during magnetic reconnection activity. The observed anti-correlation is interpreted as evidence for presence of MHD sausage wave modes in the magnetic loops and strands along which these electron beams are propagating. Our observations suggest that the sites of these small scale reconnections are distributed along the magnetic flux tube. We hypothesise that small scale reconnections produces electron beams which quickly get collisionally damped. Hence, the plasma emission produced by them span only a narrow bandwidth and the features seen even a few MHz apart must arise from independent electron beams.

Cunshi Wang, Yu Bai, Haibo Yuan, Song Wang, Jifeng Liu

19 pages,31 figures, accepted by AJ

In modern astronomy, machine learning as an raising realm for data analysis, has proved to be efficient and effective to mine the big data from the newest telescopes. By using support vector machine (SVM), we construct a supervised machine learning algorithm, to classify the objects in the Javalambre-Photometric Local Universe Survey (J-Plus). The sample is featured with 12-waveband, and magnitudes is labeled with spectrum-based catalogs, including Sloan Digital Sky Survey spectroscopic data, Large Sky Area Multi-Object Fiber Spectroscopic Telescope, and VERONCAT - Veron Catalog of Quasars & AGN. The performance of the classifier is presented with the applications of blind test validations based on RAdial Velocity Extension, Kepler Input Catalog, 2 MASS Redshift Survey, and UV-bright Quasar Survey. The accuracies of the classifier are 96.5% in blind test and 97.0\% in training cross validation. The F_1-scores are 95.0% for STAR, 92.9% for GALAXY and 87.0% for QSO. In the classification for J-Plus catalog, we develop a new method to constrain the potential extrapolation.

D. A. Ladeyschikov, M. S. Kirsanova, A. M. Sobolev, M. Thomasson, V. Ossenkopf-Okada, M. Juvela, S. A. Khaibrakhmanov, E. A. Popova

Accepted for publication in MNRAS

The paper aims to study relation between the distributions of the young stellar objects (YSOs) of different ages and the gas-dust constituents of the S254-S258 star-formation complex. This is necessary to study the time evolution of the YSO distribution with respect to the gas and dust compounds which are responsible for the birth of the young stars. For this purpose we use correlation analysis between different gas, dust and YSOs tracers. We compared the large-scale CO, HCO$^+$, near-IR extinction, and far-IR {\it Herschel} maps with the density of YSOs of the different evolutionary Classes. The direct correlation analysis between these maps was used together with the wavelet-based spatial correlation analysis. This analysis reveals a much tighter correlation of the gas-dust tracers with the distribution of Class I YSOs than with that of Class II YSOs. We argue that Class I YSOs which were initially born in the central bright cluster S255-IR (both N and S parts) during their evolution to Class II stage ($\sim$2 Myr) had enough time to travel through the whole S254-S258 star-formation region. Given that the region contains several isolated YSO clusters, the evolutionary link between these clusters and the bright central S255-IR (N and S) cluster can be considered. Despite the complexity of the YSO cluster formation in the non-uniform medium, the clusters of Class II YSOs in the S254-258 star-formation region can contain objects born in the different locations of the complex.

Yuhu Miao, Libo Fu, Xian Du, Ding Yuan, Chaowei Jiang, Jiangtao Su, Mingyu Zhao, Sergey Anfinogentov

MNRAS, 5 pages,3 figures

A light bridge is a magnetic intrusion into a sunspot, it interacts with the main magnetic field and excites a variety of dynamical processes. In the letter, we studied magnetic connectivity between a light bridge and coronal loops rooted at the sunspot. We used the data of the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory (SDO) to study the features of sunspots with light bridges. It is found that if a light bridge anchors at the umbra-penumbra boundary, the coronal loops could not be formed around the anchoring point. If the a light bridge become detached from the penumbra, the coronal loop starts to form again. The vector magnetogram provided by the Helioseismic Magnetic Imager onboard SDO shows that the anchoring region of a light bridge usually have an accompanying opposite minor-polarities. We conjugate that the magnetic field line could connect to these opposite polarities and form short-range magnetic loops, and therefore, coronal loops that extend to long-range could not be formed. A model of light bridge is proposed to explain the magnetic connectivity between a light bridge and the coronal loops. This model could explain many physical processes associated with light bridges.

Ze-Cheng Zou, Bin-Bin Zhang, Yong-Feng Huang, Xiao-Hong Zhao

10 pages, 6 figures, 1 table, submitted to ApJL

Typical models of gamma-ray bursts (GRBs) assume that they happen in isolated central engine systems, with little attention being paid to the fact that most stars (including those evolving as GRB progenitors) are in binary or multi-body systems. This paper will study in detail the observational effects when a GRB occurs in a binary system. We show that the blockage by the companion star becomes nonnegligible when it is located within a typical GRB jet opening angle (e.g., 10 degrees) and at an orbit of 67 au, assuming the GRB emission site is at $10^{15}$ cm radius. In such a case, an on-axis observer will see a GRB with a similar temporal behavior but 25% dimmer. On the other hand, an off-axis observer outside the jet opening angle (hence missed the original GRB) can see a delayed "reflected" GRB, which is much fainter in brightness, much wider in the temporal profile and slightly softer in energy. Our study can naturally explain the origin of some low-luminosity GRBs. Moreover, we also point out that the companion star may be shocked if it is located inside the GRB emission site, which can give rise to an X-ray transient or a GRB followed by a delayed X-ray bump on top of X-ray afterglows.

We study the redshift-space fluctuations induced by a stochastic gravitational wave background via the Sachs-Wolfe effect. The redshift-space fluctuations can be encapsulated in a line-of-sight integral that is useful for studying the imprint of short-wavelength gravitational waves on the CMB anisotropy, as well as providing us with a precise redshift fluctuation correlation between a pair of pulsars in pulsar timing measurements.

Lami Suleiman, Morgane Fortin, Julian Leszek Zdunik, Pawel Haensel

19 pages, 12 figures, 3 tables. Accepted for publication in PRC

Measurements of neutron-star macrophysical properties thanks to multimessenger observations offer the possibility to constrain the properties of nuclear matter. Indeed cold and dense matter as found inside neutron stars, in particular in their core, is not accessible to terrestrial laboratories.We investigate the consequences of using equations of state that employ models for the core and the crust that are not calculated consistently on the neutron-star macrophysical properties, on some of the so-called universal relations and on the constraints obtained from gravitational-wave observations. We use various treatments found in the literature to connect together nonconsistent core and crust equations of state. We then compute the mass, the radius, the tidal deformability, and the moment of inertia for each model. Finally, we assess the discrepancies in the neutron-star macrophysical properties obtained when consistent models for the whole star and nonconsistent ones are employed. The use of crust models nonconsistent with the core introduces an error on the macrophysical parameters which can be as large as the estimated accuracy of current and next generation telescopes. The precision of some of the universal relations reported in the literature is found to be overestimated. We confirm that the equation of the crust has limited influence on the macrophysical properties. The discrepancy between results obtained for a fully consistent equation of state and a nonconsistent one can be reduced if one connects the core and the crust models at baryon densities around 0.08-0.1 fm$^{-3}$. The equation of the crust cannot be probed with current multimessenger observations and near-future ones.

Ogata Erika, Ohsuga Ken, Yajima Hidenobu

14 pages, 8 figures, accepted for publication in PASJ

We investigate the Hoyle-Lyttleton accretion of dusty-gas for the case where the central source is the black hole accretion disk. By solving the equation of motion taking into account the radiation force which is attenuated by the dust absorption, we reveal the steady structure of the flow around the central object. We find that the mass accretion rate tends to increase with an increase of the optical thickness of the flow and the gas can accrete even if the disk luminosity exceeds the Eddington luminosity for the dusty-gas, since the radiation force is weakened by the attenuation via the dust absorption. When the gas flows in from the direction of the rotation axis for the disk with ${\Gamma}^{'}=3.0$, the accretion rate is about 93% of the Hoyle-Lyttleton accretion rate if ${\tau}_{\rm{HL}}=3.3$ and zero for ${\tau}_{\rm{HL}}=1.0$, where ${\Gamma}^{'}$ is the Eddington ratio for the dusty-gas and ${\tau}_{\rm{HL}}$ is the typical optical thickness of the Hoyle-Lyttleton radius. Since the radiation flux in the direction of disk plane is small, the radiation force tends not to prevent gas accretion from the direction near the disk plane. For ${\tau}_{\rm{HL}}=3.3$ and ${\Gamma}^{'}=3.4$, although the accretion is impossible in the case of ${\Theta}=0$, the accretion rate is 28% of the Hoyle-Lyttleton one in the case of ${\Theta}=90$, where ${\Theta}$ is the angle between the direction the gas is coming from and the rotation axis of the disk. We also obtain relatively high accretion luminosity that is realized when the accretion rate of the disk onto the BH is consistent with that via the Hoyle-Lyttleton mechanism taking into account the effect of radiation. This implies the intermediate-mass black holes moving in the dense dusty-gas are identified as luminous objects in the infrared band.

The r-band of the Sloan Digital Sky Survey (SDSS) for 17,924 brightest cluster galaxies (BCGs) in clusters and groups within 0.02 z 0.20 are used to study possible environmental relations affecting the nature of these galaxies. We find a correlation between BCGs physical properties (the effective radius (Re ), absolute magnitude and central velocity dispersion ({\sigma}0 )) and their host groups and clusters velocity dispersion ({\sigma}cl ). This type of relations suggests that the most massive groups or clusters host larger central galaxies. On the other hand, the {\sigma}0 /{\sigma}cl ratio as a function of {\sigma}cl is consistent with [10].

It is shown with several concrete examples that the Dyer-Roeder approximation is valid in spacetimes which fulfill the condition that fluctuations in the expansion rate along a light ray locally cancels with the shear contribution to the redshift. This is the case for standard cosmological scenarios including perturbed FLRW spacetimes, N-body simulations and Swiss-cheese models. With another concrete example it is then illustrated that it is possible to construct statistically homogeneous and effectively statistically isotropic cosmological models which do not fulfill the condition. In this case, the Dyer-Roeder approximation is invalid. Instead, the mean redshift-distance relation can be described using a relation based on the spatial averages of the transparent part of the spacetime.

Jackson M. Norris, Sowgat Muzahid, Jane C. Charlton, Glenn G. Kacprzak, Bart P. Wakker, Christopher W. Churchill

18 pages, 9 figures

We have obtained new observations of the absorption system at $z_\mathrm{abs}=0.48$ toward QSO Q0454-220, which we use to constrain its chemical and physical conditions. The system features metal-enriched gas and previously unknown low-metallicity gas detected $\sim 200 \, \mathrm{km \, s^{-1}}$ blueward of the metal-enriched gas. The low-metallicity gas is detected in multiple Lyman series lines but is not detected in any metal lines. Our analysis includes low-ionization (e.g., Fe II, Mg II) metal lines, high-ionization (e.g., C IV, O VI, N V) metal lines, and several Lyman series lines. We use new UV spectra taken with HST/COS along with data taken from HST/STIS, Keck/HIRES, and VLT/UVES. We find that the absorption system can be explained with a photoionized low-ionization phase with $\mathrm{[Fe/H]} \sim -0.5$ and $n_\mathrm{H} \sim 10^{-2.3} \, \mathrm{cm}^{-3}$, a photoionized high-ionization phase with a conservative lower limit of $-3.3 < \mathrm{[Fe/H]}$ and $n_\mathrm{H} \sim 10^{-3.8} \, \mathrm{cm}^{-3}$, and a low-metallicity component with a conservative upper limit of $\mathrm{[Fe/H]} < -2.5$ that may be photoionized or collisionally ionized. We suggest that the low-ionization phase may be due to cold-flow accretion via large-scale filamentary structure or due to recycled accretion while the high-ionization phase is the result of ancient outflowing material from a nearby galaxy. The low-metallicity component may come from pristine accretion. The velocity spread and disparate conditions among the absorption system's components suggest a combination of gas arising near galaxies along with gas arising from intergroup material.

E.Kokoulina, A. Matter, B. Lopez, E.Pantin, N.Ysard, G.Weigelt, E.Habart, J.Varga, A.Jones, A.Meilland, E.Dartois, L.Klarmann, J.-C. Augereau, R.van Boekel, M.Hogerheijde, G.Yoffe, L.B.F.M.Waters, C.Dominik, W.Jaffe, F.Millour, Th.Henning, K.-H.Hofmann, D.Schertl, S.Lagarde, R.G.Petrov, P.Antonelli, F.Allouche, P.Berio, S.Robbe-Dubois, P.Abraham, U.Beckmann, A.Bensberg, F.Bettonvil, P.Bristow, P.Cruzalèbes, W.C.Danchi, M.Dannhoff, U.Graser, M.Heininger, L.Labadie, M.Lehmitz, C.Leinert, K.Meisenheimer, C.Paladini, I.Percheron, Ph.Stee, J.Woillez, S.Wolf, G.Zins, M.Delbo, J.Drevon, J.Duprat, V.Gámez Rosas, V.Hocdé, J.Hron, C.A.Hummel, J.W.Isbell, J.Leftley, A.Soulain, F.Vakili, M.Wittkowski

Carbon is one of the most abundant components in the Universe. While silicates have been the main focus of solid phase studies in protoplanetary discs (PPDs), little is known about the solid carbon content especially in the planet-forming regions ($\sim $0.1 to 10 au). Fortunately, several refractory carbonaceous species present C-H bonds (such as hydrogenated nano-diamond and amorphous carbon as well as polycyclic aromatic hydrocarbons (PAHs)), which generate infrared (IR) features that can be used to trace the solid carbon reservoirs. The new mid-IR instrument MATISSE, installed at the Very Large Telescope Interferometer (VLTI), can spatially resolve the inner regions ($\sim$ 1 to 10 au) of PPDs and locate, down to the au-scale, the emission coming from carbon grains. Our aim is to provide a consistent view on the radial structure, down to the au-scale, as well as basic physical properties and the nature of the material responsible for the IR continuum emission in the inner disk region around HD 179218. We implemented a temperature-gradient model to interpret the disk IR continuum emission, based on a multiwavelength dataset comprising a broadband spectral energy distribution (SED) and VLTI H-, L-, and N-bands interferometric data obtained in low spectral resolution. Then, we added a ring-like component, representing the carbonaceous L-band features-emitting region, to assess its detectability in future higher spectral resolution observations employing mid-IR interferometry.

Weiss W.W, Zwintz K., Kuschnig R., Handler G., Moffat A.F.J., Baade D., Bowman D.M., Granzer T., Kallinger T., Koudelka O.F., Lovekin C.C., Neiner C., Pablo H., Pigulski A., Popowicz A., Ramiaramanantsoa T., Rucinski S.M., Strassmeier K.G., Wade G.A

BRITE-Constellation was designed, built, launched, and is operated and supported by the Austrian Research Promotion Agency, the University of Vienna, the Technical University of Graz, the University of Innsbruck, the Canadian Space Agency, the University of Toronto Institute for Aerospace Studies, the Foundation for Polish Science and Technology and the National Science Centre. 23 figures

BRITE-Constellation is devoted to high-precision optical photometric monitoring of bright stars, distributed all over the Milky Way, in red and/or blue passbands. Photometry from space avoids the turbulent and absorbing terrestrial atmosphere and allows for very long and continuous observing runs with high time resolution and thus provides the data necessary for understanding various processes inside stars (e.g., asteroseismology) and in their immediate environment. While the first astronomical observations from space focused on the spectral regions not accessible from the ground it soon became obvious around 1970 that avoiding the turbulent terrestrial atmosphere significantly improved the accuracy of photometry and satellites explicitly dedicated to high-quality photometry were launched. A perfect example is BRITE-Constellation, which is the result of a very successful cooperation between Austria, Canada and Poland. Research highlights for targets distributed nearly over the entire HRD are presented, but focus primarily on massive and hot stars.

Ambient radiation background contributed by the penetrating cosmic ray particles and the radionuclides present in the rock materials have been measured at an underground laboratory located inside a mine at 555 m depth. The laboratory is being set up to explore rare event search processes, such as direct dark matter search, neutrinoless double beta decay, axion search, supernova neutrino detection, etc., that require specific knowledge of the nature and extent of the radiation environment in order to assess the sensitivity reach and also to plan for its reduction for the targeted experiment. The gamma ray background, which is mostly contributed by the primordial radionuclides and their decay chain products, have been measured inside the laboratory and found to be dominated by rock radioactivity for $E_\gamma \lesssim 3 \,{\rm MeV}$. Shielding of these residual gamma rays for the experiment was also evaluated. The cosmic muon flux, measured inside the laboratory using large area plastic scintillator telescope, was found to be: $(2.051 \pm 0.142 \pm 0.009) \times 10^{-7}\, {\rm cm}^{-2}.{\rm sec}^{-1}$, which agrees reasonably well with simulation results. The neutron background flux has been measured for the radiogenic neutrons and found to be: $(1.61 \pm 0.03) \times 10^{-4} \, {\rm cm}^{-2}.{\rm sec}^{-1}$ for no threshold cut. Detailed GEANT4 simulation for the radiogenic neutrons and the cosmogenic neutrons have been carried out. Effects of multiple scattering of both the types of neutrons within the surrounding rock and the cavern walls were studied and the results for the radiogenic neutrons are found to be in reasonable agreement with experimental results. Neutron fluxes contributed by those neutrons of cosmogenic origin have been reported as function of the energy threshold.

Francesco Guarneri, Giorgio Calderone, Stefano Cristiani, Fabio Fontanot, Konstantina Boutsia, Guido Cupani, Andrea Grazian, Valentina D'Odorico

Accepted for publication in MNRAS, 12 pages, 11 figures, 4 tables

The number of known, bright ($i<18$), high-redshift ($z>2.5$) QSOs in the Southern Hemisphere is considerably lower than the corresponding number in the Northern Hemisphere due to the lack of multi-wavelength surveys at $\delta<0$. Recent works, such as the QUBRICS survey, successfully identified new, high-redshift QSOs in the South by means of a machine learning approach applied on a large photometric dataset. Building on the success of QUBRICS, we present a new QSO selection method based on the Probabilistic Random Forest (PRF), an improvement of the classic Random Forest algorithm. The PRF takes into account measurement errors, treating input data as probability distribution functions: this allows us to obtain better accuracy and a robust predictive model. We applied the PRF to the same photometric dataset used in QUBRICS, based on the SkyMapper DR1, Gaia DR2, 2MASS, WISE and GALEX databases. The resulting candidate list includes $626$ sources with $i<18$. We estimate for our proposed algorithm a completeness of $\sim84\%$ and a purity of $\sim78\%$ on the test datasets. Preliminary spectroscopic campaigns allowed us to observe 41 candidates, of which 29 turned out to be $z>2.5$ QSOs. The performances of the PRF, currently comparable to those of the CCA, are expected to improve as the number of high-z QSOs available for the training sample grows: results are however already promising, despite this being one of the first applications of this method to an astrophysical context.

This study is the direct continuation of the work performed in Hillman et al. (2020) where they used their feedback dominated numerical simulations to model the evolution of four initial models with white dwarf (WD) masses of 0.7 and 1.0M_Solar and red dwarf (RD) masses of 0.45 and 0.7M_Solar from first Roche-lobe contact of the donor RD, over a few times 10^9 years, until the RD was eroded down to below 0.1M_Solar. This study presents an in-depth analysis of their four models complimented by three models with a higher WD mass of 1.25M_Solar, one of which comprises an oxygen-neon (ONe) core. Common features were found for all seven models on a secular time scale as well as on a cyclic time scale. On the other hand, certain features were found that are strongly dependent either on the WD or the RD mass but are indifferent to the other of the two. Additionally, a model with a WD composed of an ONe core was compared with its corresponding carbon oxygen (CO) core WD model and found to have a significant impact on the heavy element abundances in the ejecta composition.

J. Healy, T. Deb, M.A.W. Verheijen, S-L. Blyth, P. Serra, M. Ramatsoku, B. Vulcani

20 pages, 11 figures, accepted for publication in Astronomy and Astrophysics

The morphology-density relation manifests the environmental dependence of the formation and evolution of galaxies as they continuously migrate through the cosmic web to ever denser environments. As gas-rich galaxies traverse the outskirts and inner regions of galaxy clusters they experience sudden and radical changes in their gas content and star formation activity. The goal of this work is to gain an H$\,$I perspective on gas depletion mechanisms acting on galaxies and galaxy groups that are being accreted by a moderately massive galaxy cluster. We aim to study the relative importance and efficiency of processes such as ram-pressure stripping and tidal interactions as well as their dependency on the local and global environment of galaxies in the cluster core and in its surroundings. We have conducted a blind radio continuum and H$\,$I spectral line imaging survey with the MeerKAT radio telescope of a 2$^\circ$ $\times$ 2$^\circ$ area centred on the galaxy cluster Abell 2626. We have used the CARAcal pipeline to reduce the data, SoFiA to detect sources within the H$\,$I data cube, and GIPSY to construct spatially resolved information on the H$\,$I morphologies and kinematics of the H$\,$I detected galaxies. We have detected H$\,$I in 219 galaxies with optical counterparts within the entire surveyed volume. We present the H$\,$I properties of each of the detected galaxies as a data catalogue and as an atlas page for each galaxy, including H$\,$I column-density maps, velocity fields, position-velocity diagrams and global H$\,$I profiles. These data will also be used for case studies of identified ``jellyfish'' galaxies and galaxy population studies by means of morphological classification of the direct H$\,$I detections as well as using the H$\,$I stacking technique.

Jaime Vargas-González, Jan Forbrich, Sergio A. Dzib, John Bally

Accepted for publication in MNRAS, 17 pages, 12 figures

We present a newly enlarged census of the compact radio population towards the Orion Nebula Cluster (ONC) using high-sensitivity continuum maps (3-10 $\mu$Jy bm$^{-1}$) from a total of $\sim30$ h centimeter-wavelength observations over an area of $\sim$20$'\times20'$ obtained in the C-band (4$-$8 GHz) with the Karl G. Jansky Very Large Array (VLA) in its high-resolution A-configuration. We thus complement our previous deep survey of the innermost areas of the ONC, now covering the field of view of the Chandra Orion Ultra-deep Project (COUP). Our catalog contains 521 compact radio sources of which 198 are new detections. Overall, we find that 17% of the (mostly stellar) COUP sources have radio counterparts, while 53% of the radio sources have COUP counterparts. Most notably, the radio detection fraction of X-ray sources is higher in the inner cluster and almost constant for $r>3'$ (0.36 pc) from $\theta^1$ Ori C suggesting a correlation between the radio emission mechanism of these sources and their distance from the most massive stars at the center of the cluster, for example due to increased photoionisation of circumstellar disks. The combination with our previous observations four years prior lead to the discovery of fast proper motions of up to $\sim$373 km s$^{-1}$ from faint radio sources associated with ejecta of the OMC1 explosion. Finally, we search for strong radio variability. We found changes in flux density by a factor of $\lesssim$5 within our observations and a few sources with changes by a factor $>$10 on long timescales of a few years.

Sara Webb, Chris Flynn, Jeff Cooke, Jielai Zhang, Ashish Mahabal, Tim Abbott, Rebecca Allen, Igor Andreoni, Sarah Bird, Simon Goode, Michelle Lochner, Tyler Pritchard

20 pages, 14 figures, 6 tables

We present our 500 pc distance-limited study of stellar fares using the Dark Energy Camera as part of the Deeper, Wider, Faster Program. The data was collected via continuous 20-second cadence g band imaging and we identify 19,914 sources with precise distances from Gaia DR2 within twelve, ~3 square-degree, fields over a range of Galactic latitudes. An average of ~74 minutes is spent on each field per visit. All light curves were accessed through a novel unsupervised machine learning technique designed for anomaly detection. We identify 96 flare events occurring across 80 stars, the majority of which are M dwarfs. Integrated are energies range from $\sim 10^{31}-10^{37}$ erg, with a proportional relationship existing between increased are energy with increased distance from the Galactic plane, representative of stellar age leading to declining yet more energetic are events. In agreement with previous studies we observe an increase in flaring fraction from M0 -> M6 spectral types. Furthermore, we find a decrease in the flaring fraction of stars as vertical distance from the galactic plane is increased, with a steep decline present around ~100 pc. We find that ~70% of identified flares occur on short timescales of ~8 minutes. Finally we present our associated are rates, finding a volumetric rate of $2.9 \pm 0.3 \times 10^{-6}$ flares pc$^{-3}$ hr$^{-1}$.

James O. Chibueze, Haruka Sakemi, Takumi Ohmura, Mami Machida, Hiroki Akamatsu, Takuya Akahori, Hiroyuki Nakanishi, Viral Parekh, Ruby van Rooyen, Tsutomu T. Takeuchi

4 figures

Galaxy clusters are known to harbour magnetic fields. The nature of the intra-cluster magnetic fields remains an unresolved question. Intra-cluster magnetic field can be observed at the density contact discontinuity formed by cool and dense plasma running into hot ambient plasma, and the discontinuity exists near the 2nd BCG MRC 0600-399 of a merging galaxy cluster Abell 3376 (z=0.0461, hereafter as A3376). Elongated X-ray image in the east-west direction with a comet-like structure reaches a Mpc-scale (Fig1.(a)). Previous radio observations detected the bent jets from MRC 0600-399, moving in same direction as the sub-cluster's motion against ram pressure.Here we report a new radio observation of a radio galaxy MRC 0600-399 which is 3.4 times and 11 times higher resolution and sensitivity than the previous results. Contrary to typical jets, the MRC 0600-399 shows a 90deg bend at the contact discontinuity and the collimated jets further extend over 100 kpc from the bend point. Diffuse, elongated emission named "double-scythe" structures were detected for the first time. The spectral index flattens downstream of the bend point, indicating cosmic-ray re-acceleration. High-resolution numerical simulations reveal that the ordered magnetic field along the discontinuity plays a significant role in the change in the jet direction. The morphology of the "double-scythe" bear remarkable similarities with the simulations, which strengthens our understanding of the interaction between relativistic electrons and intra-cluster magnetic field.

Fernando Moreno, Javier Licandro, Antonio Cabrera-Lavers, David Morate, Daniel Guirado

Accepted by MNRAS (June 24th, 2021)

We report on the characterisation of the dust activity and dynamical evolution of two faint active asteroids, P/2019 A4, and P/2021 A5, observed with the 10.4m GTC using both imaging and spectroscopy. Asteroid P/2019 A4 activity is found to be linked to an impulsive event occurring some $\pm$10 days around perihelion, probably due to a collision or a rotational disruption. Its orbit is stable over 100 Myr timescales. Dust tail models reveal a short-term burst producing (2.0$\pm$0.7)$\times$10$^6$ kg of dust for maximum particle radius rmax=1 cm. The spectrum of P/2019 A4 is featureless, and slightly redder than the Sun. P/2021 A5 was active $\sim$50 days after perihelion, lasting $\sim$5 to $\sim$60 days, and ejecting (8$\pm$2)$\times$10$^6$ kg of dust for rmax=1 cm. The orbital simulations show that a few percent of dynamical clones of P/2021 A5 are unstable on 20-50 Myr timescales. Thus, P/2021 A5 might be an implanted object from the JFC region or beyond. These facts point to water ice sublimation as the activation mechanism. This object also displays a featureless spectrum, but slightly bluer than the Sun. Nuclei sizes are estimated in the few hundred meters range for both asteroids. Particle ejection speeds ($\sim$0.2 m/s) are consistent with escape speeds from those small-sized objects.

Simón Díaz-García, Ute Lisenfeld, Isabel Pérez, Almudena Zurita, Simon Verley, Françoise Combes, Daniel Espada, Stéphane Leon, Vicent Martínez-Badenes, José Sabater, Lourdes Verdes-Montenegro

Accepted for publication in A&A (June 22, 2021). 18 pages, 9 figures, 2 tables, plus appendices (16 pages, 18 figures, 3 tables)

While some galactic bars show recent massive star formation (SF) along them, some others present a lack of it. Whether bars with low level of SF are a consequence of low star formation efficiency (SFE), low gas inflow rate, or dynamical effects, remains a matter of debate. We perform a multi-wavelength analysis of 12 strongly barred massive galaxies, chosen to host different degrees of SF along the bar major axis without any prior condition on gas content. We observe the CO(1-0) and CO(2-1) emission within bars with the IRAM-30m telescope, which we use to estimate molecular gas masses. SF rates (SFR) are calculated from GALEX near- and far- ultraviolet (UV) and WISE 12 and 22 micron images within the beam pointings, covering the full bar extent. We detect molecular gas along the bars of all probed galaxies. The SFE in bars varies between galaxies by up to an order of magnitude. On average, SFEs are roughly constant along bars. SFEs are not significantly different from the mean value in spiral galaxies reported in the literature. Interestingly, the higher the total stellar mass of the host galaxy, the lower the SFE within their bars. In particular, the two galaxies in our sample with lowest SFEs and SFR surface densities (NGC 4548 and NGC 5850) are also the ones hosting massive bulges and signs of past interactions with nearby companions. The SFE in strong bars is not systematically inhibited (either in the central, mid- or end-parts of the bar). Both environmental and internal quenching are likely responsible for the lowest SFEs reported in this work (Abridged).

Florian Sarron, Christopher J Conselice

Accepted in MNRAS. 21 pages (18 + 3 in appendix). 9 figures. Include online only appendix of the published version

We use a large K-selected sample of 299,961 galaxies from the REFINE survey, consisting of a combination of data from three of the deepest near-infrared surveys: UKIDSS UDS, COSMOS/UltraVISTA and CFHTLS-D1/VIDEO, that were homogeneously reduced to obtain photometric redshifts and stellar masses. We detect 2588 candidate galaxy groups up to $z=3.15$ at $S/N>1.5$. We build a very pure ($>90\%$) sample of 448 candidate groups up to $z=2.5$ and study some of their properties. Cluster detection is done using the DElaunay TEssellation ClusTer IdentiFication with photo-z (DETECTIFz) algorithm that we describe. This new group finder algorithm uses the joint probability distribution functions (PDF) of redshift and stellar-mass of galaxies to detect groups as stellar-mass overdensities in overlapping redshift slices, where density is traced using Monte Carlo realisation of the Delaunay Tessellation Field Estimator (DTFE). We compute the algorithm selection function using mock galaxy catalogues taken from cosmological N-body simulation lightcones. Based on these simulations, we reach a completeness of $\sim80\%$ for clusters ($M_{200}>10^{14} M_{\odot}$) at a purity of $\sim90\%$ at $z<2.5$. Using our 403 most massive candidate groups, we constrain the redshift evolution of the group galaxy quenched fraction at $0.12\le z<2.32$, for galaxies with $10.25 < \log M_\star/M_{\odot} < 11$ in $0.5\times R_{200}$. We find that the quenched fraction in group cores is higher than in the field in the full redshift range considered, the difference growing with decreasing redshift. This indicates either more efficient quenching mechanisms in group cores at lower redshift or pre-processing by cosmic filaments.

Itamar Reis, Rennan Barkana, Anastasia Fialkov

Submitted to MNRAS on June 7th

The 21-cm signal from cosmic dawn and the epoch of reionization (EoR) probes the characteristics of the high redshift galaxy population. Many of the astrophysical properties of galaxies at high redshifts are currently unconstrained due to the lack of observations. This creates a vast space of possible astrophysical scenarios where the 21-cm signal needs to be modeled in order to plan for, and eventually fit, future observations. This is done with fast numerical methods which make simplifying approximations for the underlying physical processes. In this work we quantify the effect of Poisson fluctuations and scatter in the star formation efficiency; while Poisson fluctuations are included in some works and not in others, scatter in the star formation efficiency is usually neglected, and all galaxies of a given mass are assumed to have the same properties. We show that both features can have a significant effect on the 21-cm power spectrum, most importantly in scenarios where the signal is dominated by massive galaxies. Scatter in the star formation efficiency does not simply enhance the effect of Poisson fluctuations; for example we show that the power spectrum shape at cosmic dawn has a feature corresponding to the width of the galaxy brightness distribution. We also discuss some of the consequences for 21-cm imaging, and the signature of reduced correlation between the density and radiation fields.

Pauline Zarrouk, Omar Ruiz-Macias, Shaun Cole, Peder Norberg, Carlton Baugh, David Brooks, Enrique Gaztañaga, Ellie Kitanidis, Robert Kehoe, Martin Landriau, John Moustakas, Francisco Prada, Gregory Tarlé

16 pages, 17 figures

We characterise the selection cuts and clustering properties of a magnitude-limited sample of bright galaxies that is part of the Bright Galaxy Survey (BGS) of the Dark Energy Spectroscopic Instrument (DESI) using the ninth data release of the Legacy Imaging Surveys (DR9). We describe changes in the DR9 selection compared to the DR8 one as explored in Ruiz-Macias et al. (2021). We also compare the DR9 selection in three distinct regions: BASS/MzLS in the north Galactic Cap (NGC), DECaLS in the NGC, and DECaLS in the south Galactic Cap (SGC). We investigate the systematics associated with the selection and assess its completeness by matching the BGS targets with the Galaxy and Mass Assembly (GAMA) survey. We measure the angular clustering for the overall bright sample (r $\leq$ 19.5) and as function of apparent magnitude and colour. This enables to determine the clustering strength and slope by fitting a power-law model that can be used to generate accurate mock catalogues for this tracer. We use a counts-in-cells technique to explore higher-order statistics and cross-correlations with external spectroscopic data sets in order to check the evolution of the clustering with redshift and the redshift distribution of the BGS targets using clustering-redshifts. While this work validates the properties of the BGS bright targets, the final target selection pipeline and clustering properties of the entire DESI BGS will be fully characterised and validated with the spectroscopic data of Survey Validation.

Francesca Matteucci

93 pages,37 figures. Accepted in Astronomy & Astrophysics Review

In this review I will discuss the comparison between model results and observational data for the Milky Way, the predictive power of such models as well as their limits. Such a comparison, known as Galactic archaeology, allows us to impose constraints on stellar nucleosynthesis and timescales of formation of the various Galactic components (halo, bulge, thick disk and thin disk).

N. Jurlin, M. Brienza, R. Morganti, Y. Wadadekar, C. H. Ishwara-Chandra, N. Maddox, V. Mahatma

24 pages, 7 figures, 7 tables, A&A accepted

Remnant radio galaxies represent an important phase in the life-cycle of radio active galactic nuclei. It is suggested that in this phase, the jets have switched off and the extended emission is fading rapidly. This phase is not well-studied due to the lack of statistical samples observed at both low and high frequencies. In this work, we study a sample of 23 candidate remnant radio galaxies previously selected using the Low Frequency Array at 150 MHz in the Lockman Hole field. We examine their morphologies and study their spectral properties to confirm their remnant nature and revise the morphological and spectral criteria used to define the initial sample. We present new observations with the Karl G. Jansky Very Large Array at 6000 MHz at both high and low resolution. These observations allowed us to observe the presence or absence of cores and study the spectral curvature and steepness of the spectra of the total emission expected at these high frequencies for the remnant candidates. We confirm 13 out of 23 candidates as remnant radio sources. This corresponds to 7% of the full sample of active, restarted, and remnant candidates from the Lockman Hole field. Surprisingly, only a minority of remnants reside in a cluster (23%). The remnant radio galaxies show a range of properties and morphologies. The majority do not show detection of the core at 6000 MHz and their extended emission often shows ultra-steep spectra (USS). However, there are also remnants with USS total emission and a detection of the core at 6000 MHz, possibly indicating a variety of evolutionary stages in the remnant phase. We confirm the importance of the combination of morphological and spectral criteria and this needs to be taken into consideration when selecting a sample of remnant radio sources.

Pedro Carrilho, Chiara Moretti, Benjamin Bose, Katarina Markovič, Alkistis Pourtsidou

30 pages, 10 figures

Interacting dark energy models have been proposed as attractive alternatives to $\Lambda$CDM. Forthcoming Stage-IV galaxy clustering surveys will constrain these models, but they require accurate modelling of the galaxy power spectrum multipoles on mildly non-linear scales. In this work we consider a dark scattering model with a simple 1-parameter extension to $w$CDM - adding only $A$, which describes a pure momentum exchange between dark energy and dark matter. We then provide a comprehensive comparison of three approaches of modeling non-linearities, while including the effects of this dark sector coupling. We base our modeling of non-linearities on the two most popular perturbation theory approaches: TNS and EFTofLSS. To test the validity and precision of the modelling, we perform an MCMC analysis using simulated data corresponding to a $\Lambda$CDM fiducial cosmology and Stage-IV surveys specifications in two redshift bins, $z=0.5$ and $z=1$. We find the most complex EFTofLSS-based model studied to be better suited at both, describing the mock data up to smaller scales, and extracting the most information. Using this model, we forecast uncertainties on the dark energy equation of state, $w$, and on the interaction parameter, $A$, finding $\sigma_w=0.06$ and $\sigma_A=1.1$ b/GeV for the analysis at $z=0.5$ and $\sigma_w=0.06$ and $\sigma_A=2.0$ b/GeV for the analysis at $z=1$. In addition, we show that a false detection of exotic dark energy up to 3$\sigma$ would occur should the non-linear modelling be incorrect, demonstrating the importance of the validation stage for accurate interpretation of measurements.

E.W. Guenther, D. Woeckel, P. Chaturvedi, V. Kumar, M.K. Srivastava, P. Muheki

14 pages, 12 figures

Flares are known to play an important role for the evolution of the atmospheres of young planets. In order to understand the evolution of planets, it is thus important to study the flare-activity of young stars. This is particularly the case for young M-stars, because they are very active. We study photometrically and spectroscopically the highly active M-star 2MASS J16111534-1757214. We show that it is a member of the Upper Sco OB association, which has an age of 5-10 Myrs. We also re-evaluate the status of other bona-fide M-stars in this region and identify 42 members. Analyzing the K2-light curves, we find that 2MASS J16111534-1757214 has, on average, one super-flare with E > 1.0E35 erg every 620 hours, and one with E >1.0E34 erg every 52 hours. Although this is the most active M-star in the Upper Sco association, the power-law index of its flare-distribution is similar to that of other M-stars in this region. 2MASS J16111534-1757214 as well as other M-stars in this region show a broken power-law distribution in the flare-frequency diagram. Flares larger than E >3E34 erg have a power-law index beta=-1.3+/-0.1 and flares smaller than that beta=-0.8+/-0.1. We furthermore conclude that the flare-energy distribution for young M-stars is not that different from solar-like stars.

Nina-Elisabet Nèmec (1, 2) ( (1) Max-Planck-Institut für Sonnensystemforschung, Göttingen, Germany, (2) Institut für Astrophysik, Georg-August-Universität Göttingen, Göttingen, Germany)

Doctoral thesis, 95 pages, 54 figures

The unprecedented precision of broadband stellar photometry achieved with the planet-hunting missions CoRoT and \textit{Kepler} initiated a new era in examining the magnetically-driven brightness variations of hundreds of thousands of stars. Such brightness variations are well studied and understood for the Sun. The plethora of data allows to accurately compare solar and stellar brightness variations. An intriguing question is whether the observed trends in the stellar photometric variability (e.g. the dependence of the variability on the stellar rotation period) can be explained by utilising the solar paradigm, in particular the physical concepts of brightness variations learnt from the Sun. The goal of this work is to find out, through comparison of observational and simulated data, if any physical concepts of solar brightness variability have to be altered to reproduce the distribution of Sun-like stars variabilities.

K. Balasubramaniam, L .Stawarz, C.C. Cheung, M. Sobolewska, V. Marchenko, R. Thimmappa, D.L. Krol, G. Migliori, A. Siemiginowska

Submitted to the Astrophysical Journal

The radio source 1146+596 is hosted by an elliptical/S0 galaxy NGC\,3894, with a low-luminosity active nucleus. The radio structure is compact, suggesting a very young age of the jets in the system. Recently, the source has been confirmed as a high-energy (HE, $>0.1$\,GeV) $\gamma$-ray emitter, in the most recent accumulation of the {\it Fermi} Large Area Telescope (LAT) data. Here we report on the analysis of the archival {\it Chandra} X-ray Observatory data for the central part of the galaxy, consisting of a single 40\,ksec-long exposure. We have found that the core spectrum is best fitted by a combination of an ionized thermal plasma with the temperature of $\simeq 0.8$\,keV, and a moderately absorbed power-law component (photon index $\Gamma \simeq 1.4\pm 0.4$, hydrogen column density $N_{\rm H}/10^{22}$\,cm$^{-2}$\,$\simeq 2.4\pm 0.7$). We have also detected the iron K$\alpha$ line at $\simeq 6.47\pm 0.07$\,keV, with a large equivalent width of EW\,$\simeq 1$\,keV. Based on the simulations of the {\it Chandra}'s Point Spread Function (PSF), we have concluded that, while the soft thermal component is extended on the scale of the galaxy host, the hard X-ray emission within the narrow photon energy range 6.0--7.0\,keV originates within the unresolved core (effectively the central $<5^{\prime\prime}\simeq 1.2$\,kpc radius). The line is therefore indicative of the X-ray reflection from a cold neutral gas in the central regions of NGC\,3894. We discuss the implications of our findings in the context of the X-ray Baldwin effect, and the overall energetic of the system. We note that NGC\,3894 becomes the first HE $\gamma$-ray source with the detected K$\alpha$ iron line.

Ralf Aurich, Thomas Buchert, Martin J. France, Frank Steiner

31 pages, 10 figures and 3 tables, submitted

In this work we investigate the standard deviation of the Cosmic Microwave Background (CMB) temperature gradient field as a signature for a multiply connected nature of the Universe. CMB simulations of a spatially infinite universe model within the paradigm of the standard cosmological model present non-zero two-point correlations at any angular scale. This is in contradiction with the extreme suppression of correlations at scales above $60^{\circ}$ in the observed CMB maps. Universe models with spatially multiply connected topology contain typically a discrete spectrum of the Laplacian with a specific wave-length cut-off and thus lead to a suppression of the correlations at large angular scales, as observed in the CMB (in general there can be also an additional continuous spectrum). Among the simplest examples are 3-dimensional tori which possess only a discrete spectrum. To date, the universe models with non-trivial topology such as the toroidal space are the only models that possess a two-point correlation function showing a similar behaviour as the one derived from the observed Planck CMB maps. In this work it is shown that the normalized standard deviation of the CMB temperature gradient field does hierarchically detect the change in size of the cubic 3-torus. It is also shown that the variance of the temperature gradient of the Planck maps is in slight anomaly with the median value of simulations within the standard cosmological model. All flat tori are globally homogeneous, but are globally anisotropic. However, this study also presents a test showing a level of homogeneity and isotropy of all the CMB map ensembles for the different torus sizes considered that are nearly at the same weak level of anisotropy revealed by the CMB in the standard cosmological model.

Soler, J.D., Beuther, H., Syed, J., Wang, Y., Henning, Th., Glover, S. C. O., Klessen, R. S., Sormani, M. C., Heyer, M., Smith, R. J., Urquhart, J. S., Yang, J., Su, Y., Zhou, X

Accepted for publication in Astronomy & Astrophysics (7 pages, 7 figures). arXiv admin note: text overlap with arXiv:2007.07285

We present a statistical study of the filamentary structure orientation in the CO emission observations obtained in the Milky Way Imaging Scroll Painting (MWISP) survey in the range $25.8\deg < l < 49.7\deg$, $|b| \leq 1.25\deg$, and $-100 < v_{\rm LSR} < 135$ km/s. We found that most of the filamentary structures in the $^{12}$CO and $^{13}$CO emission do not show a global preferential orientation either parallel or perpendicular to the Galactic plane. However, we found ranges in Galactic longitude and radial velocity where the $^{12}$CO and $^{13}$CO filamentary structures are parallel to the Galactic plane. These preferential orientations are different from those found for the HI emission. We consider this an indication that the molecular structures do not simply inherit these properties from parental atomic clouds. Instead, they are shaped by local physical conditions, such as stellar feedback, magnetic fields, and Galactic spiral shocks.

Pierre Burger, Oliver Friedrich, Joachim Harnois-Déraps, Peter Schneider

21 pages, 13 figures

Studying the statistical properties of the large-scale structure in the Universe with weak gravitational lensing is a prime goal of several current and forthcoming galaxy surveys. The power that weak lensing has to constrain cosmological parameters can be enhanced by considering statistics beyond second-order shear correlation functions or power spectra. One such higher-order probe that has proven successful in observational data is the density split statistics (DSS), in which one analyses the mean shear profiles around points that are classified according to their foreground galaxy density. In this paper, we generalise the most accurate DSS model to allow for a broad class of angular filter functions used for the classification of the different local density regions. This approach is motivated by earlier findings showing that an optimised filter can provide tighter constraints on model parameters compared to the standard top-hat case. We build on large deviation theory approaches and approximations thereof to model the matter density PDF, and on perturbative calculations of higher-order moments of the density field. The novel addition relies on the generalisation of these previously employed calculations to allow for general filter functions and is validated on several sets of numerical simulations. The revised model fits well the simulation measurements, with a residual systematic offset that is small compared to the statistical accuracy of current weak lensing surveys. The accuracy of the model is slightly lower for a compensated filter than for a non-negative filter function, and that it increases with the filter size. Using a Fisher matrix approach, we find constraints comparable to the commonly used two-point cosmic shear measures. Hence, our DSS model can be used in competitive analyses of current cosmic shear data, while it may need refinements for forthcoming lensing surveys.

Mohammad Ansari Fard, Zahra Baghkhani, Laya Ghodsi, Sina Taamoli, Farbod Hassani, Shant Baghram

17 pages, 11 figures, 8 tables. Comments are welcomed!

In non-linear scales, the matter density distribution is not Gaussian. Consequently, the widely used two-point correlation function is not adequate anymore to capture the matter density field's entire behaviour. Among all statistics beyond correlation functions, the spherical contact and nearest neighbour distribution function seem promising tools to probe matter distribution in a non-linear regime. In this work, we use halos from cosmological $N$-body simulations and galaxy groups from the volume-limited galaxy group catalogues to compare the spherical contact distribution function with the nearest neighbour distribution function. We also calculate the J-function for different samples. Moreover, we consider the redshift evolution and mass-scale dependence of statistics in the simulations and dependence on the magnitude of volume-limited samples in group catalogues. The shape of the spherical contact probability distribution function is nearly skew-normal, with skewness and kurtosis being approximately 0.5 and 3, respectively. On the other hand, the nearest neighbour probability distribution function is nearly log-normal, with logarithmic skewness and kurtosis being approximately 0.1 and 2.5, respectively. Accordingly, the spherical contact distribution function probes larger scales compared to the nearest neighbour distribution function, which is influenced by details of structures. We also find a linear relation between the first and second moment of the spherical contact probability distribution function in simulations, which could be used as a distinguishing probe of cosmological models.

Yang Bai, Joshua Berger, Mrunal Korwar, Nicholas Orlofsky

15 pages, 3 figures

A novel mechanism, "catalyzed baryogenesis," is proposed to explain the observed baryon asymmetry in our universe. In this mechanism, the motion of a ball-like catalyst provides the necessary out-of-equilibrium condition, its outer wall has CP-violating interactions with the Standard Model particles, and its interior has baryon number violating interactions. We use the electroweak-symmetric ball model as an example of such a catalyst. In this model, electroweak sphalerons inside the ball are active and convert baryons into leptons. The observed baryon number asymmetry can be produced for a light ball mass and a large ball radius. Due to direct detection constraints on relic balls, we consider a scenario in which the balls evaporate, leading to dark radiation at testable levels.

Maximilian Dax, Stephen R. Green, Jonathan Gair, Jakob H. Macke, Alessandra Buonanno, Bernhard Schölkopf

7+12 pages, 4+11 figures

We demonstrate unprecedented accuracy for rapid gravitational-wave parameter estimation with deep learning. Using neural networks as surrogates for Bayesian posterior distributions, we analyze eight gravitational-wave events from the first LIGO-Virgo Gravitational-Wave Transient Catalog and find very close quantitative agreement with standard inference codes, but with inference times reduced from O(day) to a minute per event. Our networks are trained using simulated data, including an estimate of the detector-noise characteristics near the event. This encodes the signal and noise models within millions of neural-network parameters, and enables inference for any observed data consistent with the training distribution, accounting for noise nonstationarity from event to event. Our algorithm -- called "DINGO" -- sets a new standard in fast-and-accurate inference of physical parameters of detected gravitational-wave events, which should enable real-time data analysis without sacrificing accuracy.

Gongjun Choi, Weikang Lin, Luca Visinelli, Tsutomu T. Yanagida

5 pages, 1 table, 1 figure

Taking the recently reported non-zero rotation angle of the cosmic microwave background (CMB) linear polarization $\beta=0.35\pm0.14{\rm\, deg}$ as the hint for a pseudo Nambu-Goldstone boson quintessence dark energy (DE), we study the electroweak (EW) axion quintessence DE model where the axion mass is generated by the EW instantons. We find that the observed value of $\beta$ implies a non-trivial $U(1)$ electromagnetic anomaly coefficient ($c_{\gamma}$), once the current constraint on the DE equation of state is also taken into account. With the aid of the hypothetical high energy structure of the model inspired by the experimentally inferred $c_{\gamma}$, the model is shown to be able to make prediction for the current equation of state ($w_{\rm DE,0}$) of the quintessence DE. This is expected to make our scenario distinguishable in comparison with the cosmological constant ($w=-1$) and testable in future when the error in the future measurement of $w_{\rm DE,0}$ is reduced to $\mathcal{O}(1)\%$ level ($\delta w=\mathcal{O}(10^{-2})$).

Davide Astesiano, Sergio L. Cacciatori, Federico Re

13 pages

We analyze the dynamics of a single spiral galaxy from a general relativistic viewpoint. We employ the known family of stationary axially-symmetric solutions to Einstein gravity coupled with dust in order to model the halo external to the bulge. In particular, we generalize the known results of Balasin and Grumiller, relaxing the condition of co-rotation, thus including non co-rotating dust. This further highlights the discrepancy between Newtonian theory of gravity and general relativity at low velocities and energy densities. We investigate the role of dragging in simulating dark matter effects. In particular, we show that non co-rotance further reduce the amount of energy density required to explain the rotation curves for spiral galaxies.

We present a short introduction to a non-standard cosmological scenario motivated by the duality symmetries of string theory, in which the big bang singularity is replaced with a "big bounce" at high but finite curvature. The bouncing epoch is prepared by a long (possibly infinitely extended) phase of cosmic evolution, starting from an initial state asymptotically approaching the string perturbative vacuum.

The spontaneous scalarization during the stellar core collapse in the massive scalar-tensor theories of gravity introduces extra polarizations (on top of the plus and cross modes) in gravitational waves, whose amplitudes are determined by several model parameters. Observations of such scalarization-induced gravitational waveforms therefore offer valuable probes into these theories of gravity. Considering a triple-scalar interactions in such theories, we find that the self-coupling effects suppress the magnitude of the scalarization and thus reduce the amplitude of the associated gravitational wave signals. In addition, the self-interacting effects in the gravitational waveform are shown to be negligible due to the dispersion throughout the astrophysically distant propagation. As a consequence, the gravitational waves observed on the Earth feature the characteristic inverse-chirp pattern. Although not with the on-going ground-based detectors, we illustrate that the scalarization-induced gravitational waves may be detectable at a signal-to-noise ratio level of ${\cal O}(100)$ with future detectors, such as Einstein Telescope and Cosmic Explorer.

The H + D_2^+(v=0,1 and 2) charge transfer reaction is studied using an accurate wave packet method, using recently proposed coupled diabatic potential energy surfaces. The state-to-state cross section is obtained for three different channels: non-reactive charge transfer, reactive charge transfer, and exchange reaction. The three processes proceed via the electronic transition from the first excited to the ground electronic state. The cross section for the three processes increases with the initial vibrational excitation. The non-reactive charge transfer process is the dominant channel, whose branching ratio increases with collision energy, and it compares well with experimental measurements at collision energies around 0.5 eV. For lower energies the experimental cross section is considerably higher, suggesting that it corresponds to higher vibrational excitation of D_2^+(v) reactants. Further experimental studies of this reaction and isotopic variants are needed, where conditions are controlled to obtain a better analysis of the vibrational effects of the D_2^+ reagents.