Abstract visibility
Abstract text size

Papers for Tuesday, Jul 11 2023

Papers with local authors

Jiazheng Zhu, Ning Jiang, Tinggui Wang, Shifeng Huang, Zheyu Lin, Yibo Wang, Jian-Guo Wang

9 pages, 6 figures; Accepted for ApJL (July, 2023)

0 votes
Paper 36 — arXiv:2307.04297
0 votes
Paper 36 — arXiv:2307.04297

We report the discovery of a faint optical tidal disruption event (TDE) in the nearby star-forming galaxy NGC 3799. Identification of the TDE is based on its position at the galaxy nucleus, a light curve declining as t^-5/3, a blue continuum with an almost constant blackbody temperature of ~12,000K, and broad (~15,000kms^-1) Balmer lines and characteristic He~II 4686A emission. The light curve of AT 2023clx peaked at an absolute magnitude of -17.16mag in the g-band and a maximum blackbody bolometric luminosity of 4.56*10^42 ergs^-1, making it the faintest TDE discovered to date. With a redshift of 0.01107 and a corresponding luminosity distance of 47.8Mpc, it is also the closest optical TDE ever discovered to our best knowledge. Furthermore, our analysis of Swift/XRT observations of AT 2023clx yields a very tight 3 sigma upper limit of 9.53*10^39 ergs^-1 in the range 0.3--10keV. AT 2023clx, together with very few other faint TDEs such as AT 2020wey, prove that there are probably a large number of faint TDEs yet to be discovered at higher redshifts, which is consistent with the prediction of luminosity functions (LFs). The upcoming deeper optical time-domain surveys, such as the Legacy Survey of Space and Time (LSST) and the Wide-Field Survey Telescope (WFST) will discover more TDEs at even lower luminosities, allowing for a more precise constraint of the low-end of the LF.

Si-Yue Yu, Cheng Cheng, Yue Pan, Fengwu Sun, Yang A. Li

21 pages, 17 figures; A&A in press

0 votes
Paper 61 — arXiv:2307.04753
0 votes
Paper 61 — arXiv:2307.04753

Observations of high-redshift galaxies with unprecedented detail have now been rendered possible with JWST. However, accurately quantifying their morphology remains uncertain due to potential biases and uncertainties. To address this issue, we used a sample of 1816 nearby DESI galaxies, with a mass range of $10^{9.75-11.25}M_{\odot}$, to compute artificial images of galaxies of the same mass located at $0.75\leq z\leq 3$ and observed at rest-frame optical wavelength in CEERS. We analyzed the effects of cosmological redshift on the measurements of Petrosian radius ($R_p$), half-light radius ($R_{50}$), asymmetry ($A$), concentration ($C$), axis ratio ($q$), and S\'ersic index ($n$). Our results show that $R_p$ and $R_{50}$, calculated using non-parametric methods, are slightly overestimated due to PSF smoothing, while $R_{50}$, $q$, and $n$ obtained through model fitting does not exhibit significant biases. We improve the computation of $A$ by incorporating a more accurate noise effect removal procedure. Due to PSF asymmetry, there is a minor overestimation of $A$ for intrinsically symmetric galaxies. However, for intrinsically asymmetric galaxies, PSF smoothing dominates and results in an underestimation of $A$, an effect that becomes more significant with higher intrinsic $A$ or at lower resolutions. Moreover, PSF smoothing also leads to an underestimation of $C$, which is notably more pronounced in galaxies with higher intrinsic $C$ or at lower resolutions. We developed functions based on resolution level, defined as $R_p/$FWHM, for correcting these biases and the associated statistical uncertainties. Applying these corrections, we measured the bias-corrected morphology for the simulated CEERS images and we find that the derived quantities are in good agreement with their intrinsic values -- except for $A$, which is robust only for angularly large galaxies where $R_p/{\rm FWHM}\geq 5$.

All other papers

Stefan H. Strub, Luigi Ferraioli, Cédric Schmelzbach, Simon C. Stähler, Domenico Giardini

13 pages, 11 figures

The Laser Interferometer Space Antenna (LISA) is a planned space-based gravitational wave telescope with the goal of measuring gravitational waves in the milli-Hertz frequency band, which is dominated by millions of Galactic binaries. While some of these binaries produce signals that are loud enough to stand out and be extracted, most of them blur into a confusion foreground. Current methods for analyzing the full frequency band recorded by LISA to extract as many Galactic binaries as possible and to obtain Bayesian posterior distributions for each of the signals are computationally expensive. We introduce a new approach to accelerate the extraction of the best fitting solutions for Galactic binaries across the entire frequency band from data with multiple overlapping signals. Furthermore, we use these best fitting solutions to omit the burn-in stage of a Markov chain Monte Carlo method and to take full advantage of GPU-accelerated signal simulation, allowing us to compute posterior distributions in 2 seconds per signal on a laptop-grade GPU.

T. Eftekhari, W. Fong, A. C. Gordon, N. Sridhar, C. D. Kilpatrick, S. Bhandari, A. T. Deller, Y. Dong, A. Rouco Escorial, K. E. Heintz, J. Leja, B. Margalit, B. D. Metzger, A. B. Pearlman, J. X. Prochaska, S. D. Ryder, P. Scholz, R. M. Shannon, N. Tejos

15 pages, 5 figures; submitted to ApJ

We present the first X-ray census of fast radio burst (FRB) host galaxies to conduct the deepest search for AGN and X-ray counterparts to date. Our sample includes seven well-localized FRBs with unambiguous host associations and existing deep Chandra observations, including two events for which we present new observations. We find evidence for AGN in two FRB host galaxies based on the presence of X-ray emission coincident with their centers, including the detection of a luminous ($L_X \approx 5 \times 10^{42} \ \rm erg \ s^{-1}$) X-ray source at the nucleus of FRB20190608B's host, for which we infer an SMBH mass of $\rm M_{\rm BH} \sim 10^{8} \ M_{\odot}$ and an Eddington ratio $L_{\rm bol}/ L_{\rm Edd} \approx 0.02$, characteristic of geometrically thin disks in Seyfert galaxies. We also report nebular emission line fluxes for 24 highly secure FRB hosts (including 10 hosts for the first time), and assess their placement on a BPT diagram, finding that FRB hosts trace the underlying galaxy population. We further find that the hosts of repeating FRBs are not confined to the star-forming locus, contrary to previous findings. Finally, we place constraints on associated X-ray counterparts to FRBs in the context of ultraluminous X-ray sources (ULXs), and find that existing X-ray limits for FRBs rule out ULXs brighter than $L_X \gtrsim 10^{40} \ \rm erg \ s^{-1}$. Leveraging the CHIME/FRB catalog and existing ULX catalogs, we search for spatially coincident ULX-FRB pairs. We identify two ULX in the galaxy NGC 2633 that are spatially coincident with the repeating FRB20180908B and for which the DM-inferred redshift is comparable to the distance of the galaxy, assuming a $\rm DM_{host}$ contribution of $150 \ \rm pc \ cm^{-3}$.

C.M. Harrison (Newcastle University), A. Girdhar, S.R. Ward

Accepted for publication in IAU Symposium 378 Conference Proceedings "Black Hole Winds at all Scales"

Cosmological simulations fail to reproduce realistic galaxy populations without energy injection from active galactic nuclei (AGN) into the interstellar medium (ISM) and circumgalactic medium (CGM); a process called `AGN feedback'. Consequently, observational work searches for evidence that luminous AGN impact their host galaxies. Here, we review some of this work. Multi-phase AGN outflows are common, some with potential for significant impact. Additionally, multiple feedback channels can be observed simultaneously; e.g., radio jets from `radio quiet' quasars can inject turbulence on ISM scales, and displace CGM-scale molecular gas. However, caution must be taken comparing outflows to simulations (e.g., kinetic coupling efficiencies) to infer feedback potential, due to a lack of comparable predictions. Furthermore, some work claims limited evidence for feedback because AGN live in gas-rich, star-forming galaxies. However, simulations do not predict instantaneous, global impact on molecular gas or star formation. The impact is expected to be cumulative, over multiple episodes.

Louise Welsh, Ryan Cooke, Michele Fumagalli, Max Pettini

15 pages, 12 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

We present new Keck/HIRES data of the most metal-poor damped Lyman-alpha (DLA) system currently known. By targeting the strongest accessible Fe II features, we have improved the upper limit of the [Fe/H] abundance determination by ~1 dex, finding [Fe/H]<-3.66 (2 sigma). We also provide the first upper limit on the relative abundance of an odd-atomic number element for this system [Al/H]<-3.82 (2 sigma). Our analysis thus confirms that this z_abs=3.07 DLA is not only the most metal-poor DLA but also the most iron-poor DLA currently known. We use the chemistry of this DLA, combined with a stochastic chemical enrichment model, to probe its enrichment history. We find that this DLA is best modelled by the yields of an individual Population III progenitor rather than multiple Population III stars. We then draw comparisons with other relic environments and, particularly, the stars within nearby ultra-faint dwarf galaxies. We identify a star within Bootes I, with a similar chemistry to that of the DLA presented here, suggesting that it may have been born in a gas cloud that had similar properties. The extremely metal-poor DLA at redshift z_abs=3.07 (i.e. ~2 Gyrs after the Big Bang) may reside in one of the least polluted environments in the early Universe.

Giorgio Laverda, Javier Rubio

24 pages, 10 figures

A Hubble-induced phase transition is a natural spontaneous symmetry breaking mechanism allowing for explosive particle production in non-oscillatory models of inflation involving non-minimally coupled spectator fields. In this work, we perform a comprehensive characterisation of this type of transitions as a tachyonic Ricci-heating mechanism, significantly extending previous results in the literature. By performing $\mathcal{O}(100)$ 3+1-dimensional classical lattice simulations, we explore the parameter space of two exemplary scenarios, numerically determining the main timescales in the process. Based on these results, we formulate a set of parametric equations that offer a practical approach for determining the efficiency of the heating process, the temperature at the onset of radiation domination, and the minimum number of e-folds of inflation needed to resolve the flatness and horizon problems in specific quintessential inflation scenarios. These parametric equations eliminate the need for additional lattice simulations, providing a convenient and efficient method for evaluating these key quantities.

Christopher Tiede, Daniel J. D'Orazio

17 pages, 14 figures

Modern numerical hydrodynamics tools have recently enabled detailed examinations of binaries accreting from prograde circumbinary disks that have re-framed the current understanding of binary-disk interactions and disk driven orbital evolution. We present the first full-domain grid-based hydrodynamics simulations of equal-mass, eccentric binaries accreting from retrograde circumbinary disks. We study binary eccentricities that span $e=0.0$ to $e = 0.8$ continuously, and explore the influence of retrograde accretion on the binary orbital response, disk morphology, and observational properties. We find that, at all eccentricities, retrograde accretion shrinks the binary semi-major axis and pumps its eccentricity leading to the previously identified possibility of highly eccentric mergers. Contrary to past studies and models, we observe gravitational forces to dominate the binary's orbital evolution as opposed to the physical accretion of mass and momentum. Retrograde accretion variability also differs strongly from prograde solutions. Preeminently, binaries with $e > 0.55$ reveal a unique two-period, double-peaked accretion signature that has not previously been identified. We additionally find evidence for the emergence of retrograde Lindblad resonances at large eccentricities in accordance with predictions from linear theory. Our results suggest that some astrophysical binaries where retrograde accretion is possible will experience factors of a few times faster orbital decay than in prograde disks and will have their eccentricities pumped beyond the limits found from prograde solutions. Such effects could lead to rapid inward migration for some young stellar binaries, the detection of highly-eccentric LISA mergers, and the tentatively observed turnover at the low-frequency end of the gravitational wave background.

Anna Niemiec, Mathilde Jauzac, Dominique Eckert, David Lagattuta, Keren Sharon, Anton M. Koekemoer, Keiichi Umetsu, Ana Acebron, Jose M. Diego, David Harvey, Eric Jullo, Vasily Kokorev, Marceau Limousin, Guillaume Mahler, Priyamvada Natarajan, Mario Nonino, Juan D. Remolina, Charles Steinhardt, Sut-Ieng Tam, Adi Zitrin

29 pages, 17 figures, 3 tables

