Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Perturbative 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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Electric polarization in inhomogeneous crystals

In this work, we derive the charge density up to second order in spatial gradient in inhomogeneous crystals using the semiclassical coarse graining procedure based on the wave packet method. It can be recast as divergence of polarization, whose first-order contribution consists of three parts, a perturbative correction to the original Berry connection expression, a topological part that can be written as an integral of the Chern-Simons 3-form, and a previously-unknown, quadrupolelike contribution. The topological part can be related to the quantized fractional charge carried by a vortex in two-dimensional systems. We then generalize our results to the multiband case and show that the quadrupolelike contribution plays an important role, as it makes the total polarization gauge-independent. Finally, we verify our theory in several model systems.

36 MATERIALS SCIENCE↗
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 OSTI2021

The origin of CEMA and its relation to CSP

Currently, there exist two methods for analysing an explosive mode introduced by chemical kinetics in a reacting process: the Computational Singular Perturbation (CSP) algorithm and the Chemical Explosive Mode Analysis (CEMA). CSP was introduced in 1989 and addressed both dissipative and explosive modes encountered in the multi-scale dynamics that characterize the process, while CEMA was introduced in 2009 and addressed only the explosive modes. It is shown that (i) the algorithmic tools incorporated in CEMA were developed previously on the basis of CSP and (ii) the examination of explosive modes has been the subject of CSP-based works, reported before the introduction of CEMA.

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

Structure−Property Relationships in a Series of Hybrid Halopyridinium Ruthenium(IV) Halides

Reported is the synthesis and characterization of a family of hybrid ruthenium halides consisting of haloruthenium octahedra (X = Cl, Br) charge balanced with halopyridinium (XPy; X= H, Cl, Br, I) organic cations which assemble via noncovalent interactions between ion pairs. Diffuse reflectance spectroscopy showed that compounds containing RuBr 6 2- octahedra displayed a lower bandgap (1.05 eV< x < 1.08 eV) compared to compounds with RuCl 6 2- octahedra (1.22 eV < x < 1.43 eV). Additionally, computational density functional theory (DFT) based natural bonding orbital (NBO) analysis and density of state (DOS) methods were used to characterize second-sphere non-covalent interaction strengths and elucidate molecular orbital perturbations to elucidate their influence on Ru-X orbital constructs and in turn rationalize band gap trends. Analyses showed ligand to metal charge transfer is responsible for the small bandgap energies and are unperturbed by second-sphere interactions. Here, this report serves as a platform for probing the relationship between the structural and photophysical properties within the haloruthenium family of low dimensional transition metal halide perovskites.

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

Microstructural impacts on ionic conductivity of oxide solid electrolytes from a combined atomistic-mesoscale approach

Abstract Although multiple oxide-based solid electrolyte materials with intrinsically high ionic conductivities have emerged, practical processing and synthesis routes introduce grain boundaries and other interfaces that can perturb primary conduction channels. To directly probe these effects, we demonstrate an efficient and general mesoscopic computational method capable of predicting effective ionic conductivity through a complex polycrystalline oxide-based solid electrolyte microstructure without relying on simplified equivalent circuit description. We parameterize the framework for Li 7- x La 3 Zr 2 O 12 (LLZO) garnet solid electrolyte by combining synthetic microstructures from phase-field simulations with diffusivities from molecular dynamics simulations of ordered and disordered systems. Systematically designed simulations reveal an interdependence between atomistic and mesoscopic microstructural impacts on the effective ionic conductivity of polycrystalline LLZO, quantified by newly defined metrics that characterize the complex ionic transport mechanism. Our results provide fundamental understanding of the physical origins of the reported variability in ionic conductivities based on an extensive analysis of literature data, while simultaneously outlining practical design guidance for achieving desired ionic transport properties based on conditions for which sensitivity to microstructural features is highest. Additional implications of our results are discussed, including a possible connection between ion conduction behavior and dendrite formation.

25 ENERGY STORAGE↗
DOE OSTI2024

Out of time order correlation of the Hubbard model with random local disorder

