Search NASA⌕ Search

SEARCH · Search NASA

Results for “perturbation test”

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 OSTI2022

Bacterial Membranes Are More Perturbed by the Asymmetric Versus Symmetric Loading of Amphiphilic Molecules

Characterizing the biophysical properties of bacterial membranes is critical for understanding the protective nature of the microbial envelope, interaction of biological membranes with exogenous materials, and designing new antibacterial agents. Presented here are molecular dynamics simulations for two cationic quaternary ammonium compounds, and the anionic and nonionic form of a fatty acid molecule interacting with a Staphylococcus aureus bacterial inner membrane. The effect of the tested materials on the properties of the model membranes are evaluated with respect to various structural properties such as the lateral pressure profile, lipid tail order parameter, and the bilayer’s electrostatic potential. Conducting asymmetric loading of molecules in only one leaflet, it was observed that anionic and cationic amphiphiles have a large impact on the Staphylococcus aureus membrane’s electrostatic potential and lateral pressure profile as compared to a symmetric distribution. Nonintuitively, we find that the cationic and anionic molecules induce a similar change in the electrostatic potential, which points to the complexity of membrane interfaces, and how asymmetry can induce biophysical consequences. Finally, we link changes in membrane structure to the rate of electroporation for the membranes, and again find a crucial impact of introducing asymmetry to the system. Understanding these physical mechanisms provides critical insights and viable pathways for the rational design of membrane-active molecules, where controlling the localization is key.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

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

Drift orbit islands of energetic particles due to 3D fields in ITER

In this study, we report a comparative study of the orbit islands formed by the guiding center drift motion of test energetic particles (EPs), including both high-energy deuterium ions and fusion-born α-particles, with respect to the magnetic islands formed by tracing the field lines, for an ITER plasma representing the 15 MA baseline scenario. The particle drift orbit is modified by the presence of static resonant magnetic perturbations (RMPs) of different toroidal mode numbers n (n = 1 and 3 in this study), forming orbit islands in the Poincaré plane by passing EPs. The key findings are as follows. (i) With an n = 1 RMP field, the size of the orbit islands combined with the total RMP field, including the plasma response, is about three times smaller than that obtained by assuming the corresponding vacuum field. (ii) The orbit island size is not sensitive to the EP energy, and is comparable to those of the magnetic islands. (iii) Passing EPs with outward radial orbital drift form orbit islands that shift inward with respect to the corresponding magnetic islands, and vice versa. This radial shift, ΔΨ, measured in the normalized equilibrium poloidal flux, is quantified as ΔΨ ≃ 0.03X for the ITER baseline plasma and assuming n = 1 RMPs, with $X\equiv {(M/{M}_{\mathrm{p}})}^{1/2}{(E[\mathrm{M}\mathrm{e}\mathrm{V}])}^{1/2}{Z}^{-1}$, where M is the particle mass, Mp is the proton mass, E is the particle energy, and Z is the particle charge number. (iv) Trapped EPs do not form drift orbit islands, even in the presence of the 3D fields produced by the 90 kAt ELM control coil current in ITER, and thus possess good confinement properties. Nevertheless, the deformation of the trapped particle drift orbits in the Poincaré plane shows that the canonical toroidal angular momentum of EPs is no longer conserved in the presence of the RMP fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Direct nonparametric multimessenger constraints on the equation of state of cold dense nuclear matter

We utilize the now substantial amount of astrophysical observations of neutron stars (NSs), along with perturbative quantum chromodynamics (pQCD) calculations at high density, to directly constrain the NS equation of state (EOS). To this end, we construct nonparametric EOS priors by using Gaussian processes trained on 75 EOSs, which include models with either hadrons, hyperons, or quarks at high densities. We create a prior using the full EOS sample (model agnostic), and one prior for each EOS family to test model discrimination. We introduce a novel inference approach, which allows the simultaneous sampling of intrinsic and extrinsic parameters of binary NS mergers, as well as a nonparametric equation of state. We showcase this method in a Bayesian updating scheme by first performing a complete analysis of the binary NS merger event GW170817 with minimal assumptions, and sequentially adding information from x-ray and radio NS observations, along with pQCD calculations. Besides providing standard constraints, such as the pressure at twice nuclear saturation density 𝑝⁡(2⁢𝜌 sat ) = 4.3$^{+0.6}_{−0.6}$ × 10 34 dyne/cm 2 , at 95% confidence level, for the model agnostic prior, our methodology shows how the choice of EOS families used in conditioning changes the inferred astrophysical properties of the EOS, namely tidal deformability and maximum supported NS mass. We find hyperonic priors predicting higher tidal deformabilities for a 1.4⁢𝑀 ⊙ NS, and hadronic priors being preferred by the considered astrophysical data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Ice Phase Classification Made Easy with Score-Based Denoising