The HST treasury program BUFFALO provides extended wide-field imaging of the six Hubble Frontier Fields galaxy clusters. Here we present the combined strong and weak-lensing analysis of Abell 370, a massive cluster at z=0.375. From the reconstructed total projected mass distribution in the 6arcmin x 6arcmin BUFFALO field-of-view, we obtain the distribution of massive substructures outside the cluster core and report the presence of a total of seven candidates, each with mass $\sim 5 \times 10^{13}M_{\odot}$. Combining the total mass distribution derived from lensing with multi-wavelength data, we evaluate the physical significance of each candidate substructure, and conclude that 5 out of the 7 substructure candidates seem reliable, and that the mass distribution in Abell 370 is extended along the North-West and South-East directions. While this finding is in general agreement with previous studies, our detailed spatial reconstruction provides new insights into the complex mass distribution at large cluster-centric radius. We explore the impact of the extended mass reconstruction on the model of the cluster core and in particular, we attempt to physically explain the presence of an important external shear component, necessary to obtain a low root-mean-square separation between the model-predicted and observed positions of the multiple images in the cluster core. The substructures can only account for up to half the amplitude of the external shear, suggesting that more effort is needed to fully replace it by more physically motivated mass components. We provide public access to all the lensing data used as well as the different lens models.

Kevin France, Nicole Arulanantham, Erin Maloney, P. Wilson Cauley, P. Abraham, Juan M. Alcala, Justyn Campbell-White, Eleonora Fiorellino, Gregory J. Herczeg, Brunella Nisini, Miguel Vioque

19 pages, 11 figures. Accepted to the Astronomical Journal

The spatial distribution and evolution of gas in the inner 10 au of protoplanetary disks form the basis for estimating the initial conditions of planet formation. Among the most important constraints derived from spectroscopic observations of the inner disk are the radial distributions of the major gas phase constituents, how the properties of the gas change with inner disk dust evolution, and how chemical abundances and excitation conditions are influenced by the high-energy radiation from the central star. We present a survey of the radial distribution, excitation, and evolution of inner disk molecular hydrogen (H$_{2}$) obtained as part of the $HST$/ULLYSES program. We analyze far-ultraviolet spectroscopy of 71 (63 accreting) pre-main sequence systems in the ULLYSES DR5 release to characterize the H$_{2}$ emission lines, H$_{2}$ dissociation continuum emission, and major photochemical/disk evolution driving UV emissions (Ly$\alpha$, UV continuum, and C IV). We use the widths of the H$_{2}$ emission lines to show that most fluorescent H$_{2}$ arises between 0.1 - 1.4 au from the parent star, and show positive correlations of the average emitting radius with the accretion luminosity and with the dust disk mass. We find a strong correlation between H$_{2}$ dissociation emission and both the accretion-dominated Ly$\alpha$ luminosity and the inner disk dust clearing, painting a picture where water molecules in the inner 3 au are exposed to and dissociated by strong Ly$\alpha$ emission as the opacity of the inner disk declines with time.

Andrea Banzatti, Klaus M. Pontoppidan, John Carr, Evan Jellison, Ilaria Pascucci, Joan Najita, Carlos E. Munoz-Romero, Karin I. Oberg, Anusha Kalyaan, Paola Pinilla, Sebastiaan Krijt, Feng Long, Michiel Lambrechts, Giovanni Rosotti, Gregory J. Herczeg, Colette Salyk, Ke Zhang, Nick Ballering, Michael R. Meyer, Simon Bruderer, the JDISCS collaboration

Posted as submitted to ApJ Letters; feedback and input from the community is welcome

Previous analyses of mid-infrared water spectra from young protoplanetary disks observed with the Spitzer-IRS found an anti-correlation between water luminosity and the millimeter dust disk radius observed with ALMA. This trend was suggested to be evidence for a fundamental process of inner disk water enrichment, used to explain properties of the Solar System 40 years ago, in which icy pebbles drift inward from the outer disk and sublimate after crossing the snowline. Previous analyses of IRS water spectra, however, were very uncertain due to the low spectral resolution that blended lines together. We present new JWST-MIRI spectra of four disks, two compact and two large with multiple radial gaps, selected to test the scenario that water vapor inside the snowline is regulated by pebble drift. The higher spectral resolving power of MIRI-MRS now yields water spectra that separate individual lines, tracing upper level energies from 900 K to 10,000 K. These spectra clearly reveal excess emission in the low-energy lines in compact disks, compared to the large disks, establishing the presence of a cooler component with $T \approx$ 170-400 K and equivalent emitting radius $R_{\rm{eq}}\approx$ 1-10 au. We interpret the cool water emission as ice sublimation and vapor diffusion near the snowline, suggesting that there is indeed a higher inwards mass flux of icy pebbles in compact disks. Observation of this process opens up multiple exciting prospects to study planet formation chemistry in inner disks with JWST.

Navin Chaurasiya, Surhud More, Shogo Ishikawa, Shogo Masaki, Daichi Kashino, Teppei Okumura

We infer the connection between the stellar mass of galaxies from the Subaru Hyper Suprime-Cam (HSC) survey, and their dark matter halo masses and its evolution in two bins of redshifts between $[0.3, 0.8]$. We use the measurements of the weak lensing signal of galaxies using background sources from the Year 1 shape catalog from the HSC survey. We bin galaxies in stellar mass with varying thresholds ranging from $8.6 \leq \log [ M_*/(h^{-2} {M_\odot})] \leq 11.2$ and use stringent cuts in the selection of source galaxies to measure the weak lensing signal. We model these measurements of the weak lensing signal together with the abundance of galaxies in the halo occupation distribution framework. We obtain constraints on the halo occupation parameters of central galaxies $M_{\rm min}$ and $\sigma_{\log M}$, which correspond to the halo mass at which central galaxies for each threshold sample reach half occupancy, and its scatter, respectively, along with parameters that describe the occupation of the satellite galaxies. The measurements of abundance and weak lensing individually constrain different degeneracy directions in the $M_{\rm min}$ and $\sigma_{\log M}$ plane, thus breaking the degeneracy in these parameters. We demonstrate that the weak lensing measurements are best able to constrain the average central halo masses, $\langle M_{\rm cen} \rangle$. We compare our measurements to those obtained using the abundance and clustering of these galaxies as well as the subhalo abundance matching measurements and demonstrate qualitative agreement. We find that the galaxy-dark matter connection does not vary significantly between redshift bins we explore in this study. Uncertainties in the photometric redshift of the lens galaxies imply that more efforts are required to understand the true underlying stellar mass-halo mass relation of galaxies and its evolution over cosmic epoch.

Morgan T. Chidester, F.X. Timmes, Ebraheem Farag

18 pages, 10 total figures, 2 online interactive figures, zenodo data set link included. Accepted to the Astrophysical Journal

We consider the combined effects that overshooting and the 12C({\alpha}, {\gamma})16O reaction rate have on variable white dwarf stellar models. We find that carbon-oxygen white dwarf models continue to yield pulsation signatures of the current experimental 12C({\alpha}, {\gamma})16O reaction rate probability distribution function when overshooting is included in the evolution. These signatures hold because the resonating mantle region, encompassing $\simeq$\,0.2\,\Msun\ in a typical $\simeq$\,0.6\,\Msun\ white dwarf model, still undergoes radiative helium burning during the evolution to a white dwarf. Our specific models show two potential low-order adiabatic g-modes, $g_2$ and $g_6$, that signalize the 12C({\alpha}, {\gamma})16O reaction rate probability distribution function. Both g-mode signatures induce average relative period shifts of $\Delta P/P = 0.44 \%$ and $\Delta P/P = 1.33\%$ for $g_2$ and $g_6$ respectively. We find that $g_6$ is a trapped mode, and the $g_2$ period signature is inversely proportional to the 12C({\alpha}, {\gamma})16O reaction rate. The $g_6$ period signature generally separates the slower and faster reaction rates, and has a maximum relative period shift of $\Delta P/P = 3.45\%$. We conclude that low-order g-mode periods from carbon-oxygen white dwarfs may still serve as viable probes for the 12C({\alpha}, {\gamma})16O reaction rate probability distribution function when overshooting is included in the evolution.

Galaxy and galaxy clusters exhibit tight robust physical scaling relations between baryons and system dynamics. One such phenomenon is mass discrepancy with two leading solution spaces occupied by LCDM and MOND. Here, we propose an alternative solution to this puzzling problem exclusively based on application of the scalar virial theorem. For these dynamically equilibrated systems, we demonstrate there is ample virially-induced kinetic energy available to modify bulk structure dynamics in apparent violation of Newtonian law. We propose the ubiquitous Baryonic Tully-Fisher Relation represents the preferred dynamic configuration that best assures long-term survivability for these thermodynamic quasi-equilibrated systems. We compare total mass estimates guided by the empirical evidence to those obtained from NFW dark matter halo fits ranging from small dwarf galaxies to massive galaxy clusters.

Somsawat Rattanasoon, Eugene Semenko, David Mkrtichian, Saran Poshyachinda

This paper was submitted as a contribution to the proceedings of the third Belgo-Indian Network for Astronomy and Astrophysics (BINA) workshop, which was held in Bhimtal, Uttarakhand (India) on 22-24 March 2023. The final, peer-reviewed version will be published in Bulletin de la Soci\'et\'e Royale des Sciences de Li\`ege. The manuscript contains 5 pages, 3 figures, and 4 references

The National Astronomical Research Institute of Thailand (NARIT) has a manifold network of small telescopes installed worldwide. These telescopes serve educational and research purposes and are equipped mainly with CCD detectors for direct imaging and photometry. To extend the possible field of applications, several telescopes were fitted with commercially available medium-resolution spectrographs eShel from Shelyak. With these devices, researchers in NARIT obtained a versatile tool for stellar spectroscopy. Here we describe the current status of available equipment, possible ways of upgrading, and briefly introduce the achieved results of the asteroseismologic study of fast-rotating stars.

Kevin Wagner, Jordan Stone, Andrew Skemer, Steve Ertel, Ruobing Dong, Dániel Apai, Eckhart Spalding, Jarron Leisenring, Michael Sitko, Kaitlin Kratter, Travis Barman, Mark Marley, Brittany Miles, Anthony Boccaletti, Korash Assani, Ammar Bayyari, Taichi Uyama, Charles E. Woodward, Phil Hinz, Zackery Briesemeister, Kellen Lawson, François Ménard, Eric Pantin, Ray W. Russell, Michael Skrutskie, John Wisniewski

Published in Nature Astronomy

Understanding the driving forces behind spiral arms in protoplanetary disks remains a challenge due to the faintness of young giant planets. MWC 758 hosts such a protoplanetary disk with a two-armed spiral pattern that is suggested to be driven by an external giant planet. We present new thermal infrared observations that are uniquely sensitive to redder (i.e., colder or more attenuated) planets than past observations at shorter wavelengths. We detect a giant protoplanet, MWC 758c, at a projected separation of ~100 au from the star. The spectrum of MWC 758c is distinct from the rest of the disk and consistent with emission from a planetary atmosphere with Teff = 500 +/- 100 K for a low level of extinction (AV<30), or a hotter object with a higher level of extinction. Both scenarios are commensurate with the predicted properties of the companion responsible for driving the spiral arms. MWC 758c provides evidence that spiral arms in protoplanetary disks can be caused by cold giant planets or by those whose optical emission is highly attenuated. MWC 758c stands out both as one of the youngest giant planets known, and also as one of the coldest and/or most attenuated. Furthermore, MWC 758c is among the first planets to be observed within a system hosting a protoplanetary disk.

A. E. Camisasca, I. A. Steele, M. Bulla, C. Guidorzi, M. Shrestha

9 pages, 8 figures, published in MNRAS

We consider the optimisation of the observing strategy (cadence, exposure time and filter choice) using medium size (2-m class) optical telescopes in the follow-up of kilonovae localised with arcminute accuracy to be able to distinguish among various kilonova models and viewing angles. To develop an efficient observation plan, we made use of the synthetic light curves obtained with the Monte Carlo radiative transfer code POSSIS for different kilonova models and as a function of different viewing angles and distances. By adding the appropriate photon counting noise to the synthetic light curves, we analysed four alternative sequences having the same total time exposure of 8 hours, with different time windows (0.5, 1, 2, 4 h), each with $i$, $r$, and $u$ filters, to determine the observing sequence that maximises the chance of a correct identification of the model parameters. We suggest to avoid $u$ filter and to avoid the use of colour curves. We also found that, if the error on distance is $\le$ 2%, 0.5, 1, 2-hour time window sequences are equivalent, so we suggest to use 2-hour one, because it has 1 day cadence, so it can be easily realised. When the distance of the source is unknown, 0.5 h time window sequence is preferable.