The out-of-time-order correlator (OTOC) serves as a powerful tool for investigating quantum information spreading and chaos in complex systems. We present a method employing non-equilibrium dynamical mean-field theory and coherent potential approximation combined with diagrammatic perturbation on the Schwinger–Keldysh contour to calculate the OTOC for correlated fermionic systems subjected to both random disorder and electron interaction. Furthermore, our key finding is that random disorder enhances the OTOC decay in the Hubbard model for the metallic phase in the weakly interacting limit. However, the current limitation of our perturbative solver restricts the applicability to weak interaction regimes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Investigating the impact of intermediate-mode perturbations on diagnosing plasma conditions in DT cryogenic implosions via synthetic x-ray Thomson scattering

The pursuit of inertial confinement fusion ignition target designs requires precise experimental validation of the conditions within imploding capsules, in particular the density and temperature of the compressed shell. Previous work has identified x-ray Thomson scattering (XRTS) as a viable diagnostic tool for inferring the in-flight compressed deuterium-tritium shell conditions during capsule implosions (Poole et al 2022 Phys. Plasmas 29 072703). However, this study focused on one-dimensional simulations, which do not account for the growth of hydrodynamic instabilities. In this work, two-dimensional DRACO simulations incorporating intermediate-mode perturbations up to Legendre mode l = 50 were used to generate synthetic XRTS spectra with the SPECT3D code. The analysis employed Markov-Chain Monte Carlo techniques to infer plasma conditions from these spectra. The results demonstrate that the XRTS diagnostic platform can effectively discern the in-flight compressed shell conditions for targets with varying adiabats, even in the presence of intermediate-mode perturbations. This work underscores the potential of XRTS for realistic inertial confinement fusion experiments, providing a robust method for probing the complex dynamics of fusion implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

First-order formalism for β functions in bosonic sigma models from supersymmetry breaking

We consider the renormalization group flow equation for the two-dimensional sigma models with the Kähler target space. The first-order formulation allows us to treat perturbations in these models as current-current deformations. We demonstrate, however, that the conventional first-order formalism misses certain anomalies in the measure, and should be amended. We reconcile beta functions obtained within the conformal perturbation theory for the current-current deformations with traditional “geometric” results obtained in the background field methods, in this way resolving the peculiarities pointed out in O. Gamayun et al. [Peculiarities of beta functions in sigma models, J. High Energy Phys. 10 (2023) 097]. The result is achieved by the supersymmetric completion of the first-order sigma model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Optimally tuned starting point for single-shot GW calculations of solids

The dependence of ab initio many-body perturbation theory within the $\textit{GW}$ approximation on the eigensystem used in calculating quasiparticle corrections limits this method's predictive power. Here, we investigate the accuracy of the recently developed Wannier-localized optimally tuned screened range-separated hybrid (WOT-SRSH) functional as a generalized Kohn-Sham starting point for single-shot $\textit{GW}$ ($G_0W_0$) calculations for a range of semiconductors and insulators. Comparison to calculations based on well-established functionals, namely, PBE, PBE0, and HSE, as well as to self-consistent GW schemes and to experiment, shows that band gaps computed via $G_0W_0$@WOT-SRSH have a level of precision and accuracy that is comparable to that of more advanced methods such as quasiparticle self-consistent $\textit{GW}$ and eigenvalue self-consistent $\textit{GW}$. In this work, we also find that $G_0W_0$@WOT-SRSH improves the description of states deeper in the valence band manifold. Finally, we show that $G_0W_0$@WOT-SRSH significantly reduces the sensitivity of computed band gaps to ambiguities in the underlying WOT-SRSH tuning procedure.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Photodynamic melting of phase-reversed charge stripes and enhanced condensation

The interplay between charge stripes and pairing has long been a subject of scrutiny in a broad class of unconventional superconductors, as in some cases it is unclear whether this interplay benefits the ensuing superfluidity. Experiments that explore the out-of-equilibrium dynamics of these systems aim to tip the balance toward one phase or the other by selectively coupling to relevant modes. Leveraging the fact that competition between stripes and pairing is not exclusive to fermionic systems, we explore the photoirradiation dynamics of interacting hardcore bosons in which density-wave phase-reversal melting leads to enhanced phase-coherent transport response, as quantified by the dynamic amplification of both the zero-momentum occupancy and the condensate fraction, as well as finite out-of-equilibrium charge stiffness and superfluid weight, for a given system size. Finally, our results, obtained using unbiased methods for an interacting system on a ladder geometry, demonstrate how one can engineer time-dependent perturbations to release suppressed orders, potentially providing insight into the underlying mechanism in related experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Radial-offset Optics in EIC Hadron Lattices: Trajectory Lengthening

