Engineering Papers⌕ Search

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.

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

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The approximate second order coupled-cluster method based on a size-consistent Brillouin–Wigner partitioning

We present a variant of the approximate second order coupled-cluster method (CC2) with a two-parameter size-consistent Brillouin–Wigner (BW-s) partitioning instead of a Møller–Plesset (MP) partitioning for the unperturbed Hamiltonian, which we refer to as BWs-CC2. The computational complexity of this model scales identically to CC2 with molecular size. Conventional CC2 and its regularized BWs-CC2 variants, as well as conventional MP2 and two of its regularized BW-s2 variants, were assessed on a 535 element database spanning thermochemistry, non-covalent interactions, barrier heights, and isomerization energies. To ensure a well-defined model chemistry, the assessment was performed using internally stable spin-polarized Hartree–Fock (HF) orbitals in the finite aug-cc-pVQZ basis without counterpoise corrections. As a result of using stable orbitals, contrary to conventional wisdom, we find that CC2 substantially outperforms MP2 on molecules with significantly spin contaminated reference orbitals without a significant increase in error on systems with a spin-pure reference, showing the value of its single substitutions. While no single choice of regularization parameters can be optimal for all datasets, we find that BWs-CC2 generally outperforms both CC2 and BW-s2 with a single judicious parameter choice. Additional tests on dipole moments and bond lengths of diatomics provide further support for the utility of this choice. Furthermore, the main outliers and poorest performing cases are associated with large amounts of spin-contamination in the HF reference, which is indicative of systems with either strong correlation or extensive artificial symmetry breaking. Overall, these findings argue that the perception of the quality of the CC2 ground state should be reevaluated and that it can be further improved upon by the soundly based BWs-CC2 variant with the recommended parameter choice.

Correlation energy↗
DOE OSTI2020

Performance of an atomic mean-field spin–orbit approach within exact two-component theory for perturbative treatment of spin–orbit coupling

The paper reports development of an atomic mean-field (AMF) spin–orbit approach within the spin-free exact two-component theory in its one-electron variant (SFX2C-1e), a new relativistic quantum-chemical method for treating spin-orbit coupling in molecules. The effective one-electron spin–orbit integral matrix in the four-component representation is assembled as a direct sum of one-centre spin–orbit integral matrices with the mean-field two-electron contributions evaluated using atomic SFX2C-1e Hartree–Fock density matrices. It is then transformed into two-component representation using analytic SFX2C-1e energy derivative formulation. The resulting two-component spin–orbit integral matrix is by design suitable for use in perturbative calculations of spin–orbit coupling, treating SFX2C-1e wavefunctions as unperturbed states. The accuracy of the present AMF approach has been demonstrated using benchmark calculations of spin–orbit splittings for representative diatomic radicals at the equation-of-motion coupled-cluster singles and doubles level. In conclusion, to demonstrate the applicability and accuracy of the present perturbative spin–orbit scheme in calculations of challenging heavy-element containing systems, a thorough computational investigation of six low-lying electronic states of ThO+ is reported.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Quasiparton distributions in massive QED2: Toward quantum computation

We analyze the quasiparton distributions of the lightest 𝜂′ meson in massive two-dimensional quantum electrodynamics (QED2) by exact diagonalization. The Hamiltonian and boost operators are mapped onto spin qubits in a spatial lattice with open boundary conditions. The lowest excited state in the exact diagonalization is shown to interpolate continuously between an anomalous 𝜂′ state at strong coupling, and a nonanomalous heavy meson at weak coupling, with a cusp at the critical point. The boosted 𝜂′ state follows relativistic kinematics but with large deviations in the luminal limit. The spatial quasiparton distribution function and amplitude for the 𝜂′ state are computed numerically for increasing rapidity both at strong and weak coupling, and compared to the exact light front results. The numerical results from the boosted form of the spatial parton distributions, compare fairly with the inverse Fourier transformation of the luminal parton distributions, derived in the lowest Fock space approximation. Our analysis points out some of the limitations facing the current lattice program for the parton distributions.

