Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “self consistent field 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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Ab initio study of the molecular structure and vibrational spectrum of nitric acid and its protonated forms

The equilibrium structures, harmonic vibrational frequencies, IR intensities, and relative energetics of HNO3 and its protonated form H2NO3+ were investigated using double-zeta plus polarization and triple-zeta plus polarization basis sets in conjunction with high-level ab initio methods. The latter include second-order Moller-Plesset perturbation theory, the single and double excitation coupled cluster (CCSD) methods, a perturbational estimate of the effects of connected triple excitations (CCSD(T)), and the self-consistent field. To determine accurate energy differences CCSD(T) energies were computed using large atomic natural orbital basis sets. Four different isomers of H2NO3+ were considered. The lowest energy form of protonated nitric acid was found to correspond to a complex between H2O and NO2+, which is consistent with earlier theoretical and experimental studies.

Lee, Timothy J.↗
DOE OSTI2022

A non-hierarchical correlation discrete variable representation

The correlation discrete variable representation (CDVR) facilitates (multi-layer) multi-configurational time-dependent Hartree (MCTDH) calculations with general potentials. It employs a layered grid representation to efficiently evaluate all potential matrix elements appearing in the MCTDH equations of motion. The original CDVR approach and its multi-layer extension show a hierarchical structure: the size of the grids employed at the different layers increases when moving from an upper layer to a lower one. In this work, a non-hierarchical CDVR approach, which uses identically structured quadratures at all layers of the MCTDH wavefunction representation, is introduced. Here, the non-hierarchical CDVR approach crucially reduces the number of grid points required, compared to the hierarchical CDVR, shows superior scaling properties, and yields identical results for all three representations showing the same topology. Numerical tests studying the photodissociation of NOCl and the vibrational states of CH 3 demonstrate the accuracy of the non-hierarchical CDVR approach.

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

A new step forward in realistic cluster lens mass modelling: analysis of Hubble Frontier Field Cluster Abell S1063 from joint lensing, X-ray, and galaxy kinematics data

We present a new method to simultaneously and self-consistently model the mass distribution of galaxy clusters that combines constraints from strong lensing features, X-ray emission, and galaxy kinematics measurements. We are able to successfully decompose clusters into their collisionless and collisional mass components thanks to the X-ray surface brightness, as well as use the dynamics of cluster members, to obtain more accurate masses exploiting the fundamental plane of elliptical galaxies. Knowledge from all observables is included through a consistent Bayesian approach in the likelihood or in physically motivated priors. We apply this method to the galaxy cluster Abell S1063 and produce a mass model that we publicly release with this paper. The resulting mass distribution presents different ellipticities for the intra-cluster gas and the other large-scale mass components as well as deviation from elliptical symmetry in the main halo. We assess the ability of our method to recover the masses of the different elements of the cluster using a mock cluster based on a simplified version of our Abell S1063 model. Thanks to the wealth of mutliwavelength information provided by the mass model and the detected X-ray emission, we also found evidence for an ongoing merger event with gas sloshing from a smaller infalling structure into the main cluster. In agreement with previous findings, the total mass, gas profile, and gas mass fraction are all consistent with small deviations from the hydrostatic equilibrium. This new mass model for Abell S1063 is publicly available, as the lenstool extension used to construct it.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Graph-based quantum response theory and shadow Born–Oppenheimer molecular dynamics

Graph-based linear scaling electronic structure theory for quantum-mechanical molecular dynamics simulations [A. M. N. Niklasson et al., J. Chem. Phys. 144, 234101 (2016)] is adapted to the most recent shadow potential formulations of extended Lagrangian Born–Oppenheimer molecular dynamics, including fractional molecular-orbital occupation numbers [A. M. N. Niklasson, J. Chem. Phys. 152, 104103 (2020) and A. M. N. Niklasson, Eur. Phys. J. B 94, 164 (2021)], which enables stable simulations of sensitive complex chemical systems with unsteady charge solutions. The proposed formulation includes a preconditioned Krylov subspace approximation for the integration of the extended electronic degrees of freedom, which requires quantum response calculations for electronic states with fractional occupation numbers. For the response calculations, we introduce a graph-based canonical quantum perturbation theory that can be performed with the same natural parallelism and linear scaling complexity as the graph-based electronic structure calculations for the unperturbed ground state. Further, the proposed techniques are particularly well-suited for semi-empirical electronic structure theory, and the methods are demonstrated using self-consistent charge density-functional tight-binding theory both for the acceleration of self-consistent field calculations and for quantum-mechanical molecular dynamics simulations. Graph-based techniques combined with the semi-empirical theory enable stable simulations of large, complex chemical systems, including tens-of-thousands of atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2012

