Search NASA⌕ Search

SEARCH · Search NASA

Results for “astronomy computing”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Problem-tailored Simulation of Energy Transport on Noisy Quantum Computers

The transport of conserved quantities like spin and charge is fundamental to characterizing the behavior of quantum many-body systems. Numerically simulating such dynamics is generically challenging, which motivates the consideration of quantum computing strategies. However, the relatively high gate errors and limited coherence times of today's quantum computers pose their own challenge, highlighting the need to be frugal with quantum resources. In this work we report simulations on quantum hardware of infinite-temperature energy transport in the mixed-field Ising chain, a paradigmatic many-body system that can exhibit a range of transport behaviors at intermediate times. We consider a chain with L = 12 sites and find results broadly consistent with those from ideal circuit simulators over 90 Trotter steps, containing up to 990 entangling gates. To obtain these results, we use two key problem-tailored insights. First, we identify a convenient basis – the Pauli Y basis – in which to sample the infinite-temperature trace and provide theoretical and numerical justifications for its efficiency relative to, e.g., the computational basis. Second, in addition to a variety of problem-agnostic error mitigation strategies, we employ a renormalization strategy that compensates for global nonconservation of energy due to device noise. We discuss the applicability of the proposed sampling approach beyond the mixed-field Ising chain and formulate a variational method to search for a sampling basis with small sample-to-sample fluctuations for an arbitrary Hamiltonian. This opens the door to applying these techniques in more general models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Exact evaluation of large-charge correlation functions in nonrelativistic conformal field theory

The large-charge master field which generates all n -point correlation functions with an insertion of large charge Q in nonrelativistic conformal field theory is obtained. This field is used to compute Schrödinger-invariant n -point correlation functions of large-charge operators via a direct evaluation of the path integral. Conformal dimensions are found to agree with calculations based on the state-operator correspondence. The master field solution exhibits an emergent harmonic trap whose frequency is a function of the Euclidean time. The large-charge effective action with operator insertions describes a droplet of superfluid matter whose spatial size scales with the time separation of sources. The solution is used to compute Schrödinger symmetry breaking corrections in the large-charge effective field theory (EFT) due to a finite scattering length in the fundamental theory of fermions near unitarity. The scaling of these effects in the large-charge power counting scheme is established, and the size of the effects is quantified using input from quantum Monte Carlo simulations of the near-unitary gas, as well as from the large- N expansion at large charge. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2026

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Heavy dark matter in white dwarfs: multiple-scattering capture and thermalization

We present an improved treatment for the scattering of heavy dark matter from the ion constituents of a white dwarf. In the heavy dark matter regime, multiple collisions are required for the dark matter to become gravitationally captured. Our treatment incorporates all relevant physical effects including the dark matter trajectories, nuclear form factors, and radial profiles for the white dwarf escape velocity and target number densities. Our capture rates differ by orders of magnitude from previous estimates, which have typically used approximations developed for dark matter scattering in the Earth. We also compute the time for the dark matter to thermalize in the center of the white dwarf, including in-medium effects such as phonon emission and absorption from the ionic lattice in the case where the star has a crystallized core. We find much shorter thermalization timescales than previously estimated, especially if the white dwarf core has crystallized. We illustrate the importance of our improved approach by determining the cross section required for accumulated asymmetric dark matter to self-gravitate.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Toward a Ginsparg-Wilson lattice Hamiltonian

To address quantum computation of quantities in quantum chromodynamics (QCD) for which chiral symmetry is important, it would be useful to have the Hamiltonian for a fermion satisfying the Ginsparg-Wilson (GW) equation. I work with a solution to the GW equation which is fractional linear in time derivatives. The resulting Hamiltonian is nonlocal and has ghosts, but is free of doublers and has the correct continuum limit. This construction works in general odd spatial dimensions, and I provide an explicit expression for the Hamiltonian in one spatial dimension. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Thermodynamics of the near-extremal Kerr spacetime