Lattice field theory↗
DOE OSTI2023

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

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

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

Unusual behavior of Cooper minima of n s subshells in high- Z atoms

Here, a study of Cooper minima (CM) arising from the photoionization of 6$\textit{s}$, 5$\textit{s}$, and 4$\textit{s}$ subshells of high-$\textit{Z}$ atoms has been performed using Dirac-Fock (DF), two-channel relativistic-random-phase approximation (RRPA), and fully coupled RRPA. The results show huge splittings between $ns → εp_{3/2}$ and $ns → εp_{1/2}$ CM which increase with $\textit{Z}$ owing primarily to the relativistic interactions (spin orbit) that are attractive for the $εp_{1/2}$ final state but repulsive for the corresponding $εp_{3/2}$. In addition, it was found that correlation in the form of interchannel coupling (essentially configuration interaction in the final continuum states) plays a huge role in determining the location of the CM. For 6$\textit{s}$ photoionization, the 6$s → εp_{3/2}$ and 6$s → εp_{1/2}$ CM behave completely different as a function of $\textit{Z}$ . It was also found that for 5$\textit{s}$ and 4$\textit{s}$ photoionization, the CM move below the threshold, with increasing $\textit{Z}$ , and, at high enough $\textit{Z}$ , the 5s → εp3/2 and 4s → εp3/2 CM re-emerge into the continuum. The calculations have been carried out for the $\textit{ns}$ subshells of Hg ($\textit{Z}$ = 80), Rn ($\textit{Z}$ = 86), Ra ($\textit{Z}$ = 88), No ($\textit{Z}$ = 102), Cn ($\textit{Z}$ = 112), and Og ($\textit{Z}$ = 118).

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1992

L-shell X-ray opacity of many-electron atoms

The study presents a formulation for obtaining the L-shell photoelectric cross section (sigma sub L) for many-electron atomic systems. A scaling law is suggested for the evaluation of the contribution to sigma sub L from 2s and 2p ejection from a general system in any ionization state. The formula is expressed in terms of the effective nuclear charge (Z sub 2s and Z sub 2p) seen by the corresponding bound electrons, with Z sub 2s and Z sub 2p, obtained from improved screening constants derived in a recent investigation. A comparison is made of Hartree-Fock and hydrogenic expressions for sigma sub L with experimental results for atomic neon at X-ray photon energies, providing the basis for the scaling law and an approximate energy dependence for the cross sections. The special case of sigma sub L for atomic Fe is considered in detail. This L edge has structure because of the different threshold for 2s and 2p photoejection; the latter has fine structure.

Jung, Young-Dae↗
DOE OSTI2024

Modeling Multi-Step Organic Reactions: Can Density Functional Theory Deliver Misleading Chemistry?

Many organic reactions are characterized by a complex mechanism with a variety of transition states and intermediates of different chemical natures. Their correct and accurate theoretical characterization critically depends on the accuracy of the computational method used. In this work, we study a complex ambimodal cycloaddition with five transition states, two intermediates, and three products, and we ask whether density functional theory (DFT) can provide a correct description of this type of complex and multifaceted reaction. Our work fills a gap in that most systematic benchmarks of DFT for chemical reactions have considered much simpler reactions. Our results show that many density functionals not only lead to seriously large errors but also differ from one another in predicting whether the reaction is ambimodal. Only a few of the available functionals provide a balanced description of the complex and multifaceted reactions. The parameters varied in the tested functionals are the ingredients, the treatment of medium-range and nonlocal correlation energy, and the inclusion of Hartree–Fock exchange. Furthermore, these results show a clear need for more benchmarks on the mechanisms of large molecules in complex reactions.

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

Breakdown of helical edge state topologically protected conductance in time-reversal-breaking excitonic insulators

