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 379 records · Page 21

Variational quantum solver employing the PDS energy functional

In our previous work (J. Chem. Phys. 2020, 153, 201102),we reported a new class of quantum algorithms that are based on the quantum computation of the connected moment expansion to find the ground and excited state energies. In particular, the Peeters-Devreese-Soldatov (PDS) formulation is found variational and bearing the potential for further combining with the existing variational quantum infrastructure. Following this direction, here we propose a variational quantum solver employing the PDS energy gradient. In comparison with the usual variational quantum eigensolver (VQE) and the original static PDS approach, the proposed variational quantum solver offers an effective approach to achieve high ac-curacy at finding the ground state and its energy through the rotation of the trial wave function of modest quality guided by the low order PDS energy gradients, thus improves the ac-curacy and efficiency of the quantum simulation. We demonstrate the performance of the proposed variational quantum solver for toy models, H2molecule, and strongly correlated planar H4system in some challenging situations. In all the case studies, the proposed variational quantum approach out-performs the usual VQE and static PDS calculations even at the lowest order.

Peng, Bo↗

Electron impact ionization in dense plasmas

Here, the distorted-wave with exchange (DWE) method is employed to evaluate electron impact ionization cross sections in dense electron–ion plasmas. Bound and continuum electron wave functions are obtained from a non-relativistic average-atom code. Plots of DWE cross sections are presented for 1s electrons in Li and Be plasmas and 2p electrons in Na and Mg plasmas. For each of these elements, cross sections are evaluated at metallic density in a range of temperatures from 10 to 100 eV and at their respective melting points. Resonances in the cross sections appear near the incident energy threshold at high temperatures. The origin of these resonances is discussed. In general, the distorted wave (DW) method without exchange is found to be a good approximation to the DWE method for electron impact ionization calculations. In the resonance region, however, exchange effects are found to be very important and cannot be neglected.

74 ATOMIC AND MOLECULAR PHYSICS↗

Resonant Inelastic X-ray Scattering Calculations of Transition Metal Complexes Within a Simplified Time-Dependent Density Functional Theory Framework

We present a time-dependent density functional theory (TDDFT) approach to compute the light-matter couplings between two different manifolds of excited states relative to a common ground state in the context of 4d transition metal systems. These quantities are the necessary ingredients to solve the Kramers-Heisenberg (KH) equation for resonant inelastic X-ray scattering (RIXS) and several other types of two-photon spectroscopies. The procedure is based on the pseudo-wavefunction approach, where the solutions of a TDDFT calculation can be used to construct excited-state wave- functions, and on the restricted energy window approach, where a manifold of excited states can be rigorously defined based on the energies of the occupied molecular orbitals involved in the excitation process. Thus, the present approach bypasses the need to solve the costly TDDFT quadratic-response equations. We illustrate the applicability of the method to 4d transition metal molecular complexes by calculating the 2p4d RIXS maps of three representative ruthenium complexes and comparing them to experimental results. The method can capture all the experimental features in all three complexes to allow the assignment of the experimental peaks, and relative energies correct to within 0.6 eV at the cost of two independent TDDFT calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Photoelectron spectra of early 3 d -transition metal dioxide molecular anions from GW calculations

Photoelectron spectra of early 3 d -transition metal dioxide anions, Sc O 2 - , Ti O 2 - , V O 2 - , Cr O 2 - , and Mn O 2 - , are calculated using semilocal and hybrid density functional theory (DFT) and many-body perturbation theory within the GW approximation using one-shot perturbative and eigenvalue self-consistent formalisms. Different levels of theory are compared with each other and with available photoelectron spectra. We show that one-shot GW with a PBE0 starting point ( G 0 W 0 @PBE0) consistently provides very good agreement for all experimentally measured binding energies (within 0.1 eV–0.2 eV or less). We attribute this to the success of PBE0 in mitigating self-interaction error and providing good quasiparticle wave functions, which renders a first-order perturbative GW correction effective. One-shot GW calculations with a Perdew–Burke–Ernzerhof (PBE) starting point do poorly in predicting electron removal energies by underbinding orbitals with typical errors near 1.5 eV. A higher exact exchange amount of 50% in the DFT starting point of one-shot GW does not provide very good agreement with experiment by overbinding orbitals with typical errors near 0.5 eV. While not as accurate as G 0 W 0 @PBE0, the G -only eigenvalue self-consistent GW scheme with W fixed to the PBE level provides a reasonably predictive level of theory (typical errors near 0.3 eV) to describe photoelectron spectra of these 3 d -transition metal dioxide anions. Adding eigenvalue self-consistency also in W , on the other hand, worsens the agreement with experiment overall. Overall, our findings on the performance of various GW methods are discussed in the context of our previous studies on other transition metal oxide molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An intriguing intermediate state as a bridge between antiferroelectric and ferroelectric perovskites

