Search NASA⌕ Search

SEARCH · Search NASA

Results for “perturbation”

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

Perturbation of Tropical Stratospheric Ozone Through Homogeneous and Heterogeneous Chemistry Due To Pinatubo

Abstract The Pinatubo eruption in 1991 injected 10–20 Tg SO 2 into the stratosphere, which formed sulfate aerosols through oxidation. Our modeling results show that volcanic heating significantly perturbs the heterogeneous and homogeneous chemistry including NO x and HO x catalytic cycles in the tropical stratosphere. The simulated tropical chemical ozone tendency is positive at 20 mb while negative at 10 mb in the tropics. The simulated ozone chemical tendency is of the same magnitude as the dynamical ozone tendency caused by the accelerated tropical upwelling, but with the opposite sign. Our study finds that the tropical ozone chemical tendency due to homogeneous chemistry becomes more important than heterogeneous chemistry 3 months after eruption. Sensitivity simulations further suggest that the tropical ozone tendency through heterogeneous chemistry is saturated when the injected amount exceeds 2 Tg.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Assessing Impacts of Environmental Perturbations on Urban Biogenic Carbon Exchange in the Chicago Region

Carbon dioxide (CO 2 ) quantification is critical for assessing city-level carbon emissions and sustainable urban development. While urban vegetation has the potential to provide environmental benefits, such as heat and carbon mitigation, the CO 2 exchange from biogenic sectors and its impact from the environmental perturbations are often overlooked. It is also challenging to simulate the plant functions in the complex urban terrain. This study presents a processed-based modeling approach to assess the biogenic carbon fluxes from the vegetated areas over the Chicago Metropolitan Area (CMA) using the Weather Research and Forecast—Urban Biogenic Carbon exchange model. We investigate the change of CO 2 sink power in CMA under heatwaves and irrigation. The results indicate that the vegetation plays a significant role in the city's carbon portfolio and the landscaping management has the potential to reduce carbon emissions significantly. Furthermore, based on the competing mechanisms in the biogenic carbon balance identified in this study, we develop a novel Environmental Benefit Score metrics framework to identify the vulnerability and mitigation measures associated with nature-based solutions (NbS) within CMA. By using the generalized portable framework and our science-policy confluence analysis presented in this study, global cities can maximize the effectiveness of NbS and accelerate carbon neutrality.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

4D Electrical Resistivity Imaging of Stress Perturbations Induced During High–Pressure Shear Stimulation Tests

Fluid flow through fractured media is typically governed by the distribution of fracture apertures, which are in turn governed by stress. Consequently, understanding subsurface stress is critical for understanding and predicting subsurface fluid flow. Although laboratory–scale studies have established a sensitive relationship between effective stress and bulk electrical conductivity in crystalline rock, that relationship has not been extensively leveraged to monitor stress evolution at the field scale using electrical or electromagnetic geophysical monitoring approaches. In this paper we demonstrate the use time–lapse 3–dimensional (4D) electrical resistivity tomography to image perturbations in the stress field generated by pressurized borehole packers deployed during shear–stimulation attempts in a 1.25 km deep metamorphic crystalline rock formation.

58 GEOSCIENCES↗
DOE OSTI2025

Relationships Between Mesoscale Convective System Properties and Midlevel Dynamic Perturbations

Abstract Past studies implicate dynamic anomalies operating on subsynoptic scales as a possible initiation source of summertime (July–August) mesoscale convective systems (MCSs) in the central United States during northwesterly flow regimes. To improve our understanding of warm season MCSs occurring over a variety of flow regimes, we track midlevel (600 hPa) vorticity perturbations (“MPs”) as 2D objects comprising wavelengths of 500–2,500 km over the central US from May–August of 2004–2021. We perform statistical analysis of relationships between metrics of MP objects (e.g., duration, size, intensity, and origin) and high‐resolution MCS precipitation characteristics (e.g., duration, total rainfall, rain coverage area, and motion) that occur while collocated with or in the absence of MPs to discern predictive capability of background dynamic features on storm precipitation potential. Although the majority of MPs collocated with MCS initiation occur during July–August, a significant number (40%) occur between May and June. Northwesterly flow MPs comprise a relative minority of our events, suggesting that MPs can affect MCSs across a variety of warm season flow regimes. MPs affecting MCSs initiated primarily over the high plains near the central Rockies. Only approximately 20% of tracked MCS initiation events were collocated with MPs, but these storms produced ∼25% greater lifetime rainfall and coverage area, and ∼29% more stratiform rain than non‐MP‐induced MCSs. In general, larger and more vigorous MPs resulted in more hydrologically impactful MCSs. The most directly attributable benefit to MCS initiation was from MP‐enhanced background vertical motion and thermodynamic instability (e.g., increased CAPE).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Anisotropic leaky-like perturbation with subwavelength gratings enables zero crosstalk

