SEARCH · Engineering Papers

Results for “Fock states”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Semiclassical approach to atomic and molecular interactions

A general approach, combining quantum and classical mechanics, is used to determine the electron position and velocity distributions in atoms and atomic ions (positive and negative). The Hartree-Fock electronic wave functions and the classical central field approximation are used for evaluation of the dynamic properties of the localized electrons. The distributions, which are of fundamental importance in applications of the binary encounter approximation to description of atomic and ionic collissions, are obtained in the form of simple analytical expressions. The quantum-classical distributions of this work are compared with several other distributions in Ne, Ar, and Al atoms in the ground state.

Kunc, Joseph A.↗

DOE OSTI2020

Distillation of maximally correlated bosonic matter from many-body quantum coherence

We construct quantum coherence resource theories in symmetrized Fock space (QCRTF), thereby providing an information-theoretic framework that connects analyses of quantum coherence in discrete-variable (DV) and continuous variable (CV) bosonic systems. Unlike traditional quantum coherence resource theories, QCRTF can be made independent of the single-particle basis and allow to quantify coherence within and between particle number sectors. For example, QCRTF can be formulated in such a way that neither Bose-Einstein condensates nor Heisenberg-Weyl coherent states are considered as quantum many-body coherence resources, whereas spin-squeezed and quadrature squeezed states are. The QCRTF framework is utilized to calculate the optimal asymptotic distillation rate of maximally correlated bosonic states both for particle number conserving resource states and resource states of indefinite particle number. In particular, we show how to generate a uniform superposition of maximally correlated bosonic states from a state of maximal bosonic coherence with asymptotically unit efficiency using only free operations in the QCRTF.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2019

K-line X-ray fluorescence from highly charged iron ions under dense astrophysical plasma conditions

In the present work, we report an investigation of plasma environment effects on the atomic parameters associated with the K-vacancy states in highly charged iron ions within the astrophysical context of accretion disks around black holes. More particularly, the sensitivity of K-line X-ray fluorescence parameters (wavelengths, radiative transition probabilities, and Auger rates) in Fe XVII–Fe XXV ions has been estimated for plasma conditions characterized by an electron temperature ranging from 105 to 107K and an electron density ranging from 1018 to1022 cm−3. In order to do this, relativistic multiconfiguration Dirac-Fock atomic structure calculations have been carried out by considering a time averaged Debye-Hückel potential for both the electron–nucleus and electron–electron interactions.

Jérôme Deprince↗

NASA NTRS1990

Polyatomic molecular Dirac-Hartree-Fock calculations with Gaussian basis sets

Numerical methods have been used successfully in atomic Dirac-Hartree-Fock (DHF) calculations for many years. Some DHF calculations using numerical methods have been done on diatomic molecules, but while these serve a useful purpose for calibration, the computational effort in extending this approach to polyatomic molecules is prohibitive. An alternative more in line with traditional quantum chemistry is to use an analytical basis set expansion of the wave function. This approach fell into disrepute in the early 1980's due to problems with variational collapse and intruder states, but has recently been put on firm theoretical foundations. In particular, the problems of variational collapse are well understood, and prescriptions for avoiding the most serious failures have been developed. Consequently, it is now possible to develop reliable molecular programs using basis set methods. This paper describes such a program and reports results of test calculations to demonstrate the convergence and stability of the method.

Dyall, Kenneth G.↗

DOE OSTI2022

Multinode Multi-GPU Two-Electron Integrals: Code Generation Using the Regent Language

