Search NASA⌕ Search

SEARCH · Search NASA

Results for “astrophysical simulations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

Numerical Relativity Simulations for Black Hole Merger Astrophysics

Massive black hole mergers are perhaps the most energetic astronomical events, establishing their importance as gravitational wave sources for LISA, and also possibly leading to observable influences on their local environments. Advances in numerical relativity over the last five years have fueled the development of a rich physical understanding of general relativity's predictions for these events. Z will overview the understanding of these event emerging from numerical simulation studies. These simulations elucidate the pre-merger dynamics of the black hole binaries, the consequent gravitational waveform signatures ' and the resulting state, including its kick velocity, for the final black hole produced by the merger. Scenarios are now being considered for observing each of these aspects of the merger, involving both gravitational-wave and electromagnetic astronomy.

Baker, John G.↗

Nucleobases and Other Prebiotic Species from the UV Irradiation of Pyrimidine in Astrophysical Ices

Nucleobases are aromatic N-heterocycles that constitute the informational subunits of DNA and RNA and are divided into two families: pyrimidine bases (uracil, cytosine, and thymine) and purine bases (adenine and guanine). Nucleobases have been detected in meteorites and their extraterrestrial origin confirmed by isotope measurement. Although no N-heterocycles have been individually identified in the ISM, the 6.2-micron interstellar emission feature seen towards many astronomical objects suggests a population of such molecules is likely present. We report on a study of the formation of pyrimidine-based molecules, including nucleobases and other species of prebiotic interest, from the ultraviolet (UV) irradiation of pyrimidine in low temperature ices containing H2O, NH3, C3OH, and CH4, to simulate the astrophysical conditions under which prebiotic species may be formed in the Solar System.

Sandford, Scott↗

Nucleobases and other Prebiotic Species from the Ultraviolet Irradiation of Pyrimidine in Astrophysical Ices

Nucleobases are N-heterocycles that are the informational subunits of DNA and RNA, and are divided into two families: pyrimidine bases (uracil, cytosine, and thymine) and purine bases (adenine and guanine). Nucleobases have been detected in meteorites and their extraterrestrial origin confirmed by isotope measurement. Although no Nheterocycles have ever been observed in the ISM, the positions of the 6.2-m interstellar emission features suggest a population of such molecules is likely to be present. In this work we study the formation of pyrimidine-based molecules, including nucleobases, as well as other species of prebiotic interest, from the ultraviolet (UV) irradiation of pyrimidine in combinations of H2O, NH3, CH3OH, and CH4 ices at low temperature, in order to simulate the astrophysical conditions under which prebiotic species may be formed in the interstellar medium and icy bodies of the Solar System. Experimental: Gas mixtures are prepared in a glass mixing line (background pressure approx. 10(exp -6)-10(exp -5) mbar). Relative proportions between mixture components are determined by their partial pressures. Gas mixtures are then deposited on an aluminum foil attached to a cold finger (15-20 K) and simultaneously irradiated with an H2 lamp emitting UV photons (Lyman and a continuum at approx.160 nm). After irradiation samples are warmed to room temperature, at which time the remaining residues are recovered to be analyzed with liquid and gas chromatographies. Results: These experiments showed that the UV irradiation of pyrimidine mixed in these ices at low temperature leads to the formation of several photoproducts derived from pyrimidine, including the nucleobases uracil and cytosine, as well as their precursors 4(3H)-pyrimidone and 4-aminopyrimidine (Fig. 1). Theoretical quantum calculations on the formation of 4(3H)-pyrimidone and uracil from the irradiation of pyrimidine in pure H2O ices are in agreement with their experimental formation pathways. In those residues, other species of prebiotic interest such as urea and the amino acids glycine and alanine could also be identified. However, no pyrimidine derivatives containing CH3 groups, including the third nucleobase thymine, could be identified, suggesting that the addition of methyl groups to pyrimidine is not an efficient process.

DNA and RNA↗

Conceptual Design of an Experiment to Study Dust Destruction by Astrophysical Shock Waves

