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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Electronic structure and exchange interactions in altermagnetic MnGeP 2 in the quasiparticle self-consistent G W approach

The quasiparticle self-consistent 𝐺⁡𝑊 method is used to study the electronic band structure, optical dielectric function, and exchange interactions in chalcopyrite, 𝐼⁢ ‾ 4 ⁢2⁢𝑑, structure MnGeP 2 . The material is found to be an antiferromagnetic semiconductor with a lowest direct gap of 2.44 eV at the Γ point and a lower indirect gap of 1.87 eV from Γ to 𝑀. The material is an altermagnet because the two magnetic atoms of opposite spin are related by a twofold rotation operation perpendicular to the main fourfold rotation inversion axis. The spin splittings along a low symmetry line like 𝑃⁢𝑁 is sizable while at k points on the diagonal mirror planes or on the twofold symmetry axes the spin splitting is zero. The exchange interactions are calculated using a linear response approach. The antiferromagnetic exchange interaction between nearest neighbors in the primitive unit cell is dominating and found to be slightly decreasing upon carrier doping but not sufficiently to change the interaction to become ferromagnetic. The bare (noninteracting) and interacting transverse spin susceptibilities, which provide interatomic site exchange interactions after averaging over the muffin-tin spheres, are calculated from the 𝐺⁡𝑊 band structure and wave functions. From these exchange interactions, the spin wave spectra are obtained along the high symmetry lines and the Néel temperature is calculated using the mean-field and Tyablikov estimations. The dielectric function and the optical absorption spectra are calculated including excitonic effects using the Bethe Salpeter equation. The exchange interactions around Mn Ge defect sites is also studied. While we find it can generate ferromagnetic interactions with neighboring spins, we did not find direct evidence of producing an overall ferromagnetic phase. First, if Mn antisites are introduced by exchanging Mn with a nearby Ge, the interactions stay largely antiferromagnetic. Second, when we add additional Mn, in other words in Mn-rich stoichiometry, the Mn Ge antisites produce a strong ferromagnetic interaction primarily with the Mn in the same basal plane but weaker ferromagnetic interaction with adjacent plane Mn. The interactions between regular lattice Mn stay antiferromagnetic as before and thus favor keeping the antiferromagnetic order along the [001] direction. Adding Mn antisites, however, does lead to a metallic band structure.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Accuracy of ghost rotationally invariant slave-boson and dynamical mean field theory as a function of the impurity-model bath size

Here, we compare the accuracy of the ghost rotationally invariant slave-boson (g-RISB) theory and dynamical mean field theory (DMFT) on the single-band Hubbard model, as a function of the number of bath sites in the embedding impurity Hamiltonian. Our benchmark calculations confirm that the accuracy of g-RISB can be systematically improved by increasing the number of bath sites, similar to DMFT. With a few bath sites, we observe that g-RISB is systematically more accurate than DMFT for the ground-state observables. On the other hand, the relative accuracy of these methods is generally comparable for the quasiparticle weight and the spectral function. As expected, we observe that g-RISB satisfies the variational principle in infinite dimensions, as the total energy decreases monotonically towards the exact value as a function of the number of bath sites, suggesting that the g-RISB wave function may approach the exact ground state in infinite dimensions. Our results suggest that the g-RISB is a promising method for first-principles simulations of strongly correlated matter, which can capture the behavior of both static and dynamical observables, at a relatively low computational cost.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Scattering phase shifts and mixing angles for an arbitrary number of coupled channels on the lattice

We present a lattice method for determining scattering phase shifts and mixing angles for the case of an arbitrary number of coupled channels. Previous nuclear lattice effective field theory simulations were restricted to mixing of up to two partial waves for scattering of two spin-1/2 particles, which is insufficient for analyzing nucleon-nucleus or nucleus-nucleus scattering processes. In the proposed method, the phase shifts and mixing angles are extracted from the radial wave functions obtained by projecting the three-dimensional lattice Hamiltonian onto the partial wave basis. In this work, we use a spherical wall potential as a boundary condition along with a channel-mixing auxiliary potential to construct the full-rank S matrix. Our method can be applied to particles with any spin, but we focus here on scattering of two spin-1 bosons involving up to four coupled channels. For a considered test potential, the phase shifts and mixing angles extracted on the lattice are shown to agree with the ones calculated by solving the Schrödinger equation in the continuum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