We examine the thermodynamics of a near-extremal Kerr black hole, and demonstrate that the geometry behaves as an ordinary quantum system with a vanishingly small degeneracy at low temperatures. This is in contrast with the classical analysis, which instead predicts a macroscopic entropy for the extremal Kerr black hole. Our results follow from a careful analysis of the gravitational path integral. Specifically, the low temperature canonical partition function behaves as Z ~ T$\frac{3}{2}$ e S 0 +c log S 0 , with S 0 the classical degeneracy and c a numerical coefficient we compute. This is in line with the general expectations for non-supersymmetric near-extremal black hole thermodynamics, as has been clarified in the recent past, although cases without spherical symmetry have not yet been fully analyzed until now. We also point out some curious features relating to the rotational zero modes of the near-extremal Kerr black hole background that affects the coefficient c. This raises a puzzle when considering similar black holes in string theory. Our results generalize to other rotating black holes, as we briefly exemplify.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Cosmic Ray Mediated Thermal Fronts in the Warm-Hot Circumgalactic Medium

We investigate the 1D plane-parallel front connecting the warm ($10^4$ K) and hot ($10^6$ K) phases of the circumgalactic medium (CGM), focusing on the influence of cosmic rays (CRs) in shaping these transition layers. We find that cosmic rays dictate the thermal balance while other fluxes (thermal conduction, radiative cooling, and gas flow) adjust to compensate. We compute column density ratios for selected transition temperature ions and compare them with observational data. While most of our models fail to reproduce the observations, a few are successful, although we make no claims for their uniqueness. Some of the discrepancies may indicate challenges in capturing the profiles in cooler, photoionized regions, as has been suggested for by previous efforts to model thermal transition layers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A-B Transition in Superfluid $^3$He and Cosmological Phase Transitions

First-order phase transitions in the very early universe are a prediction of many extensions of the Standard Model of particle physics and could provide the departure from equilibrium needed for a dynamical explanation of the baryon asymmetry of the Universe. They could also produce gravitational waves of a frequency observable by future space-based detectors such as the Laser Interferometer Space Antenna. All calculations of the gravitational wave power spectrum rely on a relativistic version of the classical nucleation theory of Cahn-Hilliard and Langer, due to Coleman and Linde. The high purity and precise control of pressure and temperature achievable in the laboratory made the first-order A to B transition of superfluid $^3$He ideal for test of classical nucleation theory. As Leggett and others have noted, the theory fails dramatically. The lifetime of the metastable A phase is measurable, typically of order minutes to hours, far faster than classical nucleation theory predicts. If the nucleation of B phase from the supercooled A phase is due to a new, rapid intrinsic mechanism that would have implications for first-order cosmological phase transitions as well as predictions for gravitational wave production in the early universe. Here we discuss studies of the A-B phase transition dynamics in $^3$He, both experimental and theoretical, and show how the computational technology for cosmological phase transition can be used to simulate the dynamics of the A-B transition, support the experimental investigations of the A-B transition in the QUEST-DMC collaboration with the goal of identifying and quantifying the mechanism(s) responsible for nucleation of stable phases in ultra-pure metastable quantum phases.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Electron and muon dynamics in neutron stars beyond chemical equilibrium

A neutron star harbors $\mathcal{O}$(10 56 ) electrons in its core, and almost the same number of muons, with muon decay prohibited by Pauli blocking. However, as macroscopic properties of the star such as its mass, rotational velocity, or magnetic field evolve over time, the equilibrium lepton abundances (dictated by the weak interactions) change as well. Scenarios where this can happen include spin-down, accretion, magnetic field decay, and tidal deformation. Here, we discuss the mechanisms by which a star disrupted in one of these ways re-establishes lepton chemical equilibrium. In most cases, the dominant processes are out-of-equilibrium Urca reactions, the rates of which we compute for the first time. If, however, the equilibrium muon abundance decreases, while the equilibrium electron abundance increases (or decreases less than the equilibrium muon abundance), outward diffusion of muons plays a crucial role as well. This is true in particular for stars older than about 10 4 yrs whose core has cooled to ≲ 20 keV. The muons decay in a region where Pauli blocking is lifted, and we argue that these decays lead to a flux of Script $\mathcal{O}$(10 MeV) neutrinos. Realistically, however, this flux will remain undetectable for the foreseeable future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Finite Bubble Statistics Constrain Late Cosmological Phase Transitions