The configuration of electric or magnetic dipoles embedded into a crystal lattice gives rise to the defining physical properties of a functional material. The prevailing antiferroelectric materials have so far been allusively assumed to be in a strict antiparallel electric polarization configuration, although they have shown huge potential applications for use in electronic devices. Herein, we have revealed for the first time the complex crystal structure of an intermediate state between the two conventional antiferroelectric (AFE) and ferroelectric (FE) states in PbZrO 3 -based perovskites by employing a combination of advanced synchrotron X-ray, neutron and three-dimensional (3D) electron diffraction techniques. We found that the intermediate state displays an intriguing flexible configuration of a sinusoidal wave functional with imbalanced polarization, different to the commonly believed antiparallel and balanced polarization. Such an intermediate state was concluded to act as a bridge between the conventional AFE and FE states. Our observations offer new insight into the AFE state and its transformation into the FE state. Such an intermediate state with a unique polarization configuration is believed to be a common feature that occurs in high-performance AFE perovskite materials.

36 MATERIALS SCIENCE↗

Magnetic Properties of Tetravalent Pu in the Perovskites BaPuO 3 and SrPuO 3

BaPuO 3 and SrPuO 3 were synthesized, and their structures were refined in the orthorhombic space group Pbnm , a common distortion from the classic $Pm\bar{3}$m over barm cubic perovskite. Magnetic-susceptibility measurements, obtained as a function of temperature over the range of 1.8-320 K, exhibit temperature-dependent behavior, with evidence of long-range magnetic order at temperatures higher than their lanthanide and actinide analogues: BaPuO 3 below 164(1) K and SrPuO 3 below 76(1) K. Effective moments of 1.66(10)μ B for BaPuO 3 and 1.84(8)μ B for SrPuO 3 were obtained by fitting their paramagnetic susceptibilities using the Curie-Weiss law. Further, both are below the free-ion value of 2.68 μ B expected for a Pu 4+ 5 $\mid$ 4 ground level. Ab initio wave function calculations, performed at the relativistic complete active space level including spin-orbit coupling and with an embedded cluster approach that neglects interactions between Pu centers, were used to generate embedded-cluster Pu 4+ magnetic susceptibilities. The calculations agree well with experimental data at higher temperatures, providing evidence that a single-ion representation is sufficient to account for the observed paramagnetic behavior without the need to invoke charge transfer, disproportionation, strong covalent bonding, or other more complex electronic behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A new generation of effective core potentials from correlated and spin–orbit calculations: Selected heavy elements

We introduce new correlation consistent effective core potentials (ccECPs) for the elements I, Te, Bi, Ag, Au, Pd, Ir, Mo, and W with 4d, 5d, 6s, and 6p valence spaces. These ccECPs are given as a sum of spin-orbit averaged relativistic effective potential (AREP) and effective spin–orbit (SO) terms. The construction involves several steps with increasing refinements from more simple to fully correlated methods. The optimizations are carried out with objective functions that include weighted many-body atomic spectra, norm-conservation criteria, and SO splittings. Transferability tests involve molecular binding curves of corresponding hydride and oxide dimers. The constructed ccECPs are systematically better and in a few cases on par with previous effective core potential (ECP) tables on all tested criteria and provide a significant increase in accuracy for valence-only calculations with these elements. Our study confirms the importance of the AREP part in determining the overall quality of the ECP even in the presence of sizable spin–orbit effects. The subsequent quantum Monte Carlo calculations point out the importance of accurate trial wave functions that, in some cases (mid-series transition elements), require treatment well beyond a single-reference.

Wang, Guangming (ORCID:0000000231327925)↗

Large-momentum effective theory

