Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “perturbation methods”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Random Phase Approximation Correlation Energy Using Real-Space Density Functional Perturbation Theory

We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn–Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic-scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method’s accuracy by comparing with plane-wave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 s on 4608 processors.

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

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

DeepAdversaries: examining the robustness of deep learning models for galaxy morphology classification

With increased adoption of supervised deep learning methods for work with cosmological survey data, the assessment of data perturbation effects (that can naturally occur in the data processing and analysis pipelines) and the development of methods that increase model robustness are increasingly important. In the context of morphological classification of galaxies, we study the effects of perturbations in imaging data. In particular, we examine the consequences of using neural networks when training on baseline data and testing on perturbed data. We consider perturbations associated with two primary sources: (a) increased observational noise as represented by higher levels of Poisson noise and (b) data processing noise incurred by steps such as image compression or telescope errors as represented by one-pixel adversarial attacks. We also test the efficacy of domain adaptation techniques in mitigating the perturbation-driven errors. We use classification accuracy, latent space visualizations, and latent space distance to assess model robustness in the face of these perturbations. For deep learning models without domain adaptation, we find that processing pixel-level errors easily flip the classification into an incorrect class and that higher observational noise makes the model trained on low-noise data unable to classify galaxy morphologies. On the other hand, we show that training with domain adaptation improves model robustness and mitigates the effects of these perturbations, improving the classification accuracy up to 23% on data with higher observational noise. Domain adaptation also increases up to a factor of ${\approx}2.3$ the latent space distance between the baseline and the incorrectly classified one-pixel perturbed image, making the model more robust to inadvertent perturbations. Successful development and implementation of methods that increase model robustness in astronomical survey pipelines will help pave the way for many more uses of deep learning for astronomy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Methods and compositions for obtaining useful plant traits

The present invention provides methods for obtaining plants that exhibit useful traits by perturbation of plastid function in plant rootstocks and grafting the rootstocks to scions. Methods for identifying genetic loci that provide for useful traits in plants and plants produced with those loci are also provided. In addition, plants that exhibit the useful traits, parts of the plants including seeds, and products of the plants are provided as well as methods of using the plants. Recombinant DNA vectors and transgenic plants comprising those vectors that provide for plastid perturbation are also provided.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Methods and compositions for obtaining useful plant traits

The present invention provides methods for obtaining plants that exhibit useful traits by perturbation of plastid function in plant rootstocks and grafting the rootstocks to scions. Methods for identifying genetic loci that provide for useful traits in plants and plants produced with those loci are also provided. In addition, plants that exhibit the useful traits, parts of the plants including seeds, and products of the plants are provided as well as methods of using the plants. Recombinant DNA vectors and transgenic plants comprising those vectors that provide for plastid perturbation are also provided.

Mackenzie, Sally↗
DOE OSTI2020

