Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variational wave functional methods”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

The scattering of low energy positrons by helium

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.

Humberston, J. W.↗
NASA NTRS1993

The quantum dynamics of electronically nonadiabatic chemical reactions

Considerable progress was achieved on the quantum mechanical treatment of electronically nonadiabatic collisions involving energy transfer and chemical reaction in the collision of an electronically excited atom with a molecule. In the first step, a new diabatic representation for the coupled potential energy surfaces was created. A two-state diabatic representation was developed which was designed to realistically reproduce the two lowest adiabatic states of the valence bond model and also to have the following three desirable features: (1) it is more economical to evaluate; (2) it is more portable; and (3) all spline fits are replaced by analytic functions. The new representation consists of a set of two coupled diabatic potential energy surfaces plus a coupling surface. It is suitable for dynamics calculations on both the electronic quenching and reaction processes in collisions of Na(3p2p) with H2. The new two-state representation was obtained by a three-step process from a modified eight-state diatomics-in-molecules (DIM) representation of Blais. The second step required the development of new dynamical methods. A formalism was developed for treating reactions with very general basis functions including electronically excited states. Our formalism is based on the generalized Newton, scattered wave, and outgoing wave variational principles that were used previously for reactive collisions on a single potential energy surface, and it incorporates three new features: (1) the basis functions include electronic degrees of freedom, as required to treat reactions involving electronic excitation and two or more coupled potential energy surfaces; (2) the primitive electronic basis is assumed to be diabatic, and it is not assumed that it diagonalizes the electronic Hamiltonian even asymptotically; and (3) contracted basis functions for vibrational-rotational-orbital degrees of freedom are included in a very general way, similar to previous prescriptions for locally adiabatic functions in various quantum scattering algorithms.

Truhlar, Donald G.↗
DOE OSTI2020

Starting-point-independent quantum Monte Carlo calculations of iron oxide

Quantum Monte Carlo (QMC) methods are useful for studies of strongly correlated materials because they are many body in nature and use the physical Hamiltonian. Typical calculations assume as a starting point a wave function constructed from single-particle orbitals obtained from one-body methods, e.g., density functional theory. However, mean-field-derived wave functions can sometimes lead to systematic QMC biases if the mean-field result poorly describes the true ground state. In this study, we examine the accuracy and flexibility of QMC trial wave functions using variational and fixed-node diffusion QMC estimates of the total spin density and lattice distortion of antiferromagnetic iron oxide (FeO) in the ground state B1 crystal structure. We found that for relatively simple wave functions the predicted lattice distortion was controlled by the choice of single-particle orbitals used to construct the wave function, rather than by subsequent wave function optimization techniques within QMC. By optimizing the orbitals with QMC, we then demonstrate starting-point independence of the trial wave function with respect to the method by which the orbitals were constructed by demonstrating convergence of the energy, spin density, and predicted lattice distortion for two qualitatively different sets of orbitals. The results suggest that orbital optimization is a promising method for accurate many-body calculations of strongly correlated condensed phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

VQE method: a short survey and recent developments

Abstract The variational quantum eigensolver (VQE) is a method that uses a hybrid quantum-classical computational approach to find eigenvalues of a Hamiltonian. VQE has been proposed as an alternative to fully quantum algorithms such as quantum phase estimation (QPE) because fully quantum algorithms require quantum hardware that will not be accessible in the near future. VQE has been successfully applied to solve the electronic Schrödinger equation for a variety of small molecules. However, the scalability of this method is limited by two factors: the complexity of the quantum circuits and the complexity of the classical optimization problem. Both of these factors are affected by the choice of the variational ansatz used to represent the trial wave function. Hence, the construction of an efficient ansatz is an active area of research. Put another way, modern quantum computers are not capable of executing deep quantum circuits produced by using currently available ansatzes for problems that map onto more than several qubits. In this review, we present recent developments in the field of designing efficient ansatzes that fall into two categories—chemistry–inspired and hardware–efficient—that produce quantum circuits that are easier to run on modern hardware. We discuss the shortfalls of ansatzes originally formulated for VQE simulations, how they are addressed in more sophisticated methods, and the potential ways for further improvements.

