SEARCH · Engineering Papers

Results for “HAMILTONIAN”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1978

Numerical integration of nearly-Hamiltonian systems

The reported investigation is concerned with the solution of systems of differential equations which are derived from a Hamiltonian function in the extended phase space. The problem selected involves a one-dimensional perturbed harmonic oscillator. The van der Pol equation considered has an exact asymptotic value for its amplitude. Comparisons are made between a numerical solution and a known analytical solution. In addition to the van der Pol problem, known solutions regarding the restricted problem of three bodies are used as examples for perturbed Keplerian motion. The extended phase space Hamiltonian discussed by Stiefel and Scheifele (1971) is considered. A description is presented of two canonical formulations of the perturbed harmonic oscillator.

Bond, V. R.↗

NASA NTRS1988

A study to evaluate STS heads-up ascent trajectory performance employing a minimum-Hamiltonian optimization strategy

A study was conducted to evaluate the performance implications of a heads-up ascent flight design for the Space Transportation System, as compared to the current heads-down flight mode. The procedure involved the use of the Minimum Hamiltonian Ascent Shuttle Trajectory Evaluation Program, which is a three-degree-of-freedom moment balance simulation of shuttle ascent. A minimum-Hamiltonian optimization strategy was employed to maximize injection weight as a function of maximum dynamic pressure constraint and Solid Rocket Motor burnrate. Performance Reference Mission Four trajectory groundrules were used for consistency. The major conclusions are that for heads-up ascent and a mission nominal design maximum dynamic pressure value of 680 psf, the optimum solid motor burnrate is 0.394 ips, which produces a performance enhancement of 4293 lbm relative to the baseline heads-down ascent, with 0.368 ips burnrate solid motors and a 680 psf dynamic pressure constraint. However, no performance advantage exists for heads-up flight if the current Solid Rocket Motor target burnrate of 0.368 ips is used. The advantage of heads-up ascent flight employing the current burnrate is that Space Shuttle Main Engine throttling for dynamic pressure control is not necessary.

Sinha, Sujit↗

NASA NTRS1990

Boson Hamiltonians and stochasticity for the vorticity equation

The evolution of the vorticity in time for two-dimensional inviscid flow and in Lagrangian time for three-dimensional viscous flow is written in Hamiltonian form by introducing Bose operators. The addition of the viscous and convective terms, respectively, leads to an interpretation of the Hamiltonian contribution to the evolution as Langevin noise.

Shen, Hubert H.↗

NASA NTRS1992

The Hamiltonian structure of Dirac's equation in tensor form and its Fermi quantization

Currently, there is some interest in studying the tensor forms of the Dirac equation to elucidate the possibility of the constrained tensor fields admitting Fermi quantization. We demonstrate that the bispinor and tensor Hamiltonian systems have equivalent Fermi quantizations. Although the tensor Hamiltonian system is noncanonical, representing the tensor Poisson brackets as commutators for the Heisenberg operators directly leads to Fermi quantization without the use of bispinors.

Reifler, Frank↗

NASA NTRS1992

Non-Hamiltonian perturbations of integrable systems and resonance trapping

This paper studies general, non-Hamiltonian perturbations of integrable systems with two degrees of freedom and derives conditions for temporary and permanent resonance trapping. The analysis involves a noncanonical transformation of variables near the resonant manifold and averaging with respect to the fast phase to investigate oscillatory behavior on the intermediate timescale. The resulting reduced system is Hamiltonian to leading order and permits, after averaging on the intermediate, or libration, timescale, a canonical transformation to action-angle variables in the oscillation zone. The final system so obtained reveals the possible existence of two- and three-dimensional invariant tori in the vicinity of the resonant manifold. An explicit divergence condition for general perturbations to be dissipative on the slow timescale follows from the analysis. An application of this approach to the problem of resonant trapping and escape is outlined for the restricted problem of three bodies subject to dissipative perturbations with a radial symmetry.