On exact-WKB analysis, resurgent structure, and quantization conditions

There are two well-known approaches to studying nonperturbative aspects of quantum mechanical systems: saddle point analysis of the partition functions in Euclidean path integral formulation and the exact-WKB analysis based on the wave functions in the Schrödinger equation. In this work, based on the quantization conditions obtained from the exact-WKB method, we determine the relations between the two formalism and in particular show how the two Stokes phenomena are connected to each other: the Stokes phenomenon leading to the ambiguous contribution of different sectors of the path integral formulation corresponds to the change of the “topology” of the Stoke curves in the exact-WKB analysis. We also clarify the equivalence of different quantization conditions including Bohr-Sommerfeld, path integral and Gutzwiller’s ones. In particular, by reorganizing the exact quantization condition, we improve Gutzwiller’s analysis in a crucial way by bion contributions (incorporating complex periodic paths) and turn it into an exact result. Furthermore, we argue the novel meaning of quasi-moduli integral and provide a relation between the Maslov index and the intersection number of Lefschetz thimbles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Theoretical Evaluation of Metal–Ligand Bonding in Neptunium Compounds in Relation to 237 Np Mössbauer Spectroscopy

The 237 Np Mössbauer isomer shift and quadrupole splitting (QS) are powerful probes for the metal–ligand bonding of neptunium, a 5f-element of vital importance in the nuclear fuel cycle. A large set of Np compounds with different oxidation states (III) to (VII) is studied to investigate, by first-principles calculations, isomer shifts and the QS trends in relation to the Np oxidation state. Natural Bond Orbital analysis reveals that in addition to donation bonding to the 5f shell, participation of the 6d and 7s neptunium shells in covalent (donation) bonding substantially impacts the isomer shifts. The isomer shift cannot be interpreted solely by the 5f shell electron count. The isomer shift for Np(II) compounds is estimated to be in the range of 31–34 mm/s, less positive than for Np(III) compounds. For the QS, density functional calculations fail to reproduce the quadrupole splitting for some Np(VI) ionic solids. A multiconfigurational wave function approach reproduces the observed QS trends. The calculations give a semiquantitative interpretation of the trends for Np oxidation states (V) to (VII). Furthermore, the contrasting QS for standard and “reverse” neptunyl(VI), at the opposite extremes of the observed QS scale, arises predominantly from the different crystal environments.

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

Automation of Active Space Selection for Multireference Methods via Machine Learning on Chemical Bond Dissociation

Predicting and understanding the chemical bond is one of the major challenges of computational quantum chemistry. Kohn–Sham density functional theory (KS-DFT) is the most common method, but approximate density functionals may not be able to describe systems where multiple electronic configurations are equally important. Multiconfigurational wave functions, on the other hand, can provide a detailed understanding of the electronic structures and chemical bonds of such systems. In the complete active space self-consistent field (CASSCF) method, one performs a full configuration interaction calculation in an active space consisting of active electrons and active orbitals. However, CASSCF and its variants require the selection of these active spaces. This choice is not black box; it requires significant experience and testing by the user, and thus active space methods are not considered particularly user-friendly and are employed only by a minority of quantum chemists. Our goal is to popularize these methods by making it easier to make good active space choices. Here, we present a machine learning protocol that performs an automated selection of active spaces for chemical bond dissociation calculations of main group diatomic molecules. The protocol shows high prediction performance for a given target system as long as a properly correlated system is chosen for training. Good active spaces are correctly predicted with a considerably better success rate than random guess (larger than 80% precision for most systems studied). Our automated machine learning protocol shows that a “black-box” mode is possible for facilitating and accelerating the large-scale calculations on multireference systems where single-reference methods such as KS-DFT cannot be applied.

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

Quantum tensor network algorithms for evaluation of spectral functions on quantum computers

