Search NASA⌕ Search

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.

At least 451 records · Page 25

Natural geometric representation for electron local observables

An existence of the quartic identities for the electron local observables that define orthogonality relations for the 3D quantities quadratic in the electron observables is found. It is shown that the joint solution of the quartic and bilinear identities for the electron observables defines a unique natural representation of the observables. In the natural representation the vector type electron local observables have well-defined fixed positions with respect to a local 3D orthogonal reference frame. It is shown that the natural representation of the electron local observables can be defined in six different forms depending on a choice of the orthogonal unit vectors. The natural representation is used to determine the functional dependence of the electron wave functions on the local observables valid for any shape of the electron wave packet. -- Highlights: •Quartic identities that define the orthogonality relations for the electron local observables are found. •Joint solution of quartic and bilinear identities defines a unique natural representation of the electron local observables. •Functional dependence of the electron wave functions on the electron local observables is determined.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Gamow shell model description of proton scattering on 18 Ne

Background: The structure of weakly bound/unbound nuclei close to particle drip lines is different from that around the valley of beta stability. A comprehensive description of these systems goes beyond the standard shell model (SM) and demands an open quantum system description of the nuclear many-body system. Purpose: For that purpose, we are using the Gamow shell model (GSM) in this study, which provides a fully microscopic description of bound and unbound nuclear states, nuclear decays, and reactions. Here, we formulate the GSM in coupled-channel (GSM-CC) representation to describe low-energy elastic and inelastic scattering of protons on $\mathrm{^{18}Ne}$. Method: The GSM-CC formalism is applied to a translationally invariant Hamiltonian with an effective finite-range two-body interaction. We discuss in detail the GSM-CC formalism in coordinate space and give the description of the novel equivalent potential method for solving the GSM-CC system of integrodifferential equations. This method is then applied for the description of ($p, p'$) reaction cross-sections. The reactions channels are built by GSM wave functions for the ground state $\mathrm{0^+}$ and the first excited $\mathrm{2^+}$ of $\mathrm{^{18}Ne}$ and a proton wave function expanded in different partial waves. The completeness of this basis is verified by comparing GSM and GSM-CC energies of low-energy resonant states in $\mathrm{^{19}Na}$. The differences between the two calculations provide a measure of the missing configurations in the GSM-CC calculation of low-energy states of $\mathrm{^{19}Na}$ due to the restriction on the number of excited states of $\mathrm{^{18}Ne}$. Results: We present the first application of the GSM-CC formalism for the calculation of excited states of $\mathrm{^{18}Ne}$ and $\mathrm{^{19}Na}$, the excitation function, and the elastic/inelastic differential cross-sections in the $\mathrm{^{18}Ne}$ ($p,p'$) reaction at different energies. This is the first unified description of the spectra and reaction cross-sections in the GSM formalism. The method is shown to be both feasible and accurate. The approximate equivalence of GSM and GSM-CC in describing spectra of $\mathrm{^{19}Na}$ has been demonstrated numerically. Conclusions: The GSM in the coupled-channel representation opens a possibility for the unified description of low-energy nuclear structure and reactions using the same Hamiltonian. While both GSM and GSM-CC can describe energies, widths, and wave functions of the many-body states, the GSM-CC can in addition yield reaction cross-sections. The combined application of GSM and GSM-CC to describe energies of resonant states allows to test the exactitude of calculated cross-sections for a given many-body Hamiltonian.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Time-dependent restricted-active-space self-consistent-field singles method for many-electron dynamics

