Astrophysical Models Research Articles

Rate coefficients for C and O2 reactive collisions relevant to interstellar clouds from QCT and machine learning.

The chemical reactions between certain interstellar molecules are exothermic in nature and barrierless in the entrance channel, allowing these reactions to occur rapidly even at low astronomical temperatures, e.g., C and O2 interaction. Obtaining detailed rovibrational transition parameters for the reaction between C and O2, such as state-selected rate coefficients, is crucial for studying the associated atmospheric and astronomical environments. Hence, this work presents an approach that combines quasi-classical trajectory calculations with machine learning techniques based on Neural Network (NN) and Gaussian Process Regression (GPR) to determine state-selected rate coefficients. Employing this approach, we significantly reduced the computational requirements while simultaneously obtaining the accurate state-selected reaction cross sections and rate coefficients for the collision of C and O2. Both the NN-based and GPR-based models established in this work accurately predict the results calculated from explicit numerical calculations in the explored temperature range of 50-1500K, achieving a coefficient of determination R2 > 0.96. Most importantly, the current work provides the most comprehensive dataset of rovibrational rate coefficients of v = 0-4, j = 0-70 → v' = 0-15 for the astrophysical modeling of the C-O2 collision system.

Measurement of the Ep=416.9\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$E_\ extrm{p}=416.9$$\\end{document} keV resonance strength in the 29\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{29}$$\\end{document}Si(p,γ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\gamma $$\\end{document})30\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{30}$$\\end{document}P reaction

Silicon isotopic ratios measured in meteoritic presolar grains can provide useful information about the nucleosynthesis origin of these isotopes if the rates of nuclear reactions responsible for their production are known. One of the key reactions determining the Si isotopic abundances is 29\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{29}$$\\end{document}Si(p,γ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\gamma $$\\end{document})30\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{30}$$\\end{document}P. Its reaction rate is not known with sufficient precision due in part to some ambiguous resonance strength values. In the present work, the strength of the Ep\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$E_\ extrm{p}$$\\end{document} = 416.9 keV resonance has been measured with high precision using the activation technique. The new strength of ωγ=219±\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\omega \\gamma = 219 \\pm $$\\end{document} 16 meV can be used in updated reaction rate estimations and astrophysical models.

Inelastic Triatom-Atom Quantum Close-Coupling Dynamics in Full Dimensionality: All Rovibrational Mode Quenching of Water Due to the H Impact on a Six-Dimensional Potential Energy Surface.

The rovibrational level populations, and subsequent emission in various astrophysical environments, are driven by inelastic collision processes. The available rovibrational rate coefficients for water have been calculated using a number of approximations. We present a numerically exact calculation for the rovibrational quenching for all water vibrational modes due to collisions with atomic hydrogen. The scattering theory implements a quantum close-coupling (CC) method on a high level ab initio six-dimensional (6D) potential energy surface (PES). Total rovibrational quenching cross sections for excited bending levels were compared with earlier results on a 4D PES with the rigid-bender close-coupling (RBCC) approximation. General agreement between 6D-CC and 4D-RBCC calculations are found, but differences are evident including the energy and amplitude of low-energy orbiting resonances. Quenching cross sections from the symmetric and asymmetric stretch modes are provided for the first time. The current 6D-CC calculation provides accurate inelastic data needed for astrophysical modeling.

VENICE: A multi-scale operator-splitting algorithm for multi-physics simulations

Context. We present VENICE, an operator-splitting algorithm to integrate a numerical model on a hierarchy of timescales. Aims. VENICE allows a wide variety of different physical processes operating on different scales to be coupled on individual and adaptive time-steps. It therewith mediates the development of complex multi-scale and multi-physics simulation environments with a wide variety of independent components. Methods. The coupling between various physical models and scales is dynamic, and realised through (Strang) operators splitting using adaptive time-steps. Results. We demonstrate the functionality and performance of this algorithm using astrophysical models of a stellar cluster, first coupling gravitational dynamics and stellar evolution, then coupling internal gravitational dynamics with dynamics within a galactic background potential, and finally combining these models while also introducing dwarf galaxy-like perturbers. These tests show numerical convergence for decreasing coupling timescales, demonstrate how VENICE can improve the performance of a simulation by shortening coupling timescales when appropriate, and provide a case study of how VENICE can be used to gradually build up and tune a complex multi-physics model. Although the examples provided here couple dedicated numerical models, VENICE can also be used to efficiently solve systems of stiff differential equations.

Comparing Performance in the Codes AREPO, GIZMO and GAMER-2