Plasma Diffusion in Self-Consistent Fluctuations

The problem of particle diffusion in position space, as a consequence ofeleclromagnetic fluctuations is addressed. Numerical results obtained with a self-consistent hybrid code are presented, and a method to calculate diffusion coefficient in the direction perpendicular to the mean magnetic field is proposed. The diffusion is estimated for two different types of fluctuations. The first type (resuiting from an agyrotropic in itiai setting)is stationary, wide band white noise, and associated to Gaussian probability distribution function for the magnetic fluctuations. The second type (result ing from a Kelvin-Helmholtz instability) is non-stationary, with a power-law spectrum, and a non-Gaussian probabi lity distribution function. The results of the study allow revisiting the question of loading particles of solar wind origin in the Earth magnetosphere.

Smets, R.↗
NASA NTRS1997

Self-Consistent Superthermal Electron Effects on Plasmaspheric Refilling

The effects of self-consistently including superthermal electrons in the definition of the ambipolar electric field are investigated for the case of plasmaspheric refilling after a geomagnetic storm. By using the total electron population in the hydrodynamic equations, a method for incorporating superthermal electron parameters in the electric field and electron temperature calculation is developed. Also, the ambipolar electric field is included in the kinetic equation for the superthermal electrons through a change of variables using the total energy and the first adiabatic invariant. Calculations based on these changes are performed by coupling time-dependent models of the thermal plasma and superthermal electrons. Results from this treatment of the electric field and the self-consistent development of the solution are discussed in detail. Specifically, there is a decreased thermal electron density in the plasmasphere during the first few minutes of refilling, a slightly accelerated proton shock front, and a decreased superthermal electron flux due to the deceleration by the electric field. The timescales of plasmaspheric refilling are discussed and determined to be somewhat shorter than previously calculated for the thermal plasma and superthermal electron population due to the effects of the field-aligned potential.

Liemohn, M. W.↗
DOE OSTI2021

Simulating dust grain-radiation coupling on a moving mesh

ABSTRACT We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code arepo-rt, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in arepo-rt, which treats UV, optical photons as single scattered and IR photons as multi scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Coupled magnetostructural continuum model for multiferroic BiFeO 3

Here, a continuum approach to study magnetoelectric multiferroic BiFeO 3 (BFO) is proposed. Our modeling effort marries the ferroelectric (FE) phase field method and micromagnetic simulations to describe the entire multiferroic order parameter sector (polarization, oxygen antiphase tilts, strain, and magnetism) self-consistently on the same time- and length scale. In this paper, we discuss our choice of FE and magnetic energy terms and demonstrate benchmarks against known behavior. We parametrize the lowest order couplings of the structural distortions against previous predictions from density functional theory calculations giving access to simulations of the FE domain wall (DW) topology. This allows us to estimate the energetic hierarchy and thicknesses of the numerous structural DWs. We then extend the model to the canted antiferromagnetic order and demonstrate how the FE domain boundaries influence the resulting magnetic DWs. We also highlight some capabilities of this model by providing two examples relevant for applications. We demonstrate spin-wave transmission through the multiferroic domain boundaries which identify rectification in qualitative agreement with recent experimental observations. As a second example of application, we model fully dynamical magnetoelectric switching, where we find a sensitivity on the Gilbert damping with respect to switching pathways. We envision that this modeling effort will set the basis for further work on properties of arbitrary 3D nanostructures of BFO (and related multiferroics) at the mesoscale.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1997

The Vibrational Frequencies of CaO2, ScO2, and TiO2: A Comparison of Theoretical Methods

The vibrational frequencies of several states of CaO2, ScO2, and TiO2 are computed at using density functional theory (DFT), the Hatree-Fock approach, second order Moller-Plesset perturbation theory (MP2), and the complete-active-space self-consistent-field theory. Three different functionals are used in the DFT calculations, including two hybrid functionals. The coupled cluster singles and doubles approach including the effect of unlinked triples, determined using perturbation theory, is applied to selected states. The Becke-Perdew 86 functional appears to be the cost effective method of choice, although even this functional does not perform well for one state of CaO2. The MP2 approach is significantly inferior to the DFT approaches.

Rosi, Marzio↗
DOE OSTI2021

Electronic structure of LaNiO 2 and CaCuO 2 from a self-consistent vertex-corrected GW approach