Abstract Electromagnetic coupling via an evanescent field or radiative wave is a primary characteristic of light, allowing optical signal/power transfer in a photonic circuit but limiting integration density. A leaky mode, which combines both evanescent field and radiative wave, causes stronger coupling and is thus considered not ideal for dense integration. Here we show that a leaky oscillation with anisotropic perturbation rather can achieve completely zero crosstalk realized by subwavelength grating (SWG) metamaterials. The oscillating fields in the SWGs enable coupling coefficients in each direction to counteract each other, resulting in completely zero crosstalk. We experimentally demonstrate such an extraordinarily low coupling between closely spaced identical leaky SWG waveguides, suppressing the crosstalk by ≈40 dB compared to conventional strip waveguides, corresponding to ≈100 times longer coupling length. This leaky-SWG suppresses the crosstalk of transverse–magnetic (TM) mode, which is challenging due to its low confinement, and marks a novel approach in electromagnetic coupling applicable to other spectral regimes and generic devices.

42 ENGINEERING↗
DOE OSTI2020

Ab initio electron-two-phonon scattering in GaAs from next-to-leading order perturbation theory

Electron-phonon (e–ph) interactions are usually treated in the lowest order of perturbation theory. Here we derive next-to-leading order e–ph interactions, and compute from first principles the associated electron-two-phonon (2ph) scattering rates. The derivations involve Matsubara sums of two-loop Feynman diagrams, and the numerical calculations are challenging as they involve Brillouin zone integrals over two crystal momenta and depend critically on the intermediate state lifetimes. Using Monte Carlo integration together with a self-consistent update of the intermediate state lifetimes, we compute and converge the 2ph scattering rates, and analyze their energy and temperature dependence. We apply our method to GaAs, a weakly polar semiconductor with dominant optical-mode long-range e–ph interactions. We find that the 2ph scattering rates are as large as nearly half the value of the one-phonon rates, and that including the 2ph processes is necessary to accurately predict the electron mobility in GaAs from first principles.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Machine learning the electric field response of condensed phase systems using perturbed neural network potentials

Abstract The interaction of condensed phase systems with external electric fields is of major importance in a myriad of processes in nature and technology, ranging from the field-directed motion of cells (galvanotaxis), to geochemistry and the formation of ice phases on planets, to field-directed chemical catalysis and energy storage and conversion systems including supercapacitors, batteries and solar cells. Molecular simulation in the presence of electric fields would give important atomistic insight into these processes but applications of the most accurate methods such as ab-initio molecular dynamics (AIMD) are limited in scope by their computational expense. Here we introduce Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) to push back the accessible time and length scales of such simulations. We demonstrate that important dielectric properties of liquid water including the field-induced relaxation dynamics, the dielectric constant and the field-dependent IR spectrum can be machine learned up to surprisingly high field strengths of about 0.2 V Å −1 without loss in accuracy when compared to ab-initio molecular dynamics. This is remarkable because, in contrast to most previous approaches, the two neural networks on which PNNP MD is based are exclusively trained on molecular configurations sampled from zero-field MD simulations, demonstrating that the networks not only interpolate but also reliably extrapolate the field response. PNNP MD is based on rigorous theory yet it is simple, general, modular, and systematically improvable allowing us to obtain atomistic insight into the interaction of a wide range of condensed phase systems with external electric fields.

Science & Technology - Other Topics↗
DOE OSTI2024

Second-harmonic generation tensors from high-throughput density-functional perturbation theory

Optical materials play a key role in enabling modern optoelectronic technologies in a wide variety of domains such as the medical or the energy sector. Among them, nonlinear optical crystals are of primary importance to achieve a broader range of electromagnetic waves in the devices. However, numerous and contradicting requirements significantly limit the discovery of new potential candidates, which, in turn, hinders the technological development. In the present work, the static nonlinear susceptibility and dielectric tensor are computed via density-functional perturbation theory for a set of 579 inorganic semiconductors. The computational methodology is discussed and the provided database is described with respect to both its data distribution and its format. Several comparisons with both experimental and ab initio results from literature allow to confirm the reliability of our data. The aim of this work is to provide a relevant dataset to foster the identification of promising nonlinear optical crystals in order to motivate their subsequent experimental investigation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Large-scale application of free energy perturbation calculations for antibody design

