SEARCH · Search NASA
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.
Calculation of e - H total and elastic partial wave cross sections is being carried out using the complex correlation variational T-matrix method. In this preliminary study, elastic partial wave phase shifts are calculated with the correlation functions which are confined to be real. In that case the method reduces to the conventional optical potential approach with projection operators.
Calculation of e-H total and elastic partial wave cross sections is being carried out using the complex correlation variational T-matrix method. In this preliminary study, elastic partial wave phase shifts are calculated with the correlation functions which are confined to be real. In that case the method reduces to the conventional optical potential approach with projection operators. The number of terms in the Hylleraas-type wave function for the S phase shifts is 95 while for the S it is 56, except for k=0.8 where it is 84. Our results, which are rigorous lower bounds, are given. They are seen to be in general agreement with those of Schwartz, but they are of 0 greater accuracy and outside of his error limits for k=0.3 and 0.4 for S. The main aim of this approach' is the application to higher energy scattering. By virtue of the complex correlation functions, the T matrix is not unitary so that elastic and total scattering cross sections are independent of each other. Our results will be compared specifically with those of Bray and Stelbovics.
Here, a generalized, projection-based transformation of the method-agnostic Fock operator in various ab initio fragment-based quantum chemistry methods has been developed for the treatment of interfragment covalent bonds. This transformation freezes the relevant localized molecular orbital associated with each interfragment bond, thereby restricting the variational subspace of the fragment wave functions, in order to maintain the proper physical characteristics of the involved covalent bonds. In addition, sets of orbitals that would lead to multiple occupancy of certain orbitals are explicitly removed from the variational space. The transformation is developed for the specific case of mutually orthonormal frozen and unfrozen orbitals within each fragment. The newly developed approach is then used to study model systems with two popular ab initio fragment-based methods, and the results of these calculations are compared to those obtained by existing methodologies. Analysis is focused on both quantitative and qualitative accuracy as well as computational scalability and stability. Other methods for which the developed formalisms are appropriate are outlined, and future extensions of the methods are discussed.
Context. Recent studies of the solar wind sunward of 0.25 AU reveal that it contains quiescent regions, with low-amplitude plasma and magnetic field fluctuations, and a magnetic field direction similar to the Parker spiral. The quiescent regions are thought to have a more direct magnetic connection to the solar corona than other types of solar wind, suggesting that waves or instabilities in the quiescent regions are indicative of the early evolution of the solar wind as it escapes the corona. The quiescent solar wind regions are highly unstable to the formation of plasma waves near the electron cyclotron frequency (fce). Aims. We examine high time resolution observations of these waves in an effort to understand their impact on electron distribution functions of the quiescent near-Sun solar wind. Methods. High time resolution waveform captures of near-fce waves were examined to determine variations of their amplitude and frequency in time as well as their polarization properties. Results. We demonstrate that the near-fce wave intervals contain several distinct wave types, including electron Bernstein waves and extremely narrowband waves that are highly sensitive to the ambient magnetic field orientation. Using the properties of these waves, we suggest possible plasma wave mode classifications and possible instabilities that generate these waves. The results of this analysis indicate that these waves may modify the cold core of the electron distribution functions in the quiescent near-Sun solar wind.
In the variational subspace valence bond (VSVB) [G. D. Fletcher, J. Chem. Phys. 142, 134112 (2015)] method, the electronic orbitals comprising the wave function correspond to chemically meaningful objects, such as bonds, lone pairs, atomic cores, and so on. Selected regions of a molecule (for example, a single chemical bond, as opposed to all the bonds) can be modeled with different levels of basis set and possible methods for modeling correlation from the other regions. The interactions between the components of a molecule (say, a bond and a neighboring orbital) can then be studied in detail for their impact on a chemical phenomenon while avoiding the expense of necessarily applying the higher levels and methods to the entire molecule. Here, this work presents the theoretical basis for modeling correlation effects between specific electron pairs by incorporating terms in the inter-electronic coordinates (“r 12 ”) into VSVB. The approach is validated with calculations on small systems using single-reference wave functions.
A numerically exact spectral method for solving the time-dependent Schroedinger equation in spherical coordinates is described. The angular dependence of the wave function is represented on a two-dimensional grid of evenly spaced points. The fast Fourier transform algorithm is used to transform between the angle space representation of the wave function and its conjugate representation in momentum space. The time propagation of the wave function is evaluated using an expansion of the time evolution operator as a series of Chebyshev polynomials. Calculations performed for a model system representing H2 scattering from a rectangular corrugated surface yield transition probabilities that are in excellent agreement with those obtained using the close-coupling wave packet (CCWP) method. However, the new method is found to require substantially more computation time than the CCWP method because of the large number of grid points needed to represent the angular dependence of the wave function and the variation in the number of terms required in the Chebyshev representation of the time evolution operator.
While the electron scatters from the target, the target no longer stays in its original form. One of the first methods to take into account the distortion of the target at low incident energies is the method of polarized orbitals. In this method the wave function for the process is written using the first-order perturbation theory and the equation for the scattering function is derived from the Schradinger equation. This method has been very successful in calculating the phase shifts and therefore the cross sections at various energies. The total wave function can be used to calculate photoionization cross sections. The disadvantage of this approach is that the method is not variational and therefore does not provide bounds on the phase shifts. These difficulties can be overcome by using the Feshbach projection operator formalism. This approach has been employed for the scattering of electrons and positrons from targets. Results of various calculations will be discussed.
Ab initio quantum-chemical methods that perform well for computing the electronic ground state are not straightforwardly transferable to electronically excited states, particularly in large molecular systems. Wave function theory offers high accuracy, but is often prohibitively expensive. Methods based on time-dependent density functional theory (TD-DFT) are crucially sensitive to the chosen exchange-correlation functional (XCF) parameterization, and system-specific tuning protocols were therefore proposed to address the method's robustness. Methods based on the variational relaxation of the excited-state electron density showcased promising results for the calculation of charge-transfer excitations, but the complex shape of the electronic hypersurface makes convergence to a specific excited state much more difficult than for the ground state when standard variational techniques are applied. We address the latter aspect by providing suitable initial guesses, which we obtain by two separate constrained algorithms. Combined with the squared-gradient minimization algorithm for all-electrons relaxation in a freeze-and-release scheme (FRZ-SGM), we demonstrate that orbital-optimized density functional theory (OO-DFT) calculations can reliably converge to the charge-transfer states of interest even for large molecular systems. We test the FRZ-SGM method on a phenothiazine-anthraquinone CT excitation in a supramolecular Pd(II) coordination cage complex as a function of the cage conformation. This compound has been studied experimentally prior to our work. We compare this freeze-and-release scheme to two XCF reparameterizations, which were recently proposed as low-cost TD-DFT-based alternatives to variational methods. Two dye-semiconductor complexes, which were previously investigated in the context of photovoltaic applications, serve as a second example to investigate the convergence and stability of the FRZ-SGM approach. Our results demonstrate that FRZ-SGM provides reliable convergence for charge-transfer excited states and avoids variational collapse to lower-lying electronic states, whereas time-dependent DFT calculations with an adequate tuning procedure for the range-separation parameter provide a computationally efficient initial estimate of the corresponding energies, with a computational cost comparable to that of configuration-interaction singles (CIS) calculations.
The quantum mechanical problem of reactive or nonreactive scattering of atoms and molecules is formulated in terms of square-integrable basis sets with variational expressions for the reactance matrix. Several formulations involving expansions of the wave function (the Schwinger variational principle) or amplitude density (a generalization of the Newton variational principle), single-channel or multichannel distortion potentials, and primitive or contracted basis functions are presented and tested. The test results, for inelastic and reactive atom-diatom collisions, suggest that the methods may be useful for a variety of collision calculations and may allow the accurate quantal treatment of systems for which other available methods would be prohibitively expensive.
Accurately describing strongly correlated systems with affordable quantum resources remains a central challenge for quantum chemistry applications on near and intermediate term quantum computers. The localized active space self-consistent field (LASSCF) approximates the complete active space self-consistent field (CASSCF) by generating active space-based wave functions within specific fragments while treating interfragment correlation with mean-field approach, hence is computationally less expensive. Hardware-efficient ansatzes (HEA) offer affordable and shallower circuits, yet they often fail to capture the necessary correlation. Previously, Jastrow-factor-inspired nonunitary qubit operators were proposed to use with HEA for variational quantum eigensolver (VQE) calculations (so-called nuVQE), as they do not increase circuit depths and recover correlation beyond the mean-field level for Hartree–Fock initial states. Here, in this study, we explore running nuVQE with LASSCF as the initial state. The method, named LAS-nuVQE, is shown to recover interfragment correlations, reach chemical accuracy with a small number of gates (<70) in both H 4 and square cyclobutadiene (C 4 H 4 ), and produces more accurate energetics than its HEA counterparts at all circuit depths. To further address the inherent symmetry-breaking in HEA, we implemented spin-constrained LAS-nuVQE to extend the capabilities of HEA further and show spin-pure results for square cyclobutadiene. We also mitigate the increased measurement overhead of nuVQE via Pauli grouping and shot-frugal sampling, reducing measurement costs by up to 2 orders of magnitude compared to ungrouped operator, and show that one can achieve better accuracy with a small number of shots (10 3–4 ) per one expectation value calculation compared to noiseless simulations with one or two orders of magnitude more shots. Finally, wall clock time estimates show that, with our measurement mitigation protocols, nuVQE becomes a cheaper and more accurate alternative than vanilla VQE with HEA. Taken together, these developments illustrate a practical pathway toward performing multireference chemical simulations with accuracy and affordable resources on today’s quantum hardware, achieving both accuracy and affordability in challenging correlated systems.
Simulating molecules using the Variational Quantum Eigensolver method is one of the promising applications for NISQ-era quantum computers. Designing an efficient ansatz to represent the electronic wave function is crucial in such simulations. Standard unitary coupled-cluster with singles and doubles (UCCSD) ansatz tends to have a large number of insignificant terms that do not lower the energy of the system. In this work, we present a unitary selective coupled-cluster method, a way to construct a unitary coupled-cluster ansatz iteratively using a selection procedure with excitations up to fourth order. This approach uses the electronic Hamiltonian matrix elements and the amplitudes for excitations already present in the ansatz to find the important excitations of higher order and to add them to the ansatz. The important feature of the method is that it systematically reduces the energy error with increasing ansatz size for a set of test molecules. {The main advantage of the proposed method is that the effort to increase the ansatz does not require any additional measurements on a quantum computer.}
The problem of acoustic wave propagation in a lined duct with non-uniform liner admittance distribution has been analyzed. Two different methods of solution are discussed. Computations are made for a simple example of symmetrically lined two-dimensional duct containing a unit pressure source distribution. The non-uniform admittance variation is assumed to be a homogeneous random function along the direction of wave propagation. The influences of the spatial scale of the admittance variation, the magnitude of the variation and their dependence on mean liner admittance on sound attenuation are presented. The predicted trends agree qualitatively with existing experimental findings.
Here, we compare the accuracy of the ghost rotationally invariant slave-boson (g-RISB) theory and dynamical mean field theory (DMFT) on the single-band Hubbard model, as a function of the number of bath sites in the embedding impurity Hamiltonian. Our benchmark calculations confirm that the accuracy of g-RISB can be systematically improved by increasing the number of bath sites, similar to DMFT. With a few bath sites, we observe that g-RISB is systematically more accurate than DMFT for the ground-state observables. On the other hand, the relative accuracy of these methods is generally comparable for the quasiparticle weight and the spectral function. As expected, we observe that g-RISB satisfies the variational principle in infinite dimensions, as the total energy decreases monotonically towards the exact value as a function of the number of bath sites, suggesting that the g-RISB wave function may approach the exact ground state in infinite dimensions. Our results suggest that the g-RISB is a promising method for first-principles simulations of strongly correlated matter, which can capture the behavior of both static and dynamical observables, at a relatively low computational cost.
Tetraquarks, bound states composed of two quarks and two antiquarks, have been the subject of intense study but are challenging to understand from first principles. We apply variational and Green’s function Monte Carlo methods to compute tetraquark ground-state energies in potential nonrelativistic QCD using a wide range of color and spatial wave functions. We find no evidence for bound tetraquarks composed of equal-mass quarks and antiquarks. Conversely, we find clear evidence for the existence of bound tetraquarks for sufficiently unequal quark/antiquark mass ratios at all overall mass scales where our effective theory results are applicable. We predict the critical mass ratios for bound state formation and study tetraquark bound states’ spatial and color structure at leading order and next-to-leading order in potential nonrelativistic QCD. Published by the American Physical Society 2024
In this work, we present a variational Monte Carlo method that solves the nuclear many-body problem in the occupation number formalism exploiting an artificial neural network representation of the groundstate wave function. A memory-efficient version of the stochastic reconfiguration algorithm is developed to train the network by minimizing the expectation value of the Hamiltonian. We benchmark this approach against widely used nuclear many-body methods by solving a model used to describe pairing in nuclei for different types of interaction and different values of the interaction strength. Despite its polynomial computational cost, our method outperforms coupled-cluster and provides energies that are in excellent agreement with the numerically-exact full configuration interaction values.
Conditions under which optimal wave functions satisfy various time-dependent Hellmann-Feynman theorems
In this study, the lateral variations of the P wave velocity as a function of depth were examined for the regions of the Rio Grande rift and the Jemez lineament, to infer spatial distributions of partial melt in the upper mantle source zones for the rift and the lineament. The method involved measurements of teleismic P wave delays at a 22-station network followed by performing a damped least-squares three-dimensional inversion for these lateral variations. Results indicate that, directly beneath the Jemez lineament (but not beneath the Rio Grande rift), there is a 100-km-wide 1-2-percent low-P-wave-velocity feature in the depth range of 50-160 km. This implies that the volcanic potential of the Jemez lineaments continues to greatly exceed that of the Rio Grande rift.
Kohn's variational method is used to calculate the positron-helium scattering length and low energy S-wave phase shifts for a quite realistic Hylleraas type of helium function containing an electron-electron correlation term. The zero energy wavefunction is used to calculate the value of the annihilation rate parameter Z sub eff. All the results are significantly different from those for Drachman's helium model B, but are in better agreement with the available experimental data.