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 415 records · Page 23

Reproducibility of calculations on Li species with correlation-consistent basis sets

Here we have found that two different basis set series, both labeled cc-pVxZ, are widely available and widely used for calculation on Li compounds such that calculations that should be the same are different, causing a reproducibility problem. Here we study the magnitude of the disagreement of these nominally-the-same basis sets for wave function and density functional calculations of Li clusters. We find deviations among the basis sets of as large as 22 millihartrees. In most cases one obtains lower energies with the noncanonical choices. Therefore, it is important to specify which version of the basis set is used for calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Convergence of variational Monte Carlo simulation and scale-invariant pre-training

We provide theoretical convergence bounds for the variational Monte Carlo (VMC) method as applied to optimize neural network wave functions for the electronic structure problem. Here, we study both the energy minimization phase and the supervised pre-training phase that is commonly used prior to energy minimization. For the energy minimization phase, the standard algorithm is scale-invariant by design, and we provide a proof of convergence for this algorithm without modifications. The pre-training stage typically does not feature such scale-invariance. We propose using a scale-invariant loss for the pretraining phase and demonstrate empirically that it leads to faster pre-training.

97 MATHEMATICS AND COMPUTING↗

A rigorous non-empirical theoretical analysis of the 2p XPS of NiO: Is it necessary to invoke nonlocal screening?

The role of ligand field effects and many-body effects arising from angular momentum coupling and shake effects for the XPS of NiO, as a representative transition metal oxide, have been investigated using rigorous, non-empirical wave functions for single site cluster models. Here it is shown that important features of the XPS arise from proper treatment of angular momentum coupling and it is unnecessary to invoke non-local screening to explain the experimental XPS. Contrary to the usual understanding, it is shown that inclusion of shake excitations in the many body wavefunctions is responsible for the high BE satellites that are observed. The contribution of covalent mixing of metal and ligand orbitals in the closed shells to screening of core-holes is demonstrated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Precision Møller polarimetry for PREX-2 and CREX

The PREX-2 and CREX experiments in Hall A at Jefferson Lab are precision measurements of parity violating elastic electron scattering from complex nuclei. One requirement was that the incident electron beam polarization, typically , be known with 1% precision. We commissioned and operated a Møller polarimeter on the beam line that exceeds this requirement, achieving a precision of 0.89% for PREX-2, and 0.85% for CREX. The uncertainty is purely systematic, accumulated from several different sources, but dominated by our knowledge of the target polarization. Our analysis also demonstrates the need for accurate atomic wave functions in order to correct for the Levchuk Effect. We describe the details of the polarimeter operation and analysis, as well as (for CREX) a comparison to results from a different polarimeter based on Compton scattering.

43 PARTICLE ACCELERATORS↗

Quantifying uncertainties in α -nucleus reaction dynamics informed from first principles

The ab initio symmetry-adapted no-core shell model is a microscopic many-body method which naturally describes challenging collective and clustering features of atomic nuclei. Wave functions and observables computed with realistic nucleon-nucleon forces in this framework are tied to first principles, and are hence well-suited for rigorous uncertainty quantification. We discuss α-deuteron and α- 12 C cluster potentials informed by symmetry-adapted calculations, and propagate uncertainties in the effective binary cluster method as well as those in the nuclear interaction to reaction observables, namely scattering phase shifts, cross sections, partial widths and resonance energies. Here, we find that the overall uncertainties are dominated by those originating in the underlying nuclear force, speaking to the need for tighter constraints on realistic nucleon-nucleon interactions.

Ab initio↗

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↗

Central vortex steady states and dynamics of Bose–Einstein condensates interacting with a microwave field

Here we study central vortex steady states and dynamics of a two-dimensional (2D), two-component, Gross–Pitaevskii equation (CGPE) system for two pseudo-spinor Bose Einstein condensates (BECs) interacting with an electromagnetic field (microwave) analytically and numerically. For the central vortex steady state at any given winding number S, we prove its existence in a reduced, single component detuning limit when contact-interaction strength β < β b , and nonexistence when β > β b , respectively, where β b is a threshold value, whose value is given in Theorem 3.1 in the paper. We extend the existence and nonexistence result to the general two pseudo-spinor case and prove that a central vortex steady state exists for any given S if β ≤ β b while it does not exist when β > 2β b . We then derive dynamical equations for some observables (expectations of the matter wave function) such as the center-of-mass, position of dispersion, the linear and angular momentum of the two pseudo-spinor CGPEs. Finally, numerical computations are brought in to validate and extend the existence of vortex steady state results to β b ≤ β < 2β b for the two pseudo-spinor case and to explore transient dynamics of the observables.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Basis light-front quantization approach to photon