In our previous work, we searched for super-flares on different types of stars, focusing on G-type dwarfs using entire Kepler data to study statistical properties of the occurrence rate of super-flares. The said study also considered how the statistics change with stellar rotation period, which in turn, had to be determined. Using such new data, as a by-product, we found 138 Kepler IDs of F and G types main sequence stars with rotation periods less than a day ($P_{\rm rot}<1$ d). On one hand, previous studies have revealed short activity cycles in F-type and G-type stars and the question investigated was whether or not short-term activity cycles are a common phenomenon in these stars. On the other hand, extensive studies exist which establish empirical connection between a star's activity cycle and rotation periods. In this study, we compile all available Kepler data with $P_{\rm rot}<1$ d and derive, as well as use plausible, established empirical relations between $P_{\rm cyc}$ and $P_{\rm rot}$ with the aim to provide predictions for very short $5.13\leq P_{\rm cyc}\leq 38.14$ d cases in a tabular form. As a result, we invite others to measure $P_{\rm cyc}$ using monitoring program of stellar activity (e.g. activity-related chromospheric emission S-index) or similar means for the Kepler IDs found in this study in order put to test the derived and/or established empirical relations between $P_{\rm cyc}$ and $P_{\rm rot}$. We also propose an alternative method for measuring very short $P_{\rm cyc}$, using flare-detection algorithms applied to future space mission data.

R. Andrassy, G. Leidi, J. Higl, P. V. F. Edelmann, F. R. N. Schneider, F. K. Roepke

10 pages, 12 figures, submitted to A&A

There is strong observational evidence that convective cores of intermediate-mass and massive main-sequence stars are substantially larger than standard stellar-evolution models predict. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-time-scale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number Seven-League Hydro (SLH) code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a $15$ M$_\odot$ zero-age main-sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor ranging from $10^3$ to $10^6$. After reaching thermal equilibrium, all of our models show a clear penetration layer. Its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. As the luminosity is decreased, the penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification. This structure corresponds to the adiabatic ,,step overshoot'' model often employed in stellar-evolution calculations. The thickness of the penetration layer slowly decreases with decreasing luminosity. Depending on how we extrapolate our 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from $0.09$ to $0.44$ pressure scale heights, which are broadly compatible with observations.

Đorđe V. Savić, Isidora Jankov, Weixiang Yu, Vincenzo Petrecca, Matthew J. Temple, Qingling Ni, Raphael Shirley, Andjelka B. Kovacevic, Mladen Nikolic, Dragana Ilic, Luka C. Popovic, Maurizio Paolillo, Swayamtrupta Panda, Aleksandra Ciprijanovic, Gordon T. Richards

Accepted by ApJ. 21 pages, 14 figures, 5 tables

Development of the Rubin Observatory Legacy Survey of Space and Time (LSST) includes a series of Data Challenges (DC) arranged by various LSST Scientific Collaborations (SC) that are taking place during the projects preoperational phase. The AGN Science Collaboration Data Challenge (AGNSCDC) is a partial prototype of the expected LSST AGN data, aimed at validating machine learning approaches for AGN selection and characterization in large surveys like LSST. The AGNSC-DC took part in 2021 focusing on accuracy, robustness, and scalability. The training and the blinded datasets were constructed to mimic the future LSST release catalogs using the data from the Sloan Digital Sky Survey Stripe 82 region and the XMM-Newton Large Scale Structure Survey region. Data features were divided into astrometry, photometry, color, morphology, redshift and class label with the addition of variability features and images. We present the results of four DC submitted solutions using both classical and machine learning methods. We systematically test the performance of supervised (support vector machine, random forest, extreme gradient boosting, artificial neural network, convolutional neural network) and unsupervised (deep embedding clustering) models when applied to the problem of classifying/clustering sources as stars, galaxies or AGNs. We obtained classification accuracy 97.5% for supervised and clustering accuracy 96.0% for unsupervised models and 95.0% with a classic approach for a blinded dataset. We find that variability features significantly improve the accuracy of the trained models and correlation analysis among different bands enables a fast and inexpensive first order selection of quasar candidates

Tiancheng Sun, Zhishuai Ge, Xunzhou Chen, Shaolan Bi, Tanda Li, Xianfei Zhang, Yaguang Li, Yaqian Wu, Sarah A. Bird, Ferguson J. W., Jianzhao Zhou, Lifei Ye, Liu Long, Jinghua Zhang

18 pages, 12 figures, accepted for publication in ApJS

Varying oxygen abundance could impact the modeling-inferred ages. This work aims to estimate the ages of dwarfs considering observed oxygen abundance. To characterize 67,503 LAMOST and 4,006 GALAH FGK-type dwarf stars, we construct a grid of stellar models which take into account oxygen abundance as an independent model input. Compared with ages determined with commonly-used $\alpha$-enhanced models, we find a difference of $\sim$9% on average when the observed oxygen abundance is considered. The age differences between the two types of models are correlated to [Fe/H] and [O/$\alpha$], and they are relatively significant on stars with [Fe/H] $\lesssim$ -0.6 dex. Generally, varying 0.2 dex in [O/$\alpha$] will alter the age estimates of metal-rich (-0.2 $<$ [Fe/H] $<$ 0.2) stars by $\sim$10%, and relatively metal-poor (-1 $<$ [Fe/H] $<$ -0.2) stars by $\sim$15%. Of the low-O stars with [Fe/H] $<$ 0.1 dex and [O/$\alpha$] $\sim$ -0.2 dex, many have fractional age differences of $\geq$ 10%, and even reach up to 27%. The fractional age difference of high-O stars with [O/$\alpha$] $\sim$ 0.4 dex reaches up to -33% to -42% at [Fe/H] $\lesssim$ -0.6 dex. We also analyze the chemical properties of these stars. We find a decreasing trend of [Fe/H] with age from 7.5-9 Gyr to 5-6.5 Gyr for the stars from the LAMOST and GALAH. The [O/Fe] of these stars increases with decreasing age from 7.5-9 Gyr to 3-4 Gyr, indicating that the younger population is more O-rich.

Germán Molpeceres, Joan Enrique-Romero, Yuri Aikawa

Accepted for publication in Astronomy and Astrophysics

CO2 is one of the dominant components of the interstellar ice. Recent observations show CO2 exists more abundantly in polar (H2O-dominated) ice than in apolar (H2O-poor) ice. CO2 ice formation is primarily attributed to the reaction between CO and OH, which has a barrier. Highly accurate quantum chemical calculations were employed to analyze the stationary points of the potential energy surfaces of the title reaction in the gas phase on a H2O and CO clusters. Microcanonical transition state theory was used as a diagnostic tool for the efficiency of the reaction under ISM conditions. We simulate the kinetics of ice chemistry, considering different scenarios involving non-thermal processes and energy dissipation. The CO + OH reaction proceeds through the remarkably stable intermediate HOCO radical. On the H2O cluster, the formation of this intermediate is efficient, but the subsequent reaction leading to CO2 formation is not. Conversely, HOCO formation on the CO cluster is inefficient without external energy input. Thus, CO2 ice cannot be formed by the title reaction alone either on the H2O cluster or CO cluster. In the polar ice, CO2 ice formation is possible via CO + OH -> HOCO, followed by HOCO + H ->CO2 + H2, as demonstrated by abundant experimental literature. In apolar ice, CO2 formation is less efficient because HOCO formation requires external energy. Our finding is consistent with the JWST observations. Further experimental work is encouraged using low-temperature OH radicals.

Mingxu Sun, Bingqiu Chen, Helong Guo, He Zhao, Ming Yang, Wenyuan Cui

16 pages, 10 figures, Accepted by The Astronomical Journal on 7/8/2023

Modern large-scale photometric surveys have provided us with multi-band photometries of billions of stars. Determining the stellar atmospheric parameters, such as the effective temperature (\teff) and metallicities (\feh), absolute magnitudes ($M_{G}$), distances ($d$) and reddening values (\ebr) is fundamental to study the stellar populations, structure, kinematics and chemistry of the Galaxy. This work constructed an empirical stellar library which maps the stellar parameters to multi-band photometries from a dataset with Gaia parallaxes, LAMOST atmospheric parameters, and optical to near-infrared photometry from several photometric surveys. Based on the stellar library, we developed a new algorithm, SPar (\textbf{S}tellar \textbf{P}arameters from multib\textbf{a}nd photomet\textbf{r}y), which fits the multi-band stellar photometries to derive the stellar parameters (\teff, \feh, $M_G$, $d$ and \ebr) of the individual stars. The algorithm is applied to the multi-band photometric measurements of a sample of stars selected from the SMSS survey, which have stellar parameters derived from the spectroscopic surveys. The stellar parameters derived from multi-band photometries by our algorithm are in good agreement with those from the spectroscopic surveys. The typical differences between our results and the literature values are 170\,K for \teff, 0.23\,dex for \feh, 0.13\,mag for $M_G$ and 0.05\,mag for \ebr. The algorithm proved to be robust and effective and will be applied to the data of future large-scale photometric surveys such as the Mephisto and CSST surveys.

Linghua Xie, Nicola R. Napolitano, Xiaotong Guo, Crescenzo Tortora, Haicheng Feng, Antonios Katsianis, Rui Li, Sirui Wu, Mario Radovich, Leslie K. Hunt, Yang Wang, Lin Tang, Baitian Tang, Zhiqi Huang

39 pages,4 tables, 14 figures, SCIENCE CHINA Physics, Mechanics & Astronomy Accepted

The Kilo Degree Survey (KiDS) is currently the only sky survey providing optical ($ugri$) plus near-infrared (NIR, $ZYHJK_S$) seeing matched photometry over an area larger than 1000 $\rm deg^2$. This is obtained by incorporating the NIR data from the VISTA Kilo Degree Infrared Galaxy (VIKING) survey, covering the same KiDS footprint. As such, the KiDS multi-wavelength photometry represents a unique dataset to test the ability of stellar population models to return robust photometric stellar mass ($M_*$) and star-formation rate (SFR) estimates. Here we use a spectroscopic sample of galaxies for which we possess $u g r i Z Y J H K_s$ ``gaussianized'' magnitudes from KiDS data release 4. We fit the spectral energy distribution from the 9-band photometry using: 1) three different popular libraries of stellar {population} templates, 2) single burst, simple and delayed exponential star-formation history models, and 3) a wide range of priors on age and metallicity. As template fitting codes we use two popular softwares: LePhare and CIGALE. We investigate the variance of the stellar masses and the star-formation rates from the different combinations of templates, star formation recipes and codes to assess the stability of these estimates and define some ``robust'' median quantities to be included in the upcoming KiDS data releases. As a science validation test, we derive the mass function, the star formation rate function, and the SFR-$M_*$ relation for a low-redshift ($z<0.5$) sample of galaxies, that result in excellent agreement with previous literature data. The final catalog, containing $\sim290\,000$ galaxies with redshift $0.01<z<0.9$, is made publicly available.

Yue Shao, Yidong Xu, Yougang Wang, Wenxiu Yang, Ran Li, Xin Zhang, Xuelei Chen

53 pages, 6 figures, 2 extended data figures, 9 supplementary figures. Published in Nature Astronomy

The absorption features in spectra of high-redshift background radio sources, caused by hyperfine structure lines of hydrogen atoms in the intervening structures, are known collectively as the 21-cm forest. They provide a unique probe of small-scale structures during the epoch of reionization, and can be used to constrain the properties of the dark matter (DM) thought to govern small-scale structure formation. However, the signals are easily suppressed by heating processes that are degenerate with a warm DM model. Here we propose a probe of both the DM particle mass and the heating history of the Universe, using the one-dimensional power spectrum of the 21-cm forest. The one-dimensional power spectrum measurement not only breaks the DM model degeneracy but also increases the sensitivity, making the probe actually feasible. Making 21-cm forest observations with the upcoming Square Kilometre Array has the potential to simultaneously determine both the DM particle mass and the heating level in the early Universe, shedding light on the nature of DM and the first galaxies.

E.F. Ocran, M. Vaccari, J.M. Stil, A.R. Taylor, C.H. Ishwara-Chandra, Jae-Woo Kim

17 pages, 13 figures, accepted for publication in MNRAS