Modern astrophysical simulation is the main aspect of human getting insight from the cosmos. Numerical simulations became irreplaceable tools for generating various phases of the universe, the formation of galaxies, and other astrophysical phenomena. As the advances of computational power, validity and diversity of astrophysical modeling have been growing. In addition to the development of the processor performance, the diversity of the software community is also rising therefore, research between various simulation codes is worthwhile. It allowed researchers to know the nature of different simulation techniques linked to their basic algorithms. This study summarizes the basics of common simulation methods used by researchers including N-body Simulations, Smoothed Particle Hydrodynamics, Monte Carlo Methods, Finite Difference and Finite Element Methods. Furthermore, this research analyzes the basics of AREPO, GIZMO, and GAMER-2 codes, and compared the performance of these codes. Last, according to the analysis, this study discusses and demonstrates the advantages and limitations of these codes. These results shed light on guiding further exploration of cosmology simulations.

AM3: An Open-source Tool for Time-dependent Lepto-hadronic Modeling of Astrophysical Sources

We present the Astrophysical Multimessenger Modeling (AM 3 ) software. AM 3 is a documented open-source software (source code at https://gitlab.desy.de/am3/am3; user guide and documentation at https://am3.readthedocs.io/en/latest/) that efficiently solves the coupled integro-differential equations describing the temporal evolution of the spectral densities of particles interacting in astrophysical environments, including photons, electrons, positrons, protons, neutrons, pions, muons, and neutrinos. The software has been extensively used to simulate the multiwavelength and neutrino emission from active galactic nuclei (including blazars), gamma-ray bursts, and tidal disruption events. The simulations include all relevant nonthermal processes, namely synchrotron emission, inverse Compton scattering, photon–photon annihilation, proton–proton and proton–photon pion production, and photo-pair production. The software self-consistently calculates the full cascade of primary and secondary particles, including nonlinear feedback processes and predictions in the time domain. It also allows the user to track separately the particle densities produced by means of each distinct interaction process, including the different hadronic channels. With its efficient hybrid solver combining analytical and numerical techniques, AM 3 combines efficiency and accuracy at a user-adjustable level. We describe the technical details of the numerical framework and present three examples of applications to different astrophysical environments.

Time delay of light in the gravitational lensing of supermassive black holes in dark matter halos

AbstractThe dark matter halo has non-negligible effects on the gravitational lensing of supermassive black hole in the galaxy center. Our work presents a study on the time-delay of light in gravitational lensing of black holes enclosed by dark matter halos. To provide a precise description on the distribution of dark matter in galaxies, we choose several famous phenomenological dark matter halo models in astrophysics, including the NFW, Beta, Burkert and Moore models, to carry out the present study. Through numerically calculating the time-delay of light in gravitational lensing, a comparative analysis of the dark matter effects within different halo models has been performed. Assuming typical length scales associated with the galactic gravitational lensing, numerical results indicate that the NFW, Beta, Burkert and Moore dark matter halos can significantly enhance the time delay of light in gravitational lenisng of central supermassive black holes. The enhancing effect becomes more pronounced with a small dark matter halo scale and an increasing dark matter halo mass.

Gravitational lens system as a long baseline detector of extremely low frequency primordial gravitational wave

ABSTRACT The effect of extremely low frequency primordial gravitational wave with arbitrary direction of propagation on a gravitational lens system in expanding universe is investigated. From the point of view of real astrophysical lens model, singular isothermal sphere lens model is adopted in the gravitational lens system. The results show that, under the perturbation from extremely low frequency primordial gravitational wave, time delay in the gravitational lens system is very sensitive to extremely low frequency primordial gravitational wave and could strongly deviate from that deduced from theoretical model. This means that the strongly deviated time delay could be the imprint of extremely low frequency primordial gravitational wave on gravitational lens system, indicating that gravitational lens system could be used as a long baseline detector to detect extremely low frequency primordial gravitational wave.

Resource Letter: Synthesis of the elements in stars

This Resource Letter provides a guide to the literature on the field of nuclear astrophysics, and particularly the origin of the elements. Nuclear astrophysics is a multidisciplinary field that aims at understanding where everything we see around us comes from and how it came to be. Astronomical observations, astrophysics modeling, and nuclear physics experiment and theory come together to answer important questions like: Where and how are the elements created? How do stars evolve? What drives the different types of stellar explosions? What is left behind after the cataclysmic death of a star? This Resource Letter presents our current understanding of the origin of the various chemical elements, together with modern research and new developments in the field, with a particular focus on the measurement of nuclear properties for astrophysical applications.

Gravitational lensing of spherically symmetric black holes in dark matter halos

The gravitational lensing of supermassive black holes surrounded by dark matter halo has attracted a great number of interests in recent years. However, many studies employed simplified dark matter density models, which makes it very hard to give a precise prediction on the dark matter effects in real astrophysical galaxies. In this work, to more accurately describe the distribution of dark matter in real astrophysical galaxies, we study the gravitational lensing of black holes in astrophysical dark matter halo models (Beta, Burkert, Brownstein, and Moore). The deflection angle is obtained using a generalized Gibbons-Werner approach. The visual angular positions and the Einstein rings are also calculated by adopting the gravitational lens equation. Specifically, we choose the supermassive black holes in Milky Way Galaxy, Andromeda galaxy (M31), Virgo galaxy (M87), and ESO138-G014 galaxy as examples, including the corresponding fitted value of dark matter halos. The results suggest that the dark matter halo described by the Beta model has non-negligible influences on the gravitational deflection angle and gravitational lensing observations. However, the Burkert, Brownstein, and Moore models have relatively small influences on angular position of images and the Einstein ring.

Cosmic-ray propagation models elucidate the prospects for antinuclei detection

Tentative observations of cosmic-ray antihelium by the AMS-02 collaboration have re-energized the quest to use antinuclei to search for physics beyond the standard model. However, our transition to a data-driven era requires more accurate models of the expected astrophysical antinuclei fluxes. We use a state-of-the-art cosmic-ray propagation model, fit to high-precision antiproton and cosmic-ray nuclei (B, Be, Li) data, to constrain the antinuclei flux from both astrophysical and dark matter annihilation models. We show that astrophysical sources are capable of producing 𝒪(1) antideuteron events and 𝒪(0.1) antihelium-3 events over 15 years of AMS-02 observations. Standard dark matter models could potentially produce higher levels of these antinuclei, but showing a different energy-dependence. Given the uncertainties in these models, dark matter annihilation is still the most promising candidate to explain preliminary AMS-02 results. Meanwhile, any robust detection of antihelium-4 events would require more novel dark matter model building or a new astrophysical production mechanism.

Tetrahedral grids in Monte Carlo radiative transfer

Context. To understand the structures of complex astrophysical objects, 3D numerical simulations of radiative transfer processes are invaluable. For Monte Carlo radiative transfer, the most common radiative transfer method in 3D, the design of a spatial grid is important and non-trivial. Common choices include hierarchical octree and unstructured Voronoi grids, each of which has advantages and limitations. Tetrahedral grids, commonly used in ray-tracing computer graphics, can be an interesting alternative option. Aims. We aim to investigate the possibilities, advantages, and limitations of tetrahedral grids in the context of Monte Carlo radiative transfer. In particular, we want to compare the performance of tetrahedral grids to other commonly used grid structures. Methods. We implemented a tetrahedral grid structure, based on the open-source library TetGen, in the generic Monte Carlo radiative transfer code SKIRT. Tetrahedral grids can be imported from external applications or they can be constructed and adaptively refined within SKIRT. We implemented an efficient grid traversal method based on Plücker coordinates and Plücker products. Results. The correct implementation of the tetrahedral grid construction and the grid traversal algorithm in SKIRT were validated using 2D radiative transfer benchmark problems. Using a simple 3D model, we compared the performance of tetrahedral, octree, and Voronoi grids. With a constant cell count, the octree grid outperforms the tetrahedral and Voronoi grids in terms of traversal speed, whereas the tetrahedral grid is poorer than the other grids in terms of grid quality. All told, we find that the performance of tetrahedral grids is relatively poor compared to octree and Voronoi grids. Conclusions. Although the adaptively constructed tetrahedral grids might not be favourable in most media representative of astrophysical simulation models, they still form an interesting unstructured alternative to Voronoi grids for specific applications. In particular, they might prove useful for radiative transfer post-processing of hydrodynamical simulations run on tetrahedral or unstructured grids.

Exploring the role of the halo mass function for inferring astrophysical parameters during reionisation

Abstract Detecting the 21-cm signal at z ≳ 6 will reveal insights into the properties of the first galaxies responsible for driving reionisation. To extract this information, we perform parameter inference with 3D simulations of the 21-cm signal embedded within a Bayesian inference pipeline. Presently, when performing inference, we must choose which sources of uncertainty to sample and which to hold fixed. Since the astrophysics of galaxies is much more uncertain than that of the underlying halo-mass function (HMF), we typically parameterise and model the former while fixing the latter. However, doing so may bias our inference of the galaxy properties. In this work, we explore the consequences of assuming an incorrect HMF and quantify the relative biases on our inferred astrophysical model parameters when considering the wrong HMF. We then relax this assumption by constructing a generalised five parameter HMF model and simultaneously recover it with our underlying astrophysical model. For this, we use 21CMFAST and perform Simulation-Based Inference using marginal neural ratio estimation to learn the likelihood-to-evidence ratio with Swyft. Using a mock 1000-hour observation of the 21-cm power spectrum from the forthcoming Square Kilometre Array, conservatively assuming foreground wedge avoidance, we find that assuming the incorrect HMF can bias the recovered astrophysical parameters by up to ∼3 − 4σ even when including independent information from observed luminosity functions. Using our generalised HMF model, although we recover our astrophysical parameters with a factor of ∼2 − 4 larger marginalised uncertainties, the constraints are unbiased, agnostic to the underlying HMF and therefore more conservative.

Debiasing with Diffusion: Probabilistic Reconstruction of Dark Matter Fields from Galaxies with CAMELS

Galaxies are biased tracers of the underlying cosmic web, which is dominated by dark matter (DM) components that cannot be directly observed. Galaxy formation simulations can be used to study the relationship between DM density fields and galaxy distributions. However, this relationship can be sensitive to assumptions in cosmology and astrophysical processes embedded in galaxy formation models, which remain uncertain in many aspects. In this work, we develop a diffusion generative model to reconstruct DM fields from galaxies. The diffusion model is trained on the CAMELS simulation suite that contains thousands of state-of-the-art galaxy formation simulations with varying cosmological parameters and subgrid astrophysics. We demonstrate that the diffusion model can predict the unbiased posterior distribution of the underlying DM fields from the given stellar density fields while being able to marginalize over uncertainties in cosmological and astrophysical models. Interestingly, the model generalizes to simulation volumes ≈500 times larger than those it was trained on and across different galaxy formation models. The code for reproducing these results can be found at https://github.com/victoriaono/variational-diffusion-cdm ✎.

Strong observational support for the [formula omitted] timeline in the early universe

JWST’s recent discovery of well-formed galaxies and supermassive black holes only a few hundred Myr after the big bang seriously challenges the timeline predicted by ΛCDM. Now, the latest identification of polycyclic aromatic hydrocarbons (PAHs) at z=6.71, together with these earlier inconsistencies, makes the time compression problem in this model quite overwhelming. We consider the timeline associated with the formation and growth of PAH grains based on current astrophysical models and argue that their appearance at z=6.71 favors the structure formation history in Rh=ct rather than that of Planck-ΛCDM. We estimate the time at which they must have started growing in each case, and then trace their history through various critical events, such as the end of the ‘dark ages’, the beginning of Pop III star formation, and the onset of reionization. Together, these three distinct discoveries by JWST, viz. high-z galaxies, high-z quasars and the surprisingly early appearance of PAHs, all paint a fully consistent picture in which the timeline in ΛCDM is overly compressed at z≳6, while strongly supporting the expansion history in the early Universe predicted by Rh=ct.

Constraining the Clustering and 21 cm Signature of Radio Galaxies at Cosmic Dawn

The efficiency of radio emission is an important unknown parameter of early galaxies at cosmic dawn, as models with high efficiency have been shown to modify the cosmological 21 cm signal substantially, deepening the absorption trough and boosting the 21 cm power spectrum. Such models have been previously directly constrained by the overall extragalactic radio background, as observed by Absolute Radiometer for Cosmology, Astrophysics, and Diffuse Emission 2 and Long Wavelength Array. In this work, we constrain the clustering of high-redshift radio sources by utilizing the observed upper limits on arcminute-scale anisotropy from the Very Large Array at 4.9 GHz and Australia Telescope Compact Array at 8.7 GHz. Using a seminumerical simulation of a plausible astrophysical model for illustration, we show that the clustering constraints on the radio efficiency are much stronger than those from the overall background intensity by a factor that varies between 18 and 55 in the redshift range of 7–22. As a result, the predicted maximum depth of the global 21 cm signal is lowered by a factor of 6 (to 1400 mK), and the maximum 21 cm power spectrum peak at cosmic dawn is lowered by a factor of 45 (to 1.3 × 105 mK2). We conclude that the observed clustering is the strongest current direct constraint on such models, but strong early radio emission from galaxies remains viable for producing a strongly enhanced 21 cm signal from cosmic dawn.

Observational appearances of hairy black holes in the framework of gravitational decoupling

Abstract We investigate the observational appearance of static and spherically symmetric hairy black holes in the framework of gravitational decoupling with the weak energy condition (WEC). Two types of thin illumination conditions are studied: spherical accretion and disk accretion. As the hairy parameter increases, the size of the photon sphere and photon rings in both models decreases, and the overall luminosity attenuation becomes more pronounced. In spherical accretion, the luminosity of infalling accretion is significantly lower than that of stationary accretion. In disk accretion the luminosity of the black hole is contributed by direct emission, the lensing ring and the photon ring. Employing four types of astrophysical disk luminosity model, we investigate the appearance of halos and note that their luminosities do not superimpose when the source is on or beyond the innermost stable circular orbit.

Atomic structure and transition parameters of the C-like ion K XIV

Energy levels, lifetimes, oscillator strengths and transition probabilities for the multicharged carbon like K XIV ion have been calculated with the configuration expansion: 2s22p2, 2p4, 2s22p3p, 2s22p4p, 2s2p23s, 2s2p24s, 2s2p23d, 2s2p24d, 2s2p23p, 2s2p24p, 2s22p3s, 2s22p4s, 2s22p3d, 2s22p4d and 2s2p3. Two methods were used in the calculations: the Hartree-Fock pseudo-relativistic approach and the Thomas-Fermi-Dirac-Amaldi potential approach using the Cowan and the AUTOSTRUCTURE atomic structure codes respectively. Results have been compared with data from National Institute of Standards and Technology Atomic Spectra Database (NIST-ASD) and from other theoretical methods. Obtained new data will be important for plasma analysis in laboratories and for astrophysical modeling.

Calibrating gravitational-wave search algorithms with conformal prediction

In astronomy, we frequently face the decision problem: does this data contain a signal? Typically, a statistical approach is used, which requires a threshold. The choice of threshold presents a common challenge in settings where signals and noise must be delineated, but their distributions overlap. Gravitational-wave astronomy, which has gone from the first discovery to catalogs of hundreds of events in less than a decade, presents a fascinating case study. For signals from colliding compact objects, the field has evolved from a frequentist to a Bayesian methodology. However, the issue of choosing a threshold and validating noise contamination in a catalog persists. Confusion and debate often arise due to the misapplication of statistical concepts, the complicated nature of the detection statistics, and the inclusion of astrophysical background models. We introduce conformal prediction (CP), a framework developed in machine learning to provide distribution-free uncertainty quantification to point predictors. We show that CP can be viewed as an extension of the traditional statistical frameworks whereby thresholds are calibrated such that the uncertainty intervals are statistically rigorous and the error rate can be validated. Moreover, we discuss how CP offers a framework to optimally build a metapipeline combining the outputs from multiple independent searches. We introduce CP with a toy cosmic-ray detector, which captures the salient features of most astrophysical search problems and allows us to demonstrate the features of CP in a simple context. We then apply the approach to a recent gravitational-wave mock data challenge using multiple search algorithms for compact binary coalescence signals in interferometric gravitational-wave data. Finally, we conclude with a discussion on the future potential of the method for gravitational-wave astronomy. Published by the American Physical Society 2024

Abundance, sizes, and major element compositions of components in CR and LL chondrites: Formation from single reservoirs

AbstractAbundances, apparent sizes, and individual chemical compositions of chondrules, refractory inclusions, other objects, and surrounding matrix have been determined for Semarkona (LL3.00) and Renazzo (CR2) using consistent methods and criteria on X‐ray element intensity maps. These represent the non‐carbonaceous (NC, Semarkona) and carbonaceous chondrite (CC, Renazzo) superclans of chondrite types. We compare object and matrix abundances with similar data for CM, CO, K, and CV chondrites. We assess, pixel‐by‐pixel, the major element abundance in each object and in the entire matrix. We determine the abundance of “metallic chondrules” in LL chondrites. Chondrules with high Mg/Si and low Fe/Si and matrix carrying opposing ratios complement each other to make whole rocks with near‐solar major element ratios in Renazzo. Similar Mg/Si and Fe/Si chondrule–matrix relationships are seen in Semarkona, which is within 11% of solar Mg/Si but significantly Fe‐depleted. These results provide a robust constraint in support of single‐reservoir models for chondrule formation and accretion, ruling out whole classes of astrophysical models and constraining processes of chondrite component formation and accretion into chondrite parent bodies.