Accurate identification of ice phases is essential for understanding various physicochemical phenomena. However, such classification for structures simulated with molecular dynamics is complicated by the complex symmetries of ice polymorphs and thermal fluctuations. For this purpose, both traditional order parameters and data-driven machine learning approaches have been employed, but they often rely on expert intuition, specific geometric information, or large training data sets. In this work, we present an unsupervised phase classification framework that combines a score-based denoiser model with a subsequent model-free classification method to accurately identify ice phases. Further, the denoiser model is trained on perturbed synthetic data of ideal reference structures, eliminating the need for large data sets and labeling efforts. The classification step utilizes the smooth overlap of atomic position (SOAP) descriptors as the atomic fingerprint, ensuring Euclidean symmetries and transferability to various structural systems. Our approach achieves a remarkable 100% accuracy in distinguishing ice phases of test trajectories using only seven ideal reference structures of ice phases as model inputs. This demonstrates the generalizability of the score-based denoiser model in facilitating phase identification for complex molecular systems. The proposed classification strategy can be broadly applied to investigate structural evolution and phase identification for a wide range of materials, offering new insights into the fundamental understanding of water and other complex systems.

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

AP-Net: An atomic-pairwise neural network for smooth and transferable interaction potentials

Intermolecular interactions are critical to many chemical phenomena, but their accurate computation using ab initio methods is often limited by computational cost. The recent emergence of machine learning (ML) potentials may be a promising alternative. Useful ML models should not only estimate accurate interaction energies but also predict smooth and asymptotically correct potential energy surfaces. However, existing ML models are not guaranteed to obey these constraints. Indeed, systemic deficiencies are apparent in the predictions of our previous hydrogen-bond model as well as the popular ANI-1X model, which we attribute to the use of an atomic energy partition. As a solution, we propose an alternative atomic-pairwise framework specifically for intermolecular ML potentials, and we introduce AP-Net—a neural network model for interaction energies. The AP-Net model is developed using this physically motivated atomic-pairwise paradigm and also exploits the interpretability of symmetry adapted perturbation theory (SAPT). We show that in contrast to other models, AP-Net produces smooth, physically meaningful intermolecular potentials exhibiting correct asymptotic behavior. Initially trained on only a limited number of mostly hydrogen-bonded dimers, AP-Net makes accurate predictions across the chemically diverse S66x8 dataset, demonstrating significant transferability. On a test set including experimental hydrogen-bonded dimers, AP-Net predicts total interaction energies with a mean absolute error of 0.37 kcal mol−1, reducing errors by a factor of 2–5 across SAPT components from previous neural network potentials. The pairwise interaction energies of the model are physically interpretable, and an investigation of predicted electrostatic energies suggests that the model “learns” the physics of hydrogen-bonded interactions.

Glick, Zachary L. (ORCID:0000000309002849)↗
DOE OSTI2023

An in-silico NMR laboratory for nuclear magnetic shieldings computed via finite fields: Exploring nucleus-specific renormalizations of MP2 and MP3