In this work, the electronic structure of one of the nickelates (LaNiO 2 ) and one of the cuprates (CaCuO 2 ) is studied with three self-consistent GW-based methods: scGW, sc(GW+vertex), and quasiparticle self-consistent GW. Low-energy features obtained in our study are in many respects similar to the features reported in previous density functional theory plus dynamical mean-field theory (DFT+DMFT) studies. Consistent with the DFT+DMFT conclusion, we find LaNiO 2 to be more correlated than CaCuO 2 . However, correlation effects included in our study change the DFT Fermi surface near the Γ point differently from that reported in DMFT studies. Features that are a few electronvolts away from the Fermi level are broader in our calculations than in the DFT+DMFT, which reflects the differences between the DFT and the GW methods. Our results are in qualitative agreement with previous G0W0 results, but the self-consistency results in quantitative differences. Generally, correlation effects are found to be sufficiently weak in both materials, which allows one to use totally ab initio diagrammatic approaches such as sc(GW+Vertex) and to avoid the methods with adjustable parameters (DFT+U or DFT+DMFT). However, the possibility of some strong correlations at low energy that cannot be captured by perturbative methods cannot be completely excluded. For instance, differences in the Fermi surface should be resolved, thus experimental studies are necessary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1986

Theoretical research program to study chemical reactions in AOTV bow shock tubes

Progress in the development of computational methods for the characterization of chemical reactions in aerobraking orbit transfer vehicle (AOTV) propulsive flows is reported. Two main areas of code development were undertaken: (1) the implementation of CASSCF (complete active space self-consistent field) and SCF (self-consistent field) analytical first derivatives on the CRAY X-MP; and (2) the installation of the complete set of electronic structure codes on the CRAY 2. In the area of application calculations the main effort was devoted to performing full configuration-interaction calculations and using these results to benchmark other methods. Preprints describing some of the systems studied are included.

Taylor, P.↗
DOE OSTI2023

Toward continuum gyrokinetic study of high-field mirrors

High-temperature superconducting (HTS) magnetic mirrors under development exploit strong fields with high mirror ratio to compress loss cones and enhance confinement and may offer cheaper, more compact fusion power plant candidates. This new class of devices could exhibit largely unexplored interchange and gradient-driven modes. Such instabilities, and methods to stabilize them, can be studied with gyrokinetics, given the strong magnetization and prevalence of kinetic effects. Our focus here is to (a) determine if oft-used gyrokinetic models for open field lines produce the electron-confining (Pastukhov) electrostatic potential and (b) examine and address challenges faced by gyrokinetic codes in studying HTS mirrors. Here, we show that a one-dimensional limit of said models self-consistently develops a potential qualitatively approaching the analytical Pastukhov level. Additionally, we describe the computational challenges of studying high mirror ratios with open field line gyrokinetic solvers and offer a force softening method to mitigate small time steps needed for time integration in colossal magnetic field gradients produced by HTS coils, providing a 19X speedup.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Hybrid algorithm for the time-dependent Hartree–Fock method using the Yang–Baxter equation on quantum computers *

Abstract The time-dependent Hartree–Fock (TDHF) method is an approach to simulate the mean field dynamics of electrons within the assumption that the electrons move independently in their self-consistent average field and within the space of single Slater determinants. One of the major advantages of performing time dynamics within Hartree–Fock theory is the free fermionic nature of the problem, which makes TDHF classically simulatable in polynomial time. Here, we present a hybrid TDHF implementation for quantum computers. This quantum circuit grows with time; but with our recent work on circuit compression via the Yang–Baxter equation (YBE), the resulting circuit is constant depth. This study provides a new way to simulate TDHF with the aid of a quantum device as well as provides a new direction for the application of YBE symmetry in quantum chemistry simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Nonunitary Variational Quantum Eigensolver with the Localized Active Space Method and Cost Mitigation

Accurately describing strongly correlated systems with affordable quantum resources remains a central challenge for quantum chemistry applications on near and intermediate term quantum computers. The localized active space self-consistent field (LASSCF) approximates the complete active space self-consistent field (CASSCF) by generating active space-based wave functions within specific fragments while treating interfragment correlation with mean-field approach, hence is computationally less expensive. Hardware-efficient ansatzes (HEA) offer affordable and shallower circuits, yet they often fail to capture the necessary correlation. Previously, Jastrow-factor-inspired nonunitary qubit operators were proposed to use with HEA for variational quantum eigensolver (VQE) calculations (so-called nuVQE), as they do not increase circuit depths and recover correlation beyond the mean-field level for Hartree–Fock initial states. Here, in this study, we explore running nuVQE with LASSCF as the initial state. The method, named LAS-nuVQE, is shown to recover interfragment correlations, reach chemical accuracy with a small number of gates (<70) in both H 4 and square cyclobutadiene (C 4 H 4 ), and produces more accurate energetics than its HEA counterparts at all circuit depths. To further address the inherent symmetry-breaking in HEA, we implemented spin-constrained LAS-nuVQE to extend the capabilities of HEA further and show spin-pure results for square cyclobutadiene. We also mitigate the increased measurement overhead of nuVQE via Pauli grouping and shot-frugal sampling, reducing measurement costs by up to 2 orders of magnitude compared to ungrouped operator, and show that one can achieve better accuracy with a small number of shots (10 3–4 ) per one expectation value calculation compared to noiseless simulations with one or two orders of magnitude more shots. Finally, wall clock time estimates show that, with our measurement mitigation protocols, nuVQE becomes a cheaper and more accurate alternative than vanilla VQE with HEA. Taken together, these developments illustrate a practical pathway toward performing multireference chemical simulations with accuracy and affordable resources on today’s quantum hardware, achieving both accuracy and affordability in challenging correlated systems.