Ghil, M.↗

NASA NTRS1994

Sqeezing generated by a nonlinear master equation and by amplifying-dissipative Hamiltonians

In the first part of this contribution we show that the master equation derived from the generalized version of the nonlinear Doebner-Goldin equation leads to the squeezing of one of the quadratures. In the second part we consider two familiar Hamiltonians, the Bateman- Caldirola-Kanai and the optical parametric oscillator; going back to their classical Lagrangian form we introduce a stochastic force and a dissipative factor. From this new Lagrangian we obtain a modified Hamiltonian that treats adequately the simultaneous amplification and dissipation phenomena, presenting squeezing, too.

Dodonov, V. V.↗

NASA NTRS1994

On The Choice of a Zeroth-Order Hamiltonian for Second-Order Perturbation Theory with A CASSCF Reference Function

A new approach to perturbation theory based on a CASSCF reference function has been developed. The key to the approach is the definition of the zeroth order Hamiltonian, H(sub 0), which includes the full CI Hamiltonian for the active space. In the the inactive and secondary spaces, operators may be chosen which reduce to the usual Moller-Plesset or Epstein-Nesbet forms in the limit of a null active space. These operators are diagonal in the orbital indices and permit the block-diagonalization of H(sub 0). The reference is an eigenfunction of Ho without N-particle projection. H(sub 0) automatically incorporates denominator shifts in the style of those appearing in recent open-shell perturbation theories. Comparative results are presented for a few test cases.

Dyall, Kenneth G.↗

NASA NTRS2017

Exponentially-Biased Ground-State Sampling of Quantum Annealing Machines with Transverse-Field Driving Hamiltonians

We study the performance of the D-Wave 2X quantum annealing machine on systems with well-controlled ground-state degeneracy. While obtaining the ground state of a spin-glass benchmark instance represents a difficult task, the gold standard for any optimization algorithm or machine is to sample all solutions that minimize the Hamiltonian with more or less equal probability. Our results show that while naive transverse-field quantum annealing on the D-Wave 2X device can find the ground-state energy of the problems, it is not well suited in identifying all degenerate ground-state configurations associated to a particular instance. Even worse, some states are exponentially suppressed, in agreement with previous studies on toy model problems [New J. Phys. 11, 073021 (2009)]. These results suggest that more complex driving Hamiltonians are needed in future quantum annealing machines to ensure a fair sampling of the ground-state manifold.

Quantum Annealing↗

NASA NTRSDate not supplied

Improved Fermion Hamiltonians for Quantum Simulations

Constructing improved hamiltonians for gauge theories coupled to fermonic matter will be important for improving continuum limit extrapolations of quantum computations. In this talk we will present a formulation for simulating ASQTAD fermions for lattice computation and provide fault tolerant resource costs in terms of primitive group operations. We additionally show that the scaling of energies with respect to the lattice spacing are better than for the unimproved Hamiltonian for toy models.

Erik Joseph Gustafson↗

NASA NTRSDate not supplied

Improved Fermion Hamiltonians for Quantum Simulations

Erik Gustafson↗

NASA NTRS2023

Improved Fermion Hamiltonians for Quantum Simulation

The Symanzik improvement program has been quite successful in classical simulations of quantum chromodynamics allowing calculations to be performed at coarser lattice spacings and with reduced computational resource costs. It is expected that improved Hamiltonians will be essential to simulate lattice field theories using quantum computers. In this work I will discuss the formulation of an ASQTAD and HISQ Hamiltonian amenable for quantum simulations. I will also show preliminary results that demonstrate significant tree-level contributions are removed in the spectrum of the 1 flavor Schwinger model.

quantum computing↗

NASA NTRS2023

Improved Fermion Hamiltonians for Quantum Simulations

Quantum Algorithms↗

DOE OSTI2026

Gradient-informed Hamiltonian Monte Carlo for multicomponent CALPHAD model optimization and uncertainty quantification