The computation of two-electron repulsion integrals (ERIs) is often the most expensive step of integral-direct self-consistent field methods. Formally it scales as O(N 4 ), where N is the number of Gaussian basis functions used to represent the molecular wave function. In practice, this scaling can be reduced to O(N 2 ) or less by neglecting small integrals with screening methods. The contributions of the ERIs to the Fock matrix are of Coulomb (J) and exchange (K) type and require separate algorithms to compute matrix elements efficiently. We previously implemented highly efficient GPU-accelerated J-matrix and K-matrix algorithms in the electronic structure code TeraChem. Although these implementations supported the use of multiple GPUs on a node, they did not support the use of multiple nodes. This presents a key bottleneck to cutting-edge ab initio simulations of large systems, e.g., excited state dynamics of photoactive proteins. We present our implementation of multinode multi-GPU J- and K-matrix algorithms in TeraChem using the Regent programming language. Regent directly supports distributed computation in a task-based model and can generate code for a variety of architectures, including NVIDIA GPUs. We demonstrate multinode scaling up to 45 GPUs (3 nodes) and benchmark against hand-coded TeraChem integral code. Finally, we also outline our metaprogrammed Regent implementation, which enables flexible code generation for integrals of different angular momenta.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1984

Potential energy surfaces of MH2 (M = Co, Fe, and Cu)

In many important catalytic processes, the dissociation of the hydrogen molecule is a fundamental step, and the chemisorption of dihydrogen by metal surfaces has been studied. However, quantitative results from theoretical treatments are not very encouraging. Blomberg and Siegbahn (1983) have found that the negative binding energy for NiH2 obtained with the aid of the commonly used Hartree-Fock approximation differs greatly from the result found on the basis of a more sophisticated calculation, including electron correlation. Another study showed the great importance of d orbitals in the dissociation process. A study, similar to the one on NH2, has been conducted for three other metals, including cobalt, iron, and copper. The present investigation is concerned with the results of this study. The potential surfaces of CoH2, FeH2, and CuH2 are discussed, taking into account the lowest atomic states, the bent low spin complex, the bent high spin saddle points, and the linear states.

Siegbahn, P. E. M.↗

NASA NTRS2016

Laboratory Measurements of the K-Shell Transition Energies in L-Shell Ions of Si and S

We have measured the energies of the strongest 1s-2 ℓ (azimuthal quantum number) (ℓ = s, p (s, p are angular momentum states)) transitions in He- through Ne-like silicon and sulfur ions to an accuracy of less than 1 electronvolt using the Lawrence Livermore National Laboratory's electron beam ion traps, EBIT-I and SuperEBIT, and the NASA/GSFC EBIT Calorimeter Spectrometer (ECS). We identify and measure the energies of 18 and 21 X-ray features from silicon and sulfur, respectively. The results are compared to new Flexible Atomic Code calculations and to semi-relativistic Hartree-Fock calculations by Palmeri et al. (2008). These results will be especially useful for wind diagnostics in high-mass X-ray binaries, such as Vela X-1 and Cygnus X-1, where high-resolution spectral measurements using Chandraʼs high-energy transmission grating has made it possible to measure Doppler shifts of 100 kilometers per second. The accuracy of our measurements is consistent with that needed to analyze Chandra observations, exceeding Chandraʼs 100 kilometers per second limit. Hence, the results presented here not only provide benchmarks for theory, but also accurate rest energies that can be used to determine the bulk motion of material in astrophysical sources. We show the usefulness of our results by applying them to redetermine Doppler shifts from Chandra observations of Vela X-1.

atomic data↗

DOE OSTI2022

The roles of three-nucleon force and continuum coupling in mirror symmetry breaking of oxygen mass region

With both three-nucleon force and continuum coupling included, we have developed a self-consistent ab initio Gamow shell model within the Gamow Hartree-Fock (GHF) basis obtained by the realistic interaction itself. With the chiral two-nucleon N 3 LO and three-nucleon N 2 LO interactions, the Gamow shell model has been applied to the mirror systems of Z = 8 neutron-rich isotopes and N = 8 proton-rich isotones, giving good agreements with data in binding energies, dripline positions and excitation spectra. The GHF calculated that the 0 ds/2 , 1 S1/2 and 1 p3/2 orbitals are resonances. The resonance states and their interplay with nonresonant continua play a crucial role in the descriptions of nuclei around driplines. Excitation spectra and Thomas-Ehrman shifts observed can be better described when both three-nucleon force and continuum coupling are considered in calculations. The three-nucleon force and continuum coupling produce a combined effect on the Thomas-Ehrman shift, e.g., for the 1/2 + resonance level of 19Na. The calculations help the understandings of related nuclear astrophysical processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