Since the parton model was introduced by Feynman more than 50 years ago, much has been learned about the partonic structure of the proton through a large body of high-energy experimental data and dedicated global fits. However, limited progress has been made in calculating partonic observables such as the parton distribution function (PDFs) from the fundamental theory of strong interactions, quantum chromodynamics (QCD). Recently some advocated for a formalism, large-momentum effective theory (LaMET), through which one can extract parton physics from the properties of the proton traveling at a moderate boost factor such as γ ~ 2 – 5 . The key observation behind this approach is that Lorentz symmetry allows the standard formalism of partons in terms of light-front operators to be replaced by an equivalent one with large-momentum states and time-independent operators of a universality class. With LaMET, the PDFs, generalized PDFs or generalized parton distributions, transverse-momentum-dependent PDFs, and light-front wave functions can all be extracted in principle from lattice simulations of QCD (or other nonperturbative methods) through standard effective field theory matching and running. Future lattice QCD calculations with exascale computational facilities could help one to understand the experimental data related to the hadronic structure, including those from the upcoming electron-ion colliders dedicated to exploring the partonic landscape of the proton. In this work, the progress made in the past few years in the development of the LaMET formalism and its applications is reviewed, with an emphasis on a demonstration of its effectiveness from initial lattice QCD simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Landau-Ginzburg theory of charge density wave formation accompanying lattice and electronic long-range ordering

Here we propose an analytical Landau-Ginzburg (LG) theory of the charge density waves coupled with lattice and electronic long-range order parameters. Examples of long-range order include the electronic wave function of superconducting Cooper pairs, structural distortions, electric polarization, and magnetization. We formulate the LG free energy density as a power expansion with respect to the charge density and other long-range order parameters as well as their spatial gradients and biquadratic coupling terms. We introduced a biquadratic coupling between the charge density gradient and long-range order parameters as well as nonlinear higher gradients of the long-range order parameters. The biquadratic gradient coupling is critical to the appearance of different spatially modulated phases in charge-ordered ferroics and high-temperature superconductors. We derived the thermodynamic conditions for the stability of the spatially modulated phases, which are the intertwined spatial waves of charge density and lattice/electronic long-range order. The analytical expressions for the energies of different phases, corresponding order parameters, charge density waves amplitudes, and modulation periods obtained in this paper can be employed to guide the comprehensive physical explanation, deconvolution, and Bayesian analysis of experimental data on quantum materials ranging from charge-ordered ferroics to high-temperature superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Light (anti)nuclei production in Pb-Pb collisions at s N N = 5.02 TeV

The measurement of the production of deuterons, tritons and 3 He and their antiparticles in Pb-Pb collisions at $\sqrt{s_{NN}}$= 5.02 TeV is presented in this article. The measurements are carried out at midrapidity (|y| < 0.5) as a function of collision centrality using the ALICE detector. The p T -integrated yields, the coalescence parameters and the ratios to protons and antiprotons are reported and compared with nucleosynthesis models. The comparison of these results in different collision systems at different center-of-mass collision energies reveals a suppression of nucleus production in small systems. In the Statistical Hadronisation Model framework, this can be explained by a small correlation volume where the baryon number is conserved, as already shown in previous fluctuation analyses. However, a different size of the correlation volume is required to describe the proton yields in the same data sets. The coalescence model can describe this suppression by the fact that the wave functions of the nuclei are large and the fireball size starts to become comparable and even much smaller than the actual nucleus at low multiplicities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nonlocal Coherent Optical Nonlinearities of a Macroscopic Quantum System

The optical responses of solids are typically understood to be local in space. Whether locality holds for the optical response of a macroscopic quantum system has, however, been largely unexplored. Here, we use multidimensional coherent spectroscopy at the optical diffraction limit to demonstrate nonlocal optical nonlinearities in a semiconductor microcavity. In conclusion, these nonlocal optical responses are both coherent and quantum in nature, deriving from the macroscopic length scale of confined exciton-polariton wave functions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nanoscale and Element-Specific Lattice Temperature Measurements Using Core-Loss Electron Energy-Loss Spectroscopy

Measuring nanoscale local temperatures, particularly in vertically integrated and multicomponent systems, remains challenging. Spectroscopic techniques like X-ray absorption and core-loss electron energy-loss spectroscopy (EELS) are sensitive to lattice temperature, but understanding thermal effects is nontrivial. This work explores the potential for nanoscale and element-specific core-loss thermometry by comparing the Si L2,3 edge’s temperature- dependent redshift against plasmon energy expansion thermometry (PEET) in a scanning TEM. Using density functional theory (DFT), time-dependent DFT, and the Bethe−Salpeter equation, we ab initio model both the Si L 2,3 and plasmon redshift. We find that the core-loss redshift occurs due to bandgap reduction from electron−phonon renormalization. Our results indicate that despite lower core-loss signal intensity compared to plasmon features, core-loss thermometry has key advantages and can be more accurate through standard spectral denoising. Specifically, we show that the Varshni equation easily interprets the core-loss redshift for semiconductors, which avoids plasmon spectral convolution for PEET in complex junctions and interfaces. We also find that core-loss thermometry is more accurate than PEET at modeling thermal lattice expansion in semiconductors, unless the specimen’s temperature-dependent dielectric properties are fully characterized. Furthermore, core-loss thermometry has the potential to measure nanoscale heating in multicomponent materials and stacked interfaces with elemental specificity at length scales smaller than the plasmon’s wave function.