Fedorov, Dmitry A. (ORCID:0000000316598580)↗
DOE OSTI2025

Overset-Grid Method with Smooth Orbital Partitioning for Molecular Scattering Calculations

To solve molecular photoionization and electron scattering problems, we use an overset-grid representation of electronic continuum functions, which has an extended central spherical grid that overlaps small spherical grids (subgrids) centered on each atom of a polyatomic molecule. Here, in this work, we present an improved algorithm that smoothly partitions the total wave function between the central grid and the atomic subgrids. The smooth partitioning allows one to use approximately one-fourth the number of partial waves on the central grid compared to our previous implementation with switching functions. The resulting numerical method for treating electron scattering and photoionization of polyatomic molecules combines the accuracy and flexibility of pure numerical grid representations with the rapid convergence of hybrid combinations of atom-centered basis-set expansions and grid methods. The overset-grid representation is implemented using the complex Kohn variational principle for scattering and photoionization amplitudes. The faster convergence with respect to the number of central grid partial waves is demonstrated and accuracy is verified by comparisons with the previous implementation and with far more computationally demanding single-center numerical expansions in electron-molecule scattering and photoionization calculations on the neon dimer (Ne 2 ) system, carbon tetrafluoride (CF 4 ) molecule, and the pyridine (C 5 H 5 N) molecule in the static-exchange approximation.

Molecules↗
DOE OSTI2020

Atmospheric condition identification in multivariate data through a metric for total variation

Identification of atmospheric conditions within a multivariable atmospheric data set is a necessary step in the validation of emerging and existing high-fidelity models used to simulate wind plant flows and operation.Atmospheric conditions relevant for wind energy research include stationary conditions, given the need for well-converged statistics for model validation, as well as conditions observed less frequently, such as extreme atmospheric events, which are used in wind turbine and wind plant design.Aggregation of observations without regard to covariance between time series discounts the dynamical nature of the atmosphere and is not sufficiently representative of atmospheric conditions.Identification and characterization of continuous time periods with atmospheric conditions that have a high value for analysis or simulation set the stage for more advanced model validation and the development of real-time control and operational strategies.The current work explores a single metric for variation in a multivariate data sample that quantifies variability within each channel as well as covariance between channels.The total variation is used to identify conditions of interest that conform to desired objective functions, such as stationary conditions, ramps or waves of wind speed, and changes in wind direction.Total variation is somewhat sensitive to the presence of outliers in the input data, and the method is best complemented by quality-control procedures to ensure reliable results.The direct detection and classification of events or conditions of interest within atmospheric data sets is vital to developing our understanding of wind plant response and to the formulation of forecasting and control models.

17 WIND ENERGY↗
NASA NTRS1990

Polyatomic molecular Dirac-Hartree-Fock calculations with Gaussian basis sets

Numerical methods have been used successfully in atomic Dirac-Hartree-Fock (DHF) calculations for many years. Some DHF calculations using numerical methods have been done on diatomic molecules, but while these serve a useful purpose for calibration, the computational effort in extending this approach to polyatomic molecules is prohibitive. An alternative more in line with traditional quantum chemistry is to use an analytical basis set expansion of the wave function. This approach fell into disrepute in the early 1980's due to problems with variational collapse and intruder states, but has recently been put on firm theoretical foundations. In particular, the problems of variational collapse are well understood, and prescriptions for avoiding the most serious failures have been developed. Consequently, it is now possible to develop reliable molecular programs using basis set methods. This paper describes such a program and reports results of test calculations to demonstrate the convergence and stability of the method.

Dyall, Kenneth G.↗
NASA NTRS2004

Electron-Hydrogen Elastic Scattering

Scattering by single-electron systems is always of interest because the wave function of the target is known exactly. Various approximations have been employed to take into account distortion produced in the target. Among them are the method of polarized orbitals and the close coupling approximation. Recently, e-H and e-He+ S-wave scattering in the elastic region has been studied using the Feshbach projection operator formalism. In this approach, the usual Hartree-Fock and exchange potentials are augmented by an optical potential and the resulting phase shifts have rigorous lower bounds. Now this method is being applied to the e-H P-wave scattering in the elastic region. The number of terms in the Hylleraas-type wave function for the 1,3 P phase shifts is 84 and the resulting phase shifts (preliminary) are given. The results have been given up to five digits because to that accuracy they are rigorous lower bounds. They are in general agreement with the variational (VAR) results of Armstead, and those obtained from the intermediate energy R-matrix method (RM) of Scholz et al., and the finite element method (FEM) of Botero and Shertzer. The later two methods do not provide any bounds on phase shifts.

