Search NASA⌕ Search

SEARCH · Search NASA

Results for “variational method”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Fast & accurate emulation of two-body scattering observables without wave functions

We combine Newton's variational method with ideas from eigenvector continuation to construct a fast & accurate emulator for two-body scattering observables. The emulator will facilitate the application of rigorous statistical methods for interactions that depend smoothly on a set of free parameters. Our approach begins with a trial K or T matrix constructed from a small number of exact solutions to the Lippmann–Schwinger equation. Subsequent emulation only requires operations on small matrices. We provide several applications to short-range potentials with and without the Coulomb interaction and partial-wave coupling. It is shown that the emulator can accurately extrapolate far from the support of the training data. When used to emulate the neutron-proton cross section with a modern chiral interaction as a function of 26 free parameters, it reproduces the exact calculation with negligible error and provides an over 300x improvement in CPU time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Understanding the origin of Tabula Rasa-induced defects in n-type Cz c-Si: The case of nitrogen atmosphere

Phosphorus-doped Czochralski-grown silicon (Cz-Si) has been gaining market share in the large-scale manufacturing of high-efficiency silicon (Si)-based photovoltaic (PV) devices thanks to higher carrier lifetimes than their boron-doped counterpart. However, the fabrication of n-type Cz-Si based solar cells often requires process steps with much higher temperatures and longer times than p-type Silicon. Defect interaction with the high temperatures during such processes tend to be detrimental to the n-type Cz-Si carrier lifetime, therefore limiting the final device efficiency. Short thermal anneals before cell processing, known as Tabula Rasa (TR), have been proposed to mitigate the thermally induced lifetime degradation during n-type Cz-Si solar cell fabrication. This work thoroughly investigates the defects responsible for the lifetime degradation after TR in a N 2 atmosphere treatment. We use temperature-injection-dependent lifetime spectroscopy and the thickness variation method to decouple the effects of TR treatment in the bulk and the surface of the n-type Cz-Si wafers. Using the defect parameter contour mapping (DPCM), we identify the defect energy level (E t ) and the capture cross-section ratio (k) of the most likely process-induced defect, which aligns with previously proposed Si vacancy-associated defects. The DPCM reveals that these vacancy-associated defects have a shallow energy level E t - E v ~0.13 eV and very efficient electron capture cross section k~600. Unexpectedly, the bulk degradation due to vacancy defects in the volume of the wafer, is accompanied by a significant increase in the surface recombination as well. Through evaluating the surface recombination velocity temperature- and injection dependence, we show that after TR, at room temperature and for an injection level of 10 15 cm -3 , in a wafer passivated with a-Si:H(i) the surface recombination dominates the overall lifetime response. Here we hypothesize that the near surface vacancy-associated bulk defects play a role in lowering the electron diffusion current into the a-Si:H(i) from the c-Si(n) reducing the field-effect passivation.

14 SOLAR ENERGY↗
DOE OSTI2023

Revealing Variable Dependences in Hexagonal Boron Nitride Synthesis via Machine Learning

Wafer-scale monolayer two-dimensional (2D) materials have been realized by epitaxial chemical vapor deposition (CVD) in recent years. To scale up the synthesis of 2D materials, a systematic analysis of how the growth dynamics depend on the growth parameters is essential to unravel its mechanisms. However, the studies of CVD-grown 2D materials mostly adopted the control variate method and considered each parameter as an independent variable, which is not comprehensive for 2D materials growth optimization. Herein, we synthesized a representative 2D material, monolayer hexagonal boron nitride (hBN), on single-crystalline Cu (111) by epitaxial chemical vapor deposition and varied the growth parameters to regulate the hBN domain sizes. Furthermore, we explored the correlation between two growth parameters and provided the growth windows for large flake sizes by the Gaussian process. Here, this new analysis approach based on machine learning provides a more comprehensive understanding of the growth mechanism for 2D materials.

36 MATERIALS SCIENCE↗
DOE OSTI2019

The Microphysics of Stratiform Precipitation During OLYMPEX: Compatibility Between Triple‐Frequency Radar and Airborne In Situ Observations

