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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2011

Wave-function functionals for the density

We extend the idea of the constrained-search variational method for the construction of wave-function functionals {psi}[{chi}] of functions {chi}. The search is constrained to those functions {chi} such that {psi}[{chi}] reproduces the density {rho}(r) while simultaneously leading to an upper bound to the energy. The functionals are thereby normalized and automatically satisfy the electron-nucleus coalescence condition. The functionals {psi}[{chi}] are also constructed to satisfy the electron-electron coalescence condition. The method is applied to the ground state of the helium atom to construct functionals {psi}[{chi}] that reproduce the density as given by the Kinoshita correlated wave function. The expectation of single-particle operators W={Sigma}{sub i}r{sub i}{sup n}, n=-2,-1,1,2, W={Sigma}{sub i}{delta}(r{sub i}) are exact, as must be the case. The expectations of the kinetic energy operator W=-(1/2){Sigma}{sub i}{nabla}{sub i}{sup 2}, the two-particle operators W={Sigma}{sub n}u{sup n}, n=-2,-1,1,2, where u=|r{sub i}-r{sub j}|, and the energy are accurate. We note that the construction of such functionals {psi}[{chi}] is an application of the Levy-Lieb constrained-search definition of density functional theory. It is thereby possible to rigorously determine which functional {psi}[{chi}] is closer to the true wave function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Superconductivity in the bilayer Hubbard model: Two Fermi surfaces are better than one

Fully occupied or unoccupied bands in a solid are often considered inert and irrelevant to a material's low-energy properties. But the discovery of enhanced superconductivity in heavily electron-doped FeSe-derived superconductors poses questions about the possible role of incipient bands (those laying close to but not crossing the Fermi level) in pairing. To answer this question, researchers have studied pairing correlations in the bilayer Hubbard model, which has an incipient band for large interlayer hopping t⊥, using many-body perturbation theory and variational methods. They have generally found that superconductivity is enhanced as one of the bands approaches the Lifshitz transition and even when it becomes incipient. Here we address this question using the nonperturbative quantum Monte Carlo (QMC) dynamical cluster approximation (DCA) to study the bilayer Hubbard model's pairing correlations. We find that the model has robust s± pairing correlations in the large t⊥ limit, which can become stronger as one band is made incipient. While this behavior is linked to changes in the effective interaction, we further find that it is counteracted by a suppression of the intrinsic pair-field susceptibility and does not translate to an increased Tc. Our results demonstrate that the highest achievable transition temperatures in the bilayer Hubbard model occur when the system has two bands crossing the Fermi level.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Composite excitonic states in doped semiconductors

We present a theoretical model of composite excitonic states in doped semiconductors. Many-body interactions between a photoexcited electron-hole pair and the electron gas are integrated into a computationally tractable few-body problem, solved by the variational method. We focus on electron-doped monolayer (ML)-MoSe 2 and ML-WSe 2 due to the contrasting character of their conduction bands. In both cases, the core of the composite is a tightly bound trion (two electrons and a valence-band hole) surrounded by a region depleted of electrons. The composite in ML-WSe 2 further includes a satellite electron with different quantum numbers. The theory is general and can be applied to semiconductors with various energy-band properties, allowing one to calculate their excitonic states and to quantify the interaction with the Fermi sea.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Vestigial nematicity from spin and/or charge order in the cuprates

Nematic order has manifested itself in a variety of materials in the cuprate family. We propose an effective field theory of a layered system with incommensurate, intertwined spin- and charge-density wave (SDW and CDW) orders, each of which consists of two components related by C4 rotations. Using a variational method (which is exact in a large N limit), we study the development of nematicity from partially melting those density waves by either increasing temperature or adding quenched disorder. As temperature decreases we first find a transition to a nematic phase, but depending on the range of parameters (e.g. doping concentration) the strongest fluctuations associated with this phase reflect either proximate SDW or CDW order. We also discuss the changes in parameters that can account for the differences in the SDW-CDW interplay between the (214) family and the other hole-doped cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Quantum Computing for Neutrino-Nucleus Scattering