Abstract Alchemical free energy perturbation (FEP) is a rigorous and powerful technique to calculate the free energy difference between distinct chemical systems. Here we report our implementation of automated large-scale FEP calculations, using the Amber software package, to facilitate antibody design and evaluation. In combination with Hamiltonian replica exchange, our FEP simulations aim to predict the effect of mutations on both the binding affinity and the structural stability. Importantly, we incorporate multiple strategies to faithfully estimate the statistical uncertainties in the FEP results. As a case study, we apply our protocols to systematically evaluate variants of the m396 antibody for their conformational stability and their binding affinity to the spike proteins of SARS-CoV-1 and SARS-CoV-2. By properly adjusting relevant parameters, the particle collapse problems in the FEP simulations are avoided. Furthermore, large statistical errors in a small fraction of the FEP calculations are effectively reduced by extending the sampling, such that acceptable statistical uncertainties are achieved for the vast majority of the cases with a modest total computational cost. Finally, our predicted conformational stability for the m396 variants is qualitatively consistent with the experimentally measured melting temperatures. Our work thus demonstrates the applicability of FEP in computational antibody design.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

A non-perturbative pairwise-additive analysis of charge transfer contributions to intermolecular interaction energies

Energy decomposition analysis (EDA) based on absolutely localized molecular orbitals (ALMOs) decomposes the interaction energy between molecules into physically interpretable components like geometry distortion, frozen interactions, polarization, and charge transfer (CT, also sometimes called charge delocalization) interactions. In this work, a numerically exact scheme to decompose the CT interaction energy into pairwise additive terms is introduced for the ALMO-EDA using density functional theory. Unlike perturbative pairwise charge-decomposition analysis, the new approach does not break down for strongly interacting systems, or show significant exchange–correlation functional dependence in the decomposed energy components. Both the energy lowering and the charge flow associated with CT can be decomposed. Complementary occupied–virtual orbital pairs (COVPs) that capture the dominant donor and acceptor CT orbitals are obtained for the new decomposition. It is applied to systems with different types of interactions including DNA base-pairs, borane-ammonia adducts, and transition metal hexacarbonyls. While consistent with most existing understanding of the nature of CT in these systems, we report the results also reveal some new insights into the origin of trends in donor–acceptor interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Improving the accuracy of the variational quantum eigensolver for molecular systems by the explicitly-correlated perturbative [2] R12 -correction

We provide an integration of the universal, perturbative explicitly correlated [2] R12 -correction in the context of the Variational Quantum Eigensolver (VQE). This approach is able to increase the accuracy of the underlying reference method significantly while requiring no additional quantum resources. The proposed approach only requires knowledge of the one- and two-particle reduced density matrices (RDMs) of the reference wavefunction; these can be measured after having reached convergence in the VQE. This computation comes at a cost that scales as the sixth power of the number of electrons. Here, we explore the performance of the VQE + [2] R12 approach using both conventional Gaussian basis sets and our recently proposed directly determined pair-natural orbitals obtained by multiresolution analysis (MRA-PNOs). Both Gaussian orbital and PNOs are investigated as a potential set of complementary basis functions in the computation of [2] R12 . In particular the combination of MRA-PNOs with [2] R12 has turned out to be very promising – persistently throughout our data, this allowed very accurate simulations at a quantum cost of a minimal basis set. Additionally, we found that the deployment of PNOs as complementary basis can greatly reduce the number of complementary basis functions that enter the computation of the correction at a complexity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Polymorphism and phase transitions in Na 2 U 2 O 7 from density functional perturbation theory

