Search NASA⌕ Search

SEARCH · Search NASA

Results for “Relativistic wavefunctions”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

X-ray absorption spectra of f-element complexes: insight from relativistic multiconfigurational wavefunction theory

X-ray absorption near edge structure (XANES) spectroscopy, coupled with ab initio calculations, has emerged as the state-of-the-art tool for elucidating the metal-ligand bonding in f-element complexes. Furthermore, this highlight presents recent efforts in calculating XANES spectra of lanthanide and actinide compounds with relativistic multiconfiguration wavefunction approaches that account for differences in donation bonding in the ground state (GS) versus a core-excited state (ES), multiplet effects, and spin-orbit-coupling. With the GS and ES wavefunctions available, including spin-orbit effects, an arsenal of chemical bonding tools that are popular among chemists can be applied to rationalize the observed intensities in terms of covalent bonding.

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

Magnus method for electronic structure calculations at extreme conditions

We present the application of Magnus based methods to the solution of first order coupled ordinary differential equations in High Energy Density (HED) physics applications. Our focus is on the application to quantum mechanical methods, specifically on the solution of the radial Dirac equation for real and complex energies. HED applications require accurate solutions across a wide range of spatial and energy domains, including regimes where the solutions exhibit pronounced oscillatory behavior. Such cases pose significant computational challenges. We demonstrate that Magnus-based integrators can efficiently and accurately address these challenges. We discuss the implementation of the Magnus method for the solution of the radial Dirac equation, including practical considerations such as the evaluation of matrix exponentials, numerical integration, error estimation, and adaptive step size control. We also discuss the application of these methods to complex energy Green’s function techniques and the efficient approximation of integrals of the solutions relevant to HED electronic structure calculations. Here, we demonstrate the accuracy and robustness of the resulting method in applications to the free-particle case, for which analytic solutions are available for comparison, as well as the challenging case of gold at HED conditions.

general physics↗
NASA NTRS1995

Kramers' Restricted Closed Shell CCSD Theory

A Kramers' restricted version of the closed shell coupled cluster singles doubles theory is presented. The theory may be used in conjunction with 2 or 4-component relativistic reference wavefunctions. The intrinsic treatment of the spin-orbit coupling doubles the number of independent quantities (amplitudes and integrals) relative to a spin-independent formalism. The number of operations required to evaluate the equations is four times larger than in the optimal spin-independent closed shell formalism.

Visscher, Lucas↗
DOE OSTI2025

A comprehensive approach for elucidating the interplay between 4f n +1 and 4f n 5d 1 configurations in Ln 2+ complexes

Lanthanides (Ln) are typically found in the +3 oxidation state. However, in recent decades, their chemistry has been expanded to include the less stable +2 oxidation state across the entire series except promethium (Pm), facilitated by the coordination of ligands such as trimethylsilylcyclopentadienyl, C 5 H 4 SiMe 3 (Cp'). The [LnCp' 3 ] complexes have been the workhorse for the synthesis and theoretical study of the fundamental aspects of divalent lanthanide chemistry, where experimental and computational evidence have suggested the existence of different ground state (GS) configurations, 4f n+1 or 4f n 5d 1 , depending on the specific metal. Standard reduction potentials and 4f n+1 to 4f n 5d 1 promotion energies have been two factors usually considered to rationalize the occurrence of these variable GS configurations, however the driving force behind this phenomenon is still not clear. In this work we present a comprehensive theoretical approach to shed light on this matter using the [LnCp 3 ] - model systems. We begin by calculating 4f n+1 to 4f n 5d 1 promotion energies and successfully correlate them with existing experimental data. Furthermore, we analyze how changes in the GS charge distribution between the Ln ions, LnCp 3 and the reduced [LnCp 3 ] - complexes (Ln = La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu) correlate with experimental trends in redox potentials and the calculated promotion energies. For this purpose, a comprehensive theoretical work that includes relativistic ligand field density functional theory (LFDFT) and relativistic ab initio wavefunction methods was performed. This study will help the rational design of suitable environments to tune the different GS configurations as well as modulating the spectroscopic properties of new Ln 2+ complexes.

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

Probing Multiconfigurational States by Spectroscopy: The Cerium XAS L 3 ‐edge Puzzle