Bhatia, A. K.↗
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↗
NASA NTRS1989

Composite structure dynamics and space applicability

Guided wave modes in cladded or uncladded fiber-reinforced composite plates and tubes have been analyzed using a stiffness method in which the displacement variation through the thickness is approximated by polynomial interpolation functions. This allows for an arbitrary number of laminations and fiber orientations different from lamina to lamina. It is shown that dispersive behavior of guided modes depends significantly on the cladding, number of laminae, and interfaces between the adjacent laminae. A hybrid modeling technique is described in which an inner region containing cracks (or other defects) is discretized by finite elements, and the field in the exterior region is represented in terms of modes that are found using the stiffness method described above. It is found that the reflected and transmitted amplitudes of modes vary significantly with the size of a transverse or longitudinal (delamination) crack and frequency. We have also studied the impact response of a unidirectional fiber-reinforced plate. Received signals at the epicentral and other locations are shown. Strong longitudinal anisotropy of the graphite/epoxy plate causes the signal to be considerably different from that in an isotropic plate.

Datta, Subhendu K.↗
NASA NTRS2003

Mass Polarization Effect in He-like Ions: First and Second Order

In a paper with a similar title Yamanaka has calculated the mass polarization effect (to first order in mu/M for several low-lying states of the two-electron atoms and ions with atomic number Z from 2 to 10. Here we improve the previous results by using Hylleraas variational wave functions with up to 560 terms and extend the calculation to include some additional states and the Z=1ground state. In addition, we compute the second-order effect using the method of pseudostate summation. In one appendix another method of computation is discussed and used as a check, while the energies of the Z=1 ions are presented and discussed in a second appendix.

Bhatia, A. K.↗
NASA NTRS1984

A variational method for the calculation of dynamic polarizabilities and two-photon transition moments - The dressed molecule approach

The solution of the time-dependent Schroedinger equation of the molecule + radiation field system is analyzed. A quantized radiation field is used. The relationship between the oscillatory wave function and the dressed state wave function is established, the oscillatory wave function being the solution if the radiation field is turned on at t = 0 and the dressed state being the stationary solution satisfying the boundary condition only if the field is present at t = negative infinity. In general, the oscillatory wave function can be expressed using a complete set of dressed states. However, in the presence of a nonresonant radiation field, the system is well represented by a single dressed state. Molecular properties such as dynamic polarizabilities and two-photon transition moments can be deduced from the dressed state wave function instead of the oscillatory wave function as in previous methods. Because of its stationary character, the dressed state is more amenable to approximate calculations. A CI method is developed for this purpose. The CI matrix is simple in form and can be readily constructed using existing computer codes. The present method can also be adapted to calculate other optical properties.

Huo, W. M.↗
NASA NTRS1998

Phase Inversion: Inferring Solar Subphotospheric Flow and Other Asphericity from the Distortion of Acoustic Waves

A method is proposed for analyzing an almost monochromatic train of waves propagating in a single direction in an inhomogeneous medium that is not otherwise changing in time. An effective phase is defined in terms of the Hilbert transform of the wave function, which is related, via the JWKB approximation, to the spatial variation of the background state against which the wave is propagating. The contaminating effect of interference between the truly monochromatic components of the train is eliminated using its propagation properties. Measurement errors, provided they are uncorrelated, are manifest as rapidly varying noise; although that noise can dominate the raw phase-processed signal, it can largely be removed by low-pass filtering. The intended purpose of the analysis is to determine the distortion of solar oscillations induced by horizontal structural variation and material flow. It should be possible to apply the method directly to sectoral modes. The horizontal phase distortion provides a measure of longitudinally averaged properties of the Sun in the vicinity of the equator, averaged also in radius down to the depth to which the modes penetrate. By combining such averages from different modes, the two-dimensional variation can be inferred by standard inversion techniques. After taking due account of horizontal refraction, it should be possible to apply the technique also to locally sectoral modes that propagate obliquely to the equator and thereby build a network of lateral averages at each radius, from which the full three-dimensional structure of the Sun can, in principle, be determined as an inverse Radon transform.

Gough, Douglas↗
NASA NTRS1958

Application of the Method of Coordinate Perturbation to Unsteady Duct Flow

The method of coordinate perturbation is applied to the unsteady flow of a compressible fluid in ducts of variable cross section. Solutions, in the form of perturbation series, are obtained for unsteady flows in ducts for which the logarithmic derivative of area variation with respect to the space coordinate is a function of the 'smallness' parameter of the perturbation series. This technique is applied to the problem of the interaction of a disturbance and a shock wave in a diffuser flow. It is found that, for a special choice of the function describing the disturbance, the path of the shock wave can be expressed in closed form to first order. The method is then applied to the determination of the flow field behind a shock wave moving on a prescribed path in the x,t-plane. Perturbation series solutions for quite general paths are developed. The perturbation series solutions are compared with the more exact solutions obtained by the application of the method of characteristics. The approximate solutions are shown to be in reasonably accurate agreement with the solutions obtained by the method of characteristics.

Himmel, Seymour C.↗
DOE OSTI2020

A Generalized Variational Principle with Applications to Excited State Mean Field Theory

We present a generalization of the variational principle that is compatible with any Hamiltonian eigenstate that can be specified uniquely by a list of properties. This variational principle appears to be compatible with a wide range of electronic structure methods, including mean field theory, density functional theory, multireference theory, and quantum Monte Carlo. Like the standard variational principle, this generalized variational principle amounts to the optimization of a nonlinear function that, in the limit of an arbitrarily flexible wave function, has the desired Hamiltonian eigenstate as its global minimum. Unlike the standard variational principle, it can target excited states and select individual states in cases of degeneracy or near-degeneracy. As an initial demonstration of how this approach can be useful in practice, we employ it to improve the optimization efficiency of excited state mean field theory by an order of magnitude. With this improved optimization, we are able to demonstrate that the accuracy of the corresponding second-order perturbation theory rivals that of singles-and-doubles equation-of-motion coupled cluster in a substantially broader set of molecules than could be explored by our previous optimization methodology.

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

Nondestructive Evaluation of Concrete: Elastic Property Imaging Through Full-Waveform Inversion

Concrete is a vital material in construction—especially in the nuclear industry, where it is used in critical structures such as containment vessels. Over time, concrete can degrade due to harsh operational and environmental conditions, necessitating that its elastic properties be accurately evaluated to ensure structural integrity and safety. Traditional nondestructive evaluation methods such as ultrasound-based techniques often rely on simplifying assumptions that may not hold true for concrete. This paper presents an advanced ultrasound-based method that uses elastic full-waveform inversion (EFWI) to create detailed images of concrete’s mechanical properties. By accurately modeling wave behaviors such as scattering and reflection, we aim to overcome the limitations of conventional ultrasonic-based methods. In this work, the imaging problem involved reconstructing the various elastic properties of a heterogenous concrete block with three steel rebars embedded in it. The ultrasonic measurements were synthetically generated from multiple sources and receivers, and the reconstruction process was performed using a gradient-based optimization algorithm. Our approach leveraged EFWI to reconstruct high-resolution images of the pressure wave speed, shear wave speed, and density. Multiple misfit functions—including L2-norm, cross-correlation (CC), and L1-norm—combined with total variation (TV) regularization and parameter constraints using a Sigmoid function—were explored for the reconstruction. The results demonstrated that using the L1-norm misfit function in conjunction with TV regularization and Sigmoid constraints significantly improved the reconstruction quality in comparison to traditional methods. This approach provided clearer images with fewer artifacts and better captured background heterogeneity. Our findings highlight that, when properly designed, EFWI carries great potential for providing comprehensive, more accurate, and more reliable assessments of concrete conditions, as is crucial for the maintenance and safety of nuclear power plant structures.

97 - MATHEMATICS AND COMPUTING↗