Zeeman effect and configuration interaction in germanium.
Zeeman effect of germanium I and II, discussing configuration-interaction, deriving intermediate- coupling wave functions
ELECTRON SPECTRUM↗
SEARCH · Search NASA
Results for “Wave Functions”
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.
Discontinuity of the exchange-correlation potential and the functional derivative of the noninteracting kinetic energy as the number of electrons crosses integer boundaries in Li, Be, and B
Accurate densities were determined from configuration interaction wave functions for atoms and ions of Li, Be, and B with up to four electrons. Exchange-correlation potentials, V{sub xc}(r), and functional derivatives of the noninteracting kinetic energy, δK[ρ]/δρ(r), obtained from these densities were used to examine their discontinuities as the number of electrons N increases across integer boundaries for N = 1, N = 2, and N = 3. These numerical results are consistent with conclusions that the discontinuities are characterized by a jump in the chemical potential while the shape of V{sub xc}(r) varies continuously as an integer boundary is crossed. The discontinuity of the V{sub xc}(r) is positive, depends on the ionization potential, electron affinity, and orbital energy differences, and the discontinuity in δK[ρ]/δρ(r) depends on the difference between the energies of the highest occupied and lowest unoccupied orbitals. The noninteracting kinetic energy and the exchange correlation energy have been computed for integer and noninteger values of N between 1 and 4.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
A pure-sampling quantum Monte Carlo algorithm
The objective of pure-sampling quantum Monte Carlo is to calculate physical properties that are independent of the importance sampling function being employed in the calculation, save for the mismatch of its nodal hypersurface with that of the exact wave function. To achieve this objective, we report a pure-sampling algorithm that combines features of forward walking methods of pure-sampling and reptation quantum Monte Carlo (RQMC). The new algorithm accurately samples properties from the mixed and pure distributions simultaneously in runs performed at a single set of time-steps, over which extrapolation to zero time-step is performed. In a detailed comparison, we found RQMC to be less efficient. It requires different sets of time-steps to accurately determine the energy and other properties, such as the dipole moment. We implement our algorithm by systematically increasing an algorithmic parameter until the properties converge to statistically equivalent values. As a proof in principle, we calculated the fixed-node energy, static α polarizability, and other one-electron expectation values for the ground-states of LiH and water molecules. These quantities are free from importance sampling bias, population control bias, time-step bias, extrapolation-model bias, and the finite-field approximation. We found excellent agreement with the accepted values for the energy and a variety of other properties for those systems.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
ALmost EXact boundary conditions for transient Schrödinger–Poisson system
For the Schrödinger–Poisson system, we propose an ALmost EXact (ALEX) boundary condition to treat accurately the numerical boundaries. Being local in both space and time, the ALEX boundary conditions are demonstrated to be effective in suppressing spurious numerical reflections. Together with the Crank–Nicolson scheme, we simulate a resonant tunneling diode. The algorithm produces numerical results in excellent agreement with those in Mennemann et al. [1], yet at a much reduced complexity. Primary peaks in wave function profile appear as a consequence of quantum resonance, and should be considered in selecting the cut-off wave number for numerical simulations.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Quantum Newtonian cosmology and the biconfluent Heun functions
We obtain the exact solution of the Schrödinger equation for a particle (galaxy) moving in a Newtonian universe with a cosmological constant, which is given in terms of the biconfluent Heun functions. The first six Heun polynomials of the biconfluent function are written explicitly. The energy spectrum which resembles the one corresponding to the isotropic harmonic oscillator is also obtained. The wave functions as well as the energy levels codify the role played by the cosmological constant.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Some remarks on the sternheimer potential
Generalization of Sternheimer potential to include wave functions with spin
POTENTIAL ENERGY↗
Perturbation theory of the constrained variational method in molecular quantum mechanics.
Energy minimization for trial wave function constrained to give expectation value of operator
PERTURBATION THEORY↗
Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance
By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Effect of spin-orbit coupling on the structure of the electron ground state in silicon nanocrystals
The effect of spin-orbit interaction on the structure of the ground state in the conduction band of spherical silicon nanocrystals is theoretically studied using the envelope-function approximation and the k {center_dot} p method. It is shown that the arising weak spin-orbit coupling of the conduction- and valence bands leads to specific asymmetric hybridization of the s- and p-type envelope functions with opposite spin orientations caused by the anisotropy of spin mixing in the silicon conduction band. As a result, the wave functions of the ground-state transform which is accompanied by an insignificant decrease in its energy. In this case, the spin-mixing parameter in nanocrystals depends strongly on their size due to the quantum-confinement effect.
36 MATERIALS SCIENCE↗
An Ab Initio Multiple Cloning Method for Non-Adiabatic Excited-State Molecular Dynamics in NWChem
The recently developed ab-initio multiple cloning (AIMC) approach based on the multicongurational Ehrenfest (MCE) method provides a powerful and accurate way of describing the excited-state dynamics of molecular systems. The AIMC method is a controlled approximation to non-adiabatic dynamics with a particular strength in the proper description of decoherence effects due to branching of vibrational wavepackets at a level crossing. In this work, we report a new implementation of the AIMC algorithm in the open source NWChem computational chemistry program. The framework combines linear-response time-dependent density functional theory with Ehrenfest mean-eld theory to determine the equations of motion for classical trajectories. The multi-dimensional wave function is decomposed into a superposition of Gaussian coherent states guided by Ehrenfest trajectories (i.e., MCE approach), which can clone with fully quantum mechanical amplitudes and phases. By using an efficient time-derivative based non-adiabatic coupling approach within the AIMC method, all observable are calculated on-the-y in the non-adiabatic molecular dynamics process. As a representative example, we apply our implementation to study the ultrafast photoinduced electronic and vibrational energy transfer in a pyridine molecule. The effects of the cloning procedure on electronic and vibrational coherence, relaxation and unidirectional energy transfer are discussed. This new AIMC implementation provides a high-level non-adiabatic molecular dynamics framework for simulating photoexcited dynamics in complex molecular systems and experimentally relevant ultrafast spectroscopic probes, such as nonlinear coherent optical and X-ray signals.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
GENERALIZATION OF THE THEORY OF PARTIAL SEPARATION OF VARIABLES TO INCLUDE THE CASE OF MULTIVALENT ATOMS
Wave function of multivalent atom from linear determinants formulated from laplace theorem
WAVE FUNCTION↗
Collective motion in prolate γ-rigid nuclei within minimal length concept via a quantum perturbation method
Based on the minimal length concept, inspired by Heisenberg algebra, a closed analytical formula is derived for the energy spectrum of the prolate γ-rigid Bohr–Mottelson Hamiltonian of nuclei, within a quantum perturbation method (QPM), by considering a scaled Davidson potential in β shape variable. In the resulting solution, called X(3)-D-ML, the ground state and the first β-band are all studied as a function of the free parameters. The fact of introducing the minimal length concept with a QPM makes the model very flexible and a powerful approach to describe nuclear collective excitations of a variety of vibrational-like nuclei. The introduction of scaling parameters in the Davidson potential enables us to get a physical minimum of this latter in comparison with previous works. The analysis of the corrected wave function, as well as the probability density distribution, shows that the minimal length parameter has a physical upper bound limit.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Resolution-of-identity stochastic time-dependent configuration interaction for dissipative electron dynamics in strong fields
In this contribution, we introduce a method for simulating dissipative, ultrafast many-electron dynamics in intense laser fields. The method is based on the norm-conserving stochastic unraveling of the dissipative Liouville-von Neumann equation in its Lindblad form. The N-electron wave functions sampling the density matrix are represented in the basis of singly excited configuration state functions. The interaction with an external laser field is treated variationally and the response of the electronic density is included to all orders in this basis. The coupling to an external environment is included via relaxation operators inducing transition between the configuration state functions. Single electron ionization is represented by irreversible transition operators from the ionizing states to an auxiliary continuum state. The method finds its efficiency in the representation of the operators in the interaction picture, where the resolution-of-identity is used to reduce the size of the Hamiltonian eigenstate basis. The zeroth-order eigenstates can be obtained either at the configuration interaction singles level or from a time-dependent density functional theory reference calculation. The latter offers an alternative to explicitly time-dependent density functional theory which has the advantage of remaining strictly valid for strong field excitations while improving the description of the correlation as compared to configuration interaction singles. The method is tested on a well-characterized toy system, the excitation of the low-lying charge transfer state in LiCN.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Some Remarks on the Use of the Variational Principle for the Second Order Energy
Approximate ground state wave function used in variational principle for second order energy
WAVE FUNCTION↗
Calculations of atomic energy levels
Vector model for wave functions of atomic energy levels
STATE VECTOR↗
Atomistic k ⋅ p theory
Pseudopotentials, tight-binding models, and k ⋅ p theory have stood for many years as the standard techniques for computing electronic states in crystalline solids. Here, we present the first new method in decades, which we call atomistic k ⋅ p theory. In its usual formulation, k ⋅ p theory has the advantage of depending on parameters that are directly related to experimentally measured quantities, however, it is insensitive to the locations of individual atoms. We construct an atomistic k ⋅ p theory by defining envelope functions on a grid matching the crystal lattice. The model parameters are matrix elements which are obtained from experimental results or ab initio wave functions in a simple way. This is in contrast to the other atomistic approaches in which parameters are fit to reproduce a desired dispersion and are not expressible in terms of fundamental quantities. This fitting is often very difficult. We illustrate our method by constructing a four-band atomistic model for a diamond/zincblende crystal and show that it is equivalent to the sp{sup 3} tight-binding model. We can thus directly derive the parameters in the sp{sup 3} tight-binding model from experimental data. We then take the atomistic limit of the widely used eight-band Kane model and compute the band structures for all III–V semiconductors not containing nitrogen or boron using parameters fit to experimental data. Our new approach extends k ⋅ p theory to problems in which atomistic precision is required, such as impurities, alloys, polytypes, and interfaces. It also provides a new approach to multiscale modeling by allowing continuum and atomistic k ⋅ p models to be combined in the same system.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Asymptotics of the Resonant Tunneling of High-Energy Electrons in Two-Dimensional Quantum Waveguides of Variable Cross-Section
A waveguide occupies a strip in ℝ{sup 2} having two identical narrows of small diameter ε. An electron wave function satisfies the Helmholtz equation with the homogeneous Dirichlet boundary condition. The energy of electrons may be rather high, i.e., any (fixed) number of waves can propagate in the strip far from the narrows. As ε → 0, a neighborhood of a narrow is assumed to transform into a neighborhood of the common vertex of two vertical angles. The part of the waveguide between the narrows as ε = 0 is called the resonator. An asymptotics of the transmission coefficient is obtained in the waveguide as ε → 0. Near a degenerate eigenvalue of the resonator, the leading term of the asymptotics has two sharp peaks. Positions and shapes of the resonant peaks are described.
97 MATHEMATICS AND COMPUTING↗
Multinode Multi-GPU Two-Electron Integrals: Code Generation Using the Regent Language
The computation of two-electron repulsion integrals (ERIs) is often the most expensive step of integral-direct self-consistent field methods. Formally it scales as O(N 4 ), where N is the number of Gaussian basis functions used to represent the molecular wave function. In practice, this scaling can be reduced to O(N 2 ) or less by neglecting small integrals with screening methods. The contributions of the ERIs to the Fock matrix are of Coulomb (J) and exchange (K) type and require separate algorithms to compute matrix elements efficiently. We previously implemented highly efficient GPU-accelerated J-matrix and K-matrix algorithms in the electronic structure code TeraChem. Although these implementations supported the use of multiple GPUs on a node, they did not support the use of multiple nodes. This presents a key bottleneck to cutting-edge ab initio simulations of large systems, e.g., excited state dynamics of photoactive proteins. We present our implementation of multinode multi-GPU J- and K-matrix algorithms in TeraChem using the Regent programming language. Regent directly supports distributed computation in a task-based model and can generate code for a variety of architectures, including NVIDIA GPUs. We demonstrate multinode scaling up to 45 GPUs (3 nodes) and benchmark against hand-coded TeraChem integral code. Finally, we also outline our metaprogrammed Regent implementation, which enables flexible code generation for integrals of different angular momenta.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