A novel laboratory experimental design is described that will investigate the processing of dust grains in astrophysical shocks. Dust is a ubiquitous ingredient in the interstellar medium (ISM) of galaxies; however, its evolutionary cycle is still poorly understood. Especially shrouded in mystery is the efficiency of grain destruction by astrophysical shocks generated by expanding supernova remnants. While the evolution of these remnants is fairly well understood, the grain destruction efficiency in these shocks is largely unknown. The experiments described herein will fill this knowledge gap by studying the dust destruction efficiencies for shock velocities in the range of approximately 10-30 kilometers per second (microns per nanosecond), at which most of the grain destruction and processing in the ISM takes place. The experiments focus on the study of grain-grain collisions by accelerating small (1 millimeter) dust particles into a large (approximately 5-10 millimeter diameter) population; this simulates the astrophysical system well in that the more numerous, small grains impact and collide with the large population. Facilities that combine the versatility of high-power optical lasers with the diagnostic capabilities of X-ray free-electron lasers, e.g., the Matter in Extreme Conditions instrument at the SLAC (originally named Stanford Linear Accelerator Center) National Accelerator Laboratory, provide an ideal laboratory environment to create and diagnose dust destruction by astrophysically relevant shocks at the micron scale.

Manuel, M. J.-E.↗

Multiple scattering theory for dense plasmas

Dense plasmas occur in stars, giant planets, and in inertial fusion experiments. Accurate modeling of the electronic structure of these plasmas allows for prediction of material properties that can in turn be used to simulate these astrophysical objects and terrestrial experiments. But modeling them remains a challenge. Here we explore the Korringa-Kohn-Rostoker Green's function (KKR-GF) method for this purpose. We find that it is able to predict equation of state in good agreement with other state-of-the-art methods, where they are accurate and viable. In addition, it is shown that the computational cost does not significantly change with temperature, in contrast with other approaches. Moreover, the method does not use pseudopotentials—core states are calculated self consistently. We conclude that KKR-GF is a very promising method for dense plasma simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fluorescence Spectroscopy of Gas-phase Polycyclic Aromatic Hydrocarbons

The purpose of this investigation was to produce fluorescence spectra of polycyclic aromatic hydrocarbon (PAH) molecules in the gas-phase for comparison with blue luminescence (BL) emission observed in astrophysical sources Vijh et al. (2004, 2005a,b). The BL occurs roughly from 350 to 450 nm, with a sharp peak near 380 nm. PAHs with three to four rings, e.g. anthracene and pyrene, were found to produce luminescence in the appropriate spectral region, based on existing studies. Relatively few studies of the gas-phase fluorescence of PAHs exist; those that do exist have dealt primarily with the same samples commonly available for purchase such as pyrene and anthracene. In an attempt to understand the chemistry of the nebular environment we also obtained several nitrogen substituted PAHs from our colleagues at NASA Ames. In order to simulate the astrophysical environment we also took spectra by heating the PAHs in a flame. The flame environment counteracts the formation of eximers and permits the spectroscopy of free-flying neutral molecules. Experiments with coal tar demonstrate that fluorescence spectroscopy reveals primarily the presence of the smallest molecules, which are most abundant and which possess the highest fluorescence efficiencies. One gas-phase PAH that seems to fit the BL spectrum most closely is phenanthridine. In view of the results from the spectroscopy of coal tar, a compound containing a mixture of PAHs ranging from small to very large PAH molecules, we can not preclude the presence of larger PAHs in interstellar sources exhibiting BL.

Thomas, J. D.↗

Formation of Nucleobases from the UV Irradiation of Pyrimidine in Astrophysical Ice Analogs

Nucleobases are the informational subunits of DNA and RNA. They consist of Nheterocycles that belong to either the pyrimidine-base group (uracil, cytosine, and thymine) or the purinebase group (adenine and guanine). Several nucleobases, mostly purine bases, have been detected in meteorites [1-3], with isotopic signatures consistent with an extraterrestrial origin [4]. Uracil is the only pyrimidine-base compound formally reported in meteorites [2], though the presence of cytosine cannot be ruled out [5,6]. However, the actual process by which the uracil was made and the reasons for the non-detection of thymine in meteorites have yet to be fully explained. Although no N-heterocycles have ever been observed in the ISM [7,8], the positions of the 6.2-μm interstellar emission features suggest a population of such molecules is likely to be present [9]. In this work we study the formation of pyrimidine-based molecules, including the three nucleobases uracil, cytosine, and thymine from the ultraviolet (UV) irradiation of pyrimidine in ices consisting of several combinations of H(sub2)O, NH(sub3), CH(sub3)OH, and CH(sub4) at low temperature, in order to simulate the astrophysical conditions under which prebiotic species may be formed in the interstellar medium, in the protosolar nebula, and on icy bodies of the Solar System.

UV Irradiation↗