Abstract The link between stratiform precipitation microphysics and multifrequency radar observables is thoroughly investigated by exploiting simultaneous airborne radar and in situ observations collected from two aircraft during the OLYMPEX/RADEX (Olympic Mountain Experiment/Radar Definition Experiment 2015) field campaign. Above the melting level, in situ images and triple‐frequency radar signatures both indicate the presence of moderately rimed aggregates. Various mass‐size relationships of ice particles and snow scattering databases are used to compute the radar reflectivity from the in situ particle size distribution. At Ku and Ka band, the best agreement with radar observations is found when using the self‐similar Rayleigh‐Gans approximation for moderately rimed aggregates. At W band, a direct comparison is challenging because of the non‐Rayleigh effects and of the probable attenuation due to ice aggregates and supercooled liquid water between the two aircraft. A variational method enables the retrieval of the full precipitation profile above and below the melting layer, by combining the observations from the three radars. Even with three radar frequencies, the retrieval of rain properties is challenging over land, where the integrated attenuation is not available. Otherwise, retrieved mean volume diameters and water contents of both solid and liquid precipitation are in agreement with in situ observations and indicate local changes of the degree of riming of ice aggregates, on the scale of 5 km. Finally, retrieval results are analyzed to explore the validity of using continuity constraints on the water mass flux and diameter within the melting layer in order to improve retrievals of ice properties.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

The phase diagram of quantum chromodynamics in one dimension on a quantum computer

The quantum chromodynamics (QCD) phase diagram, which reveals the state of strongly interacting matter at different temperatures and densities, is key to answering open questions in physics, ranging from the behaviour of particles in neutron stars to the conditions of the early universe. However, classical simulations of QCD face significant computational barriers, such as the sign problem at finite matter densities. Quantum computing offers a promising solution to overcome these challenges. Here, we take an important step toward exploring the QCD phase diagram with quantum devices by preparing thermal states in one-dimensional non-Abelian gauge theories. We experimentally simulate the thermal states of SU(2) and SU(3) gauge theories at finite densities on a trapped-ion quantum computer using a variational method. This is achieved by introducing two features: Firstly, we add motional ancillae to the existing qubit register to efficiently prepare thermal probability distributions. Secondly, we introduce charge-singlet measurements to enforce colour-neutrality constraints. This work pioneers the quantum simulation of QCD at finite density and temperature for two and three colours, laying the foundation to explore QCD phenomena on quantum platforms.

Quantum information↗
DOE OSTI2022

Preparation of metrological states in dipolar-interacting spin systems

Abstract Spin systems are an attractive candidate for quantum-enhanced metrology. Here we develop a variational method to generate metrological states in small dipolar-interacting spin ensembles with limited qubit control. For both regular and disordered spatial spin configurations the generated states enable sensing beyond the standard quantum limit (SQL) and, for small spin numbers, approach the Heisenberg limit (HL). Depending on the circuit depth and the level of readout noise, the resulting states resemble Greenberger-Horne-Zeilinger (GHZ) states or Spin Squeezed States (SSS). Sensing beyond the SQL holds in the presence of finite spin polarization and a non-Markovian noise environment. The developed black-box optimization techniques for small spin numbers ( N ≤ 10) are directly applicable to diamond-based nanoscale field sensing, where the sensor size limits N and conventional squeezing approaches fail.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Nonthermal effects on the electron-atom bremsstrahlung process in generalized Lorentzian plasmas

The nonthermal effects on the electron-atom bremsstrahlung process are investigated in generalized Lorentzian plasmas. The screened atomic wave function and energy for the hydrogen atom in Lorentzian plasmas are obtained by the Ritz variational method. The electron-atom bremsstrahlung radiation cross section in Lorentzian plasmas is also obtained as a function of the spectral index and plasma parameters. It is shown that the nonthermal character of the Lorentzian plasma strongly suppresses the electron-atom bremsstrahlung cross section. In addition, it is shown that the nonthermal effect on the bremsstrahlung process decreases with increasing radiation photon energy. It is also shown that the nonthermal effect on the bremsstrahlung cross section is more significant in plasmas with higher densities.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2012

Shape coexistence, shape evolution and Gamow-Teller {beta}-decay of neutron-rich A Asymptotically-Equal-To 100 nuclei

The structure of neutron-rich nuclei in the A Asymptotically-Equal-To 100 mass region relevant for the astrophysical r process manifests drastic changes in some isotopic chains and often sudden variations of particular nuclear properties have been identified. For a realistic description of the evolution in structure with increasing energy, spin, and isospin determined by shape coexistence and mixing beyond-mean-field approaches are required. Our recent studies represent an attempt to the self-consistent description of the shape coexistence phenomena in neutron-rich A Asymptotically-Equal-To 100 nuclei within the complex Excited Vampir variational model with symmetry projection before variation using a realistic effective interaction based on the Bonn A potential in a large model space. Results concerning the triple shape coexistence and the shape evolution in the N=58 Sr and Zr isotopes, the shape evolution in a chain of Zr nuclei, as well as the Gamow-Teller {beta}-decay properties of neutron-rich Zr and Tc nuclei are presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