We consider first order cosmological phase transitions (PTs) happening at late times below standard model temperatures T PT ≲ GeV. The inherently stochastic nature of bubble nucleation and the finite number of bubbles associated with a late-time PT lead to superhorizon fluctuations in the PT completion time. We compute how such fluctuations eventually source curvature fluctuations with universal properties, independent of the microphysics of the PT dynamics. Using cosmic microwave background (CMB) and large scale structure measurements, we constrain the energy released in a dark-sector PT. For 0.1 eV ≲ T PT ≲ keV this constraint is stronger than both the current bound from additional neutrino species Δ⁢N eff , and in some cases, even CMB-S4 projections. Future measurements of CMB spectral distortions and pulsar timing arrays will also provide competitive sensitivity for keV ≲ T PT ≲ GeV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Discovering Strong Gravitational Lenses in the Dark Energy Survey with Interactive Machine Learning and Crowd-sourced Inspection with Space Warps

We conduct a search for strong gravitational lenses in the Dark Energy Survey (DES) Year 6 imaging data. We implement a pre-trained Vision Transformer (ViT) for our machine learning (ML) architecture and adopt interactive machine learning to construct a training sample with multiple classes to address common types of false positives. Our ML model reduces ∼236 million DES cutout images to 22,564 targets of interest, including ∼85% of previously reported galaxy–galaxy lens candidates discovered in DES. These targets were visually inspected by citizen scientists, who ruled out ∼90% as false positives. Of the remaining 2618 candidates, 149 were expert-classified as “definite” lenses and 516 as “probable” lenses, for a total of 665 systems, with 147 of these candidates being newly identified. Additionally, we trained a second ViT to find double-source plane lens systems, finding at least one double-source system. Our main ViT excels at identifying galaxy–galaxy lenses, consistently assigning high scores to candidates with high expert assessments. The top 800 ViT-scored images include ∼100 of our “definite” lens candidates. This selection is an order of magnitude higher in purity than previous convolutional neural-network-based lens searches and demonstrates the feasibility of applying our methodology for discovering large samples of lenses in future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Kilonova light-curve interpolation with neural networks

Kilonovae are the electromagnetic transients created by the radioactive decay of freshly synthesized elements in the environment surrounding a neutron star merger. To study the fundamental physics in these complex environments, kilonova modeling requires, in part, the use of radiative transfer simulations. The microphysics involved in these simulations results in high computational cost, prompting the use of emulators for parameter inference applications. Utilizing a training set of 22 248 high-fidelity simulations (composed of 412 unique ejecta parameter combinations evaluated at 54 viewing angles), we use a neural network to efficiently train on existing radiative transfer simulations and predict light curves for new parameters in a fast and computationally efficient manner. Our neural network can generate millions of new light curves in under a minute. We discuss our emulator's degree of off-sample reliability and parameter inference of the AT2017gfo observational data. Finally, we discuss tension introduced by multiband inference in the parameter inference results, particularly with regard to the neural network's recovery of viewing angle. Published by the American Physical Society 2024

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Energy loss of a heavy fermion in a collisional QED plasma

We compute the energy loss of heavy fermions moving in a plasma, taking into account the modification of the photon collective modes induced by collisions using a Bhatnagar-Gross–Krook collisional kernel. We include contributions from both hard and soft scatterings of the heavy fermion using a collisionally modified hard-thermal-loop resummed propagator. Using this method, one does not need to introduce a separation scale between hard- and soft-momentum exchanges. To place our calculation in context, we review other theoretical approaches to computing the collisional energy loss of fermions and discuss the systematics and results obtained in each approach compared to using a resummed propagator for both hard and soft momentum exchanges. Our final results indicate that self-consistently including the effect of collisions in the self-energies of the resummed propagator results in an increased energy loss compared to using collisionless hard-thermal-loop propagators. The effect becomes larger as the magnitude of the coupling constant and the velocity of the fermion increase. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Boundary of the Distribution of Solar Wind Proton Beta versus Temperature Anisotropy

Abstract The frequency distribution of solar wind protons, measured in the vicinity of Earth’s orbit, is customarily plotted in ( β ∥ , T ⊥ / T ∥ ) phase space. Here, T ⊥ / T ∥ is the ratio of perpendicular and parallel temperatures, and β ∥ = 8 π nT ∥ / B 2 is the ratio of parallel thermal energy to background magnetic field energy, the so-called “parallel beta,” with ⊥ and ∥ denoting directions with respect to the ambient magnetic field. Such a frequency distribution, plotted as a two-dimensional histogram, forms a peculiar rhombic shape defined with an outer boundary in the said phase space. Past studies reveal that the threshold conditions for temperature anisotropy–driven plasma instability partially account for the boundary on the high- β ∥ side. The low- β ∥ side remains largely unexplained despite some efforts. Work by Vafin et al. recently showed that certain contours of collisional relaxation frequency, ν pp , when parameterized by T ⊥ / T ∥ and β ∥ , could match the overall shape of the left-hand boundary, thus suggesting that the collisional relaxation process might be closely related to the formation of the left-hand boundary. The present paper extends the analysis by Vafin et al. and carries out the dynamical computation of the collisional relaxation process for an ensemble of initial proton states with varying degrees of anisotropic temperatures. The final states of the relaxed protons are shown to closely match the observed boundary to the left of the ( β ∥ , T ⊥ / T ∥ ) phase space. When coupled with a similar set of calculations for the ensemble in the collective instability regime, it is found that the combined collisional/collective effects provide the baseline explanation for the observation.

