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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Few-electron atomic ions in non-relativistic QED: The ground state

Following detailed analysis of relativistic, QED and mass corrections for helium-like and lithium-like ions with static nuclei for Z≤20 the domain of applicability of Non-Relativistic QED (NRQED) is localized for ground state energy. It is demonstrated that for both helium-like and lithium-like ions with Z≤20 the finite nuclear mass effects do not change 4–5 significant digits (s.d.), and the leading relativistic and QED effects leave unchanged 3–4 s.d. in the ground state energy. It is shown that the non-relativistic ground state energy can be interpolated with accuracy not less than 13 s.d. for Z≤12, and not less than 12 s.d. for Z≤50 for helium-like as well as for Z≤20 for lithium-like ions by a compact meromorphic function in variable λ=√(Z−Z{sub B}) (here Z{sub B} is the 2nd critical charge (Turbiner et al., 2016)), P{sub 9}(λ)∕Q{sub 5}(λ). It is found that both the Majorana formula – a second degree polynomial in Z with two free parameters – and a fourth degree polynomial in λ (a generalization of the Majorana formula) reproduce the ground state energy of the helium-like and lithium-like ions for Z≤20 in the domain of applicability of NRQED, thus, at least, 3 s.d. It is noted that ≳99.9% of the ground state energy is given by the variational energy for properly optimized trial function of the form of (anti)-symmetrized product of three (six) screened Coulomb orbitals for two-(three) electron system with 3 (7) free parameters for Z≤20, respectively. It may imply that these trial functions are, in fact, exact wavefunctions in non-relativistic QED, thus, the NRQED effective potential can be derived. It is shown that the sum of relativistic and QED effects in leading approximation - 3 s.d. - for both 2 and 3 electron systems is interpolated by 4th degree polynomial in Z for Z≤20.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Distribution of solutions of the fastest apparent convergence condition in optimized perturbation theory and its relation to anti-Stokes lines

We discuss fundamental properties of the fastest apparent convergence (FAC) condition which is used as a variational criterion in optimized perturbation theory (OPT). We examine an integral representation of the FAC condition and a distribution of the zeros of the integral in a complex artificial parameter space on the basis of theory of Lefschetz thimbles. We find that the zeros accumulate on a certain line segment known as an anti-Stokes line in the limit K→∞, where K is a truncation order of a perturbation series. This phenomenon gives an underlying mechanism that physical quantities calculated by OPT can be insensitive to the choice of the artificial parameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

A covariant simultaneous action for branes

Highlights: • A novel variational approach is presented. • It provides a shortcut towards the second variation of geometrical models. • This should be of interest, for its generality, to various areas of theoretical physics. - Abstract: A covariant simultaneous action for branes in an arbitrary curved background spacetime is considered. The term ‘simultaneous’ is imported from variational calculus and refers to the fact that extremization of the action produces at once both the first and second variations of a given geometrical action for the brane. The action depends on a pair of independent field variables, the brane embedding functions, through the canonical momentum of a reparametrization invariant geometric model for the brane, and an auxiliary vector field. The form of the action is analogous to a symplectic potential. Extremization of the simultaneous action produces at once the equations of motion and the Jacobi equations for the brane geometric model, and it also provides a convenient shortcut towards its second variation. In this note, we consider geometric models depending only on the intrinsic geometry of the brane worldvolume, and discuss briefly the generalization to extrinsic geometry dependent models. The approach is illustrated for Dirac–Nambu–Goto [DNG] branes. For a relativistic particle, a simultaneous action was introduced by Ba z-dot ański, that served as an inspiration for the present work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Helicity, spin, and infra-zilch of light: A Lorentz covariant formulation

