Spin-Orbit Quantum Monte Carlo for heavy element Equations of state
Explore the source record for details and available documents.
Clay III, Raymond C. [Sandia National Laboratories↗
SEARCH · Engineering Papers
Results for “ORBIT EQUATION”
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.
One-Body Properties and Their Perturbative Accuracy with Aufbau Suppressed Coupled Cluster Theory
In this work we derived and implemented the calculation of the one-body reduced density matrix for Aufbau suppressed coupled cluster theory, from which excited state natural orbitals and one-body properties, like atomic populations and dipole moments, are obtained. We utilized the natural orbitals to refine the ASCC solution for simple valence and Rydberg systems, exploring the process of repeatedly solving the ASCC equations in successive natural orbital bases to achieve independence from the starting molecular orbitals. For dipole moments in small molecules where high-level comparison data is available, we find that the accuracy of ASCC essentially matches that of linear response and equation-of-motion coupled cluster as long as care is taken to preserve the response's perturbative completeness.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Linear-Scaling Local Natural Orbital-Based Full Triples Treatment in Coupled-Cluster Theory
We present an efficient, asymptotically linear-scaling implementation of the canonically O(N 8 ) coupled-cluster method with singles, doubles, and full triples excitations (CCSDT) method. We apply the domain-based local pair natural orbital (DLPNO) approach for computing CCSDT amplitudes. Our method, called DLPNO–CCSDT, uses the converged coupled-cluster amplitudes from a preceding DLPNO–CCSD(T) computation as a starting point for the solution of the CCSDT equations in the local natural orbital basis. To simplify the working equations, we t1-dress our two-electron integrals and Fock matrices, allowing our equations to take on the form of CCDT. With appropriate parameters, our method can recover more than 99.99% of the total canonical CCSDT correlation energy. In addition, we demonstrate that our method consistently yields sub-kJ mol –1 errors in relative energies when compared to canonical CCSDT, and, likewise, when computing the difference between CCSDT and CCSD(T). Finally, to highlight the low scaling of our algorithm, we present timings on linear alkanes (up to 30 carbons and 730 basis functions) and water clusters (up to 131 water molecules and 3144 basis functions).
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Impact of solvation on the electronic resonances in uracil
Interactions of low-energy electrons with the DNA and RNA nucleobases are known to form metastable states, known as electronic resonances. In this work, we study electron attachment to solvated uracil, an RNA nucleobase, using the orbital stabilization method at the Equation of Motion-Coupled Cluster for Electron Affinities with Singles and Doubles (EOM-EA-CCSD) level of theory with the Effective Fragment Potential (EFP) solvation method. We benchmarked the approach using multireference methods, as well as by comparing EFP and full quantum calculations. The impact of solvation on the first one particle (1p) shape resonance, formed by electron attachment to the π* LUMO orbital, as well as the first two particle one hole (2p1h) resonance, formed by electron attachment to neutral uracil's π–π* excited state, was investigated. We used molecular dynamics simulations for solvent configurations and applied charge stabilization technique-based biased sampling to procure configurations adequate to cover the entire range of the electron attachment energy distribution. The electron attachment energy in solution is found to be distributed over a wide range of energies, between 4.6 eV to 6.8 eV for the 2p1h resonance, and between −0.1 eV to 2 eV for the 1p resonance. The solvent effects were similar for the two resonances, indicating that the exact electron density of the state is not as important as the solvent configurations. Multireference calculations extended the findings showing that solvation effects are similar for the lowest four resonances, further indicating that the specific solute electron density is not as important, but rather the water configurations play the most important role in solvation effects. Lastly, by comparing bulk solvation to clusters of uracil with a few water molecules around it, we find that the impact of microsolvation is very different from that of bulk solvation.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Transverse orbital angular momentum in the proton
Using the equations of motion of QCD and Lorentz invariance relations, we show a new, more direct way of obtaining the decomposition of the proton's quark transverse angular momentum into its spin and orbital components. The new decomposition can be understood as a twist-three relation, with the quark transverse spin described by the parton distribution g T . We find that the orbital component can be expressed in terms of a moment in the quark transverse momentum, k T , of the generalized transverse momentum-dependent distribution F 12 , or alternatively in terms of the twist-three generalized parton distributions H 2T and $\tilde{E}$ 2T .
Lorentz invariance relations↗
Equation-of-motion internally contracted multireference unitary coupled-cluster theory
The accurate computation of excited states remains a challenge in electronic structure theory, especially for systems with a ground state that requires a multireference treatment. In this work, we introduce a novel equation-of-motion (EOM) extension of the internally contracted multireference unitary coupled-cluster framework (ic-MRUCC), termed EOM-ic-MRUCC. EOM-ic-MRUCC follows the transform-then-diagonalize approach, in analogy to its non-unitary counterpart. By employing a projective approach to optimize the ground state, the method retains additive separability and proper scaling with system size. We show that excitation energies are size-intensive if the EOM operator satisfies the “killer” and the projective conditions. Furthermore, we propose to represent changes in the reference state upon electron excitation via projected many-body operators that span the active orbitals and show that the EOM equations formulated in this way are invariant with respect to active orbital rotations. We test the EOM-ic-MRUCC method truncated to single and double excitations by computing the potential energy curves for several excited states of a BeH2 model system, the HF molecule, and water undergoing symmetric dissociation. Across these systems, our method delivers accurate excitation energies and potential energy curves within 5 mE h (∼0.14 eV) from full configuration interaction. Here, we find that truncating the Baker–Campbell–Hausdorff series to fourfold commutators contributes negligible errors (on the order of 10 −5 E h or less), offering a practical route to highly accurate excited-state calculations with reduced computational overhead.
74 ATOMIC AND MOLECULAR PHYSICS↗
Linear-Scaling Asymmetric Triples Correction through the Solution of the DLPNO–CCSD Lambda Equations: DLPNO–CCSD(T) Λ
In this research, we derive equations for solving for the stationary points of the DLPNO–CCSD Lagrangian, in the t 1 -transformed formalism introduced earlier and as currently implemented in the P SI 4 quantum chemistry software package. These lambda equations in the local pair natural orbital basis allow for the evaluation of CCSD(T) Λ energetics with linear-scaling computational effort, also known as the asymmetric triples correction. This DLPNO–CCSD(T) Λ method allows for accurate triples contributions to be computed for larger molecules, especially in cases that CCSD(T) is known to be insufficient, such as with multireference systems and bond-breaking systems. We showcase the accuracy of our code on reaction energies, barrier heights, and noncovalent interaction energies. Also showcased are the capabilities of our code by evaluating DLPNO–CCSD(T) Λ energetics on large noncovalent dimers up to 112 atoms, as well as a rhodium catalyst complex containing 66 atoms.
Cluster chemistry↗
Understanding plasma turbulence through exact coherent structures
Plasma turbulence is a key challenge in understanding transport phenomena in magnetically confined plasmas. This work presents a generalized framework to analyze plasma turbulence that utilizes periodic orbit theory. In periodic orbit theory, doubly periodic solutions (coherent structures) of the governing equation(s) serve as building blocks of the considered turbulent dynamics. To illustrate the concept and method, the particularly simple Kuramoto–Sivashinsky (referred to here as LMRT for the original authors: LaQuey, Mahajan, Rutherford, and Tang) trapped-ion mode toy model is used. By applying numerical optimization techniques to the LMRT equation, we extract coherent spacetime patterns that represent the library of allowable fundamental structures of the equation. These structures provide a framework to systematically describe turbulence as a composition of recurrent solutions, revealing an underlying order within chaotic plasma motion. Although illustrated here using the simplified LMRT model for clarity, this framework provides a general strategy that can be extended to more complex and realistic models of plasma turbulence, including gyrokinetic systems. This offers a new method for predicting and potentially controlling transport processes in fusion plasmas by providing a bridge between nonlinear dynamical systems theory and plasma physics in the form of a generalized framework with which to analyze and understand spatially extended nonlinear partial differential equations.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Quasiparticle band structure and excitonic optical response in V 2 O 5 bulk and monolayer
Not Available
Garcia, Claudio (ORCID:0009000077503536)↗
Quantum stress and torsion distributions in the deuteron
Stress distributions in the deuteron are related to form factors of the asymmetric energy-momentum tensor through three-dimensional Fourier transforms. There are eleven such form factors, which we calculate in an impulse approximation. We compare the obtained form factors to prior results for the six form factors that have been previously calculated. We then elaborate on the formalism for relating the form factors to internal distributions of mass, mass flux, momentum, stresses, and forces, and obtain results for all of these distributions. We obtain the principal stresses for the symmetric part of the stress tensor, and show that the antisymmetric part describes the reorientation of fermion spin by torsion stress when the nucleon moves between the S and D waves. Force distributions in the nucleons depend on the so-called nonconserved form factors through the Cauchy momentum equation, and are nonradial owing to the presence of tensor forces and spin-orbit coupling.
QCD phenomenology↗
Formation and Redshift Evolution of Dark Matter Spikes
Dark matter density spikes forming around adiabatically growing black holes can dramatically enhance indirect and direct detection signals. Canonical predictions, however, assume a zero-mass seed in a purely dark matter environment and do not track the long-term dynamical impact of surrounding stars. We present a semi-analytic framework that first generalizes adiabatic spike formation to include finite seed masses, stellar cusps, and non-circular orbits, and then studies the subsequent cosmic evolution by solving coupled Fokker-Planck equations for the dark matter and stellar phase-space distributions, with a heating rate modulated by the cosmic star formation rate. Starting conservatively from canonical Gondolo-Silk spikes and marginalizing over astrophysical uncertainties, we find that stellar gravitational heating drives the inner slope towards $γ_χ\simeq 1.5$ within a few Gyrs (e.g by $z \lesssim 2$ for spikes formed at $z\simeq 10$), yielding overdensities two to four orders of magnitude below canonical expectations but still well above an NFW-like cusp. We provide redshift-dependent benchmarks for the column density and $J$-factor relevant to scattering, decay and annihilation signatures. Any robust interpretation of indirect dark matter signals from galactic nuclei must account for this evolution.
Herrera, Gonzalo [MIT, LNS; MIT, MKI; Harvard U.]↗
Time inversion symmetry in the Dirac and Schrödinger-Pauli theories
The Schrödinger-Pauli theory is generally believed to give a faithful representation of the nonrelativistic and weakly relativistic limit of the Dirac theory. However, the Schrödinger-Pauli theory is fundamentally incomplete in its account of broken time inversion symmetry, e.g., in magnetically ordered systems. Here, in the Dirac theory of the electron, magnetic order breaks time inversion symmetry even in the nonrelativistic limit, whereas time inversion symmetry is effectively preserved in the Schrödinger-Pauli theory in the absence of spin-orbit coupling. In the Dirac theory, the Berry curvature $1/(2m^2 c^2)$ is thus an intrinsic property of nonrelativistic electrons similar to the well-known spin magnetic moment $e\hbar/(2m)$, while this result is missed by the nonrelativistic or weakly relativistic Schrödinger-Pauli equation. In ferromagnetically ordered systems, the intrinsic Berry curvature yields a contribution to the anomalous Hall conductivity independent of spin-orbit coupling.
Winkler, R. [Northern Illinois Univ., DeKalb, IL (↗
Local Pair Natural Orbital-Based Coupled-Cluster Theory through Full Quadruples (DLPNO–CCSDTQ)
In this work, we implement a local pair natural orbitalbased coupled-cluster method through the full treatment of quadruple excitations (CCSDTQ). The domain-based local pair natural orbital (DLPNO) approach, which has successfully been applied to lower levels of coupled-cluster theory, is utilized in our algorithm, and thus our algorithm is called DLPNO-CCSDTQ. For simplicity in the working equations and in the implementation, we t 1 -dress the twoelectron integrals as well as Fock matrix elements. Our method can recover CCSDTQ-CCSDT and CCSDTQ-CCSDT(Q) energy differences on the order of 0.01−0.05 kcal mol −1 , even at a loose quadruples natural orbital (QNO) occupation number cutoff of 3.33 × 10 −6 . To highlight the capabilities of our code and its potential future applications, we showcase computations that would be intractable with canonical CCSDTQ, such as the benzene dimer, (H 2 O) 17 , and adamantane. With sufficient computing resources, computations up to 15 heavy atoms (40 atoms overall) may be feasible for fully bonded 3D systems.
Cluster chemistry↗
Direct and indirect spin current generation and spin-orbit torques in ferromagnet/nonmagnet/ferromagnet trilayers
Spin-orbit torques in ferromagnet/nonmagnet/ferromagnet trilayers are studied using a combination of symmetry analysis, circuit theory, semiclassical simulations, and first-principles calculations using the nonequilibrium Green's function method with supercell disorder averaging. Here, we focus on unconventional processes involving the interplay between the two ferromagnetic layers, which are classified into direct and indirect mechanisms. The direct mechanism involves spin current generation by one ferromagnetic layer and its subsequent absorption by the other. In the indirect mechanism, the in-plane spin-polarized current from one ferromagnetic layer “leaks” into the other layer, where it is converted into an out-of-plane spin current and reabsorbed by the original layer. The direct mechanism results in a predominantly dampinglike torque, which damps the magnetization towards a certain direction 𝐬 𝑑 . The indirect mechanism results in a predominantly fieldlike torque with respect to a generally different direction 𝐬 𝑓 . Similarly to the current-in-plane giant magnetoresistance, the indirect mechanism is only active if the thickness of the nonmagnetic spacer is smaller than or comparable to the mean free path. Numerical calculations for a semiclassical model based on the Boltzmann equation confirm the presence of both direct and indirect mechanisms of spin current generation. First-principles calculations reveal sizable unconventional spin-orbit torques in Co/Cu/Co, Py/Cu/Py, and Co/Pt/Co trilayers and provide strong evidence of indirect spin current generation.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Two-component atomic Fermi superfluid with spin-orbital coupling in thin-spherical-shell geometry
We present a theory of two-component atomic Fermi superfluid with tunable pairing interaction in a thin spherical shell subject to spin-orbit coupling (SOC). By incorporating SOC into the Fermi superfluid in the BCS-Bose Einstein condensation (BEC) crossover, we obtain the energy spectrum and equations of state. While the order parameter and chemical potential are suppressed by SOC on the BCS side, the former remains positive but the latter may be pushed to negative values by SOC.Meanwhile, the compressibility exhibits kinks as the pairing interaction or SOC varies, indicating singular behavior of higher derivatives of the free energy despite the presence of the superfluid order parameter. The minimum of the energy dispersion indicates a decoupling of the energy gap from the order parameter, and the compressibility kinks occur when the energy gap approaches 0. We found the particle-hole mixing in the Fermi superfluid essential for the compressibility kinks since a Fermi gas with the same SOC but no pairing interaction only exhibits complicated dispersions but no singular behavior. Therefore, our results show that a combination of superfluid, SOC, and geometry can lead to interesting phenomena. Here, we also discuss possible experimental realizations and implications.
Fermi gases↗
EMPDF : inferring the Milky Way mass with data-driven distribution function in phase space
We introduce the emPDF (empirical distribution function), a novel dynamical modelling method that infers the gravitational potential from kinematic tracers with optimal statistical efficiency under the minimal assumption of steady state. emPDF determines the best-fitting potential by maximizing the similarity between instantaneous kinematics and the time-averaged phase-space distribution function (DF), which is empirically constructed from observation upon the theoretical foundation of oPDF (Han et al. 2016). This approach eliminates the need for presumed functional forms of DFs or orbit libraries required by conventional DF- or orbit-based methods. emPDF stands out for its flexibility, efficiency, and capability in handling observational effects, making it preferable to the popular Jeans equation or other minimal assumption methods, especially for the Milky Way (MW) outer halo where tracers often have limited sample size and poor data quality. We apply emPDF to infer the MW mass profile using Gaia DR3 data of satellite galaxies and globular clusters, obtaining enclosed masses of M (,r) = 26±8, 46±8, 90±13, and 149±40 x 10 10 M ⊙ at r = 30, 50, 100, and 200 kpc, respectively. These are consistent with the updated constraints from simulation-informed DF fitting (Li et al. 2020). While the simulation-informed DF offers superior precision owing to the additional information extracted from simulations, emPDF is independent of such supplementary knowledge and applicable to general tracer populations. emPDF is currently implemented for tracers with complete 6D kinematics within spherical potentials, but it can potentially be extended to address more general problems.
Astrophysics of Galaxies (astro-ph.GA)↗
An–imidophosphorane (An = U–Pu) bond covalency and proton-coupled electron transfer thermodynamics driven by orbital energy matching
A series of mid-actinide (An = U–Pu) tetrahomoleptic complexes supported by highly electron-donating imidophosphorane ligands, NPC ([NP t Bu(pyrr) 2 ] − , where t Bu = C(CH 3 ) 3 ; pyrr = pyrrolidinyl = N(C 4 H 8 )), are systematically investigated computationally and experimentally to elucidate the nature of actinide–ligand (An–L) covalency across the An 3+/4+/5+ oxidation states. Trends in An–L bonding and redox properties for these complexes, together with their protonated counterparts, are examined using orbital-, electron density-, and energy-decomposition-based methods. This integrated approach reveals progressively improved energy matching between α-spin An 5f and N im 2p orbitals with increasing atomic number and oxidation state, becoming particularly pronounced in the ligand-dominant π-bonding orbitals of An 4+ and An 5+ . In contrast to the An 3+ species, the enhanced An 5f π contributions in the higher-valent counterparts drive the increase in An–N im covalency for later An, thereby inverting the covalency trend to U < Np < Pu. Redistribution of electron density towards the An and N im atomic basins due to the growing energy-matching assisted covalency correlates with higher pKa values and increased N im –H bond dissociation free energies in protonated An 4+ complexes. Electron density at Nim in An 4+ shows a linear correlation with the p K a values calculated via the Bordwell equation. Calculations predict a cathodic shift of 0.84–1.00 V in the redox couples upon protonation, a trend validated when experimentally accessible. These findings demonstrate an increasing role of covalency driven by orbital energy matching from U to Pu in tuning the thermodynamic driving force for proton-coupled electron transfer in the An 5+ species.
38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
Thermal quasiparticle theory
The widely used thermal Hartree–Fock (HF) theory is generalized to include the effect of electron correlation while maintaining its quasi-independent-particle framework. An electron-correlated internal energy (or grand potential) is postulated in consultation with the second-order finite-temperature many-body perturbation theory (MBPT), which then dictates the corresponding thermal orbital (quasiparticle) energies in such a way that all fundamental thermodynamic relations are obeyed. The associated density matrix is of a one-electron type, whose diagonal elements take the form of the Fermi–Dirac distribution functions, when the grand potential is minimized. The formulas for the entropy and chemical potential are unchanged from those of Fermi–Dirac or thermal HF theory. The theory thus stipulates a finite-temperature extension of the second-order Dyson self-energy of one-particle many-body Green’s function theory and can be viewed as a second-order, diagonal, frequency-independent, thermal inverse Dyson equation. At low temperatures, the theory approaches finite-temperature MBPT of the same order, but it may outperform the latter at intermediate temperatures by including additional electron-correlation effects through orbital energies. Here, a physical meaning of these thermal orbital energies is proposed (encompassing that of thermal HF orbital energies, which has been elusive) as a finite-temperature version of Janak’s theorem.
74 ATOMIC AND MOLECULAR PHYSICS↗