The time-dependent restricted-active-space self-consistent-field singles (TD-RASSCF-S) method is presented for investigating TD many-electron dynamics in atoms and molecules. Adopting the SCF notion from the muticonfigurational TD Hartree-Fock (MCTDHF) method and the RAS scheme (single-orbital excitation concept) from the TD configuration-interaction singles (TDCIS) method, the TD-RASSCF-S method can be regarded as a hybrid of them. We prove that, for closed-shell N{sub e}-electron systems, the TD-RASSCF-S wave function can be fully converged using only N{sub e}/2 + 1 ⩽ M ⩽ N{sub e} spatial orbitals. Importantly, based on the TD variational principle, the converged TD-RASSCF-S wave function with M = N{sub e} is more accurate than the TDCIS wave function. The accuracy of the TD-RASSCF-S approach over the TDCIS is illustrated by the calculation of high-order harmonic generation spectra for one-dimensional models of atomic helium, beryllium, and carbon in an intense laser pulse. The electronic dynamics during the process is investigated by analyzing the behavior of electron density and orbitals. The TD-RASSCF-S method is accurate, numerically tractable, and applicable for large systems beyond the capability of the MCTDHF method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Two-photon excitation cross section in light and intermediate atoms in frozen-core LS-coupling approximation

Using the method of explicit summation over the intermediate states two-photon absorption cross sections in light and intermediate atoms based on the simplistic frozen-core approximation and LS coupling have been formulated. Formulas for the cross section in terms of integrals over radial wave functions are given. Two selection rules, one exact and one approximate, valid within the stated approximations are derived. The formulas are applied to two-photon absorptions in nitrogen, oxygen, and chlorine. In evaluating the radial integrals, for low-lying levels, the Hartree-Fock wave functions, and for high-lying levels, hydrogenic wave functions obtained by the quantum-defect method have been used. A relationship between the cross section and the oscillator strengths is derived.

Omidvar, K.↗

Comparison of the canonical transformation and energy functional formalisms for ab initio calculations of self-localized polarons

In materials with strong electron-phonon (e-ph) interactions, charge carriers can distort the surrounding lattice and become trapped, forming self-localized (small) polarons. We recently developed an ab initio approach based on canonical transformations to efficiently compute the formation and energetics of small polarons A different approach based on a Landau-Pekar energy functional has been proposed in the recent literature. Here in this work, we analyze and compare these two methods in detail. We show that the small polaron energy is identical in the two formalisms when using the same polaron wave function. We also show that our canonical transformation formalism can predict polaron band structures and can properly treat zero- and finite-temperature lattice vibration effects, although at present using a fixed polaron wave function. Conversely, the energy functional approach can compute the polaron wave function, but as we show here, it neglects lattice vibrations and cannot address polaron self-localization and thermal band narrowing. Taken together, this work relates two different methods developed recently to study polarons from first-principles, highlighting their merits and shortcomings and discussing them both in a unified formalism.

36 MATERIALS SCIENCE↗

Spin-orbit couplings within the equation-of-motion coupled-cluster framework: Theory, implementation, and benchmark calculations

We present a formalism and an implementation for calculating spin-orbit couplings (SOCs) within the EOM-CCSD (equation-of-motion coupled-cluster with single and double substitutions) approach. The following variants of EOM-CCSD are considered: EOM-CCSD for excitation energies (EOM-EE-CCSD), EOM-CCSD with spin-flip (EOM-SF-CCSD), EOM-CCSD for ionization potentials (EOM-IP-CCSD) and electron attachment (EOM-EA-CCSD). We employ a perturbative approach in which the SOCs are computed as matrix elements of the respective part of the Breit-Pauli Hamiltonian using zeroth-order non-relativistic wave functions. We follow the expectation-value approach rather than the response-theory formulation for property calculations. Both the full two-electron treatment and the mean-field approximation (a partial account of the two-electron contributions) have been implemented and benchmarked using several small molecules containing elements up to the fourth row of the periodic table. The benchmark results show the excellent performance of the perturbative treatment and the mean-field approximation. When used with an appropriate basis set, the errors with respect to experiment are below 5% for the considered examples. The findings regarding basis-set requirements are in agreement with previous studies. The impact of different correlation treatment in zeroth-order wave functions is analyzed. Overall, the EOM-IP-CCSD, EOM-EA-CCSD, EOM-EE-CCSD, and EOM-SF-CCSD wave functions yield SOCs that agree well with each other (and with the experimental values when available). Using an EOM-CCSD approach that provides a more balanced description of the target states yields more accurate results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Generalization of the Kohn-Sham system that can represent arbitrary one-electron density matrices

