Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hamiltonian principles”

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

Nature of Hyperfine Interactions in TbPc 2 Single-Molecule Magnets: Multiconfigurational Ab Initio Study

Lanthanide-based single-ion magnetic molecules can have large magnetic hyperfine interactions as well as large magnetic anisotropy. Recent experimental studies reported tunability of these properties by changes of chemical environments or by application of external stimuli for device applications. In order to provide insight onto the origin and mechanism of such tunability, here we investigate the magnetic hyperfine and nuclear quadrupole interactions for a 159 Tb nucleus in TbPc 2 (Pc = phthalocyanine) single-molecule magnets using multiconfigurational ab initio methods including spin–orbit interaction. Since the electronic ground and first-excited (quasi)doublets are well separated in energy, the microscopic Hamiltonian can be mapped onto an effective Hamiltonian with an electronic pseudospin S = 1/2. From the ab initio calculated parameters, we find that the magnetic hyperfine coupling is dominated by the interaction of the Tb nuclear spin with electronic orbital angular momentum. The asymmetric 4f-like electronic charge distribution leads to a strong nuclear quadrupole interaction with significant transverse terms for the molecule with low symmetry. The ab initio calculated electronic–nuclear spectrum including the magnetic hyperfine and quadrupole interactions is in excellent agreement with the experiment. We further find that the transverse quadrupole interactions significantly influence the avoided level crossings in magnetization dynamics and that the molecular distortions affect mostly the Fermi contact terms as well as the transverse quadrupole interactions. Synopsis Hyperfine interactions for 159 Tb nucleus in TbPc 2 single-molecule magnets are investigated from first-principles using multiconfigurational calculations. Strong nuclear quadrupole and magnetic hyperfine coupling are found, with the latter being dominated by the paramagnetic spin-orbital mechanism. Furthermore, we construct an ab initio pseudospin Hamiltonian and obtain the electronic-nuclear spectrum that is in excellent agreement with experiment. The Zeeman diagram is calculated, and the magnetization dynamics is discussed. The effects of molecular distortions are researched.

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

An attempt to simulate laser-induced all-optical spin switching in a crystalline ferrimagnet

Interest in all-optical spin switching (AOS) is growing rapidly. The recent discovery of AOS in Mn 2 RuGa provides a much needed clean case of crystalline ferrimagnets for theoretical simulations. Here, we attempt to simulate it using the state-of-the-art first-principles method combined with the Heisenberg exchange model. We first compute the spin moments at two inequivalent manganese sites and then feed them into our model Hamiltonian. We employ an ultrafast laser pulse to switch the spins. We find that there is a similar optimal laser field amplitude to switch spins. However, we find that the exchange interaction has a significant effect on the system switchability. Weakening the exchange interaction could make the system unswitchable. This provides a crucial insight into the switching mechanism in ferrimagnets.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Passivity-Based Control of Grid Forming and Grid Following Converters in Microgrids

With the integration of more and more power electronics devices into the grid, it is not easy to guarantee the stability of the whole system due to the complexity of the control structure of various power converters. In this paper, we consider grid-forming and grid-following converters in one microgrid, in which the grid-forming converters support the voltage and frequency of the microgrid and the grid-following converters inject their maximum power into the microgrid. To handle the stability issue, the passivity principle is applied to guarantee the stability of the whole microgrid. We ensure every grid-following converter satisfies the passivity via the port-controlled Hamiltonian method. In addition, we apply the passivity-based proportional-resonant controller to the grid-forming converter. One of the advantages of the proposed method is that the phase lock loop system is eliminated, which may cause stability issues. Simulation results show that the proposed method can control the microgrid effectively.

Gui, Yonghao (Young)↗
DOE OSTI2026

$\mathscr{PT}$ symmetry enforced twin exchange as the origin of chirality-induced spin selectivity

Chiral molecules, ubiquitous in chemistry and biology, can differentiate electrons by their spin, a phenomenon known as chirality-induced spin selectivity (CISS). Despite its robustness and technological relevance, CISS has resisted conventional explanation: Spin-orbit coupling (SOC) models cannot fully account for the observed magnitude, room-temperature persistence, or equilibrium signatures. Here, we argue that structural chirality enforces a twin-pair exchange mechanism via the indistinguishability principle, which intrinsically couples spin and spatial degrees of freedom such that wave functions cannot be factorized into spin and spatial components. We derive an effective Hamiltonian that describes both transport and equilibrium CISS phenomena and is non-Hermitian. However, the inherent pseudo-Hermiticity, with $\mathscr{PT}$ symmetry as a special case, ensures real eigenvalues and thermodynamic consistency. We demonstrate that our framework is a step toward resolving long-standing anomalies of CISS. It situates CISS alongside equilibrium symmetry-breaking phenomena such as ferromagnetism and superconductivity, with implications for spintronics, catalysis, and the origins of biological homochirality.

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