Abstract The Ce L 3 edge XAS spectra of CeO 2 and cerocene [Ce(C 8 H 8 ) 2 ] were calculated with relativistic ab‐initio multireference wavefunction approaches capable of reproducing the observed spectra accurately. The study aims to resolve the decades‐long puzzle regarding the relationship between the number and relative intensities of the XAS peaks and the 4f electron occupation in the ground state (GS) versus the core‐excited states (ESs). CeO 2 and cerocene exemplify the different roles of covalent bonding and wavefunction configurational composition in the observed intensity patterns. Good agreement is found between the calculated GS 4f‐shell occupations and the value derived from XAS measurements using peak areas ( n f ). The identity of the two‐peaked Ce L 3 edge is fully rationalized from the perspective of the relaxed wavefunctions for the GS and core ESs. The states underlying the different peaks differ from each other in a surprisingly simple way that can be associated with 4f 1 vs. 4f 0 sub‐configurations. Furthermore, part of one of the cerocene spectral peaks is associated with 4f 2 sub‐configurations. The pattern therefore reveals excited states that can be interpreted in terms of Ce IV and III oxidation numbers, as long assumed, with Ce II states additionally appearing in the cerocene spectrum. While this work demonstrates the rough accuracy of the conventional approach to determining n f from Ce L 3 ‐edge XAS, limitations are highlighted in terms of the ultimate accuracy of this approach and the potential of observing new types of excited states. The need to determine the sources of n f by calculations, is stressed.

Sergentu, Dumitru‐Claudiu↗
DOE OSTI2021

Electron–Nucleus Hyperfine Coupling Calculated from Restricted Active Space Wavefunctions and an Exact Two-Component Hamiltonian

Exact two-component (X2C) relativistic nuclear hyperfine magnetic field operators were incorporated in X2C ab-initio wavefunction calculations at the multi-reference restricted active space (RAS) level for calculations of nuclear hyperfine magnetic properties. Spin-orbit coupling was treated via RAS state interaction (SO-RASSI). The method was tested by calculations of electron – nucleus hyperfine coupling constants. The approach, implemented in the OpenMolcas program, overcomes a major limitation of a previous SO-RASSI implementation for hyperfine coupling that relied on non-relativistic hyperfine operators [J. Chem. Theor. Comput. 2015, 11, 538–549] and therefore had only limited applicability. Furthermore, results from calculations on systems with light and heavy main group elements, transition metals, lanthanides, and one actinide complex, demonstrate reasonably good agreement with experimental data, where available, as long as the active space can generate sufficient spin polarization.

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

Covalency of Trivalent Actinide Ions with Different Donor Ligands: Do Density Functional and Multiconfigurational Wavefunction Calculations Corroborate the Observed “Breaks”?

A comprehensive ab initio study of periodic actinide–ligand bonding trends for trivalent actinides is performed. Relativistic density functional theory (DFT) and complete active-space (CAS) self-consistent field wavefunction calculations are used to dissect the chemical bonding in the [AnCl 6 ] 3– , [An(CN) 6 ] 3– , [An(NCS) 6 ] 3– , [An(S 2 PMe 2 ) 3 ], [An(DPA) 3 ] 3– , and [An(HOPO)] – series of actinide (An = U–Es) complexes. Except for some differences for the early actinide complexes with DPA, bond orders and excess 5f-shell populations from donation bonding show qualitatively similar trends in 5f n active-space CAS vs DFT calculations. The influence of spin–orbit coupling on donation bonding is small for the tested systems. Along the actinide series, chemically soft vs chemically harder ligands exhibit clear differences in bonding trends. There are pronounced changes in the 5f populations when moving from Pu to Am or Cm, which correlate with previously noted “breaks” in chemical trends. As a result, bonding involving 5f becomes very weak beyond Cm/Bk. We propose that Cm(III) is a borderline case among the trivalent actinides that can be meaningfully considered to be involved in ground-state 5f covalent bonding.

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

Ab Initio Study of Vibronic and Magnetic 5f-to-5f and Dipole-Allowed 5f-to-6d and Charge-Transfer Transitions in [UX 6 ] n – (X = Cl, Br; n = 1, 2)