Gapless helical edge modes are a hallmark of the quantum spin Hall effect. Protected by time-reversal symmetry, each edge contributes a quantized zero-temperature conductance quantum G 0 ≡ e 2 /h. However, the experimentally observed conductance in WT e2 decreases below G 0 per edge already at edge lengths around 100 nm, even in the absence of explicit time-reversal breaking due to an external field or magnetic impurities. Here we show how a time-reversal breaking excitonic condensate with a spin-spiral order that can form in WT e2 leads to the breakdown of conductance quantization. We perform Hartree-Fock calculations to compare time-reversal breaking and preserving excitonic insulators. Using these mean-field models we demonstrate via quantum transport simulations that weak nonmagnetic disorder reproduces the edge length scaling of resistance observed in the experiments. We complement this by analysis in the Luttinger liquid picture, shedding additional light on the mechanism behind the quantization breakdown.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Physics-Driven Construction of Compact Primitive Gaussian Density Fitting Basis Sets

We present a model-assisted density fitting (MADF) basis set generator, an algorithm for generating primitive atomic Gaussian density fitting (DF) basis sets (DFBSs) from a contracted Gaussian orbital basis set (OBS). The MADF algorithm produces DFBSs suitable for accurate robust DF approximation of 2-particle interactions in mean-field and correlated electronic structures. The algorithm is designed to (a) saturate the OBS product space by a large regularized set of primitive solid-harmonic Gaussian shells with nonuniform distribution of exponents, followed by (b) pruning of the shells according to their contributions to the 2- body energy of a correlated atomic ensemble. Building the DFBS generator model almost exclusively on mathematical and physical principles allows one to limit the number of parameters that control the density fitting error to three, with a single set of parameters sufficient for computations with all basis cardinal numbers, with and without correlation of core electrons, with and without scalar and spin-dependent relativistic effects, spanning almost all of the Periodic Table. Performance assessment included basis sets up to quadruple-ζ quality from several major basis set families, using molecules composed of main-group, d-block, and f-block elements. The resulting DF errors in Hartree−Fock and second-order MP2 energies (with relativistic all-electron treatments, when appropriate) were on the order of 20 and 10 μE h per electron, respectively.

Approximation↗
DOE OSTI2023

Shape evolution in even-mass 98-104 Zr isotopes via lifetime measurements using the $\gamma \gamma $-coincidence technique

The Zirconium (Z = 40) isotopic chain has attracted interest for more than four decades. The abrupt lowering of the energy of the first 2 + state and the increase in the transition strength B(E2; 2$^+_1$ →0$^+_1$) going from 98 Zr to 100 Zr has been the first example of “quantum phase transition” in nuclear shapes, which has few equivalents in the nuclear chart. Although a multitude of experiments have been performed to measure nuclear properties related to nuclear shapes and collectivity in the region, none of the measured lifetimes were obtained using the Recoil Distance Doppler Shift method in the γγ-coincidence mode where a gate on the direct feeding transition of the state of interest allows a strict control of systematical errors. Here this work reports the results of lifetime measurements for the first yrast excited states in 98-104 Zr carried out to extract reduced transition probabilities. The new lifetime values in γγ-coincidence and γ-single mode are compared with the results of former experiments. Recent predictions of the Interacting Boson Model with Configuration Mixing, the Symmetry Conserving Configuration Mixing model based on the Hartree–Fock–Bogoliubov approach and the Monte Carlo Shell Model are presented and compared with the experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Axionic dark matter halos in the gravitational field of baryonic matter

We consider a dark matter halo (DMH) of a spherical galaxy as a Bose–Einstein condensate (BEC) of the ultra-light axions (ULA) interacting with the baryonic matter. In the mean-field (MF) limit, we have derived the integro-differential equation of the Hartree–Fock type for the spherically symmetrical wave function of the DMH component. This equation includes two independent dimensionless parameters: (i) β is the ratio of baryon and axion total mases and (ii) ξ is the ratio of characteristic baryon and axion spatial parameters. We extended our “dissipation algorithm” for studying numerically the ground state of the axion halo in the gravitational field produced by the baryonic component. We calculated the characteristic size, x c of DMH as a function of β and ξ and obtained an analytical approximation for x c .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Exact-two-component block-localized wave function: A simple scheme for the automatic computation of relativistic ΔSCF