Here we developed and implemented a method-independent, fully numerical, finite difference approach to calculating nuclear magnetic resonance shieldings, using gauge-including atomic orbitals. The resulting capability can be used to explore non-standard methods, given only the energy as a function of finite-applied magnetic fields and nuclear spins. For example, standard second-order Møller-Plesset theory (MP2) has well-known efficacy for 1H and 13C shieldings and known limitations for other nuclei such as 15N and 17O. It is, therefore, interesting to seek methods that offer good accuracy for 15N and 17O shieldings without greatly increased compute costs, as well as exploring whether such methods can further improve 1H and 13C shieldings. Using a small molecule test set of 28 species, we assessed two alternatives: κ regularized MP2 (κ-MP2), which provides energy-dependent damping of large amplitudes, and MP2.X, which includes a variable fraction, X, of third-order correlation (MP3). The aug-cc-pVTZ basis was used, and coupled cluster with singles and doubles and perturbative triples [CCSD(T)] results were taken as reference values. Our κ-MP2 results reveal significant improvements over MP2 for 13C and 15N, with the optimal κ value being element-specific. κ-MP2 with κ = 2 offers a 30% rms error reduction over MP2. For 15N, κ-MP2 with κ = 1.1 provides a 90% error reduction vs MP2 and a 60% error reduction vs CCSD. On the other hand, MP2.X with a scaling factor of 0.6 outperformed CCSD for all heavy nuclei. These results can be understood as providing renormalization of doubles amplitudes to partially account for neglected triple and higher substitutions and offer promising opportunities for future applications.

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

Bootstrapping more QM systems

Abstract We test the bootstrap approach for determining the spectrum of one dimensional Hamiltonians. In this paper we focus on problems that have a two parameter search space in the bootstrap approach: the double well and a periodic potential associated with the Mathieu equation. For the double well, we compare the bootstrapped energies with contributions from perturbative and non-perturbative results, finding good agreement. For the periodic potentials, we notice that the bootstrap approach gives the band structure of the periodic potential, but it has trouble finding the quasi-momentum of the system. To make further progress on determining the dispersion relation of the bands, new techniques are needed.

Physics↗
DOE OSTI2020

Hadronic structure in high-energy collisions

Here, Parton distribution functions (PDFs) describe the structure of hadrons as composed of quarks and gluons. They are needed to make predictions for short-distance processes in high-energy collisions and are determined by fitting to cross-section data. Definitions of the PDFs and their relations to high-energy cross sections are reviewed. The focus is on the PDFs in protons, but PDFs in nuclei are also discussed. The standard statistical treatment needed to fit the PDFs to data using the Hessian method is reviewed in some detail. Tests are discussed that critically examine whether the needed assumptions are indeed valid. Also presented are some ideas of what one can do in case tests indicate that the assumptions fail.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Impact parameter dependence of color charge correlations in the proton

The impact parameter dependence of color charge correlators in the proton is obtained from the light front formalism in light cone gauge. We include NLO corrections due to the |qqqg\rangle | q q q g ⟩ Fock state via light-cone perturbation theory. Near the center of the proton, the b b -dependence of the correlations is very different from a ``transverse profile function’’. The resulting t t -dependence of exclusive J/\Psi J / Ψ photoproduction transitions from exponential to power law at |t| \approx 1 | t | ≈ 1 ~GeV ^2 2 . This prediction could be tested at upcoming DIS facilities or in nucleus-proton ultraperipheral collisions (UPCs).

Dumitru, Adrian↗
DOE OSTI2022

A novel flow cell for optical particle analyzers—application to measurements of Malvern Insitec under high pressure and temperature

A flow cell is a necessary measurement interface for some important optical analyzers. In our application of the flow cell, we utilize a state-of-the-art optical analyzer (Malvern Insitec) to measure particle size distribution and concentration in-situ in a sampled flow from a novel pressurized oxygen-fuel combustion process. However, since this sampling flow is a flow of moist flue gas under high temperature and pressure, and the flow contains particles and corrosive acid gases, it is an extreme challenge to obtain a flow cell with a high optical quality that does not perturb the measurement. To address this challenge, we propose a new design for an optical flow cell. By using a unique flow field in the proposed flow cell, the measurement zone can be well defined by the sampling flow, minimizing the influence of purge flow. To demonstrate this flow cell, we have built a test system, and conduct measurements utilizing polydisperse-particle standards (10-100 µm and 1-10 µm). The results reveal that the optical windows are well protected by the purge flow field, without risk of deposition from the sampling flow, and the Malvern Insitec can measure the particle size distribution by using this flow cell, without generating sample bias.

Cheng, Mao↗
DOE OSTI2020

Applying Noether's theorem to matter in the Milky Way: evidence for external perturbations and non-steady-state effects from Gaia Data Release 2