State Preparation of Antisymmetrized Geminal Power on a Quantum Computer without Number Projection

The antisymmetrized geminal power (AGP) is equivalent to the number projected Bardeen–Cooper–Schrieffer (PBCS) wave function. It is also an elementary symmetric polynomial (ESP) state. We generalize previous research on deterministically implementing the Dicke state to a state preparation algorithm for an ESP state, or equivalently AGP, on a quantum computer. Our method is deterministic and has polynomial cost, and it does not rely on number symmetry breaking and restoration. We also show that our circuit is equivalent to a disentangled unitary paired coupled cluster operator and a layer of unitary Jastrow operator acting on a single Slater determinant. Here, the method presented herein highlights the ability of disentangled unitary coupled cluster to capture nontrivial entanglement properties that are hardly accessible with traditional Hartree–Fock based electronic structure methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Utilizing Essential Symmetry Breaking in Auxiliary-Field Quantum Monte Carlo: Application to the Spin Gaps of the C 36 Fullerene and an Iron Porphyrin Model Complex

We present three distinct examples where phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) can be reliably performed with a single-determinant trial wave function with essential symmetry breaking. Essential symmetry breaking was first introduced by Lee and Head-Gordon [ Phys. Chem. Chem. Phys. 2019, 21, 4763-4778, 10.1039/C8CP07613H]. We utilized essential complex and time-reversal symmetry breaking with ph-AFQMC to compute the triplet–singlet energy gap in the TS12 set. We found statistically better performance of ph-AFQMC with complex-restricted orbitals than with spin-unrestricted orbitals. We then showed the utilization of essential spin symmetry breaking when computing the singlet–triplet gap of a known biradicaloid, C 36 . ph-AFQMC with spin-unrestricted Hartree–Fock (ph-AFQMC+UHF) fails catastrophically even with spin-projection and predicts no biradicaloid character. With approximate Brueckner orbitals obtained from regularized orbital-optimized second-order Møller–Plesset perturbation theory (κ-OOMP2), ph-AFQMC quantitatively captures strong biradicaloid character of C 36 . Lastly, we applied ph-AFQMC to the computation of the quintet–triplet gap in a model iron porphyrin complex where brute-force methods with a small active space fail to capture the triplet ground state. We show unambiguously that neither triplet nor quintet is strongly correlated using UHF, κ-OOMP2, and coupled-cluster with singles and doubles (CCSD) performed on UHF and κ-OOMP2 orbitals. There is no essential symmetry breaking in this problem. By virtue of this, we were able to perform UHF+ph-AFQMC reliably with a cc-pVTZ basis set and predicted a triplet ground state for this model geometry. The largest ph-AFQMC in this work correlated 186 electrons in 956 orbitals. Our work highlights the utility, scalability, and accuracy of ph-AFQMC with a single-determinant trial wave function with essential symmetry breaking for systems mainly dominated by dynamical correlation with little static correlation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Shadow energy functionals and potentials in Born–Oppenheimer molecular dynamics

