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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Exchange correlation potentials from full configuration interaction in a Slater orbital basis

Ryabinkin–Kohut–Staroverov (RKS) theory builds a bridge between wave function theory and density functional theory by using quantities from the former to produce accurate exchange-correlation potentials needed by the latter. In this work, the RKS method is developed and tested alongside Slater atomic orbital basis functions for the first time. To evaluate this approach, full configuration interaction computations in the Slater orbital basis are employed to give quality input to RKS, allowing full correlation to be present along with correct nuclei cusps and asymptotic decay of the wavefunction. SlaterRKS is shown to be an efficient algorithm to arrive at exchange-correlation potentials without unphysical artifacts in moderately-sized basis sets. Furthermore, enforcement of the nuclear cusp conditions will be shown to be vital for the success of the Slater-basis RKS method. As a result, examples of weakly and strongly correlated molecular systems will demonstrate the main features of SlaterRKS.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Tuning perovskite nanocrystal superlattices for superradiance in the presence of disorder

The cooperative emission of interacting nanocrystals is an exciting topic fueled by recent reports of superfluorescence and superradiance in assemblies of perovskite nanocubes. Several studies estimated that coherent coupling is localized to a small fraction of nanocrystals (10 -7 –10 -3 ) within the assembly, raising questions about the origins of localization and ways to overcome it. In this work, we examine single-excitation superradiance by calculating radiative decays and the distribution of superradiant wave function in two-dimensional CsPbBr 3 nanocube superlattices. The calculations reveal that the energy disorder caused by size distribution and large interparticle separations reduces radiative coupling and leads to the excitation localization, with the energy disorder being the dominant factor. The single-excitation model clearly predicts that, in the pursuit of cooperative effects, having identical nanocubes in the superlattice is more important than achieving a perfect spatial order. The monolayers of large CsPbBr 3 nanocubes (L NC = 10–20 nm) are proposed as model systems for experimental tests of superradiance under conditions of non-negligible size dispersion, while small nanocubes (L NC = 5–10 nm) are preferred for realizing the Dicke state under ideal conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Design of AlGaN-Zn(Si,Ge)N2 quantum wells for high-efficiency ultraviolet light emitters

The effect of inserting a nm-scale layer of Zn(Si,Ge)N2 into an AlGaN quantum well structure designed for light emission in the wavelength range from 255 to 305 nm is investigated here. The enhanced confinement of the hole within the quantum well results in an enhancement of the overlap of the hole and electron wave functions, resulting in an enhancement of the radiative recombination rate. In this theoretical calculation, for emission at a 270 nm wavelength, the enhancement in the wavefunction overlap can reach a factor of 7 when compared to an AlGaN quantum well device specifically engineered for optimal emission at the identical wavelength. Increases of almost an order of magnitude in both the peak spontaneous emission intensity and the radiative recombination rate are predicted. The peak emission wavelength can be tuned from 255 to 305 nm by adjusting the width and/or the composition of the inserted layer. The proposed structures provide a route to higher efficiency ultraviolet practical light emitting diodes and lasers.

Physics↗
DOE OSTI2024

Hartree–Fock–Bogoliubov theory for number-parity-violating fermionic Hamiltonians

It is usually asserted that physical Hamiltonians for fermions must contain an even number of fermion operators. This is indeed true in electronic structure theory. However, when the Jordan–Wigner (JW) transformation is used to map physical spin Hamiltonians to Hamiltonians of spinless fermions, terms that contain an odd number of fermion operators may appear. The resulting fermionic Hamiltonian thus does not have number parity symmetry and requires wave functions that do not have this symmetry either. In this work, we discuss the extension of standard Hartree–Fock–Bogoliubov (HFB) theory to the number-parity-nonconserving case. These ideas had appeared in the literature before but, perhaps for lack of practical applications, had, to the best of our knowledge, never been employed. We here present a useful application for this more general HFB theory based on coherent states of the SO(2M + 1) Lie group, where M is the number of orbitals. Here, we also show how using these unusual mean-field states can provide significant improvements when studying the JW transformation of chemically relevant spin Hamiltonians.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Method-independent cusps for atomic orbitals in quantum Monte Carlo