Highlights: • Spin and orbital parts of angular momentum. • The helicity tensor contains the helicity, spin and infra-zilch currents. • The helicity tensor is a Noether current for a duality-symmetric action. • Symmetry generator for the helicity tensor. • The helicity tensor is conserved in Lorenz gauge. In this paper, a novel conserved Lorentz covariant tensor, termed the helicity tensor, is introduced in Maxwell theory. The conservation of the helicity tensor expresses the conservation laws contained in the helicity array, introduced by Cameron et al. (2012), including helicity, spin, as well as the spin-flux or infra-zilch. The Lorentz covariance of the helicity tensor is in contrast to previous formulations of the helicity hierarchy of conservation laws, which required the non-Lorentz covariant transverse gauge. The helicity tensor is shown to arise as a Noether current for a variational symmetry of a duality-symmetric Lagrangian for Maxwell theory. This symmetry transformation generalizes the duality symmetry and includes the symmetry underlying the conservation of the spin part of the angular momentum.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Attaining regularization length insensitivity in phase-field models of ductile failure

In this work, a cohesive phase-field model of ductile fracture in a finite-deformation setting is presented. The model is based on a free-energy function in which both elastic and plastic work contributions are coupled to damage. Using a strictly variational framework, the field evolution equations, damage kinetics, and flow rule are jointly derived from a scalar least-action principle. Particular emphasis is placed on the use of a rational function for the stress degradation that maintains a fixed effective strength with decreasing regularization length. The model is employed to examine crack growth in pure mode-I problems through the generation of crack growth resistance (J-R) curves. In contrast to alternative models, the current formulation gives rise to J-R curves that are insensitive to the regularization length. Numerical evidence suggests convergence of local fields with respect to diminishing regularization length as well.

42 ENGINEERING↗
DOE OSTI2024

Scalable semidefinite programming approach to variational embedding for quantum many-body problems

In quantum embedding theories, a quantum many-body system is divided into localized clusters of sites which are treated with an accurate ‘high-level’ theory and glued together self-consistently by a less accurate ‘low-level’ theory at the global scale. The recently introduced variational embedding approach for quantum many-body problems combines the insights of semidefinite relaxation and quantum embedding theory to provide a lower bound on the ground-state energy that improves as the cluster size is increased. The variational embedding method is formulated as a semidefinite program (SDP), which can suffer from poor computational scaling when treated with black-box solvers. Here, we exploit the interpretation of this SDP as an embedding method to develop an algorithm which alternates parallelizable local updates of the high-level quantities with updates that enforce the low-level global constraints. Moreover, we show how translation invariance in lattice systems can be exploited to reduce the complexity of projecting a key matrix to the positive semidefinite cone.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Noether currents for Eulerian variational principles in non-barotropic magnetohydrodynamics and topological conservations laws

We derive a Noether current for the Eulerian variational principle of ideal non-barotropic magnetohydrodynamics (MHD). It was shown previously that ideal non-barotropic MHD is mathematically equivalent to a five function field theory with an induced geometrical structure in the case that field lines cover surfaces and this theory can be described using a variational principle. Here we use various symmetries of the flow to derive topological constants of motion through the derived Noether current and discuss their implication for non-barotropic MHD.

42 ENGINEERING↗
DOE OSTI2020

Cross-Sections of Nanocellulose from Wood Analyzed by Quantized Polydispersity of Elementary Microfibrils

Bio-based nanocellulose has been shown to possess impressive mechanical properties and simplicity for chemical modifications. The chemical properties are largely influenced by the surface area and functionality of the nanoscale materials. However, finding the typical cross-sections of nanocellulose, such as cellulose nanofibers (CNFs), has been a long-standing puzzle, where subtle changes in extraction methods seem to yield different shapes and dimensions. Here, we extracted CNFs from wood with two different oxidation methods and variations in degree of oxidation and high-pressure homogenization. Additionally, the cross-sections of CNFs were characterized by small-angle X-ray scattering and wide-angle X-ray diffraction in dispersed and freeze-dried states, respectively, where the results were analyzed by assuming that the cross-sectional distribution was quantized with an 18-chain elementary microfibril, the building block of the cell wall. We find that the results agree well with a pseudosquare unit having a size of about 2.4 nm regardless of sample, while the aggregate level strongly depends on the extraction conditions. Furthermore, we find that aggregates have a preferred cohesion of phase boundaries parallel to the (110)-plane of the cellulose fibril, leading to a ribbon shape on average.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Emulation With Uncertainty Quantification of Regional Sea‐Level Change Caused by the Antarctic Ice Sheet