Solving gyrokinetic systems with higher-order time dependence

We discuss theoretical and numerical aspects of gyrokinetics as a Lagrangian field theory when the field perturbation is introduced into the symplectic part. A consequence is that the field equations and particle equations of motion in general depend on the time derivatives of the field. The most well-known example is when the parallel vector potential is introduced as a perturbation, where a time derivative of the field arises only in the equations of motion, so an explicit equation for the fields may still be written. We will consider the conceptually more problematic case where the time-dependent fields appear in both the field equations and equations of motion, but where the additional term in the field equations is formally small. The conceptual issues were described by Burby (J. Plasma Phys., vol. 82 (3), 2016, 905820304): these terms lead to apparent additional degrees of freedom to the problem, so that the electric field now requires an initial condition, which is not required in low-frequency (Darwin) Vlasov–Maxwell equations. Also, the small terms in the Euler–Lagrange equations are a singular perturbation, and these two issues are interlinked. For well-behaved problems the apparent additional degrees of freedom are spurious, and the physically relevant solution may be directly identified. Because we needed to assume that the system is well behaved for small perturbations when deriving gyrokinetic theory, we must continue to assume that when solving it, and the physical solutions are thus the regular ones. The spurious nature of the singular degrees of freedom may also be seen by changing coordinate systems so the varying field appears only in the Hamiltonian. We then describe how methods appropriate for singular perturbation theory may be used to solve these asymptotic equations numerically. We then describe a proof-of-principle implementation of these methods for an electrostatic strong-flow gyrokinetic system; two basic test cases are presented to illustrate code functionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

On a variational formulation of the weakly nonlinear magnetic Rayleigh–Taylor instability

The magnetic-Rayleigh–Taylor (MRT) instability is a ubiquitous phenomenon that occurs in magnetically-driven Z-pinch implosions. It is important to understand this instability since it can decrease the performance of such implosions. In this work, I present a theoretical model for the weakly nonlinear MRT instability. I obtain such a model by asymptotically expanding an action principle, whose Lagrangian leads to the fully nonlinear MRT equations. After introducing a suitable choice of coordinates, I show that the theory can be cast as a Hamiltonian system, whose Hamiltonian is calculated up to the sixth order in a perturbation parameter. The resulting theory captures the harmonic generation of MRT modes. It is shown that the amplitude at which the linear magnetic-Rayleigh–Taylor instability exponential growth saturates depends on the stabilization effect of the magnetic-field tension. Overall, the theory provides an intuitive interpretation of the weakly nonlinear MRT instability and provides a systematic approach for studying this instability in more complex settings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Leakage Reduces Device Coherence Demands for Pulse-Level Molecular Simulations

Quantum simulation on noisy intermediate-scale quantum devices is severely limited by short qubit coherence times. A variational pulse-shaping algorithm known as ctrl-VQE was recently proposed to address this issue by eliminating the need for parameterized quantum circuits, which lead to long state preparation times. Here, in this study, we find the fastest possible pulses that prepare target molecular wave functions for a given device Hamiltonian describing coupled transmon qubits using simulations. We find that the resulting time-optimal pulses develop a bang-bang form consistent with Pontryagin’s maximum principle. We further investigate how the minimal state preparation time is impacted by the number of energy levels active in the transmon simulations. We find that leakage outside the computational subspace (something that is usually considered problematic) speeds up the state preparation, further reducing device coherence-time demands. Our analysis reveals that this speedup is due to both an enlarged solution space of target wave functions and the appearance of additional channels connecting initial and target states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Temperature dependence of magnetic anisotropy and magnetoelasticity from classical spin-lattice calculations