Here, we present an approach for augmenting Gaussian atomic orbitals with correct nuclear cusps. Like the atomic orbital basis set itself and unlike previous cusp corrections, this approach is independent of the many-body method used to prepare wave functions for quantum Monte Carlo. Once the basis set and molecular geometry are specified, the cusp-corrected atomic orbitals are uniquely specified, regardless of which density functionals, quantum chemistry methods, or subsequent variational Monte Carlo optimizations are employed. We analyze the statistical improvement offered by these cusps in a number of molecules and find them to offer similar advantages as molecular-orbital-based approaches while remaining independent of the choice of many-body method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Symmetry-projected spin-AGP methods applied to spin systems

Symmetry-projected wave function methods capture static correlation by breaking and restoring the symmetries of a system. In this article, we present the symmetry-projected spin antisymmetrized geminal power (spin-AGP) state projected onto space group symmetry as well as complex conjugation, spin-flip, and time-reversal symmetries. The method is benchmarked on the 1D XXZ model and the 2D J 1 − J 2 model with square and triangular lattices. Our results indicate that symmetry projection methods provide a powerful tool for frustrated spin systems.

Antisymmetrized geminal power↗
DOE OSTI2026

Degenerate coupled-cluster theory

A size-extensive, converging, black-box, ab initio coupled-cluster (ΔCC) ansatz is introduced that computes the energies and wave functions of states from any degenerate or nondegenerate Slater-determinant references with any numbers of α- and β-spin electrons, any patterns of orbital occupancy, any spin multiplicities, and any spatial symmetries. For a nondegenerate reference, it reduces to the single-reference coupled-cluster ansatz. For a degenerate multireference, it is a natural coupled-cluster extension of degenerate Møller–Plesset perturbation (ΔMP) theory. For ionized and electron-attached references, it is a coupled-cluster Green’s function, although the present theory is convergent toward the full-configuration-interaction limits, while the Feynman–Dyson many-body Green’s function (MBGF) theory generally is not. Its single-excitation instance is a projection Hartree–Fock theory as per the Thouless theorem, which may be useful for core ionizations, high-spin states, and possibly electron affinities. Additionally, a new multireference coupled-cluster theory for a general model space is developed. This quasidegenerate coupled-cluster (QCC) theory is exactly converging, but not black-box, and intended for strong correlation. Determinant-based, general-order algorithms of ΔCC and QCC theories are implemented and compared with configuration-interaction (CI) and equation-of-motion coupled-cluster (EOM-CC) theories through octuple excitations and with ΔMP and MBGF theories up to the nineteenth order. An algebraic, optimal-scaling algorithm of the ΔCC theory is computer-synthesized at the levels of single excitations (ΔCCS) and of single and double excitations (ΔCCSD). As a result, the order of performance is QCC ≈ ΔCC > EOM-CC > CI at the same order or QCC ≈ ΔCC > ΔMP > MBGF at the same cost scaling.

Hirata, So [University of Illinois at Urbana-Champ↗
DOE OSTI2022

The cosine-sine decomposition with different-orbitals-for-different-spins determinants

We report the spin decomposition of a spin-contaminated different-orbitals-for-different-spins (DODS) wave function is formulated in terms of the Krylov space of the $\hat{S}^2$ ||operator. The number of determinants that contribute to the various projected spin eigenfunctions vertical | $\psi ; S, M \rangle$ depends on the sparsity of the orbital overlap matrix X between the α and β spatial orbitals. The cosine-sine decomposition (CSD) procedure may be applied to this overlap matrix, and the resulting redundant orbital transformations may be applied to the α and β spatial orbitals. This produces a sparse X' matrix in the transformed basis in which each row or column has either one or two nonzero elements. This sparse X' matrix simplifies the spin-decomposition procedure in three ways:1) it reduces the number of contributing determinants within the Krylov basis functions and within the projected spin functions, 2) it reduces the effective orbital dimension through elimination of the frozen core and frozen virtual orbitals, and 3) it simplifies the spin-decomposition procedure in both the original and transformed bases by limiting the Krylov space dimension. These simplifications all reduce the computational effort for the spin-decomposition process. This procedure is implemented within a string-based DODS determinant formulation and applied to spin projection of unrestricted Hartree-Fock determinants.