In Born–Oppenheimer molecular dynamics (BOMD) simulations based on the density functional theory (DFT), the potential energy and the interatomic forces are calculated from an electronic ground state density that is determined by an iterative self-consistent field optimization procedure, which, in practice, never is fully converged. The calculated energies and forces are, therefore, only approximate, which may lead to an unphysical energy drift and instabilities. Here, we discuss an alternative shadow BOMD approach that is based on backward error analysis. Instead of calculating approximate solutions for an underlying exact regular Born–Oppenheimer potential, we do the opposite. Instead, we calculate the exact electron density, energies, and forces, but for an underlying approximate shadow Born–Oppenheimer potential energy surface. In this way, the calculated forces are conservative with respect to the approximate shadow potential and generate accurate molecular trajectories with long-term energy stabilities. We show how such shadow Born–Oppenheimer potentials can be constructed at different levels of accuracy as a function of the integration time step, δt, from the constrained minimization of a sequence of systematically improvable, but approximate, shadow energy density functionals. For each energy functional, there is a corresponding ground state Born–Oppenheimer potential. These pairs of shadow energy functionals and potentials are higher-level generalizations of the original “zeroth-level” shadow energy functionals and potentials used in extended Lagrangian BOMD. The proposed shadow energy functionals and potentials are useful only within this extended dynamical framework, where also the electronic degrees of freedom are propagated as dynamical field variables together with the atomic positions and velocities. The theory is quite general and can be applied to MD simulations using approximate DFT, Hartree–Fock, or semi-empirical methods, as well as to coarse-grained flexible charge models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Zero-pairing and zero-temperature limits of finite-temperature Hartree-Fock-Bogoliubov theory

Recently, variational Hartree-Fock-Bogoliubov (HFB) mean-field equations were shown to possess a mathematically well-defined zero-pairing limit, independently of the closed- or open-shell character of the system under consideration. This limit is non-trivial for open-shell systems such that HFB theory does {\it not} reduce to the Hartree-Fock (HF) formalism in all cases. The present work extends this analysis to finite-temperature HFB (FTHFB) theory by investigating the behavior of this more general formalism in the combined zero-temperature and zero-pairing limits. The zero-pairing and zero-temperature limits of the FTHFB statistical density operator constrained to carry an arbitrary (integer) number of particles A on average is worked out analytically and realized numerically using a two-nucleon interaction. While the FTHFB density operator reduces to the projector corresponding to a pure HF Slater determinant for closed-shell nuclei, the FTHFB formalism does not reduce to the HF theory in all cases in the zero-temperature and zero-pairing limits, i.e. for open-shell nuclei. However, the fact that a nucleus can be of open-shell character in these joint limits is necessarily the result of some symmetry restrictions. Whenever it is the case, the non-trivial description obtained for open-shell systems is shown to depend on the order with which both limits are taken, i.e. the two limits do not commute for these systems. When the zero-temperature limit is performed first, the FTHFB density operator is demoted to a projector corresponding to a pure state made out of a linear combination of a finite number of Slater determinants with different (even) numbers of particles. When the zero-pairing limit is performed first, the FTHFB density operator remains a statistical mixture of a finite number of Slater determinants with both even and odd particle numbers. While the entropy (pairing density) is zero in the first (second) case, it does not vanish in the second (first) case in spite of the temperature (pairing) tending towards zero. The difference between both limits can have striking consequences for the (thermal) expectation values of observables. For instance, the particle-number variance does not vanish in either case and has limiting values that differ by a factor of two in both cases. In conclusion, while in the textbook situation associated with closed-shell nuclei Hartree-Fock-Bogoliubov (finite-temperature Hartree-Fock) theory reduces to Hartree-Fock theory in the zero-pairing (zero-temperature) limit, the present analysis demonstrates that a non trivial and unexpected limit is obtained for this formalism in open-shell systems. This result sheds a new light on certain aspects of this otherwise very well-studied many-body formalism.

20 ≤ A ≤ 38↗

DOE OSTI2024

Thermal mean-field theories

Several closely related ab initio thermal mean-field theories for fermions, both well-established and new ones, are compared with one another at the formalism level and numerically. The theories considered are Fermi–Dirac theory; thermal Hartree–Fock (HF) theory; two modifications of the thermal single-determinant and the first-order finite-temperature many-body perturbation theory based on a zero-temperature or thermal HF reference. Furthermore, thermal full-configuration-interaction theory is used as the benchmark.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2025