CALPHAD model parameter optimization is inherently challenging due to non-smooth objective functions, high-dimensional parameter spaces, and the need for uncertainty quantification (UQ). Traditional weighted nonlinear least squares approaches are computationally efficient but local, whereas black-box global optimizers and ensemble Markov Chain Monte Carlo (MCMC) methods provide broader exploration at substantial computational cost. The objective of this work is to combine the global exploration capability of gradient-informed Hamiltonian Monte Carlo – specifically the No-U-Turn Sampler (NUTS) – with local deterministic refinement using BFGS to efficiently optimize multicomponent CALPHAD models with minimal manual intervention. Analytic gradients are computed via the Jansson derivative framework. The methodology is demonstrated on the Cr—Fe binary system and extended to the Cr—Fe—Ni ternary system with 32 degrees of freedom. For Cr—Fe, NUTS achieves comparable or superior optimality relative to ensemble MCMC while requiring over an order-of-magnitude fewer likelihood evaluations. Parameter uncertainties are quantified through NUTS sampling and propagated to thermodynamic observables using local expansion, demonstrating a novel modular approach that combines binary and ternary parameter subsets without requiring global relaxation. These results establish gradient-informed exploration as a scalable strategy for multicomponent CALPHAD optimization and provide a practical route towards efficient higher-order database development with quantified uncertainty.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Accelerating Hamiltonian Monte Carlo for Bayesian inference in neural networks and neural operators

Hamiltonian Monte Carlo (HMC) is a powerful and accurate method to sample from the posterior distribution in Bayesian inference. However, HMC techniques are computationally demanding for Bayesian neural networks due to the high dimensionality of the network’s parameter space and the non-convexity of their posterior distributions. Therefore, various approximation techniques, such as variational inference (VI) or stochastic gradient MCMC, are often employed to infer the posterior distribution of the network parameters. Such approximations introduce inaccuracies in the inferred distributions, resulting in unreliable uncertainty estimates. In this work, we propose a hybrid approach that combines inexpensive VI and accurate HMC methods to efficiently and accurately quantify uncertainties in neural networks and neural operators. The proposed approach leverages an initial VI training on the full network. We examine the influence of individual parameters on the prediction uncertainty, which shows that a large proportion of the parameters do not contribute substantially to uncertainty in the network predictions. This information is then used to significantly reduce the dimension of the parameter space, and HMC is performed only for the subset of network parameters that strongly influence prediction uncertainties. This yields a framework for accelerating the full batch HMC for posterior inference in neural networks. We demonstrate the efficiency and accuracy of the proposed framework on deep neural networks and operator networks, showing that inference can be performed for large networks with tens to hundreds of thousands of parameters. Finally, we show that this method can effectively learn surrogates for complex physical systems by modeling the operator that maps from upstream conditions to wall-pressure data on a cone in hypersonic flow.

Bayesian inference↗

DOE OSTI2025

CLPNets: Coupled Lie–Poisson neural networks for multi-part Hamiltonian systems with symmetries

To accurately compute data-based prediction of Hamiltonian systems, it is essential to utilize methods that preserve the structure of the equations over time. We consider a particularly challenging case of systems with interacting parts that do not reduce to pure momentum evolution. Such systems are essential in scientific computations, such as discretization of a continuum elastic rod, which can be viewed as the group of rotations and translations $SE(3)$. The evolution involves not only the momenta but also the relative positions and orientations of the particles. The presence of Lie group-valued elements, such as relative positions and orientations, poses a problem for applying previously derived methods for data-based computing. We develop a novel method of data-based computation and complete phase space learning of such systems. We follow the original framework of SympNets (Jin et al., 2020) and LPNets (Eldred et al., 2024), building the neural network from phase space mappings that preserve the Lie–Poisson structure. We derive a novel system of mappings that are built into neural networks describing the evolution of such systems. We call such networks Coupled Lie–Poisson Neural Networks, or CLPNets. We consider increasingly complex examples for the applications of CLPNets, starting with the rotation of two rigid bodies about a common axis, progressing to the free rotation of two rigid bodies, and finally to the evolution of two connected and interacting $SE(3)$ components, describing the discretization of an elastic rod into two elements. Our method preserves all Casimir invariants to machine precision, preserves energy to high accuracy, and shows good resistance to the curse of dimensionality, requiring only a few thousand data points for all cases studied (three to eighteen dimensions). Additionally, the method is highly economical in memory requirements, requiring only about 200 parameters for the most complex case considered.