Bethe–Salpeter equation↗

Exploring the Computational Aspects of Propylene Oligomerization Catalysis Using M′ 2 M Type Trimetallic MOF Nodes

Metal-organic frameworks have emerged as promising materials in the field of catalysis. They offer an optimal ground for screening catalysts and tailoring their catalytic properties. Here, in this work, via density functional theory (DFT) calculations, we investigated the catalytic activity of the trimetallic MOF nodes, M' 2 M for propylene oligomerization, by varying the active metal M from Sc to Cu with M' being Fe, aiming to grasp the impact of altering the active atoms on the catalyst's activity. Additionally, we examined how substituting the spectator atom, M', with other transition metals, i.e., from Sc to Cu, affects these energy barriers, keeping Ni as the active metal. We proposed several cases with lower or comparable energy barriers to the experimentally reported Fe 2 Ni trimetallic MOF node. In addition, we found a correlative relationship between spin-density from natural population analysis and energy barriers in the realm of C-C bond formation, whereby an elevation in spin-density is found to be inversely proportional to the magnitude of the energy barriers. Moreover, we calculated the energy barriers for C-C coupling and beta-hydride elimination using multireference NEVPT2 calculations on top of the CASSCF wave function to validate the rate-determining step that is predicted by DFT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A parallel strategy for density functional theory computations on accelerated nodes

Using the Löwdin orthonormalization of tall-skinny matrices as a proxy-app for wavefunction-based Density Functional Theory solvers, we investigate a distributed memory parallel strategy focusing on Graphics Processing Unit (GPU)-accelerated nodes as available on some of the top ranked supercomputers at the present time. Here we present numerical results in the strong limit regime, as it is particularly relevant for First-Principles Molecular Dynamics. We also examine how matrix product-based iterative solvers provide a competitive alternative to dense eigensolvers on GPUs, allowing to push the strong scaling limit of these computations to a larger number of distributed tasks. Our strategy, which relies on replicated Gram matrices and efficient collective communications using the NCCL library, leads to a time-to-solution under 0.5 s for the Löwdin orthonormalization of a tall-skinny matrix of 3000 columns on Summit at Oak Ridge Leadership Facility (OLCF). Given the similarity in computational operations between one iteration of a DFT solver and this proxy-app, this shows the possibility of solving accurately the DFT equations well under a minute for 3000 electronic wave functions, and thus perform First-Principles molecular dynamics of physical systems much larger than traditionally solved on CPU systems.

97 MATHEMATICS AND COMPUTING↗

Twisted symmetric trilayer graphene: Single-particle and many-body Hamiltonians and hidden nonlocal symmetries of trilayer moiré systems with and without displacement field

Here, we derive the Hamiltonian for trilayer moiré systems with the Coulomb interaction projected onto the bands near the charge neutrality point. Motivated by the latest experimental results, we focus on the twisted symmetric trilayer graphene (TSTG) with a mirror symmetry with respect to the middle layer. We provide a full symmetry analysis of the noninteracting Hamiltonian with a perpendicular displacement field coupling the band structure made otherwise of the twisted bilayer graphene (TBG) and the high-velocity Dirac fermions, and we identify a hidden nonlocal symmetry of the problem. In the presence of this displacement field, we construct an approximate single-particle model, akin to the tripod model for TBG, capturing the essence of noninteracting TSTG. We also derive more quantitative perturbation schemes for the low-energy physics of TSTG with displacement field, obtaining the corresponding eigenstates. This allows us to obtain the Coulomb interaction Hamiltonian projected in the active band TSTG wave functions and derive the full many-body Hamiltonian of the system. We also provide an efficient parametrization of the interacting Hamiltonian. Finally, we show that the discrete symmetries at the single-particle level promote the U (2) × U (2) spin-valley symmetry to enlarged symmetry groups of the interacting problem under different limits. The interacting part of the Hamiltonian exhibits a large U (4) × U (4) × U (4) × U (4) symmetry in the chiral limit. Moreover, by identifying a symmetry which we dub spatial many-body charge conjugation, we show that the physics of TSTG is symmetric around charge neutrality.

2-dimensional systems↗

Interplay of broken symmetry and delocalized excitations in the insulating state of 1 T –Ta⁢S 2