We investigate quantum algorithms derived from tensor networks to simulate the static and dynamic properties of quantum many-body systems. Using a sequentially prepared quantum circuit representation of a matrix product state (MPS) that we call a quantum tensor network (QTN), we demonstrate algorithms to prepare ground and excited states on a quantum computer and apply them to molecular nanomagnets (MNMs) as a paradigmatic example. In this setting, we develop two approaches for extracting the spectral correlation functions measured in neutron-scattering experiments: (a) a generalization of the SWAP test for computing wave function overlaps and, (b) a generalization of the notion of matrix product operators to the QTN setting which generates a linear combination of unitaries. The latter method is discussed in detail for translationally invariant spin-half systems, where it is shown to reduce the qubit resource requirements compared with the SWAP method and may be generalized to other systems. We demonstrate the versatility of our approaches by simulating spin-1/2 and spin-3/2 MNMs, with the latter being an experimentally relevant model of a Cr$^{3+}_{8}$ ring. Here, our approach has qubit requirements that are independent of the number of constituents of the many-body system and scale only logarithmically with the bond dimension of the MPS representation, making them appealing for implementation on near-term quantum hardware with mid-circuit measurement and reset.

Neutron scattering↗
DOE OSTI2020

Stochastic Inversion of Gaussian Random Media Using Transverse Coherence Functions for Reflected Waves: Theory and Method

The transverse coherence functions (TCFs) of phase and amplitude fluctuations of a seismic wave are powerful to estimate the spatial distribution, length scales, and strength of random heterogeneities. However, TCFs have been formulated for transmitted waves only, not for reflected waves. In this paper, we derive reflection TCFs for Gaussian random media. Furthermore, we propose to invert for Gaussian random media using the reflection TCFs based on the grid search. We validate the new reflection TCF formulas using 2D finite-difference numerical experiments. The numerical example also illustrates the feasibility and efficiency of the inversion. The stochastic inversion using reflected waves can be used in both exploration and global seismology.

58 GEOSCIENCES↗
DOE OSTI2020

Semi-supervised Machine Learning Enables the Robust Detection of Multireference Character at Low Cost

Multireference (MR) diagnostics are common tools for identifying strongly correlated electronic structure that makes single-reference (SR) methods (e.g., density functional theory or DFT) insufficient for accurate property prediction. However, MR diagnostics typically require computationally demanding correlated wave function theory (WFT) calculations, and diagnostics often disagree or fail to predict MR effects on properties. To overcome these challenges, we introduce a semi-supervised machine learning (ML) approach with virtual adversarial training (VAT) of an MR classifier using 15 WFT and DFT MR diagnostics as inputs. In semi-supervised learning, only the most extreme SR or MR points are labeled, and the remaining point labels are learned. The resulting VAT model outperforms the alternatives, as quantified by the distinct property distributions of SR- and MR-classified molecules. Additionally, to reduce the cost of generating inputs to the VAT model, we leverage the VAT model’s robustness to noisy inputs by replacing WFT MR diagnostics with regression predictions in an MR decision engine workflow that preserves excellent performance. We demonstrate the transferability of our approach to larger molecules and those with distinct chemical composition from the training set. This MR decision engine demonstrates promise as a low-cost, high-accuracy approach to the automatic detection of strong correlation for predictive high-throughput screening.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Electronic structure and optical properties of quantum crystals from first principles calculations in the Born–Oppenheimer approximation

In this work, we develop a formalism to accurately account for the renormalization of the electronic structure due to quantum and thermal nuclear motions within the Born–Oppenheimer approximation. We focus on the fundamental energy gap obtained from electronic addition and removal energies from quantum Monte Carlo calculations in either the canonical or grand-canonical ensembles. The formalism applies as well to effective single electron theories such as those based on density functional theory. We show that the electronic (Bloch) crystal momentum can be restored by marginalizing the total electron–ion wave function with respect to the nuclear equilibrium distribution, and we describe an explicit procedure to establish the band structure of electronic excitations for quantum crystals within the Born–Oppenheimer approximation. Based on the Kubo–Greenwood equation, we discuss the effects of nuclear motion on optical conductivity. Our methodology applies to the low temperature regime where nuclear motion is quantized and, in general, differs from the semi-classical approximation. We apply our method to study the electronic structure of C2/c-24 crystalline hydrogen at 200 K and 250 GPa and discuss the optical absorption profile of hydrogen crystals at 200 K and carbon diamond at 297 K.

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

Eikonal spin-dependent odderon and gluon Sivers function of a proton and its small-𝑥 evolution