unrestricted Hartree-Fock↗
DOE OSTI2021

Non-relativistic geometry and the equivalence principle

Abstract We describe a geometric and symmetry-based formulation of the equivalence principle in non-relativistic physics. It applies both on the classical and quantum levels and states that the Newtonian potential can be eliminated in favor of a curved and time-dependent spatial metric. It is this requirement that forces the gravitational mass to be equal to the inertial mass. We identify the symmetry responsible for the equivalence principle as the remnant of time-reparameterization symmetry of the relativistic theory. We also clarify the transformation properties of the Schrödinger wave-function under arbitrary frame changes.

Astronomy & Astrophysics↗
DOE OSTI2020

One-center close-coupling approach to two-center rearrangement collisions

Calculations of ionization and electron-capture cross sections in ion-atom collisions usually require solving the Schrödinger equation governing the collision system by expanding the total scattering wave function in a basis of target- and projectile-centered pseudostates. This approach leads to the two-center close-coupling equations which in some cases may become ill-conditioned due to the non-orthogonality of the underlying combined basis. Here we develop a technique which allows to accurately extract the necessary collision information, including that for the rearrangement channels, from the computationally more convenient one-center close-coupling equations which are built from only target-centered pseudostates. The robustness of the method is demonstrated by considering the proton-hydrogen scattering problem across a wide incident energy range. The developed method is then applied to study proton scattering on multielectron target of lithium. Here, the obtained results for the 2s → 2p excitation and the total electron-capture cross sections are in good agreement with corresponding experimental data. It is concluded that the presented technique could be a simpler alternative when integrated cross sections are required.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

R-matrix calculations for opacities: I. Methodology and computations

Abstract An extended version of the R -matrix methodology is presented for calculation of radiative parameters for improved plasma opacities. Contrast and comparisons with existing methods primarily relying on the distorted wave approximation are discussed to verify accuracy and resolve outstanding issues, particularly with reference to the opacity project (OP). Among the improvements incorporated are: (i) large-scale Breit–Pauli R -matrix calculations for complex atomic systems including fine structure, (ii) convergent close coupling wave function expansions for the ( e + ion) system to compute oscillator strengths and photoionization cross sections, (iii) open and closed shell iron ions of interest in astrophysics and experiments, (iv) a treatment for plasma broadening of autoionizing resonances as function of energy-temperature-density dependent cross sections, (v) a ‘top-up’ procedure to compare convergence with R -matrix calculations for highly excited levels, and (vi) spectroscopic identification of resonances and bound ( e + ion) levels. The present R -matrix monochromatic opacity spectra are fundamentally different from OP and lead to enhanced Rosseland and Planck mean opacities. An outline of the work reported in other papers in this series and those in progress is presented. Based on the present re-examination of the OP work, opacities of heavy elements might require revisions in high temperature-density plasma sources.

Pradhan, A. K. (ORCID:0000000187753643)↗
DOE OSTI2020

S -factor and scattering-parameter extractions from ${}^{3}\mathrm{He}+{}^{4}\mathrm{He}{ \rightarrow }^{7}\mathrm{Be}+\gamma $