Abstract Projecting regional sea‐level change under various climate‐change scenarios typically involves running forward simulations of the Earth's gravitational, rotational and deformational (GRD) response to ice‐mass change, which requires substantial computational cost if applied to probabilistic frameworks requiring thousands to millions of samples. Here we build emulators of regional sea‐level change at 27 coastal locations, due to the GRD effects associated with future Antarctic Ice Sheet mass change over the 21st century. The emulators are evaluated against a numerical sea‐level model applied to an ensemble of ice‐sheet model simulations of the Antarctic Ice Sheet through 2100. We build a physics‐based emulator using a recent sensitivity kernel approach and compare it to machine learning based emulators (neural network and conditional variational autoencoder methods). In order to quantify uncertainty, we derive well‐calibrated prediction intervals for regional sea‐level change via split‐conformal inference and linear regression, and show that Monte Carlo dropout does not yield well‐calibrated uncertainties in this instance. We also demonstrate substantial gains in computational efficiency using both the physics‐based emulator and neural networks in comparison to the numerical model for the complete regional sea‐level solution. Overall, we find the physics‐based emulator modestly outperforms the machine learning emulators for this problem.

58 GEOSCIENCES↗
DOE OSTI2022

Generalization in quantum machine learning from few training data

Modern quantum machine learning (QML) methods involve variationally optimizing a parameterized quantum circuit on a training data set, and subsequently making predictions on a testing data set (i.e., generalizing). In this work, we provide a comprehensive study of generalization performance in QML after training on a limited number N of training data points. We show that the generalization error of a quantum machine learning model with T trainable gates scales at worst as $\sqrt{T/N}$. When only K$\ll$T gates have undergone substantial change in the optimization process, we prove that the generalization error improves to $\sqrt{K/N}$. Our results imply that the compiling of unitaries into a polynomial number of native gates, a crucial application for the quantum computing industry that typically uses exponential-size training data, can be sped up significantly. We also show that classification of quantum states across a phase transition with a quantum convolutional neural network requires only a very small training data set. Other potential applications include learning quantum error correcting codes or quantum dynamical simulation. Our work injects new hope into the field of QML, as good generalization is guaranteed from few training data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Evaluating the Limits of QAOA Parameter Transfer at High-Rounds on Sparse Ising Models With Geometrically Local Cubic Terms

The emergent practical applicability of the Quantum Approximate Optimization Algorithm (QAOA) for approximate combinatorial optimization is a subject of considerable interest. One of the primary limitations of QAOA is the task of finding a set of good parameters, which is usually done using a variational optimization loop. Parameter transfer, or parameter concentration, is a phenomenon where QAOA angles trained on problem instances that are self-similar tend to perform well for other problem instances from that similar class. This suggests a potentially highly efficient and scalable non-variational learning method for QAOA angle finding. In this work, we systematically study QAOA parameter transferability from small problem sizes (16 and 27 decision variables) onto large problem instances (up to 156 qubits) for heavy-hex graph Ising models with geometrically local higher order terms using the Julia based QAOA simulation tool \texttt{JuliQAOA} to perform classical angle finding for up to $49$ QAOA layers ($p$). Parameter transfer of the fixed angles is validated using a combination of full statevector, Projected Entangled Pair States (PEPS), Matrix Product State (MPS), and LOWESA numerical simulations. We find that the QAOA parameter transfer from single instances applied to other (unseen) problem instances does not in general provide monotonically improving performance as a function of $p$ - there are many cases where the performance temporarily decreases as a function of $p$ - but despite this the transferred angles have a general trend of improved expectation value as the QAOA depth increases, in many cases converging close to the true ground-state energy of the $100+$ qubit instances. We also sample the hardware-compatible Ising models using the ensemble of transfer-learned QAOA parameters on several superconducting qubit IBM Quantum processors with 127, 133, and 156 qubits. We find continuous solution quality improvement of the hardware-compatible QAOA circuits run on the IBM NISQ processors up to $p=5$ on \texttt{ibm\_fez}, up to $p=9$ on \texttt{ibm\_torino}, and up to $p=10$ on \texttt{ibm\_pittsburgh}.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Neuralized fermionic tensor networks for quantum many-body systems