Here we present a classical molecular-spin dynamics (MSD) methodology that enables accurate computations of the temperature dependence of the magnetocrystalline anisotropy as well as magnetoelastic properties of magnetic materials. The nonmagnetic interactions are accounted for by a spectral neighbor analysis potential (SNAP) machine-learned interatomic potential, whereas the magnetoelastic contributions are accounted for using a combination of an extended Heisenberg Hamiltonian and a Néel pair interaction model, representing both the exchange interaction and spin-orbit-coupling effects, respectively. All magnetoelastic potential components are parameterized using a combination of first-principles and experimental data. Our framework is applied to the α phase of iron. Initial testing of our MSD model is done using a 0 K parametrization of the Néel interaction model. After this, we examine how individual Néel parameters impact the $B$ 1 and $B$ 2 magnetostrictive coefficients using a moment-independent δ sensitivity analysis. The results from this study are then used to initialize a genetic algorithm optimization which explores the Néel parameter phase space and tries to minimize the error in the B 1 and B 2 magnetostrictive coefficients in the range of 0–1200 K. Our results show that while both the 0 K and genetic algorithm optimized parametrization provide good experimental agreement for $B$ 1 and $B$ 2 , only the genetic algorithm optimized results can capture the second peak in the $B$ 1 magnetostrictive coefficient which occurs near approximately 800 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Skyrmion formation in Ni-based Janus dihalide monolayers: Interplay between magnetic frustration and Dzyaloshinskii-Moriya interaction

Here we present a comprehensive theory of the magnetic phases in monolayer nickel-based Janus dihalides (NiIBr, NiICl, NiBrCl) through a combination of first-principles calculations and atomistic simulations. Phonon band structure calculations and finite-temperature molecular dynamics simulations show that Ni-dihalide Janus monolayers are stable. The parameters of the interacting spin Hamiltonian extracted from ab initio calculations show that nickel-dihalide Janus monolayers exhibit varying degrees of magnetic frustration together with a large Dzyaloshinskii-Moriya interaction due to their inherent inversion symmetry breaking. The atomistic simulations reveal a delicate interplay between these two competing magnetic interactions in giving rise to skyrmion formation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Promise of Graph Sparsification and Decomposition for Noise Reduction in QAOA: Analysis for Trapped-Ion Compilations

We develop new approximate compilation schemes that significantly reduce the expense of compiling the Quantum Approximate Optimization Algorithm (QAOA) for solving the Max-Cut problem. Our main focus is on compilation with trapped-ion simulators using Pauli-X operations and all-to-all Ising Hamiltonian HIsing evolution generated by Molmer-Sorensen or optical dipole force interactions, though some of our results also apply to standard gate-based compilations. Our results are based on principles of graph sparsification and decomposition; the former reduces the number of edges in a graph while maintaining its cut structure, while the latter breaks a weighted graph into a small number of unweighted graphs. Though these techniques have been used as heuristics in various hybrid quantum algorithms, there have been no guarantees on their performance, to the best of our knowledge. This work provides the first provable guarantees using sparsification and decomposition to improve quantum noise resilience and reduce quantum circuit complexity. For quantum hardware that uses edge-by-edge QAOA compilations, sparsification leads to a direct reduction in circuit complexity. For trapped-ion quantum simulators implementing all-to-all HIsing pulses, we show that for a (1−ϵ) factor loss in the Max-Cut approximation (ϵ>0), our compilations improve the (worst-case) number of HIsing pulses from O(n2) to O(nlog(n/ϵ)) and the (worst-case) number of Pauli-X bit flips from O(n2) to O(nlog(n/ϵ)ϵ2) for n-node graphs. This is an asymptotic improvement for any constant ϵ>0. We demonstrate that significant improvements to the approximation ratio are obtained using decomposition in simulated trapped-ion experiments with dephasing noise. We further present a generic argument showing that sparsification results in an exponentially improved circuit fidelity lower bound in digital computing schemes based on one- and two-qubit gates, which are relevant to a wide variety of hardwares such as superconducting qubits and certain neutral atom or trapped ion setups, and more sophisticated noise models. We anticipate these approximate compilation techniques will be useful tools in a variety of future quantum computing experiments.

Moondra, Jai [Georgia Institute of Technology]↗
DOE OSTI2025

DFT-based insight into finite-temperature properties of ferroelectric perovskites with lone-pair: the case of CsGeX 3 (X = Cl, Br, I)