Polymorphism and phase transitions in sodium diuranate, Na 2 U 2 O 7 , are investigated with density functional perturbation theory (DFPT). Thermal properties of crystalline α-, β- and γ-Na 2 U 2 O 7 polymorphs are predicted from DFPT phonon calculations, i.e., the first time for the high-temperature γ-Na 2 U 2 O 7 phase (R$\bar{3}$ with combining macron]m symmetry). The standard molar isochoric heat capacities predicted within the quasi-harmonic approximation are C p (298.15K)= 219.4 and 220.9 J K -1 mol -1 mfor P21/a α-Na 2 U 2 O 7 and C2/m β-Na 2 U 2 O 7 , respectively. Gibbs free energy calculations reveal that α-Na 2 U 2 O 7 (P2 1 /a) and β-Na 2 U 2 O 7 (C2/m) are almost energetically degenerate at low temperature, with β-Na 2 U 2 O 7 becoming slightly more stable than α-Na 2 U 2 O 7 as temperature increases. Here these findings are consistent with XRD data showing a mixture of α and β phases after cooling of γ-Na 2 U 2 O 7 to room temperature and the observation of a sluggish α → β phase transition above ca. 600 K. A recently observed α-Na 2 U 2 O 7 structure with P21 symmetry is also shown to be metastable at low temperature. Based on Gibbs free energy, no direct β → γ solid-solid phase transition is predicted at high temperature, although some experiments reported the existence of such phase transition around 1348 K. This, along with recent experiments, suggests the occurrence of a multi-step process consisting of initial β-phase decomposition, followed by recrystallization into γ-phase as temperature increases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Low-cost calculation and analysis of 2D IR spectra of model diiron trinitrosyl complexes in the NO stretch region with vibrational perturbation theory

Two-dimensional infrared spectroscopy offers unique capabilities for probing vibrational coupling in complex metal–ligand systems. In this paper, we combine two-dimensional infrared spectroscopy with vibrational perturbation theory to investigate vibrational coupling in a diiron trinitrosyl complex across three stable redox states. Although these systems are challenging for electronic structure methods, we demonstrate that key features of experimental 2D IR spectra can be accurately reproduced using reduced-dimensional anharmonic calculations with a small harmonic frequency scaling. Analysis reveals that N–O stretching modes maintain high locality across all redox states, with coupling patterns that directly reflect variations in Fe–N bond strength. Using curvilinear coordinate analysis, we demonstrate these differences result from systematic changes in cubic anharmonic force constants rather than mode delocalization. Our results establish N–O stretches as sensitive probes of metal–ligand bonding strength, expanding the toolkit for studying biologically relevant nitrosyl complexes.