In this work, we describe a class of neuralized fermionic tensor network states (NN-fTNSs) that introduce nonlinearity into fermionic tensor networks through configuration-dependent neural network transformations of the local tensors. The construction uses the fTNS algebra to implement a natural fermionic sign structure and is compatible with standard tensor network algorithms but gains enhanced expressivity through the neural network parametrization. Using the 1D and 2D Fermi-Hubbard models as benchmarks, we demonstrate that NN-fTNSs achieve order of magnitude improvements in the ground-state energy compared to pure fTNSs with the same bond dimension and can be systematically improved through both the tensor network bond dimension and the neural network parametrization. Compared to existing fermionic neural quantum states based on Slater determinants and Pfaffians, NN-fTNSs offer a physically motivated alternative fermionic structure. Furthermore, compared to such states, NN-fTNSs naturally exhibit improved computational scaling and we demonstrate a construction that achieves linear scaling with the lattice size.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Electromagnetic Analysis of ITER Electron Cyclotron Emission Model

The ITER electron cyclotron emission (ECE) diagnostic system is located at diagnostics shielding module 2 (DSM2), equatorial port 9 (EP9), to measure electron temperature profile and electron temperature fluctuations, and also assess nonthermal electron distributions via the oblique view. Therefore, ECE has both radial and oblique views with two hot sources for calibration and two couples of mirrors for two different optical views. It shall receive that the large electromagnetic (EM) loads up to 100-MN/m3 force density due to the eddy currents on the copper mirrors generated by the short transient plasma disruption. The simplified EM analysis model using magnetic field ( B ) data and flux variations (dB/dt) method based on the worst case of plasma disruption, MD_DW_EXP16MS_CATIII, has been developed to calculate the force and moments for each ECE component and support structure for EP9 DSM2 Bay2 and Bay3. The B and dB/dt methods use the constant B and dB/dt assumed for each local bay space, but the real B and dB/dt are basically varied along the radial direction of each bay. The complex global model with the local ECE components and support structure using Maxwell transient has been developed to do EM analysis to compare the results of the B and dB/dt models. The volumetric EM force density of the whole ECE components and support structure can be used for the subsequent structural analysis. In conclusion, combined with the thermal and nuclear loads and seismic and initial loads, the ECE integration analysis can be finalized for the operation case.