Previous studies of the reaction ${}^{3}\mathrm{He}+{}^{4}\mathrm{He}\to {}^{7}\mathrm{Be}+\gamma $ have mainly focused on providing the best central value and error bar for the S factor at solar energies. Experimental measurements of this capture reaction at higher energies, the ${}^{3}\mathrm{He}$–${}^{4}\mathrm{He}$ scattering phase shifts, as well as properties of ${}^{7}\mathrm{Be}$ and its excited state, have been used to constrain the theoretical models employed for this purpose. Here we show that much more information than was previously appreciated can be extracted from angle-integrated capture data alone. We use the next-to-leading-order (NLO) amplitude in an effective field theory (EFT) for ${}^{3}\mathrm{He}+{}^{4}\mathrm{He}\to {}^{7}\mathrm{Be}+\gamma $ to perform the extrapolation. At this order the EFT describes the capture process using an s-wave scattering length and effective range, the asymptotic properties of ${}^{7}\mathrm{Be}$ and its excited state, and short-distance contributions to the E1 capture amplitude. We extract the multi-dimensional posterior of all these parameters via a Bayesian analysis that uses capture data below 2 MeV. We find that properties of the ${}^{7}\mathrm{Be}$ ground and excited states are well constrained. The total S factor $S(0)\,={0.577}_{-0.016}^{+0.015}$ keV b, while the branching ratio for excited- to ground-state capture at zero energy, ${Br}(0)={0.406}_{-0.011}^{+0.013}$, both at 68% degree of belief. This S(0) is broadly consistent with other recent evaluations, and agrees with the previously recommended value $S(0)=0.56\pm 0.03\,\,\mathrm{keV}$ b, but has a smaller error bar. We also find significant constraints on ${}^{3}\mathrm{He}$–${}^{4}\mathrm{He}$ scattering parameters, and we obtain constraints on the angular distribution of capture gamma rays, which is important for interpreting experiments. The path forward for this reaction seems to lie with better measurements of the scattering phase shifts and S(E)'s angular dependence away from zero energy, together with better understanding of the asymptotic normalization coefficients of the ${}^{7}$Be bound states' wave functions. Lastly, data on these could further reduce the uncertainty on S(0).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Rotational bands beyond the Elliott model

Rotational bands are commonplace in the spectra of atomic nuclei. Inspired by early descriptions of these bands by quadrupole deformations of a liquid drop, Elliott constructed discrete nucleon representations of SU(3) from fermionic creation and annihilation operators. Ever since, Elliott's model has been foundational to descriptions of rotation in nuclei. Later work, however, suggested the symplectic extension Sp(3, R) provides a more unified picture. We decompose no-core shell-model nuclear wave functions into symmetry-defined subspaces for several beryllium isotopes, as well as 20 Ne, using the quadratic Casimirs of both Elliott's SU(3) and Sp(3, R). The band structure, delineated by strong B(E2) values, has a more consistent description in Sp(3, R) rather than SU(3). In particular, we confirm previous work finding in some nuclides strongly connected upper and lower bands with the same underlying symplectic structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

The charge and mass symmetry breaking in the KK K¯ system

Abstract In the framework of the Faddeev equations in configuration space, we investigate theK(1460) meson as a resonant state of the KK K ¯ kaonic system. We perform calculations for the particle configurationsK 0 K + K − and K 0 K + K 0 ¯ within two models: theABCmodel, in which all three particles are distinguishable, and theAACmodel when two particles are identical. The models differ in their treatment of the kaon mass difference and the attractive Coulomb force between theK + K − pair. We found that the Coulomb shift adds over 1 MeV to the three-body binding energy. The expected correction to the binding energy due to mass redistribution fromAAtoABis found to be negligible, up to a maximum of 6% of the relative mass correction. At the same time, the symmetry of the wave function is distorted depending on the mass ratio value. We found that the repulsiveKKinteraction plays an essential role in the binding energy of the KK K ¯ system and reports the mass of 1461.8 or 1464.1 MeV for the neutralK 0 (1460) and 1466.5 or 1468.8 MeV for the chargedK + (1460) resonances, respectively, depending on the parameter sets forKKand K K ¯ interactions.

Physics↗
DOE OSTI2021

Full configuration interaction simulations of exchange-coupled donors in silicon using multi-valley effective mass theory