The absorption spectra of the octahedral [UX 6 ] n– (X = Cl, Br; n = 1, 2) complexes in the near-infrared (NIR) and UV–vis spectral regions were studied theoretically, using a relativistic restricted active space second-order perturbation theory (RASPT2) wavefunction framework, with the spin–orbit (SO) coupling treated by state interaction, in conjunction with Kohn–Sham density functional theory calculations for determining the vibrational normal modes. The electric-dipole-allowed ligand-to-metal charge-transfer (LMCT) and 5f-to-6d transitions, and the electric-dipole-forbidden 5f-to-5f ligand field (LF) transitions, are thereby obtained within the same theoretical framework. For the 5f-to-5f LF transitions, the observed absorption intensity is mostly due to vibronic coupling with low-energy electric-dipole-allowed transitions, but in some cases, the magnetic dipole intensity of the purely electronic transition has comparable intensity to the vibronic transitions. Here, experimental LF spectra of 5f 2 open-shell systems have been reported decades back, but ab initio calculations of their vibronic intensity have not yet been reported in the literature. Although the LF spectra for the 5f 2 systems can be assigned in detail, based on the calculations, the spectra are very complex and the underlying electronic states are strongly multiconfigurational. Therefore, the usefulness of the LF spectra beyond serving as a “fingerprint” of the LF and the metal oxidation state appears to be limited.

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

Finite nucleus effects on relativistic energy corrections

The effect of using a finite nucleus model in quantum-chemical calculations is examined. Relativistic corrections from the first order Foldy-Wouthuysen terms are affected indirectly by the change in wavefunction, but also directly as a result of revised expressions for the Darwin and spin-orbit terms due to the change in nuclear potential. A calculation for the Rn atom indicates that the mass-velocity and Darwin corrections are much more sensitive to the finite nucleus than the non-relativistic total energy, but that the total contribution for these two terms is quite stable provided the revised form of the Darwin term is used. The spin-orbit interaction is not greatly affected by the choice of nuclear model.

Dyall, Kenneth G.↗
DOE OSTI2022

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). Here, in this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H 2 Z 2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

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

Shake-up and shake-off spectra in the electron capture decay of atomic $^7$Be

The most stringent laboratory-based experimental limits on the existence of sub-MeV sterile neutrinos are currently set by decay spectroscopy of radioactive $^7$Be embedded into superconducting sensors. The systematic uncertainties are dominated by the modeling of the electron shake-up and shake-off spectra that are not based on state-of-the-art atomic theory and do not include electron correlations or relativistic effects. We have used the multiconfiguration Dirac-Fock formalism to obtain correlated wavefunctions ab initio and compute all single and double shake processes in the electron capture decay of atomic $^7$Be. The simulations can explain some but not all of the observed spectral features, likely because the wave functions are modified by the Ta sensor material that the $^7$Be is embedded into. The new models also show that the L/K electron capture ratio of $^7$Be in Ta has previously been slightly underestimated revising the previous value of 0.070(7) to a new value of 0.0756(20).

Atomic Physics (physics.atom-ph)↗
NASA NTRS1984

Theoretical electric quadrupole transition probabilities for Ca, Sr and Ba

The 1D-1S quadrupole transition probabilities for Ca, Sr and Ba have been computed using extended GTO and STO valence basis sets and configuration-interaction wavefunctions that include the important core-valence correlation effects. For Ba and Sr, the relativistic contraction of the core orbitals was accounted for in the GTO calculations by a relativistic effective-core potential. The computed Einstein coefficient for Ca of 39.6/s is in excellent agreement with the recent experimental value of 40 + or - 8/s. The best Einstein coefficients for Sr (44.7/s) and Ba (2.98/s) imply increasing quadrupole line strengths down the column. Relativistic effects substantially increase the quadrupole Einstein coefficient for Ba.

Bauschlicher, C. W., Jr.↗
NASA NTRS1978

Relativistic scattered-wave theory. II - Normalization and symmetrization

Formalisms for normalization and symmetrization of one-electron Dirac scattered-wave wavefunctions are presented. The normalization integral consists of one-dimensional radial integrals for the spherical regions and an analytic expression for the intersphere region. Symmetrization drastically reduces the size of the secular matrix to be solved. Examples for planar Pb2Se2 and tetrahedral Pd4 are discussed.

Yang, C. Y.↗
DOE OSTI2025

Numerically exact configuration interaction at quadrillion-determinant scale