We solve for the light-front wave functions (LFWFs) of the physical photon from the eigenvectors of the light-front quantum electrodynamics (QED) Hamiltonian with the aim to determine its bare photon and electron-positron Fock components. We then employ the resulting LFWFs to compute the transverse momentum dependent parton distributions (TMDs) and the generalized parton distributions (GPDs) at zero skewness of the photon. The TMDs are found to be in excellent agreement with the lowest-order perturbative results calculated using the electron-positron quantum fluctuation of the photon. The GPDs are also consistent with the perturbative calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Glueball-meson molecules

Experimental searches for pure glueball states have proven challenging and so far yielded no results. This is believed to occur because glueballs mix with the ordinary $q\bar{q}$ states with the same quantum numbers. We will discuss an alternative mechanism, the formation of the glueball-meson molecular states. We will argue that the wave functions of already observed excited meson states may contain a significant part due to such molecular states. We discuss the phenomenology of glueball-meson molecules and comment on a possible charmless component of the XYZ states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Intruder band mixing in an ab initio description of 12 Be

The spectrum of 12 Be exhibits exotic features, e.g., an intruder ground state and shape coexistence, normally associated with the breakdown of a shell closure. While previous phenomenological treatments indicated the ground state has substantial contributions from intruder configurations, it is only with advances in computational abilities and improved interactions that this intruder mixing is observed in ab initio no-core shell model (NCSM) predictions. In this work, we extract electromagnetic observables and symmetry decompositions from the NCSM wave functions to demonstrate that the low-lying positive parity spectrum can be explained in terms of mixing of rotational bands with very different intrinsic structure coexisting within the low-lying spectrum. These observed bands exhibit an approximate SU(3) symmetry and are qualitatively consistent with Elliott model predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quenching of the quantum p -strength in light exotic nuclei

One-nucleon spectroscopic overlaps, their strength or spectroscopic factors (SFs), and nucleon removal cross sections for light nuclei in the mass range A≤12 were evaluated using the fully correlated Quantum Monte Carlo (QMC) wave functions (WFs). Harder (AV18+UX) and softer (NV2+3) bare interactions were used providing consistent results. The SFs were also taken from simple Shell Model (SM) calculations. This structure information was incorporated in the standard three-body Faddeev/Alt-Grassberger-Sandhas reaction formalism to evaluate the (p,pN) cross sections. The results further our understanding of the quenching of the quantum p-shell strength obtained from structure and reactions. We have found the ratios of the total QMC sums of SFs to the SM ones to be uniform and ∼ 3/4. The corresponding ratios of the sums Below Particle Threshold (BPT) of SFs, as well as of total cross sections, deviate strikingly from the uniform trend in some special cases. We find these ratios to be close to unit for 9 Li → 8 Li + n and 9 C → 8 B + p. In contrast, they are strongly reduced for 11 C → 10 C + n and 11 B → 10 Be + p mirror transitions, resulting from the quenching of the strength for low-lying 2 + final state transitions. The QMC theoretical cross section BPT for 11 C(p,pn) is about two times smaller than the experimental data.

(p,pN) reactions↗

Coalescence formation of muonic atoms at RHIC

The discovery of exotic muonic atoms, including antimatter hydrogen muonic atoms and kaon muonic atoms, constitutes a milestone in our ability to make and study new forms of matter. Relativistic heavy-ion collisions provide the only likely condition for production and detection of these exotic atoms. Taking Coulomb correlations into account from the time of the fireball freeze-out, which refers to the cessation of hadronic interactions in the high-temperature, high-density region formed after the collision, until the formation of a stable atom has dramatic consequences on the expected yields of these atoms. When the coalescence model with the assumption of quantum wave function localization is applied to the formation of muonic atoms, we find that the atom yields are about two orders of magnitude higher than previously predicted.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Teaching Nonradiative Transitions with MATLAB and Python

Nonradiative transitions are changes in energy states of atoms, ions, or molecules that do not involve the emission or absorption of photons. Despite their importance in understanding luminescent properties and photochemical reaction mechanisms, nonradiative transitions are rarely given more than a qualitative overview in undergraduate and even graduate physical chemistry curricula. To supplement the coverage of nonradiative transition topics, we provide here a set of active learning exercises to help students develop an intuitive understanding of the factors that determine the rate of nonradiative transitions. Here, we start by outlining the theoretical background through the formulations of the Franck–Condon factor and its relation to the rate of nonradiative transition. We then introduce three teaching modules, with associated MATLAB and Python codes, to explore how (1) the excited state nuclear displacement, (2) the electronic energy gap between excited and ground state and (3) the excited/ground state vibrational mode frequencies affect the magnitude of the Franck–Condon factor and thereby the rate of nonradiative transitions. The wave function overlap plots that accompany all teaching modules provide direct visualization of the effect of input parameters on the magnitude of Franck–Condon overlap integral.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spin–Orbit Matrix Elements for a Combined Spin-Flip and IP/EA approach