We measure the specific star formation rates of \textit{K}-band selected galaxies from the ELAIS-N1 by stacking GMRT data at 610 MHz. We identify a sample of SFGs, spanning $\rm{0.1\leq\,\textit{z}\,\leq\,1.5}$ and $\rm{10^{8.5}<\,{\textit{M}_{\star}}/{\textit{M}_{\odot}}<10^{12.4}}$, using a combination of multi-wavelength diagnostics obtained from the deep LoTSS multi-wavelength catalogue. We measure the flux densities in the radio map and estimate the radio SFR in order to probe the nature of the galaxies below the noise and confusion limits. The massive galaxies in our sample have the lowest sSFRs which is in agreement with previous studies. For the different populations, we show that the sSFR-mass relation steepens with redshift, with an average slope of $\rm{\langle \beta_{All} \rangle\,=\, -0.49\pm0.01}$ for the whole sample, and $\rm{\langle \beta_{SFG} \rangle\,=\, -0.42\pm0.02}$ for the SFGs. Our results indicate that galaxy populations undergo 'downsizing', whereby most massive galaxies form their stars earlier and more rapidly than low mass galaxies. Both populations show a strong decrease in their sSFR toward the present epoch. The sSFR evolution with redshift is best described by a power law $\rm{(1\,+\,\textit{z})^\textit{n}}$, where $\rm{\langle \textit{n}_{ALL}\rangle\sim4.94\pm0.53}$ for all galaxies, and $\rm{\langle \textit{n}_{SFG}\rangle \sim3.51\pm0.52}$ for SFGs. Comparing our measured sSFRs to results from literature, we find a general agreement in the \textit{sSFR-M$_{\star}$} plane.

Jordan Conrad-Burton, Alon Shabi, Sivan Ginzburg

8 pages, submitted to MNRAS, comments welcome

Black widows and redbacks are binary millisecond pulsars with close low-mass companions that are irradiated and gradually ablated by the pulsar's high-energy luminosity $L_{\rm irr}$. These binaries evolve primarily through magnetic braking, which extracts orbital angular momentum and pushes the companion to overflow its Roche lobe. Here, we use the stellar evolution code MESA to examine how the irradiation modifies the companion's structure. Strong $L_{\rm irr}$ inhibits convection to the extent that otherwise fully convective stars become almost fully radiative. By computing the convective velocities and assuming a dynamo mechanism, we find that the thin convective envelopes of such strongly irradiated companions ($L_{\rm irr}\gtrsim 3\,{\rm L}_\odot$) generate much weaker magnetic fields than previously thought - halting binary evolution. With our improved magnetic braking model, we explain most observed black widow and redback companions as remnants of main-sequence stars. We also apply our model (with $L_{\rm irr}$) to evolved companions that overflow their Roche lobe close to the end of their main-sequence phase. The evolutionary tracks of such companions bifurcate, explaining the shortest period systems (which are potential gravitational wave sources) as well as the longest period ones (which are the progenitors of common pulsar-white dwarf binaries). The variety of black widow structures and evolutionary trajectories may be utilized to calibrate the dependence of magnetic braking on the size of the convective layer and on the existence of a radiative-convective boundary, with implications for single stars as well as other binaries, such as cataclysmic variables and AM Canum Venaticorum stars.

L.E. Kristensen, B. Godard, P. Guillard, A. Gusdorf, G. Pineau des Forets

Published in A&A

(Abridged) H2 is the most abundant molecule in the Universe. Thanks to its widely spaced energy levels, it predominantly lights up in warm gas, T > 100 K, such as shocked regions, and it is one of the key targets of JWST observations. These include shocks from protostellar outflows, all the way up to starburst galaxies and AGN. Shock models are able to simulate H2 emission. We aim to explore H2 excitation using such models, and to test over which parameter space distinct signatures are produced in H2 emission. We present simulated H2 emission using the Paris-Durham shock code over an extensive grid of 14,000 plane-parallel stationary shock models, a large subset of which are exposed to an external UV radiation field. The grid samples 6 input parameters: preshock density, shock velocity, transverse magnetic field strength, UV radiation field strength, cosmic-ray-ionization rate, and PAH abundance. Physical quantities, such as temperature, density, and width, have been extracted along with H2 integrated line intensities. The strength of the transverse magnetic field, set by the scaling factor, b, plays a key role in the excitation of H2. At low values of b (<~ 0.3, J-type shocks), H2 excitation is dominated by vibrationally excited lines; at higher values (b >~ 1, C-type shocks), rotational lines dominate the spectrum for shocks with an external radiation field comparable to (or lower than) the solar neighborhood. Shocks with b >= 1 can be spatially resolved with JWST for nearby objects. When the input kinetic energy flux increases, the excitation and integrated intensity of H2 increases similarly. An external UV field mainly serves to increase the excitation, particularly for shocks where the input radiation energy is comparable to the input kinetic energy flux. These results provide an overview of the energetic reprocessing of input kinetic energy flux and the resulting H2 line emission.

Donald V. Reames

19 pages, 9 figures, submitted to Frontiers of Space Science

The earliest evidence on spatial distributions of solar energetic particles (SEPs) compared events from many different source longitudes on the Sun, but the early Pioneers provided the first evidence of the large areas of equal SEP intensities across the magnetically-confined "reservoirs" late in the events. More-detailed measurements of the importance of self-generated waves and trapping structures around the shock waves that accelerate SEPs were obtained from the Helios mission plus IMP 8, especially during the year when the two Voyager spacecraft also happened by. The extent of the dozen widest SEP events in a solar cycle, that effectively wrap around the Sun, was revealed by the widely separated STEREO spacecraft with three-point intensities fit to Gaussians. Element abundances of the broadest SEP events favor average coronal element abundances with little evidence of heavy-element-enhanced "impulsive suprathermal" ions that often dominate the seed population of the shocks, even in extremely energetic local events. However, it is hard to define a distribution with two or three points. Advancing the physics of SEPs may require a return to the closer spacing of the Helios era with coverage mapped by a half-dozen spacecraft to help disentangle the distribution of the SEPs from the underlying structure of the magnetic field and the accelerating shock.

K. McDonough, K. Hughes, D. Smith, A. G. Vieregg

18 pages, 7 figures

The origin of the diffuse astrophysical neutrino flux measured by the IceCube Observatory remains largely unknown. Although NGC 1068 and TXS 0506+056 have been identified as potential neutrino sources, the diffuse flux of neutrinos must have additional sources that have not yet been identified. Here we investigate potential correlations between IceCube's neutrino events and the Fermi and MOJAVE source catalogs, using the publicly-available IceCube data set. We perform three separate spatially-dependent, energy-dependent, and time-dependent searches, and find no statistically significant sources outside of NGC 1068. We find that no more than 13% of IceCube's neutrino flux originates from blazars over the whole sky. Then, using an energy-dependent likelihood analysis, the limit on neutrinos originating from blazars reduces to 9% in the Northern hemisphere. Finally, we set limits on individual sources from the MOJAVE radio catalog after finding no statistically significant time-flaring sources.

Athulya M. P., Anuj Nandi

19 pages, 9 figures, 6 tables. Accepted for publication in MNRAS