Selected Negative Linear Compressibilities in the Metal–Organic Framework of [Cu(4,4'-bpy) 2 (H 2 O) 2 ]·SiF 6

Negative linear compressibility (NLC) is a rare high-pressure observation that lattice contraction is accompanied by the structural expansion along a specific direction. Generally, this counterintuitive phenomenon only derives from the intrinsic structural characteristics of materials and cannot be tuned by external perturbations. Searching for an effective method to control NLC effect is still a challenge in both chemical and material science. In this work, we successfully discover and select the NLC behaviors in the metal-organic framework (MOF) of [Cu(4,4'-bpy) 2 (H 2 O) 2 ]·SiF 6 (Cu(bpy)·SiF) through controlling the pressure conditions therein. The NLC effect of Cu(bpy)·SiF originates from the wine-rack mechanism that quasi-square grids transfer to rhombic ones with the structural expansion along the diagonal direction at high pressure. Meanwhile, both of the pressure range and magnitude of the NLC responses are enlarged with optimized pressure conditions. This study not only presents the intriguing selected NLC behaviors of a MOF but also proves the effects of pressure conditions on NLC, which offers promising strategies for further design and applications of NLC materials.

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

Zero-Field Splitting Calculations by Multiconfiguration Pair-Density Functional Theory

Zero-field splitting (ZFS) is a fundamental molecular property that is especially relevant for single-molecule magnets (SMMs), electron paramagnetic resonance spectra, and quantum computing. Developing a method that can accurately predict ZFS parameters can be very powerful for designing new SMMs. One of the challenges is to include external correlation in an inherently multiconfigurational open-shell species for the accurate prediction of magnetic properties. Previously available methods depend on expensive multireference perturbation theory calculations to include external correlation. In this paper, we present spin-orbit-inclusive multiconfiguration and multistate pair-density functional theory (MC-PDFT) calculations of ZFSs; these calculations have a cost comparable to complete-active-space self-consistent field (CASSCF) theory, but they include correlation external to the active space. We found that combining a multistate formulation of MC-PDFT, namely, compressed-state multistate pair-density functional theory, with orbitals optimized by weighted-state-averaged CASSCF, yields reasonably accurate ZFS results

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

GPU Acceleration of Large-Scale Full-Frequency GW Calculations

Many-body perturbation theory is a powerful method to simulate electronic excitations in molecules and materials starting from the output of density functional theory calculations. By implementing the theory efficiently so as to run at scale on the latest leadership high-performance computing systems it is possible to extend the scope of GW calculations. Here, we present a GPU acceleration study of the full-frequency GW method as implemented in the WEST code. Excellent performance is achieved through the use of (i) optimized GPU libraries, e.g., cuFFT and cuBLAS, (ii) a hierarchical parallelization strategy that minimizes CPU-CPU, CPU-GPU, and GPU-GPU data transfer operations, (iii) nonblocking MPI communications that overlap with GPU computations, and (iv) mixed precision in selected portions of the code. A series of performance benchmarks has been carried out on leadership high-performance computing systems, showing a substantial speedup of the GPU-accelerated version of WEST with respect to its CPU version. Good strong and weak scaling is demonstrated using up to 25 920 GPUs. Finally, we showcase the capability of the GPU version of WEST for large-scale, full-frequency GW calculations of realistic systems, e.g., a nanostructure, an interface, and a defect, comprising up to 10 368 valence electrons.

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

Electron and ion spectroscopy of azobenzene in the valence and core shells

Azobenzene is a prototype and a building block of a class of molecules of extreme technological interest as molecular photo-switches. We present a joint experimental and theoretical study of its response to irradiation with light across the UV to x-ray spectrum. The study of valence and inner shell photo-ionization and excitation processes combined with measurement of valence photoelectron-photoion coincidence and mass spectra across the core thresholds provides a detailed insight into the site- and state-selected photo-induced processes. Photo-ionization and excitation measurements are interpreted via the multi-configurational restricted active space self-consistent field method corrected by second order perturbation theory. Using static modeling, we demonstrate that the carbon and nitrogen K edges of azobenzene are suitable candidates for exploring its photoinduced dynamics thanks to the transient signals appearing in background-free regions of the NEXAFS and XPS.

Chemistry↗
DOE OSTI2026

Understanding polaronic transport in complex oxides by combining precise synthesis and first-principles many-body theory

In complex oxides, charge carriers often couple strongly with lattice vibrations to form polarons–entangled electron–phonon quasiparticles whose transport properties remain difficult to characterize. Experimental access to intrinsic polaronic transport requires ultraclean samples, while theoretical description demands methods beyond low-order perturbation theory. Here, we show a predictive theory–experiment workflow to study polaron transport in complex oxides. Focusing on a prototypical polaronic oxide, anatase TiO 2 , we combine growth of high-quality oxygen-vacancy-doped films using hybrid molecular beam epitaxy with a first-principles electron–phonon diagrammatic Monte-Carlo (FEP-DMC) framework recently developed for accurate polaron predictions. Our films exhibit record-high electron mobility for anatase TiO 2 , in excellent agreement with FEP-DMC calculations conducted prior to experiment, which predict a room-temperature mobility of 45 ± 15 cm −2 V −1 s −1 and a mobility-temperature scaling of μ ∝ T −1.9 ± 0.077 . Microscopic analysis using scanning transmission electron microscopy and x-ray photoelectron spectroscopy reveals the role of oxygen vacancies in modulating transport at lower temperatures. FEP-DMC further provides quantitative insight into polaron formation energy, phonon cloud distribution, lattice distortion around the polaron, and the polaronic contribution to mobility. Together, these results provide a deeper microscopic understanding of large-polaron transport in a complex oxide and provide the blueprint to characterize other polaronic materials.

anatase TiO2↗
DOE OSTI2021

Density functional study of the variants of inter-Coulombic decay resonances in the photoionization of Cl@C 60

Inter-Coulombic decay (ICD) resonances in the photoionization of Cl@C 60 endofullerene molecule are calculated using a perturbative density functional theory (DFT) method. This is the first ICD study of an open shell atom in a fullerene cage. Here, three classes of resonances are probed: (i) Cl inner vacancies decaying through C 60 outer continua, (ii) C 60 inner vacancies decaying through Cl outer continua, and (iii) inner vacancies of either system decaying through the continua of Cl-C 60 hybrid levels, the hybrid Auger-ICD resonances. Comparisons with Ar@C 60 results reveal that the properties of hybrid Auger-ICD resonances are affected by the extent of level hybridization.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Lieb-Robinson Bounds with Exponential-in-Volume Tails

Lieb-Robinson bounds demonstrate the emergence of locality in many-body quantum systems. Intuitively, Lieb-Robinson bounds state that, with local or exponentially decaying interactions, the correlation that can be built up between two sites separated by distance 𝑟 after a time 𝑡 decays as exp (𝑣⁢𝑡 −𝑟), where 𝑣 is the emergent Lieb-Robinson velocity. In many problems, it is important to also capture how much of an operator grows to act on 𝑟 𝑑 sites in 𝑑 spatial dimensions. Perturbation theory and cluster expansion methods suggest that, at short times, these volume-filling operators are suppressed as exp (−𝑟 𝑑 ). We confirm this intuition, showing that, for 𝑟 >𝑣⁢𝑡, the volume-filling operator is suppressed by exp (−(𝑟−𝑣⁢𝑡) 𝑑 /(𝑣⁢𝑡) 𝑑−1 ). This closes a conceptual and practical gap between the cluster expansion and the Lieb-Robinson bound. We then present two very different applications of this new bound. Firstly, we obtain improved bounds on the classical computational resources necessary to simulate many-body dynamics with error tolerance 𝜀 for any finite time 𝑡: as 𝜀 becomes sufficiently small, only 𝜀 −O⁡(𝑡 𝑑−1 ) resources are needed. A protocol that likely saturates this bound is given. Secondly, we prove that disorder operators have volume-law suppression near the “solvable (Ising) point” in quantum phases with spontaneous symmetry breaking, which implies a new diagnostic for distinguishing many-body phases of quantum matter.

computational complexity↗
DOE OSTI2020

Nonrenormalization and Operator Mixing via On-Shell Methods

Using on-shell methods, we present a new perturbative nonrenormalization theorem for operator mixing in massless four-dimensional quantum field theories. By examining how unitarity cuts of form factors encode anomalous dimensions, we show that longer operators are often restricted from renormalizing shorter operators at the first order where Feynman diagrams exist. The theorem applies quite generally and depends only on the field content of the operators involved. We apply our theorem to operators of dimension five through seven in the standard model effective field theory, including examples of nontrivial zeros in the anomalous-dimension matrix at one through four loops. The zeros at two and higher loops go beyond those previously explained using helicity selection rules. We also include explicit sample calculations at two loops.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Disorder, interactions, and their interplay in novel narrow-gap Dirac materials and Weyl semimetals

Progress of the modern day condensed matter physics is to a large extent driven by the synthesis of new materials, advances in their experimental characterization and theoretical description. Recent discoveries of novel gapless Weyl semimetals, such as NaBi,CdAs, and BiTe-based films, in which magnetic dopants essentially suppress the gap, have added to the family of graphene and topological insulators actively investigated over the past decade. With the field of novel semimetals rapidly maturing, its focus necessarily shifts from demonstrations of the feasibility of such materials to their quantitative characterization. While the transport and optical properties of graphene and topological insulators are well captured within the picture of free non-interacting electrons, gapless 3D Weyl semimetals and narrow-gap 2D semiconductors with Dirac spectrum are known to be extremely susceptible to disorder and electron-electron interactions. This susceptibility obscures the manifestations of nontrivial band structure -- like quantum anomalous Hall effect -- of the new topological materials. Among particular projects to be addressed are: 1) optical conductivity of 3D gapless Dirac fermions in the presence of smooth disorder, 2) interplay of disorder and Coulomb interactions in the spectral properties of such fermions, 3) formation and structure of the impurity band with Coulomb supercritical clusters, 4) Coulomb interaction-driven renormalization of the electron spectrum and of the transport response in the presence of strong magnetic field, 5) instanton approach to the disorder-induced fluctuation states in zero-gap 3D materials, and 6) the role of disorder in quantum anomalous Hall effect. The proposal relies upon the investigators' previous broad expertise in interacting and disordered electron systems. The methods to be employed include perturbative diagrammatic technique, non-perturbative instanton and self-consistent approximations, hydrodynamics of electron liquid. Both analytical as well as numerical approaches are to be employed. The anticipated broader outcome of the proposal includes gaining an in-depth understanding of the interplay of the disorder and interactions under the conditions when this interplay has the most dramatic impact on observables. Traditionally, interaction effects are among the most challenging and interesting problems of condensed matter physics. Similarly, disordered systems typically present very difficult but extremely rich problems in the description of various materials. Importantly, understanding the spectral and transport properties of such materials not only presents the fundamental objective, but is also of particular interest for many applications, such as computation, memory, optics, plasmonics. In particular realization of the quantum anomalous Hall effect may lead to the development of low-power-consumption electronics. Indeed, a major constraint for practical use of the quantum Hall effect is limited by the requirement of the quantizing magnetic field. At the same time, the quantum anomalous Hall effect samples exhibit non-dissipative edge quantum transport in a zero magnetic field.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Scattering wave packets of hadrons in gauge theories: Preparation on a quantum computer