Ferroelectrics remain in the focus of scientific attention for decades owing to their fundamental and practical appeal. Recently, ferroelectricity has been demonstrated in semiconducting halide perovskites (Zhang et al 2022 Sci. Adv. 8 eabj5881), offering both a rare combination of ferroelectricity and semiconductivity in the same material and a possible alternative to the prevailing perovskite oxide ferroelectrics. We propose a route to simulating such materials at finite temperatures capable of reproducing key experimental and first-principle data, such as Curie temperature, phase transition sequence, spontaneous polarization, and soft mode frequencies. The key methodological finding is the superior performance of hybrid exchange correlation functionals in parametrization of effective Hamiltonians for ferroelectrics with lone pair. The parametrization for effective Hamiltonians for CsGeX 3 (X = Cl, Br, I) is reported. The application of methodology to study polarization reversal in CsGeX 3 allows for the development of a ‘minimalistic’ model for polarization reversal in ferroelectrics that provides an insight into the mechanisms of polarization reversal and its key features, such as the relationship between the coercive field, temperature, and AC field frequency. Importantly, the model reveals the origin of the well-known and ever-puzzling overestimation of coercive fields in computations. Furthermore, we report a variety of finite-temperature properties of CsGeX 3 ferroelectrics, such as dielectric susceptibility, pyroelectric coefficients, and energy storage density, which reveal that these halide perovskites possess properties comparable to their oxide counterparts. Here, we believe that our work provides significant methodological advancements, deepens fundamental understanding of ferroelectrics, and reveals the potential of halide perovskite ferroelectrics.

effective Hamiltonian↗
DOE OSTI2025

Systematic many-fermion Hamiltonian input scheme and spectral calculations on quantum computers

We present a novel input scheme for general second-quantized Hamiltonians of relativistic or non-relativistic many-fermion systems. This input scheme incorporates the fermionic anticommutation relations, particle number variations, and respects the symmetries of the Hamiltonian. Based on our input scheme, we propose a hybrid quantum-classical framework for spectral calculations on future quantum hardwares. We provide explicit circuit designs and the associated gate cost. We demonstrate our hybrid framework by solving the low-lying spectra of 42 Ca and 46 Ca. Our input scheme provides new pathways to solving the spectra and dynamics of the relativistic and nonrelativistic many-fermion systems via first-principles approaches.

Hybrid spectral calculation framework↗
DOE OSTI2020

27 AlNMR chemical shift of Al(OH)$^{-}_{4}$ calculated from first principles: Assessment of error cancellation in chemically distinct reference and target systems

Predicting accurate nuclear magnetic resonance chemical shieldings relies upon cancellation of different types of errors between the theoretically calculated shielding constant of the analyte of interest and the reference. Often, the intrinsic error in computed shieldings due to basis sets, approximations in the Hamiltonian, description of the wave function, and dynamic effects is nearly identical between the analyte and reference, yet if the electronic structure or sensitivity to local environment differs dramatically, this cannot be taken for granted. Detailed prior work has examined the octahedral trivalent cation Al(H 2 O)$^{3+}_{6}$, accounting for ab initio intrinsic errors. However, the use of this species as a reference for the chemically distinct tetrahedral anion Al(OH)$^{-}_{4}$ requires an understanding of how these errors cancel, in order to define the limits of accurately predicting 27 Al chemical shielding in Al(OH)$^{-}_{4}$. In this work, we estimate the absolute shielding of the 27 Al nucleus in Al(OH)$^{-}_{4}$ at the coupled cluster level (515.1± 5.3 ppm). Shielding sensitivity to the choice of method approximation and atomic basis sets used has been evaluated. Solvent and thermal effects are assessed through ensemble averaging techniques using ab-initio molecular dynamics. The contribution of each type of intrinsic error is assessed for the Al(H 2 O)$^{3+}_{6}$ and Al(OH)$^{-}_{4}$ ions, revealing significant differences that fundamentally hamper the ability to accurately calculate the 27 Al chemical shift of Al(OH)$^{-}_{4}$ from first principles.

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

Benchmarking machine learning interatomic potentials via phonon anharmonicity

Abstract Machine learning approaches have recently emerged as powerful tools to probe structure-property relationships in crystals and molecules. Specifically, machine learning interatomic potentials (MLIPs) can accurately reproduce first-principles data at a cost similar to that of conventional interatomic potential approaches. While MLIPs have been extensively tested across various classes of materials and molecules, a clear characterization of the anharmonic terms encoded in the MLIPs is lacking. Here, we benchmark popular MLIPs using the anharmonic vibrational Hamiltonian of ThO 2 in the fluorite crystal structure, which was constructed from density functional theory (DFT) using our highly accurate and efficient irreducible derivative methods. The anharmonic Hamiltonian was used to generate molecular dynamics (MD) trajectories, which were used to train three classes of MLIPs: Gaussian approximation potentials, artificial neural networks (ANN), and graph neural networks (GNN). The results were assessed by directly comparing phonons and their interactions, as well as phonon linewidths, phonon lineshifts, and thermal conductivity. The models were also trained on a DFT MD dataset, demonstrating good agreement up to fifth-order for the ANN and GNN. Our analysis demonstrates that MLIPs have great potential for accurately characterizing anharmonicity in materials systems at a fraction of the cost of conventional first principles-based approaches.

