SEARCH · Search NASA
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.
Gravitational confinement is clearly very effective but what can we do on Earth? Magnetic confinement has been studied since around 1950. Inertial confinement has been associated with lasers for over 50 years. We have made steady progress in both MCF and ICF over the last half-century.
Cislunar space, encompassing the region from geosynchronous orbit to beyond the Moon, is poised to become a cornerstone for future exploration, scientific discovery, and national security. Missions in this region, spanning durations from weeks to decades, require robust infrastructure and reliable transit capabilities. The complex gravitational influences of the Moon, Sun, and planets, along with thermal radiation from Earth and the Sun, lead to significant trajectory deviations, resulting in kilometer-scale errors within days. Leveraging the high-performance computing resources at Lawrence Livermore National Laboratory (LLNL), we have simulated one million high-fidelity cislunar trajectories, now publicly available via LLNL’s Green Data Oasis and the Unified Data Library. Generated using the open-source Space Situational Awareness Python package, these trajectories match the precision of commercial tools such as AGI’s Systems Tool Kit and NASA’s General Mission Analysis Tool. This data set is a valuable resource for reference, statistical analysis of cislunar orbit populations, and training machine learning models for rapid orbit classification with minimal observational input. Preliminary analysis reveals stable bands in Keplerian element space, particularly around five geosynchronous radii across a range of inclinations and eccentricities. Beyond this threshold, the Moon’s influence disrupts most unassisted orbits, though co-orbiting L4/L5 Lunar Trojans persist throughout the six-year simulation.
SUMMARY Building on global adjoint tomography model GLAD-M15, we present transversely isotropic global model GLAD-M25, which is the result of 10 quasi-Newton tomographic iterations with an earthquake database consisting of 1480 events in the magnitude range 5.5 ≤ Mw ≤ 7.2, an almost sixfold increase over the first-generation model. We calculated fully 3-D synthetic seismograms with a shortest period of 17 s based on a GPU-accelerated spectral-element wave propagation solver which accommodates effects due to 3-D anelastic crust and mantle structure, topography and bathymetry, the ocean load, ellipticity, rotation and self-gravitation. We used an adjoint-state method to calculate Fréchet derivatives in 3-D anelastic Earth models facilitated by a parsimonious storage algorithm. The simulations were performed on the Cray XK7 ‘Titan’ and the IBM Power 9 ‘Summit’ at the Oak Ridge Leadership Computing Facility. We quantitatively evaluated GLAD-M25 by assessing misfit reductions and traveltime anomaly histograms in 12 measurement categories. We performed similar assessments for a held-out data set consisting of 360 earthquakes, with results comparable to the actual inversion. We highlight the new model for a variety of plumes and subduction zones.
Axisymmetric disks of high-eccentricity, low-mass bodies on near-Keplerian orbits are unstable to an out-of-plane buckling. This “inclination instability” exponentially grows the orbital inclinations, raises perihelion distances, and clusters in the argument of perihelion. Here we examine the instability in a massive primordial scattered disk including the orbit-averaged gravitational influence of the giant planets. We show that differential apsidal precession induced by the giant planets will suppress the inclination instability unless the primordial mass is ≳20 Earth masses. We also show that the instability should produce a “perihelion gap” at semimajor axes of hundreds of astronomical units, as the orbits of the remnant population are more likely to have extremely large perihelion distances (O(100 au)) than intermediate values.
Mutual gravitation between a pulsating star and an orbital companion leads to a time-dependent variation in path length for starlight traveling to Earth. These variations can be used for coherently pulsating stars, such as the δ Scuti variables, to constrain the masses and orbits of their companions. Observing these variations for δ Scuti stars has previously relied on subdividing the light curve and measuring the average pulsation phase in equally sized subdivisions, which leads to undersampling near periapsis. We introduce a new approach that simultaneously forward models each sample in the light curve and show that this method improves upon current sensitivity limits—especially in the case of highly eccentric and short-period binaries. We find that this approach is sensitive enough to observe Jupiter mass planets around δ Scuti stars under ideal conditions, and use gravity-mode pulsations in the subdwarf B star KIC 7668647 to detect its companion without radial velocity data. We further provide robust detection limits as a function of the signal-to-noise ratio of the pulsation mode and determine that the minimum detectable light travel time amplitude for a typical Kepler δ Scuti is around 2 s. This new method significantly enhances the application of light travel time variations to detecting short-period binaries with pulsating components, and pulsating A-type exoplanet host stars, especially as a tool for eliminating false positives.
We answer the question of whether long-period comet C/1961 T1 (Seki) is the parent body of a meteor shower observed in Earth’s atmosphere. We model discrete parts of a meteoroid stream associated with this comet and follow the orbital evolution of its particles via a numerical integration of their orbits. The gravitational perturbations of major planets and the Poynting–Robertson effect, which influences the motion of the particles, are taken into account. We predict two showers which are related to C/1961 T1. One of these showers was identified with the December ρ-Virginids, #502. This identification confirms an earlier suggested relationship between the shower and the comet. The second predicted shower was identified with the γ-Sagittariids, #657. The radiant areas of both showers are symmetrical in respect to the apex of Earth’s heliocentric motion.
The observation of electromagnetic radiation from radio to γ-ray wavelengths has provided a wealth of information about the Universe. However, at PeV (10 15 eV) energies and above, most of the Universe is impenetrable to photons. New messengers, namely cosmic neutrinos, are needed to explore the most extreme environments of the Universe where black holes, neutron stars, and stellar explosions transform gravitational energy into non-thermal cosmic rays. These energetic particles havemillions of times higher energies than those produced in the most powerful particle accelerators on Earth. As neutrinos can escape from regions otherwise opaque to radiation, they allow an unique view deep into exploding stars and the vicinity of the event horizons of black holes. The discovery of cosmic neutrinos with IceCube has opened this new window on the Universe. IceCube has been successful in finding first evidence for cosmic particle acceleration in the jet of an active galactic nucleus. Yet, ultimately, its sensitivity is too limited to detect even the brightest neutrino sources with high significance, or to detect populations of less luminous sources. In thiswhite paper, we present an overview of a next-generation instrument, IceCube-Gen2, which will sharpen our understanding of the processes and environments that govern the Universe at the highest energies. IceCube-Gen2 is designed to: (a) Resolve the high-energy neutrino sky from TeV to EeV energies (b) Investigate cosmic particle acceleration through multi-messenger observations (c) Reveal the sources and propagation of the highest energy particles in the Universe (d) Probe fundamental physics with high-energy neutrinos IceCube-Gen2 will enhance the existing IceCube detector at the South Pole. It will increase the annual rate of observed cosmic neutrinos by a factor of ten compared to IceCube, and will be able to detect sources five times fainter than its predecessor. Furthermore, through the addition of a radio array, IceCube- Gen2 will extend the energy range by several orders of magnitude compared to IceCube. Construction will take 8 years and cost about $350M. The goal is to have IceCube-Gen2 fully operational by 2033. IceCube-Gen2 will play an essential role in shaping the new era of multimessenger astronomy, fundamentally advancing our knowledge of the highenergy Universe. Overall, this challenging mission can be fully addressed only through the combination of the information from the neutrino, electromagnetic, and gravitational wave emission of high-energy sources, in concert with the new survey instruments across the electromagnetic spectrum and gravitational wave detectors which will be available in the coming years.
We propose the first method for water Cherenkov detectors to constrain GeV-scale dark matter (DM) below the solar evaporation mass. While previous efforts have highlighted the Sun and Earth as DM capture targets, we demonstrate that Jupiter is a viable target. Jupiter’s unique characteristics, such as its lower core temperature and significant gravitational potential, allow it to capture and retain light DM more effectively than the Sun, particularly in the mass range below 4 GeV where direct detection sensitivity diminishes. Our calculations provide the first sensitivity estimates to GeV-scale annihilating DM within Jupiter, predicting Hyper-K can reach spin dependent cross sections as low as $𝜎^{SD}_{𝑝𝜒}$ = 2×10 −35 cm 2 for DM masses below 2 GeV. This surpasses current solar limits and direct detection results. We additionally provide estimates for Super-K ORCA, and the IceCube-Upgrade, showing that these experiments could provide complimentary bounds to direct detection experiments.
The Sun represents a promising target for indirect dark matter searches, as dark matter particles from the Galactic halo can be gravitationally trapped in its core or in external orbits, and their annihilations can lead to final states with standard model particles that are able to reach the Earth. In this work we have considered a scenario in which dark matter particles can annihilate into pairs of long-lived mediators, which in turn can escape from the Sun and decay into pairs of gamma rays or into the bb̅, τ+τ-, μ+μ- channels, with the production of gamma rays in the final states. All these processes are expected to yield an excess in the energy spectrum of gamma rays towards the Sun. We have therefore analyzed the data collected by the Fermi Large Area Telescope during its first 13.5 years of operation, searching for possible excesses in the solar gamma-ray spectrum. Since no statistically significant excess is found, we have set constraints on the dark matter-nucleon scattering cross sections in both the spin-dependent and spin-independent cases. For all the mediator decay channels explored and for dark matter masses between a few GeV/c 2 and 1 TeV/c 2 , we have found that the upper limits on the spin-dependent and spin-independent cross sections are in the ranges from 10 -45 to 10 -39 cm 2 and from 10 -47 up to 10 -42 cm 2 , respectively.
The cosmic neutrino background ( C ν B ) constitutes the last observable prediction of the standard cosmological model, which has yet to be detected directly. In this work, we show how the coherent scattering of neutrinos off dense neutron matter can lead to an additional cooling channel in neutron stars (NSs). We also include the effects of gravitational capture and boosting, but find that the cooling is efficient only in the presence of large overdensities. We further discuss the prediction of a boosted C ν B flux on Earth from nearby NSs and the potential detection prospects in the case of a future nearby galactic supernova. Although currently these ideas do not offer any detection prospects, they can be used to constrain overdensities η ≲ 10 11 – 10 14 on short length scales O ( 10 km ) . We also discuss the impact of new physics scenarios, such as long-range forces, on NS cooling through the C ν B . Published by the American Physical Society 2025
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.
Supercooled first-order phase transition (FOPT) can lead to the formation of primordial black holes (PBHs). This scenario imposes stringent requirements on the profile of the effective potential. In this work, we use the singlet extended Standard Model (xSM) as a benchmark model to investigate this possibility at the electroweak scale. The PBHs formed during a supercooled FOPT have a narrow mass distribution around the mass of Earth. This distribution is closely tied to the temperature at which the PBHs form, corresponding to the FOPT at the electroweak scale. This scenario can be probed with microlensing experiments, space-based gravitational wave detectors, and collider experiments. Remarkably, the future space-based gravitational wave detector LISA will hold the potential to either confirm this scenario that leads to PBH formation across the observable Universe in the xSM or completely rule it out. Interestingly, our findings suggest that PBHs within the xSM framework may align with observations of the six ultrashort timescale events reported by the OGLE microlensing experiment. Published by the American Physical Society 2025
Neutron star mergers are an important source of information for nuclear physics. Except for black holes, neutron stars are the densest macroscopic objects known to exist in the Universe. They provide us with a remarkable laboratory to study the poorly understood high-density states of matter, complementing the study of heavy nuclei performed by nuclear physics experiments on Earth. Collisions of neutron stars are of particular interest. They are among the most energetic events observable in the Universe, powering a broad range of signals across the electromagnetic spectrum as well as gravitational wave signals. These signals contain important information about the properties of extremely dense matter. Colliding neutron stars additionally eject large amounts of neutron-rich material into the surrounding interstellar medium -- material that then undergoes rapid neutron-capture (r-process) nucleosynthesis, the mechanism thought to be responsible for the production of about half of the heavy nuclei. In order to understand current and future observations of neutron star mergers, we need reliable models for the signals that they power. An important component in the construction of such models are numerical simulations of colliding neutron stars. These simulations are costly, running for multiple months on supercomputers, and require the inclusion of complex physics (general relativity, magnetohydrodynamics, neutrino physics, nuclear reactions). In this document, we describe the development of new methods for the treatment of neutrinos in merger simulations, as well as the first simulations capable of evolving the equations of neutrino radiation transport directly and their impact on our understanding of neutron star mergers and more broadly in nuclear astrophysics. We also discuss additional work performed in the study of nucleosynthesis and neutrino physics in neutron star mergers as part of Early Career Award DE-SC0020435.
Axion dark matter can resonantly convert to photons in the magnetosphere of neutron stars, possibly giving rise to radio signals observable on Earth. This method for the indirect detection of axion dark matter has recently received significant attention in the literature. The calculation of the radio signal is complicated by a number of effects; most importantly, the gravitational infall of the axions onto the neutron star accelerates them to semi-relativistic speed, and the neutron star magnetosphere is highly anisotropic. Both of these factors complicate the calculation of the conversion of axions to photons. In this work, we present the first fully three-dimensional calculation of the axion-photon conversion in highly magnetised anisotropic media. Depending on the axion trajectory, this calculation leads to orders-of-magnitude differences in the conversion compared to the simplified one-dimensional calculation used so far in the literature, altering the directionality of the produced photons. Our results will have important implications for the radio signal one would observe in a telescope.
With ground-based interferometers detecting hundreds of gravitational-wave (GW) events, GW astronomy has continued to blossom. U.S. and European scientists are developing plans to construct third-generation ground-based interferometers. ESA has proceeded through the mission formulation phase of space-based interferometer in a lower-frequency band, 10 –4 –0.1 Hz. Despite all these exciting developments, there is still a missing frequency band, 0.1–10 Hz. This mid-frequency band is rich with interesting astrophysical events. Coalescence and merger of intermediate-mass black holes (IMBHs) will occur in this frequency band. Coalescing stellar-mass BHs will pass through this frequency band days before they reach the frequency band of LIGO and Virgo. Detection of such signals would enable a mid-frequency detector to issue an advance notice to the high-frequency GW detectors, as well as to optical, x-ray and γ-ray telescopes. We propose Lunar Accelerometer Network Gravitational Observatory (LANGO) to detect GWs in this frequency band. In the first phase (LANGO 1), we propose to deploy four ambient-temperature (250 K) accelerometers in the tetrahedral or in a square configuration on the hemisphere facing the Earth. After successful operation at 250 K, LANGO would be upgraded to a cryogenic (4 K) version with over two orders of magnitude increased sensitivity (LANGO 2), with coherent rejection of seismic noise implemented. LANGO is a full-tensor detector, capable of determining the source direction and wave polarization. Each test mass (TM) is suspended as a pendulum with resonance frequency ∼ 0.01 Hz, thus is only weakly coupled to the lunar surface horizontally. LANGO is designed to detect the relative motion of globally separated, nearly free, TMs by using the Moon as a large quiet platform. LANGO 1 and 2 accelerometers aim at sensitivities ⩽ 10 –11 m s –2 Hz –1/2 and ⩽ 10 –13 m s –2 Hz –1/2 in the horizontal axes over the frequency band of 1 mHz–10 Hz, which yield GW sensitivities 1.1 x 10 -21 Hz -1/2 and 3.8 x 10 -24 Hz -1/2 at 1 Hz, respectively. The LANGO accelerometers will be 10 3 –10 5 times more sensitive than Apollo seismometers. With such sensitivity, LANGO will also make great contribution to the advancement of lunar geophysics.
Viscous and gravitational flow instabilities cause a displacement front to break up into finger-like fluids. The detection and evolutionary analysis of these fingering instabilities are critical in multiple scientific disciplines such as fluid mechanics and hydrogeology. However, previous detection methods of the viscous and gravitational fingers are based on density thresholding, which provides limited geometric information of the fingers. The geometric structures of fingers and their evolution are important yet little studied in the literature. In this work, we explore the geometric detection and evolution of the fingers in detail to elucidate the dynamics of the instability. Additionally, we propose a ridge voxel detection method to guide the extraction of finger cores from three-dimensional (3D) scalar fields. After skeletonizing finger cores into skeletons, we design a spanning tree based approach to capture how fingers branch spatially from the finger skeletons. Finally, we devise a novel geometric-glyph augmented tracking graph to study how the fingers and their branches grow, merge, and split over time. Feedback from earth scientists demonstrates the usefulness of our approach to performing spatio-temporal geometric analyses of fingers.
Here, diverting hazardous small bodies on impact trajectories with the Earth can in some circumstances be impossible without risking disrupting them. Disruption is a much more difficult planetary defense scenario to assess, being linked both to the response of the body to shock loading and the much more complicated gravitational dynamics of the fragments in the solar system relative to pure deflection scenarios. In this work we present a new simulation suite built on N-body gravitational methods that solves fragment orbits in the full gravitational system without recourse to more approximate methods. We assess the accuracy of our simulations and the simplifying assumptions we adopt to make the system tractable, and then discuss in more detail several specific, plausible planetary defense scenarios based on real close encounters. We find that disruption can be a very effective planetary defense strategy even for very late (sub-year) interventions, and should be considered an effective backup strategy should preferred methods, which require long warning times, fail.
Phase separation, interfacial tension gradients, and molecular self-assembly are all phenomena frequently observed when manipulating soft condensed matter systems on Earth. Depending on the system these processes can also be quite predictable. For instance, on Earth water droplets placed on hydrophobic surfaces will have contact angles > 90° with limited spread, contrary to their wetting behavior on hydrophilic surfaces. However, in an environment that is in a constant state of free-fall, where the gravitational field is reduced by 10%, such as on the International Space Station (ISS), buoyancy and its effect on fluid convection also become negligible. Fluid flow and phase behavior become less predictable. Thus, continuing the prior example, even if water droplets still remain in contact with a surface in microgravity (µg), the interface between the air and droplets will still remain significantly curved, appearing almost as if they were placed over completely nonwetting surfaces (Fig. 1a). Moreover, the negligible opposing buoyant force would also mean that if immiscible liquid droplets of a less dense character, or air bubbles even (Fig. 1b), were to be suspended within a larger water drop, these immiscible droplets or air bubbles would not migrate toward any one direction unless either artificially propelled by an external force or internal reaction inducing flow.