Abstract visibility
Abstract text size

Papers for Wednesday, Jul 12 2023

We conduct two-dimensional particle-in-cell simulations to investigate the scattering of electron heat flux by self-generated oblique electromagnetic waves. The heat flux is modeled as a bi-kappa distribution with a T_parallel > T_perp temperature anisotropy maintained by continuous injection at the boundaries. The anisotropic distribution excites oblique whistler waves and filamentary-like Weibel instabilities. Electron velocity distributions taken after the system has reached a steady state show that these in stabilities inhibit the heat flux and drive the total distributions towards isotropy. Electron trajectories in velocity space show a circular-like diffusion along constant energy surfaces in the wave frame. The key parameter controlling the scattering rate is the average speed, or drift speed vd, of the heat flux compared with the electron Alfven speed vAe, with higher drift speeds producing stronger fluctua tions and a more significant reduction of the heat flux. Reducing the density of the electrons carrying the heat flux by 50% does not significantly affect the scattering rate. A scaling law for the electron scattering rate versus vd/vAe is deduced from the simulations. The implications of these results for understanding energetic electron transport during solar flare energy release are discussed.

Stephanie M. Urbano Stawinski, M. C. Cooper, Steven L. Finkelstein, Intae Jung, Pablo G. Pérez-González, Caitlin M. Casey, Olivia R. Cooper, Nimish P. Hathi, Benne W. Holwerda, Anton M. Koekemoer, Vital Fernández, Rebecca L. Larson, Ray A. Lucas, L. Y. Aaron Yung

8 pages, 6 figures, 2 tables, submitted to MNRAS

We present a spectroscopic survey of Ly$\alpha$ emitters in the Extended Groth Strip (EGS) field, targeting the regime near the Epoch of Reionization. Using Keck/DEIMOS, we observed 947 high-$z$ candidates with photometric redshifts from 3 $< z_\text{phot} <$ 7 and down to an $H$-band (HST/WFC3 F160W) magnitude limit of < 27.5. Observations were taken over the course of 8 nights, with integration times ranging from 4 to 7.8 hours. Our survey secured 137 unique redshifts, 126 of which are Ly$\alpha$ emitters at 2.8 $< z <$ 6.5 with a mean redshift of $\overline{z} = 4.3$. We provide a comprehensive redshift catalog for our targets, as well as the reduced one- and two- dimensional spectra for each object. These observations will provide an important auxiliary dataset for the JWST Directors Discretionary Early Release Science (DD-ERS) program the Cosmic Evolution Early Release Science Survey (CEERS), which recently completed near- and mid-IR imaging and spectroscopy of galaxies in the EGS field.

David Berardo, Julien de Wit, Benjamin V. Rackham

10 pages, 3 figures

Transmission spectroscopy is currently the most powerful technique to study a wide range of planetary atmospheres, leveraging the filtering of a stars light by a planets atmosphere rather than its own emission. However, both a planet and its star contribute to the information encoded in a transmission spectrum and a particular challenge relate to disentangling their contributions. As measurements improve, the lack of fidelity of stellar spectra models present a bottleneck for accurate disentanglement. Considering JWST and future high-precision spectroscopy missions, we investigate the ability to derive empirical constraints on the emission spectra of stellar surface heterogeneities (i.e., spots and faculae) using the same facility as used to acquire the transmission spectra intended to characterize a given atmosphere. Using TRAPPIST-1 as a test case, we demonstrate that it is possible to constrain the photospheric spectrum to 0.2% and the spectra of stellar heterogeneities to within 1-5%, which will be valuable benchmarks to inform the new generation of theoretical stellar models. Long baseline of observations (90% of the stellar rotation period) are necessary to ensure the photon-limited (i.e., instrument-limited) exploration of exoplanetary atmospheres via transmission spectroscopy.

Katerina Slavicinska, Marina Gomes Rachid, Will Robson Monteiro Rocha, Ko-Ju Chuang, Ewine Fleur van Dishoeck, Harold Linnartz

Accepted for publication in A&A. 25 pages, 11 figures, 9 tables

This work aims at characterizing the mid-IR spectra of formamide in its pure form as well as in mixtures of the most abundant interstellar ices via laboratory simulation of such ices, as well as demonstrating how these laboratory spectra can be used to search for formamide in ice observations. Mid-IR spectra (4000 - 500 cm$^{-1}$, 2.5 - 20 $\mu$m) of formamide, both in its pure form as well as in binary and tertiary mixtures with H$_2$O, CO$_2$, CO, NH$_3$, CH$_3$OH, H$_2$O:CO$_2$, H$_2$O:NH$_3$, CO:NH$_3$, and CO:CH$_3$OH, are collected at temperatures ranging from 15 - 212 K. Apparent band strengths and positions of eight IR bands of pure amorphous and crystalline formamide at various temperatures are provided. Three bands are identified as potential formamide tracers in observational ice spectra: the overlapping C=O stretch and NH$_2$ scissor bands at 1700.3 and 1630.4 cm$^{-1}$ (5.881 and 6.133 $\mu$m), the CH bend at 1388.1 cm$^{-1}$ (7.204 $\mu$m), and the CN stretch at 1328.1 cm$^{-1}$ (7.529 $\mu$m). The relative apparent band strengths, positions, and FWHM of these features in mixtures at various temperatures are also determined. Finally, the laboratory spectra are compared to observational spectra of low- and high-mass young stellar objects as well as pre-stellar cores observed with the Infrared Space Observatory, the Spitzer Space Telescope, and the JWST. A comparison between the formamide CH bend in laboratory data and the 7.24 $\mu$m band in the observations tentatively indicates that, if formamide ice is contributing significantly to the observed absorption, it is more likely in a polar matrix. Upper limits ranging from 0.35-5.1\% with respect to H$_{2}$O are calculated. These upper limits are in agreement with gas-phase formamide abundances and take into account the effect of a H$_{2}$O matrix on formamide's band strengths.

Andris Dorozsmai, Silvia Toonen, Alejandro Vigna-Gómez, Selma E. de Mink, Floris Kummer

27 pages, 17 figures, 3 tables. Submitted to MNRAS

Observations suggest that massive stellar triples are common. However, their evolution is not yet fully understood. We investigate the evolution of hierarchical triples in which the stars of the inner binary experience chemically homogeneous evolution (CHE), particularly to understand the role of the tertiary star in the formation of gravitational-wave (GW) sources. We use the triple-star rapid population synthesis code TRES to determine the evolution of these systems at two representative metallicities: $Z = 0.005$ and $Z = 0.0005$. About half of all triples harbouring a CHE inner binary (CHE triples) experience tertiary mass transfer (TMT) episodes, an event which is rare for classically evolving stars. In the majority of TMT episodes, the inner binary consists of two main-sequence stars (58-60 per cent) or two black holes (BHs, 24-31 per cent). Additionally, we explore the role of von Zeipel-Lidov-Kozai (ZLK) oscillations for CHE triples. ZLK oscillations can result in eccentric stellar mergers or lead to the formation of eccentric compact binaries in systems with initial outer pericenters smaller than $\sim$ 1200 $R_{\odot}$. Approximately 24-30 per cent of CHE triples form GW sources, and in 31 per cent of these, the tertiary star plays a significant role and leads to configurations that are not predicted for isolated binaries. We conclude that the evolution of CHE binaries can be affected by a close tertiary companion, resulting in astronomical transients such as BH-BH binaries that merge via GW emission orders of magnitude faster than their isolated binary counterparts and tertiary-driven massive stellar mergers.

Junhan Kim, Jack Sayers, Mauro Sereno, Iacopo Bartalucci, Loris Chappuis, Sabrina De Grandi, Federico De Luca, Marco De Petris, Megan E. Donahue, Dominique Eckert, Stefano Ettori, Massimo Gaspari, Fabio Gastaldello, Raphael Gavazzi, Adriana Gavidia, Simona Ghizzardi, Asif Iqbal, Scott Kay, Lorenzo Lovisari, Ben J. Maughan, Pasquale Mazzotta, Nobuhiro Okabe, Etienne Pointecouteau, Gabriel W. Pratt, Mariachiara Rossetti, Keiichi Umetsu

submitted to A&A, comments welcome

Galaxy clusters are the products of structure formation through myriad physical processes that affect their growth and evolution throughout cosmic history. As a result, the matter distribution within galaxy clusters, or their shape, is influenced by cosmology and astrophysical processes, in particular the accretion of new material due to gravity. We introduce an analysis method to investigate the 3D triaxial shapes of galaxy clusters from the Cluster HEritage project with XMM-Newton -- Mass Assembly and Thermodynamics at the Endpoint of structure formation (CHEX-MATE). In this work, the first paper of a CHEX-MATE triaxial analysis series, we focus on utilizing X-ray data from XMM and Sunyaev-Zel'dovich (SZ) effect maps from Planck and ACT to obtain a three dimensional triaxial description of the intracluster medium (ICM) gas. We present the forward modeling formalism of our technique, which projects a triaxial ellipsoidal model for the gas density and pressure to compare directly with the observed two dimensional distributions in X-rays and the SZ effect. A Markov chain Monte Carlo is used to estimate the posterior distributions of the model parameters. Using mock X-ray and SZ observations of a smooth model, we demonstrate that the method can reliably recover the true parameter values. In addition, we apply the analysis to reconstruct the gas shape from the observed data of one CHEX-MATE galaxy cluster, Abell 1689, to illustrate the technique. The inferred parameters are in agreement with previous analyses for that cluster, and our results indicate that the geometrical properties, including the axial ratios of the ICM distribution, are constrained to within a few percent. With much better precision than previous studies, we thus further establish that Abell 1689 is significantly elongated along the line of sight, resulting in its exceptional gravitational lensing properties.

Huanqing Chen

4 pages, 4 figures, submitted to MNRAS Letter

Spectroscopic analysis of Ly$\alpha$ damping wings of bright sources at $z>6$ is a promising way to measure the reionization history of the universe. However, the theoretical interpretation of the damping wings is challenging due to the inhomogeneous nature of the reionization process and the proximity effect of bright sources. In this Letter, we analyze the damping wings arising from the neutral patches in the radiative transfer cosmological simulation suite Cosmic Reionization on Computers (CROC). We find that the damping wing profile remains a tight function of volume-weighted neutral fraction $\left< x_{\rm HI} \right>_{\rm v}$, especially when $\left< x_{\rm HI} \right>_{\rm v}>0.5$, despite the patchy nature of reionization and the proximity effect. This small scatter indicates that with a well-measured damping wing profile, we could constrain the volume-weighted neutral fraction as precise as $\Delta \left< x_{\rm HI} \right>_{\rm v} \lesssim 0.1$ in the first half of reionization.

Matthew A. Price, Matthijs Mars, Matthew M. Docherty, Alessio Spurio Mancini, Augustin Marignier, Jason. D. McEwen

Cosmic strings are linear topological defects that may have been produced during symmetry-breaking phase transitions in the very early Universe. In an expanding Universe the existence of causally separate regions prevents such symmetries from being broken uniformly, with a network of cosmic string inevitably forming as a result. To faithfully generate observables of such processes requires computationally expensive numerical simulations, which prohibits many types of analyses. We propose a technique to instead rapidly emulate observables, thus circumventing simulation. Emulation is a form of generative modelling, often built upon a machine learning backbone. End-to-end emulation often fails due to high dimensionality and insufficient training data. Consequently, it is common to instead emulate a latent representation from which observables may readily be synthesised. Wavelet phase harmonics are an excellent latent representations for cosmological fields, both as a summary statistic and for emulation, since they do not require training and are highly sensitive to non-Gaussian information. Leveraging wavelet phase harmonics as a latent representation, we develop techniques to emulate string induced CMB anisotropies over a 7.2 degree field of view, with sub-arcminute resolution, in under a minute on a single GPU. Beyond generating high fidelity emulations, we provide a technique to ensure these observables are distributed correctly, providing a more representative ensemble of samples. The statistics of our emulations are commensurate with those calculated on comprehensive Nambu-Goto simulations. Our findings indicate these fast emulation approaches may be suitable for wide use in, e.g., simulation based inference pipelines. We make our code available to the community so that researchers may rapidly emulate cosmic string induced CMB anisotropies for their own analysis.