interatomic potentials↗
DOE OSTI2022

Analysis of the classical trajectory treatment of photon dynamics for polaritonic phenomena

Simulating photon dynamics in strong light–matter coupling situations via classical trajectories is proving to be powerful and practical. Here, in this work, we analyze the performance of the approach through the lens of the exact factorization approach. Since the exact factorization enables a rigorous definition of the potentials driving the photonic motion, it allows us to identify that the underestimation of photon number and intensities observed in earlier work is primarily due to an inadequate accounting of light–matter correlation in the classical Ehrenfest force rather than errors from treating the photons quasiclassically per se. The latter becomes problematic when the number of photons per mode begins to exceed a half.

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

Degradation of performance in ICF implosions due to Rayleigh–Taylor instabilities: A Hamiltonian perspective

The Rayleigh–Taylor instability (RTI) is an ubiquitous phenomenon that occurs in inertial-confinement-fusion (ICF) implosions and is recognized as an important limiting factor of ICF performance. To analytically understand the RTI dynamics and its impact on ICF capsule implosions, we develop a first-principle variational theory that describes an imploding spherical shell undergoing RTI. The model is based on a thin-shell approximation and includes the dynamical coupling between the imploding spherical shell and an adiabatically compressed fluid within its interior. Using a quasilinear analysis, we study the degradation trends of key ICF performance metrics (e.g., stagnation pressure, residual kinetic energy, and areal density) as functions of initial RTI parameters (e.g., the initial amplitude and Legendre mode), as well as the 1D implosion characteristics (e.g., the convergence ratio). We compare analytical results from the theory against nonlinear results obtained by numerically integrating the governing equations of this reduced model. Our findings emphasize the need to incorporate polar flows in the calculation of residual kinetic energy and demonstrate that higher convergence ratios in ICF implosions lead to significantly greater degradation of key performance metrics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Feedback-based quantum algorithms for ground state preparation

The ground state properties of quantum many-body systems are a subject of interest across chemistry, materials science, and physics. Thus, algorithms for finding ground states can have broad impacts. Variational quantum algorithms are one class of ground state algorithms that has received significant attention in recent years. These algorithms utilize a hybrid quantum-classical computing framework to prepare ground states on quantum computers. However, this requires solving a classical optimization problem that can become prohibitively expensive in high dimensions. Here, we develop formulations of feedback-based quantum algorithms for ground state preparation that can be used to address this challenge for two broad classes of Hamiltonians: Fermi-Hubbard Hamiltonians, and molecular Hamiltonians represented in second quantization. Feedback-based quantum algorithms are optimization-free; in place of classical optimization, quantum circuit parameters are set according to a deterministic feedback law derived from quantum Lyapunov control principles. This feedback law guarantees a monotonic improvement in solution quality with respect to the depth of the quantum circuit. A variety of numerical illustrations are provided that analyze the convergence and robustness of feedback-based quantum algorithms for these problem classes. Published by the American Physical Society 2024

Larsen, James B. (ORCID:000000020777440X)↗
DOE OSTI2025

Mechanical cosmology: Simulating scalar fluctuations in expanding universes using synthetic mechanical lattices

Inspired by recent advances in observational astrophysics and continued explorations in the field of analog gravity, we discuss the prospect of simulating models of cosmology within the context of synthetic mechanical lattice experiments. We focus on the physics of expanding universe scenarios described by the Friedmann-Lemaître-Robertson-Walker (FLRW) metric. Specifically, quantizing scalar fluctuations in a background FLRW spacetime leads to a quadratic bosonic Hamiltonian with temporally varying pair production terms. Here we present a mapping that provides a one-to-one correspondence between these classes of cosmology models and feedback-coupled mechanical oscillators. As proof of principle, we then perform experiments on a synthetic mechanical lattice composed of such oscillators. We simulate two different FLRW expansion scenarios with universes dominated by vacuum energy and matter and discuss our experimental results. Published by the American Physical Society 2025

Rhyno, Brendan (ORCID:0000000235266998)↗