Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variational approach”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A variational Monte Carlo approach for core excitations

In this work, we present a systematically improvable approach to core excitations in variational Monte Carlo. Building on recent work in excited-state-specific Monte Carlo, we show how a straightforward protocol, starting from a quantum chemistry guess, is able to capture core state's strong orbital relaxations, maintain accuracy in the near-nuclear region during these relaxations, and explicitly balance accuracy between ground and core excited states. In water, ammonia, and methane, which serve as prototypical representatives for oxygen, nitrogen, and carbon core states, respectively, this approach predicts core excitation energies within 0.3 eV of experiment and core excitation peak separations within 0.1 eV of experiment.

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

Constructing nested coordinates inside strongly shaped toroids using an action principle

A new approach for constructing polar-like boundary-conforming coordinates inside a toroid with strongly shaped cross-sections is presented. A coordinate mapping is obtained through a variational approach, which involves identifying extremal points of a proposed action in the mapping space from [0,2π] 2 ×[0,1] to a toroidal domain in $\mathbb{R}$ 3 . This approach employs an action built on the squared Jacobian and radial length. Extensive testing is conducted on general toroidal boundaries using a global Fourier–Zernike basis via action minimisation. The results demonstrate successful coordinate construction capable of accurately describing strongly shaped toroidal domains. The coordinate construction is successfully applied to the computation of three-dimensional magnetohydrodynamic equilibria in the GVEC code where the use of traditional coordinate construction by interpolation from the boundary failed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

A new approach for magnetic curves in 3D Riemannian manifolds

A magnetic field is defined by the property that its divergence is zero in a three-dimensional oriented Riemannian manifold. Each magnetic field generates a magnetic flow whose trajectories are curves called as magnetic curves. In this paper, we give a new variational approach to study the magnetic flow associated with the Killing magnetic field in a three-dimensional oriented Riemann manifold, (M{sup 3}, g). And then, we investigate the trajectories of the magnetic fields called as N-magnetic and B-magnetic curves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Renormalization of excitonic properties by polar phonons

Here we employ quasiparticle path integral molecular dynamics to study how the excitonic properties of model semiconductors are altered by electron–phonon coupling. We describe ways within a path integral representation of the system to evaluate the renormalized mass, binding energy, and radiative recombination rate of excitons in the presence of a fluctuating lattice. To illustrate this approach, we consider Fröhlich-type electron–phonon interactions and employ an imaginary time influence functional to incorporate phonon-induced effects nonperturbatively. The effective mass and binding energies are compared with perturbative and variational approaches, which provide qualitatively consistent trends. We evaluate electron-hole recombination rates as mediated through both trap-assisted and bimolecular processes, developing a consistent statistical mechanical approach valid in the reaction limited regime. These calculations demonstrate how phonons screen electron–hole interactions, generically reducing exciton binding energies and increasing their radiative lifetimes.

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

Radiative decay rate of excitons in square quantum wells: Microscopic modeling and experiment

The binding energy and the corresponding wave function of excitons in GaAs-based finite square quantum wells (QWs) are calculated by the direct numerical solution of the three-dimensional Schrödinger equation. The precise results for the lowest exciton state are obtained by the Hamiltonian discretization using the high-order finite-difference scheme. The microscopic calculations are compared with the results obtained by the standard variational approach. The exciton binding energies found by two methods coincide within 0.1 meV for the wide range of QW widths. The radiative decay rate is calculated for QWs of various widths using the exciton wave functions obtained by direct and variational methods. The radiative decay rates are confronted with the experimental data measured for high-quality GaAs/AlGaAs and InGaAs/GaAs QW heterostructures grown by molecular beam epitaxy. The calculated and measured values are in good agreement, though slight differences with earlier calculations of the radiative decay rate are observed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Nonlinear saturation of ballooning modes in stellarators