This Tech. Note describes a method for a radial orbit offset in the arcs of the EIC hadron lattices, based on the perturbation dB/B of the magnetic field of the arc main bend magnets. This can be done independently for each arc (allowing multiple rigidities) and requires the use of pairs of orbit correctors at the extremities of the arc to ensure that the radially shifted orbit is on the zero-orbit (i.e. with the horizontal coordinates x = 0, p x = 0) at the exit and entry point of the respectively upstream and downstream interaction region (IR). The goal is to produce appropriate orbit lengthening as required for time of flight adjustments between the electron and hadron beams as part of EIC operations. An ad hoc parameter is introduced: a “geometric compaction factor” α B = (dC/C)/(dB/B), to quantify the orbit lengthening, or equivalently the radial orbit shift dR = R · α B · dB/B, from bending perturbation. Detailed numerical results are produced and discussed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Impact of Coating Defects on Performance of Coated Zirconium Cladding

Research on accident tolerant fuels (ATF) has started after the Fukushima accident [1–3]. While efforts have been expended on both fuel and cladding ATF concepts, the bulk of work has been devoted to improved cladding. The overarching goal of these approaches is to extend the coping time available during a severe accident before any event would result in release of radioactivity to the public. The most basic ATF cladding concept is obtained by applying a thin coating of highly corrosion-resistant material on the surface of a licensed zirconium cladding alloy. This thin coating is intended to not interfere with the neutronic or mechanical performance of the base cladding under normal operating conditions. Different coating materials, thicknesses, coating processes, process parameters, and testing methods have an impact on the microstructure and mechanical properties and therefore on the results of the applied investigation methods. These challenges have motivated an initial focus on demonstrating that the presence of coatings do not perturb the critical performance benchmarks of uncoated material. Ongoing lead test assembly irradiations of coated zirconium concepts in commercial reactions is intended to establish baseline performance in this regard in the coming years.

36 MATERIALS SCIENCE↗
DOE OSTI2022

GANISP: A GAN-Assisted Importance Splitting Probability Estimator: Preprint