The matrix element in the proton of the eikonal Odderon operator, with a helicity flip, has been shown to correspond to the dipole gluon Sivers function. We employ a three quark light-front model of the proton to determine the Sivers function at moderately small 𝑥 0 ∼ 0.1 and transverse momentum 𝑘 ⊥ ≲ 1 GeV. The model light-cone (LC) wave function predicts the properties of 𝑥⁢𝑓$^{⊥⁢𝑔}_{1⁢𝑇⁡}$(𝑥,𝑘 ⊥ ) such as its overall magnitude, the position of its peak in 𝑘 ⊥ , and its behavior at small 𝑘 ⊥ . We then compute numerically the Balitsky-Fadin-Kuraev-Lipatov anomalous dimension characterizing the power-law tail at 𝑘 ⊥ ≳ 1.5 GeV of the gluon Sivers function at small (but preasymptotic) LC momentum fractions 𝑥 ∼ 𝑥 0 ⁢𝑒 −1/𝛼 𝑠 : 𝑥⁢𝑓$^{⊥⁢𝑔}_{1⁢𝑇⁡}$⁡(𝑥,𝑘 ⊥ ) ∼ 𝑘$^{−3.3}_{⊥}$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Collins-Soper kernel and reduced soft function in lattice QCD

We evaluate the Collins-Soper kernel and the reduced soft function in lattice QCD, incorporating 𝒪⁡(𝛼 𝑠 ) matching corrections. The calculation relies on the evaluation of the quasitransverse momentum–dependent wave function with asymmetric staple-shaped quark bilinear operators and four-point meson form factors. These quantities are computed nonperturbatively using two 𝑁 𝑓 =2 + 1 + 1 twisted-mass fermion ensembles with the same lattice spacing of 𝑎 = 0.093 fm: the first ensemble has a lattice size of 24 3 × 48 and a pion mass of 346 MeV, and the second one has a lattice size of 32 3 × 64 and a pion mass of 261 MeV. The Collins-Soper kernel and the soft function are needed for the determination of the transverse momentum–dependent parton distribution functions.

Lattice QCD↗
DOE OSTI2022

Crust and Upper Mantle Structure Beneath the Eastern United States

Abstract The Eastern United States (EUS) has a complex geological history and hosts several seismic active regions. We investigate the subsurface structure beneath the broader EUS. To produce reliable images of the subsurface, we simultaneously invert smoothed P‐wave receiver functions, Rayleigh‐wave phase and group velocity measurements, and Bouguer gravity observations for the 3D shear‐wave speed. Using surface‐wave observations (3–250 s) and spatially smoothed receiver functions, our velocity models are robust, reliable, and rich in detail. The shear‐wave velocity models fit all three types of observations well. The resulting velocity model for the eastern U.S. shows thinner crust beneath New England, the east coast, and the Mississippi Embayment (ME). A relatively thicker crust was found beneath the stable North America craton. A relatively slower upper mantle was imaged beneath New England, the east coast, and western ME. A comparison of crust thickness derived from our model against four recent published models shows first‐order consistency. A relatively small upper mantle low‐speed region correlates with a published P‐wave analysis that has associated the anomaly with a 75 Ma kimberlite volcanic site in Kentucky. We also explored the relationship between the subsurface structure and seismicity in the eastern U.S. We found that earthquakes often locate near regions with seismic velocity variations, but not universally. Not all regions of significant subsurface wave speed changes are loci of seismicity. A weak correlation between upper mantle shear velocity and earthquake focal mechanism has been observed.

58 GEOSCIENCES↗
DOE OSTI2026

Majorana Edge Modes in Isolated Wires

Topological superconductors are believed to host exotic quasiparticle excitations known as Majorana zero modes (MZMs), with much of the evidence based on BCS mean-field theory. The direct application of mean-field arguments is tenuous in finite, isolated systems relevant in some experiments. Here, we develop a new correlation-based method for identifying MZMs in interacting, number-conserving systems. Using the density matrix renormalization group, we study fermion number-conserving models with long-range interactions, which under periodic boundary conditions exhibit robust topological and nontopological superconductivity, tuned by the strength of interaction [Ortiz et al., Phys. Rev. Lett. 113, 267002 (2014)]. We find evidence that, on the topological side, Majorana edge modes appear in open chains, manifesting as the vanishing of the energy splitting between odd- and even-parity ground states with increasing system size. Additionally, off-diagonal two-point correlation functions show nonlocal, parity-dependent edge effects. These correlations reveal the spatial structure of Majorana modes in the many-body wave function. We show that the correlation diagnostic applies broadly, including to short-range interacting models, where topological superconductivity is more fragile due to the absence of a bulk excitation gap.

Thomas-Markarian, Jaden↗
DOE OSTI2025