Abstract Donor spins in silicon have achieved record values of coherence times and single-qubit gate fidelities. The next stage of development involves demonstrating high-fidelity two-qubit logic gates, where the most natural coupling is the exchange interaction. To aid the efficient design of scalable donor-based quantum processors, we model the two-electron wave function using a full configuration interaction method within a multi-valley effective mass theory. We exploit the high computational efficiency of our code to investigate the exchange interaction, valley population, and electron densities for two phosphorus donors in a wide range of lattice positions, orientations, and as a function of applied electric fields. The outcomes are visualized with interactive images where donor positions can be swept while watching the valley and orbital components evolve accordingly. Our results provide a physically intuitive and quantitatively accurate understanding of the placement and tuning criteria necessary to achieve high-fidelity two-qubit gates with donors in silicon.

Joecker, Benjamin (ORCID:0000000302635440)↗
DOE OSTI2023

Quantum simulation of bosons with the contracted quantum eigensolver

Abstract Quantum computers are promising tools for simulating many-body quantum systems due to their potential scaling advantage over classical computers. While significant effort has been expended on many-fermion systems, here we simulate a model entangled many-boson system with the contracted quantum eigensolver (CQE). We generalize the CQE to many-boson systems by encoding the bosonic wavefunction on qubits. The CQE provides a compact ansatz for the bosonic wave function whose gradient is proportional to the residual of a contracted Schrödinger equation. We apply the CQE to a bosonic system, where N quantum harmonic oscillators are coupled through a pairwise quadratic repulsion. The model is relevant to the study of coupled vibrations in molecular systems on quantum devices. Results demonstrate the potential efficiency of the CQE in simulating bosonic processes such as molecular vibrations with good accuracy and convergence even in the presence of noise.

Physics↗
DOE OSTI2022

Investigating color screening in proton-nucleus collisions with complex potentials*

Abstract Color screening and parton inelastic scattering modify the heavy-quark antiquark potential in mediums consisting of particles from quantum chromodynamics (QCD), leading to the suppression of quarkonium production in relativistic heavy-ion collisions. Owing to the small charm/anti-charm ( ) pair production number in proton-nucleus ( pA ) collisions, the correlation between different pairs is negligible, which makes the Schrödinger equation viable for tracking the evolution of only one pair. We employ the time-dependent Schrödinger equation with an in-medium potential to study the evolution of charmonium wave functions in a hydrodynamic-like QCD medium produced in pA collisions. We explore different parametrizations of real and imaginary parts of the potential and calculate the nuclear modification factors ( ) of and in TeV energy p -Pb collisions at the Large Hadron Collider (LHC). Comparing strong and weak screening scenarios with experimental data in this approach, we arrive at the conclusion that color screening is weak at temperatures close to the deconfined phase transition. Moreover, the imaginary part of the potential is crucial in describing the experimental data, which is consistent with widely studied semi-classical approaches, where dissociation rates are essential.

Physics↗
DOE OSTI2021

Stability of trapped solutions of a nonlinear Schrödinger equation with a nonlocal nonlinear self-interaction potential

This work focuses on the study of the stability of trapped soliton-like solutions of a (1 + 1)-dimensional nonlinear Schrödinger equation (NLSE) in a nonlocal, nonlinear, self-interaction potential of the form [|Ψ(x,t)| 2 +|Ψ(-x,t)| 2 ] κ where κ is an arbitrary nonlinearity parameter. Although the system with κ = 1 (i.e. fully integrable case) was first reported by Yang (2018 Phys. Rev. E 98 042202), here in the present work, we extend this model to the one in which κ is arbitrary. This allows us to compare the stability properties of the now trapped solutions to previously found solutions of the more usual NLSE with κ ≠ 1 which are moving soliton solutions. We show that there is a simple, one-component, nonlocal Lagrangian and corresponding action governing the dynamics of the system. Using a collective coordinate method derived from the action as well as assuming the validity of Derrick's theorem, we find that these trapped solutions are stable for 0 < κ < 2 and unstable when κ > 2. At the critical value of κ, i.e. κ = 2, the solution can either collapse or blowup linearly in time when q 0 = 0, where q 0 is the center of the initial density ρ(x, t = 0) = ψ*ψ of the solution. For q 0 ≠ 0 the displaced solution collapses. When κ > 2 initial small displacements from the origin also lead to collapse of the wave function. This phenomenon is not seen in the usual NLSE.

collective coordinates↗