The choice of trial functions in the perturbation-variation method.
Nonlinear variational trial function for atomic hydrogen dynamic polarizability, discussing first order perturbed wave function and perturbation- variation method
HYDROGEN ATOM↗
SEARCH · Search NASA
Results for “Variational wave functional methods”
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.
Variational method for the ground state of liquid He super 4.
Liquid He 4 ground state studied by variation method, deriving radial distribution function
HELIUM 4↗
Gauge constrained algorithm of variational discrete action theory at N = 3 for the multiorbital Hubbard model
The recently developed variational discrete action theory (VDAT) provides a systematic variational approach to the ground state of the quantum many-body problem, where the quality of the solution is controlled by an integer N, and increasing N monotonically approaches the exact solution. VDAT can be exactly evaluated in the d = ∞ multiorbital Hubbard model using the self-consistent canonical discrete action theory (SCDA), which requires a self-consistency condition for the integer time Green's functions. Previous work demonstrates that N = 3 accurately captures multiorbital Mott/Hund physics at a cost similar to the Gutzwiller approximation. Here we employ a gauge constraint to automatically satisfy the self-consistency condition of the SCDA at N = 3, yielding an even more efficient algorithm with enhanced numerical stability. We derive closed form expressions of the gauge constrained algorithm for the multiorbital Hubbard model with general density-density interactions, allowing VDAT at N = 3 to be straightforwardly applied to the seven-orbital Hubbard model. We present results and a performance analysis using N = 2 and N = 3 for the SU(2N orb ) Hubbard model in d = ∞ with N orb = 2–8, and compare to numerically exact dynamical mean-field theory solutions where available. Finally, the developments in this work will greatly facilitate the application of VDAT at N = 3 to strongly correlated electron materials.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Revisiting excitation gaps in the fractional quantum Hall effect
Recent systematic measurements of the quantum well width dependence of the excitation gaps of fractional quantum Hall states in high mobility samples open the possibility of a better quantitative understanding of this important issue. Here we present what we believe to be accurate theoretical gaps including the effects of finite width and Landau level (LL) mixing. While theory captures the width dependence, there still remains a deviation between the calculated and the measured gaps, presumably caused by disorder. It is customary to model the experimental gaps of the n/(2n ± 1) states as Δ n/(2n+1) = Ce 2 /[(2n ± 1)ϵl]-Γ, where ϵ is the dielectric constant of the background semiconductor and l is the magnetic length; the first term is interpreted as the cyclotron energy of composite fermions and Γ as a disorder-induced broadening of composite-fermion LLs. Fitting the gaps for various fractional quantum Hall states, we find that Γ can be nonzero even in the absence of disorder.
36 MATERIALS SCIENCE↗
Variational study of the Kitaev-Heisenberg-Gamma model
Here we compute the low-energy excitation spectrum and the dynamical spin structure factor of the Kitaev-Heisenberg-Gamma model through a variational approach based on the exact fractionalized excitations of the pure Kitaev honeycomb model. This novel approach reveals the physical reason for the asymmetric stability of the Kitaev spin liquid phases around the ferromagnetic and antiferromagnetic Kitaev limits. Moreover, we demonstrate that the fractionalized excitations form bound states in specific regions of each Kitaev spin liquid phase and that certain phase transitions induced by Heisenberg and Gamma interactions are driven by the condensation of such a bound state. Remarkably, this bound state appears as a sharp mode in the dynamical spin structure factor, while its condensation patterns at the appropriate phase transitions provide a simple explanation for the magnetically ordered phases surrounding each Kitaev spin liquid phase.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Phase transitions of the ferroelectric (ND 4 ) 2 FeCl 5 ·D 2 O under magnetic field
Due to the strong coupling between magnetism and ferroelectricity, (ND 4 ) 2 FeCl 5 ·D 2 O exhibits several intriguing magnetic and electric phases. In this Letter, we include higher-order onsite anisotropic spin interactions to explain the ferroelectric phase transitions of (ND 4 ) 2 FeCl 5 ·D 2 O in a magnetic field and to produce the large weights of high-order harmonic components in the cycloid structure that are observed from neutron diffraction experiments. Moreover, we predict a new ferroelectric phase sandwiched between the FE II and FE III phases in a magnetic field. Furthermore, by emphasizing the importance of the higher-order spin anisotropic interactions, our work provides a framework to understand multiferroic materials with rich phase diagrams.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Spin and orbital metallic magnetism in rhombohedral trilayer graphene
We provide a theoretical interpretation of the metallic broken spin-valley (flavor) symmetry states recently discovered in hole-doped rhombohedral trilayer (ABC) graphene in large electric displacement fields. Our conclusions about the phase diagram and phase transitions combine insights from ABC graphene electronic structure models and mean-field theory, and are guided by the precise magneto-oscillation Fermi-surface-area measurements of recent experiments. We find that the principle of momentum-space condensation plays a key role in determining Fermi-surface reconstructions enabled by broken flavor symmetries when the single-particle bands imply thin annular Fermi seas. The reconstructed Fermi sea consists of one large outer Fermi-surface-enclosed majority-flavor states in reciprocal-space area A maj and one or more small inner holelike Fermi-surface-enclosed minority-flavor states in A min that are primarily responsible for nematic order. The competing ground states (valley-Ising, valley-XY, and spin-polarized state) have different A maj /A min and exchange energy maximizes this ratio and selects valley-XY nematic metal as the lowest-energy state. Here, we discuss how the nematic pockets explain the observed fractionalization of quantum oscillation frequencies, and propose anisotropic transport and the nonlinear Hall effect as additional observables.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Cascades between Light and Heavy Fermions in the Normal State of Magic-Angle Twisted Bilayer Graphene
We present a framework for understanding the cascade transitions and the Landau level degeneracies of twisted bilayer graphene. The Coulomb interaction projected onto narrow bands causes the charged excitations at an integer filling to disperse, forming new bands. If the excitation moves the filling away from the charge neutrality point, then it has a band minimum at the moiré Brillouin zone center with a small mass that compares well with the experiment; if towards the charge neutrality point, then it has a much larger mass and a higher degeneracy. At a nonzero density away from an integer filling the excitations interact. The system on the small mass side has a large bandwidth and forms a Fermi liquid. On the large mass side the bandwidth is narrow, the compressibility is negative and the Fermi liquid is likely unstable. This explains the observed sawtooth features in compressibility, the Landau fans pointing away from charge neutrality and their degeneracies. The framework sets the stage for superconductivity at lower temperatures.
36 MATERIALS SCIENCE↗
Quantum transitions of nematic phases in a spin-1 bilinear-biquadratic model and their implications for FeSe
Not Available
2-dimenional systems↗
Dirac Node Pinning from Dzyaloshinskii-Moriya Interactions in a Kagome Spin Liquid
Recent experiments on the Kagome spin liquid candidate material YCu 3 (OH) 6 Br 2 [Br 1−𝑦 (OH) 𝑦 ] suggest the presence of Dirac fermionic spinons near the magnetization plateau at 1/9. Theories suggest that the spinons are charge neutral spin-1/2 excitations, in a 2𝜋/3 flux, which triples the unit cell. Generally a gap is expected, and there is no symmetry protection for the Dirac nodes in this system. The question arises as to what causes the nodes and stabilizes them. In this work, we propose a node-creation and node-pinning mechanism driven by Dzyaloshinskii-Moriya (DM) interactions. Employing Gutzwiller-projected variational Monte Carlo calculations, we demonstrate that DM interactions induce a band closing phase transition in the spinon spectrum. There is a change in the Chern number when the bands are inverted. Together with the DM-generated internal gauge flux, the coupling to the spinon orbital magnetization counteracts the band reopening. Furthermore, this interplay energetically pins the Dirac nodes over a range of parameters, resulting in a pinning mechanism distinct from the usual one from symmetry protection.
Kagome lattice↗
Variational Monte Carlo calculations of n + H 3 scattering
A paramount goal in nuclear physics is to unify ab initio treatments of bound and unbound states. The position-space quantum Monte Carlo (QMC) methods have a long history of successful bound-state calculations in light systems but have seen minimal implementation in unbound systems. Here, we introduce a numerical method to improve the efficiency and accuracy of unbound-state calculations in QMC. As an initial application, we compute scattering observables for the smallest system available to probe three-body forces, the neutron-triton system, using variational Monte Carlo (VMC) wave functions. The method involves inferring long-range amplitudes in the wave function from integrals over the short-range region where all the particles interact. This approach using integral relations is well established in the literature; here, we develop it for the QMC framework. We validate our approach with a consistency check between short-range spectroscopic overlap functions computed from direct evaluation and from the integral relations; scattering amplitudes are long-range asymptotics of those overlaps. Comparison against published benchmark calculations using the same potential demonstrates that when applied to the current VMC wave functions, the integral method produces more accurate scattering observables than direct evaluation from the same variational wave function. However, it still differs noticeably from the exact results. Using additional interactions, we then present phase shifts and mixing parameters for the n + 3 H system. In particular, we present one of the first applications of the Norfolk family of local coordinate-space chiral potentials in unbound systems of A > 2. The Norfolk results accurately describe s-wave scattering but predict p-wave cross sections too large. Compared with previous QMC scattering calculations, the integral method avoids difficulties associated with the precise computation of energy differences and with convergence outside the interaction region, which is particularly severe in the variational calculation. Application of the integral method here paves the way for its use in Green's function Monte Carlo (GFMC) calculations. In GFMC, the wave functions are more accurate, but the high-precision convergence of their tails is slow, and there are additional difficulties in reading out amplitudes. Here, the integral methods will address both of those remaining problems.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Plasma-screening effects on the electron-impact excitation of hydrogenic ions in dense plasmas
Plasma-screening effects are investigated on electron-impact excitation of hydrogenic ions in dense plasmas. Scaled cross sections Z(exp 4) sigma for 1s yields 2s and 1s yields 2p are obtained for a Debye-Hueckel model of the screened Coulomb interaction. Ground and excited bound wave functions are modified in the screened Coulomb potential (Debye-Hueckel model) using the Ritz variation method. The resulting atomic wave functions and their eigenenergies agree well with the numerical and high-order perturbation theory calculations for the interesting domain of the Debye length not less than 10. The Born approximation is used to describe the continuum states of the projectile electron. Plasma screening effects on the atomic electrons cannot be neglected in the high-density cases. Including these effects, the cross sections are appreciably increased for 1s yields 2s transitions and decreased for 1s yields 2p transitions.
Jung, Young-Dae↗
Fast & accurate emulation of two-body scattering observables without wave functions
We combine Newton's variational method with ideas from eigenvector continuation to construct a fast & accurate emulator for two-body scattering observables. The emulator will facilitate the application of rigorous statistical methods for interactions that depend smoothly on a set of free parameters. Our approach begins with a trial K or T matrix constructed from a small number of exact solutions to the Lippmann–Schwinger equation. Subsequent emulation only requires operations on small matrices. We provide several applications to short-range potentials with and without the Coulomb interaction and partial-wave coupling. It is shown that the emulator can accurately extrapolate far from the support of the training data. When used to emulate the neutron-proton cross section with a modern chiral interaction as a function of 26 free parameters, it reproduces the exact calculation with negligible error and provides an over 300x improvement in CPU time.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Reexamining the variational two-particle reduced density matrix for nuclear systems
We report that in most nuclear many-body methods, observables are calculated using many-body wave functions explicitly. The variational two-particle reduced density matrix method is one of the few exceptions to the rule. Ground-state energies of both closed-shell and open-shell nuclear systems can indeed be evaluated by minimizing a constrained linear functional of the two-particle reduced density matrix. However, it has virtually never been used in nuclear theory, because nuclear ground states were found to be well overbound, contrary to those of atoms and molecules. Consequently, we introduced new constraints in the nuclear variational two-particle reduced density matrix method, developed recently for atomic and molecular systems. Our calculations then show that this approach can provide a proper description of nuclear systems where only valence neutrons are included. For the nuclear systems where both neutrons and protons are active, however, the energies obtained with the variational two-particle reduced density matrix method are still overbound. The possible reasons for the noticed discrepancies and solutions to this problem will be discussed.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Theoretical calculation of ozone vibrational infrared intensities
An ab initio dipole moment function for ozone has been computed using the CASSCF (complete active space self-consistent field) method, and forms the basis for a calculation of ozone infrared band intensities. Vibrational wave functions were generated using the variational method with potential energy surfaces derived from experimental force constants. Computed values of the permanent dipole moment, dipole moment derivatives, and infrared band strengths are all found to be in remarkably good agreement with experiment. Intensities are predicted for hot bands for which experimental values are unavailable, and implications for atmospheric ozone spectroscopy are discussed. As the dipole moment matrix element signs are now established for nearly all of the observed bands, further refinement of the dipole moment function is possible.
Adler-Golden, S. M.↗
Constrained-variation method in molecular quantum mechanics - Comparison of different approaches.
Constrained variation method in molecular quantum mechanics, presenting wave function, degree of constraint and energy loss values for lithium hydride molecule
QUANTUM MECHANICS↗
A light weight regularization for wave function parameter gradients in quantum Monte Carlo
The parameter derivative of the expectation value of the energy, ∂E/∂p, is a key ingredient in variational Monte Carlo (VMC) wave function optimization methods. In some cases, a naïve estimate of this derivative suffers from an infinite variance, which inhibits the efficiency of optimization methods that rely on a stable estimate of the derivative. In this work, we derive a simple regularization of the naïve estimator, which is trivial to implement in existing VMC codes, has finite variance, and a negligible bias, which can be extrapolated to zero bias with no extra cost. We use this estimator to construct an unbiased, finite variance estimation of ∂E/∂p for a multi-Slater–Jastrow trial wave function on the LiH molecule and in the optimization of a multi-Slater–Jastrow trial wave function on the CuO molecule. This regularized estimator is a simple and efficient estimator of ∂E/∂p for VMC optimization techniques.
36 MATERIALS SCIENCE↗
SchrödingerNet: A Universal Neural Network Solver for the Schrödinger Equation
Recent advances in machine learning have facilitated numerically accurate solution of the electronic Schrödinger equation (SE) by integrating various neural network (NN)-based wave function ansatzes with variational Monte Carlo methods. Nevertheless, such NN-based methods are all based on the Born–Oppenheimer approximation (BOA) and require computationally expensive training for each nuclear configuration. In this work, we propose a novel NN architecture, SchrödingerNet, to solve the full electronic-nuclear SE by defining a loss function designed to equalize local energies across the system. This approach is based on a translationally, rotationally and permutationally symmetry-adapted total wave function ansatz that includes both nuclear and electronic coordinates. Furthermore, this strategy not only allows for an efficient and accurate generation of a continuous potential energy surface at any geometry within the well-sampled nuclear configuration space, but also incorporates non-BOA corrections, through a single training process. Comparison with benchmarks of atomic and small molecular systems demonstrates its accuracy and efficiency.
Chemical calculations↗