Ballooning mode saturation is investigated in realistic stellarator configurations using the flux tube approach of Ham et al (2018 Plasma Phys. Control. Fusion 60 075017), Ham et al (2016 Phys. Rev. Lett. 116 235001). The method is adapted to account for the lack of exact force balance in stellarator equilibrium solvers that assume existence of nested flux surfaces. A variational approach for calculating flux tube energy is developed to overcome this force error problem in stellarator numerical equilibria. Saturated (equilibrium) flux tube states that cross 10%–20% of the plasma minor radius are shown to exist for linearly ballooning unstable profiles. It is shown that several features of the displaced flux tube structure in a full nonlinear MHD simulation of Wendelstein 7X are reproduced by our model. Saturated states are found in a compact stellarator equilibrium close but below the marginal ballooning linear instability, i.e. the unperturbed equilibrium is metastable. This suggests that edge-localized-mode-like explosive MHD behavior may be possible in stellarators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Solitary waves in the nonlinear Dirac equation in the presence of external driving forces

In this paper, we consider the nonlinear Dirac (NLD) equation in (1 + 1) dimensions with scalar–scalar self interaction g 2 /κ + 1 (Ψ¯Ψ) κ + 1 in the presence of external forces as well as damping of the form f(x) - iμγ 0 Ψ, where both f and Ψ are two-component spinors. We develop an approximate variational approach using collective coordinates (CC) for studying the time dependent response of the solitary waves to these external forces. This approach predicts intrinsic oscillations of the solitary waves, i.e. the amplitude, width and phase all oscillate with the same frequency. The translational motion is also affected, because the soliton position oscillates around a mean trajectory. For κ = 1 we solve explicitly the CC equations of the variational approximation for slow moving solitary waves in a constant external force without damping and find reasonable agreement with solving numerically the CC equations. Finally, we then compare the results of the variational approximation with no damping with numerical simulations of the NLD equation for κ = 1, when the components of the external force are of the form f j = r j exp(–iΚx) and again find agreement if we take into account a certain linear excitation with specific wavenumber that is excited together with the intrinsic oscillations such that the momentum in a transformed NLD equation is conserved.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

A velocity/stress mixed stabilized nodal finite element for elastodynamics: Analysis and computations with strongly and weakly enforced boundary conditions

A new nodal mixed finite element is proposed for the simulation of linear elastodynamics and wave propagation problems in time domain. Our method is based on equal-order interpolation discrete spaces for both the velocity (or displacement) and stress (or strain) tensor variables. The mixed form is derived using either the velocity/stress or velocity/strain pair of unknowns, the latter being instrumental in extensions of the method to nonlinear mechanics. The proposed approach works equally well on hexahedral or tetrahedral grids and, for this reason, it is suitable for time-domain engineering applications in complex geometry. The peculiarity of the proposed approach is the use of the rate form of the stress update equation, which yields a set of governing equations with the structure of a non-dissipative space/time Friedrichs’ system. We complement standard traction boundary conditions for the stress with strongly and weakly enforced boundary conditions for the velocity (or displacement). Weakly enforced boundary conditions are particularly suitable when considering complex geometrical shapes, because they do not require dedicated data structures for the imposition of the boundary degrees of freedom, but, rather, they utilize the structure of the variational formulation. Here we also show how the framework of weakly enforced boundary conditions can be used to develop variational forms for multi-domain simulations of heterogeneous media. A complete analysis including stability and convergence proofs is included, in the case of a space–time variational approach. A series of computational tests are used to demonstrate and verify the performance of the proposed approach.

42 ENGINEERING↗
NASA NTRS1997

An MHD variational principle that admits reconnection

The variational approach of Pfirsch and Sudan's averaged magnetohydrodynamics (MHD) to the stability of a line-tied current layer is summarized. The effect of line-tying on current sheets that might arise in line-tied magnetic flux tubes by estimating the growth rates of a resistive instability using a variational method. The results show that this method provides a potentially new technique to gauge the stability of nearly ideal magnetohydrodynamic systems. The primary implication for the stability of solar coronal structures is that tearing modes are probably constant at work removing magnetic shear from the solar corona.

Rilee, M. L.↗
DOE OSTI2015

Quantum mechanical derivation of the Wallis formula for π

A famous pre-Newtonian formula for π is obtained directly from the variational approach to the spectrum of the hydrogen atom in spaces of arbitrary dimensions greater than one, including the physical three dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Compact submanifolds supporting singular interactions