Coexistence of localized and extended excitations is central to the macroscopic properties of correlated materials. For 5⁢d transition-metal compounds, electron correlations alone generally do not lead to a metal-insulator (Mott) transition, with insulating behavior usually resulting from their coupling with magnetic ordering and/or structural distortions. 1⁢T–Ta⁢S 2 is a prototypical example of such correlated insulating behavior, with a high-symmetry metallic phase transforming into a distorted, charge-density wave (CDW) insulating state at low temperatures. The nature of the insulating phase as well as the existence and relevance of the localized electron physics remains debated. We resolve this standing controversy in 1⁢T–Ta⁢S 2 combining resonant inelastic x-ray spectroscopy and first-principles calculations. We observe five electronic excitations arising from the interband transitions of the Ta 5d orbitals and the S 3⁢p ligand state, with none of the excitations on the order of the Mott gap. These excitations cannot be explained within the framework of standard multiplet calculations that assume a localized wave function, but instead, are captured by a band-theory framework accounting for the low symmetry of the crystal field in the CDW state. Our findings suggest that the electronic properties of 1T–Ta⁢S 2 are dominated in the visible by both plasmonic quasiparticles and interband transitions, with no resonance associated with a putative Mott transition observed in the 0–3 eV energy range. Finally, our discovery provides insights into the electron localization and the Mott vs band-insulator debate in 1T–Ta⁢S 2 and other transition-metal materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Search for the 1/2 + intruder state in 35 P

The excitation energy of deformed intruder states (specifically the 2p2h bandhead) as a function of proton number Z along N = 20 is of interest both in terms of better understanding the evolution of nuclear structure between spherical Ca 40 and the Island of Inversion nuclei, and for benchmarking theoretical descriptions in this region. At the center of the N = 20 Island of Inversion, the n p n h (where n = 2 , 4 , 6 ) neutron excitations across a diminished N = 20 gap result in deformed and collective ground states, as observed in Mg 32 . In heavier isotones, n p n h excitations do not dominate in the ground states but are present in the relatively low-lying level schemes. With the aim of identifying the expected 2 p 2 h x s 1 / 2 + state in P 35 , the only N = 20 isotone for which the neutron 2p2h excitation bandhead has not yet been identified, the S 36 ( d , He 3 ) P 35 reaction has been revisited in inverse kinematics with the HELical Orbit Spectrometer (HELIOS) at the Argonne Tandem Linac Accelerator System (ATLAS). While a candidate state has not been located, an upper limit for the transfer reaction cross section to populate such a configuration within a 2.5 to 3.6 MeV energy range provides a stringent constraint on the wave function compositions in both S 36 and P 35 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Pb-Apatite Framework as a Generator of Novel Flat-Band CuO-Based Physics

Based upon density functional theory (DFT) calculations, we present the basic electronic structure of CuPb 9 (PO 4 ) 6 O (Cu-doped lead apatite, aka LK-99), in two scenarios: (1) where the structure is constrained to the P3 symmetry and (2) where no symmetry is imposed. At the DFT level, the former is predicted to be metallic while the latter is found to be a charge-transfer insulator. In both cases the filling of these states is nominally d 9 , consistent with the standard Cu 2+ valence state, and Cu with a local magnetic moment of order 0.7 μ B . In the metallic case we find these states to be unusually flat (∼0.2 eV dispersion), giving a very high density of electronic states (DOS) at the Fermi level that we argue can be a host for novel electronic physics. The flatness of the bands is the likely origin of symmetry-lowering gapping possibilities that would remove the spectral weight from E F . Motivated by some initial experimental observations of metallic or semiconducting behavior, we propose that disorder (likely structural) is responsible for closing the gap. Here, we consider a variety of possibilities that could possibly close the charge-transfer gap but limit consideration to kinds of disorder that preserve electron count. Of the possible kinds we considered (spin disorder, O populating vacancy sites, and Cu on less energetically favorable Pb sites), the local Cu moment, and consequently the charge-transfer gap, remains robust. We conclude that disorder responsible for metallic behavior entails some kind of doping where the electron count changes. Further, we claim that the emergence of the flat bands should be due to weak wave function overlap between the orbitals on Cu and O sites, owing to the directional character of the constituent orbitals. Therefore, finding an appropriate host structure for minimizing hybridization between Cu and O while allowing them to still weakly interact should be a promising route for generating flat bands at E F which can lead to interesting electronic phenomena, regardless of whether LK-99 is a superconductor.

36 MATERIALS SCIENCE↗