Cross sections for neutron-induced reactions from surrogate data: Reexamining the Weisskopf-Ewing approximation for ( n , n ' ) and ( n , 2 n ) reactions

Background: Modeling nuclear reaction networks for nuclear science applications and for simulations of astrophysical environments relies on cross section data for a vast number of reactions, many of which have never been measured. Cross sections for neutron-induced reactions on unstable nuclei are particularly scarce, since they are the most difficult to measure. Consequently, we must rely on theoretical predictions or indirect measurements to obtain the requisite reaction data. For compound nuclear reactions, the surrogate reaction method can be used to determine many cross sections of interest. Purpose: Earlier work has demonstrated that cross sections for neutron-induced fission and radiative neutron capture can be determined from a combination of surrogate reaction data and theory. For the fission case, it was shown that the Weisskopf-Ewing approximation, which significantly simplifies the implementation of the surrogate method, can be employed. Capture cross sections cannot be obtained, and require a detailed description of the surrogate reaction process. Here, we examine the validity of the Weisskopf-Ewing approximation for determining unknown (n, n') and (n, 2n) cross sections from surrogate data. Methods: Using statistical reaction calculations with realistic parametrizations, we investigate first whether the assumptions underlying the Weisskopf-Ewing approximation are valid for (n, n') and (n, 2n) reactions on representative target nuclei. We then produce simulated surrogate reaction data and assess the impact of applying the Weisskopf-Ewing approximation when extracting (n, n') and (n, 2n) cross sections in situations where the approximation is not strictly justified. Results: We find that peak cross sections can be estimated using the Weisskopf-Ewing approximation, but the shape of the (n, n') and (n, 2n) cross sections, especially for low neutron energies, cannot be reliably determined without accounting for the angular-momentum differences between the neutron-induced and surrogate reaction. Conclusions: To obtain reliable (n, n') and (n, 2n) cross sections from surrogate reaction data, a detailed description of the surrogate reaction mechanisms is required. To do so for the compound-nucleus energies and decay channels relevant to these reactions, it becomes necessary to extend current modeling capabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

On the statistical theory of self-gravitating collisionless dark matter flow

Dark matter, if it exists, accounts for five times as much as the ordinary baryonic matter. Compared to hydrodynamic turbulence, the flow of dark matter might possess the widest presence in our universe. This paper presents a statistical theory for the flow of dark matter that is compared with N-body simulations. By contrast to hydrodynamics of normal fluids, dark matter flow is self-gravitating, long-range, and collisionless with a scale-dependent flow behavior. The peculiar velocity field is of constant divergence nature on small scale and irrotational on large scale. The statistical measures, i.e., correlation, structure, dispersion, and spectrum functions, are modeled on both small and large scales, respectively. Kinematic relations between statistical measures are fully developed for incompressible, constant divergence, and irrotational flow. Incompressible and constant divergence flow share the same kinematic relations for even order correlations. The limiting correlation of velocity $\mathrm{ρ_{L}=1/2}$ on the smallest scale ( r = 0) is a unique feature of collisionless flow (⁠$\mathrm{ρ_{L}=1}$ for incompressible flow). On large scale, transverse velocity correlation has an exponential form $T_{2}∝e^{–r/r_2}$ with a constant comoving scale r 2 =21.3 Mpc/h that may be related to the horizon size at matter–radiation equality. All other correlation, structure, dispersion, and spectrum functions for velocity, density, and potential fields are derived analytically from kinematic relations for irrotational flow. On small scale, longitudinal structure function follows one-fourth law of ${S}_{2}^{1}∝r^{1/4}$. All other statistical measures can be obtained from kinematic relations for constant divergence flow. Vorticity is negatively correlated for scale r between 1 and 7 Mpc/h. Divergence is negatively correlated for r > 30 Mpc/h that leads to a negative density correlation.

79 ASTRONOMY AND ASTROPHYSICS↗

Association between optically identified galaxy clusters and the underlying dark matter halos

Clusters of galaxies trace massive dark matter halos in the Universe, but they can include multiple halos projected along lines of sight. Here, we study the halos contributing to clusters using the Cardinal simulation, which mimics the Dark Energy Survey data. We use the red-sequence-based cluster finding algorithm redMaPPer as a case study. For each cluster, we identify the halos hosting its member galaxies, and we define the main halo as the one contributing the most to the cluster's richness ($λ$, the estimated number of member galaxies). At $z=0.3$, for clusters with $λ> 60$, the main halo typically contributes to $92\%$ of the richness, and this fraction drops to $67\%$ for $λ\approx 20$. Defining "clean" clusters as those with $\geq50\%$ of the richness contributed by the main halo, we find that $100\%$ of the $λ> 60$ clusters are clean, while $73\%$ of the $λ\approx 20$ clusters are clean. Three halos can usually account for more than $80\%$ of the richness of a cluster. The main halos associated with redMaPPer clusters have a completeness ranging from $98\%$ at virial mass $10^{14.6}~h^{-1}M_{\odot}$ to $64\%$ at $10^{14}~h^{-1}M_{\odot}$. In addition, we compare the inferred cluster centers with true halo centers, finding that $30\%$ of the clusters are miscentered with a mean offset $40\%$ of the cluster radii, in agreement with recent X-ray studies. These systematics worsen as redshift increases, but we expect that upcoming surveys extending to longer wavelengths will improve the cluster finding at high redshifts. Our results affirm the robustness of the redMaPPer algorithm and provide a framework for benchmarking other cluster-finding strategies.

79 ASTRONOMY AND ASTROPHYSICS↗

Dark energy survey: Modeling strategy for multiprobe cluster cosmology and validation for the full six-year dataset

Here, we introduce an updated To&Krause2021 model for joint analyses of cluster abundances and large-scale two-point correlations of weak lensing and galaxy and cluster clustering (termed CL+3×2 pt analysis) and validate that this model meets the systematic accuracy requirements of analyses with the statistical precision of the final Dark Energy Survey (DES) Year 6 (Y6) dataset. The validation program consists of two distinct approaches, (i) identification of modeling and parametrization choices and impact studies using simulated analyses with each possible model misspecification and (ii) end-to-end validation using mock catalogs from customized Cardinal simulations that incorporate realistic galaxy populations and DES-Y6-specific galaxy and cluster selection and photometric redshift modeling, which are the key observational systematics. In combination, these validation tests indicate that the model presented here meets the accuracy requirements of DES-Y6 for CL+3×2 pt based on a large list of tests for known systematics. In addition, we also validate that the model is sufficient for several other data combinations: the CL+GC subset of this data vector (excluding galaxy–galaxy lensing and cosmic shear two-point statistics) and the CL+3×2 pt+BAO+SN (combination of CL+3×2 pt with the previously published Y6 DES baryonic acoustic oscillation and Y5 supernovae data).

79 ASTRONOMY AND ASTROPHYSICS↗

DG-IMEX method for a two-moment model for radiation transport in the $\mathscr{O}$($v$/$c$) limit

Here, we consider neutral particle systems described by moments of a phase-space density and propose a realizability-preserving numerical method to evolve a spectral two-moment model for particles interacting with a background fluid moving with nonrelativistic velocities. The system of nonlinear moment equations, with special relativistic corrections to $\mathscr{O}$($v$/$c$), expresses a balance between phase-space advection and collisions and includes velocity-dependent terms that account for spatial advection, Doppler shift, and angular aberration. The model is conservative for the correct $\mathscr{O}$($v$/$c$) Eulerian-frame number density and is consistent, to $\mathscr{O}$($v$/$c$), with Eulerian-frame energy and momentum conservation. This model is closely related to the one promoted by Lowrie et al. and similar to models currently used to study transport phenomena in large-scale simulations of astrophysical environments. The proposed numerical method is designed to preserve moment realizability, which guarantees that the moments correspond to a nonnegative phase-space density. The realizability-preserving scheme consists of the following key components: (i) a strong stability-preserving implicit-explicit (IMEX) time-integration method; (ii) a discontinuous Galerkin (DG) phase-space discretization with carefully constructed numerical uxes; (iii) a realizability-preserving implicit collision update; and(iv) a realizability-enforcing limiter. In time integration, nonlinearity of the moment model necessitates solution of nonlinear equations, which we formulate as fixed-point problems and solve with tailored iterative solvers that preserve moment realizability with guaranteed global convergence. We also analyze the simultaneous Eulerian-frame number and energy conservation properties of the semi-discrete DG scheme and propose a "spectral redistribution" scheme that promotes Eulerian-frame energy conservation. Through numerical experiments, we demonstrate the accuracy and robustness of this DG-IMEX method and investigate its Eulerian-frame energy conservation properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The upgraded summing NaI(Tl) (SuN++) absorption spectrometer

Simulations of astrophysical processes require a plethora of nuclear physics input. In particular, models of neutron-capture nucleosynthesis like the s, i, and r processes require β-decay information and experimentally constrained neutron-capture reaction rates. Past experiments with the 4π Summing NaI(Tl) (SuN) total absorption spectrometer have provided these physics quantities. Here, we outline an upgrade of SuN to SuN++, where 20 new segments (12 NaI(Tl) and 8 CeBr 3 ) have been integrated into the pre-existing SuN total absorption spectrometer to provide increased energy and time resolution in β-decay experiments. The details of the newly upgraded SuN++ total absorption spectrometer are discussed with results from the commissioning experiment at the Facility for Rare Isotope Beams (FRIB) utilizing a 70 Cu beam.

CeBr3↗

Quantified limits of the nuclear landscape

The chart of the nuclides is limited by particle drip lines beyond which nuclear stability to proton or neutron emission is lost. Predicting the range of particle-bound isotopes poses an appreciable challenge for nuclear theory as it involves extreme extrapolations of nuclear masses well beyond the regions where experimental information is available. Still, quantified extrapolations are crucial for a wide variety of applications, including the modeling of stellar nucleosynthesis. We use microscopic nuclear global mass models, current mass data, and Bayesian methodology to provide quantified predictions of proton and neutron separation energies as well as Bayesian probabilities of existence throughout the nuclear landscape all the way to the particle drip lines. Here, we apply nuclear density functional theory with several energy density functionals. We also consider two global mass models often used in astrophysical nucleosynthesis simulations. To account for uncertainties, Bayesian Gaussian processes are trained on the separation-energy residuals for each individual model, and the resulting predictions are combined via Bayesian model averaging. This framework allows to account for systematic and statistical uncertainties and propagate them to extrapolative predictions. We establish and characterize the drip-line regions where the probability that the nucleus is particle- bound decreases from 1 to 0. In these regions, we provide quantified predictions for one- and two-nucleon separation energies. According to our Bayesian model averaging analysis, 7759 nuclei with Z ≤ 119 have a probability of existence ≥ 0.5. The extrapolation results obtained in this study will be put through stringent tests when new experimental information on existence and masses of exotic nuclei becomes available. In this respect, the quantified landscape of nuclear existence obtained in this study should be viewed as a dynamical prediction that will be fine-tuned when new experimental information and improved global mass models become available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimal multifrequency weighting for CMB lensing

We report that extragalactic foregrounds in cosmic microwave background (CMB) temperature maps lead to significant biases in CMB lensing reconstruction if not properly accounted for. Combinations of multifrequency data have been used to minimize the overall map variance (internal linear combination, or ILC), or specifically null a given foreground, but these are not tailored to CMB lensing. In this paper, we derive an optimal multifrequency combination to jointly minimize CMB lensing noise and bias. We focus on the standard lensing quadratic estimator, as well as the "shear-only"and source-hardened estimators, whose responses to foregrounds differ. We show that an optimal multifrequency combination is a compromise between the ILC and joint deprojection, which nulls the thermal Sunyaev-Zel'dovich (tSZ) and cosmic infrared background (CIB) contributions. In particular, for a Simons Observatory-like experiment with ℓ max,T =3000, we find that profile hardening alone (with the standard ILC) reduces the bias to the lensing power amplitude by 40%, at a 20% cost in noise, while the bias to the cross-correlation with a LSST-like sample is reduced by nearly an order of magnitude at a 10% noise cost, relative to the standard quadratic estimator. With a small amount of joint deprojection the bias to the profile hardened estimator can be further reduced to less than half the statistical uncertainty on the respective amplitudes, at a 20% and 5% noise cost for the auto- and cross-correlation respectively, relative to the profile hardened estimator with the standard ILC weights. Finally, we explore possible improvements with more aggressive masking and varying ℓ max,T .

79 ASTRONOMY AND ASTROPHYSICS↗

Dark branches of immortal stars at the Galactic Center

We show that stars in the inner parsec of the Milky Way can be significantly affected by dark matter annihilation, producing population-level effects that are visible in a Hertzsprung-Russell (HR) diagram. We establish the dark HR diagram, where stars lie on a new stable dark main sequence with similar luminosities, but lower temperatures, than the standard main sequence. The dark matter density in these stars continuously replenishes, granting these stars immortality and solving multiple stellar anomalies. Upcoming telescopes could detect the dark main sequence, offering a new dark matter discovery avenue.

astrophysical simulations↗