Quantum simulation holds promise of enabling a complete description of high-energy scattering processes rooted in gauge theories of the Standard Model. A first step in such simulations is preparation of interacting hadronic wave packets. To create the wave packets, one typically resorts to adiabatic evolution to bridge between wave packets in the free theory and those in the interacting theory, rendering the simulation resource intensive. In this work, we construct a wave-packet creation operator directly in the interacting theory to circumvent adiabatic evolution, taking advantage of resource-efficient schemes for ground-state preparation, such as variational quantum eigensolvers. By means of an ansatz for bound mesonic excitations in confining gauge theories, which is subsequently optimized using classical or quantum methods, we show that interacting mesonic wave packets can be created efficiently and accurately using digital quantum algorithms that we develop. Specifically, we obtain high-fidelity mesonic wave packets in the Z 2 and U(1) lattice gauge theories coupled to fermionic matter in 1+1 dimensions. Our method is applicable to both perturbative and non-perturbative regimes of couplings. The wave-packet creation circuit for the case of the Z 2 lattice gauge theory is built and implemented on the Quantinuum H1-1 trapped-ion quantum computer using 13 qubits and up to 308 entangling gates. The fidelities agree well with classical benchmark calculations after employing a simple symmetry-based noise-mitigation technique. This work serves as a step toward quantum computing scattering processes in quantum chromodynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Sensitivity Coefficients Calculated for the Prompt Neutron Decay Constant at or Near Delayed Critical [Abstract]