Block-localized wave function is a useful method for optimizing constrained determinants. In this article, we extend the generalized block-localized wave function technique to relativistic two-component framework. Optimization of excited state determinants for two-component wave functions presents a unique challenge because the excited state manifold is often quite dense with degenerate states. Furthermore, we test the degree to which certain symmetries result naturally from the ΔSCF optimization such as time reversal symmetry and symmetry with respect to the total angular momentum operator on a series of atomic systems. Variational optimizations may often break symmetry in order to lower the overall energy, just as unrestricted Hartree-Fock breaks spin symmetry. Overall, we demonstrate that time reversal symmetry is roughly maintained when using Hartree-Fock, but less so when using Kohn-Sham density functional theory. Additionally, maintaining total angular momentum symmetry appears to be system dependent and not guaranteed. Finally, we were able to trace the breaking of total angular momentum symmetry to the relaxation of core electrons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1980

Two-photon excitation cross-section in light and intermediate atoms

The method of explicit summation over the intermediate states is used along with LS coupling to derive an expression for two-photon absorption cross section in light and intermediate atoms in terms of integrals over radial wave functions. Two selection rules, one exact and one approximate, are also derived. In evaluating the radial integrals, for low-lying levels, the Hartree-Fock wave functions, and for high-lying levels, hydrogenic wave functions obtained by the quantum defect method are used. A relationship between the cross section and the oscillator strengths is derived. Cross sections due to selected transitions in nitrogen, oxygen, and chlorine are given. The expression for the cross section is useful in calculating the two-photon absorption in light and intermediate atoms.

Omidvar, K.↗
DOE OSTI2025

Measurement of 𝐷-line energies in sodiumlike Ir

We report measurements of the 𝐷 1 and 𝐷 2 transition energies in Na-like Ir (Ir 66+ ). The 𝐷 1 3⁢𝑠−3⁢𝑝 1/2 transition energy, measured in the extreme ultraviolet (EUV) region at 169.977 ± 0.005 eV with a relative uncertainty of 29 ppm, represents one of the most precise 𝐷-line energy measurements to date for high atomic numbers 𝑍 ≥ 30. The 𝐷 2 3⁢𝑠−3⁢𝑝 3/2 transition energy, measured in the x-ray region, is 621.11 ± 0.06 eV. Therefore, there is a fine-structure separation of 451.13 ± 0.06 eV for the Na-like Ir 3⁢𝑝 levels. We present theoretical calculations of these energies using quantum electrodynamic-extended relativistic many-body perturbation theory and multiple-configuration Dirac-Hartree-Fock methods, and discuss the discrepancies between experimental and theoretical values. The prospect of determining the absolute nuclear charge radii of heavy elements, including rare isotopes, from these measurements is explored.

Atomic spectra↗
DOE OSTI2023

SQMBox: Interfacing a semiempirical integral library to modular ab initio electronic structure enables new semiempirical methods