Binding energy levels of a slowly moving ion in dusty plasmas

The near field electric potential of a slowly moving ion in complex plasmas is studied. We find that the potential consists of the Debye-Hueckel potential, the wake potential, and the potential associated with charge fluctuations. The binding energy levels of the ion are calculated by use of the Ritz variation method. The results show that the binding energy levels are related to the magnetic quantum number m. The binding energy levels are affected by speed of the ion and dust grain number density. In contract to isolated ion or static ion in plasmas, the binding energy levels of the ion are pushed up and even become unbounded.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Hamiltonian magnetohydrodynamics: Lagrangian, Eulerian, and dynamically accessible stability—Theory

Stability conditions of magnetized plasma flows are obtained by exploiting the Hamiltonian structure of the magnetohydrodynamics (MHD) equations and, in particular, by using three kinds of energy principles. First, the Lagrangian variable energy principle is described and sufficient stability conditions are presented. Next, plasma flows are described in terms of Eulerian variables and the noncanonical Hamiltonian formulation of MHD is exploited. For symmetric equilibria, the energy-Casimir principle is expanded to second order and sufficient conditions for stability to symmetric perturbation are obtained. Then, dynamically accessible variations, i.e., variations that explicitly preserve invariants of the system, are introduced and the respective energy principle is considered. General criteria for stability are obtained, along with comparisons between the three different approaches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Erratum: “A finite difference Davidson procedure to sidestep full ab initio hessian calculation: Application to characterization of stationary points and transition state searches”

Certain key contributions to the field were not presented in the Introduction section of our recent paper. In addition to the variational method described in Ref. 14 and iterative diagonalization methods in Refs. 15–17, we wish to direct the attention of the reader to other significant efforts that have utilized similar principles for various applications. Reiher and Neugebauer developed a mode-tracking algorithm using the Davidson method for calculating only the modes relevant for vibrational analysis. Similarly, a block-Davidson approach via finite differentiation of the gradient was implemented for calculating hessian-vector products for normal mode analysis. Deglmann and Furche developed a method for stationary point characterization with density functional theory where full solution of all coupled-perturbed Kohn-Sham equations, the bottleneck in hessian calculation, could be avoided through iterative calculation of hessian-vector products to determine only the lowest eigenvalues.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2015

Thermodynamical properties of liquid lanthanides-A variational approach

Thermodynamical properties like Entropy (S), Internal energy (E) and Helmholtz free energy (F) of liquid lanthanides using a variation principle based on the Gibbs-Bogoliubuv (GB) inequality with Percus Yevick hard sphere reference system have been reported in the present investigation. To describe electron-ion interaction we have used our newly constructed parameter free model potential along with Sarkar et al. local field correction function. Lastly, we conclude that our newly constructed model potential is capable to explain the thermodynamical properties of liquid lanthanides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Effects of the built-in electric field on polaron effects in a wurtzite ZnO/Mg{sub x}Zn{sub 1−x}O quantum well

The effects of the built-in electric field on the polarons in wurtzite ZnO/Mg{sub x}Zn{sub 1−x}O quantum wells are studied numerically via the improved Lee-Low-Pines intermediate coupling variational method. The contributions of the different branches of the optical phonons to the polaron energies of the ZnO/Mg{sub x}Zn{sub 1−x}O quantum wells are calculated as functions of well width d and composition x. The anisotropy effects of the electronic effective masses, the dielectric constants, and the frequencies of the different branches of the phonon modes (including both the longitudinal-like and transverse-like confined optical phonon modes, the interface optical phonon modes and the half-space phonon modes) on the polaron energies are considered in the calculations. Comparisons between the cases with and without the built-in electric field (F ≠ 0 and F = 0) are made for the optical phonon contributions to the polaron energies for the different branches of the phonon modes. The results show that the built-in electric field has marked effects on the contributions of the phonons with the different modes; in detail, it makes positive contributions to the interface and the half-space phonons, but negative contributions to the confined phonons, and thus its effect on the total phonon contribution is not obvious. Detailed comparisons of the contributions of the symmetric and antisymmetric phonon modes to the polaron energies as functions of d and x are also presented and a heuristic argument is provided to explain the numerical results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Theory of torsional control for G 16 -type molecules