Neutrino-nucleus cross section uncertainties are expected to be a dominant systematic in future accelerator neutrino experiments. The cross sections are determined by the linear response of the nucleus to the weak interactions of the neutrino, and are dominated by energy and distance scales of the order of the separation between nucleons in the nucleus. These response functions are potentially an important early physics application of quantum computers. Here we present an analysis of the resources required and their expected scaling for scattering cross section calculations. The current estimates of Trotter steps needed to achieve an energy resolution of 10 MeV and the number of CNOT gates for analyzing Ar40 highlights the need for significant improvements in algorithms. We also examine simple small-scale neutrino-nucleus models on modern quantum hardware. In this paper, we use variational methods to obtain the ground state of a three nucleon system (the triton) and then implement the relevant time evolution. To tame the errors in present-day NISQ devices, we explore the use of different error-mitigation techniques to increase the fidelity of the calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Finite-volume pionless effective field theory for few-nucleon systems with differentiable programming

Finite-volume pionless effective field theory provides an efficient framework for the extrapolation of nuclear spectra and matrix elements calculated at finite volume in lattice QCD to infinite volume, and to nuclei with larger atomic number. In this work, it is demonstrated how this framework may be implemented via a set of correlated Gaussian wave functions optimized using differentiable programming and via solution of a generalized eigenvalue problem. This approach is shown to be significantly more efficient than a stochastic implementation of the variational method based on the same form of correlated Gaussian wave functions, yielding comparably accurate representations of the ground-state wave functions with an order of magnitude fewer terms. The efficiency of representation allows such calculations to be extended to larger systems than in previous work. Further, the method is demonstrated through calculations of the binding energies of nuclei with atomic number A ϵ {2,3,4} in finite volume, matched to lattice QCD calculations at quark masses corresponding to m π = 806 MeV, and infinite-volume effective field theory calculations of A ϵ {2,3,4,5,6} systems based on this matching.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Automated label flows for excited states of correlation functions in lattice gauge theory

Extracting excited states from lattice gauge theory correlation functions can be achieved through $χ^2$ minimization fits or algebraic approaches such as the variational method and Prony's method. Performing any kind of error analysis often leads to overlapping confidence regions of model parameters, even when the spectrum is not particularly dense. To correctly estimate errors, one must beware of mislabeling the states. In this work, we provide an algorithm that we call automated label flows which consistently and systematically identifies a deterministic labeling of states. This is a black-box approach in the sense that it gives a sensible set of labels without user guidance. Further, as an example, we pair one black box method with another, analyzing a lattice correlation function from real data using automated label flows in the context of Prony's method.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2012

Comparison of Spectroscopic Data with Cluster Calculations of Plutonium, Plutonium Dioxide and Uranium Dioxide

Using spectroscopic data produced in the experimental investigations of bulk systems, including X-Ray Absorption Spectroscopy (XAS), Photoelectron Spectroscopy (PES) and Bremstrahlung Isochromat Spectroscopy (BIS), the theoretical results within for UO{sub 2}{sup 6}, PuO{sub 2}{sup 6} and Pu{sup 7} clusters have been evaluated. The calculations of the electronic structure of the clusters have been performed within the framework of the Relativistic Discrete-Variational Method (RDV). The comparisons between the LLNL experimental data and the Russian calculations are quite favorable. The cluster calculations may represent a new and useful avenue to address unresolved questions within the field of actinide electron structure, particularly that of Pu. Observation of the changes in the Pu electronic structure as a function of size suggests interesting implications for bulk Pu electronic structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

SU-E-T-147: Beam Specific Planning Target Volumes Incorporating 4DCT for Pencil Beam Scanning Proton Therapy of Thoracic Tumors