Data-Driven Recommendation of Optimal Tuning Scheme for Range-Separated Hybrid Functionals in Solution-Phase UV/Vis Absorption Energy Prediction

Time-dependent density functional theory (TDDFT) combined with range-separated hybrid (RSH) functionals and a tuned range-separation parameter γ offers a computationally economical approach for high-throughput excited- state property predictions. The γ-tuning procedure in the gas phase is well established. However, no agreement on the best γ- tuning procedure has been made when considering the solvent effect with implicit solvent models like the polarizable continuum model (PCM). To answer that question, this study created a diverse dataset with 937 molecules with experimental solutionphase UV/vis absorption spectra. Three γ-tuning methods, the gasphase γ-tuning (GPγT), the partial vertical γ-tuning (PVγT), and the strict vertical γ-tuning (SVγT), were evaluated for the ωPBEh functional over the entire dataset. Additional benchmarks are done for the optimally tuned screened range-separated hybrid combined with the PCM approach (SRSH-PCM) and the solvation-mediated tuning procedure (sol-med-OT). Our findings revealed that the optimal γ-values obtained by the PVγT and the SVγT are significantly smaller than the GPγT. This trend holds consistently across all molecules in our dataset, and we explained the origin of this phenomenon. TDDFT calculations with PVγTand SVγT-tuned γ-values and default global Fock exchange fraction achieve superior performance compared to those using GPγTtuned or default γ and slightly outperform SRSH-PCM and sol-med-OT with similar or lesser computational cost. Furthermore, we found that the smaller γ-values from SVγT captured the expected 1/(εR) asymptotic behavior in the solution phase, resulting in accurate prediction of solution-phase CT excitations, consistent with the screened asymptote behavior encoded in SRSH-PCM. These results show that SVγT is the best scheme for high-throughput UV/vis absorption spectrum calculations using the ωPBEh functional from a data-driven perspective.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Real-Time Time-Dependent Nuclear–Electronic Orbital Approach: Dynamics beyond the Born–Oppenheimer Approximation

The quantum mechanical treatment of both electrons and nuclei is crucial in nona- diabatic dynamical processes such as proton-coupled electron transfer. The nuclear electronic orbital (NEO) method provides an elegant framework for including nuclear quantum effects beyond the Born-Oppenheimer approximation. To enable the study of nonequilibrium properties, we derive and implement a real-time NEO (RT-NEO) approach based on time-dependent Hatree-Fock or density functional theory, in which both the electronic and nuclear degrees of freedom are propagated in a time-dependent variational framework. Nuclear and electronic spectral features can be resolved from the time-dependent dipole moment computed using the RT-NEO method. The test cases demonstrate the dynamical interplay between the quantum nuclei and the electrons through vibronic coupling. Moreover, vibrational excitation in the RT-NEO approach is demonstrated by applying a resonant driving eld, and electronic excitation is demonstrated by simulating excited state intramolecular proton transfer. This work shows that the RT-NEO approach is a promising tool to study nonadiabatic quantum dynamical processes within a time-dependent variational description for the coupled electronic and nuclear degrees of freedom.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2006

New Accurate Oscillator Strengths and Electron Excitation Collision Strengths for N1

