Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Wave Function”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Effective field theory analysis of the Coulomb breakup of the one-neutron halo nucleus $$^{19}$$C

Abstract We analyse the Coulomb breakup of $$^{19}$$ 19 C measured at 67 A MeV at RIKEN. We use the Coulomb-Corrected Eikonal (CCE) approximation to model the reaction and describe the one-neutron halo nucleus $$^{19}$$ 19 C within Halo Effective Field Theory (Halo EFT). At leading order we obtain a fair reproduction of the measured cross section as a function of energy and angle. The description is insensitive to the choice of optical potential, as long as it accurately represents the size of $${}^{18}$$ 18 C. It is also insensitive to the interior of the $${}^{19}$$ 19 C wave function. Comparison between theory and experiment thus enables us to infer asymptotic properties of the ground state of $$^{19}$$ 19 C: these data put constraints on the one-neutron separation energy of this nucleus and, for a given binding energy, can be used to extract an asymptotic normalisation coefficient (ANC). These results are confirmed by CCE calculations employing next-to-leading order Halo EFT descriptions of $${}^{19}$$ 19 C: at this order the results for the Coulomb breakup cross section are completely insensitive to the choice of the regulator. Accordingly, this reaction can be used to constrain the one-neutron separation energy and ANC of $${}^{19}$$ 19 C.

Physics↗
DOE OSTI2024

Anisotropic excitonic magnetism from discrete C 4 symmetry in CeRhIn 5

Anisotropy in strongly correlated materials is a central parameter in determining the electronic ground state and is tuned through the local crystalline electric field. This is notably the case in the CeCo x Rh 1 − x In 5 system where the ground-state wave function can provide the basis for antiferromagnetism and/or unconventional superconductivity. We develop a methodology to understand the local magnetic anisotropy and experimentally investigate with neutron spectroscopy applied to antiferromagnetic ( T N = 3.8 K ) CeRhIn 5 , which is isostructural to d -wave superconducting ( T c = 2.3 K ) CeCoIn 5 . Through diagonalizing the local crystal field Hamiltonian with discrete tetragonal C 4 point group symmetry and coupling these states with the random phase approximation, we find two distinct modes polarized along the crystallographic c and a − b planes, agreeing with experiment. The anisotropy and bandwidth, underlying the energy scale of these modes, are tuneable with a magnetic field which we use experimentally to separate in energy single and multiparticle excitations thereby demonstrating the instability of excitations polarized within the crystallographic a − b plane in CeRhIn 5 . We compare this approach to a S eff = 1 2 parametrizations and argue for the need to extend conventional SU(2) theories of magnetic excitations to utilize the multilevel nature of the underlying crystal-field basis states constrained by the local point-group C 4 symmetry. Published by the American Physical Society 2024

Brener, D. J. (ORCID:0000000333400177)↗
DOE OSTI2025

Hydrogen Bond Benchmark: Focal‐Point Analysis and Assessment of DFT Functionals

We performed a hierarchical, convergent ab initio benchmark study and systematically analyzed the performance of density functional approximations for describing hydrogen bonds in small neutral, cationic, and anionic complexes, as well as in larger systems involving amide, urea, deltamide, and squaramide moieties. Focal point analyses (FPA), extrapolating to the ab initio limit, were carried out using correlated wave function methods up to CCSDT(Q) for the small complexes and CCSD(T) for the larger systems, together with correlation-consistent Gaussian basis sets up to the complete basis set limit. Optimized geometries and vibrational frequencies were obtained at the CCSD(T) level. The resulting FPA hydrogen-bond energies converge within a few tenths of a kcal mol −1 . These reference data were used to evaluate 60 density functionals (including 12 dispersion-corrected), spanning the local-density approximation (LDA), generalized gradient approximations (GGAs), meta-GGAs, hybrids, meta-hybrids, double-hybrids, and range-separated hybrids. Overall, the meta-hybrid M06-2X provides the best performance for both hydrogen bond energies and geometries, while the dispersion-corrected GGAs BLYP-D3(BJ) and BLYP-D4 also yield accurate hydrogen-bond data and can serve as cost-effective options for studying large and complex systems.

coupled cluster theory↗
DOE OSTI2022

Multireference Methods are Realistic and Useful Tools for Modeling Catalysis

Abstract Highly correlated systems, in particular those that include transition metals, are ubiquitous in catalysis. The significant static correlation found in such systems is often poorly accounted for using Kohn Sham density functional theory methods, as they are single determinantal in nature. Applications to catalysis of more rigorous and appropriate multiconfigurational methods have been reported in select instances, but their use remains rare. We discuss obstacles that hinder the routine application of multireference (MR) wave function theoretical calculations to catalytic systems and the current state of the art with respect to removing those obstacles.

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