We introduce a practical approach for computing the Breit–Pauli spin–orbit matrix elements of multiconfigurational systems with both spin and spatial degeneracies based on our recently developed RAS-nSF-IP/EA method (Houck, S. E.; et al. J. Chem. Theory Comput. 2019, 15, 2278). The spin–orbit matrix elements over all the multiplet components are computed using a single one-particle reduced density matrix as a result of the Wigner–Eckart theorem. A mean field spin–orbit approximation was used to account for the two-electron contributions. Basis set dependence as well as the effect of including additional excitations is presented. The effect of correlating the core and semicore orbitals is also examined. Surprisingly accurate results are obtained for spin–orbit coupling constants, despite the fact that the efficient wave function approximations we research neglect the bulk of dynamical correlation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

N 5 -Scaling Excited-State-Specific Perturbation Theory

In this paper, we show that by working in a basis similar to that of the natural transition orbitals and using a modified zeroth-order Hamiltonian, the cost of a recently introduced perturbative correction to excited-state mean field theory can be reduced from seventh to fifth order in the system size. The (occupied) 2 (virtual) 3 asymptotic scaling matches that of ground-state second-order Møller-Plesset theory but with a significantly higher prefactor because the bottleneck is iterative: it appears in the Krylov-subspace-based solution of the linear equation that yields the first-order wave function. Here, we discuss the details of the modified zeroth-order Hamiltonian we use to reduce the cost and the automatic code generation process we used to derive and verify the cost scaling of the different terms. Overall, we find that our modifications have little impact on the method's accuracy, which remains competitive with singles and doubles equation-of-motion coupled cluster.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Ranked-Orbital Approach to Select Active Spaces for High-Throughput Multireference Computation

The past decade has seen a great increase in the application of high-throughput computation to a variety of important problems in chemistry. However, one area which has been resistant to the high-throughput approach is multireference wave function methods, in large part due to the technicalities of setting up these calculations and in particular the not always intuitive challenge of active space selection. As we look toward a future of applying high-throughput computation to all areas of chemistry, it is important to prepare these methods for large-scale automation. Here, we propose a ranked-orbital approach to select active spaces with the goal of standardizing multireference methods for high-throughput computation. This method allows for the meaningful comparison of different active space selection schemes and orbital localizations, and we demonstrate the utility of this approach across 1120 multireference calculations for the excitation energies of small molecules. Our results reveal that it is helpful to distinguish the method used to generate orbitals from the method of ranking orbitals in terms of importance for the active space. Additionally, we propose our own orbital ranking scheme that estimates the importance of an orbital for the active space through a pair-interaction framework from orbital energies and features of the Hartree–Fock exchange matrix. Here, we call this new scheme the “approximate pair coefficient” (APC) method and we show that it performs quite well for the test systems presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spin–Orbit Natural Transition Orbitals and Spin-Forbidden Transitions

Natural transition orbitals (NTOs) are in widespread use for visualizing and analyzing electronic transitions. Here, the present work introduces the analysis of formally spin-forbidden transitions with the help of complex-valued spin-orbit (SO) NTOs. The analysis specifically focuses on the components in such transitions that cause their intensity to be non-zero because of SO coupling. Transition properties such as transition dipole moments are partitioned into SO-NTO hole-particle pairs, such that contributions to the intensity from specific occupied and unoccupied orbitals are obtained. The method has been implemented within the restricted active space (RAS) self-consistent field wave function theory framework, with SO coupling treated by RAS state interaction. SO-NTOs have a broad range of potential applications, which is illustrated by the T 2 -S 1 state mixing in pyrazine, spin-forbidden vs. spin-allowed 4f–5d transitions in the Tb 3+ ion, and the phosphorescence of tris(2-phenylpyridine) iridium [Ir(ppy) 3 ].

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An Orbital Basis Set for Double Photoionization of Atoms and Molecules

The ab initio theoretical treatment of one-photon double photoionization processes has been limited to atoms and diatomic molecules by the challenges posed by large grid-based representations of the double ionized continuum wave function. To provide a path for extensions to polyatomics, an energy-adapted orbital basis approach is demonstrated that reduces the dimensions of such representations and simultaneously allows larger time steps in time-dependent computational descriptions of double ionization. Additionally, an algorithm that exploits the diagonal nature of the two-electron integrals in the grid basis and dramatically accelerates the transformation between grid and orbital representations is presented. Excellent agreement between the present results and benchmark theoretical calculations is found for H – and Be atoms, as well as the hydrogen molecule, including for the triply differential cross sections that relate the angular distribution and energy sharing of all of the particles in the molecular frame.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