The nonorthogonal orbitals technique in a multiconfiguration Hartree-Fock approach is used to calculate oscillator strengths and transition probabilities of N(I) lines. The relativistic effects are allowed by means of Breit-Pauli operators. The length and velocity forms of oscillator strengths show good agreement for most transitions. The B-spline R-matrix with pseudostates approach has been used to calculate electron excitation collision strengths and rates. The nonorthogonal orbitals are used for an accurate description of both target wave functions and the R-matrix basis functions. The 24 spectroscopic bound and autoionizing states together with 15 pseudostates are included in the close-coupling expansion. The collision strengths for transitions between fine-structure levels are calculated by transforming the LS-coupled K-matrices to K-matrices in an intermediate coupling scheme. Thermally averaged collision strengths have been determined by integrating collision strength over a Maxwellian distribution of electron energies over a temperature range suitable for the modeling of astrophysical plasmas. The oscillator strengths and thermally averaged collision strengths are presented for transitions between the fine-structure levels of the 2s(sup 2)p(sup 3) (sup 4)S(sup 0), (sup 2)D(sup 0), (sup 2)P(sup 0), 2s2p(sup 4) (sup 4)P, 2s(sup 2)2p(sup 2)3s (sup 4)P, and (sup 2)P terms and from these levels to the levels of the 2s(sup 2)2p(sup 2)3p (sup 2)S(sup 0), (sup 4)D(sup 0), (sup 4)P(sup 0), (sup 4)S(sup 0), (sup 2)D(sup 0), (sup 2)P(sup 0),2s(sup 2)2p(sup 2)3s(sup 2)D, 2s(sup 2)2p(sup 2)4s(sup 4)P, (sup 2)P, 2s(sup 2)2p(sup 2)3d(sup 2)P, (sup 4)F,(sup 2)F,(sup 4)P, (sup 4)D, and (sup 2)D terms. Thermally averaged collision strengths are tabulated over a temperature range from 500 to 50,000 K.

Tayal, S. S.↗

DOE OSTI2025

A Meta-Generalized Gradient Approximation for the Cavity-Dependent Exchange-Correlation Interaction in Strongly Coupled Light–Matter Systems

Strong light–matter coupling in optical cavities enables the manipulation of chemical and physical properties without altering molecular composition. Theoretical modeling of such phenomena requires exchange-correlation (XC) functionals that account for both electron–electron and electron–photon (ep) interactions within quantum electrodynamical density functional theory (QEDFT). In this work, we develop a meta-generalized gradient approximation (meta-GGA) specifically targeting the cavity-dependent XC interaction in strongly coupled light–matter systems. This novel approximation is built upon a new semilocal polarizability approximation, which draws from the jellium-with-a-gap model, and can be extended to a “global hybrid” variant that goes beyond the isotropic model from previous approximations. The polarizability model yields significantly improved dispersion coefficients and benchmark calculations with the cavity-dependent XC functional demonstrate improved agreement with QED Hartree–Fock (QED-HF) reference energies. Application to the regioselectivity of brominated nitrobenzene intermediates reveals the functional’s capacity to capture cavity-induced energetic shifts. In conclusion, our results advance the Jacob’s ladder of functionals for QEDFT and provide a practical tool for modeling polaritonic chemistry.

Approximation↗

DOE OSTI2023

Extracting model-independent nuclear level densities away from stability

The nuclear level density (NLD) is a fundamental measure of the complex structure of atomic nuclei at relatively high energies. Here, in this study, we present the first model-independent measurement of the absolute partial NLD for a short-lived nucleus. For this purpose we adapt the recently introduced “shape method” for β-decay experiments, providing the shape of the γ-ray strength function for exotic nuclei. In this work, we show that combining the shape method with the β-Oslo technique allows for the extraction of the NLD of the populated states without the need for theoretical input. This development opens the way for the extraction of experimental NLDs far from stability with major implications in astrophysical and other applications. We benchmark our approach using data for the stable 76 Ge nucleus, finding excellent agreement with previous experimental results. In addition, we present new experimental data and determine the absolute partial level density for the short-lived 88 Kr nucleus. Our results suggest a fivefold increase in the NLD for the case of 88 Kr, compared to the recommended values from semimicroscopic Hartree-Fock Bogoliubov calculations recommended by the RIPL3 nuclear data library. However, our results are in good agreement with other semimicroscopic level density models. We demonstrate the impact of our method on the 87 Kr(n, γ) neutron capture rate and show that our experimental uncertainties for NLDs fulfill the requirements needed for astrophysical calculations predicting r-process abundances.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