Density functional theory is currently the most widely applied method in electronic structure theory. The Kohn-Sham method, based on a fictitious system of noninteracting particles, is the workhorse of the theory. The particular form of the Kohn-Sham wave function admits only idempotent one-electron density matrices whereas wave functions of correlated electrons in post-Hartree-Fock methods invariably have fractional occupation numbers. Here we show that by generalizing the orbital concept and introducing a suitable dot product as well as a probability density, a noninteracting system can be chosen that can represent the one-electron density matrix of any system, even one with fractional occupation numbers. This fictitious system ensures that the exact electron density is accessible within density functional theory. It can also serve as the basis for reduced density matrix functional theory. Moreover, to aid the analysis of the results the orbitals may be assigned energies from a mean-field Hamiltonian. This produces energy levels that are akin to Hartree-Fock orbital energies such that conventional analyses based on Koopmans' theorem are available. Lastly, this system is convenient in formalisms that depend on creation and annihilation operators as they are trivially applied to single-determinant wave functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Finite-volume pionless effective field theory for few-nucleon systems with differentiable programming

Finite-volume pionless effective field theory provides an efficient framework for the extrapolation of nuclear spectra and matrix elements calculated at finite volume in lattice QCD to infinite volume, and to nuclei with larger atomic number. In this work, it is demonstrated how this framework may be implemented via a set of correlated Gaussian wave functions optimized using differentiable programming and via solution of a generalized eigenvalue problem. This approach is shown to be significantly more efficient than a stochastic implementation of the variational method based on the same form of correlated Gaussian wave functions, yielding comparably accurate representations of the ground-state wave functions with an order of magnitude fewer terms. The efficiency of representation allows such calculations to be extended to larger systems than in previous work. Further, the method is demonstrated through calculations of the binding energies of nuclei with atomic number A ϵ {2,3,4} in finite volume, matched to lattice QCD calculations at quark masses corresponding to m π = 806 MeV, and infinite-volume effective field theory calculations of A ϵ {2,3,4,5,6} systems based on this matching.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Revisiting the Hole Size in Double Helical DNA with Localized Orbital Scaling Corrections

The extent of electronic wave function delocalization for the charge carrier (electron or hole) in double helical DNA plays an important role in determining the DNA charge transfer mechanism and kinetics. The size of the charge carrier’s wave function delocalization is regulated by the solvation induced localization and the quantum delocalization among the π stacked base pairs at any instant of time. Using a newly developed localized orbital scaling correction (LOSC) density functional theory method, we accurately characterized the quantum delocalization of the hole wave function in double helical B-DNA. This approach can be used to diagnose the extent of delocalization in fluctuating DNA structures. Our studies indicate that the hole state tends to delocalize among 4 guanine–cytosine (GC) base pairs and among 3 adenine–thymine (AT) base pairs when these adjacent bases fluctuate into degeneracy. The relatively small delocalization in AT base pairs is caused by the weaker π–π interaction. This extent of delocalization has significant implications for assessing the role of coherent, incoherent, or flickering coherent carrier transport in DNA.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The nature of the chemical bond and the role of non-dynamical and dynamical correlation in Be2