Data-based modeling↗

DOE OSTI2025

Scalable Implementation of Mean-Field and Correlation Methods Based on Lie-Algebraic Similarity Transformation of Spin Hamiltonians in the Jordan–Wigner Representation

Recent work has highlighted that the strong correlation inherent in spin Hamiltonians can be effectively reduced by mapping spins to Fermions via the Jordan−Wigner transformation (JW). The Hartree−Fock method is straightforward in the Fermionic domain and may provide a reasonable approximation to the ground state. Correlation with respect to the Fermionic mean field can be recovered based on Lie-algebraic similarity transformation (LAST) with two-body correlators. Specifically, a unitary LAST variant eliminates the dependence on site ordering, while a nonunitary LAST yields size-extensive correlation energies. Whereas the first recent demonstration of such methods was restricted to small spin systems, we present efficient implementations using analytical gradients for the optimization with respect to the mean-field reference and the LAST parameters, thereby enabling the treatment of larger clusters, including systems with local spins s > $\frac{1}{2}$.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

SAP-X2C: Optimally-Simple Two-Component Relativistic Hamiltonian with Size-Intensive Picture Change

We present a simple relativistic exact 2-component (X2C) Hamiltonian that models two-electron picture-change effects using Lehtola’s superposition of atomic potentials (SAP) [S. Lehtola, J. Chem. Theory Comput. 15, 1593−1604 (2019)]. The SAP-X2C approach retains the low cost and technical simplicity of the popular 1-electron X2C (1eX2C) predecessor but is significantly more accurate and has a well-defined thermodynamic limit, making it applicable to extended systems (such as large molecules and periodic crystals). The assessment of the SAP-X2C-based Hartree−Fock total and spinor energies, spin−orbit splittings, equilibrium bond distances, and harmonic vibrational frequencies suggests that SAP-X2C is similar to the more complex atomic meanfield (AMF) X2C counterparts in its ability to approximate the 4-component Dirac−Hartree−Fock reference.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Laplace Transform–Based Quantum Eigenvalue Transformation via Linear Combination of Hamiltonian Simulation

Eigenvalue transformations, which include solving time-dependent differential equations as a special case, have a wide range of applications in scientific and engineering computation. While quantum algorithms for singular value transformations are well studied, eigenvalue transformations are distinct, especially for nonnormal matrices. Here, we propose an efficient quantum algorithm for performing a class of eigenvalue transformations that can be expressed as a certain type of matrix Laplace transformation. This allows us to significantly extend the recently developed linear combination of Hamiltonian simulation method [D. An, J.-P. Liu, and L. Lin, Phys. Rev. Lett., 131 (2023), 150603; D. An, A. M. Childs, and L. Lin, Commun. Math. Phys. 407, 19 (2026)] to represent a wider class of eigenvalue transformations, such as powers of the matrix inverse, 𝐴 −𝑘 , and the exponential of the matrix inverse, 𝑒 −𝐴 −1 . The latter can be interpreted as the solution of a mass-matrix differential equation of the form form 𝐴⁢𝑢′⁡⁡(𝑡) =−𝑢⁡(𝑡). We demonstrate that our eigenvalue transformation approach can solve this problem without explicitly inverting 𝐴, thereby reducing the computational complexity.

Laplace transform↗