Experimenters at Los Alamos National Laboratory (LANL) measure the prompt neutron decay constant for many experiments at the National Criticality Experiments Research Center (NCERC) to infer reactivity and the effective neutron multiplication factor. These quantities are very important for nuclear criticality safety and validating nuclear data. Uncertainty in measures of criticality of an experimental configuration can be determined prior to physically performing the experiment by applying first order perturbation theory to Monte Carlo codes, such as MCNP®. The first order perturbation theory produces first derivatives of some nuclear parameter to nuclear data (e.g., cross section data). This first derivative is commonly referred to as a sensitivity coefficient. Currently, the MCNP® software has the capability of computing effective neutron multiplication factor sensitivity coefficients to cross section data. This work builds off of this MCNP® capability and the first order perturbation theory to provide a method of calculating sensitivity coefficients for the prompt neutron decay constant at or near delayed critical to cross section data. The prompt neutron decay constant sensitivity coefficient calculated in this work does not depend on any modification of the MCNP® source code. Prompt neutron decay constant sensitivity coefficient calculations can be used to infer reactivity and effective neutron multiplication factor sensitivity coefficient values as well. By investigating the trends of prompt neutron decay constant sensitivity coefficients for nuclide-reaction pairs across energy spectra, experiments can be designed to maximize or minimize the uncertainty in the prompt neutron decay constant in a particular energy region, which can lead to further optimization studies. Prompt neutron decay constant sensitivity coefficients will be calculated for the Jezebel benchmark. Subsequently, these sensitivity coefficients will be used in a data assimilation process to determine if there is or are optimal experiments that can be performed to provide insight into adjustments of uncertain/inaccurate cross section data. Specifically, the effect of the prompt neutron decay constant sensitivity coefficients on the nuclear data-induced uncertainty in the effective neutron multiplication factor will be examined

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