Moran, Hayden A. [Texas A&M University, College St↗
DOE OSTI2026

Perturbation of nanoplastics on biomembranes: molecular insights from neutron scattering

Plastic waste is now pervasive in the environment, breaking down into microplastics and nanoplastics under many environmental conditions. These particles have been found in various ecosystems and even in human tissues, raising significant environmental and health concerns. In this study, we investigated the interaction of polystyrene nanoplastics, with and without surface modifications, on biomembrane structures using contrast-matching small-angle neutron scattering and neutron spin echo spectroscopy. The neutron contrast matching enabled the selective study of biomembranes in the presence of nanoplastics. Two model membranes were employed: a simple zwitterionic bilayer (i.e., dimyristoylphosphatidylcholine [DMPC]) and an Escherichia coli lipid extract as a bacterial membrane model. The results show profound membrane disruptions, including possible thinning, vesicle fragmentation, lipid monolayer formation, and inter-vesicle aggregation, with the more severe effects observed in DMPC membranes. Notably, E. coli membranes exhibited greater resilience, suggesting that natural membranes with diverse lipid compositions may reduce susceptibility to perturbation by extracellular nanoplastics. Here, these findings highlight potential risks posed by environmental nanoplastic particles to biological membranes, with insights into molecular-level interactions and the environmental toxicity of nanoplastics. This work provides a foundation for future studies into nanoplastic–biomembrane interactions and their broader implications for health and environment using neutrons.

Qian, Shuo [Oak Ridge National Laboratory (ORNL), ↗
DOE OSTI2025

Topological perturbation to a standard dehydrogenation catalyst, Pt 3 Sn

Topological materials, which exhibit protected topological surface states (TSS) near the Fermi level, have been proposed to be good catalysts. Topological catalysis may be more prevalent than we suspect, and not limited to exotic new materials. Here we study a known dehydrogenation catalyst, Pt 3 Sn alloy, which happens to be a topological semimetal, and probe the participation of TSSs in catalytic dehydrogenation of methane catalyzed by this material. Through first principle modeling and detailed analysis of the electronic structure for topological and non-topological surfaces of Pt 3 Sn, we find that TSS get significantly altered by the binding of reaction intermediates, particularly H. However, this effect of TSS on the binding of the reagents is merely perturbative, as the majority of the adsorbate binding is achieved by not-surface-focused electronic states, located much deeper below the Fermi level. Therefore, the reaction energetics and selectivity are predominantly determined by electronic states other than TSS. The fact that TSS are available for the reagent binding does not alone guarantee that the catalysis is strongly driven by TSS. However, TSS are not to be ignored, as small changes in the energetics along the reaction profile can translate into substantial differences in the reaction rate. Hence, in our view, Pt 3 Sn – a topological material – is first and foremost a standard catalyst, with added topological features, and not purely a topological catalyst. Our results point at the need to carefully consider all the bonding effects at the topological material interface.

Chemistry↗
DOE OSTI2025

Beyond the two-point correlation: Constraining primordial non-Gaussianity with density-perturbation moments

Constraining primordial non-Gaussianity (PNG) on the large-scale cosmic structure (LSS) is an important step in understanding properties of the early Universe, specifically in distinguishing between different inflationary models. Measuring PNG relies on evaluating the scale-dependent correlations in the density field. New summary statistics beyond the two- and three-point correlation functions in configuration space and their Fourier-space counterparts, the power- and bispectrum may provide increased sensitivity. We introduce a new method for extracting the PNG signal imprinted on the LSS by using the first three Gaussian moments of the normalized correlation in density perturbations, evaluated on varying distance scales. We aim to assess this method’s sensitivity to local PNG, parameterized by f NL . We performed spherical convolutions on a range of scales on dark-matter-halo simulations to measure the scale-dependent correlations in the density field. From these, we computed the first three moments and compared them to a model expectation vector, parameterized to the second power in f NL . Our method provides about 21% improvement in sensitivity to f NL with respect to using the two-point correlation function alone. Notably, we find that the second moment alone carries nearly as much constraining power as the mean, highlighting the potential of higher order statistics. Given its simplicity and efficiency, this framework is well suited for application to current and upcoming large-scale surveys such as the Dark Energy Spectroscopic Instrument (DESI).

early universe↗
DOE OSTI2020

Gyrokinetic understanding of the edge pedestal transport driven by resonant magnetic perturbations in a realistic divertor geometry

Self-consistent simulations of neoclassical and electrostatic turbulent transport in a DIII-D H-mode edge plasma under resonant magnetic perturbations (RMPs) have been performed using the global total-f gyrokinetic particle-in-cell code x-point gyrokinetic code (XGC), in order to study density pump-out and electron heat confinement. The RMP field is imported from the extended magneto-hydrodynamics code M3D-C1, taking into account the linear two-fluid plasma response. With both neoclassical and turbulence physics considered together, the XGC simulation reproduces two key features of experimentally observed edge transport under RMPs: increased radial particle transport in the pedestal region that is sufficient to account for the experimental pump-out rate and suppression of the electron heat flux in the steepest part of the edge pedestal. In the simulation, the density fluctuation amplitude of modes moving in the electron diamagnetic direction increases due to interaction with RMPs in the pedestal shoulder and outward, while the electron temperature fluctuation amplitude decreases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Theory of edge localized mode suppression by static resonant magnetic perturbations in the DIII-D tokamak

According to a recent paper [Hu et al ., Phys. Plasmas 26, 120702 (2019)], mode penetration at the top of the pedestal is a necessary and sufficient condition for the suppression of edge localized modes (ELMs) by a resonant magnetic perturbation (RMPs) in an H-mode tokamak discharge. This paper employs asymptotic matching theory to model a particular DIII-D discharge in which ELMs were suppressed by an externally generated, static, n ¼ 2, RMP whose amplitude was modulated at a frequency of 1 Hz. It is demonstrated that the response of the plasma to the applied RMP, in the immediate vicinities of the rational (i.e., resonant) surfaces, is governed by nonlinear, rather than by linear, physics. This is the case because the magnetic island widths associated with driven reconnection exceed the linear layer widths, even in cases where driven reconnection is strongly suppressed by plasma rotation. The natural frequency at a given rational surface (i.e., the helical frequency at which the locally resonant component of the RMP would need to propagate in order to maximize driven reconnection) is found to be offset from the local E B frame in the ion diamagnetic direction. The size of the offset is mostly determined by neoclassical poloidal rotation. Finally, the predictions of a fully nonlinear plasma response model are found to be broadly consistent with the DIII-D experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