Weighted nodal domain averages of eigenstates for quantum Monte Carlo and beyond

In this report we study the nodal properties of many-body eigenstates of stationary Schrödinger equation that affect the accuracy of real-space quantum Monte Carlo calculations. In particular, we introduce weighted nodal domain averages that provide a new probe of nodal surfaces beyond the usual expectations. Particular choices for the weight function reveal, for example, that the difference between two arbitrary fermionic eigenvalues is given by the nodal hypersurface integrals normalized by overlaps with the bosonic ground state of the given Hamiltonian. Noninteracting and fully interacting Be atom with corresponding almost exact and approximate wave functions are used to illustrate several aspects of these concepts. Variational formulations that employ different weights are proposed for prospective improvement of nodes in variational and fixed-node diffusion Monte Carlo calculations.

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

Shadow Lagrangian dynamics for superfluidity

Motivated by a similar approach for Born-Oppenheimer molecular dynamics, this paper proposes an extended "shadow" Lagrangian density for quantum states of superfluids. The extended Lagrangian contains an additional field variable that is forced to follow the wave function of the quantum state through a rapidly oscillating extended harmonic oscillator. By considering the adiabatic limit for large frequencies of the harmonic oscillator, we can derive the two equations of motions, a Schrödinger-type equation for the quantum state and a wave equation for the extended field variable. The equations are coupled in a nonlinear way, but each equation individually is linear with respect to the variable that it defines. The computational advantage of this new system is that it can be easily discretized using linear time stepping methods, where we propose to use a Crank-Nicolson-type approach for the Schrödinger equation and an extended leapfrog scheme for the wave equation. Furthermore, the difference between the quantum state and the extended field variable defines a consistency error that should go to zero if the frequency tends to infinity. By coupling the time-step size in our discretization to the frequency of the harmonic oscillator we can extract an easily computable consistency error indicator that can be used to estimate the numerical error without additional costs. The findings are illustrated in numerical experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

The graviton four-point function in de Sitter space

We compute the tree-level late-time graviton four-point correlation function, and the related quartic wavefunction coefficient, for Einstein gravity in de Sitter spacetime. We derive this result in several ways: by direct calculation, using the in-in formalism and the wavefunction of the universe; by a heuristic derivation leveraging the flat space wave-function coefficient; and by using the boostless cosmological bootstrap, in particular the combination of the cosmological optical theorem, the amplitude limit, and the manifestly local test. We find agreement among the different methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

First principles study of the Fermi surface topology of CeCu 2 ⁢Si 2

Since the discovery of heavy-fermion superconductivity in CeCu 2 ⁢Si 2 , the material has attracted great interest, particularly with regard to the nature of the superconducting pairing and its mechanism. Consequently, it is essential to better understand the electronic Fermi surface topology and its role in strong antiferromagnetic fluctuations. The standard density functional theory method is insufficient to model the interplay of strong on-site Coulomb repulsion in localized 4⁢𝑓 electrons and their hybridization with itinerant ligand-orbital electrons. We have performed electronic ground-state calculations on CeCu 2 ⁢Si 2 using the Gutzwiller wave function approximation. The Gutzwiller approximation captures the quasiparticle band renormalization from the strong on-site Coulomb repulsion. We have performed an analysis of this effect on the electronic structure and the Fermi surface topology by varying the interaction strength and taking into account the crystal-field splitting. Using the de Haas-van Alphen effect, the extremal Fermi surface cross-sectional areas were calculated to quantify the effects of quasiparticle mass renormalization on the Fermi surface. Our results confirm the presence of two Fermi surface sheets corresponding to the heavy (488⁢𝑚 𝑒 ) and light (4.35⁢𝑚 𝑒 ) quasiparticles when the crystal-field splitting is accounted for on equal footing with the electronic correlations. This method gives the best agreement with experimental measurements as well as the renormalized band method.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Transverse structure of the pion beyond leading twist with basis light-front quantization

We investigate the twist-3transverse-momentum-dependent parton distribution functions (TMDs) of the pion with basis light-front quantization. The twist-3 TMDs are not independent and can be decomposed into twist-2and genuine twist-3 terms from the equations of motion (EOM). We compute the TMDs from the resulting light-front wave functions obtained by diagonalizing the light-front QCD Hamiltonian, determined for pion’s constituent quark-antiquark and quark-antiquark-gluon Fock sectors, with three-dimensional confinement. We also obtain the twist-3 parton distribution functions (PDFs) and show that they preserve the sum rule, which affirms the robustness of our approach. This is the first time that theoretical predictions are made for subleading twist structures of the pion containing interference terms between two light-front Fock sectors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Nonadiabatic dynamics of photoexcited thiopyridone isomers: An interplay between El-Sayed’s conditions and energy gap law