In the spin-coupled generalized valence bond (SCGVB) description of Be2, there is a pair of electrons in highly overlapping “inner” orbitals corresponding to a traditional σ bond, but this bond is compromised by Pauli repulsion arising from its overlap with a second “outer” pair. The presence of this outer pair of electrons leads to a repulsive potential energy curve at long range and a bound, but metastable molecule at short range. To obtain further insights into the nature of the bond in Be2, we determined the non-dynamical and dynamical correlation contributions to the potential energy curve of Be2 using four different choices for the zero-order wave function: Restricted Hartree–Fock (RHF), SCGVB, valence-CASSCF(4,4), and valence-CASSCF(4,8). The SCGVB and valence-CASSCF(4,4) wave functions yield similar breakdowns of the total correlation energy, with non-dynamical correlation being the more important contribution. For the RHF and valence-CASSCF(4,8) wave functions, dynamical correlation is critical, without which the potential energy curve is purely repulsive. High accuracy calculations on the HBen−1Be–BeBen−1H molecule as a function of n (n = 1–6) suggest that the intrinsic strength of a Be–Be σ bond uncompromised by Pauli repulsion is on the order of 62–63 kcal/mol, and its length is 2.13–2.14 Å, ∼60 kcal/mol stronger and ∼0.35 Å shorter than in Be2.

Xu, Lu. T. (ORCID:0000000182064902)↗

Communication: Active space decomposition with multiple sites: Density matrix renormalization group algorithm

We extend the active space decomposition method, recently developed by us, to more than two active sites using the density matrix renormalization group algorithm. The fragment wave functions are described by complete or restricted active-space wave functions. Numerical results are shown on a benzene pentamer and a perylene diimide trimer. It is found that the truncation errors in our method decrease almost exponentially with respect to the number of renormalization states M, allowing for numerically exact calculations (to a few μE{sub h} or less) with M = 128 in both cases. This rapid convergence is because the renormalization steps are used only for the interfragment electron correlation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transition probabilities and radiative lifetimes of levels in F I

The electric dipole transition probabilities and the lifetimes of excited levels have been calculated using the weakest bound electron potential model theory (WBEPMT) and the quantum defect orbital theory (QDOT) in atomic fluorine. In the calculations, many of transition arrays included both multiplet and fine-structure transitions are considered. We employed Numerical Coulomb Approximation (NCA) wave functions and numerical non-relativistic Hartree-Fock (NRHF) wave functions for expectation values of radii in determination of parameters. The necessary energy values have been taken from experimental energy data in the literature. The calculated transition probabilities and lifetimes have been compared with available theoretical and experimental results. A good agreement with results in literature has been obtained. Moreover, some transition probability and the lifetime values not existing in the literature for some highly excited levels have been obtained using these methods.

74 ATOMIC AND MOLECULAR PHYSICS↗

Stochastic approach to correlations beyond the mean field with the Skyrme interaction

Large-scale calculation based on the multi-configuration Skyrme density functional theory is performed for the light N=Z even-even nucleus, {sup 12}C. Stochastic procedures and the imaginary-time evolution are utilized to prepare many Slater determinants. Each state is projected on eigenstates of parity and angular momentum. Then, performing the configuration mixing calculation with the Skyrme Hamiltonian, we obtain low-lying energy-eigenstates and their explicit wave functions. The generated wave functions are completely free from any assumption and symmetry restriction. Excitation spectra and transition probabilities are well reproduced, not only for the ground-state band, but for negative-parity excited states and the Hoyle state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Diffusion Monte Carlo approaches for studying nuclear quantum effects in fluxional molecules

Abstract Diffusion quantum Monte Carlo (DMC) provides a powerful approach for obtaining the ground state energy and wave function of molecules, ions, and molecular clusters. The approach is uniquely well suited for studies of fluxional molecules, which undergo large amplitude vibrational motions even in their ground state. In contrast to the electronic structure problem, where the wave function must be antisymmetric with respect to exchange of any pair of electrons, the wave function for the ground vibrational state is nodeless. This greatly simplifies the application of DMC for vibrational problems. Because there is not a single potential function that can be used to describe the intramolecular and intermolecular interactions in all molecular systems, most methods that are used to describe nuclear quantum effects rely on a carefully chosen zero‐order description of the molecular vibrations. In contrast, DMC calculations can be performed in Cartesian coordinates, making the DMC algorithm easily transferable between different chemical systems. In this contribution, the theory that underlies DMC will be discussed along with important considerations for performing DMC calculations. Extensions for evaluating vibrationally excited states and molecular properties are also discussed. Insights that can be obtained from DMC calculations are illustrated in the context of the protonated water clusters. This article is categorized under: Molecular and Statistical Mechanics > Molecular Dynamics and Monte‐Carlo Methods Theoretical and Physical Chemistry > Spectroscopy