Fang, J. [Princeton Plasma Physics Laboratory (PPP↗
DOE OSTI2018

Metastable Bound States of the Two-Dimensional Bimagnetoexcitons in the Lowest Landau Levels Approximation

The possible existence of the bound states of the interacting two-dimensional (2D) magnetoexcitons in the lowest Landau levels (LLLs) approximation was investigated using the Landau gauge description. The magnetoexcitons taking part in the formation of the bound state with resultant wave vector k=0 have opposite in-plane wave vectors k and −k and look as two electric dipoles with the arms oriented in-plane perpendicularly to the corresponding wave vectors. The bound state of two antiparallel dipoles moving with equal probability in any direction of the plane with equal but antiparallel wave vectors is characterized by the variational wave function of the relative motion φ{sub n}(k) depending on the modulus |k|. The spins of two electrons and the effective spins of two holes forming the bound states were combined separately in the symmetric or in the antisymmetric forms (↑↓ +η ↓↑) with the same parameter η=±1 for electrons and holes. In the case of the variational wave function φ{sub 2}(k)=(8α{sup 3}){sup 1/2}k{sup 2}l{sub 0}{sup 2}exp [−αk{sup 2}l{sub 0}{sup 2}] the maximum density of the magnetoexcitons in the momentum space representation is concentrated on the in-plane ring with the radius k{sub r}=1/(√αl{sub 0}). The stable bound states of the bimagnetoexciton molecule do not exist for both spin orientations. Instead of them, a deep metastable bound state with the activation barrier comparable with the ionization potential I{sub l} of the magnetoexciton with k=0 was revealed in the case η=1 and α=0.5. In the case η=−1 and α=3.4 only a shallow metastable bound state can appear.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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 OSTI2025

Data-Driven Clustering and Classification of Outage Patterns with Insights into their Links to Extreme Events

At a global level extreme events have increased in both scale and impact. These events have the potential to affect the electrical grid infrastructure and cause a wide range of outages, which can lead to a disruption in daily patterns, cost millions of dollars and also the loss of life. Currently, to track these outage events there have been various approaches developed ranging from regional to national level quantifications for what defines an outage. However, this variation in methods can potentially lead to subjective decision-making and a lack of proper management in relation to the event. While previous work has made strides in determining spatio-temporal patterns, minimal attention has been given to the type and number of outages an area may be exposed to. The differences in incurred cost and the overall severity of an event between a transformer box malfunction and a hurricane are drastic, and by finding historical signals, we can allow for more efficient management, potentially saving lives and millions of dollars. Here, we leverage unsupervised machine learning techniques to delineate outage patterns among 22 counties within the United States and find that there are clear, segregated clusters (0.93 silhouette) of data which are related by event behavior and underlying cause. This finding will allow for energy stakeholders, policy makers, and researchers to gain a deeper understanding of the extent and severity of historic events and to better prepare for electrical grid infrastructure planning and management.

Koob, Benjamin [ORNL]↗
DOE OSTI2022

The Impact of Constrained Data Assimilation on the Forecasts of Three Convection Systems During the ARM MC3E Field Campaign

A constrained data assimilation (CDA) system based on the ensemble variational (EnVar) method and physical constraints of mass and water conservations is evaluated through three convective cases during the Midlatitude Continental Convective Clouds Experiment (MC3E) of the Atmospheric Radiation Measurement (ARM) program. Compared to the original data assimilation (ODA), the CDA is shown to perform better in the forecasted state variables and simulated precipitation. The CDA is also shown to greatly mitigate the loss of forecast skills in observation denial experiments when radar radial winds are withheld in the assimilation. Modifications to the algorithm and sensitivities of the CDA to the calculation of the time tendencies in the constraints are described.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Dispersion Validation for Flow Involving a Large Structure Revisited: 45 Degree Rotation

The atmospheric dispersion of contaminants in the wake of a large urban structure is a challenging fluid mechanics problem of interest to the scientific and engineering communities. Magnetic Resonance Velocimetry (MRV) and Magnetic Resonance Concentration (MRC) are relatively new techniques that leverage diagnostic equipment used primarily by the medical field to make 3D engineering measurements of flow and contaminant dispersal. SIERRA/Fuego, a computational fluid dynamics (CFD) code at Sandia National Labs is employed to make detailed comparisons to the dataset to evaluate the quantitative and qualitative accuracy of the model. This work is the second in a series of scenarios. In the prior work, a single large building in an array of similar buildings was considered with the wind perpendicular to a building face. In this work, the geometry is rotated by 45 degrees and improved studies are performed for simulation credibility. The comparison exercise shows conditionally good comparisons between the model and experiment. Model uncertainties are assessed through parametric variations. Various methods of quantifying the accuracy between experiments and data are examined Three-dimensional analysis of accuracy is performed. The effort helped identify deficiencies in the techniques used to make these comparisons, and further methods development therefore becomes one of the main recommendations for follow-on work.

42 ENGINEERING↗