Search NASA⌕ Search

SEARCH · Search NASA

Results for “continuum”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Modeling laser interactions with aluminum and tantalum targets using a hybrid atomistic-continuum model

A hybrid atomistic-continuum method can model the microstructure evolution of metals subjected to laser irradiation. This method combines classical molecular dynamics (MD) simulations with the two-temperature model (TTM) to account for the laser energy absorption and heat diffusion behavior. Accurate prediction of the temperature evolution in the combined MD-TTM method requires reliable accuracy in electron heat capacity, electron thermal conductivity, and electron–phonon coupling factor across the temperatures generated. This study uses the electronic density of states (DOS) obtained from first-principle calculations. The calculated electron temperature-dependent parameters are used in MD-TTM simulations to study the laser metal interactions in FCC and BCC metals and the phenomenon of laser shock loading and melting. This study uses FCC Al and BCC Ta as model systems to demonstrate this capability. When subjected to short pulsed laser shocks, the dynamic failure behavior predicted using temperature-dependent parameters is compared with the experimentally reported single-crystal and nanocrystalline Al and Ta systems. The MD-TTM simulations also investigate laser ablation and melting behavior of Ta to compare with the ablation threshold reported experimentally. In conclusion, this manuscript demonstrates that integrating the temperature-dependent parameters into MD-TTM simulations leads to the accurate modeling of the laser–metal interaction and allows the prediction of the kinetics of the solid–liquid interface.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Polarization consistent dielectric screening in polarizable continuum model calculations of solvation energies

A polarization consistent framework, where dielectric screening is affected consistently in polarizable continuum model (PCM) calculations, is employed for the study of solvation energies. The computational framework combines a screened range-separated-hybrid functional (SRSH) with PCM calculations, SRSH-PCM, where dielectric screening is imposed in both PCM self-consistent reaction field (SCRF) iterations and the electronic structure Hamiltonian. We begin by demonstrating the impact of modifying the Hamiltonian to include such dielectric screening in SCRF iterations by considering the solutions of electrostatically embedded Hartree–Fock (HF) exact exchange equations. Long-range screened HF-PCM calculations are shown to capture properly the linear dependence of gap energy of frontier orbitals on the inverse of the dielectric constant, whereas unscreened HF-PCM orbital energies are fallaciously semi-constant with respect to the dielectric constant and, therefore, inconsistent with the ionization energy gaps. Similar trends affect density functional theory (DFT) calculations that aim to achieve predictive quality. Importantly, the dielectric screened calculations are shown to significantly affect DFT- and HF PCM-based solvation energies, where screened solvation energies are smaller compared to the unscreened values. Importantly, SRSH-PCM, therefore, appears to reduce the tendency of DFT-PCM to overestimate solvation energies, where we find the effect to increase with the dielectric constant and the polarity of the molecular solute, trends that enhance the quality of DFT-PCM calculations of solvation energy. Understanding the relationship of dielectric screening in the Hamiltonian and DFT-PCM calculations can ultimately benefit on-going efforts for the design of predictive and parameter free descriptions of solvation energies.

Chemistry↗
DOE OSTI2023

X-ray continuum spectroscopy of inertial confinement fusion implosions at the National Ignition Facility

A methodology for measuring x-ray continuum spectra of inertial confinement fusion (ICF) implosions is described. The method relies on the use of ConSpec, a high-throughput spectrometer using a highly annealed pyrolytic graphite crystal , which measures the spectra in the ≃20–30 keV range. Due to its conical shape, the crystal is sagittally focusing a Bragg-reflected x-ray spectrum into a line, which enhances the recorded x-ray emission signal above the high neutron-induced background accompanying ICF implosions at the National Ignition Facility. To improve the overall measurement accuracy, the sensitivity of the spectrometer measured in an off-line x-ray laboratory setting was revised. The error analysis was expanded to include the accuracy of the off-line measurements, the effect of the neutron-induced background, as well as the influence of possible errors in alignment of the instrument to the ICF target. Here we demonstrate how the improved methodology is applied in the analysis of ConSpec data with examples of a relatively low-neutron-yield implosion using a tritium–hydrogen–deuterium mix as a fuel and a high-yield deuterium–tritium (DT) implosion producing high level of the background. In both cases, the shape of the measured spectrum agrees with the exponentially decaying spectral shape of bremsstrahlung emission to within ±10%. In the case of the high-yield DT experiment, non-monotonic deviations slightly exceeding the measurement uncertainties are observed and discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Calibration of a broadband x-ray crystal spectrometer using a continuum x-ray source and photon-counting detectors