We investigated the non-adiabatic dynamics of photoexcited thiopyridone systems across their ortho-, meta-, and para-isomeric forms. The relaxation pathways of the three isomers in both gas phase and solvent environments are mapped using surface hopping dynamics based on time-dependent density functional theory. Our analysis highlights the influence of isomeric structures on photophysical behavior, offering insights into design principles to control photochemical phenomena. The simulations suggest a systematic reduction in the rate of intersystem crossing (ISC) from ortho- to meta- to para-isomer. Comparisons with multiconfigurational wave function methods in the gas phase further demonstrate how electronic structure influences the predicted dynamical pathways. The simulated dynamics demonstrates that the spin–orbit coupling strength alone does not determine the rate of ISC, as both state energetics and underlying electronic and structural features play decisive roles. These aspects explain the much slower ISC in the para-isomer, as well as the non-negligible role of the El-Sayed forbidden pathway in the computed ISC dynamics.

Complete-active space self-consistent field↗
DOE OSTI2025

Diverse electronic topography in a distorted kagome metal LaTi 3 ⁢Bi 4

Recent reports on a family of kagome metals of the form 𝐿⁢𝑛⁢Ti 3 ⁢Bi 4 (𝐿⁢𝑛 = Lanthanide) have stoked interest due to the combination of highly anisotropic magnetism and a rich electronic structure. The electronic structure near the Fermi level is proposed to exhibit Dirac points and van Hove singularities (VHSs). Here, in this manuscript, we use angle-resolved photoemission spectroscopy measurements in combination with density functional theory calculations to investigate the electronic structure of an interesting kagome metal LaTi 3 ⁢Bi 4 . Our results reveal multiple VHSs with one VHS located in the vicinity of the Fermi level. We clearly observe two flat bands, which originate from the destructive interference of wave functions within the Ti kagome motif. These flat bands and VHSs originate from Ti 𝑑 orbitals and are very responsive to the polarization of the incident beam. We notice a significant anisotropy in the electronic structure, resulting from the breaking of sixfold rotational symmetry in this material. Our findings demonstrate this member of Ti based kagome materials as a promising platform to explore novel emerging phenomena in the wider 𝐿⁢𝑛⁢Ti 3⁢ Bi 4 (𝐿⁢𝑛 = lanthanide) family.

electronic structure↗
DOE OSTI2021

Fermion Pair Dynamics in Open Quantum Systems

Three-body decay is a rare decay mode observed in a handful of unbound rare isotopes. The angular and energy correlations between emitted nucleons are of particular interest, as they provide invaluable information on the interplay between structure and reaction aspects of the nuclear open quantum system. To study the mechanism of two-nucleon emission, we developed a time-dependent approach that allows us to probe emitted nucleons at long times and large distances. We successfully benchmarked the new method against the Green’s function approach and applied it to low-energy two-proton and two-neutron decays. In particular, we studied the interplay between initial-state nucleon-nucleon correlations and final-state interaction. We demonstrated that the time evolution of the two-nucleon wave function is strongly impacted by the diproton/dineutron dynamics and that the correlations between emitted nucleons provide invaluable information on the dinucleon structure in the initial state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

6Li as a three-body system in the (p, 3He) reaction at astrophysical energies

Several astrophysical processes are governed by the occurrence of nuclear reactions involving light nuclei at energies below the Coulomb barrier. Among other effects, their understanding is challenged by the appearance of clustered structures in the ground-state configuration of some of these nuclei, which may have a significant impact on the reaction cross section. In this contribution, we focus on the Li 6 (p, He 3 ) He 4 reaction, to probe the role of clustered configurations of Li 6 . In particular, we consider a three-body ab-initio calculation, based on the hyperspherical harmonics (HH) method, of the Li 6 wave function (WF), together with a more phenomenological three-body model. We observe that the HH WF entails a degree of clustering much larger than obtained from the phenomenological WFs. However, the corresponding reaction cross section, evaluated as a direct two-nucleon transfer in distorted-wave Born approximation, still follows the scaling with the clustering strength already pointed out in a previous work [1] and exhibits an energy trend very similar to that obtained with realistic phenomenological WFs. This opens up interesting perspectives towards constraining

Perrotta, S.↗
DOE OSTI2023

Threshold resummation for computing large-x parton distribution through large-momentum effective theory

