Unoccupied natural orbitals of two-electron wave-functions.
Functions complementing set of occupied natural orbitals on two-electron wave functions
Nazaroff, G. V.↗
SEARCH · Search NASA
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.
A new-old approach to composite scalars with chiral fermion constituents
We develop a dynamical, Lorentz invariant theory of composite scalars in configuration space consisting of chiral fermions, interacting by the perturbative exchange of a massive “gluon” of coupling g 0 and mass M 0 2 (the coloron model). The formalism is inspired by, but goes beyond, old ideas of Yukawa and the Nambu-Jona-Lasinio (NJL) model. It yields a non-pointlike internal wave-function of the bound state, ϕ ( r ) , which satisfies a Schrödinger-Klein-Gordon (SKG) equation with eigenvalue μ 2 . For super-critical coupling, g 0 2 > g 0 c 2 , we have μ 2 < 0 leading to spontaneous symmetry breaking. The binding of chiral fermions is semiclassical, not loop-level as in NJL. The mass scale is determined by the interaction as in NJL. We mainly focus on the short-distance, large M 0 2 limit, yielding an NJL pointlike interaction, but the bound state internal wave-function, ϕ ( r → ) , remains spatially extended and dilutes ϕ ( 0 ) . This leads to power-law suppression of the induced Yukawa and quartic couplings and requires radically less fine-tuning of a hierarchy than does the NJL model. We include a discussion of loop corrections of the theory. A realistic top condensation model appears possible.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Unexpected origin of the quantum efficiency reduction in long-wavelength (In,Ga)N light-emitting diodes
Differential carrier lifetime (DCL) measurements were performed on c-plane In 𝑥 Ga 1−𝑥 N/GaN single-quantum-well (QW) light-emitting diodes (LEDs) with varying indium-content QWs (x = 13.5%, 16%, 22%), emitting with violet, blue, and green wavelengths. The recombination lifetimes of LEDs were found to increase with increasing indium composition, resulting in increased carrier densities n measured by DCL. Extraction of the A coefficients, which are assumed to not vary with n, and of the effective B(n) and C(n) coefficients of the ABC model of the internal quantum efficiency (IQE) of QWs showed no significant changes in the A coefficient with increasing indium content [In] in the In 𝑥 Ga 1−𝑥 N QW, and a reduction in the B(n) and C(n) coefficients with increasing [In]. When looking at the Shockley-Read-Hall (SRH) recombination rate An, the radiative recombination rate B(n)𝑛 2 , and the Auger-Meitner (AM) recombination rate C(n)𝑛 3 , we observed that at any given n, the SRH rate is the same for the three [In] measured, while the radiative rate decreases (by up to approximately 9 times) and the AM rate decreases (by up to approximately 7 times) with increasing indium content. This shows that the larger reduction in radiative recombination rates relative to nonradiative recombination rates with increasing [In] is the largest contributor to the decreased quantum efficiencies of LEDs (the “green gap”) at any given n. While some of the reduction in the effective recombination coefficients B(n) and C(n) can be explained by the reduced wave-function overlaps of the QW with increasing indium content, the larger reduction of B(n) relative to C(n) motivates further study of the intrinsic recombination coefficients 𝐵 0 and 𝐶 0 of bulk In 𝑥 Ga 1−𝑥 N alloys and how they are related to the effective recombination coefficients B(n) and C(n) in QWs. The relative contributions of nonradiative recombination are enhanced at any given current density due to the sublinear relationship between the carrier density n and the current density. Thus, the solution to the green gap, up to any intrinsic limitations set by the indium content of the QW, from an IQE perspective, requires the design of long-wavelength (In,Ga)N LEDs with improved wave-function overlaps that can operate with a lower carrier density at any given current density.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Response Functions of Correlated Systems within Green's Function Theory
Why use Green's functions as the fundamental variable? Wave-function (Psi) methods are king for high-fidelity and Density-functional (Rho) methods are very efficient (Kohn-Sham). Goldilocks principle: Green's function (G) methods straddle the Rho and Psi methods, intermediate in both accuracy and efficiency. Also, when interest lies in excitations & 2-particle properties: G-methods are natural - intrinsic to the theory.
DMFT↗
A physical basis for cosmological correlators from cuts
Significant progress has been made in our understanding of the analytic structure of FRW wavefunction coefficients, facilitated by the development of efficient algorithms to derive the differential equations they satisfy. Moreover, recent findings indicate that the twisted cohomology of the associated hyperplane arrangement defining FRW integrals overestimates the number of integrals required to define differential equations for the wave-function coefficient. We demonstrate that the associated dual cohomology is automatically organized in a way that is ideal for understanding and exploiting the cut/residue structure of FRW integrals. Utilizing this understanding, we develop a systematic approach to organize compatible sequential residues, which dictates the physical subspace of FRW integrals for any n -site, ℓ-loop graph. In particular, the physical subspace of tree-level FRW wavefunction coefficients is populated by differential forms associated to cuts/residues that factorize the integrand of the wavefunction coefficient into only flat space amplitudes. After demonstrating the validity of our construction using intersection theory, we develop simple graphical rules for cut tubings that enumerate the space of physical cuts and, consequently, differential forms without any calculation.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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↗
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↗
Out-of-time-ordered-correlators for the pure inverted quartic oscillator: classical chaos meets quantum stability
Out-of-time-ordered-correlators (OTOCs) have been suggested as a means to diagnose chaotic behavior in quantum mechanical systems. Recently, it was found that OTOCs display exponential growth for the inverted quantum harmonic oscillator, mirroring the fact that this system is classically and quantum mechanically unstable. In this work, I study OTOCs for the inverted anharmonic (pure quartic) oscillator in quantum mechanics, finding only oscillatory behavior despite the classically unstable nature of the system. For higher temperature, OTOCs seem to exhibit saturation consistent with a value of –2< x 2 > T < p 2 > T at late times. I provide analytic evidence from the spectral zeta-function and the WKB method as well as direct numerical solutions of the Schrödinger equation that the inverted quartic oscillator possesses a real and positive energy eigenspectrum, and normalizable wave-functions.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Natural top quark condensation (a redux)
The Nambu--Jona-Lasinio (NJL) model involves a pointlike 4-fermion interaction. While it gives a useful description of chiral dynamics (mainly in QCD), it nonetheless omits the crucially important internal wave-function of a two-body bound state, $\phi(r)$. This becomes significant near critical coupling where $\phi(r)$ extends to large distance, leading to dilution and suppression of induced couplings $\propto \phi(0)$, such as the Yukawa and quartic couplings, as well as reduced fine-tuning of a hierarchy. In top quark condensation, where the BEH boson is a $\bar{t}t$ bound state and we have a UV completion such as topcolor, we must go beyond the NJL model and include effects of $\phi(r)$. We provide a formulation of this for the Brout-Englert-Higgs boson, and find that it leads to an extended $\phi(r)$, a significantly reduced and natural composite scale of $M_0 \sim 6$ TeV, a successful prediction for the quartic coupling, $\lambda$, and fine tuning that is reduced to a few percent, providing a compelling candidate solution to the naturalness problem of the BEH boson. The theory is testable and new physics should begin to emerge on the multi-TeV mass scales and possibly accessible to the LHC.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Majorana zero modes in Y -shape interacting Kitaev wires
Motivated by the recent experimental realization of minimal Kitaev chains using quantum dots, we investigate the Majorana zero modes (MZM) in Y-shape Kitaev wires. We solve the associated Kitaev models analytically at the sweet spot (t h = Δ) and derive the exact form of MZM wave-functions in this geometry. We observed exotic multi-site MZMs located near the junction center, on the nearby edge sites of each leg. This result is important for potential braiding of Majoranas and the performance of Y-junctions made from arrays of quantum dots. Furthermore, we study the stability of local (single-site) and multi-site MZMs modes in the presence of Coulomb repulsion, using density matrix renormalization group theory. Our local density-of-states calculation shows that these multi-site MZMs are as equally topologically protected as the single-site MZMs when in the presence of Coulomb repulsion or when away from the sweet-spot.
36 MATERIALS SCIENCE↗
Simulating topological quantum gates in two-dimensional magnet-superconductor hybrid structures
The creation of topological quantum gates using Majorana zero modes—an outstanding problem in the field of topological quantum computing—relies on our ability to control the braiding process in time and space. Here, we propose two-dimensional magnet-superconductor hybrid structures as a new platformfor the successful implementation of topologically protected √σ z -, σ z - and σ x -quantum gates using Majorana zero modes. Employing a novel theoretical formalism to compute the full timedependent many-body wave-function and utilizing a braiding protocol motivated by recent advances in electron-spin-resonance techniques we simulate quantum gates in 2D systems up to 600 sites, on timescales from a few femto- to nanoseconds. We demonstrate that the braiding process can be visualized in time and space by computing the non-equilibrium local density of states, which is proportional to the time-dependent differential conductance measured in scanning tunneling spectroscopy experiments, allowing us to directly image Majorana world lines.
Superconducting properties and materials↗
GPAW: An open Python package for electronic structure calculations
We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe–Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn–Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Roadmap on methods and software for electronic structure based simulations in chemistry and materials
This Roadmap article provides a succinct, comprehensive overview of the state of electronic structure methods and software for molecular and materials simulations. Seventeen distinct sections collect insights by 51 leading scientists in the field. Each contribution addresses the status of a particular area, as well as current challenges and anticipated future advances, with a particular eye towards software related aspects and providing key references for further reading. Foundational sections cover density functional theory and its implementation in real-world simulation frameworks, Green's function based many-body perturbation theory, wave-function based and stochastic electronic structure approaches, relativistic effects and semiempirical electronic structure theory approaches. Subsequent sections cover nuclear quantum effects, real-time propagation of the electronic structure, challenges for computational spectroscopy simulations, and exploration of complex potential energy surfaces. The final sections summarize practical aspects, including computational workflows for complex simulation tasks, the impact of current and future high-performance computing architectures, software engineering practices, education and training to maintain and broaden the community, as well as the status of and needs for electronic structure based modeling from the vantage point of industry environments. Overall, the field of electronic structure software and method development continues to unlock immense opportunities for future scientific discovery, based on the growing ability of computations to reveal complex phenomena, processes and properties that are determined by the make-up of matter at the atomic scale, with high precision.
36 MATERIALS SCIENCE↗
Random insights into the complexity of two-dimensional tensor network calculations
Projected entangled pair states (PEPS) offer memory-efficient representations of some quantum many-body states that obey an entanglement area law and are the basis for classical simulations of ground states in two-dimensional (2d) condensed matter systems. However, rigorous results show that exactly computing observables from a 2d PEPS state is generically a computationally hard problem. Yet approximation schemes for computing properties of 2d PEPS are regularly used, and empirically seen to succeed, for a large subclass of (“not too entangled”) condensed matter ground states. Adopting the philosophy of random matrix theory, in this work, we analyze the complexity of approximately contracting a 2d random PEPS by exploiting an analytic mapping to an effective replicated statistical mechanics model that permits a controlled analysis at a large bond dimension. Through this statistical-mechanics lens, we argue that (i) although approximately sampling wave-function amplitudes of random PEPS faces a computational-complexity phase transition above a critical bond dimension, and (ii) one can generically efficiently estimate the norm and correlation functions for any finite bond dimension. Furthermore, these results are supported numerically for various bond-dimension regimes. It is an important open question whether the above results for random PEPS apply more generally also to PEPS representing physically relevant ground states.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Alternative to perturbative renormalization in ( 3 + 1 )-dimensional field theories
Perturbative renormalization provides the bedrock of understanding quantum field theories. In this work, I point out an alternative way of renormalizing quantum field theories, which is naturally encountered and well-known for the case of large N scalar field theories. In terms of bare parameters, this nonperturbative alternative renormalization differs qualitatively from its perturbative cousin: in the continuum limit, the bare coupling constant goes to zero instead of infinity, and there is no wave-function counterterm. Despite these differences, the resulting n -point functions of the theory are finite. I provide explicit results for alternative renormalization for the O(N) model and QCD with N f = 12 flavors in 3 + 1 dimensions. Published by the American Physical Society 2024
Astronomy & Astrophysics↗
Exact block encoding of imaginary time evolution with universal quantum neural networks
We develop a constructive approach to generate quantum neural networks capable of representing the exact thermal states of all many-body qubit Hamiltonians. The Trotter expansion of the imaginary time propagator is implemented through an exact block encoding by means of a unitary, restricted Boltzmann machine architecture. Marginalization over the hidden-layer neurons (auxiliary qubits) creates the nonunitary action on the visible layer. Then, we introduce a unitary deep Boltzmann machine architecture in which the hidden-layer qubits are allowed to couple laterally to other hidden qubits. We prove that this wave-function is closed under the action of the imaginary time propagator and, more generally, can represent the action of a universal set of quantum gate operations. We provide analytic expressions for the coefficients for both architectures, thus enabling exact network representations of thermal states without stochastic optimization of the network parameters. In the limit of large imaginary time, the yields the ground state of the system. The number of qubits grows linearly with the number of interactions and total imaginary time for a fixed interaction order. Both networks can be readily implemented on quantum hardware via midcircuit measurements of auxiliary qubits. If only one auxiliary qubit is measured and reset, the circuit depth scales linearly with imaginary time and number of interactions, while the width is constant. Alternatively, one can employ a number of auxiliary qubits linearly proportional to the number of interactions, and circuit depth grows linearly with imaginary time only. Every midcircuit measurement has a postselection success probability, and the overall success probability is equal to the product of the probabilities of the midcircuit measurements.
97 MATHEMATICS AND COMPUTING↗
Kolmogorov-Arnold wavefunctions
Here, this work investigates Kolmogorov-Arnold network-based (KAN) wave-function Ansätz as viable representations for quantum Monte Carlo simulations. Through systematic analysis of one-dimensional model systems, we evaluate their computational efficiency and representational power against established methods. Our numerical experiments suggest some efficient training methods and we explore how the computational cost scales with desired precision, particle number, and system parameters. Roughly speaking, KANs seem to be 10 times cheaper computationally than other neural-network-based Ansätz . We also introduce a novel approach for handling strong short-range potentials—a persistent challenge for many numerical techniques—which generalizes efficiently to higher-dimensional, physically relevant systems with short-ranged strong potentials common in atomic and nuclear physics.
1-dimensional systems↗
Breaking the Million-Electron and 1 EFLOP/s Barriers: Biomolecular-Scale Ab Initio Molecular Dynamics Using MP2 Potentials
The accurate simulation of complex biochemical phenomena has historically been hampered by the computational requirements of high-fidelity molecular-modeling techniques. Quantum mechanical methods, such as ab initio wave-function (WF) theory, deliver the desired accuracy, but have impractical scaling for modeling biosystems with thousands of atoms. Combining molecular fragmentation with MP2 perturbation theory, this study presents an innovative approach that enables biomolecular-scale ab initio molecular dynamics (AIMD) simulations at WF theory level. Leveraging the resolution-of-the-identity approximation for Hartree-Fock and MP2 gradients, our approach eliminates computationally intensive four-center integrals and their gradients, while achieving near-peak performance on modern GPU architectures. The introduction of asynchronous time steps minimizes time step latency, overlapping computational phases and effectively mitigating load imbalances. Utilizing up to 9,400 nodes of Frontier and achieving 59% (1006.7 PFLOP/s) of its double-precision floating-point peak, our method enables us to break the million-electron and 1EFLOP/s barriers for AIMD simulations with quantum accuracy.
Kurzak, Jakub↗