The combinatorial growth of configuration interaction (CI) has long limited this formally exact quantum chemistry method to only the smallest molecules. Here, we report a numerically exact CI calculation exceeding one quadrillion (10 15 ) determinants, made possible by a lossless categorical compression strategy within the small-tensor-product distributed active space (STP-DAS) framework. This approach overcomes the traditional memory bottlenecks of CI by a numerically exact compression of the wavefunction representation and reformulating the most computationally demanding matrix–vector operations. Using this method, we performed a fully relativistic CI calculation of the ground state of HBrTe with over 10 15 complex-valued determinants in just 34.5 h on 1000 computing nodes—the largest CI calculation ever reported. We further achieved fast computation for systems with hundreds of billions of determinants on only a few compute nodes. Extensive benchmarks confirm that the method retains full numerical exactness while cutting memory and computational cost by orders of magnitude. Compared to previous state-of-the-art CI calculations, this work achieves a 1000 times increase in CI space, a 10 6 -fold increase in floating-point operations performed, and a 10 6 -fold improvement in computational speed.

Computational chemistry↗
NASA NTRS2004

Accurate Calculation of Oscillator Strengths for CI II Lines Using Non-orthogonal Wavefunctions

Non-orthogonal orbitals technique in the multiconfiguration Hartree-Fock approach is used to calculate oscillator strengths and transition probabilities for allowed and intercombination lines in Cl II. The relativistic corrections are included through the Breit-Pauli Hamiltonian. The Cl II wave functions show strong term dependence. The non-orthogonal orbitals are used to describe the term dependence of radial functions. Large sets of spectroscopic and correlation functions are chosen to describe adequately strong interactions in the 3s(sup 2)3p(sup 3)nl (sup 3)Po, (sup 1)Po and (sup 3)Do Rydberg series and to properly account for the important correlation and relaxation effects. The length and velocity forms of oscillator strength show good agreement for most transitions. The calculated radiative lifetime for the 3s3p(sup 5) (sup 3)Po state is in good agreement with experiment.

Tayal, S. S.↗
DOE OSTI2021

Theoretical Prediction and Interpretation of 237 Np Mössbauer Isomer Shifts

Different theoretical approaches for the calculation of 237 Np Mössbauer isomer shifts are investigated. The traditional contact density route is compared to a previously proposed alternative approach that uses energy derivatives with respect to the nuclear radius. Both approaches yield similar results as long as suitable basis sets augmented with large exponents and relativistic methods are used. Density functional theory (DFT) calculations do not show a strong dependency of the 237 Np isomer shift on the chosen functional. Wavefunction calculations show that dynamic electron correlation can be important when covalent bonding influences the isomer shift. Effects from spin–orbit coupling are small. The isomer shifts of ionic solids and Np(III) organometallic complexes are largely governed by the oxidation state of Np. Isomer shifts of organometallic Np(IV) complexes are strongly affected by donation bonding. Furthermore, detailed analysis of the wavefunction results with different active spaces demonstrates that correlation among the outer core Np and occupied ligand frontier orbitals contributes significantly to isomer shifts of Np(IV) compounds.

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

Simulating Effective QED on Quantum Computers

In recent years simulations of chemistry and condensed materials has emerged as one of the preeminent applications of quantum computing, offering an exponential speedup for the solution of the electronic structure for certain strongly correlated electronic systems. To date, most treatments have ignored the question of whether relativistic effects, which are described most generally by quantum electrodynamics (QED), can also be simulated on a quantum computer in polynomial time. Here we show that effective QED, which is equivalent to QED to second order in perturbation theory, can be simulated in polynomial time under reasonable assumptions while properly treating all four components of the wavefunction of the fermionic field. In particular, we provide a detailed analysis of such simulations in position and momentum basis using Trotter-Suzuki formulas. We find that the number of T -gates needed to perform such simulations on a 3 D lattice of n s sites scales at worst as O ( n s 3 / ϵ ) 1 + o ( 1 ) in the thermodynamic limit for position basis simulations and O ( n s 4 + 2 / 3 / ϵ ) 1 + o ( 1 ) in momentum basis. We also find that qubitization scales slightly better with a worst case scaling of O ~ ( n s 2 + 2 / 3 / ϵ ) for lattice eQED and complications in the prepare circuit leads to a slightly worse scaling in momentum basis of O ~ ( n s 5 + 2 / 3 / ϵ ) . We further provide concrete gate counts for simulating a relativistic version of the uniform electron gas that show challenging problems can be simulated using fewer than 10 13 non-Clifford operations and also provide a detailed discussion of how to prepare multi-reference configuration interaction states in effective QED which can provide a reasonable initial guess for the ground state. Finally, we estimate the planewave cutoffs needed to accurately simulate heavy elements such as gold.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