A quantum particle moving under the influence of singular interactions on embedded surfaces furnish an interesting example from the spectral point of view. In these problems, the possible occurrence of a bound-state is perhaps the most important aspect. Such systems can be introduced as quadratic forms and generically they do not require renormalization. Yet an alternative path through the resolvent is also beneficial to study various properties. In the present work, we address these issues for compact surfaces embedded in a class of ambient manifolds. We discover that there is an exact bound state solution written in terms of the heat kernel of the ambient manifold for a range of coupling strengths. Moreover, we develop techniques to estimate bounds on the ground state energy when several surfaces, each of which admits a bound state solution, coexist. -- Highlights: •Schrödinger operator with singular interactions supported on compact submanifolds. •Exact bound-state solution in terms of the heat kernel of the ambient manifold. •Generalization of the variational approach to a collection of submanifolds. •Existence of a lower bound for a unique ground state energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Classical and quantum simulations of 1+1-dimensional ${\mathbb{Z}}_{2}$ gauge theory at finite temperature and density

Simulating strongly coupled gauge theories at finite temperature and density is a longstanding challenge in nuclear and high-energy physics with fundamental implications for condensed matter physics. Here, we simulate such systems using minimally entangled typical thermal state (METTS) approaches, which combine classical random sampling with imaginary-time evolution, implementable on either classical or quantum computers, to estimate thermal averages of observables. We study 1+1-dimensional ${\mathbb{Z}}_{2}$ gauge theory coupled to spinless fermionic matter, which maps onto a local quantum spin chain. We benchmark both a classical matrix-product-state implementation of METTS and a recently proposed adaptive variational approach for near-term quantum devices, focusing on the equation of state and measures of fermion confinement. Of particular importance is the choice of basis for METTS sampling, which impacts both the sampling overhead and quantum circuit complexity. Our work sets the stage for future studies of strongly coupled gauge theories using classical and quantum hardware.