Parton distribution functions (PDFs) at large x are poorly constrained by high-energy experimental data, but extremely important for probing physics beyond standard model at colliders. We study the calculation of PDFs at large-x through large-momentum P z expansion of the lattice quasi PDFs. Similar to deep-inelastic scattering, there are two distinct perturbative scales in the threshold limit where the matching coefficient can be factorized into a space-like jet function at scale P z |1 – y| and a pair of heavy-light Sudakov form factors at scale P z . The matching formula allows us to derive a full renormalization group resummation of large threshold logarithms, and the result is consistent with the known calculation to the next-to-next to leading order (NNLO). This paves the way for direct large-x PDFs calculations in lattice QCD. As by-products, we find that the space-like jet function is related to a time-like version calculated previously through analytic continuation, and the heavy-light Sudakov form factor, calculated here to NNLO, is a universal object appearing as well in the large momentum expansion of quasi transverse-momentum-dependent PDFs and quasi wave-function amplitudes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Crystal field splittings and magnetic ground state of the square-lattice antiferromagnets YbBi 2 ⁢ClO 4 and YbBi 2 ⁢IO 4 with J eff = $\frac{1}{2}$

Here, we report on the crystal field level splitting and magnetic ground state of the J eff = $\frac{1}{2}$ square lattice antiferromagnets YbBi 2 ⁢ClO 4 and YbBi 2 ⁢IO 4 using powder inelastic neutron scattering (INS) and neutron diffraction measurements. Both compounds exhibit a well-isolated Γ 7 doublet ground state under a tetragonal crystal field environment, confirming a robust J eff = $\frac{1}{2}$ picture with slight XY-type anisotropic character in the g-tensor. Notably, the ground state wave functions closely resemble the Γ 7 doublet expected in the perfect cubic limit, consistent with the nearly cubic ligand configuration of eight O 2- ions surrounding Yb 3+ . Below T N = 0.21 K, YbBi 2 ⁢IO 4 exhibits a stripe long-range magnetic order characterized by an ordering wave vector q m = (1/2, 0, 0) or its symmetry-equivalent (0, 1/2, 0), with magnetic moments aligned along q m . The ordered moment is approximately 79% of the classical prediction, significantly larger than expected from the isotropic J 1 -J 2 model, suggesting the possible involvement of exchange anisotropy in explaining this observation. We show that symmetry-allowed XXZ and bond-dependent anisotropic exchange terms in a square lattice can play a critical role in stabilizing the stripe order and suppressing the moment reduction as observed. These findings establish YbBi 2 ⁢ClO 4 and YbBi 2 ⁢IO 4 as unique platforms for exploring rich J eff = $\frac{1}{2}$ magnetism from two less investigated perspectives: (i) on a square lattice and (ii) within a (nearly) cubic ligand environment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Quantum embedding for material chemistry based on domain separation and open subsystems

Our perspective considers two theories we recently proposed to perform quantum embedding calculations for chemical systems: domain-separated density functional theory (DSDFT) and locally coupled open subsystems (LCOS). The development includes both the fundamentals of each theory as well as potential applications, some technical aspects, and related challenges. DS-DFT is suited to study intramolecular effects, where one can apply a high level of theory (based on DFT or wave function theory) to a region of interest inside a molecule or solid, and lower level theory elsewhere, with smooth switching between the regions. LCOS, in contrast, is a fragment-based embedding which offers computational advantages to study intermolecular behavior such as electron hopping, spin-environment interaction, and charge transfer excitations. Yet, both theories can exchange roles when appropriate. In addition, these theories allow for control of computational scaling of their algorithms. We investigate paths to determine the charge-transfer operator used in LCOS, and suggest an auxiliary energy minimization that can provide a practical estimate to this operator. We also briefly discuss how to implement density fitting techniques in domain separation, and how domain separation can be used for pure wavefunction-based embedding.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Modern quantum chemistry with [Open]Molcas

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree–Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

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

Projector‐Based Quantum Embedding Study of Iron Complexes

Projection‐based embedding theory (PBET) is used to calculate and assess the challenging spin‐crossover energies for a selection of small Fe‐containing systems by embedding the metal center into the frozen potential of the ligands. MP2, CCSD, and CCSD(T) are embedded in potentials from the SCAN and r 2 SCAN functionals and compared with the canonical values for the constituent methods and previously reported reference values. Considering the PBET calculations as a correction for the underlying DFT, the embedding calculations are able to provided improvement for most cases. In some cases, the PBET methods are able to compensate for limitations in the wave function methods and produce results similar to more rigorous calculations from the literature. For the systems with spin‐crossover energies near zero, the current methodology fails to provide consistent improvement. In conclusion, the isolated recalculation of the electronic structure around the metal center when embedded into a DFT treatment of the ligand field shows promise as a pragmatic and lower cost treatment compared to the canonical treatment of the whole system of the difficult class of spin‐crossover complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