Wang, Qiaohong [Univ. of Chicago, IL (United State↗
NASA NTRS1987

The role and behavior of spin in gravitational physics

A self-consistent method of introducing spin into any Lagrangian based theory of gravitation was developed. The metric variation of the Lagrangian in the theory leads to an improved energy-momentum tensor which represents the source term in the gravitational field equations. The goal of the research is the construction of a theory general enough to be used to investigate spin effects in astrophysical objects and cosmology, and also to serve as a basis for discussion of the theoretical ideas tested by the NASA Gyroscope Experiment (aboard Gravity Probe B). Specific accomplishments in the following areas are summarized: the inclusion of electromagnetism into the variational principle for spinning matter, formulation of a self-consistent theory for the case of a fluid in which particle production processes occur, and the derivation of the Raychaudhuri equation in the case of spinning matter.

Ray, John R.↗
DOE OSTI2023

Evaluation of the critical behavior near ferromagnetic to paramagnetic phase transition in CrTe 1-x Se x alloys: An experimental study

For this study, we used the conventional solid-state reaction method to prepare stoichiometric samples of CrTe 1-x Se x , where 0 ≤ x ≤ 0.10, and investigated the structural and critical behavior of the prepared samples. Room temperature powder X-ray diffraction, along with Rietveld refinement, revealed the emergence of the NiAs structure with P6 3 /mmc (194) space group with increasing Se substitution. The high-temperature linear fit to inverse of the dc-susceptibility versus temperature for the mother sample resulted in an effective moment of 3.65μ B Cr with Curie-Weiss temperature near 335K, which is slightly higher than the Tc of ~332K obtained from the inflection point of magnetization versus temperature. Magnetization isotherms were employed to investigate the critical behavior of ferromagnetic CrTe 1-x Se x with 0.0 ≤ x ≤ 0.10 near their Curie temperatures (Tc). The magnetic behavior near Tc was found to follow 3D mean-field critical exponents with a second-order phase transition in all samples investigated. We fine-tuned the critical exponents (β, γ, and δ) using a combination of an iteration technique, the Kouvel-Fisher method, and modified Arrott plots. All samples follow a mean field behavior, with Tc ranging from 298 to 340K. The acquired values exhibit self-consistency, as indicated by the results from the Widom scaling relation. Furthermore, the magnetization isotherms exhibit a universal scaling behavior, providing additional credence to the calculated critical exponents.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2024

Strategies to Obtain Reliable Energy Landscapes from Embedded Multireference Correlated Wavefunction Methods for Surface Reactions

Embedded correlated wavefunction (ECW) theory is a powerful tool for studying ground- and excited-state reaction mechanisms and associated energetics in heterogeneous catalysis. Several factors are important to obtaining reliable ECW energies, critically the construction of consistent active spaces (ASs) along reaction pathways when using a multireference correlated wavefunction (CW) method that relies on a subset of orbital spaces in the configuration interaction expansion to account for static electron correlation, e.g., complete AS self-consistent field theory, in addition to the adequate partitioning of the system into a cluster and environment, as well as the choice of a suitable basis set and number of states included in excited-state simulations. Here, in this work, we conducted a series of systematic studies to develop best-practice guidelines for ground- and excited-state ECW theory simulations, utilizing the decomposition of NH 3 on Pd(111) as an example. We determine that ECW theory results are relatively insensitive to cluster size, the aug-cc-pVDZ basis set provides an adequate compromise between computational complexity and accuracy, and that a fixed-clean-surface approximation holds well for the derivation of the embedding potential. Additionally, we demonstrate that a merging approach, which involves generating ASs from the molecular fragments at each configuration, is preferable to a creeping approach, which utilizes ASs from adjacent structures as an initial guess, for the generation of consistent potential energy curves involving open-d-shell metal surfaces, and, finally, we show that it is essential to include bands of excited states in their entirety when simulating excited-state reaction pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