Manuel Arca Sedda, Albrecht W. H. Kamlah, Rainer Spurzem, Mirek Giersz, Peter Berczik, Sara Rastello, Giuliano Iorio, Michela Mapelli, Massimiliano Gatto, Eva K. Grebel

22 pages, 21 figures, 4 tables. Comments welcome. Submitted to MNRAS

We present the first results of the \textsc{Dragon-II} simulations, a suite of 19 $N$-body simulations of star clusters with up to $10^6$ stars, with up to $33\%$ of them initially paired in binaries. In this work, we describe the main evolution of the clusters and their compact objects (COs). All \textsc{Dragon-II} clusters form in their centre a black hole (BH) subsystem with a density $10-100$ times larger than the stellar density, with the cluster core containing $50-80\%$ of the whole BH population. In all models, the BH average mass steeply decreases as a consequence of BH burning, reaching values $\langle m_{\rm BH}\rangle < 15$ M$_\odot$ within $10-30$ relaxation times. Generally, our clusters retain only BHs lighter than $30$ M$_\odot$ over $30$ relaxation times. Looser clusters retain a higher binary fraction, because in such environments binaries are less likely disrupted by dynamical encounters. We find that BH-main sequence star binaries have properties similar to recently observed systems. Double CO binaries (DCOBs) ejected from the cluster exhibit larger mass ratios and heavier primary masses than ejected binaries hosting a single CO (SCOBs). Ejected SCOBs have BH masses $m_{\rm BH} = 3-20$ M$_\odot$, definitely lower than those in DCOBs ($m_{\rm BH} = 10-100$ M$_\odot$).

Manuel Arca Sedda, Albrecht W. H. Kamlah, Rainer Spurzem, Francesco Paolo Rizzuto, Mirek Giersz, Thorsten Naab, Peter Berczik

15 pages, 6 figures, 2 tables, 1 appendix. Comments welcome. Submitted to MNRAS

The processes that govern the formation of intermediate-mass black holes (IMBHs) in dense stellar clusters are still unclear. Here, we discuss the role of stellar mergers, star-BH interactions and accretion, as well as BH binary (BBH) mergers in seeding and growing IMBHs in the \textsc{Dragon-II} simulation database, a suite of 19 direct $N$-body models representing dense clusters with up to $10^6$ stars. \textsc{Dragon-II} IMBHs have typical masses of $m_{\rm IMBH} = (100-380)$ M$_\odot$ and relatively large spins $\chi_{\rm IMBH} > 0.6$. We find a link between the IMBH formation mechanism and the cluster structure. In clusters denser than $3\times 10^5$ M$_\odot$ pc$^{-3}$, the collapse of massive star collision products represents the dominant IMBH formation process, leading to the formation of heavy IMBHs ($m_{\rm IMBH} > 200$ M$_\odot$), possibly slowly rotating, that form over times $<5$ Myr and grow further via stellar accretion and mergers in just $<30$ Myr. BBH mergers are the dominant IMBH formation channel in less dense clusters, for which we find that the looser the cluster, the longer the formation time ($10-300$ Myr) and the larger the IMBH mass, although remaining within $200$ M$_\odot$. Strong dynamical scatterings and relativistic recoil efficiently eject all IMBHs in \textsc{Dragon-II} clusters, suggesting that IMBHs in this type of cluster are unlikely to grow beyond a few $10^2$ M$_\odot$.

Manuel Arca Sedda, Albrecht W. H. Kamlah, Rainer Spurzem, Francesco Paolo Rizzuto, Mirek Giersz, Thorsten Naab, Peter Berczik

22 pages, 14 figures, 3 tables. Comments welcome. Submitted to MNRAS

Compact binary mergers forming in star clusters may exhibit distinctive features that can be used to identify them among observed gravitational-wave (GW) sources. Such features likely depend on the host cluster structure and the physics of massive star evolution. Here, we dissect the population of compact binary mergers in the \textsc{Dragon-II} simulation database, a suite of 19 direct $N$-body models representing dense star clusters with up to $10^6$ stars and $<33\%$ of stars in primordial binaries. We find a substantial population of black hole binary (BBH) mergers, some of them involving an intermediate-mass BH (IMBH), and a handful mergers involving a stellar BH and either a neutron star (NS) or a white dwarf (WD). Primordial binary mergers, $\sim 30\%$ of the whole population, dominate ejected mergers. Dynamical mergers, instead, dominate the population of in-cluster mergers and are systematically heavier than primordial ones. Around $20\%$ of \textsc{Dragon-II} mergers are eccentric in the LISA band and $5\%$ in the LIGO band. We infer a mean cosmic merger rate of $\mathcal{R}\sim 12(4.4)(1.2)$ yr$^{-1}$ Gpc$^3$ for BBHs, NS-BH, and WD-BH binary mergers, respectively, and discuss the prospects for multimessenger detection of WD-BH binaries with LISA. We model the rate of pair-instability supernovae (PISNe) in star clusters and find that surveys with a limiting magnitude $m_{\rm bol}=25$ can detect $\sim 1-15$ yr$^{-1}$ PISNe. Comparing these estimates with future observations could help to pin down the impact of massive star evolution on the mass spectrum of compact stellar objects in star clusters.

Vishnu Varma, Bernhard Mueller

Submitted to MNRAS (14 pages, 11 Figures)

We present a first 3D magnetohydrodynamic (MHD) simulation of oxygen, neon and carbon shell burning in a rapidly rotating 16 M_sun core-collapse supernova progenitor. We also run a purely hydrodynamic simulation for comparison. After 180s (15 and 7 convective turnovers respectively), the magnetic fields in the oxygen and neon shells achieve saturation at 10^{11}G and 5 x 10^{10}G. The strong Maxwell stresses become comparable to the radial Reynolds stresses and eventually suppress convection. The suppression of mixing by convection and shear instabilities results in the depletion of fuel at the base of the burning regions, so that the burning shell eventually move outward to cooler regions, thus reducing the energy generation rate. The strong magnetic fields efficiently transport angular momentum outwards, quickly spinning down the rapidly rotating convective oxygen and neon shells and forcing them into rigid rotation. The hydrodynamic model shows complicated redistribution of angular momentum and develops regions of retrograde rotation at the base of the convective shells. We discuss implications of our results for stellar evolution and for the subsequent core-collapse supernova. The rapid redistribution of angular momentum in the MHD model casts some doubt on the possibility of retaining significant core angular momentum for explosions driven by millisecond magnetars. However, findings from multi-D models remain tentative until stellar evolution calculations can provide more consistent rotation profiles and estimates of magnetic field strengths to initialise multi-D simulations without substantial numerical transients. We also stress the need for longer simulations, resolution studies, and an investigation of non-ideal effects.

Ricard Aguilera-Miret, Carlos Palenzuela, Federico Carrasco, Daniele Viganò

18 pages, 14 figures. arXiv admin note: text overlap with arXiv:2112.08413

We perform a long and accurate Large-Eddy Simulation of a binary neutron star merger, following the newly formed remnant up to 110 milliseconds. The combination of high-order schemes, high-resolution and the gradient subgrid-scale model allow us to have among the highest effective resolutions ever achieved. Our results show that, although the magnetic fields are strongly amplified by the Kelvin-Helmholtz instability, they are coherent only over very short spatial scales until t \gtrsim 30 ms. Around that time, magnetic winding becomes more efficient leading to a linear growth of the toroidal component and slowly ordering the field to more axisymmetric, large scales. The poloidal component only starts to grow at small scales at much later times t \gtrsim 90 ms, in a way compatible with the magneto-rotational instability. No strong large-scale poloidal field or jet is produced in the timescales spanned by our simulation, although there is an helicoidal structure gradually developing at late times. We highlight that soon after the merger the topology is always strongly dominated by toroidal structures, with a complex distribution in the meridional plane and highly turbulent perturbations. Thus, starting with strong purely dipolar fields before the merger is largely inconsistent with the outcomes of a realistic evolution. Finally, we confirm the universality of the evolved topology, even when starting with very different magnetic fields confined to the outermost layers of each neutron star.

Michael E. Brown, Bryan J. Butler

Planetary Science Journal, in press

We have used the Atacama Large Millimeter Array (ALMA) to measure precise absolute astrometric positions and detect the astrometric wobble of dwarf planet Orcus and its satellite Vanth over a complete orbit. We also place upper limits to the astrometric wobble induced by Dysnomia on dwarf planet Eris around its orbit. From the Vanth-Orcus barycentric motion, we find a Vanth-Orcus mass ratio of 0.16$\pm$0.02 -- the highest of any known planet or dwarf planet. This large ratio is consistent with the hypothesis that Vanth is a largely-intact impactor from a giant collision in the system, and that the system has likely evolved to a double synchronous state. We find only an upper limit of the barycenter motion of Eris, which implies a one sigma upper limit to the Dysnomia-Eris mass ratio of 0.0085, close to the modeled transition region between giant impact generated satellites which are largely intact remnants of the original impactor and those which form out of reaccreted disk material left over post-impact. The low albedo of Dysnomia leads us to marginally favor the intact impactor scenario. We find that Dysnomia has density of <1.2 g cm$^{-3}$, significantly lower than the 2.4 g cm$^{-3}$ of Eris.

Debby Tran, Benjamin Williams, Emily Levesque, Margaret Lazzarini, Julianne Dalcanton, Andrew Dolphin, Brad Koplitz, Adam Smercina, O. Grace Telford

19 pages, 11 figures, accepted for publication in ApJ

The nearby face-on star forming spiral galaxy NGC 6946 is known as the Fireworks Galaxy due to its hosting an unusually large number of supernova. We analyze its resolved near-ultraviolet (NUV) stellar photometry measured from images taken with the Hubble Space Telescope's (HST) Wide Field Camera 3 (WFC3) with F275W and F336W filters. We model the color-magnitude diagrams (CMD) of the UV photometry to derive the spatially-resolved star formation history (SFH) of NGC 6946 over the last 25 Myr. From this analysis, we produce maps of the spatial distribution of young stellar populations and measure the total recent star formation rate (SFR) of nearly the entire young stellar disk. We find the global SFR(age$\leq$25 Myr)=$13.17 \substack{+0.91 \\-0.79} M_\odot/\rm yr$. Over this period, the SFR is initially very high ($23.39\substack{+2.43\\-2.11} M_\odot/\rm yr$ between 16-25 Myr ago), then monotonically decreases to a recent SFR of $5.31\substack{+0.19\\-0.17} M_\odot/\rm yr$ in the last 10 Myr. This decrease in global star formation rate over the last 25 Myr is consistent with measurements made with other SFR indicators. We discuss in detail two of the most active regions of the galaxy, which we find are responsible for 3% and 5% of the total star formation over the past 6.3 Myr.

Branislav Rabatin, David C. Collins

Submitted to MNRAS

In star-forming clouds, high velocity flow gives rise to large fluctuations of density. In this work we explore the correlation between velocity magnitude (speed) and density. We develop an analytic formula for the joint probability distribution (PDF) of density and speed, and discuss its properties. In order to develop an accurate model for the joint PDF, we first develop improved models of the marginalized distributions of density and speed. We confront our results with a suite of 12 supersonic isothermal simulations with resolution of $1024^3$ cells in which the turbulence is driven by 3 different forcing modes (solenoidal, mixed and compressive) and 4 r.m.s. Mach numbers (1, 2, 4, 8). We show, that for transsonic turbulence, density and speed are correlated to a considerable degree and the simple assumption of independence fails to accurately describe their statistics. In the supersonic regime, the correlations tend to weaken with growing Mach number. Our new model of the joint and marginalized PDFs are a factor of 3 better than uncorrelated, and provides insight into this important process.

We examine the effect of supermassive black hole (SMBH) mass scaling relation choice on the inferred SMBH mass population since redshift $z \sim 3$. To make robust predictions for the gravitational wave background (GWB) we must have a solid understanding of the underlying SMBH demographics. Using the SDSS and 3D HST+CANDELS surveys for $0 < z < 3$ we evaluate the inferred SMBH masses from two SMBH-galaxy scaling relations: $\mathrm{M_{BH}}$-$\mathrm{M_{bulge}}$ and $\mathrm{M_{BH}}$-$\sigma$. Our SMBH mass functions come directly from stellar mass measurements for $\mathrm{M_{BH}}$-$\mathrm{M_{bulge}}$, and indirectly from stellar mass and galaxy radius measurements along with the galaxy mass fundamental plane for $\mathrm{M_{BH}}$-$\sigma$. We find that there is a substantial difference in predictions especially for $z > 1$, and this difference increases out to $z = 3$. In particular we find that using velocity dispersion predicts a greater number of SMBHs with masses greater than $10^9 \mathrm{M}_\odot$. The GWB that pulsar timing arrays find evidence for is higher in amplitude than expected from GWB predictions which rely on high redshift extrapolations of local SMBH mass-galaxy scaling relations. The difference in SMBH demographics resulting from different scaling relations may be the origin for the mismatch between the signal amplitude and predictions. Generally, our results suggest that a deeper understanding of the potential redshift evolution of these relations is needed if we are to draw significant insight from their predictions at $z > 1$