Datasets for Custom-trained Machine-learning Interatomic Potentials: Nitric Acid Aqueous Solution

This dataset was generated using an iterative active learning strategy with the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials (MLIPs) for aqueous nitric acid. Each active-learning cycle consisted of three stages: (1) training, (2) exploration, and (3) labeling. The initial training set comprised approximately 800 randomly selected configurations from a previous study by Lewis et al. (https://doi.org/10.1021/jp205510q), which investigated nitric acid solutions at 2, 3, 4, and 5 mol/L. For all configurations, single-point calculations of atomic forces and total energies were performed at the quantum density functional theory BLYP-D2 and PBE-D3 levels of theory using the CP2K Quickstep module. Valence electrons were treated explicitly, while core electrons on all atoms were represented by norm-conserving Goedecker–Teter–Hutter (GTH) pseudopotentials. Long-range dispersion interactions were accounted for using Grimme dispersion corrections. Wave functions were expanded in a mixed Gaussian-and-plane-wave scheme using TZV2P-MOLOPT basis sets for all elements and an 800 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent field convergence was accelerated using orbital transformation and Direct Inversion in the Iterative Subspace, with a convergence threshold of 10^{-6}. All single-point calculations were carried out in periodic orthorhombic cells whose dimensions match those of the molecular configurations sampled from earlier trajectories. The CELL_REF keyword in CP2K was used to define a fixed reference cell, ensuring consistency in the reference data used for MLIP training, particularly when cell fluctuations are present in NpT simulations. The resulting high-fidelity energies and forces constitute the ground-truth labels used to train the MLIPs contained in this dataset.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2024

Tetraquarks made of sufficiently unequal-mass heavy quarks are bound in QCD

Tetraquarks, bound states composed of two quarks and two antiquarks, have been the subject of intense study but are challenging to understand from first principles. We apply variational and Green’s function Monte Carlo methods to compute tetraquark ground-state energies in potential nonrelativistic QCD using a wide range of color and spatial wave functions. We find no evidence for bound tetraquarks composed of equal-mass quarks and antiquarks. Conversely, we find clear evidence for the existence of bound tetraquarks for sufficiently unequal quark/antiquark mass ratios at all overall mass scales where our effective theory results are applicable. We predict the critical mass ratios for bound state formation and study tetraquark bound states’ spatial and color structure at leading order and next-to-leading order in potential nonrelativistic QCD. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Mixed excitonic nature in water-oxidized BiVO 4 surfaces with defects

BiVO 4 is a promising photocatalyst for efficient water oxidation, with surface reactivity determined by the structure of active catalytic sites. Surface oxidation in the presence of oxygen vacancies induces electron localization, suggesting an atomistic route to improve the charge transfer efficiency within the catalytic cycle. In this paper, we study the effect of oxygen vacancies on the electronic and optical properties at BiVO 4 surfaces upon water oxidation. We use density functional theory and many-body perturbation theory to explore the change in the electronic and quasiparticle energy levels and to evaluate the electron-hole coupling as a function of the underlying structure. Further we show that while the presence of defects alters the atomic structure and largely modifies the wave-function nature, leading to defect-localized states at the quasiparticle gap region, the optical excitations remain largely unchanged due to the substantial hybridization of defect and nondefect electron-hole transitions. Our findings suggest that defect-induced surface oxidation supports improved electron transport, both through bound and tunable electronic states and via a mixed nature of the optical transitions, expected to reduce electron-hole defect trapping.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Mass of a weakly measured photon

Bohmian mechanics has garnered significant attention as an interpretation of quantum theory since the paradigmatic experiments which inferred the average trajectories of photons in the nonrelativistic regime. These experiments were largely motivated by Wiseman's formulation of Bohmian mechanics, which grounded these trajectories in weak measurements. Recently, Wiseman's framework was extended to the relativistic regime by expressing the velocity field of single photons in terms of weak values of the photon energy and momentum. Here, we propose an operational, weak value-based definition for the Bohmian “local mass” of relativistic single particles. For relativistic wave functions satisfying the scalar Klein-Gordon equation, this mass coincides with the effective mass defined by de Broglie in his relativistic pilot-wave theory, a quantity closely connected with the quantum potential that is responsible for Bohmian trajectory self-bending and the anomalous photoelectric effect. Here, we demonstrate the relationship between the photon trajectories and the mass in an interferometric setup.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