Purpose: The purpose of this study is to determine whether organ sparing and target coverage can be simultaneously maintained for pencil beam scanning (PBS) proton therapy treatment of thoracic tumors in the presence of motion, stopping power uncertainties and patient setup variations. Methods: Ten consecutive patients that were previously treated with proton therapy to 66.6/1.8 Gy (RBE) using double scattering (DS) were replanned with PBS. Minimum and maximum intensity images from 4DCT were used to introduce flexible smearing in the determination of the beam specific PTV (BSPTV). Datasets from eight 4DCT phases, using ±3% uncertainty in stopping power, and ±3 mm uncertainty in patient setup in each direction were used to create 8*12*10=960 PBS plans for the evaluation of ten patients. Plans were normalized to provide identical coverage between DS and PBS. Results: The average lung V20, V5, and mean doses were reduced from 29.0%, 35.0%, and 16.4 Gy with DS to 24.6%, 30.6%, and 14.1 Gy with PBS, respectively. The average heart V30 and V45 were reduced from 10.4% and 7.5% in DS to 8.1% and 5.4% for PBS, respectively. Furthermore, the maximum spinal cord, esophagus and heart dose were decreased from 37.1 Gy, 71.7 Gy and 69.2 Gy with DS to 31.3 Gy, 67.9 Gy and 64.6 Gy with PBS. The conformity index (CI), homogeneity index (HI), and global maximal dose were improved from 3.2, 0.08, 77.4 Gy with DS to 2.8, 0.04 and 72.1 Gy with PBS. All differences are statistically significant, with p values <0.05, with the exception of the heart V45 (p= 0.146). Conclusion: PBS with BSPTV achieves better organ sparing and improves target coverage using a repainting method for the treatment of thoracic tumors. Incorporating motion-related uncertainties is essential This work was supported by the US Army Medical Research and Materiel Command under Contract Agreement No. DAMD17-W81XWH-07-2-0121 and W81XWH-09-2-0174.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Witnessing eigenstates for quantum simulation of Hamiltonian spectra

The efficient calculation of Hamiltonian spectra, a problem often intractable on classical machines, can find application in many fields, from physics to chemistry. We introduce the concept of an “eigenstate witness” and, through it, provide a new quantum approach that combines variational methods and phase estimation to approximate eigenvalues for both ground and excited states. This protocol is experimentally verified on a programmable silicon quantum photonic chip, a mass-manufacturable platform, which embeds entangled state generation, arbitrary controlled unitary operations, and projective measurements. Both ground and excited states are experimentally found with fidelities >99%, and their eigenvalues are estimated with 32 bits of precision. We also investigate and discuss the scalability of the approach and study its performance through numerical simulations of more complex Hamiltonians. This result shows promising progress toward quantum chemistry on quantum computers.

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

Detection of Sodium in Highly Pure Graphite via High-Resolution Electrothermal Atomic Absorption Spectrometry with a Continuous Spectrum Source

This study is dedicated to elucidating the analytical abilities of electrothermal atomic absorption spectrometry with a continuous spectrum source (ETAAS-CSS) to detect the presence of sodium in highly pure graphite powder. In order to plot the calibration dependence for sodium detection, a special technique based on argon dilution of the atomic vapor obtained by introducing the aqueous reference solution in the atomizer is applied. The conditions for thermal pretreatment of highly pure graphite powder and atomization of sodium are established, and the sodium detection threshold is evaluated via ETAAS-CSS (2.6 × 10{sup –4} ng) as well. The validity of the results is confirmed via the sample variation method.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Isotropic and anisotropic dark energy models