Here, we introduce a four-dimensional quantum model for describing the torsional control with moderately strong, non-resonant laser pulses of G 16 -type molecules in the electronic ground state, based on the symmetry-adapted variational method. We define conditions for which lower-dimensional models, commonly used to simulate the strong-field control of molecular torsions, are reliable approximations to a four-dimensional treatment. In particular, we study the role of different types of rotational-torsional couplings—the field-free coupling and the field-induced coupling—and show that the conclusions recently drawn on the role of rotational-torsional couplings in the process of torsional alignment are not correct in general. Furthermore, we demonstrate how important an adequate description of the molecular polarizability is for reliably predicting the torsional alignment.

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

ADM Analysis of gravity models within the framework of bimetric variational formalism

Bimetric variational formalism was recently employed to construct novel bimetric gravity models. In these models an affine connection is generated by an additional tensor field which is independent of the physical metric. In this work we demonstrate how the ADM decomposition can be applied to study such models and provide some technical intermediate details. Using ADM decomposition we are able to prove that a linear model is unstable as has previously been indicated by perturbative analysis. Moreover, we show that it is also very difficult if not impossible to construct a non-linear model which is ghost-free within the framework of bimetric variational formalism. However, we demonstrate that viable models are possible along similar lines of thought. To this end, we consider a set up in which the affine connection is a variation of the Levi-Civita one. As a proof of principle we construct a gravity model with a massless scalar field obtained this way.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Accessing the gluon momentum fraction of nucleons through the gradient flow

We calculate the gluon momentum fraction of the nucleon using lattice QCD, with a nonperturbative renormalization technique based on the gradient flow. The gluon momentum fraction is determined on a single Wilson-clover ensemble using 𝑁 𝑓 =2 +1 flavors with pion mass 358 MeV and lattice spacing 0.094 fm. We employ the variational method to reduce excited-state contamination and apply the distillation framework to ensure a large operator basis. To reduce systematic uncertainties, we apply Bayesian model averaging to all fit procedures. We apply matching coefficients to the flow-time dependent lattice results to recover the gluon momentum fraction in the $\overline{MS}$-scheme at 2 GeV. Our final result is ⟨𝑥⟩ 𝑔 ⁢(𝜇 =2 GeV) =0.482⁢(35), where we quote only statistical uncertainties.

Lattice QCD↗
DOE OSTI2020

Quantum state reduction: Generalized bipartitions from algebras of observables

Reduced density matrices are a powerful tool in the analysis of entanglement structure, approximate or coarse-grained dynamics, decoherence, and the emergence of classicality. It is straightforward to produce a reduced density matrix with the partial-trace map by tracing out part of the quantum state, but in many natural situations this reduction may not be achievable. We investigate the general problem of identifying how the quantum state is reduced given a restriction on the observables. For example, in an experimental setting, the set of observables that can actually be measured is usually modest (compared to the set of all possible observables) and their resolution is limited. In such situations, the appropriate state-reduction map can be defined via a generalized bipartition, which is associated with the structure of irreducible representations of the algebra generated by the restricted set of observables. One of our main technical results is a general, not inherently numeric, algorithm for finding irreducible representations of matrix algebras. In our work, we demonstrate the viability of this approach with two examples of limited-resolution observables. The definition of quantum state reductions can also be extended beyond algebras of observables. To accomplish this task we introduce a more flexible notion of bipartition, the partial bipartition, which describes coarse grainings preserving information about a limited set (not necessarily algebra) of observables. We describe a variational method to choose the coarse grainings most compatible with a specified Hamiltonian, which exhibit emergent classicality in the reduced state space. We apply this construction to the concrete example of the one-dimensional Ising model. Our results have relevance for quantum information, bulk reconstruction in holography, and quantum gravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Validity of the static-exchange approximation for inner-shell photoionization of polyatomic molecules

The simple single-channel static-exchange approximation completely ignores correlation between the continuum and molecular ion electrons. In molecular systems with symmetry equivalent atoms, the single-channel approximation can seriously fail in core ionization when using delocalized orbitals to represent the core hole states. We present cross sections and molecular frame photoelectron angular distributions with both localized and delocalized core orbitals in CF 4 F (1 s ) ionization. We show that only a full coupled-channel calculation can recover an accurate description of the physics of inner-shell photoionization when using delocalized orbitals, whereas nearly the same result can be obtained from independent single-channel static-exchange calculations when localized core orbitals are used. A grid-based variational method described here makes such single-channel calculations possible on larger systems without local-exchange approximations. Illustrative calculations on the core ionization of SF 6 are presented to illustrate the power of the grid-based method.

74 ATOMIC AND MOLECULAR PHYSICS↗