We report measurements of sensitivity of a broadband (≃20–30 keV) x-ray crystal spectrometer using a bremsstrahlung continuum x-ray source and several detectors, including an energy-discriminating photon-counting point detector (Si drift diode, Amptek Inc.), imaging hybrid photon-counting detectors (Eiger2-Si and Eiger2-CdTe, DECTRIS Ltd.), and image plates (SR type). Sensitivity is defined as a ratio of the crystal-reflected energy-position-dispersed spectrum measured by a given detector to the incident non-dispersed spectrum measured using the energy-discriminating point detector. The sensitivity derived from data measured exclusively by the point detector is considered detector-independent and serves as a reference. The sensitivities derived from data collected with the hybrid photon-counting detectors were matched to the reference using a scale factor of 1.05. The image plate-derived sensitivities required a scale factor of 1.16 to match the reference. Furthermore, the resulting mismatch in the shapes of all measured scaled sensitivities in the range of the spectrometer was ≲±10%, while the mismatch between the shapes of the sensitivities corresponding to the Si-based detectors was ≲±2.5%.

Crystal optics↗
DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Inert mesoscale simulations

Mesoscale modeling of shock waves in Ni+Al multilayers poses significant challenges that are due, in part, to shock-induced chemical reactions. Current modeling approaches utilize reactive molecular dynamics (MD), but they are limited to resolving domains of only a few hundred nanometers. In contrast, actual multilayer superlattices can be tens of micrometers thick, and they exhibit non-ideal (i.e., wavy) interfaces. The second part of our research builds upon previous work developing physically based, thermodynamically complete equations of state for various Ni and Al intermetallic compositions. Here, we introduce a novel workflow for high-fidelity mesoscale simulations of Ni+Al multilayers using a continuum hydrocode. By increasing the simulation domain size beyond MD limitations (e.g., 2 × 6 μm 2 ) and incorporating explicit interfacial roughness, we investigate the shock response of Ni+Al multilayers at previously unexplored scales. Our experimental design encompasses nine multilayer geometries with varying roughness amplitudes and tilt angles (θ = 15°, 30°, and 45°), alongside 19 flyer impact velocities ranging from 0.3 to 3.0 km/s, resulting in a total of 171 high-fidelity simulations. The bulk shock state from inert 2D mesoscale simulations aligns with the law of mixtures, while temperature and pressure fluctuations strongly correlate with multilayer geometry types. A new metric dubbed the “hot spot probability integral” shows a greater dependence on a tilt angle than interfacial roughness.