We apply Noether's theorem to observations of main-sequence stars from the Gaia Data Release 2 archive to probe the matter distribution function of the Galaxy. That is, we examine the axial symmetry of stars at vertical heights $z$, $0.2 \le |z| \le 3$ kpc, to probe the quality of the angular momentum $L_z$ as an integral of motion. The failure of this symmetry test would speak to a Milky Way, in both its visible and dark matter, that is not isolated and/or not in steady state. The left-right symmetry-breaking pattern we have observed, north and south, reveals both effects, with a measured deviation from symmetry of typically 0.5%. We show that a prolate form of the gravitational distortion of the Milky Way by the Large Magellanic Cloud, determined from fits to the Orphan stream by Erkal et al., 2019, is compatible with the size and sign of the axial-symmetry-breaking effects we have discovered in our sample of up to 14.4 million main-sequence stars, speaking to a distortion of an emergent, rather than static, nature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Quantifying heavy quark transport coefficients with an improved transport model

We report that the heavy-flavor transport coefficients contain important information on the strong interaction at finite temperatures. The extraction of these numbers from experimental data requires dynamical modeling of heavy-flavor transport that is coupled to realistic medium evolution. Furthermore, meaningful extractions necessitate both a faithful implementation of the physical inputs to be tested and the quantification of model uncertainty. For these purposes, we have developed a partonic transport model LIDO. It has an improved treatment of in-medium parton bremsstrahlung, which has been calibrated to theoretical calculations in a simple medium to reduce modeling uncertainty. Regarding the interaction between heavy quark and the medium, few-body perturbative scatterings are applied to large-momentum transfer (q) processes, while a diffusion equation models the dynamics of small-q processes. Such a separation restricts the explicit use of medium quasi-particles to large-q processes only. Another advantage is that deviations from the leading-order probe-medium coupling can be parametrized as an additional contribution to the diffusion constant. The heavy quark transport coefficients are then extracted with uncertainty estimation from a Bayesian analysis including both the RHIC and the LHC data. The results are found to be consistent with earlier extraction of the light-quark transport coefficients at high momentum and be comparable with lattice calculations of the heavy-flavor diffusion constant in the static limit at low momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Electron beam characterization via quantum coherent optical magnetometry

We present a quantum optics-based detection method for determining the position and current of an electron beam. As electrons pass through a dilute vapor of rubidium atoms, their magnetic field perturbs the atomic spin's quantum state and causes polarization rotation of a laser resonant with an optical transition of the atoms. By measuring the polarization rotation angle across the laser beam, we recreate a 2D projection of the magnetic field and use it to determine the e-beam position, size, and total current. We tested this method for an e-beam with currents ranging from 30 to 110 μA. Our approach is insensitive to electron kinetic energy, and we confirmed that experimentally between 10 and 20 keV. In conclusion, this technique offers a unique platform for noninvasive characterization of charged particle beams used in accelerators for particle and nuclear physics research.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Development Of Non-Invasive Beam Diagnostics By Quantum Optics-Based Detection

We report our effort and latest progress developing diagnostics using quantum optics-based detection method for determining the spatial properties and current of electron beams. As electrons pass through a dilute vapor of rubidium atoms, their magnetic field perturb the atomic spin’s quantum state and causes polarization rotation of a laser resonant with an optical transition of the atoms. By measuring the polarization rotation angle across the laser beam, we recreate a 2D projection of the magnetic field and use it to determine the e-beam position, size and total current. We tested this method for an e-beam with currents ranging from 30 to 110 µA. Our study shows this approach is insensitive to electron kinetic energy, and is initially verified experimentally between 10 to 20 keV. A different approach was also performed with highly excited rubidium atoms in Rydberg states through a coherent two photon process. We use this detection scheme to measure the electric field generated by a