Astronomy & Astrophysics↗
DOE OSTI2025

Fast, differentiable, and extensible big bang nucleosynthesis package

Here, we introduce light isotope nucleosynthesis with JAX (LINX), a new differentiable public big bang nucleosynthesis code designed for fast parameter estimation. By leveraging JAX, LINX achieves both speed and differentiability, enabling the use of Bayesian inference, including gradient-based methods. We discuss the formalism used in LINX for rapid primordial elemental abundance predictions and give examples of how LINX can be used. When combined with differentiable cosmic microwave background power spectrum emulators, LINX can be used for joint cosmic microwave background and big bang nucleosynthesis analyses without requiring extensive computational resources, including on personal hardware.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Bare mass effects on the reheating process after inflation

We consider the effects of a bare mass term for the inflaton, when the inflationary potential takes the form V(Φ)=λΦ k about its minimum with k ≥ 4. We concentrate on k =4, but discuss general cases as well. Further, we assume $λΦ$$^{2}_{end}$ >> $m$$^{2}_{Φ}$, where Φ end is the inflaton field value when the inflationary expansion ends. We show that the presence of a mass term (which may be present due to radiative corrections or supersymmetry breaking) can significantly alter the reheating process, as the equation of state of the inflaton condensate changes from w Φ = $\frac{1}{3}$ to w Φ = 0 when λ⁢Φ 2 drops below $m$$^{2}_{Φ}$. We show that, for a mass mΦ ≳ 3⁢λ $\frac{1}{4}$ T RH , the mass term will dominate at reheating. The value of λ is relatively model independent as it is normalized by the cosmic microwave background perturbation spectrum. For T models of inflation, this leads to m Φ ≳ T RH /250. We compute the effects on the reheating temperature for cases where reheating is due to inflaton decay (to fermions, scalars, or vectors) or to inflaton scattering (to scalars or vectors). For scattering to scalars and in the absence of a decay, there is always a residual inflaton background that acts as cold dark matter. In this case, we derive a strong upper limit to the inflaton bare mass which for T models is m Φ < 350 MeV⁢(T RH /10 10 GeV) 3/5 . We also consider the effect of the bare mass term on the fragmentation of the inflaton condensate.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Millimeter-Wave Superconducting Qubit

Manipulating the electromagnetic spectrum at the single-photon level is fundamental for quantum experiments. In the visible and infrared ranges, this can be accomplished with atomic quantum emitters, and with superconducting qubits such control is extended to the microwave range (below 10 GHz). Meanwhile, the region between these two energy ranges presents an unexplored opportunity for innovation. We bridge this gap by scaling up a superconducting qubit to the millimeter-wave range (near 100 GHz). Working in this energy range greatly reduces sensitivity to thermal noise compared to microwave devices, enabling operation at significantly higher temperatures, up to 1 K. This has many advantages by removing the dependence on rare 3⁢ He for refrigeration, simplifying cryogenic systems, and providing orders-of-magnitude higher cooling power, lending the flexibility needed for novel quantum sensing and hybrid experiments. Using low-loss niobium trilayer junctions, we realize a qubit at 72 GHz cooled to 0.87 K using only 4 ⁢He. We perform Rabi oscillations to establish control over the qubit state, and measure relaxation and dephasing times of 15.8 and 17.4 ns, respectively. This demonstration of a millimeter-wave quantum emitter offers exciting prospects for enhanced sensitivity thresholds in high-frequency photon detection, provides new options for quantum transduction and for scaling up and speeding up quantum computing, enables integration of quantum systems where 3 ⁢He refrigeration units are impractical, and, importantly, paves the way for quantum experiments exploring a novel energy range.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