Chemistry↗

Quantum dynamics of electronic transitions with Gauss-Hermite wave packets

A new methodology based on the superposition of time-dependent Gauss-Hermite wave packets is developed to describe the wave function of a system in which several interacting electronic states are coupled to a bath of harmonic oscillators. The equations of motion for the wave function parameters are obtained by employing the Dirac-Frenkel time-dependent variational principle. The methodology is applied to study the quantum dynamical behaviour of model systems with two interacting electronic states characterized by a relatively large reorganization energy and a range of energy biases. The favourable scaling properties make it a promising tool for the study of the dynamics of chemico-physical processes in molecular systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A T Matrix Method Based upon Scalar Basis Functions

A surface integral formulation is developed for the T matrix of a homogenous and isotropic particle of arbitrary shape, which employs scalar basis functions represented by the translation matrix elements of the vector spherical wave functions. The formulation begins with the volume integral equation for scattering by the particle, which is transformed so that the vector and dyadic components in the equation are replaced with associated dipole and multipole level scalar harmonic wave functions. The approach leads to a volume integral formulation for the T matrix, which can be extended, by use of Green's identities, to the surface integral formulation. The result is shown to be equivalent to the traditional surface integral formulas based on the VSWF basis.

wave functions↗

A new look at the position operator in quantum theory

The postulate that coordinate and momentum representations are related to each other by the Fourier transform has been accepted from the beginning of quantum theory by analogy with classical electrodynamics. As a consequence, an inevitable effect in standard theory is the wave packet spreading (WPS) of the photon coordinate wave function in directions perpendicular to the photon momentum. This leads to the following paradoxes: if the major part of photons emitted by stars are in wave packet states (what is the most probable scenario) then we should see not separate stars but only an almost continuous background from all stars; no anisotropy of the CMB radiation should be observable; data on gamma-ray bursts, signals from directional radio antennas (in particular, in experiments on Shapiro delay) and signals from pulsars show no signs of WPS. In addition, a problem arises why there are no signs of WPS for protons in the LHC ring. We argue that the above postulate is based neither on strong theoretical arguments nor on experimental data and propose a new consistent definition of the position operator. Then WPS in directions perpendicular to the particle momentum is absent and the paradoxes are resolved. Different components of the new position operator do not commute with each other and, as a consequence, there is no wave function in coordinate representation. Implications of the results for entanglement, quantum locality and the problem of time in quantum theory are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Huzinaga projection embedding for efficient and accurate energies of systems with localized spin-densities

We demonstrate the accuracy and efficiency of the restricted open-shell and unrestricted formulation of the absolutely localized Huzinaga projection operator embedding method. Restricted open-shell and unrestricted Huzinaga projection embedding in the full system basis is formally exact to restricted open-shell and unrestricted Kohn–Sham density functional theory, respectively. By utilizing the absolutely localized basis, we significantly improve the efficiency of the method while maintaining high accuracy. Furthermore, the absolutely localized basis allows for high accuracy open-shell wave function methods to be embedded into a closed-shell density functional theory environment. The open-shell embedding method is shown to calculate electronic energies of a variety of systems to within 1 kcal/mol accuracy of the full system wave function result. For certain highly localized reactions, such as spin transition energies on transition metals, we find that very few atoms are necessary to include in the wave function region in order to achieve the desired accuracy. This extension further broadens the applicability of our absolutely localized Huzinaga level-shift projection operator method to include open-shell species. Here, we apply our method to several representative examples, such as spin splitting energies, catalysis on transition metals, and radical reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