In this review we discuss the evolution of the universe filled with dark energy with or without perfect fluid. In doing so we consider a number of cosmological models, namely Bianchi type I, III, V, VI{sub 0}, VI and FRW ones. For the anisotropic cosmological models we have used proportionality condition as an additional constrain. The exact solutions to the field equations in quadrature are found in case of a BVI model. It was found that the proportionality condition used here imposed severe restriction on the energy-momentum tensor, namely it leads to isotropic distribution of matter. Anisotropic BVI{sub 0}, BV, BIII and BIDE models with variable EoS parameter ω have been investigated by using a law of variation for the Hubble parameter. In this case the matter distribution remains anisotropic, though depending on the concrete model there appear different restrictions on the components of energy-momentum tensor. That is why we need an extra assumption such as variational a law for the Hubble parameter. It is observed that, at the early stage, the EoS parameter v is positive i.e. the universe was matter dominated at the early stage but at later time, the universe is evolving with negative values, i.e., the present epoch. DE model presents the dynamics of EoS parameter ω whose range is in good agreement with the acceptable range by the recent observations. A spatially homogeneous and anisotropic locally rotationally symmetric Bianchi-I space time filled with perfect fluid and anisotropic DE possessing dynamical energy density is studied. In the derived model, the EoS parameter of DE (ω{sup (de)}) is obtained as time varying and it is evolving with negative sign which may be attributed to the current accelerated expansion of Universe. The distance modulus curve of derived model is in good agreement with SNLS type Ia supernovae for high redshift value which in turn implies that the derived model is physically realistic. A system of two fluids within the scope of a spatially flat and isotropic FRW model is studied. The role of the two fluids, either minimally or directly coupled in the evolution of the dark energy parameter, has been investigated. In doing so we have used three different ansatzs regarding the scale factor that gives rise to a variable decelerating parameter. It is observed that, in the non-interacting case, both the open and flat universes can cross the phantom region whereas in the interacting case only the open universe can cross the phantom region. The stability and acceptability of the obtained solution are also investigated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Biexcitons formed from spatially separated electrons and holes in quasi-zero-dimensional semiconductor nanosystems

A theory of biexcitons (formed from spatially separated electron and holes) in nanosystems that consist of zinc-selenide quantum dots synthesized in borosilicate glassy matrices is developed. The dependences of the total energy and the binding energy of the singlet ground biexciton state in such a system on the spacing between the quantum-dot surfaces and the quantum-dot radius are derived by the variational method. It is shown that biexciton formation is of the threshold character and possible in nanosystems, in which the spacing between the quantum-dot surfaces is larger than a certain critical spacing.

36 MATERIALS SCIENCE↗
DOE OSTI2012

High-Order Curvilinear Finite Element Methods for Lagrangian Hydrodynamics [High Order Curvilinear Finite Elements for Lagrangian Hydrodynamics]

The numerical approximation of the Euler equations of gas dynamics in a movingLagrangian frame is at the heart of many multiphysics simulation algorithms. Here, we present a general framework for high-order Lagrangian discretization of these compressible shock hydrodynamics equations using curvilinear finite elements. This method is an extension of the approach outlined in [Dobrev et al., Internat. J. Numer. Methods Fluids, 65 (2010), pp. 1295--1310] and can be formulated for any finite dimensional approximation of the kinematic and thermodynamic fields, including generic finite elements on two- and three-dimensional meshes with triangular, quadrilateral, tetrahedral, or hexahedral zones. We discretize the kinematic variables of position and velocity using a continuous high-order basis function expansion of arbitrary polynomial degree which is obtained via a corresponding high-order parametric mapping from a standard reference element. This enables the use of curvilinear zone geometry, higher-order approximations for fields within a zone, and a pointwise definition of mass conservation which we refer to as strong mass conservation. Moreover, we discretize the internal energy using a piecewise discontinuous high-order basis function expansion which is also of arbitrary polynomial degree. This facilitates multimaterial hydrodynamics by treating material properties, such as equations of state and constitutive models, as piecewise discontinuous functions which vary within a zone. To satisfy the Rankine--Hugoniot jump conditions at a shock boundary and generate the appropriate entropy, we introduce a general tensor artificial viscosity which takes advantage of the high-order kinematic and thermodynamic information available in each zone. Finally, we apply a generic high-order time discretization process to the semidiscrete equations to develop the fully discrete numerical algorithm. Our method can be viewed as the high-order generalization of the so-called staggered-grid hydrodynamics (SGH) approach and we show that under specific low-order assumptions, we exactly recover the classical SGH method. We also present numerical results from an extensive series of verification tests that demonstrate several important practical advantages of using high-order finite elements in this context.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