Zhang, S. [Thomas Jefferson National Accelerator F↗
DOE OSTI2020

Vacuum stability in inert higgs doublet model with right-handed neutrinos

We analyze the vacuum stability in the inert Higgs doublet extension of the Standard Model (SM), augmented by right-handed neutrinos (RHNs) to explain neutrino masses at tree level by the seesaw mechanism. We make a comparative study of the high- and low-scale seesaw scenarios and the effect of the Dirac neutrino Yukawa couplings on the stability of the Higgs potential. Bounds on the scalar quartic couplings and Dirac Yukawa couplings are obtained from vacuum stability and perturbativity considerations. These bounds are found to be relevant only for low-scale seesaw scenarios with relatively large Yukawa couplings. The regions corresponding to stability, metastability and instability of the electroweak vacuum are identified. These theoretical constraints give a very predictive parameter space for the couplings and masses of the new scalars and RHNs which can be tested at the LHC and future colliders. The lightest non-SM neutral CP-even/odd scalar can be a good dark matter candidate and the corresponding collider signatures are also predicted for the model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Spatiotemporal 4D Whole-cell Modeling of a Minimal Autotroph Reveals Central Carbon Metabolism Regulated Locally by Protein Megacomplexes via Post-translational Modifications under Light Disturbance

Photosynthetic microorganisms rely on multiple pathways in central carbon metabolism to adapt to fluctuating light and energy availability across diel cycles. Mechanistic insight into the regulatory dynamics of this adaptation requires integrating processes spanning disparate timescales, from rapid redox-dependent post-translational modifications (PTMs) to slower changes in protein expression and metabolic pathway usage. To address this complexity beyond genome-based inference and traditional modeling, we develop a whole-cell four-dimensional (3D + time) model of the marine cyanobacterium Prochlorococcus marinus MED4 that explicitly represents the spatial organization of enzymatic and molecular processes in central carbon metabolism under light perturbation. We employ a perturbation-based research design to experimentally generate time-series, multi-omics measurements that provide molecular descriptors and cryo-ET derived 3D segmented volumes as constraints for this dynamic 4D framework. The integration of experiments and modeling across defined light regimes enables quantitative validation of system-level responses and forecasting under distinct light disturbances. We test the hypothesis that light-dependent redox PTMs regulating the structural assembly of a protein megacomplex, the “dark complex,” modulate metabolic flux at a conserved regulatory node of the Calvin–Benson cycle (CBC) in cyanobacteria. Our model shows that subcellular spatial organization buffers rapid light-induced changes in thylakoid reaction rates, which are followed by redox-PTM-mediated sequestration or release of CBC enzymes in the dark complex, ultimately impacting carbon fixation dynamics within carboxysomes. Comparison with an equivalently parameterized well-mixed stochastic model demonstrates that post-translational regulation not only buffers transcriptional noise and diffusion-driven fluctuations but also stabilizes phenotypic outcomes, underscoring the importance of spatial heterogeneity in phenotypic robustness. This ability to probe adaptive, spatiotemporally resolved mechanisms in photosynthetic machinery and central carbon metabolism addresses a critical gap in genotype-to-phenotype inference and expands modeling and design capabilities for understudied or genetically intractable autotrophs such as P. marinus MED4.

Johnson, Connah G.↗
DOE OSTI2020

Experimental Study on the Effect of Localized Blockages on the Friction Factor of a 61-Pin Wire-Wrapped Bundle

Abstract The thermal-hydraulic behavior of the flow in rod bundles has motivated numerous experimental and computational investigations. Previous studies have identified potential for accumulation of debris within the small subchannels of typical wire-wrapped assemblies with subsequent total or partial blockage of subchannel coolant flow. A test campaign was conducted to study the effects of localized blockages on the bundle averaged friction factor of a tightly packed wire-wrapped rod bundle. Blockages were installed within the bundle, and fluid pressure drop was measured across one wire pitch for a Reynolds number range of 500–17,200. The Darcy–Weisbach friction factor of the perturbed rod bundle geometry was compared with that of the unblocked bundle, as well as with the predictions of a well-established friction factor correlation. Differing effects based on blockage size and location for various flow regimes were studied. A number of conclusions can be made about the effects of the blockages on the friction factor, such as an increasing effect of the blockage on friction factor with an increase in Reynolds number, a change in flow behavior in the turbulent transition flow regime near Reynolds number 3000, differences in effect on friction factor for different types of subchannel blockage, and a nonlinear trend in friction factor variation with flow area impeded for edge subchannels. To this end, all data and quantified uncertainty produced in this study are made available for comparison and validation of advanced computational tools.

Engineering↗