Kittell, David E. [Sandia National Laboratories (S↗
DOE OSTI2025

Semi-implicit continuum kinetic modeling of weakly collisional parallel transport in a magnetic mirror

We present implicit-explicit (IMEX) kinetic simulations of weakly collisional parallel plasma transport in magnetic mirror configurations using the continuum code COGENT. The numerical scheme employs a Jacobian-free Newton–Krylov method with algebraic multigrid preconditioning to overcome the severe time step limitations imposed by strong mirror forces in fully explicit schemes. Applied to parameters relevant to the Wisconsin HTS Axisymmetric Mirror experiment, the IMEX approach enables time steps up to 2.5×10 4 times larger than those permitted by explicit methods, resulting in a 2500× speedup in 1D–2V simulations of parallel transport with kinetic ions and Boltzmann electrons. Additionally, a reduced bounce-averaged model for a square mirror is implemented to support the computationally intensive fully kinetic simulations. The bounce-averaged formulation is used to evaluate the numerical convergence of the velocity-space discretization algorithms and to assess the role of the collision model by comparing simulations employing the nonlinear Fokker–Planck and the simplified Lenard–Bernstein–Dougherty collision operators.

Collision theories↗
DOE OSTI2026

Investigating shock-induced chemical reactions in Ni+Al multilayers: A continuum-based mesoscale approach with Arrhenius kinetics and artificial thermal conduction

A new continuum-based mesoscale modeling approach for shock-induced chemical reactions (SICRs) in Ni+Al multilayers is demonstrated in Sandia’s shock physics hydrocode, CTH. The approach utilizes Arrhenius-type kinetics and artificial thermal conduction. Our work builds upon previous efforts to parameterize equations of state for Ni x Al y [J. Appl. Phys. 137, 075102 (2025)], as well as simulations of inert shocks in realistic 2D microstructures [J. Appl. Phys. 137, 225301 (2025)]. To calibrate the reaction kinetics, pairs of the reaction coordinate, R′, vs time are extracted from the molecular dynamics (MD) literature. Here, the MD-informed kinetics are used to simulate the dynamic evolution of pressure and temperature in 2D mesoscale simulations. Overall, the MD-informed kinetics obtained for planar interfaces are too slow, as initial reaction is not observed on a nanosecond time scale. Even with quasi-periodic shock focusing leading to the formation of so-called “hot-spots,” the hot spots are unable to grow and coalesce using the fitted Arrhenius rate constants. However, by increasing the rate constants by two orders of magnitude, SICRs are observed at a shock pressure near 30 GPa, which is supported by experiments. Consequently, these mesoscale simulations suggest that unresolved shear-based mechanical mixing might possibly account for the discrepancies in kinetic rates, with shock-generated intense perturbations, interfacial vortical flows, and elevated temperatures serving as favorable reaction conditions. Future work will calibrate a shear-dependent reaction rate from the MD simulations with realistic interfaces that are wavy, diffuse, and disordered.

Kittell, David E. [Sandia National Laboratories (S↗
DOE OSTI2023

Creating large Fock states and massively squeezed states in optics using systems with nonlinear bound states in the continuum

The quantization of the electromagnetic field leads directly to the existence of quantum mechanical states, called Fock states, with an exact integer number of photons. Despite these fundamental states being long-understood, and despite their many potential applications, generating them is largely an open problem. For example, at optical frequencies, it is challenging to deterministically generate Fock states of order two and beyond. Here, we predict the existence of an effect in nonlinear optics, which enables the deterministic generation of large Fock states at arbitrary frequencies. The effect, which we call an n-photon bound state in the continuum, is one in which a photonic resonance (such as a cavity mode) becomes lossless when a precise number of photons n is inside the resonance. Based on analytical theory and numerical simulations, we show that these bound states enable a remarkable phenomenon in which a coherent state of light, when injected into a system supporting this bound state, can spontaneously evolve into a Fock state of a controllable photon number. This effect is also directly applicable for creating (highly) squeezed states of light, whose photon number fluctuations are (far) below the value expected from classical physics (i.e., shot noise). In conclusion, we suggest several examples of systems to experimentally realize the effects predicted here in nonlinear nanophotonic systems, showing examples of generating both optical Fock states with large n (n > 10), as well as more macroscopic photonic states with very large squeezing, with over 90% less noise (10 dB) than the classical value associated with shot noise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Reaction-induced departures from continuum Navier–Stokes turbulence

Reactive hydrodynamic turbulence is an inherently multiscale phenomenon, characterized by the separation between energy-containing, viscous, and molecular length and time scales. The separation between the viscous scale (the Kolmogorov scale) and the molecular mean free path ostensibly justifies a macroscopic description of reactive turbulence via the Navier–Stokes (NS) equations. However, here we use molecular-level simulations to demonstrate that exothermic bimolecular reactions can cause the NS description of turbulence to break down in the near-continuum regime. Sufficiently energetic heat-releasing reactive collisions strongly distort the Maxwell–Boltzmann velocity distribution function, modifying not only the macroscopic chemical rate law but the kinetic-energy-transfer processes as well. This translational nonequilibrium ultimately introduces significant departures from the NS description in the kinetic energy spectra at scales orders of magnitude larger than both the molecular mean free path and the Kolmogorov length scale. These departures prove substantial enough to meaningfully alter integrated quantities, including the overall turbulence kinetic energy itself.

molecular gas dynamics↗
DOE OSTI2022

Self-consistent description of the halo nature of 31 Ne with continuum and pairing correlations

Here a relativistic structure model has previously been used to predict a halo structure for 31 Ne, consistent with halo signatures from measured reaction cross sections of Ne isotopes bombarding Carbon targets. However, previous attempts to calculate those cross sections with reaction models were missing contributions from resonances and pairing correlations in their structure input. Use a reaction model with our relativistic fully microscopic structure model input to predict these cross sections and momentum distributions and analyze for possible halo signatures. Structure input for exotic Ne isotopes were obtained via the analytical continuation of the coupling constant (ACCC) method based on the relativistic mean field (RMF) theory with Bardeen - Cooper - Schrieffer (BCS) pairing approximation, the RAB approach. Total reaction cross sections, one-neutron removal cross sections, and momentum distributions of breakup reaction products were calculated with a Glauber model using our relativistic structure input. Our predictions of total reaction and one-neutron removal cross sections of 31 Ne on a Carbon target were significantly enhanced compared with those of neighboring Neon isotopes, agreeing with measurements at 240 MeV/nucleon and consistent with a single neutron halo. Furthermore, our calculations of the inclusive longitudinal momentum distribution of the 30 Ne and valence neutron residues from the 31 Ne breakup reaction indicate a dilute density distribution in coordinate space, another halo signature. We give a full description of the halo nature of 31 Ne that includes a self-consistent use of pairing and continuum contributions that makes predictions consistent with reaction cross section measurements. This approach can be utilized to determine the halo structure of other exotic nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Lyman continuum leakage from low-mass galaxies with M ⋆ < 108 M⊙

ABSTRACT We present observations with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope of nine low-mass star-forming galaxies at redshifts, z, in the range 0.3179–0.4524, with stellar masses $M_\star \, \lt $ 108 M⊙ and very high specific star-formation rates sSFR ∼150–630 Gyr−1, aiming to study the dependence of leaking Lyman continuum (LyC) emission on stellar mass and some other characteristics of the galaxy. We detect LyC emission in four out of nine galaxies with escape fractions, fesc(LyC), in the range of 11–35 per cent, and establish upper limits for fesc(LyC) in the remaining five galaxies. We observe a narrow Ly α emission line with two peaks in seven galaxies and likely more complex Ly α profiles in the two remaining galaxies. The velocity separation between the peaks Vsep varies in the range from ∼229 to ∼512 km s−1. Our additional data on low-mass galaxies confirm and strengthen the tight anticorrelation between fesc(LyC) and Vsep found for previous low-redshift galaxy samples with higher stellar masses. Vsep remains the best indirect indicator of LyC leakage. It is better than O32 on which fesc(LyC) depends weakly, with a large scatter. Finally, contrary to expectations, we find no increase of fesc(LyC) with decreasing galaxy stellar mass M⋆.

Izotov, Y. I.↗
DOE OSTI2022

Lyman alpha and Lyman continuum emission of Mg Ⅱ -selected star-forming galaxies

We present observations with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope of seven compact low-mass star-forming galaxies at redshifts, z, in the range 0.3161–0.4276, with various O 3 Mg 2 = [O iii] λ5007/Mg ii λ2796+2803 and Mg 2 = Mg ii λ2796/Mg ii λ2803 emission-line ratios. We aim to study the dependence of leaking Lyman continuum (LyC) emission on the characteristics of Mg ii emission together with the dependencies on other indirect indicators of escaping ionizing radiation. LyC emission with escape fractions f esc (LyC) = 3.1–4.6 percent is detected in four galaxies, whereas only 1σ upper limits of f esc (LyC) in the remaining three galaxies were derived. A strong narrow Lyα emission line with two peaks separated by V sep ~ 298–592 km s -1 was observed in four galaxies with detected LyC emission and very weak Lyα emission is observed in galaxies with LyC non-detections. Our new data confirm the tight anticorrelation between f esc (LyC) and V sep found for previous low-redshift galaxy samples. V sep remains the best indirect indicator of LyC leakage among all considered indicators. It is found that escaping LyC emission is detected predominantly in galaxies with Mg 2 ≳ 1.3. A tendency of an increase of f esc (LyC) with increasing of both the O 3 Mg 2 and Mg 2 is possibly present. However, there is substantial scatter in these relations not allowing their use for reliable prediction of f esc (LyC).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Abundances of CNO elements in z ∼ 0.3–0.4 Lyman continuum leaking galaxies

We present observations with the Space Telescope Imaging Spectrograph onboard the Hubble Space Telescope (HST) of 11 Lyman continuum (LyC) leaking galaxies at redshifts, z , in the range 0.29–0.43, with oxygen abundances 12+log(O/H) = 7.64–8.16, stellar masses M ⋆ ∼ 10 7.8 –10 9.8 M ⊙ , and O 32 = [O III ] λ5007/[O II ] λ3727 of ∼ 5–20 aiming to detect the C III ] λ1908 emission line. We combine these observations with the optical Sloan Digital Sky Survey (SDSS) spectra for the determination of carbon, nitrogen, and oxygen abundances. Our sample was supplemented by 31 galaxies from the literature, for which carbon, nitrogen, and oxygen abundances can be derived from the HST and SDSS spectra. These additional galaxies, however, do not have LyC observations. We find that log(C/O) for the entire sample at 12+log(O/H) < 8.1 does not depend on metallicity, with a small dispersion of ∼0.13 dex around the average value of ∼−0.75 dex. On the other hand, the log(N/O) in galaxies at z > 0.1, including LyC leakers, is systematically higher compared to the rest of the sample with lower metallicity. We find that log(C/O) slightly decreases with increasing M ⋆ from ∼ −0.65 at M ⋆ = 10 6 M ⊙ to ∼ −0.80 at M ⋆ = 10 9 –10 10 M ⊙ , whereas log(N/O) is considerably enhanced at M ⋆ > 10 8 M ⊙ . The origin of these trends remains basically unknown. A possible solution would be to assume that the upper mass limit of the stellar initial mass function in more massive galaxies is higher. This would result in a higher production of oxygen and a larger fraction of massive stars with stellar wind polluting the interstellar medium with nitrogen.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Disorder-induced spin excitation continuum and spin-glass ground state in the inverse spinel CuGa 2 O 4

Spinel-structured compounds serve as prototypical examples of highly frustrated systems and are promising candidates for realizing the long-sought quantum spin liquid (QSL) state. However, structural disorder is inevitable in many real QSL candidates and its impact remains a topic of intense debate. In this work, we conduct comprehensive investigations on CuGa 2 ⁢ O 4 , a spinel compound with significant structural disorder, focusing on its thermodynamic properties and spectroscopic behaviors. No long-range magnetic order is observed down to ∼80 mK, as evidenced by magnetic susceptibility, specific-heat, and elastic neutron scattering measurements. More intriguingly, inelastic neutron scattering experiments reveal a broad gapless continuum of magnetic excitations around the Brillouin zone boundary, resembling the magnetic excitation spectra expected for a QSL. Nevertheless, a spin-freezing transition at 𝑇 f ≈ 0.88 K is identified from the cusp in the dc susceptibility curves, where a bifurcation between zero-field-cooling and field-cooling curves occurs. Furthermore, ac susceptibility measurements show a peak close to 𝑇 f at low frequency, which shifts to higher temperature with increasing frequency. These results show that CuGa 2 ⁢ O 4 has a spin-glass ground state, consistent with the establishment of short-range order inferred from the specific-heat measurements. Collectively, these results illustrate the crucial role of disorder in defining the excitation spectrum out of the disordered ground state. Furthermore, our findings shed light onto the broader class of 𝐴⁢𝐵 2 ⁢O 4 spinels and advance our understanding of the spin dynamics in magnetically disordered systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Bound States in the Continuum on a Silicon Chip with Dynamic Tuning

Expanding bound states in the continuum (BIC) beyond photonic crystal systems may enable broader applications benefiting from the unique properties of BIC states. We use photonic integrated circuit to realize a Fabry-Pérot BIC on a silicon chip. The devices consist of cascaded ring resonators with tunable resonance frequencies and phase delays. As a result, the BIC state is dynamically tuned with electrical I/O. We analyze the mechanism of the formation of BIC states in this waveguide system and point out the fundamental differences between the BIC states and electromagnetically induced transparency states in integrated photonics. The high transmission protected by the BIC state enables versatile optical filters, which are capable of independent control over switching, peak position, and quality factor. We also demonstrate the scalability of this platform. As a result, this integrated silicon photonic platform brings opportunities for practical BIC applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Nanostructured $\mathrm{GaAs/(Al,Ga)As}$ Waveguide for Low-Density Polariton Condensation from a Bound State in the Continuum

Exciton-polaritons are hybrid light-matter states that arise from strong coupling between an exciton resonance and a photonic cavity mode. As bosonic excitations, they can undergo a phase transition to a condensed state that can emit coherent light without a population inversion. This aspect makes them good candidates for thresholdless lasers, yet short exciton-polariton lifetime has made it difficult to achieve condensation at very low power densities. In this sense, long-lived symmetry-protected states are excellent candidates to overcome the limitations that arise from the finite mirror reflectivity of monolithic microcavities. In this work we use a photonic symmetry-protected bound state in the continuum coupled to an excitonic resonance to achieve state-of-the-art polariton condensation threshold in a GaAs/(Al,Ga)As waveguide. Most important, we show the influence of fabrication control and how surface passivation via atomic layer deposition provides a way to reduce exciton quenching at the grating sidewalls.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Lattice collective modes from a continuum model of magic-angle twisted bilayer graphene

Here, we show that the insulating states of magic-angle twisted bilayer graphene support a series of collective modes corresponding to local particle-hole excitations on triangular lattice sites. Our theory is based on a continuum model of the magic angle flat bands. When the system is insulating at moiré band filling ν = -3, our calculations show that the ground state supports seven low-energy modes that lie well below the charge gap throughout the moiré Brillouin zone, one of which couples strongly to THz photons. The low-energy collective modes are faithfully described by a model with a local SU(8) degree of freedom in each moiré unit cell that we identify as the direct product of spin, valley, and an orbital pseudospin. Apart from spin and valley-wave modes, the collective mode spectrum includes a low-energy intraflavor exciton mode associated with transitions between flat valence and conduction band orbitals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Boundary theory of the X-cube model in the continuum

We study the boundary theory of the Z N X-cube model using a continuum perspective, from which the exchange statistics of a subset of bulk excitations can be recovered. We discuss various gapped boundary conditions that either preserve or break the translation/rotation symmetries on the boundary, and further present the corresponding ground state degeneracies on T 2 × I. The low-energy physics is highly sensitive to the boundary conditions: even the extensive part of the ground state degeneracy can vary when different sets of boundary conditions are chosen on the two boundaries. Here, we also examine the anomaly inflow of the boundary theory and find that the X-cube model is not the unique (3 + 1)-dimensional theory that cancels the 't Hooft anomaly of the boundary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