Chen, I-Chi [Iowa State Univ., Ames, IA (United St↗
DOE OSTI2015

A theoretical method to compute sequence dependent configurational properties in charged polymers and proteins

A general formalism to compute configurational properties of proteins and other heteropolymers with an arbitrary sequence of charges and non-uniform excluded volume interaction is presented. A variational approach is utilized to predict average distance between any two monomers in the chain. The presented analytical model, for the first time, explicitly incorporates the role of sequence charge distribution to determine relative sizes between two sequences that vary not only in total charge composition but also in charge decoration (even when charge composition is fixed). Furthermore, the formalism is general enough to allow variation in excluded volume interactions between two monomers. Model predictions are benchmarked against the all-atom Monte Carlo studies of Das and Pappu [Proc. Natl. Acad. Sci. U. S. A. 110, 13392 (2013)] for 30 different synthetic sequences of polyampholytes. These sequences possess an equal number of glutamic acid (E) and lysine (K) residues but differ in the patterning within the sequence. Without any fit parameter, the model captures the strong sequence dependence of the simulated values of the radius of gyration with a correlation coefficient of R{sup 2} = 0.9. The model is then applied to real proteins to compare the unfolded state dimensions of 540 orthologous pairs of thermophilic and mesophilic proteins. The excluded volume parameters are assumed similar under denatured conditions, and only electrostatic effects encoded in the sequence are accounted for. With these assumptions, thermophilic proteins are found—with high statistical significance—to have more compact disordered ensemble compared to their mesophilic counterparts. The method presented here, due to its analytical nature, is capable of making such high throughput analysis of multiple proteins and will have broad applications in proteomic studies as well as in other heteropolymeric systems.

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

Dynamics of the sub-Ohmic spin-boson model: A time-dependent variational study

The Dirac-Frenkel time-dependent variation is employed to probe the dynamics of the zero temperature sub-Ohmic spin-boson model with strong friction utilizing the Davydov D{sub 1} ansatz. It is shown that initial conditions of the phonon bath have considerable influence on the dynamics. Counterintuitively, even in the very strong coupling regime, quantum coherence features still manage to survive under the polarized bath initial condition, while such features are absent under the factorized bath initial condition. In addition, a coherent-incoherent transition is found at a critical coupling strength {alpha} Almost-Equal-To 0.1 for s= 0.25 under the factorized bath initial condition. We quantify how faithfully our ansatz follows the Schroedinger equation, finding that the time-dependent variational approach is robust for strong dissipation and deep sub-Ohmic baths (s Much-Less-Than 1).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Excitons and Trions in Bilayer van der Waals Heterostructures

Excitons and trions in bilayer structures based on transition metal dichalcogenides are theoretically studied. Expressions for the effective potential of interparticle interaction in such a system with allowance for a significant dielectric contrast in the heterostructure are obtained. Simple and physically justified variational functions for electron–hole complexes are proposed. A variational calculation of the exciton and trion binding energies as functions of the interlayer distance, taking into account the specific features of the screening of the Coulomb interaction, is performed. The accuracy of the variational approach is confirmed by direct numerical calculation of the exciton binding energy as a function of the interlayer distance.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Dynamical stability of a many-body Kapitza pendulum

We consider a many-body generalization of the Kapitza pendulum: the periodically-driven sine–Gordon model. We show that this interacting system is dynamically stable to periodic drives with finite frequency and amplitude. This finding is in contrast to the common belief that periodically-driven unbounded interacting systems should always tend to an absorbing infinite-temperature state. The transition to an unstable absorbing state is described by a change in the sign of the kinetic term in the Floquet Hamiltonian and controlled by the short-wavelength degrees of freedom. We investigate the stability phase diagram through an analytic high-frequency expansion, a self-consistent variational approach, and a numeric semiclassical calculation. Classical and quantum experiments are proposed to verify the validity of our results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

High-Resolution, Rapid-Scan Dual-Doppler Retrievals of Vertical Velocity in a Simulated Supercell

Observation system simulation experiments are used to evaluate different dual-Doppler analysis (DDA) methods for retrieving vertical velocity w at grid spacings on the order of 100 m within a simulated tornadic supercell. Variational approaches with and without a vertical vorticity equation constraint are tested, along with a typical (traditional) method involving vertical integration of the mass conservation equation. The analyses employ emulated radar data from dual-Doppler placements 15, 30, and 45 km east of the mesocyclone, with volume scan intervals ranging from 10 to 150 s. The effect of near-surface data loss is examined by denying observations below 1 km in some of the analyses. At the longer radar ranges and when no data denial is imposed, the “traditional” method produces results similar to those of the variational method and is much less expensive to implement. However, at close range and/or with data denial, the variational method is much more accurate, confirming results from previous studies. The vorticity constraint shows the potential to improve the variational analysis substantially, reducing errors in the w retrieval by up to 30% for rapid-scan observations (≤30 s) at close range when the local vorticity tendency is estimated using spatially variable advection correction. However, the vorticity constraint also degrades the analysis for longer scan intervals, and the impact diminishes with increased range. Furthermore, analyses using 30-s data also frequently outperform analyses using 10-s data, suggesting a limit to the benefit of increasing the radar scan rate for variational DDA employing the vorticity constraint.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

Distorting general relativity: gravity's rainbow and f(R) theories at work

We compute the Zero Point Energy in a spherically symmetric background combining the high energy distortion of Gravity's Rainbow with the modification induced by a f(R) theory. Here f(R) is a generic analytic function of the Ricci curvature scalar R in 4D and in 3D. The explicit calculation is performed for a Schwarzschild metric. Due to the spherically symmetric property of the Schwarzschild metric we can compare the effects of the modification induced by a f(R) theory in 4D and in 3D. We find that the final effect of the combined theory is to have finite quantities that shift the Zero Point Energy. In this context we setup a Sturm-Liouville problem with the cosmological constant considered as the associated eigenvalue. The eigenvalue equation is a reformulation of the Wheeler-DeWitt equation which is analyzed by means of a variational approach based on gaussian trial functionals. With the help of a canonical decomposition, we find that the relevant contribution to one loop is given by the graviton quantum fluctuations around the given background. A final discussion on the connection of our result with the observed cosmological constant is also reported.

79 ASTRONOMY AND ASTROPHYSICS↗