To reduce the variance of rare event probability estimator, genealogical importance splitting marches towards a rare event by iteratively selecting and replicating realizations that are headed towards a rare event. The replication step is made difficult when applied to deterministic systems as the initial conditions of the offspring realizations need to be modified. Typically, a random perturbation is applied to the offspring to differentiate their trajectory from the parent realization. It is shown that a random perturbation strategy may be effective for some systems but may also fail for others, thereby preventing variance reduction in the probability estimate. To address this limitation, it is proposed to use a generative model such as a Generative Adversarial Network (GAN) to generate perturbations that are consistent with the attractor of the dynamical system. The GAN-assisted Importance SPlitting method (GANISP) improves the variance reduction for the system targeted. An implementation of the method is available in a companion repository (https://github.com/NREL/GANISP).

generative adversarial network↗
DOE OSTI2022

Error analysis of numerical methods for thick diffusive neutron transport problems on Shishkin mesh

A thin layer will develop at the boundary if the incoming angular flux is anisotropic in thick diffusive neutron transport problems. Solving such singularly perturbed problems, which have non-smooth solutions with singularity near the boundary, is computationally challenging. Standard finite difference schemes on a uniform mesh cannot yield ε-uniform convergence, where ε is a small parameter, while it can be achieved on a suitable piecewise-uniform Shishkin mesh. We present a formal error analysis of the diamond difference (DD) method and step difference (SD) method for solving the S{sub N} neutron transport equation. The analysis can be extended to other finite difference methods. Numerical results are presented to confirm the error estimates and the advantages of the Shishkin mesh. (author)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

On High-Order/Low-Order and Micro-Macro Methods for Implicit Time-Stepping of the BGK Model

In this paper, a high-order/low-order (HOLO) method is combined with a micro-macro (MM) decomposition to accelerate iterative solvers in fully implicit time-stepping of the Bhatnagar–Gross–Krook (BGK) equation for gas dynamics. The MM formulation represents a kinetic distribution as the sum of a local Maxwellian and a perturbation. In highly collisional regimes, the perturbation away from initial and boundary layers is small and can be compressed to reduce the overall storage cost of the distribution. The convergence behavior of the MM methods, the usual HOLO method, and the standard source iteration method is analyzed on a linear BGK model. Both the HOLO and MM methods are implemented using a discontinuous Galerkin (DG) discretization in phase space, which naturally preserves the consistency between high- and low-order models required by the HOLO approach. Furthermore, the accuracy and performance of these methods are compared on the Sod shock tube problem and a sudden wall heating boundary layer problem. Overall, the results demonstrate the robustness of the MM and HOLO approaches and illustrate the compression benefits enabled by the MM formulation when the kinetic distribution is near equilibrium.

BGK model↗
DOE OSTI2020

Benchmarking an Embedded Adaptive Sampling Configuration Interaction Method for Surface Reactions: H 2 Desorption from and CH 4 Dissociation on Cu(111)

Embedded (emb-) correlated wavefunction (CW) theory enables accurate assessments of both ground- and excited-state reaction mechanisms involved in heterogeneous catalysis. Embedded multireference second-order perturbation theory (emb-MRPT2) based on reference wavefunctions generated via embedded complete active space self-consistent field (emb-CASSCF) theory is currently state-of-the-art. However, the factorial scaling of CASSCF limits the size of active space and the complexity of systems that can be studied. In this work, we assess the efficacy of an alternative CW method, adaptive sampling configuration interaction (ASCI)–which enables large active spaces to be used–for studying surface reactions. We couple ASCI with density functional embedding theory (DFET) and benchmark its performance for two reactions: H 2 desorption from and CH 4 dissociation on the Cu(111) surface. Unlike embedded complete active space second-order perturbation theory (emb-CASPT2) that accurately reproduces a measured H 2 desorption barrier, embedded ASCI, using a very large active space (though one that still comprises a small portion of the full set of orbitals) fails to do so. Adding an extra correlation term from embedded Møller–Plesset second-order perturbation theory (emb-MP2) improves the desorption barrier and endothermicity predictions. Thus, the inaccuracy of embedded ASCI comes from the missing dynamic correlation from the many other electrons and orbitals not included in the active space. For CH 4 dissociation, again embedded ASCI overestimates the dissociation barrier compared to emb-CASPT2 predictions. Adding dynamic correlation from emb-MP2 helps correct the barrier. However, this composite approach suffers from double counting of correlation within embedded ASCI followed by emb-MP2 calculations. We therefore conclude that the state-of-the-art emb-MRPT2 based on reference wavefunctions generated via emb-CASSCF remains the method of choice for studying surface reactions. emb-ASCI is useful when large active spaces beyond the limit of emb-CASSCF are essential, such as to study complex surface reactions with significant multiconfigurational character (static correlation) but weak dynamic correlation.

08 HYDROGEN↗
DOE OSTI2021

Quantum Solvers for Plane-Wave Hamiltonians: Abridging Virtual Spaces Through the Optimization of Pairwise Correlations

For many-body methods such as MCSCF and CASSCF, in which the number of one-electron orbitals is optimized and independent of the basis set used, there are no problems with using plane-wave basis sets. However, for methods currently used in quantum computing such as select configuration interaction (CI) and coupled cluster (CC) methods, it is necessary to have a virtual space that is able to capture a significant amount of electron-electron correlation in the system. The virtual orbitals in a pseudopotential plane-wave Hartree–Fock calculation, because of Coulomb repulsion, are often scattering states that interact very weakly with the filled orbitals. As a result, very little correlation energy is captured from them. The use of virtual spaces derived from the one-electron operators has also been tried, and while some correlations are captured, the amount is quite low. To overcome these limitations, we have been developing new classes of algorithms to define virtual spaces by optimizing orbitals from small pairwise CI Hamiltonians, which we term as correlation optimized virtual orbitals with the abbreviation COVOs. With these procedures, we have been able to derive virtual spaces, containing only a few orbitals, which are able to capture a significant amount of correlation. The focus in this manuscript is on using these derived basis sets to target full CI (FCI) quality results for H2 on near-term quantum computers. However, the initial results for this approach were promising. We were able to obtain good agreement with FCI/cc-pVTZ results for this system with just 4 virtual orbitals, using both FCI and quantum simulations. The quality of the results using COVOs suggests that it may be possible to use them in other many-body approaches, including coupled cluster and Møller–Plesset perturbation theories, and open up the door to many-body calculations for pseudopotential plane-wave basis set methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