John Moustakas, Dustin Lang, Arjun Dey, Stéphanie Juneau, Aaron Meisner, Adam D. Myers, Edward F. Schlafly, David J. Schlegel, Francisco Valdes, Benjamin A. Weaver, Rongpu Zhou

47 pages, 13 figures, 5 tables; submitted to ApJS. All catalogs and data are publicly available through the SGA web-portal at this https URL

We present the 2020 version of the Siena Galaxy Atlas (SGA-2020), a multi-wavelength optical and infrared imaging atlas of 383,620 nearby galaxies. The SGA-2020 uses optical $grz$ imaging over $\approx20,000$ deg$^{2}$ from the DESI Legacy Imaging Surveys Data Release 9 and infrared imaging in four bands (spanning 3.4-22 $\mu$m) from the six-year unWISE coadds; it is more than 95% complete for galaxies larger than $R(26)\approx25$ arcsec and $r<18$ measured at the 26 mag arcsec$^{-2}$ isophote in the $r$-band. The atlas delivers precise coordinates, multi-wavelength mosaics, azimuthally averaged optical surface brightness profiles, model images and photometry, and additional ancillary metadata for the full sample. Coupled with existing and forthcoming optical spectroscopy from the Dark Energy Spectroscopic Instrument (DESI), the SGA-2020 will facilitate new detailed studies of the star formation and mass assembly histories of nearby galaxies; enable precise measurements of the local velocity field via the Tully-Fisher and Fundamental Plane relations; serve as a reference sample of lasting legacy value for time-domain and multi-messenger astronomical events; and more.

Jérémy Neveu, Vincent Brémaud, Pierre Antilogus, Florent Barret, Sébastien Bongard, Yannick Copin, Sylvie Dagoret-Campagne, Claire Juramy, Laurent Le-Guillou, Marc Moniez, Eduardo Sepulveda, The LSST Dark Energy Science Collaboration

30 pages, 36 figures, submitted to Astronomy and Astrophysics