Ab initio and semiempirical electronic structure methods are usually implemented in separate software packages or use entirely different code paths. As a result, it can be time-consuming to transfer an established ab initio electronic structure scheme to a semiempirical Hamiltonian. Here we present an approach to unify ab initio and semiempirical electronic structure code paths based on a separation of the wavefunction ansatz and the needed matrix representations of operators. With this separation, the Hamiltonian can refer to either an ab initio or semiempirical treatment of the resulting integrals. We built a semiempirical integral library and interfaced it to the GPU-accelerated electronic structure code TeraChem. Equivalency between ab initio and semiempirical tight-binding Hamiltonian terms is assigned according to their dependence on the one-electron density matrix. The new library provides semiempirical equivalents of the Hamiltonian matrix and gradient intermediates, corresponding to those provided by the ab initio integral library. This enables the straightforward combination of semiempirical Hamiltonians with the full pre-existing ground and excited state functionality of the ab initio electronic structure code. We demonstrate the capability of this approach by combining the extended tight-binding method GFN1-xTB with both spin-restricted ensemble-referenced Kohn–Sham and complete active space methods. We also present a highly efficient GPU implementation of the semiempirical Mulliken-approximated Fock exchange. The additional computational cost for this term becomes negligible even on consumer-grade GPUs, enabling Mulliken-approximated exchange in tight-binding methods for essentially no additional cost.

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

Accurate localization of Kosterlitz-Thouless-type quantum phase transitions for one-dimensional spinless fermions

We investigate the charge-density wave (CDW) transition for one-dimensional spinless fermions at half band filling with nearest-neighbor electron transfer amplitude t and interaction V. The model is equivalent to the anisotropic XXZ Heisenberg model for which the Bethe Ansatz provides an exact solution. For V>V c =2t, the CDW order parameter and the single-particle gap are finite but exponentially small, as is characteristic for a Kosterlitz-Thouless transition. It is notoriously difficult to locate such infinite-order phase transitions in the phase diagram using approximate analytical and numerical approaches. Second-order Hartree-Fock theory is qualitatively applicable for all interaction strengths, and predicts the CDW transition to occur at V$^{(2)}_{c,2}$≈1.5t. Second-order Hartree Fock theory is almost variational because the density of quasiparticle excitations is small. We apply the density-matrix renormalization group (DMRG) for periodic boundary conditions for system sizes up to 514 sites, which permits a reliable extrapolation of all physical quantities to the thermodynamic limit, apart from the critical region. We investigate the ground-state energy, the gap, the order parameter, the momentum distribution, the quasiparticle density, and the density-density correlation function to locate V c from the DMRG data. In conclusion, tracing the breakdown of the Luttinger liquid and the peak in the quasiparticle density at the band edge permits us to reproduce V c with an accuracy of one percent.

1-dimensional spin chains↗
DOE OSTI2020

Exact sum rules with approximate ground states

Electromagnetic and weak transitions tell us a great deal about the structure of atomic nuclei. Yet modeling transitions can be difficult: it is often easier to compute the ground state, if only as an approximation, than excited states. One alternative is through transition sum rules, in particular the non-energy-weighted and energy-weighted sum rules, which can be computed as expectation values of operators. We investigate by computing sum rules for a variety of nuclei, comparing the numerically exact full configuration-interaction shell model, as a reference, to Hartree-Fock, projected Hartree-Fock, and the nucleon pair approximation. Furthermore, these approximations yield reasonable agreement, which we explain by prior work on the systematics of transition moments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

CP2K: An Electronic Structure and Molecular Dynamics Software Package - Quickstep: Efficient and Accurate Electronic Structure Calculations

CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular and biological systems. It is especially aimed at massively-parallel and linear-scaling electronic structure methods and state-of-the-art ab-initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension. TDK has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No. 716142). VRR has been supported by the Swiss National Science Foundation in the form of Ambizione grant No. PZ00P2 174227 and RZK by the Natural Sciences and Engineering Research Council of Canada (NSERC) through Discovery Grants (RGPIN-2016-0505). GKS and CJM are supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. UK based work was funded under the embedded CSE programme of the ARCHER UK National Supercomputing Service (http://www.archer.ac.uk), grants eCSE03-011, eCSE06-6, eCSE08-9, eCSE13-17 and the EPSRC (EP/P022235/1) grant “Surface and Interface Toolkit for the Materials Chemistry Community". Computational resources were provided by the Swiss National Supercomputing Centre (CSCS) and Compute Canada. The generous allocation of computing time on the FPGA-based supercomputer “Noctua" at PC2 is kindly acknowledged.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