GRS 1915+105 is observed in an 'obscured' phase since May 2019, exhibiting steady and low X-ray luminosities, while being intervened by sporadic re-brightenings. In this work, we perform a comprehensive and wide-band analysis of the spectral and timing properties of the source during the period $2019-2021$ using AstroSat (SXT: $0.5-8$ keV; LAXPC: $3-60$ keV), NICER ($0.5-12$ keV), and NuSTAR ($3-60$ keV) observations. Spectral analysis reveals the presence of a highly variable obscurer (N$_{H_{1}}\sim~10^{22} - 10^{24}$ atoms cm$^{-2}$) throughout the observation period. Source is detected in the Low/Hard state for most of the time, with the spectra being described by a Comptonised component ($\Gamma \sim 1.16 - 1.79$, kT$_{e}\sim 2-31$ keV). The source spectra steepen ($\Gamma\sim2.5$) indicating softening of the spectrum during the rise of the re-brightenings. Various emission and absorption lines corresponding to the neutral Fe-K$\alpha$, Fe-XXV K$\alpha$, Fe-XXVI K$\alpha$, and the Ni-XXVIII K$\alpha$ were detected with equivalent widths varying between 70 eV $-$ 3.5 keV. The column density of the absorbing plasma varied between $10^{16} - 10^{18}$ atoms cm$^{-2}$ at a distance $\leq2\times$10$^{10}$ cm. Interestingly, the source is also seen exhibiting various variability classes ($\rho, \lambda, \delta, \chi$) at relatively low luminosities ($\sim$0.01L$_{Edd}$) during the re-brightening phases. Different variability classes show signature of QPOs ($\nu_{QPO}$: 20--180 mHz, rms$_{QPO}$: 7.5% - 16%). The source showed a maximum bolometric luminosity {(L$_{bol}$)} of $\sim$0.01L$_{Edd}$ (Re-brightening phases) and a minimum L$_{bol}$ of 0.004L$_{Edd}$ (Quiet phase) during the period. We discuss the possible disc dynamics around the black hole during this low-luminosity `obscured' phase.

Federico Di Vruno, Vincenza Tornatore

10 pages, 9 figures

Large LEO satellite constellations (or so-called Mega-constellations) will significantly change the view of the sky in some radio frequency bands. For VGOS telescopes it is important to understand the potential impact these constellations will have in their operations, what is the risk of its receivers going into non-linear behaviour and how much additional power would a telescope receive if observing in the same frequencies where satellites are transmitting. This work describes three of these new constellations (as they would look fully deployed) and summarizes the results of a particular study considering two VGOS telescopes (Onsala and Wettzell).

R. Marques-Chaves, D. Schaerer, A. Kuruvanthodi, D. Korber, N. Prantzos, C. Charbonnel, A. Weibel, Y. I. Izotov, M. Messa, G. Brammer, M. Dessauges-Zavadsky, P. Oesch

Submitted to A&A, 19 pages, 8 figures, 4 tables

[Abridged] Using the JWST/NIRSpec observations from CEERS we found an extreme N-emitter, CEERS-1019 at z=8.6782 showing intense NIV and NIII emission. From the observed rest-UV and optical lines we conclude that it is compatible with photoionization from stars and we determine accurate abundances for C, N, O, and Ne, relative to H, finding a highly supersolar ratio log(N/O) = -0.18+/-0.11, and normal log(C/O) = -0.75+/-0.11 and log(Ne/O) = -0.63+/-0.07, for its low metallicity, 12+log(O/H)= 7.70+/-0.18. We also analyze other N-emitters from the literature. All show strongly enhanced N/O ratios and two of them normal C/O. Massive star ejecta from WR stars are needed to explain the galaxies with enhanced C/O (Lynx arc and Mrk 996). On the other hand, supermassive stars (>1000 Msun, SMS) in the ``conveyer-belt model'' put forward to explain globular clusters (GCs), predict a high N/O and small changes in C/O, compatible with CEERS-1019, the Sunburst cluster, SMACS2031, and GN-z11. Based on the chemical abundances, possible enrichment scenarios, compactness, and high ISM density, we suggest that CEERS-1019, SMACS2031, and the Sunburst cluster could contain proto-GCs. Finally, we propose that some N-emitters enriched by SMS could also have formed intermediate-mass black holes, and we suggest that this might be the case for GN-z11. Our observations and analysis reinforce the suggested link between some N-emitters and proto-GC formation, which is supported both by empirical evidence and quantitative models. Furthermore, the observations provide possible evidence for the presence of supermassive stars in the early Universe (z>8) and at z~2-3. Our analysis also suggests that the origin and nature of the N-emitters is diverse, including also objects like GN-z11 which possibly host an AGN.

Hancheng Li, Roland Walter, Nicolas Produit, Fiona Hubert

Synthesis models of the diffuse Cosmic X-ray Background (CXB) suggest that it can be resolved into discrete sources, primarily Active Galactic Nuclei (AGNs). Measuring the CXB accurately offers a unique probe to study the AGN population in the nearby Universe. Current hard X-ray instruments suffer from the time-dependent background and cross-calibration issues. As a result, their measurements of the CXB normalization have an uncertainty of the order of $\sim$15%. In this paper, we present the concept and simulated performances of a CXB detector, which could be operated on different platforms. With a 16-U CubeSat mission running for more than two years in space, such a detector could measure the CXB normalization with $\sim$1% uncertainty.

P. Mitra, O. Scholten, T. N. G. Trinh, S. Buitink, J. Bhavani, A. Corstanje, M. Desmet, H. Falcke, B. M. Hare, J. R. Hörandel, T. Huege, N. Karastathis, G. K. Krampah, K. Mulrey, A. Nelles, H. Pandya, S. Thoudam, K. D. de Vries, S. ter Veen

Measuring the radio emission from cosmic ray particle cascades has proven to be a very efficient method to determine their properties such as the mass composition. Efficient modeling of the radio emission from air showers is crucial in order to extract the cosmic ray physics parameters from the measured radio emission. MGMR3D is a fast semi-analytic code that calculates the complete radio footprint, i.e.\ intensity, polarization, and pulse shapes, for a parametrized shower-current density and can be used in a chi-square optimization to fit a given radio data. It is many orders of magnitude faster than its Monte Carlo counterparts. We provide a detailed comparative study of MGMR3D to Monte Carlo simulations, where, with improved parametrizations, the shower maximum $\Xmax$ is found to have very strong agreement with a small dependency on the incoming zenith angle of the shower. Another interesting feature we observe with MGMR3D is sensitivity to the shape of the longitudinal profile in addition to $\Xmax$. This is achieved by probing the distinguishable radio footprint produced by a shower having a different longitudinal profile than usual. Furthermore, for the first time, we show the results of reconstructing shower parameters for LOFAR data using MGMR3D, and obtaining a $\Xmax$ resolution of 22 g/cm$^2$ and energy resolution of 19\%.

Jhonathan O. Murcia Piñeros, Antônio F. Bertachini de Almeida Prado, Walter Abrahão dos Santos, Rodolpho Vilhena de Moraes

The perturbative integral method was applied to quantify the contribution of external forces during a specific interval of time in trajectories of spacecraft around asteroids and under the Luni-solar influence. However, this method has not been used to quantify the contributions of drag in aerocapture and aerobraking. For this reason, the planet Mars is selected to apply this method during an aerogravity-assisted maneuver. Several trajectories are analyzed, making use of a drag device with area to mass ratios varying from 0.0 to 20.0 m2/kg, simulating solar sails or de-orbit devices. The mathematical model is based in the restricted three-body problem. The use of this maneuver makes it possible to obtain the variations of energy in the trajectory, replacing expensive maneuvers based on fuel consumption. To observe the effects of the maneuvers, different values of pericenter velocity and altitude were selected for prograde and retrograde orbits. The innovation of this research is the application of an integral method to quantify the delta-V of the aero gravity maneuver, comparing the cost of the maneuver with the traditional methods of space propulsion. The results allow the identification of orbits with conditions to capture, and the perturbative maps show the velocity variations.

The Bondi accretion rate of black holes in our and nearby galaxies Messier 87, NGC 3115, NGC 1600, and Cygnus A have been determined or constrained using Chandra or other observations. It, however, remains unknown how much mass from the Bondi radius reaches each black hole and how much is accreted. We determine the accretion rate and radiative efficiency for each black hole using two well-tested empirical relations: one relates a black hole's accretion rate to its mass and redshift, and the other relates the radiative efficiency to the Eddington ratio and redshift. We get an accretion rate of ~0.00002 solar mass per year and radiative efficiency of ~0.9 for Sagittarius A* and an accretion rate of ~0.09 solar masses per year and radiative efficiency of ~0.68 for NGC 1600; and values in between these extremes for the rest. The derived mass inflow rate onto each black hole (not the accretion rate) essentially matches the reported Bondi accretion rate. Thus, the results do not support the ADIOS and CDAF models, but whether the dissipated energy not reflected in a black hole's observed luminosity is advected as in the ADAF model remains uncertain. Furthermore, contrary to current model expectations, the derived radiative efficiencies are orders of magnitude higher and radiative efficiency increases as the accretion rate decreases and a BH ages. A physical basis is found relating the empirical formulation of accretion rate to Bondi accretion.

C.J. Fluke, D. Vohl, V.A. Kilborn, C. Murugeshan

21 pages, 10 figures, Accepted for publication in the Publications of the Astronomical Society of Australia

Next generation astronomical surveys naturally pose challenges for human-centred visualisation and analysis workflows that currently rely on the use of standard desktop display environments. While a significant fraction of the data preparation and analysis will be taken care of by automated pipelines, crucial steps of knowledge discovery can still only be achieved through various level of human interpretation. As the number of sources in a survey grows, there is need to both modify and simplify repetitive visualisation processes that need to be completed for each source. As tasks such as per-source quality control, candidate rejection, and morphological classification all share a single instruction, multiple data (SIMD) work pattern, they are amenable to a parallel solution. Selecting extragalactic neutral hydrogen (HI) surveys as a representative example, we use system performance benchmarking and the visual data and reasoning (VDAR) methodology from the field of information visualisation to evaluate a bespoke comparative visualisation environment: the encube visual analytics framework deployed on the 83 Megapixel Swinburne Discovery Wall. Through benchmarking using spectral cube data from existing HI surveys, we are able to perform interactive comparative visualisation via texture-based volume rendering of 180 three-dimensional (3D) data cubes at a time. The time to load a configuration of spectral cubes scale linearly with the number of voxels, with independent samples of 180 cubes (8.4 Gigavoxels or 34 Gigabytes) each loading in under 5 minutes. We show that parallel comparative inspection is a productive and time-saving technique which can reduce the time taken to complete SIMD-style visual tasks currently performed at the desktop by at least two orders of magnitude, potentially rendering some labour-intensive desktop-based workflows obsolete.

Recent observations of weak gravitational lensing surveys indicate a smoother Universe compared to the predictions of the Cosmic Microwave Background (CMB). This is known as $\sigma_8$ tension or $S_8$ tension, where $\sigma_8$ represents the present root-mean-square matter fluctuation averaged over a sphere of radius $8 h^{-1} \mathrm{Mpc}$ and $S_8 \equiv \sigma_8\sqrt{\Omega_m/0.3}$. In this letter, we investigate a kind of general Dirac-Born-Infeld (DBI) Lagrangian referred as surface-type DBI (s-DBI) model. We have found that, up to the linear order, the constraints on the s-DBI model with CMB from Planck2018 and low-redshift probes (WL and GC) yield $S_8= 0.7685_{-0.0066}^{+0.0077}$ and $S_8=0.766_{-0.0376}^{+0.0471}$, respectively, which are not only self-consistent but also consistent with the values derived from most low-redshift probes. Furthermore, we provide an outlook for searching the non-linear effects of this model, which could be helpful to resolve other issues by Cold Dark Matter on small scales.

N. Shagatova, A. Skopal, E. Kundra, R. Komžík, S. Yu. Shugarov, T. Pribulla, V. Krushevska

12 pages, 10 figures

Context. Non-dusty late-type giants without a corona and large-scale pulsations represent objects that do not fulfil the conditions under which standard mass-loss mechanisms can be applied efficiently. The driving mechanism of their winds is still unknown. Aims. The main goal of this work is to match the radial velocities of absorbing matter with a depth in the red giant (RG) atmosphere in the S-type symbiotic star EG And. Methods. We measured fluxes and radial velocities of ten FeI absorption lines from spectroscopic observations with a resolution of ~30 000. At selected orbital phases, we modelled their broadened profiles, including all significant broadening mechanisms. Results. The selected FeI absorption lines at 5151 - 6469A, originate at a radial distance ~1.03 RG radii from its centre. The corresponding radial velocity is typically ~1 km/s , which represents a few percent of the terminal velocity of the RG wind. The high scatter of the radial velocities of several km/s in the narrow layer of the stellar atmosphere points to the complex nature of the near-surface wind mass flow. The average rotational velocity of 11 km/s implies that the rotation of the donor star can contribute to observed focusing the wind towards the orbital plane. The orbital variability of the absorbed flux indicates the highest column densities of the wind in the area between the binary components, even though the absorbing neutral material is geometrically more extended from the opposite side of the giant. This wind density asymmetry in the orbital plane region can be ascribed to gravitational focusing by the white dwarf companion. Conclusions. Our results suggest that both gravitational and rotational focusing contribute to the observed enhancement of the RG wind towards the orbital plane, which makes mass transfer by the stellar wind highly efficient.

R. Abbasi, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, M. Ahrens, J.M. Alameddine, A. A. Alves Jr., N. M. Amin, K. Andeen, T. Anderson, G. Anton, C. Argüelles, Y. Ashida, S. Athanasiadou, S. Axani, X. Bai, A. Balagopal V., S. W. Barwick, V. Basu, S. Baur, R. Bay, J. J. Beatty, K.-H. Becker, J. Becker Tjus, J. Beise, C. Bellenghi, S. Benda, S. BenZvi, D. Berley, E. Bernardini, D. Z. Besson, G. Binder, D. Bindig, E. Blaufuss, S. Blot, M. Boddenberg, F. Bontempo, J. Y. Book, J. Borowka, S. Böser, O. Botner, J. Böttcher, E. Bourbeau, F. Bradascio, J. Braun, B. Brinson, S. Bron, J. Brostean-Kaiser, R. T. Burley, R. S. Busse, M. A. Campana, E. G. Carnie-Bronca, C. Chen, Z. Chen, D. Chirkin, K. Choi, B. A. Clark, K. Clark, L. Classen, A. Coleman, G. H. Collin, A. Connolly, J. M. Conrad, et al. (325 additional authors not shown)

Submitted on May 12th, 2022; Accepted on May 4th, 2023

The origin of high-energy cosmic rays, atomic nuclei that continuously impact Earth's atmosphere, has been a mystery for over a century. Due to deflection in interstellar magnetic fields, cosmic rays from the Milky Way arrive at Earth from random directions. However, near their sources and during propagation, cosmic rays interact with matter and produce high-energy neutrinos. We search for neutrino emission using machine learning techniques applied to ten years of data from the IceCube Neutrino Observatory. We identify neutrino emission from the Galactic plane at the 4.5$\sigma$ level of significance, by comparing diffuse emission models to a background-only hypothesis. The signal is consistent with modeled diffuse emission from the Galactic plane, but could also arise from a population of unresolved point sources.

Adriana Pisarčíková, Pavol Matlovič, Juraj Tóth, Stefan Loehle, Ludovic Ferrière, David Leiser, Felix Grigat, Jérémie Vaubaillon

Fragments of small solar system bodies entering Earth's atmosphere have possibly been important contributors of organic compounds to the early Earth. The cyano radical (CN) emission from meteors is considered as potentially one of the most suitable markers of organic compounds in meteoroids, however, its detection in meteor spectra has been thus far unsuccessful. With the aim to improve our abilities to identify CN emission in meteor observations and use its spectral features to characterize the composition of incoming asteroidal meteoroids, we present a detailed analysis of CN emission from high-resolution spectra of 22 laboratory simulated meteors including ordinary, carbonaceous, and enstatite chondrites, as well as a large diversity of achondrites (i.e., ureilite, aubrite, lunar, martian, howardite, eucrite, and diogenite), mesosiderite, and iron meteorites. We describe the variations of CN emission from different classes of asteroidal meteor analogues, its correlation and time evolution relative to other major meteoroid components. We demonstrate that CN can be used as a diagnostic spectral feature of carbonaceous and carbon-rich meteoroids, while most ordinary chondrites show no signs of CN. Our results point out strong correlation between CN and H emission and suggest both volatile features are suitable to trace contents of organic matter and water molecules present within meteoroids. For the application in lower resolution meteor observations, we demonstrate that CN can be best recognized in the early stages of ablation and for carbon-rich materials by measuring relative intensity ratio of CN band peak to the nearby Fe I-4 lines.

Yang Huang, Timothy C. Beers, Hai-Bo Yuan, Ke-Feng Tan, Wei Wang, Jie Zheng, Chun Li, Young Sun Lee, Hai-Ning Li, Jing-Kun Zhao, Xiang-Xiang Xue, Yu-Juan Liu, Hua-Wei Zhang, Xue-Ang Sun, Ji Li, Hong-Rui Gu, Christian Wolf, Christopher A. Onken, Ji-Feng Liu, Zhou Fan, Gang Zhao

14 pages, 14 figures, 3 tables, accepted by ApJ. arXiv admin note: text overlap with arXiv:2104.14154

We present precise photometric estimates of stellar parameters, including effective temperature, metallicity, luminosity classification, distance, and stellar age, for nearly 26 million stars using the methodology developed in the first paper of this series, based on the stellar colors from the Stellar Abundances and Galactic Evolution Survey (SAGES) DR1 and Gaia EDR3. The optimal design of stellar-parameter sensitive $uv$ filters by SAGES has enabled us to determine photometric-metallicity estimates down to $-3.5$, similar to our previous results with the SkyMapper Southern Survey (SMSS), yielding a large sample of over five million metal-poor (MP; [Fe/H]$\le -1.0$) stars and nearly one million very metal-poor (VMP; [Fe/H]$\le -2.0$) stars. The typical precision is around $0.1$ dex for both dwarf and giant stars with [Fe/H]$>-1.0$, and 0.15-0.25/0.3-0.4 dex for dwarf/giant stars with [Fe/H]$<-1.0$. Using the precise parallax measurements and stellar colors from Gaia, effective temperature, luminosity classification, distance and stellar age are further derived for our sample stars. This huge data set in the Northern sky from SAGES, together with similar data in the Southern sky from SMSS, will greatly advance our understanding of the Milky Way, in particular its formation and evolution.

M. Agundez, N. Marcelino, B. Tercero, I. Jimenez-Serra, J. Cernicharo

Accepted for publication in A&A

We report new observations of molecular anions with the Yebes 40m and IRAM 30m telescopes toward the cold dense clouds TMC-1 CP, Lupus-1A, L1527, L483, L1495B, and L1544. We detected for the first time C3N- and C5N- in Lupus-1A and C4H- and C6H- in L483. In addition, we report new lines of C6H- toward the six targeted sources, of C4H- toward TMC-1 CP, Lupus-1A, and L1527, and of C8H- and C3N- in TMC-1 CP. Excitation calculations indicate that the lines of anions accessible to radiotelescopes run from subthermally excited to thermalized as the size of the anion increases, with the degree of departure from thermalization depending on the H2 volume density and the line frequency. We noticed that the collision rate coefficients available for the radical C6H cannot explain various observational facts, which advises for a revisitation of the collision data for this species. The observations presented here, together with observational data from the literature, are used to model the excitation of interstellar anions and to constrain their abundances. In general, the anion-to-neutral ratios derived here agree within 50 % (a factor of two at most) with literature values, when available, except for the C4H-/C4H ratio, which shows higher differences due to a revision of the dipole moment of C4H. From the set of anion-to-neutral abundance ratios derived two conclusions can be drawn. First, the C6H-/C6H ratio shows a tentative trend in which it increases with increasing H2 density, as expected from theoretical grounds. And second, it is incontestable that the higher the molecular size the higher the anion-to-neutral ratio, which supports a formation mechanism based on radiative electron attachment. Nonetheless, calculated rate coefficients for electron attachment to the medium size species C4H and C3N are probably too high and too low, respectively, by more than one order of magnitude.

Mikhail Kovalev, Olivier R. Hainaut, Xuefei Chen, Zhanwen Han

accepted in MNRAS letters 2023 July 07

We present the detection of false positive double-lined spectroscopic binaries candidates (SB2) using medium-resolution survey (MRS) spectra from the one time-domain field of LAMOST data release 10 (DR10). The secondary component in all these binaries has near zero radial velocity and solar-like spectral lines. Highly likely this is light from the semi-transparent clouds illuminated by the full Moon. However we also suspect that partially this contamination can be caused by a solar light reflected from the surface of low-orbital artificial satellites launched in the beginning of 2022. We found several possible contaminant candidates using archival orbital data. We propose measures to reduce risk of such contamination for the future observations and methods to find it in archived ones.

Sujoy Kumar Nath, Dipak Debnath, Kaushik Chatterjee, Riya Bhowmick, Hsiang-Kuang Chang, Sandip K. Chakrabarti

15 pages, 8 Figures, 2 Tables (In Communication ApJ)

We study the recent outburst of the black hole candidate EXO 1846-031 which went into an outburst in 2019 after almost 34 years in quiescence. We use archival data from Swift/XRT, MAXI/GSC, NICER/XTI and NuSTAR/FPM satellites/instruments to study the evolution of the spectral and temporal properties of the source during the outburst. Low energy X-ray flux of the outburst shows multiple peaks making it a multipeak outburst. Evolving type-C quasi-periodic oscillations (QPOs) are observed in the NICER data in the hard, hard intermediate and soft intermediate states. We use the physical Two Component Advective Flow (TCAF) model to analyze the combined spectra of multiple satellite instruments. According to the TCAF model, the accreting matter is divided into Keplerian and sub-Keplerian parts, and the variation in the observed spectra in different spectral states arises out of the variable contributions of these two types of accreting matter in the total accretion rate. Studying the evolution of the accretion rates and other properties of the accretion flow obtained from the spectral analysis, we show how the multiple peaks in the outburst flux arises out of discontinuous supply and different radial velocities of two types of accreting matter from the pile-up radius. We detect an Fe emission line at $\sim6.6$ keV in the hard and the intermediate states in the NICER spectra. We determine the probable mass of the black hole to be $12.43^{+0.14}_{-0.03}~M_\odot$ from the spectral analysis with the TCAF model. We also estimate viscous time scale of the source in this outburst to be $\sim 8$ days from the peak difference of the Keplerian and sub-Keplerian mass accretion rates.

L. Mashonkina, Yu. Pakhomov, T. Sitnova, A. Smogorzhevskii, P. Jablonka, V. Hill

11 pages, 10 figures, 3 tables, revised version after the referee's positive report, MNRAS

Detailed chemical abundances of very metal-poor (VMP, [Fe/H] < -2) stars are important for better understanding the First Stars, early star formation and chemical enrichment of galaxies. Big on-going and coming high-resolution spectroscopic surveys provide a wealth of material that needs to be carefully analysed. For VMP stars, their elemental abundances should be derived based on the non-local thermodynamic equilibrium (non-LTE = NLTE) line formation because low metal abundances and low electron number density in the atmosphere produce the physical conditions favorable for the departures from LTE. The galactic archaeology research requires homogeneous determinations of chemical abundances. For this purpose, we present grids of the 1D-NLTE abundance corrections for the Na I, Mg I, Ca I, Ca II, Ti II, Fe I, Zn I, Zn II, Sr II, and Ba II lines, which are used in the galactic archaeology research. The range of atmospheric parameters represents VMP stars on various evolutionary stages and covers effective temperatures from 4000 to 6500~K, surface gravities from log g = 0.5 to log g = 5.0, and metallicities $-5.0 \le$ [Fe/H] $\le -2.0$. The data is publicly available, and we provide the tools for interpolating in the grids online.

S. Zamora, A. I. Díaz

23 pages, 21 figures, accepted for publication in MNRAS;

This work aims to derive the physical properties of the CNSFRs in the ring of the face-on spiral NGC 7742 using IFS observations. We have selected 88 individual ionising clusters that power HII regions populating the ring of the galaxy that may have originated in a minor merger event. For the HII regions the rate of Lyman continuum photon emission is between 0.025 and 1.5 10$^{51}$ which points to these regions being ionised by star clusters. Their electron density, ionisation parameter, filling factor and ionised hydrogen mass show values consistent with those found in other studies of similar regions and their metal abundances as traced by sulphur have been found to be between 0.25 and 2.4 times solar, with most regions showing values slightly below solar. The equivalent temperature of the ionising clusters is relatively low, below 40000 K which is consistent with the high elemental abundances derived. The young stellar population of the clusters has contributions of ionising and non-ionising populations with ages around 5 Ma and 300 Ma respectively. The masses of ionising clusters once corrected for the contribution of underlying non-ionising populations were found to have a mean value of 3.5 $\times$ 10$^4$ M$_{\odot}$, comparable to the mass of ionised gas and about 20 \% of the corrected photometric mass.

David Hendriks, Robert Izzard

19 pages, 21 figures includings appendices. Accepted for publication in MNRAS

Mass-transfer interactions in binary stars can lead to accretion disk formation, mass loss from the system and spin-up of the accretor. To determine the trajectory of the mass-transfer stream, and whether it directly impacts the accretor, or forms an accretion disk, requires numerical simulations. The mass-transfer stream is approximately ballistic, and analytic approximations based on such trajectories are used in many binary population synthesis codes as well as in detailed stellar evolution codes. We use binary population synthesis to explore the conditions under which mass transfer takes place. We then solve the reduced three-body equations to compute the trajectory of a particle in the stream for systems with varying system mass ratio, donor synchronicity and initial stream velocity. Our results show that on average both more mass and more time is spent during mass transfer from a sub-synchronous donor than from a synchronous donor. Moreover, we find that at low initial stream velocity the asynchronous rotation of the donor leads to self-accretion over a large range of mass ratios, especially for super-synchronous donors. The stream (self-)intersects in a narrow region of parameter space where it transitions between accreting onto the donor or the accretor. Increasing the initial stream velocity leads to larger areas of the parameter space where the stream accretes onto the accretor, but also more (self-)intersection. The radii of closest approach generally increase, but the range of specific angular momenta that these trajectories carry at the radius of closest approach gets broader. Our results are made publicly available.

Judy Chebly, Julián D. Alvarado-Gómez, Katja Poppenhäger, Cecilia Garraffo

21 pages, 10 figures

As a cool star evolves, it loses mass and angular momentum due to magnetized stellar winds which affect its rotational evolution. This change has consequences that range from the alteration of its activity to influences over the atmosphere of any orbiting planet. Despite their importance, observations constraining the properties of stellar winds in cool stars are extremely limited. Therefore, numerical simulations provide a valuable way to understand the structure and properties of these winds. In this work, we simulate the magnetized winds of 21 cool main-sequence stars (F-type to M-dwarfs), using a state-of-the-art 3D MHD code driven by observed large-scale magnetic field distributions. We perform a qualitative and quantitative characterization of our solutions, analyzing the dependencies between the driving conditions (e.g., spectral type, rotation, magnetic field strength) and the resulting stellar wind parameters (e.g., Alfv\'en surface size, mass loss rate, angular momentum loss rate, stellar wind speeds). We compare our models with the current observational knowledge on stellar winds in cool stars and explore the behaviour of the mass loss rate as a function of the Rossby number. Furthermore, our 3D models encompass the entire classical Habitable Zones (HZ) of all the stars in our sample. This allows us to provide the stellar wind dynamic pressure at both edges of the HZ and analyze the variations of this parameter across spectral type and orbital inclination. The results here presented could serve to inform future studies of stellar wind-magnetosphere interactions and stellar wind erosion of planetary atmospheres via ion escape processe.

Angelo Ricarte, Ramesh Narayan, Brandon Curd

19 pages, 7 figures, 1 table. Currently in review in ApJL

A spinning black hole accreting from a disk of strongly magnetized plasma via a magnetically arrested disk is known to produce an efficient electromagnetic jet powered by the black hole's spin energy. We present general relativistic radiative magnetohydrodynamic simulations of magnetically arrested systems covering a range of sub- to super-Eddington accretion rates. Using the numerical results from these simulation, we develop formulae to describe the magnetization, jet efficiency, and spin evolution of an accreting black hole as a function of its spin and accretion rate. A black hole with near-Eddington accretion experiences a mild degree of spin-down because of angular momentum loss through the jet, leading to an equilibrium spin of 0.8 rather than 1.0 at the Eddington limit. As the accretion rate increases above Eddington, the spin-down effect becomes progressively stronger, ultimately converging on previous predictions based on non-radiative simulations. In particular, spin evolution drives highly super-Eddington systems toward a black hole spin near zero. The formulae developed in this letter may be applied to galaxy and cosmological scale simulations that include black holes. If magnetically arrested disk accretion is common among supermassive black holes, the present results have broad implications for active galactic nucleus feedback and cosmological spin evolution.

In this work, we study the correlation between Quasi-periodic Oscillation (QPO) frequency and the spectral parameters during various X-ray states in the black hole binary GRS 1915+105 which matches well with the predicted relativistic dynamic frequency (i.e. the inverse of the sound crossing time) at the truncated radii. We have used broadband data of LAXPC and SXT instruments onboard AstroSat. Spectral fitting shows that the accretion rate varies from $\sim 0.1$ to $\sim 5.0 \times 10^{18}$ gm/s and the truncated radius changing from the last stable orbit of an almost maximally spinning black hole, $\sim$ 1.2 to $\sim$ 19 Gravitational radii. For this wide range, the frequencies of the C-type QPO (2 - 6 Hz) follow the trend predicted by the relativistic dynamical frequency model and interestingly, the high-frequency QPO at $\sim$ 70 Hz also follows the same trend, suggesting they originate from the innermost stable circular orbit with the same mechanism as the more commonly observed C-type QPO. While the qualitative trend is as predicted, there are quantitative deviations between the data and the theory, and the possible reasons for these deviations are discussed.

Amanda Pagul, F. Javier Sánchez, Iary Davidzon, Anton M. Koekemoer, Hakim Atek, Renyue Cen, Lukas J. Furtak, Mathilde Jauzac, Guillaume Mahler, Bahram Mobasher, Mireia Montes, Mario Nonino, Keren Sharon, Charles L. Steinhardt, John R. Weaver

16 pages, to be submitted to ApJS

This manuscript presents new astronomical source catalogs using data from the BUFFALO Survey. These catalogs contain detailed information for over 100,000 astronomical sources in the 6 BUFFALO clusters: Abell 370, Abell 2744, Abell S1063, MACS 0416, MACS 0717, and MACS 1149 spanning a total 240 arcmin^2. The catalogs include positions and forced photometry measurements of these objects in the F275W, F336W, F435W, F606W, F814W, F105W, F125W, F140W, and F160W HST -bands, Keck-NIRC2/VLT-HAWKI Ks band, and IRAC Channel 1 and 2 bands. Additionally, we include photometry measurements in the F475W, F625W, and F110W bands for Abell 370. This catalog also includes photometric redshift estimates computed via template fitting using LePhare. When comparing to spectroscopic reference, we obtain an outlier fraction of 9.2% and scatter, normalized median absolute deviation (NMAD), of 0.062. The catalogs are publicly available for their use by the community.

Francesco Benetti, Andrea Lapi, Giovanni Gandolfi, Balakrishna S. Haridasu, Luigi Danese

17 pages, 8 figures, Accepted by Universe

[abridged] Recently, in Benetti et al. (Astrophys. J. 2023, 949, 65), we suggested that the dark matter (DM) component in galaxies may originate fractional gravity. In such a framework, the DM component exists, but the gravitational potential associated to its density distribution is determined by a modified Poisson equation including fractional derivatives, which are meant to describe nonlocal effects. In Benetti et al., we showed that fractional gravity worked very well for reproducing the kinematics of disk-dominated galaxies, especially dwarfs; there is also preliminary evidence that the strength of fractional effects tends to weaken toward more massive systems. Here, we aim to test fractional gravity in galaxy clusters, with a twofold aim: (i) perform an independent sanity check that it can accurately describe such large and massive structures; (ii) derive a clear-cut trend for its strength in systems with different DM masses. To this purpose, we forward model the density and pressure distributions of the intracluster medium (ICM), working out the hydrostatic equilibrium equation in fractional gravity. Then, we perform a Bayesian analysis of the X-COP galaxy cluster sample and infer constraints on the fractional gravity parameters, for individual clusters as well as stacked clusters. We find that fractional gravity performs remarkably well in modeling the ICM profiles for the X-COP sample. We also confirm the weakening of the fractional gravity effects toward more massive systems and derive the overall scaling of the fractional gravity parameters from dwarf galaxies to massive clusters, spanning six orders of magnitude in DM mass. Such an overall trend implies that fractional gravity can substantially alleviate the small-scale issues of the standard DM paradigm, while remaining successful on large cosmological scales.

Strong evidence for the Helling-Downs correlations have been reported by several pulsar timing array collaborations in middle 2023. In this work, we study the state-of-the-art graviton mass bounds by analyzing the observational data of overlap ruduction functions from NANOGrav 15-year data release and CPTA first data release. The data analysis places upper limits on the graviton mass at 95\% confidence level, namely, $m_{g}\lesssim0.43\times10^{-23}\mathrm{eV}$ for NANOGrav and $m_{g}\lesssim0.57\times10^{-23}\mathrm{eV}$ for CPTA. In addition, we discuss implications of these results for scenarios of ultralight tensor dark matter.

Using data from orbits one to eleven of the Parker Solar Probe (PSP) mission, the magnetic field spectral index was measured across a range of heliocentric distances. The previously observed transition between a value of $-5/3$ far from the Sun and a value of $-3/2$ close to the Sun was recovered, with the transition occurring at around $50 \, R_{\odot}$ and the index saturating at $-3/2$ as the Sun is approached. A statistical analysis was performed to separate the variation of the index on distance from its dependence on other parameters of the solar wind that are plausibly responsible for the transition; including the cross helicity, residual energy, turbulence age and the magnitude of magnetic fluctuations. Of all parameters considered the cross helicity was found to be by far the strongest candidate for the underlying variable responsible. The velocity spectral index was also measured and found to be consistent with $-3/2$ over the range of values of cross helicity measured. Possible explanations for the behaviour of the indices are discussed, including the theorised different behaviour of imbalanced, compared to balanced, turbulence.

Giovanni Cabass, Marko Simonović, Matias Zaldarriaga

33 pages, 1 figure

We study how well perturbative forward modeling can constrain cosmological parameters compared to conventional analyses. We exploit the fact that in perturbation theory the field-level posterior can be computed analytically in the limit of small noise. In the idealized case where the only relevant parameter for the nonlinear evolution is the nonlinear scale, we argue that information content in this posterior is the same as in the $n$-point correlation functions computed at the same perturbative order. In the real universe other parameters can be important, and there are possibly enhanced effects due to nonlinear interactions of long and short wavelength fluctuations that can either degrade the signal or increase covariance matrices. We identify several different parameters that control these enhancements and show that for some shapes of the linear power spectrum they can be large. This leads to degradation of constraints in the standard analyses, even though the effects are not dramatic for a $\Lambda$CDM-like cosmology. The aforementioned long-short couplings do not affect the field-level inference which remains optimal. Finally, we show how in these examples calculation of the perturbative posterior motivates new estimators that are easier to implement in practice than the full forward modelling but lead to nearly optimal constraints on cosmological parameters.

Jéferson A. S. Fortunato, Wiliam S. Hipólito-Ricaldi, Marcelo V. dos Santos

Fast Radio Bursts (FRBs) are millisecond-duration pulses occurring at cosmological distances that have emerged as prominent cosmological probes due to their dispersion measure (DM) evolution with redshift. In this work, we use cosmography, a model-independent approach to describe the evolution of the universe, to introduce the cosmographic expansion of the DM-z relation. By fitting two different models for the intergalactic medium and host contributions to a sample of 23 well-localized FRBs, we estimate the kinematic parameters $q_0=-0.59 \substack{+0.20 \\ -0.17}$, $j_0=1.08 \substack{+0.62 \\ -0.56}$, $s_0=-2.1\pm7.0$, and $H_0=69.4\pm4.7$ achieving a precision of $6\%$ and $7\%$ for the Hubble constant depending on the models used for contributions. Furthermore, we demonstrate that this approach can be used to address the long-standing "Missing Baryons" problem in astrophysics by estimating that $82\%$ of the baryonic content of the universe resides in the intergalactic medium, within $7\%$ and $8\%$ precision, according to the contribution models considered here. Our findings highlight the potential of FRBs as a valuable tool in cosmological research and underscore the importance of ongoing efforts to improve our understanding of these enigmatic events.

Kaile Wang, Andreia Carrillo, Melissa K. Ness, Tobias Buck

12 pages, 9 figures. Submitted to MNRAS

Individual abundances in the Milky Way disc record stellar birth properties (e.g. age, birth radius ($R_{\rm birth}$)) and capture the diversity of the star-forming environments over time. Assuming an analytical relationship between ([Fe/H], [$\alpha$/Fe]) and $R_{\rm birth}$, we examine the distributions of individual abundances [X/Fe] of elements C, O, Mg, Si, Ca ($\alpha$), Al (odd-z), Mn (iron-peak), Y, and Ba (neutron-capture) for stars in the Milky Way. We want to understand how these elements might differentiate environments across the disc. We assign tracks of $R_{\rm birth}$ in the [$\alpha$/Fe] vs. [Fe/H] plane as informed by expectations from simulations for $\sim 59,000$ GALAH stars in the solar neighborhood ($R\sim7-9$ kpc) which also have inferred ages. Our formalism for $R_{\rm birth}$ shows that older stars ($\sim$10 Gyrs) have a $R_{\rm birth}$ distribution with smaller mean values (i.e., $\bar{R}_{\mbox{birth}}$$\sim5\pm0.8$ kpc) compared to younger stars ($\sim6$ Gyrs; $\bar{R}_{\mbox{birth}}$$\sim10\pm1.5$ kpc), for a given [Fe/H], consistent with inside-out growth. The $\alpha$-, odd-z, and iron-peak element abundances decrease as a function of $R_{\rm birth}$, whereas the neutron-capture abundances increase. The $R_{\rm birth}$-[Fe/H] gradient we measure is steeper compared to the present-day gradient (-0.067 dex/kpc vs -0.058 dex/kpc), which we also find true for $R_{\rm birth}$-[X/Fe] gradients. These results (i) showcase the feasibility of relating the birth radius of stars to their element abundances, (ii) the abundance gradients across $R_{\rm birth}$ are steeper than those over current radius, and (iii) offer an observational comparison to expectations on element abundance distributions from hydrodynamical simulations.

Lilly Pyras, Christian Glaser, Steffen Hallmann, Anna Nelles

19 pages, 14 figures plus references and appendix

Experiments seeking to detect radio emission stemming from neutrino interactions will soon reach sensitivities that bring a detection within reach. Since experiments like RNO-G or the future IceCube-Gen2 target more than an order of magnitude more effective volume than existing experiments, the renewed and detailed study of rare backgrounds is needed. In this paper, we study the potential background from energy losses of highly energetic atmospheric muons. Due to both limited experimental measurements and limited modeling in hadronic interaction models, the expected event rate is subject to large uncertainties. Here, we estimate rate predictions and their uncertainties for different models and instrumental parameters. We also study possible routes towards mitigation, such as parent air shower detection, and illustrate what is needed to make the first measurement of the prompt muon flux at energies above 10 PeV.

Nicole R. Crumpler (William H. Miller III Department of Physics and Astronomy, Johns Hopkins University Baltimore, MD USA)

15 pages, 3 figures, 2 tables

We constrain a broad class of "hairy" black hole models capable of directly sourcing electromagnetic radiation during a binary black hole merger. This signal is generic and model-independent since it is characterized by the black hole mass ($M$) and the fraction of that mass released as radiation ($\epsilon$). For field energy densities surpassing the Schwinger limit, this mechanism triggers pair-production to produce a gamma-ray burst. By cross-referencing gravitational wave events with gamma-ray observations, we place upper bounds of $\epsilon<10^{-5}-10^{-4}$ for $10-50$ $M_\odot$ black holes depending on the black hole mass. We discuss the weak detection of a gamma-ray burst following GW150914 and show that this event is consistent with rapid electromagnetic emission directly from a "hairy" black hole with $\epsilon\sim10^{-7}-10^{-6}$. Below the Schwinger limit, ambient charged particles are rapidly accelerated to nearly the speed of light by the strong electromagnetic field. For 1-50 $M_\odot$ black holes and $\epsilon$ ranging from $10^{-20}$ to $10^{-7}$, the typical proton energies are $\sim20$ GeV-20 TeV and electron energies are $\sim0.01-10$ GeV. At these energies, cosmic ray protons and electrons quickly diffuse into the Milky Way's background magnetic field, making it difficult to identify a point source producing them. Overall, constraining $\epsilon$ in this less energetic regime becomes difficult and future constraints may need to consider specific models of "hairy" black holes.

Peak sidelobe level reduction (PSLR) is crucial in the application of large-scale array antenna, which directly determines the radiation performance of array antenna. We study the PSLR of subarray level aperiodic arrays and propose three array structures: dislocated subarrays with uniform elements (DSUE), uniform subarrays with random elements (USRE), dislocated subarrays with random elements (DSRE). To optimize the dislocation position of subarrays and random position of elements, the improved Bat algorithm (IBA) is applied. To draw the comparison of PSLR effect among these three array structures, we take three size of array antennas from small to large as examples to simulate and calculate the redundancy and peak sidelobe level (PSLL) of them. The results show that DSRE is the optimal array structure by analyzing the dislocation distance of subarray, scanning angle and applicable frequency. The proposed design method is a universal and scalable method, which is of great application value to the design of large-scale aperiodic array antenna.

We investigate the nuclear symmetry energy and neutron star properties using a Bayesian analysis based on constraints from different chiral effective field theory calculations using new energy density functionals that allow for large variations at high densities. Constraints at high densities are included from observations of GW170817 and NICER. In particular, we show that both NICER analyses lead to very similar posterior results for the symmetry energy and neutron star properties when folded into our equation of state framework. Using the posteriors, we provide results for the symmetry energy and the slope parameter, as well as for the proton fraction, the speed of sound, and the central density in neutron stars. Moreover, we explore correlations of neutron star radii with the pressure and the speed of sound in neutron stars. Our 95\% credibility ranges for the symmetry energy $S_v$, the slope parameter $L$, and the radius of a 1.4\,$M_\odot$ neutron star $R_{1.4}$ are $S_v=(30.6-33.9)$\,MeV, $L=(43.7-70.0)$\,MeV, and $R_{1.4}=(11.6-13.2)$\,km. Our analysis for the proton fraction shows that larger and-or heavier neutron stars are more likely to cool rapidly via the direct Urca process. Within our equation of state framework a maximum mass of neutron stars $M_{\rm max}>2.1\,M_\odot$ indicates that the speed of sound needs to exceed the conformal limit.

We perform a Bayesian analysis of NANOGrav 15yr and IPTA DR2 pulsar timing residuals and show that the recently detected stochastic gravitational-wave background (SGWB) is compatible with a SGWB produced by bubble dynamics during a cosmological first-order phase transition. The timing data suggests that the phase transition would occur around QCD confinement temperature and would have a slow rate of completion. This scenario can naturally lead to the abundant production of primordial black holes (PBHs) with solar masses. These PBHs can potentially be detected by current and advanced gravitational wave detectors LIGO-Virgo-Kagra, Einstein Telescope, Cosmic Explorer, by astrometry with GAIA and by 21-cm survey.

The stochastic signal detected by NANOGrav, PPTA, EPTA, and CPTA can be explained by the scalar-induced gravitational waves. In order to determine the scalar-induced gravitational waves model that best fits the stochastic signal, we employ both single- and double-peak parameterizations for the power spectrum of the primordial curvature perturbations, where the single-peak scenarios include the $\delta$-function, box, lognormal, and broken power law model, and the double-peak scenario is described by the double lognormal form. Using Bayesian inference, we find that there is no significant evidence for or against the single-peak scenario over the double-peak model, with $\log$ (Bayes factors) among these models $\ln \mathcal{B} < 1$. Therefore, we are not able to distinguish the different shapes of the power spectrum of the primordial curvature perturbation with the current sensitivity of pulsar timing arrays.

From prolonged X-ray and optical data of the ultra-compact binary HM Cancri, two groups recently measured the second derivative of its orbital frequency. The space gravitational wave (GW) detector LISA will detect $\sim10^4$ Galactic binaries and their second frequency derivatives will be interesting observational targets for LISA. Here, we forecast the GW signal analysis for HM Cancri, as an ideal reference system for these numerous binaries. We find that, in its nominal operation period $T\sim4$yr, LISA is unlikely to realize a sufficient measurement precision for the reported second frequency derivative of this binary. However, because of a strong dependence on the time baseline, the precision will be drastically improved by extending the operation period of LISA or combining it with other missions (e.g., Taiji and TianQin) in a sequential order.

We investigate the effects of a Non-Standard Interaction (NSI) extension of the standard model of particle physics on solar neutrino flavour oscillations. This NSI model introduces a $U_{Z^\prime}(1)$ gauge symmetry through a $Z^\prime$ boson that mixes with the photon, creating a neutral current between active neutrinos and matter fields via a unique coupling to up and down quarks. The interaction is defined by a single parameter, $\zeta_o$, which is related to the $Z^\prime$ boson's mass $m_{Z^\prime}$ and coupling constant $g_{Z^\prime}$. Notably, this model relaxes the bounds on Coherent Elastic Neutrino-Nucleus Scattering experiments and fits the experimental values of the anomalous magnetic dipole moment of the muon. In this study, we use solar neutrino measurements and an up-to-date standard solar model to evaluate the neutrino flavour oscillations and assess the constraints on $\zeta_o$. Our study indicates that the NSI model aligns with the current solar neutrino data when $\zeta_o$ is between $-0.7$ and $0.002$. These models have $\chi^2_{\nu}$ values equal to or better than the standard neutrino flavor oscillation model, which stands at a $\chi^2_{\nu}$ of 3.12. The best NSI model comes with a $\zeta_o$ value of -0.2 and a $\chi^2_{\nu}$ of 2.96. Including extra data from the Darwin experiment in our analysis refines the range of $\zeta_o$ values from $-0.7$ to $0.002$, down to $-0.5$ to $-0.002$. These results hint at the possible existence of novel interactions, given that NSI models achieve a comparable or superior fit to the solar neutrino data when contrasted with the prevailing standard model of neutrino flavour oscillation.

G. Fanizza, M. Gasperini, G. Marozzi

7 pages, no figures. Published in the Special Issue on "Universe: Feature Papers 2023 -- Cosmology", edited by K. Bamba (Universe, MDPI, 2023). Available at: this https URL

In this article, we compare different formulations of the number count prescription using the convenient formalism of the Geodesic-Light-Cone gauge. We then find a simple, exact, and very general expression of such a prescription which is suitable for generalised applications.

Federico Schianchi, Henrique Gieg, Vsevolod Nedora, Anna Neuweiler, Maximiliano Ujevic, Mattia Bulla, Tim Dietrich

29 pages, 27 figures

Neutrino interactions are essential for an accurate understanding of the binary neutron star merger process. In this article, we extend the code infrastructure of the well-established numerical-relativity code BAM that until recently neglected neutrino-driven interactions. In fact, while previous work allowed already the usage of nuclear-tabulated equations of state and employing a neutrino leakage scheme, we are moving forward by implementing a first-order multipolar radiation transport scheme (M1) for the advection of neutrinos. After testing our implementation on a set of standard scenarios, we apply it to the evolution of four low-mass binary systems, and we perform an analysis of ejecta properties. We also show that our new ejecta analysis infrastructure is able to provide numerical relativity-informed inputs for the codes $\texttt{POSSIS}$ and $\texttt{Skynet}$, for the computation of kilonova lightcurves and nucleosynthesis yields, respectively.

E. Babichev, D. Gorbunov, S. Ramazanov, R. Samanta, A. Vikman

13 pages, 1 figure

We discuss cosmic domain walls described by a tension red-shifting with the expansion of the Universe. These melting domain walls emit gravitational waves (GW) with the low-frequency spectral shape $\Omega_{gw}\propto f^{2}$ corresponding to the spectral index $\gamma=3$ favoured by the recent NANOGrav 15 yrs data. We discuss a concrete high-energy physics scenario proposed in Refs. [1,2] which leads to such a melting domain wall network in the early Universe. This scenario involves a feebly coupled scalar field $\chi$, which can serve as a promising dark matter candidate. We identify parameters of the model matching the GW characteristics observed in the NANOGrav data. The dark matter mass is pushed to the ultra-light range below $10^{-11}-10^{-12}\,\text{eV}$ which is accessible through planned observations thanks to the effects of superradiance of rotating black holes.

Henri-Claude Nataf, Nathanaël Schaeffer

Supplementary Material available at: this https URL with python programs at: this https URL

Planetary cores are the seat of rich and complex fluid dynamics, in which the effects of rotation and magnetic field combine. The equilibria governing the strength of the magnetic field produced by the dynamo effect, the organisation and amplitude of the flow, and those of the density field, remain debated despite remarkable progress made in their numerical simulation. This paper proposes a new approach based on the explicit consideration of the variation of time scales $\tau$ with spatial scales $\ell$ for the different physical phenomena involved. The $\tau$-$\ell$ diagrams thus constructed constitute a very complete graphic summary of the dynamics of the object under study. They reveal the domains of validity of the different possible dynamic regimes. Several scenarios are thus constructed and discussed for the Earth's core, shedding new light on the width of convective columns and on the force equilibria to be considered. A QG-MAC scenario adapted from Aubert [2019] gives a good account of the observations. A diversion to Venus reveals the subtlety and relativity of the notion of 'fast rotator'. A complete toolbox is provided, allowing everyone to construct a $\tau$-$\ell$ diagram of a numerical simulation, a laboratory experiment, a theory, or a natural object. Supplementary material for this article is supplied as a separate archive Nataf_Schaeffer_SupMat.zip, the related data is displayed in document Nataf_Schaeffer_SupMat.pdf.

Anirban Chatterjee, Saddam Hussain, Kaushik Bhattacharya

4 pages and 4 figures; Conference proceeding for XXV DAE-BRNS HEP Symposium 2022

We examine the scenario of non-minimally coupled relativistic fluid and $k$-essence scalar field in a flat Friedmann-Lemaitre-Robertson-Walker universe. By adding a non-minimal coupling term in the Lagrangian level, we study the variation of Lagrangian with respect to independent variables, which produces modified scalar field and Friedmann equations. Using dynamical stability approach in different types of interaction models with two types of scalar field potential, we explore this coupled framework. Implementing detailed analysis, we can conclude our models can able to produce stable late-time cosmic acceleration.

Oscar López-Aguayo, J. C. López-Domínguez, M. Sabido

9 pages, 2 figures

In this work we study the effects of the generalized uncertainty principle (GUP) in cosmology. We start with the Friedmann-Robertson-Walker (FRW) model endowed with a scalar field. After introducing the GUP modification to the model, we solve for the quantum and classical cases. Finally we find the GUP modified Friedmann equations.

Wilfried Buchmuller, Valerie Domcke, Kai Schmitz

28 pages, 3 figures

Many symmetry breaking patterns in grand unified theories (GUTs) give rise to cosmic strings that eventually decay when pairs of GUT monopoles spontaneously nucleate along the string cores. These strings are known as metastable cosmic strings and have intriguing implications for particle physics and cosmology. In this article, we discuss the current status of metastable cosmic strings, with a focus on possible GUT embeddings and connections to inflation, neutrinos, and gravitational waves (GWs). The GW signal emitted by a network of metastable cosmic strings in the early universe differs, in particular, from the signal emitted by topologically stable strings by a suppression at low frequencies. Therefore, if the underlying symmetry breaking scale is close to the GUT scale, the resulting GW spectrum can be accessible at current ground-based interferometers as well as at future space-based interferometers, such as LISA, and at the same time account for the signal in the most recent pulsar timing data sets. Metastable cosmic strings thus nourish the hope that future GW observations might shed light on fundamental physics close to the GUT scale.

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.

We report an evidence that self-similar and anomalous scalings coexist in a turbulent medium, particularly in fluctuations of the magnetic field flux density in magnetized plasma of the solar photosphere. The structure function scaling exponents in the inertial range have been analyzed for fluctuations grouped according to the sign of the path-dependent stochastic entropy production. It is found that the scaling exponents for fluctuations with the positive entropy production follow the phenological linear dependence for the magnetohydrodynamic turbulence. For fluctuations with the negative entropy production, the scaling is anomalous.

Jiahang Hu, Bingrong Yu, Shun Zhou

34 pages, 7 figures, 1 table

The Higgs triplet model (HTM) extends the Standard Model (SM) by one complex triplet scalar (also known as the type-II seesaw model), offering a simple and viable way to account for nonzero neutrino masses. On the other hand, the nontrivial couplings of the triplet to the gauge fields and to the SM Higgs field are expected to influence the topological vacuum structure of the SM, and consequently, the energy and the field configuration of the electroweak sphaleron. The sphaleron process plays a crucial role in dynamically generating the baryon asymmetry of the Universe. In this work, we study the vacuum structure of the gauge and Higgs fields and calculate the saddle-point sphaleron configuration in the HTM. The coupled nonlinear equations of motion of the sphaleron are solved using the spectral method. We find the inclusion of the triplet scalar could in principle significantly change the sphaleron energy compared with the SM. Nevertheless, at zero temperature, the current stringent experimental constraint on the vacuum expectation value of the triplet suppresses the difference. Interestingly, we find that there still exists some narrow parameter space where the sphaleron energy can be enhanced up to $30\%$ compared with the SM case.

We propose an inflationary scenario where the inflaton field is non-minimally coupled to spacetime curvature and inflation is driven by a vacuum energy symmetry breaking potential without specifying \emph{a priori} whether the inflaton field is small or large. As we incorporate vacuum energy into our analysis, we further explore the implications of a non-zero potential offset in relation to the emergence of inflationary dynamics. Thus, we propose that vacuum energy can transform into particles as a result of the transition triggered by spontaneous symmetry breaking. This entails a vacuum energy cancellation that yields an effective cosmological constant during inflation by virtue of a quasi-de Sitter evolution and shows that the vacuum energy contribution can manifest as \emph{geometric particles} produced by inflaton fluctuations, with particular emphasis on super-Hubble modes. We conjecture these particles as \emph{quasi-particles} arising from interaction between the inflaton and spacetime geometry, enhanced by non-minimal coupling. Specifically, we propose that dark matter arises from a pure geometric quasi-particle contribution, and we quantify the corresponding dark matter candidate ranges of mass. In this scenario, we further find that a zero potential offset leads to a bare cosmological constant at the end of inflation, while a negative offset would require an additional kinetic (or potential) contribution in order to be fully-canceled. In this regard, we conclude that the scenario of large field inflaton is preferred since it necessitates a more appropriate selection of the offset. Our conclusion is reinforced as small field inflaton would lead to a significant screening of the Newtonian gravitational constant as inflation ends.