In the next decade, many optical surveys will aim to tackle the question of dark energy nature, measuring its equation of state parameter at the permil level. This requires trusting the photometric calibration of the survey with a precision never reached so far, controlling many sources of systematic uncertainties. The measurement of the on-site atmospheric transmission for each exposure, or on average for each season or for the full survey, can help reach the permil precision for magnitudes. This work aims at proving the ability to use slitless spectroscopy for standard star spectrophotometry and its use to monitor on-site atmospheric transmission as needed, for example, by the Vera C. Rubin Observatory Legacy Survey of Space and Time supernova cosmology program. We fully deal with the case of a disperser in the filter wheel, which is the configuration chosen in the Rubin Auxiliary Telescope. The theoretical basis of slitless spectrophotometry is at the heart of our forward model approach to extract spectroscopic information from slitless data. We developed a publicly available software called Spectractor (https://github.com/LSSTDESC/Spectractor) that implements each ingredient of the model and finally performs a fit of a spectrogram model directly on image data to get the spectrum. We show on simulations that our model allows us to understand the structure of spectrophotometric exposures. We also demonstrate its use on real data, solving specific issues and illustrating how our procedure allows the improvement of the model describing the data. Finally, we discuss how this approach can be used to directly extract atmospheric transmission parameters from data and thus provide the base for on-site atmosphere monitoring. We show the efficiency of the procedure on simulations and test it on the limited data set available.

Weida Hu, Crystal L. Martin, Max Gronke, Simon Gazagnes, Matthew Hayes, John Chisholm, Timothy Heckman, Matilde Mingozzi, Namrata Roy, Peter Senchyna, Xinfeng Xu, Danielle A. Berg, Bethan L. James, Daniel P. Stark, Karla Z. Arellano-Córdova, Alaina Henry, Anne E. Jaskot, Nimisha Kumari, Kaelee S. Parker, Claudia Scarlata, Aida Wofford, Ricardo O. Amorín, Naunet Leonhardes-Barboza, Jarle Brinchmann, Cody Carr

40 pages, 19 figures, 5 tables, submitted to ApJ, comments welcome

Lyman-alpha line profiles are a powerful probe of ISM structure, outflow speed, and Lyman continuum escape fraction. In this paper, we present the Ly$\alpha$ line profiles of the COS Legacy Archive Spectroscopic SurveY, a sample rich in spectroscopic analogs of reionization-era galaxies. A large fraction of the spectra show a complex profile, consisting of a double-peaked Ly$\alpha$ emission profile in the bottom of a damped, Ly$\alpha$ absorption trough. Such profiles reveal an inhomogeneous interstellar medium (ISM). We successfully fit the damped Ly$\alpha$ absorption (DLA) and the Ly$\alpha$ emission profiles separately, but with complementary covering factors, a surprising result because this approach requires no Ly$\alpha$ exchange between high-$N_\mathrm{HI}$ and low-$N_\mathrm{HI}$ paths. The combined distribution of column densities is qualitatively similar to the bimodal distributions observed in numerical simulations. We find an inverse relation between Ly$\alpha$ peak separation and the [O III]/[O II] flux ratio, confirming that the covering fraction of Lyman-continuum-thin sightlines increases as the Ly$\alpha$ peak separation decreases. We combine measurements of Ly$\alpha$ peak separation and Ly$\alpha$ red peak asymmetry in a diagnostic diagram which identifies six Lyman continuum leakers in the CLASSY sample. We find a strong correlation between the Ly$\alpha$ trough velocity and the outflow velocity measured from interstellar absorption lines. We argue that greater vignetting of the blueshifted Ly$\alpha$ peak, relative to the redshifted peak, is the source of the well-known discrepancy between shell-model parameters and directly measured outflow properties. The CLASSY sample illustrates how scattering of Ly$\alpha$ photons outside the spectroscopic aperture reshapes Ly$\alpha$ profiles as the distances to these compact starbursts span a large range.

Cheongho Han, Youn Kil Jung, Ian A. Bond, Sun-Ju Chung, Michael D. Albrow, Andrew Gould, Kyu-Ha Hwang, Chung-Uk Lee, Yoon-Hyun Ryu, In-Gu Shin, Yossi Shvartzvald, Hongjing Yang, Jennifer C. Yee, Weicheng Zang, Sang-Mok Cha, Doeon Kim, Dong-Jin Kim, Seung-Lee Kim, Dong-Joo Lee, Yongseok Lee, Byeong-Gon Park, Richard W. Pogge, Fumio Abe, Richard Barry, David P. Bennett, Aparna Bhattacharya, Hirosame Fujii, Akihiko Fukui, Ryusei Hamada, Yuki Hirao, Stela Ishitani Silva, Yoshitaka Itow, Rintaro Kirikawa, Naoki Koshimoto, Yutaka Matsubara, Shota Miyazaki, Yasushi Muraki, Greg Olmschenk, Clément Ranc, Nicholas J. Rattenbury, Yuki Satoh, Takahiro Sumi, Daisuke Suzuki, Mio Tomoyoshi, Paul J. Tristram, Aikaterini Vandorou, Hibiki Yama, Kansuke Yamashita

11 pages, 10 tables, 8 figures

As a part of the project aiming to build a homogeneous sample of binary-lens (2L1S) events containing brown-dwarf (BD) companions, we investigate the 2021 season microlensing data collected by the Korea Microlensing Telescope Network (KMTNet) survey. For this purpose, we first identify 2L1S events by conducting systematic analyses of anomalous lensing events. We then select candidate BD-companion events by applying the criterion that the mass ratio between the lens components is less than $q_{\rm th}\sim 0.1$. From this procedure, we find four binary-lens events including KMT-2021-BLG-0588, KMT-2021-BLG-1110, KMT-2021-BLG-1643, and KMT-2021-BLG-1770, for which the estimated mass ratios are $q\sim 0.10$, 0.07, 0.08, and 0.15, respectively. The event KMT-2021-BLG-1770 is selected as a candidate despite the fact that the mass ratio is slightly greater than $q_{\rm th}$ because the lens mass expected from the measured short time scale of the event, $t_{\rm E}\sim 7.6$~days, is small. From the Bayesian analyses, we estimate that the primary and companion masses are $(M_1/M_\odot, M_2/M_\odot)= (0.54^{+0.31}_{-0.24}, 0.053^{+0.031}_{-0.023})$ for KMT-2021-BLG-0588L, $(0.74^{+0.27}_{-0.35}, 0.055^{+0.020}_{-0.026})$ for KMT-2021-BLG-1110L, $(0.73^{+0.24}_{-0.17}, 0.061^{+0.020}_{-0.014})$ for KMT-2021-BLG-1643L, and $(0.13^{+0.18}_{-0.07}, 0.020^{+0.028}_{-0.011})$ for KMT-2021-BLG-1770L. It is estimated that the probabilities of the lens companions being in the BD mass range are 82\%, 85\%, 91\%, and 59\% for the individual events. For confirming the BD nature of the lens companions found in this and previous works by directly imaging the lenses from future high-resolution adaptive-optics (AO) followup observations, we provide the lens-source separations expected in 2030, which is an approximate year of the first AO light on 30~m class telescopes.

N. Anugu, F. Baron, D. R. Gies, C. Lanthermann, G. H. Schaefer, K. A. Shepard, T. ten Brummelaar, J. D. Monnier, S. Kraus, J.-B. Le Bouquin, C. L. Davies, J. Ennis, T. Gardner, A. Labdon, R. M. Roettenbacher, B. R. Setterholm, W. Vollmann, C. Sigismondi

20 pages, accepted for AJ

The cool hypergiant star RW Cephei is currently in a deep photometric minimum that began several years ago. This event bears a strong similarity to the Great Dimming of the red supergiant Betelgeuse that occurred in 2019-2020. We present the first resolved images of RW Cephei that we obtained with the CHARA Array interferometer. The angular diameter and Gaia distance estimates indicate a stellar radius of 900 - 1760 R_sun which makes RW Cep one of the largest stars known in the Milky Way. The reconstructed, near-infrared images show a striking asymmetry in the disk illumination with a bright patch offset from center and a darker zone to the west. The imaging results depend on assumptions made about the extended flux, and we present two cases with and without allowing extended emission. We also present a recent near-infrared spectrum of RW Cep that demonstrates that the fading is much larger at visual wavelengths compared to that at near-infrared wavelengths as expected for extinction by dust. We suggest that the star's dimming is the result of a recent surface mass ejection event that created a dust cloud that now partially blocks the stellar photosphere.

Joost P. Wardenier, Vivien Parmentier, Michael R. Line, Elspeth K. H. Lee

21 pages, 14 figures, resubmitted to MNRAS after minor revisions

Ultra-hot Jupiters are perfect targets for transmission spectroscopy. However, their atmospheres feature strong spatial variations in temperature, chemistry, dynamics, cloud coverage, and scale height. This makes transit observations at high spectral resolution challenging to interpret. In this work, we model the cross-correlation signal of five chemical species (Fe, CO, H$_\text{2}$O, OH, and TiO) on WASP-76b, a benchmark ultra-hot Jupiter. We compute phase-dependent high-resolution transmission spectra of 3D SPARC/MITgcm models. The spectra are obtained with gCMCRT, a 3D Monte-Carlo radiative-transfer code. We find that, on top of atmospheric dynamics, the phase-dependent Doppler shift of the absorption lines in the planetary rest frame is shaped by the combined effect of planetary rotation and the unique 3D spatial distribution of chemical species. For species probing the dayside (e.g., refractories or molecules like CO and OH), the two effects act in tandem, leading to increasing blueshifts with orbital phase. For species that are depleted on the dayside (e.g., H$_\text{2}$O and TiO), the two effects act in an opposite manner, and could lead to increasing redshifts during the transit. This behaviour yields species-dependent offsets from a planet's expected $K_\text{p}$ value that can be much larger than planetary wind speeds. The offsets are usually negative for refractory species. We provide an analytical formula to estimate the size of a planet's $K_\text{p}$ offsets, which can serve as a prior for atmospheric retrievals. We conclude that observing the phase-resolved absorption signal of multiple species is key to constraining the 3D thermochemical structure and dynamics of ultra-hot Jupiters.

Alice Desmons, Sarah Brough, Francois Lanusse

Accepted at the ICML 2023 Workshop on Machine Learning for Astrophysics

Low surface brightness substructures around galaxies, known as tidal features, are a valuable tool in the detection of past or ongoing galaxy mergers. Their properties can answer questions about the progenitor galaxies involved in the interactions. This paper presents promising results from a self-supervised machine learning model, trained on data from the Ultradeep layer of the Hyper Suprime-Cam Subaru Strategic Program optical imaging survey, designed to automate the detection of tidal features. We find that self-supervised models are capable of detecting tidal features and that our model outperforms previous automated tidal feature detection methods, including a fully supervised model. The previous state of the art method achieved 76% completeness for 22% contamination, while our model achieves considerably higher (96%) completeness for the same level of contamination.

Ruixuan E. Chen, Jonathan H. Jiang, Philip E. Rosen, Kristen A. Fahy, Yanbei Chen

As the search for exoplanets continues, more are being discovered orbiting Red Giant stars. We use current data from the NASA Exoplanet Archive to investigate planet distribution around Red Giant stars and their presence in the host's habitable zone. As well, we update the power law relation between planet mass and stellar radius found in previous studies and provide more detailed investigations on this topic. Ten Red Giant-hosted exoplanets are found to be in the optimistically calculated habitable zone, five of which are in a more conservatively calculated habitable zone. We believe additional exoplanets can be found in habitable zones around Red Giants using the direct imaging and other methods, along with more powerful detection instrumentation.

Kangning Diao, Yi Mao

Accepted at the ICML 2023 Workshop on Machine Learning for Astrophysics

Large-scale numerical simulations ($\gtrsim 500\rm{Mpc}$) of cosmic reionization are required to match the large survey volume of the upcoming Square Kilometre Array (SKA). We present a multi-fidelity emulation technique for generating large-scale lightcone images of cosmic reionization. We first train generative adversarial networks (GAN) on small-scale simulations and transfer that knowledge to large-scale simulations with hundreds of training images. Our method achieves high accuracy in generating lightcone images, as measured by various statistics with mostly percentage errors. This approach saves computational resources by 90% compared to conventional training methods. Our technique enables efficient and accurate emulation of large-scale images of the Universe.

In this paper, we discuss the inflationary magnetogenesis scenario, in which the coupling function is introduced to break the conformal invariance of electromagnetic action. Unlike in conventional models, we deduce the Maxwell's equations under the perturbed FRW metric. We found that, the self-consistency of the action depends on the form of the coupling function when the scalar mode perturbations have been considered. Therefore, this self-consistency can be seen as a restriction on the coupling function. In this paper, we give the restrictive equation for coupling function then obtain the specific form of the coupling function in a simple model. We found that the coupling function depends on the potential of the inflaton and thus is model dependent. We obtain the power spectrum of electric field and magnetic field in large-field inflation model. We also found that the coupling function is a incresing function of time during slow-roll era as most of inflationary magnetogenesis models, it will lead to strong coupling problem. This issue is discussed qualitatively by introducing a correction function during the preheating.

Anna C. Childs, Cody Shakespeare, David R. Rice, Chao-Chin Yang, Jason H. Steffen

15 pages, 8 figures, 4 tables, and an appendix. Accepted for publication in MNRAS

We study the formation of the TRAPPIST-1 (T1) planets starting shortly after Moon-sized bodies form just exterior to the ice line. Our model includes mass growth from pebble accretion and mergers, fragmentation, type-I migration, and eccentricity and inclination dampening from gas drag. We follow the composition evolution of the planets fed by a dust condensation code that tracks how various dust species condense out of the disc as it cools. We use the final planet compositions to calculate the resulting radii of the planets using a new planet interior structure code and explore various interior structure models. Our model reproduces the broader architecture of the T1 system and constrains the initial water mass fraction of the early embryos and the final relative abundances of the major refractory elements. We find that the inner two planets likely experienced giant impacts and fragments from collisions between planetary embryos often seed the small planets that subsequently grow through pebble accretion. Using our composition constraints we find solutions for a two-layer model, a planet comprised of only a core and mantle, that match observed bulk densities for the two inner planets b and c. This, along with the high number of giant impacts the inner planets experienced, is consistent with recent observations that these planets are likely dessicated. However, two-layer models seem unlikely for most of the remaining outer planets which suggests that these planets have a primordial hydrosphere. Our composition constraints also indicate that no planets are consistent with a core-free interior structure.

Diogo Bragança, Yaniv Donath, Leonardo Senatore, Henry Zheng

38+14 pages, 12 figures, 8 tables. arXiv admin note: text overlap with arXiv:2111.05739 by other authors

After the successful full-shape analyses of BOSS data using the Effective Field Theory of Large-Scale Structure, we investigate what upcoming galaxy surveys might achieve. We introduce a ``perturbativity prior" that ensures that loop terms are as large as theoretically expected, which is effective in the case of a large number of EFT parameters. After validating our technique by comparison with already-performed analyses of BOSS data, we provide Fisher forecasts using the one-loop prediction for power spectrum and bispectrum for two benchmark surveys: DESI and MegaMapper. We find overall great improvements on the cosmological parameters. In particular, we find that MegaMapper (DESI) should obtain at least a 12$\sigma$ ($2\sigma$) evidence for non-vanishing neutrino masses, bound the curvature $\Omega_k$ to 0.0012 (0.012), and primordial inflationary non-Gaussianities as follows: $f_{\text{NL}}^{\text{loc.}}$ to $\pm 0.26$ (3.3), $f_{\text{NL}}^{\text{eq.}}$ to $\pm16$ (92), $f_{\text{NL}}^{\text{orth.}}$ to $\pm 4.2$ (27). Such measurements would provide much insight on the theory of Inflation. We investigate the limiting factor of shot noise and ignorance of the EFT parameters.

Suk Yee Yong, Cheng Soon Ong

Accepted for publication in MNRAS. 15 pages, 11 figures, 2 tables

Precise measurements of the black hole mass are essential to gain insight on the black hole and host galaxy co-evolution. A direct measure of the black hole mass is often restricted to nearest galaxies and instead, an indirect method using the single-epoch virial black hole mass estimation is used for objects at high redshifts. However, this method is subjected to biases and uncertainties as it is reliant on the scaling relation from a small sample of local active galactic nuclei. In this study, we propose the application of conformalised quantile regression (CQR) to quantify the uncertainties of the black hole predictions in a machine learning setting. We compare CQR with various prediction interval techniques and demonstrated that CQR can provide a more useful prediction interval indicator. In contrast to baseline approaches for prediction interval estimation, we show that the CQR method provides prediction intervals that adjust to the black hole mass and its related properties. That is it yields a tighter constraint on the prediction interval (hence more certain) for a larger black hole mass, and accordingly, bright and broad spectral line width source. Using a combination of neural network model and CQR framework, the recovered virial black hole mass predictions and uncertainties are comparable to those measured from the Sloan Digital Sky Survey. The code is publicly available at https://github.com/yongsukyee/uncertain_blackholemass.

Ce Sui, Xiaosheng Zhao, Tao Jing, Yi Mao

Accepted at the ICML 2023 Workshop on Machine Learning for Astrophysics, comments welcome

The ability to compress observational data and accurately estimate physical parameters relies heavily on informative summary statistics. In this paper, we introduce the use of mutual information (MI) as a means of evaluating the quality of summary statistics in inference tasks. MI can assess the sufficiency of summaries, and provide a quantitative basis for comparison. We propose to estimate MI using the Barber-Agakov lower bound and normalizing flow based variational distributions. To demonstrate the effectiveness of our method, we compare three different summary statistics (namely the power spectrum, bispectrum, and scattering transform) in the context of inferring reionization parameters from mock images of 21~cm observations with Square Kilometre Array. We find that this approach is able to correctly assess the informativeness of different summary statistics and allows us to select the optimal set of statistics for inference tasks.

Yue Huang, Dongli Shi, Xiaolu Zhang, Xiang Ma, Peng Zhang, Shijie Zheng, Liming Song, Xiaoyun Zhao, Wei Chen, Rui Qiao, Xinying Song, Jin Wang, Ce Cai, Shuo Xiao, Yanqiu Zhang, Shaolin Xiong

17 pages, 10 figures; Accepted for publication in RAA

Gravitational Wave High-energy Electromagnetic Counterpart All-sky Monitor (GECAM), consisting of two micro-satellites, is designed to detect gamma-ray bursts associated with gravitational-wave events. Here, we introduce the real-time burst alert system of GECAM, with the adoption of the BeiDou-3 short message communication service. We present the post-trigger operations, the detailed ground-based analysis, and the performance of the system. In the first year of the in-flight operation, GECAM was triggered by 42 GRBs. GECAM real-time burst alert system has the ability to distribute the alert within $\sim$1 minute after being triggered, which enables timely follow-up observations.

Alex Goater, Justin I. Read, Noelia E. D. Noël, Matthew D. A. Orkney, Stacy Y. Kim, Martin P. Rey, Eric P. Andersson, Oscar Agertz, Andrew Pontzen, Roberta Vieliute, Dhairya Kataria, Kiah Jeneway

10 pages, 4 figures; submitted to MNRAS

Ultra-faint dwarf galaxies (UFDs) are commonly found in close proximity to the Milky Way and other massive spiral galaxies. As such, their projected stellar ellipticity and extended light distributions are often thought to owe to tidal forces. In this paper, we study the projected stellar ellipticities and faint stellar outskirts of tidally isolated ultra-faints drawn from the 'Engineering Dwarfs at Galaxy Formation's Edge' (EDGE) cosmological simulation suite. Despite their tidal isolation, our simulated dwarfs exhibit a wide range of projected ellipticities ($0.03 < \varepsilon < 0.85$), with many possessing anisotropic extended stellar haloes that mimic tidal tails, but owe instead to late-time accretion of lower mass companions. Furthermore, we find a strong causal relationship between ellipticity and formation time of an UFD, which is robust to a wide variation in the feedback model. We show that the distribution of projected ellipticities in our suite of simulated EDGE dwarfs matches well with that of 21 Local Group dwarf galaxies. Given the ellipticity in EDGE arises from an ex-situ accretion origin, the agreement in shape indicates the ellipticities of some observed dwarfs may also originate from a similar non-tidal scenario. The orbital parameters of these observed dwarfs further support that they are not currently tidally disrupting. If the baryonic content in these galaxies is still tidally intact, then the same may be true for their dark matter content, making these galaxies in our Local Group pristine laboratories for testing dark matter and galaxy formation models.

Shunhao Ji, Zhongxiang Wang (1), Yi Xing (2), Dahai Yan (1), Jintao Zheng (1) (1. Yunnan University, 2. Shanghai Astronomical Observatory)

9 pages, 6 figures, 2 tables, submitted to ApJ

Intrigued by recent high-energy study results for nearby galaxies with gamma-ray emission and in particular NGC~1068 that has been detected as a neutrino-emitting source by the IceCube Neutrino Observatory, we conduct detailed analysis of the $\gamma$-ray data for the galaxies NGC~1068 and NGC~253, obtained with the Large Area Telescope onboard {\it the Fermi Gamma-ray Space Telescope}. By checking for their possible spectral features and then constructing light curves in corresponding energy ranges, we identify flare-like activity from NGC ~1068 in $\geq$2\,GeV energy range and significant long-term variations of NGC~253 in $\geq$5\,GeV energy range. In the former, the emission appears harder in the two half-year flare-like events than that in the otherwise `quiescent' state. In the latter, there is a 2-times decrease in the flux before and after MJD~57023, which is clearly revealed by the test-statistic maps we obtain. Considering studies carried out and models proposed for the $\gamma$-ray emissions of the two sources, we discuss the implications of our findings. The jet in NGC~1068 may contribute to the \gr\ emission. The nature of the long-term variations in NGC~253 is not clear, but the variation part of the emission may be connected to the very-high-energy (VHE) emission of the galaxy and could be verified by VHE observations.

Santosh Kumar Yadav, Anil Kumar Yadav

10 pages, 4 Figure panels. arXiv admin note: text overlap with arXiv:1901.07549

In this paper, we investigate an extension of standard $\Lambda$CDM model by allowing: a temporal evolution in the equation of state (EoS) of DM via Chevallier-Polarski-Linder parametrization, and the constant non-null sound speed. We also consider the properties of neutrinos, such as the effective neutrino mass and the effective number of neutrino species as free parameters. We derive the constraints on this scenario by using the data from the Planck-2018 cosmic microwave background (CMB), baryonic acoustic oscillations (BAO), the local value of the Hubble constant from the Hubble Space Telescope (HST), and some large scale structure (LSS) information from the abundance of galaxy clusters. We find constraints on the EoS and sound speed of DM very close to the null value in all cases thus no significant evidence can be stated beyond the standard CDM paradigm, and we conclude that the present observational data favor DM as a pressureless fluid. In all cases, we find the tighter upper bounds on the sum of neutrino masses. The most tight upper bound is $\sum m_\nu < 0.17$ eV at 95\% CL, imposed with the addition of HST prior in the analysis. We also observe the effects of neutrino properties on the extended DM parameters and also on other parameters. We find significantly lower mean values of $\sigma_8$ in all cases which are in agreement with the LSS measurements. Thus, the well-known $\sigma_8$ tension is reconciled in the considered model.

D. Depaoli, A. Chiavassa, D. Corti, F. Di Pierro, M. Mariotti, R. Rando

Accepted for publication in Nuclear Instruments and Methods in Physics Research A

The Cherenkov Telescope Array (CTA) will be the next generation ground-based gamma-ray observatory. CTA consists of different telescope types of which the largest ones (Large-Sized Telescopes, LSTs) cover the lower energy range, between 20 GeV and 200 GeV. The first LST is currently being commissioned at the Roque de los Muchachos Observatory, La Palma, Canary Islands. Its camera has 1855 photomultipliers (PMTs) with 1.5 inch cathodes. Silicon Photomultipliers (SiPMs) are increasingly becoming valid alternatives to PMTs also in gamma-ray astronomy. In the context of the LST project, there is an effort to study a novel Advanced Camera, equipped with SiPMs and a completely redesigned electronics based on a fully digital approach. To study and develop solutions on the sensors of these camera, we built a prototype camera module with a fully re-designed pre-amplifying stage and sensor bias control while re-using the digitizing and triggering stages of the existing LST camera module. We report on the design choices made to achieve the highest performance in terms of timing and charge resolution and the laboratory measurements validating those choices.

J. Orell-Miquel, M. Lampón, M. López-Puertas, M. Mallorquín, F. Murgas, A. Peláez-Torres, E. Pallé, E. Esparza-Borge, J. Sanz-Forcada, H. M. Tabernero, L. Nortmann, E. Nagel, H. Parviainen, M. R. Zapatero Osorio, J. A. Caballero, S. Czesla, C. Cifuentes, G. Morello, A. Quirrenbach, P. J. Amado, A. Fernández-Martín, A. Fukui, Th. Henning, K. Kawauchi, J. P. de Leon, K. Molaverdikhani, D. Montes, N. Narita, A. Reiners, I. Ribas, A. Sánchez-López, A. Schweitzer, M. Stangret, F. Yan

Accepted for publication in A&A. 17 pages, 18 figures

HD235088 (TOI-1430) is a young star known to host a sub-Neptune-sized planet candidate. We validated the planetary nature of HD235088 b with multiband photometry, refined its planetary parameters, and obtained a new age estimate of the host star, placing it at 600-800 Myr. Previous spectroscopic observations of a single transit detected an excess absorption of He I coincident in time with the planet candidate transit. Here, we confirm the presence of He I in the atmosphere of HD235088 b with one transit observed with CARMENES. We also detected hints of variability in the strength of the helium signal, with an absorption of $-$0.91$\pm$0.11%, which is slightly deeper (2$\sigma$) than the previous measurement. Furthermore, we simulated the He I signal with a spherically symmetric 1D hydrodynamic model, finding that the upper atmosphere of HD235088 b escapes hydrodynamically with a significant mass loss rate of (1.5-5) $\times$10$^{10}$g s$^{-1}$, in a relatively cold outflow, with $T$=3125$\pm$375 K, in the photon-limited escape regime. HD235088 b ($R_{p}$ = 2.045$\pm$0.075 R$_{\oplus}$) is the smallest planet found to date with a solid atmospheric detection - not just of He I but any other atom or molecule. This positions it a benchmark planet for further analyses of evolving young sub-Neptune atmospheres.

Can Xu, Tao Wang, Qiusheng Gu, Anita Zanella, Ke Xu, Hanwen Sun, Veronica Strazzullo, Francesco Valentino, Raphael Gobat, Emanuele Daddi, David Elbaz, Mengyuan Xiao, Shiying Lu, Luwenjia Zhou

13 pages, 10 figures, 1 table

Structural properties of cluster galaxies during their peak formation epoch, $z \sim 2-4$ provide key information on whether and how environment affects galaxy formation and evolution. Based on deep HST/WFC3 imaging towards the z=2.51 cluster, J1001, we explore environmental effects on the structure, color gradients, and stellar populations of a statistical sample of cluster SFGs. We find that the cluster SFGs are on average smaller than their field counterparts. This difference is most pronounced at the high-mass end ($M_{\star} > 10^{10.5} M_{\odot}$) with nearly all of them lying below the mass-size relation of field galaxies. The high-mass cluster SFGs are also generally old with a steep negative color gradient, indicating an early formation time likely associated with strong dissipative collapse. For low-mass cluster SFGs, we unveil a population of compact galaxies with steep positive color gradients that are not seen in the field. This suggests that the low-mass compact cluster SFGs may have already experienced strong environmental effects, e.g., tidal/ram pressure stripping, in this young cluster. These results provide evidence on the environmental effects at work in the earliest formed clusters with different roles in the formation of low and high-mass galaxies.

Sumedha Biswas, Zuzanna Kostrzewa-Rutkowska, Peter G. Jonker, Paul Vreeswijk, Deepak Eappachen, Paul J. Groot, Simon Hodgkin, Abdullah Yoldas, Guy Rixon, Diana Harrison, M. van Leeuwen, Dafydd Evans

13 pages, 9 figures; Submitted to MNRAS

The discovery of gravitational wave (GW) events and the detection of electromagnetic counterparts from GW170817 has started the era of multimessenger GW astronomy.The field has been developing rapidly and in this paper,we discuss the preparation for detecting these events with the ESA Gaia satellite,during the 4th observing run of the LIGO-Virgo-KAGRA (LVK) collaboration that has started on May 24,2023. Gaia is contributing to the search for GW counterparts by a new transient detection pipeline called GaiaX. In GaiaX, a new source appearing in the field of view of only one of the two telescopes on-board Gaia is sufficient to send out an alert on the possible detection of a new transient. Ahead of O4, an experiment was conducted over a period of about two months. During the two weeks around New Moon in this period of time, the MeerLICHT (ML) telescope located in South Africa tried (weather permitting) to observe the same region of the sky as Gaia within 10 minutes. Any GaiaX detected transient was published publicly. ML and Gaia have similar limiting magnitudes for typical seeing conditions at ML. At the end of the experiment, we had 11861 GaiaX candidate transients and 15806 ML candidate transients, which we further analysed and the results of which are presented in this paper. Finally, we discuss the possibility and capabilities of Gaia contributing to the search for electromagnetic counterparts of gravitational wave events during O4 through the GaiaX detection and alert procedure.

We investigate the migration of Mars- to super-Earth-sized planets in the vicinity of a pressure bump in a 3D radiative protoplanetary disc while accounting for the effect of accretion heat release. Pressure bumps have often been assumed to act as efficient migration traps, but we show that the situation changes when the thermal forces are taken into account. Our simulations reveal that for planetary masses $\lesssim$$2\,M_{\oplus}$, once their luminosity exceeds the critical value predicted by linear theory, thermal driving causes their orbits to become eccentric, quenching the positive corotation torque responsible for the migration trap. As a result, planets continue migrating inwards past the pressure bump. Additionally, we find that planets that remain circular and evolve in the super-Keplerian region of the bump exhibit a reversed asymmetry of their thermal lobes, with the heating torque having an opposite (negative) sign compared to the standard circular case, thus leading to inward migration as well. We also demonstrate that the super-critical luminosities of planets in question can be reached through the accretion of pebbles accumulating in the bump. Our findings have implications for planet formation scenarios that rely on the existence of migration traps at pressure bumps, as the bumps may repeatedly spawn inward-migrating low-mass embryos rather than harbouring newborn planets until they become massive.

Daria Zakharova, Benedetta Vulcani, Gabriella De Lucia, Lizhi Xie, Michaela Hirschmann, Fabio Fontanot

resubmitted to MNRAS after minor comments

Filaments are elongated structures that connect groups and clusters of galaxies and are visually the striking feature in cosmological maps. In the literature, typically filaments are defined only using galaxies, assuming that these are good tracers of the dark matter distribution, despite the fact that galaxies are a biased indicator. Here we apply the topological filament extractor DisPerSE to the predictions of the semi-analytic code GAEA to investigate the correspondence between the properties of $z=0$ filaments extracted using the distribution of dark matter and the distribution of model galaxies evolving within the same large-scale structure. We focus on filaments around massive clusters with a mass comparable to Virgo and Coma, with the intent of investigating the influence of massive systems and their feeding filamentary structure on the physical properties of galaxies. We apply different methods to compare the properties of filaments based on the different tracers and study how the sample selection impacts the extraction. Overall, filaments extracted using different tracers agree, although they never coincide totally. We also find that the number of filaments ending up in the massive clusters identified using galaxies distribution is typically underestimated with respect to the corresponding dark matter filament extraction.

Hsiu-Hsien Lin, Paul Scholz, Cherry Ng, Ue-Li Pen, Mohit Bhardwaj, Pragya Chawla, Alice P. Curtin, Ketan R. Sand, Shriharsh P. Tendulkar, Bridget Andersen, Kevin Bandura, Tomas Cassanelli, Amanda M. Cook, Matt Dobbs, Fengqiu Adam Dong, Gwendolyn Eadie, Emmanuel Fonseca, Bryan M. Gaensler, Utkarsh Giri, Antonio Herrera-Martin, Jane Kaczmarek, Joseph Kania, Victoria Kaspi, Kholoud Khairy, Adam E. Lanman, Calvin Leung, Dongzi Li, Kiyoshi W. Masui, Juan Mena-Parra, Bradley W. Meyers, Daniele Michilli, Nikola Milutinovic, Aaron B. Pearlman, Ziggy Pleunis, Masoud Rafiei-Ravandi, Mubdi Rahman, Pranav Sanghavi, Kaitlyn Shin, Kendrick Smith, Ingrid Stairs, David C. Stenning, Keith Vanderlinde, Dallas Wulf

20 pages, 16 figures. Submitted to AAS Journals

We report ten fast radio bursts (FRBs) detected in the far side-lobe region (i.e., $\geq 5^\circ$ off-meridian) of the Canadian Hydrogen Intensity Mapping Experiment (CHIME) from 2018 August 28 to 2021 August 31. We localize the bursts by fitting their spectra with a model of the CHIME/FRB synthesized beam response. CHIME/FRB did not observe repetition of similar brightness from the uniform sample of 10 side-lobe FRBs in a total exposure time of 35580 hours. Under the assumption of Poisson-distributed bursts, we infer that the mean repetition interval above the detecting threshold of the far side-lobe events is longer than 11880 hours, which is at least 2380 times larger than the interval from known CHIME/FRB detected repeating sources, with some caveats, notably that very narrow-band events could have been missed. Our results from these far side-lobe events suggest one of two scenarios: either (1) all FRBs repeat and the repetition intervals span a wide range, with high-rate repeaters being a rare subpopulation, or (2) non-repeating FRBs are a distinct population different from known repeaters.

Hsiu-Hsien Lin, Paul Scholz, Cherry Ng, Ue-Li Pen, D. Z. Li, Laura Newburgh, Alex Reda, Bridget Andersen, Kevin Bandura, Mohit Bhardwaj, Charanjot Brar, Tomas Cassanelli, Pragya Chawla, Amanda M. Cook, Alice P. Curtin, Matt Dobbs, Fengqiu Adam Dong, Emmanuel Fonseca, Bryan M. Gaensler, Utkarsh Giri, Alex S. Hill, Jane Kaczmarek, Joseph Kania, Victoria Kaspi, Kholoud Khairy, Calvin Leung, Kiyoshi W. Masui, Juan Mena-Parra, Bradley W. Meyers, Anna Ordog, Aaron B. Pearlman, Ziggy Pleunis, Masoud Rafiei-Ravandi, Mubdi Rahman, Scott Ransom, Ketan R. Sand, Pranav Sanghavi, Kaitlyn Shin, Kendrick Smith, Ingrid Stairs, Shriharsh P. Tendulkar, Keith Vanderlinde, Dallas Wulf

16 pages, 8 figures. Submitted to AAS Journals

We study the 10 fast radio bursts (FRBs) detected in the far side-lobe region of the CHIME telescope from 2018 August 28 to 2021 August 31. We find that the far side-lobe events have on average $\sim$500 times greater fluxes than events detected in CHIME's main lobe. We show that the side-lobe sample is therefore statistically $\sim$20 times closer than the main-lobe sample. The median dispersion measure (DM) excess, after removing the Galactic disk component using the NE2001 for the free electron density distribution of the Milky Way, of the 10 far side-lobe and 471 non-repeating main-lobe FRBs in the first CHIME/FRB catalog is 183.0 and 433.9 pc\;cm$^{-3}$, respectively. By comparing the DM excesses of the two populations under reasonable assumptions, we statistically constrain that the local degenerate contributions (from the Milky Way halo and the host galaxy) and the intergalactic contribution to the excess DM of the 471 non-repeating main-lobe FRBs for the NE2001 model are 131.2$-$158.3 and 302.7$-$275.6 pc cm$^{-3}$, respectively, which corresponds to a median redshift for the main-lobe FRB sample of $\sim$0.3. These constraints are useful for population studies of FRBs, and in particular for constraining the location of the missing baryons.

Martin Houde, Fereshteh Rajabi, Gordon C. MacLeod, Sharmila Goedhart, Yoshihiro Tanabe, Stefanus P. van den Heever, Christopher M. Wyenberg, Yoshinori Yonekura

15 pages, 11 figures

We discuss the role that coherence phenomena can have on the intensity variability of spectral lines associated with maser radiation. We do so by introducing the fundamental cooperative radiation phenomenon of (Dicke's) superradiance and discuss its complementary nature to the maser action, as well as its role in the flaring behaviour of some maser sources. We will consider examples of observational diagnostics that can help discriminate between the two, and identify superradiance as the source of the latter. More precisely, we show how superradiance readily accounts for the different time-scales observed in the multi-wavelength monitoring of the periodic flaring in G9.62+0.20E.

Mariko Kato, Izumi Hachisu, Hideyuki Saio

4 pages, 4 figures, accepted in MNRAS (Letters)

We propose a helium nova model for the Large Magellanic Cloud (LMC) supersoft X-ray source (SSS) [HP99]159. This object has long been detected as a faint and persistent SSS for about 30 years, and recently been interpreted to be a source of steady helium-shell burning, because no hydrogen lines are observed. We find that the object can also be interpreted as in a decaying phase of a helium nova. The helium nova is slowly decaying toward the quiescent phase, during which the observed temperature, luminosity, and SSS lifetime ($\gtrsim 30$ years) are consistent with a massive white dwarf model of $\sim$ 1.2 $M_\odot$. If it is the case, this is the second discovery of a helium nova outburst after V445 Pup in our Galaxy and also the first identified helium nova in the LMC. We also discuss the nature of the companion helium star in relation to Type Ia supernova progenitors.

Qing-Zheng Li, Yang Huang, Xiao-Bo Dong

10 pages, 7 figures, 2 tables, accepted by RAA; full table can be accessed from this https URL

We present an updated catalog of 46,753 radial velocity (RV) standard stars selected from the APOGEE DR17. These stars cover the Northern and Southern Hemispheres almost evenly, with 62% being red giants and 38% being main-sequence stars. These RV standard stars are stable on a baseline longer than 200 days (54% longer than one year and 10% longer than five years) with a median stability better than 215 m s$^{-1}$. The average observation number of those stars are 5 and each observation is required to have spectral-to-noise-ratio (SNR) greater than 50 and RV measurement error smaller than 500 m s$^{-1}$. Based on the new APOGEE RV standard star catalog, we have checked the RV zero points (RVZPs) for current large-scale stellar spectroscopic surveys including RAVE, LAMOST, GALAH and Gaia. By carefully analysis, we estimate their mean RVZP to be $+0.149$ km s$^{-1}$, $+4.574$ km s$^{-1}$ (for LRS), $-0.031$ km s$^{-1}$ and $+0.014$ km s$^{-1}$, respectively, for the four surveys. In the RAVE, LAMOST (for MRS), GALAH and Gaia surveys, RVZP exhibits systematic trend with stellar parameters (mainly [Fe/H], $T_{\rm{eff}}$, log $g$, $G_{\rm{BP}}-G_{\rm{RP}}$ and $G_{\rm{RVS}}$). The corrections of those small but clear RVZPs are of vital importances for these massive spectroscopic surveys in various studies that require extremely high radial velocity accuracies.

Daniel C. M. Palumbo, George N. Wong, Andrew A. Chael, Michael D. Johnson

14 pages, 7 figures, Accepted to ApJL

Images of supermassive black hole accretion flows contain features of both curved spacetime and plasma structure. Inferring properties of the spacetime from images requires modeling the plasma properties, and vice versa. The Event Horizon Telescope Collaboration has imaged near-horizon millimeter emission from both Messier 87* (M87*) and Sagittarius A* (Sgr A*) with very-long-baseline interferometry (VLBI) and has found a preference for magnetically arrested disk (MAD) accretion in each case. MAD accretion enables spacetime measurements through future observations of the photon ring, the image feature composed of near-orbiting photons. The ordered fields and relatively weak Faraday rotation of MADs yield rotationally symmetric polarization when viewed at modest inclination. In this letter, we utilize this symmetry along with parallel transport symmetries to construct a gain-robust interferometric quantity that detects the transition between the weakly lensed accretion flow image and the strongly lensed photon ring. We predict a shift in polarimetric phases on long baselines and demonstrate that the photon rings in M87* and Sgr A* can be unambiguously detected {with sensitive, long-baseline measurements. For M87* we find that photon ring detection in snapshot observations requires $\sim1$ mJy sensitivity on $>15$ G$\lambda$ baselines at 230 GHz and above, which could be achieved with space-VLBI or higher-frequency ground-based VLBI. For Sgr A*, we find that interstellar scattering inhibits photon ring detectability at 230 GHz, but $\sim10$ mJy sensitivity on $>12$ G$\lambda$ baselines at 345 GHz is sufficient, which is accessible from the ground. For both sources, these sensitivity requirements may be relaxed by repeated observations and averaging.

Nicolas Poirier (1, 2), Victor Réville (3), Alexis P. Rouillard (3), Athanasios Kouloumvakos (4), Emeline Valette (3) ((1) Rosseland Centre for Solar Physics - University of Oslo, (2) Institute of Theoretical Astrophysics - University of Oslo, (3) Research Institute for Astrophysics and Planetology - University Toulouse III Paul Sabatier, (4) The Johns Hopkins University Applied Physics Laboratory)

23 pages, 14 figures, to appear in Astronomy & Astrophysics, associated movies available at this https URL

We analyse the signature and origin of transient structures embedded in the slow solar wind, and observed by the Wide-Field Imager for Parker Solar Probe (WISPR) during its first 10 passages close to the Sun. WISPR provides a new in-depth vision on these structures, which have long been speculated to be a remnant of the pinch-off magnetic reconnection occurring at the tip of helmet streamers. We pursue the previous modelling works of Reville (2020b, 2022) that simulate the dynamic release of quasi-periodic density structures into the slow wind through a tearing-induced magnetic reconnection at the tip of helmet streamers. Synthetic WISPR white-light (WL) images are produced using a newly developed advanced forward modelling algorithm, that includes an adaptive grid refinement to resolve the smallest transient structures in the simulations. We analyse the aspect and properties of the simulated WL signatures in several case studies, typical of solar minimum and near-maximum configurations. Quasi-periodic density structures associated with small-scale magnetic flux ropes are formed by tearing-induced magnetic reconnection at the heliospheric current sheet and within 3-7Rs. Their appearance in WL images is greatly affected by the shape of the streamer belt and the presence of pseudo-streamers. The simulations show periodicities on the ~90-180min, ~7-10hr and ~25-50hr timescales, which are compatible with WISPR and past observations. This work shows strong evidence for a tearing-induced magnetic reconnection contributing to the long-observed high variability of the slow solar wind.

J. Aschersleben, T. T. H. Arnesen, R. F. Peletier, M. Vecchi, C. Vlasakidis, M. H. F. Wilkinson

10 pages, 9 figures, submitted to Nuclear Instruments and Methods in Physics Research - section A

Imaging Atmospheric Cherenkov Telescopes (IACTs) detect very high-energy gamma rays from ground level by capturing the Cherenkov light of the induced particle showers. Convolutional neural networks (CNNs) can be trained on IACT camera images of such events to differentiate the signal from the background and to reconstruct the energy of the initial gamma ray. Pattern spectra provide a 2-dimensional histogram of the sizes and shapes of features comprising an image and they can be used as an input for a CNN to significantly reduce the computational power required to train it. In this work, we generate pattern spectra from simulated gamma-ray and proton images to train a CNN for signal-background separation and energy reconstruction for the Small-Sized Telescopes (SSTs) of the Cherenkov Telescope Array (CTA). A comparison of our results with a CNN directly trained on CTA images shows that the pattern spectra-based analysis is about a factor of three less computationally expensive but not able to compete with the performance of the CTA images-based analysis. Thus, we conclude that the CTA images must be comprised of additional information not represented by the pattern spectra.

M. Ghachoui, A. Soubkiou, R.D. Wells, B.V. Rackham, A. H.M.J. Triaud, D. Sebastian, S. Giacalone, K.G. Stassun, D.R. Ciardi, K.A. Collins, A. Liu, Y. Gómez Maqueo Chew, M. Gillon, Z. Benkhaldoun, L. Delrez, J.D. Eastman, O. Demangeon, K. Barkaoui, A. Burdanov, B.-O. Demory, J. de Wit, G. Dransfield, E. Ducrot, L. Garcia, Y. Gómez Maqueo Chew, M.A. Gómez-Muñoz, M.J. Hooton, E. Jehin, C.A. Murray, P.P. Pedersen, F.J. Pozuelos, D. Queloz, L. Sabin, N. Schanche, M. Timmermans, E.J. Gonzales, C.D. Dressing, C. Aganze, A.J. Burgasser, R. Gerasimov, C. Hsu, C.A. Theissen, D. Charbonneau, J.M. Jenkins, D.W. Latham, G. Ricker, S. Seager, A. Shporer, J.D. Twicken, R. Vanderspek, J.N. Winn, K.I. Collins, A. Fukui, T. Gan, N. Narita, R.P. Schwarz

Accepted in A&A

We report the discovery by the TESS mission of a super-Earth on a 4.8-d orbit around an inactive M4.5 dwarf (TOI-1680) and its validation with ground-based facilities. The host star is located 37.14 pc away, and it has a radius of 0.2100+/-0.0064 R_sun, a mass of 0.1800+/-0.0044 M_sun and an effective temperature of 3211+/-100 K. We validate and characterize the planet using TESS data, ground-based multi-wavelength photometry from TRAPPIST, SPECULOOS and LCO, and high-resolution AO observations from Keck/NIRC2 and Shane. Our analyses determine the planet to have a radius of 1.466+0.063/-0.049 R_earth and an equilibrium temperature of 404+/-14 K, assuming no albedo and perfect heat redistribution. Assuming a mass based on mass-radius relationships, this planet is a promising target for atmospheric characterization with the James Webb Space Telescope (JWST).

Yixiao Zhou, Anish M. Amarsi, Victor Aguirre Børsen-Koch, Klara G. Karlsmose, Remo Collet, Thomas Nordlander

24 pages, 20 figures, accepted for publication in A&A

Three-dimensional radiation-hydrodynamics (3D RHD) simulations of stellar surface convection provide valuable insights into many problems in solar and stellar physics. However, almost all 3D near-surface convection simulations to date are based on solar-scaled chemical compositions, which limit their application on stars with peculiar abundance patterns. To overcome this difficulty, we implement the robust and widely-used FreeEOS equation of state and our Blue opacity package into the Stagger 3D radiation-magnetohydrodynamics code. We present a new 3D RHD model of the solar atmosphere, and demonstrate that the mean stratification as well as the distributions of key physical quantities are in good agreement with those of the latest Stagger solar model atmosphere. The new model is further validated by comparing against solar observations. The new model atmospheres reproduce the observed flux spectrum, continuum centre-to-limb variation, and hydrogen line profiles at a satisfactory level, thereby confirming the realism of the model and the underlying input physics. These implementations open the prospect for studying other stars with different $\alpha$-element abundance, carbon-enhanced metal-poor stars and population II stars with peculiar chemical compositions using 3D Stagger model atmospheres.

Pritha Bari, Nicola Bartolo, Guillem Domènech, Sabino Matarrese

34 pages, 7 figures

This paper explores the physics of second-order gravitational waves (GWs) induced by scalar-tensor perturbation interactions in the radiation-dominated Universe. We investigate the distinctive signatures of these GWs and their detectability compared to scalar-induced GWs. Unlike scalar-scalar induced GWs, scalar-tensor induced GWs do not present resonances or a logarithmic running in the low frequency tail in the case of peaked primordial spectra. But, interestingly, they partly inherit any primordial parity violation of tensor modes. We find that chirality in primordial GWs can lead to distinguishing effects in scalar-tensor induced GWs in the ultraviolet (UV) region. We also address a potential divergence in our GWs and explore possible solutions. This study contributes to our understanding of GWs in the early Universe and their implications for cosmology and GWs detection.

Matheus Bernini-Peron, Wagner L.F. Marcolino, Andreas A.C. Sander, Jean-Claude Bouret, Varsha Ramachandran, Julian Saling, Fabian R.N. Schneider, Lidia M. Oskinova, Francisco Najarro

27 pages, 22 figures, accepted for publication in A&A

B supergiants (BSGs) are evolved stars with effective temperatures between 10 to 30 kK and are important to understand massive star evolution. Located on the edge of the line-driven wind regime, the study of their atmospheres is helpful to understand phenomena such as the bi-stability jump. Key UV features of their spectra have so far not been reproduced by models for types later than B1. Here, we aim to remedy this situation via spectral analysis that accounts for wind clumping and X-rays. In addition, we investigate the evolutionary status of our sample stars based on the obtained stellar parameters. We determined parameters via quantitative spectroscopy using CMFGEN and PoWR codes. The models were compared to UV and optical data of four BSGs: HD206165, HD198478, HD53138, and HD164353. We also study the evolutionary status of our sample using GENEC and MESA tracks. When including clumping and X-rays, we find good agreements between synthetic and observed spectra for our sample stars. For the first time, we reproduced key lines in the UV. For that, we require a moderately clumped wind (f_infty > ~0.5). We also infer relative X-ray luminosities of ~10^-7.5 to 10^-8 -- lower than the typical ratio of 10^-7. Moreover, we find a possible mismatch between evolutionary and spectroscopic masses, which could be related to the mass-discrepancy problem present in other OB stars. Our results provide evidence that X-rays and clumping are needed to describe the winds of cool BSGs. However, their winds seem less structured than in earlier type stars. This aligns with observational X-rays and clumping constraints as well as recent hydrodynamical simulations. The BSGs' evolutionary status appears diverse: some objects are potentially post-red supergiants or merger products. The wind parameters provide evidence for a moderate mass-loss rate increase around the bi-stability jump. Abstract abridged

In this paper we refine a previously developed acoustic-source filter (Bahauddin & Rast 2021), improving its reliability and extending its capabilities. We demonstrate how to fine-tune the filter to meet observational constraints and to focus on specific wavefront speeds. This refinement enables discrimination of acoustic-source depths and tracking of local-source wavefronts, thereby facilitating ultra-local helioseismology on very small scales. By utilizing the photospheric Doppler signal from a subsurface source in a MURaM simulation, we demonstrate that robust ultra-local three-dimensional helioseismic inversions for the granular flows and sound speed to depths of at least 80 km below the photosphere are possible. The capabilities of the National Science Foundation's new Daniel K. Inouye Solar Telescope (DKIST) will enable such measurements of the real Sun.

Stefan Reissl, Ralf S. Klessen, Eric W. Pellegrini, Daniel Rahner, Rüdiger Pakmor, Robert Grand, Facundo Gomez, Federico Marinacci, Volker Springel

14 pages, 3 figures

Magnetic fields are of critical importance for our understanding of the origin and long-term evolution of the Milky Way. This is due to their decisive role in the dynamical evolution of the interstellar medium (ISM) and their influence on the star-formation process. Faraday rotation measures (RM) along many different sightlines across the Galaxy are a primary means to infer the magnetic field topology and strength from observations. However, the interpretation of the data has been hampered by the failure of previous attempts to explain the observations in theoretical models and to synthesize a realistic multi-scale all-sky RM map. We here utilize a cosmological magnetohydrodynamic (MHD) simulation of the formation of the Milky Way, augment it with a novel star cluster population synthesis model for a more realistic structure of the local interstellar medium, and perform detailed polarized radiative transfer calculations on the resulting model. This yields a faithful first principles prediction of the Faraday sky as observed on Earth. The results reproduce the observations of the Galaxy not only on global scales, but also on local scales of individual star-forming clouds. They also imply that the Local Bubble containing our Sun dominates the RM signal over large regions of the sky. Modern cosmological MHD simulations of the Milky Way's formation, combined with a simple and plausible model for the fraction of free electrons in the ISM, explain the RM observations remarkably well, thus indicating the emergence of a firm theoretical understanding of the genesis of magnetic fields in our Universe across cosmic time.

Gen Ye, Alessandra Silvestri

10 pages plut appendix, 4 figures

Recently the international pulsar timing array collaboration has announced the first strong evidence for an isotropic gravitational wave background (GWB). We propose that rapid small oscillations (wiggles) in the Hubble parameter would trigger a resonance with the propagating gravitational waves, leaving novel signature in the GWB spectrum in the form of sharp resonance peaks. The proposed signal can appear at all frequency ranges and is common to continuous spectrum GWBs with arbitrary origin. Due to its resonant nature, the signal strength differs by a perturbation order depending on whether the GWB is primordial or not, which makes it a smoking gun for the primordial origin of the observed GWB. We show that a large part of the parameter space of such signal can be constrained by near future PTA observations, while fitting the signal template to the current NanoGrav 15yr data already hints an interesting feature near 15 nHz.

Karsten Jedamzik, Levon Pogosian

Invited chapter for the edited book "The Hubble Constant Tension" (Eds. E. Di Valentino and D. Brout, Springer Singapore, expected in 2024); 25 pages, 8 figures

Magnetic fields appear to be present in essentially all astrophysical environments, including galaxies, clusters of galaxies and voids. There are both observational and theoretical motives for considering the possibility of their origin tracing back to the events in the very early universe, such as the electroweak phase transition or Inflation. Such a primordial magnetic field (PMF) would remain embedded in the plasma and evolve to persist through the radiation and matter eras, and to the present day. As described in this Chapter, a PMF present in the primordial plasma prior to recombination could help relieve the Hubble tension. A stochastic magnetic field would induce inhomogeneities, pushing the baryons into regions of lower magnetic energy density and speeding up the recombination process. As a consequence, the sound horizon at last scattering would be smaller, which is a necessary ingredient for relieving the Hubble tension. Intriguingly, the strength of the magnetic field required to alleviate the tension is of the right order to also explain the observed magnetic fields in galaxies, clusters of galaxies and voids. These findings motivate further detailed studies of recombination in the presence of PMFs and observational tests of this hypothesis.

Eugenio Carretta (1), Angela Bragaglia (1) ((1) INAF-Osservatorio di Astrofisica e Scienza dello Spazio di Bologna)

19 pages plus 6 pages of Appendix, 23 figures, 5 tables; accepted for publication on Astronomy and Astrophysics

We present the homogeneous abundance analysis for a combined sample of 185 giants in the bulge globular cluster (GC) NGC 6388. Our results are used to describe the multiple stellar populations and differences or analogies with bulge field stars. Proton-capture elements indicate that a single class of first-generation polluters is sufficient to reproduce both the extreme and intermediate parts of the anti-correlations among light elements O, Na, Mg, and Al, which is at odds with our previous results based on a much smaller sample. The abundance pattern of other species in NGC 6388 closely tracks the trends observed in bulge field stars. In particular, the alpha-elements, including Si, rule out an accreted origin for NGC 6388, confirming our previous results based on iron-peak elements, chemo-dynamical analysis, and the age-metallicity relation. The neutron-capture elements are generally uniform, although the [Zr/Fe] ratio shows an intrinsic scatter, correlated to Na and Al abundances. Instead, we do not find enhancement in neutron-capture elements for stars whose photometric properties would classify NGC 6388 as a type II GC. Together with the homogeneity in [Fe/H] we found in a previous paper, this indicates we need to better understand the criteria to separate classes of GCs, coupling photometry, and spectroscopy. These results are based on abundances of 22 species (O, Na, Mg, Al, Si, Ca, Ti, Sc, V, Cr, Mn, Fe, Co, Ni, Zn, Y, Zr, Ba, La, Ce, Nd, and Eu) from UVES spectra sampling proton-, alpha-, neutron-capture elements, and Fe-peak elements. For 12 species, we also obtain abundances in a large number of giants (up to 150) from GIRAFFE spectra.

Ning Xie, Fa Peng Huang

accepted version at SCIENCE CHINA Physics, Mechanics & Astronomy

The axion or axion-like particle motivated from a natural solution of strong CP problem or string theory is a promising dark matter candidate. We study the new observational effects of ultralight axion-like particles by the space-borne gravitational wave detector and the radio telescope. Taking the neutron star-black hole binary as an example, we demonstrate that the gravitational waveform could be obviously modified by the slow depletion of the axion cloud around the black hole formed through the superradiance process. We compare these new effects on the binary with the well-studied effects from dynamical friction with dark matter and dipole radiation in model-independent ways. Finally, we discuss the constraints from LIGO/Virgo and study the detectability of the ultralight axion particles at LISA and TianQin.

The strong CP problem can be solved if the laws of nature are invariant under a space-time parity exchanging the Standard Model with its mirror copy. We review and extend different realizations of this idea with the aim of discussing Dark Matter, neutrino physics, leptogenesis and collider physics within the same context. In the minimal realization of Ref. [1] the mirror world contains a massless dark photon, which leads to a rather interesting cosmology. Mirror electrons reproduce the dark matter abundance for masses between 500-1000 GeV with traces of strongly interacting dark matter. This scenario also predicts deviations from cold dark matter, sizable $\Delta N_{\rm eff}$ and colored states in the TeV range that will be tested in a variety of upcoming experiments. We also explore scenarios where the mirror photon is massive and the mirror particles are charged under ordinary electro-magnetism with very different phenomenology. We also show that, for the measured values of the SM parameters, the Higgs effective potential can give rise to a second minimum at large field value as required to break spontaneously the parity symmetry.

First-order phase transitions, which take place when the symmetries are predominantly broken (and masses are then generated) through radiative corrections, produce observable gravitational waves and primordial black holes. We provide a model-independent approach that is valid for large-enough supercooling to quantitatively describe these phenomena in terms of few parameters, which are computable once the model is specified. The validity of a previously-proposed approach of this sort is extended here to a larger class of theories. Among other things, we identify regions of the parameter space that correspond to the background of gravitational waves recently detected by pulsar timing arrays (NANOGrav, CPTA, EPTA, PPTA) and others that are either excluded by the observing runs of LIGO and Virgo or within the reach of future gravitational wave detectors. Furthermore, we find regions of the parameter space where primordial black holes produced by large over-densities due to such phase transitions can account for dark matter. Finally, it is shown how this model-independent approach can be applied to specific cases, including a phenomenological completion of the Standard Model with right-handed neutrinos and gauged $B-L$ undergoing radiative symmetry breaking.

With focus on the cosmological evolution of linear perturbations of matter and geometry, we calculate the equivalent expressions to that of the Newtonian and Synchronous gauges within the framework of Unimodular Gravity, being these two gauges commonly used and implemented in Boltzmann codes. An important aspect of our analysis is the inclusion of the energy-momentum current violation, as well as its perturbations. Moreover, for the first time we demonstrate that it is possible to fix both gauges consistently, although as it has been already noticed in previous literature, neither of them is recovered in the sense of the dynamics given in General Relativity for matter and metric fluctuations. Specifically, we show that since the unimodular constraint at the level of linear perturbations lead to only one degree of freedom of scalar modes of metric fluctuations, the dynamics in Unimodular Gravity forces to keep the anisotropic stress in the Newtonian gauge, whereas the cold dark matter comoving frame can not be set in the Synchronous gauge. The physical implications on the density contrast of cold dark matter is reviewed, and the Sachs-Wolfe effect is obtained and compared with previous results in the literature of cosmological perturbations in Unimodular Gravity.

Constança Providência, Tuhin Malik, Milena Bastos Albino, Márcio Ferreira

16 pages, 11 figures

The general behavior of the nuclear equation of state (EOS), relevant for the description of neutron stars (NS), is studied within a relativistic mean field description of nuclear matter. Different formulations, both with density dependent couplings and with non-linear mesonic terms, are considered and their predictions compared and discussed. A special attention is drawn to the effect on the neutron star properties of the inclusion of exotic degrees of freedom as hyperons. Properties such as the speed of sound, the trace anomaly, the proton fraction and the onset of direct Urca processes inside neutron stars are discussed. The knowledge of the general behavior of the hadronic equation of state and the implication it has on the neutron star properties will allow to identify signatures of a deconfinement phase transition discussed in other studies.

David Alonso-González, Dorian W. P. Amaral, Adriana Bariego-Quintana, David Cerdeno, Martín de los Rios

22 pages, 10 figures

We explore the complementarity of direct detection (DD) and spallation source (SS) experiments for the study of sterile neutrino physics. We focus on the sterile baryonic neutrino model: an extension of the Standard Model that introduces a massive sterile neutrino with couplings to the quark sector via a new gauge boson. In this scenario, the inelastic scattering of an active neutrino with the target material in both DD and SS experiments gives rise to a characteristic nuclear recoil energy spectrum that can allow for the reconstruction of the neutrino mass in the event of a positive detection. We first derive new bounds on this model based on the data from the COHERENT collaboration on CsI and LAr targets, which we find do not yet probe new areas of the parameter space. We then assess how well future SS experiments will be able to measure the sterile neutrino mass and mixings, showing that masses in the range 15-50 MeV can be reconstructed. We show that there is a degeneracy in the measurement of the sterile neutrino mixing that substantially affects the reconstruction of parameters for masses of the order of 40 MeV. Thanks to their lower energy threshold and sensitivity to the solar tau neutrino flux, DD experiments allow us to partially lift the degeneracy in the sterile neutrino mixings and considerably improve its mass reconstruction down to 9 MeV. Our results demonstrate the excellent complementarity between DD and SS experiments in measuring the sterile neutrino mass and highlight the power of DD experiments in searching for new physics in the neutrino sector.

Ultralight bosons are well-motivated particles from various physical and cosmological theories, and can be spontaneously produced during the superradiant process, forming a dense hydrogen-like cloud around the spinning black hole. After the growth saturates, the cloud slowly depletes its mass through gravitational-wave emission. In this work we study the orbit dynamics of a binary system containing such a gravitational atom saturated in various spin-0,1,2 superradiant states, taking into account both the effects of dynamical friction and the cloud mass depletion. We estimate the significance of mass depletion, finding that although dynamical friction could dominate the inspiral phase, it typically does not affect the outspiral phase driven by the mass depletion. Focusing on the large orbit radius, we investigate the condition to observe the outspiral, and the detectability of the cloud via pulsar-timing signal in the case of black hole-pulsar binary.

Jan Niklas Reinhardt, Martin Staab, Kohei Yamamoto, Jean-Baptiste Bayle, Aurélien Hees, Olaf Hartwig, Karsten Wiesner, Gerhard Heinzel

Interspacecraft ranging is crucial for the suppression of laser frequency noise via time-delay interferometry (TDI). So far, the effect of on-board delays and ambiguities in the LISA ranging observables was neglected in LISA modelling and data processing investigations. In reality, on-board delays cause offsets and timestamping delays in the LISA measurements, and PRN ranging is ambiguous, as it only determines the range up to an integer multiple of the pseudo-random noise (PRN) code length. In this article, we identify the four LISA ranging observables: PRN ranging, the sideband beatnotes at the interspacecraft interferometer, TDI ranging, and ground-based observations. We derive their observation equations in the presence of on-board delays, noise, and ambiguities. We then propose a three-stage ranging sensor fusion to combine these observables in order to gain optimal ranging estimates. We propose to calibrate the on-board delays on ground and to compensate the associated offsets and timestamping delays in an initial data treatment (stage 1). We identify the ranging-related routines, which need to run continuously during operation (stage 2), and implement them numerically. Essentially, this involves the reduction of ranging noise, for which we develop a Kalman filter combining the PRN ranging and the sideband beatnotes. We further implement crosschecks for the PRN ranging ambiguities and offsets (stage 3). We show that both ground-based observations and TDI ranging can be used to resolve the PRN ranging ambiguities. Moreover, we apply TDI ranging to estimate the PRN ranging offsets.

Geatano Lambiase, Tanmay Kumar Poddar

23 pages, 1 figure, comments are welcome

The interaction of neutrinos with ultralight scalar and vector dark matter backgrounds induce a modification of the neutrino dispersion relation. The effects of this modification are reviewed in the framework of asymmetric emission of neutrinos from the supernova core, and, in turn, of pulsar kicks. We consider the neutrino oscillations, focusing in particular to active-sterile conversion. The ultralight dark matter induced neutrino dispersion relation contains a term of the form $\delta {\bf \Omega}\cdot \hat{{\bf{p}}}$, where $\delta {\bf \Omega}$ is related to the ultralight dark matter field and $\hat{{\bf p}}$ is the unit vector along the direction of neutrino momentum. The relative orientation of ${\bf p}$ with respect to $\delta {\bf \Omega}$ affects the mechanism for the generation of the observed pulsar velocities. We obtain the resonance condition for the active-sterile neutrino oscillation in ultralight dark matter background and calculate the star parameters in the resonance surface so that both ultralight scalar and vector dark matter backgrounds can explain the observed pulsar kicks. The asymmetric emission of neutrinos in presence of ultralight dark matter background results gravitational memory signal which can be probed from the gravitational wave detectors. We also establish a connection between the ultralight dark matter parameters and the standard model extension parameter.

We study the recent Physical Review Letter [1] which presents a new mechanism for black hole evaporation through a spatially dependent temperature. This new temperature is comparable to the Hawking result near the black hole, but is very small far away, and therefore could be a small correction. Here we apply the proposed reasoning to the case of de Sitter space, finding that it over predicts the de Sitter temperature of a minimally coupled scalar by factor of $\approx 4.3$ and therefore cannot be ignored in any limit. This indicates an inconsistency in the proposed formalism.

The proliferation of topological defects is a common out-of-equilibrium phenomenon when a system is driven into a phase of broken symmetry. The Kibble-Zurek mechanism (KZM) provides a theoretical framework for the critical dynamics and generation of topological defects in such scenarios. One of the early applications of KZM is the estimation of heavy magnetic monopoles left behind by the cosmological phase transitions in the early universe. The scarcity of such relic monopoles, which contradicts the prediction of KZM, is one of the main motivations for cosmological inflationary theories. On the other hand, magnetic monopoles as emergent quasi-particles have been observed in spin ices, a peculiar class of frustrated magnets that remain disordered at temperatures well below the energy scale of exchange interaction. Here we study the annihilation dynamics of magnetic monopoles when spin ice is cooled to zero temperature in a finite time. Through extensive Glauber dynamics simulations, we find that the density of residual monopole follows a power law dependence on the annealing rate. A kinetic reaction theory that precisely captures the annihilation process from Monte Carlo simulations is developed. We further show that the KZM can be generalized to describe the critical dynamics of spin ice, where the exponent of the power-law behavior is determined by the dynamic critical exponent $z$ and the cooling protocol.

In this paper, we explore the prospect for improving the measurement accuracy of masses and radii of neutron stars. We consider imminent and long-term upgrades of the Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo, as well as next-generation observatories -- the Cosmic Explorer and Einstein Telescope. We find that neutron star radius with single events will be constrained to within roughly 500m with the current generation of detectors and their upgrades. This will improve to 200m, 100m and 50m with a network of observatories that contain one, two or three next-generation observatories, respectively. Combining events in bins of 0.05 solar masses we find that for stiffer (softer) equations-of-state like ALF2 (APR4), a network of three XG observatories will determine the radius to within 30m (100m) over the entire mass range of neutron stars from 1 to 2.0 solar masses (2.2 solar masses), allowed by the respective equations-of-state. Neutron star masses will be measured to within 0.5 percent with three XG observatories irrespective of the actual equation-of-state. Measurement accuracies will be a factor of 4 or 2 worse if the network contains only one or two XG observatories, respectively, and a factor of 10 worse in the case of networks consisting of Advanced LIGO, Virgo KAGRA and their upgrades. Tens to hundreds of high-fidelity events detected by future observatories will allow us to accurately measure the mass-radius curve and hence determine the dense matter equation-of-state to exquisite precision.

Reconstructing urban areas in 3D out of satellite raster images has been a long-standing and challenging goal of both academical and industrial research. The rare methods today achieving this objective at a Level Of Details $2$ rely on procedural approaches based on geometry, and need stereo images and/or LIDAR data as input. We here propose a method for urban 3D reconstruction named KIBS(\textit{Keypoints Inference By Segmentation}), which comprises two novel features: i) a full deep learning approach for the 3D detection of the roof sections, and ii) only one single (non-orthogonal) satellite raster image as model input. This is achieved in two steps: i) by a Mask R-CNN model performing a 2D segmentation of the buildings' roof sections, and after blending these latter segmented pixels within the RGB satellite raster image, ii) by another identical Mask R-CNN model inferring the heights-to-ground of the roof sections' corners via panoptic segmentation, unto full 3D reconstruction of the buildings and city. We demonstrate the potential of the KIBS method by reconstructing different urban areas in a few minutes, with a Jaccard index for the 2D segmentation of individual roof sections of $88.55\%$ and $75.21\%$ on our two data sets resp., and a height's mean error of such correctly segmented pixels for the 3D reconstruction of $1.60$ m and $2.06$ m on our two data sets resp., hence within the LOD2 precision range.

We argue that if axions are the dark matter, their coupling to electromagnetism results in exponential growth of a helical magnetic field when the axion field first rolls down its potential. After an inverse cascade, the relevant length scales to day are of order 10-100 kpc, of astrophysical interest. Our mechanism for allowing the field to grow relies on a nuance of MHD. Faraday's Law says that an electric field is needed to create a magnetic field. Previous authors relied on conventional Ohm's law to calculate E, but the resistivity is negligible and therefore they assume E is as well. We use a modified Ohm's Law that includes the effects of self-induction in limiting the current driven by a given E, which allows a magnetic field to grow.

Giovanni Amelino-Camelia, Domenico Frattulillo, Giulia Gubitosi, Giacomo Rosati, Suzana Bedic

Studies of in-vacuo dispersion are the most active area of quantum-gravity phenomenology. The way in which in-vacuo dispersion produces redshift-dependent corrections to the time of flight of astrophysics particles depends on the model-dependent interplay between Planck-scale effects and spacetime curvature/expansion, and we here derive the most general formula for the leading order redshift-dependent correction to the time of flight for the scenario in which relativistic symmetries are deformed at the Planck scale (DSR). We find that, contrary to the broken symmetries scenario (LIV), where in principle any arbitrary form of redshift dependence could be allowed, for the DSR scenario only linear combinations of three possible forms of redshift dependence are allowed. We also discuss some specific combinations of these three terms whose investigation might deserve priority from the quantum-gravity perspective.