Search NASA⌕ Search

SEARCH · Search NASA

Results for “Effective Hamiltonian models”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Fast inversion, preconditioned quantum linear system solvers, fast Green's-function computation, and fast evaluation of matrix functions

Preconditioning is the most widely used and effective way for treating ill-conditioned linear systems in the context of classical iterative linear system solvers. We introduce a quantum primitive called fast inversion, which can be used as a preconditioner for solving quantum linear systems. The key idea of fast inversion is to directly block encode a matrix inverse through a quantum circuit implementing the inversion of eigenvalues via classical arithmetics. We demonstrate the application of preconditioned linear system solvers for computing single-particle Green's functions of quantum many-body systems, which are widely used in quantum physics, chemistry, and materials science. We analyze the complexities in three scenarios: the Hubbard model, the quantum many-body Hamiltonian in the plane-wave-dual basis, and the Schwinger model. We also provide a method for performing Green's function calculation in second quantization within a fixed-particle manifold and note that this approach may be valuable for simulation more broadly. Aside from solving linear systems, fast inversion also allows us to develop fast algorithms for computing matrix functions, such as the efficient preparation of Gibbs states. Furthermore, we introduce two efficient approaches for such a task, based on the contour-integral formulation and the inverse transform, respectively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Real-time chiral dynamics at finite temperature from quantum simulation

In this study, we explore the real-time dynamics of the chiral magnetic effect (CME) at a finite temperature in the (1+1)-dimensional QED, the massive Schwinger model. By introducing a chiral chemical potential μ 5 through a quench process, we drive the system out of equilibrium and analyze the induced vector currents and their evolution over time. The Hamiltonian is modified to include the time-dependent chiral chemical potential, thus allowing the investigation of the CME within a quantum computing framework. We employ the quantum imaginary time evolution (QITE) algorithm to study the thermal states, and utilize the Suzuki-Trotter decomposition for the real-time evolution. This study provides insights into the quantum simulation capabilities for modeling the CME and offers a pathway for studying chiral dynamics in low-dimensional quantum field theories.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Gapless Floquet topology

Symmetry-protected topological (SPT) phases in insulators and superconductors are known for their robust edge modes, linked to bulk invariants through the bulk-boundary correspondence. While this principle traditionally applies to gapped phases, recent advances have extended it to gapless systems, where topological edge-states persist even in the absence of a bulk gap. We extend this framework to periodically driven chains with chiral symmetry, revealing the existence of topological edge zero and π modes despite the lack of bulk gaps in the quasienergy spectrum. By examining the half-period decomposition of chiral evolutions, we construct topological invariants that circumvent the need to define the Floquet Hamiltonian, making them more suitable for generalization to the gapless regime. We provide explicit examples, including generalizations of the Kitaev chain and related spin models, where localized π modes emerge even when the bulk is gapless at the same quasienergy as the edge modes. Here, we numerically study the effect of interactions, which give a finite lifetime to the edge modes in the thermodynamic limit with the decay rate consistent with Fermi's golden rule.

1-dimensional systems↗
DOE OSTI2025

Fermi’s Golden Rule Rate Expression for Transitions Due to Nonadiabatic Derivative Couplings in the Adiabatic Basis

Starting from a general molecular Hamiltonian expressed in the basis of adiabatic electronic and nuclear position states, where a compact and complete expression for the nonadiabatic derivative coupling (NDC) Hamiltonian term is obtained, we provide a general analysis of the Fermi’s golden rule (FGR) rate expression for nonadiabatic transitions between adiabatic states. We then consider a quasi-adiabatic approximation that uses crude adiabatic states and NDC couplings, both evaluated at the minimum potential energy configuration of the initial adiabatic state, for the definition of the zeroth and first-order terms of the Hamiltonian. Although the application of this approximation is rather limited, it allows deriving a general FGR rate expression without further approximation while accounting for non-Condon contribution to the FGR rate arising from momentum operators of NDC terms and its coupling with vibronic displacements. For a generic and widely used model where all nuclear degrees of freedom and environmental effects are represented as linearly coupled harmonic oscillators, we derive a closed-form FGR rate expression that requires only Fourier transform. The resulting rate expression includes quadratic contributions of NDC terms and their couplings to Franck–Condon modes, which require evaluation of two additional bath spectral densities in addition to the conventional one that appears in a typical FGR rate theory based on the Condon approximation. Model calculations for the case where nuclear vibrations consist of both a sharp high-frequency mode and an Ohmic bath spectral density illustrate new features and implications of the rate expression. We then apply our theoretical expression to the nonradiative decay from the first excited singlet state of azulene, which illustrates the utility and implications of our theoretical results.

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

Constant-depth circuits for dynamic simulations of materials on quantum computers

Abstract Dynamic simulation of materials is a promising application for near-term quantum computers. Current algorithms for Hamiltonian simulation, however, produce circuits that grow in depth with increasing simulation time, limiting feasible simulations to short-time dynamics. Here, we present a method for generating circuits that are constant in depth with increasing simulation time for a specific subset of one-dimensional (1D) materials Hamiltonians, thereby enabling simulations out to arbitrarily long times. Furthermore, by removing the effective limit on the number of feasibly simulatable time-steps, the constant-depth circuits enable Trotter error to be made negligibly small by allowing simulations to be broken into arbitrarily many time-steps. For an N -spin system, the constant-depth circuit contains only $\mathcal {O}(N^{2})$ O ( N 2 ) CNOT gates. Such compact circuits enable us to successfully execute long-time dynamic simulation of ubiquitous models, such as the transverse field Ising and XY models, on current quantum hardware for systems of up to 5 qubits without the need for complex error mitigation techniques. Aside from enabling long-time dynamic simulations with minimal Trotter error for a specific subset of 1D Hamiltonians, our constant-depth circuits can advance materials simulations on quantum computers more broadly in a number of indirect ways.

Materials simulation↗
NASA NTRS2012

Preliminary Modeling of CH3D From 4000 to 4550 cm (exp-1)

A new study of 12CH3D line positions and intensities was performed for the upper portion of the Enneadecad polyad between 4000 and 4550/cm. For this, FTIR spectra were recorded with D-enriched methane samples (at 80 K with a Bruker 125 IFS at 0.005/cm resolution and at 291 K with the McMath-Pierce FTS at 0.011/cm resolution, respectively). Line positions and intensities were retrieved by least square curve-fitting procedures and analyzed using the effective Hamiltonian and the effective dipole moment expressed in terms of irreducible tensor operators adapted to symmetric top molecules. Initially, only the cold spectrum was used to identify quantum assignments and predict 12CH3D relative intensities in this region. To assign higher quanta up to J equal 14, additional line positions and intensities were obtained from two room temperature spectra. In the final stage, measured intensities from both the cold and the room temperature data were normalized to corresponding values at 296 K and then averaged. Combining the two temperature datasets confirmed the assumed quantum assignments and also demonstrated the relative accuracies to be better than 70.0002/cm for line positions and at least +/-6% for intensities so that approx.1160 features were selected. Including additional assignments from the room temperature spectra alone permitted 1362 line intensities of 11 bands (involving 23 vibrational symmetry components) to be reproduced with an RMS of 9%. Over 4085 selected positions for 12 bands were modeled to 0.008/cm. Nevertheless a number of known assignments could not be modeled to within our experimental precisions. More work is needed to obtain a complete characterization of this complex polyad.

Nikitin, A. V.↗
DOE OSTI2020

Bayesian differential programming for robust systems identification under uncertainty

This paper presents a machine learning framework for Bayesian systems identification from noisy, sparse and irregular observations of nonlinear dynamical systems. The proposed method takes advantage of recent developments in differentiable programming to propagate gradient information through ordinary differential equation solvers and perform Bayesian inference with respect to unknown model parameters using Hamiltonian Monte Carlo sampling. This allows an efficient inference of the posterior distributions over plausible models with quantified uncertainty, while the use of sparsity-promoting priors enables the discovery of interpretable and parsimonious representations for the underlying latent dynamics. A series of numerical studies is presented to demonstrate the effectiveness of the proposed methods, including nonlinear oscillators, predator–prey systems and examples from systems biology. Taken together, our findings put forth a flexible and robust workflow for data-driven model discovery under uncertainty. All codes and data accompanying this article are available at https://bit.ly/34FOJMj .

Science & Technology - Other Topics↗
DOE OSTI2021

Quantum Monte Carlo studies of a trimer scaling function with microscopic two- and three-body interactions

Here, we present an energy scaling function to predict, in a specific range, the energy of bosonic trimers with large scattering lengths and finite range interactions, which is validated by quantum Monte Carlo calculations using microscopic Hamiltonians with two- and three-body potentials. The proposed scaling function depends on the scattering length, effective range, and a reference energy, which we chose as the trimer energy at unitarity. We obtained the scaling function as a limit cycle from the solution of the renormalized zero-range model with effective range corrections. We proposed a simple parametrization of the energy scaling function. Besides the intrinsic interest in theoretical and experimental investigations, this scaling function allows one to probe Efimov physics with only the trimer ground states, which may open opportunities to identify Efimov trimers whenever access to excited states is limited.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Spin-orbital correlations from complex orbital order in MgV 2⁢ O 4

MgV 2 ⁢O 4 is a spinel based on magnetic V 3+ ions, which host both spin (𝑆 = 1) and orbital (𝑙 eff = 1) moments. Owing to the underlying pyrochlore coordination of the magnetic sites, the spins in MgV 2 ⁢O 4 only antiferromagnetically order once the frustrating interactions imposed by the 𝐹⁡𝑑⁢$\overline{3}$⁢𝑚 lattice are broken through an orbitally-driven structural distortion at T 𝑆 ≃ ⁢60 K. Consequently, a Néel transition occurs at T 𝑁 ≃ ⁢40 K. Low-temperature spatial ordering of the electronic orbitals is fundamental to both the structural and magnetic properties; however, considerable discussion on whether it can be described by complex or real orbital ordering is ambiguous. We apply neutron spectroscopy to resolve the nature of the orbital ground state and characterize hysteretic spin-orbital correlations using x-ray and neutron diffraction. Neutron spectroscopy finds multiple excitation bands and we parametrize these in terms of a multilevel (or excitonic) theory based on the orbitally degenerate ground state. Meaningful for the orbital ground state, we report an “optical-like” mode at high energies that we attribute to a crystal-field-like excitation from the spin-orbital 𝑗 eff = 2 ground-state manifold to an excited 𝑗 eff =1 energy level. We parametrize the magnetic excitations in terms of a Hamiltonian with spin-orbit coupling and local crystalline electric field distortions resulting from deviations from perfect octahedra surrounding the V 3+ ions. We suggest that this provides compelling evidence for complex orbital order in MgV 2⁢ O 4 . We then apply the consequences of this model to understand hysteretic effects in the magnetic diffuse scattering where we propose that MgV 2⁢ O 4 displays a high-temperature orbital memory of the low-temperature spin order.

Frustrated magnetism↗
NASA NTRS1982

Forces acting on free carriers in semiconductors of inhomogeneous composition I, II

An investigation shows that the extended Wannier-Slater theorem with its position-dependent effective mass m cannot be considered a viable description of a semiconductor with compositional changes. The correct description of the behavior of the electron or hole flow in semiconductors in which the chemical composition changes gradually from position to position is more difficult. It is necessary to find a model Hamiltonian. This model Hamiltonian is to some extent arbitry, but it must be representative of a homogeneous sample. The difference between the true and the model Hamiltonian should be small in the sense that this difference may be considered a perturbation. However, a Hamiltonian with the required specifications can always be found provided that the compositional inhomogeneities do not change drastically over the distances of interest.

Von Roos, O.↗
DOE OSTI2025

Hamiltonian parameter inference from resonant inelastic x-ray scattering with active learning

Identifying model Hamiltonians is a vital step toward creating predictive models of materials. Here, in this study, we combine Bayesian optimization with the EDRIXS numerical package to infer Hamiltonian parameters from resonant inelastic x-ray scattering (RIXS) spectra within the single atom approximation. To evaluate the efficacy of our method, we test it on experimental RIXS spectra of NiPS 3 , NiCl 2 , Ca 3 ⁢LiOsO 6 , and Fe 2⁢ O 3 , and demonstrate that it can reproduce results obtained from hand-fitted parameters to a precision similar to expert human analysis while providing a more systematic mapping of parameter space. Our work provides a key first step toward solving the inverse scattering problem to extract effective multi-orbital models from information-dense RIXS measurements, which can be applied to a host of quantum materials. We also propose atomic model parameter sets for two materials, Ca 3⁢ LiOsO 6 and Fe 2⁢ O 3 , that were previously missing from the literature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Three-dimensional imaging of proton in basis light-front quantization

Here, we employ an effective Hamiltonian that includes the transverse and longitudinal confinement and the one-gluon exchange interaction with fixed coupling constant. By solving the eigenvalue equation in basis light-front quantization (BLFQ), we generate the light-front wavefunctions (LFWFs) for the nucleon in the valence quark Fock space. Fitting the model parameters, we obtain high quality descriptions of electromagnetic form factors and radius for proton while the results deviate somewhat from experimental data for neutron.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Illuminating Ligand Field Contributions to Molecular Qubit Spin Relaxation via T 1 Anisotropy

Electron spin relaxation in paramagnetic transition metal complexes constitutes a key limitation on the growth of molecular quantum information science. However, there exist very few experimental observables for probing spin relaxation mechanisms, leading to a proliferation of inconsistent theoretical models. Here we demonstrate that spin relaxation anisotropy in pulsed electron paramagnetic resonance is a powerful spectroscopic probe for molecular spin dynamics across a library of highly coherent Cu(II) and V(IV) complexes. Here, neither the static spin Hamiltonian anisotropy nor contemporary computational models of spin relaxation are able to account for the experimental T 1 anisotropy. Through analysis of the spin-orbit coupled wavefunctions, we derive an analytical theory for the T 1 anisotropy that accurately reproduces the average experimental anisotropy of 2.5. Furthermore, compound-by-compound deviations from the average anisotropy provide a promising approach for observing specific ligand field and vibronic excited state coupling effects on spin relaxation. Finally, we present a simple density functional theory workflow for computationally predicting T 1 anisotropy. Analysis of spin relaxation anisotropy leads to deeper fundamental understanding of spin-phonon coupling and relaxation mechanisms, promising to complement temperature-dependent relaxation rates as a key metric for understanding molecular spin qubits.

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

Modified Fermi’s golden rule rate expressions

Fermi’s golden rule (FGR) serves as the basis for many expressions of spectroscopic observables and quantum transition rates. The utility of FGR has been demonstrated through decades of experimental confirmation. However, there still remain important cases where the evaluation of a FGR rate is ambiguous or ill-defined. Examples are cases where the rate has divergent terms due to the sparsity in the density of final states or time dependent fluctuations of system Hamiltonians. Strictly speaking, assumptions of FGR are no longer valid for such cases. However, it is still possible to define modified FGR rate expressions that are useful as effective rates. The resulting modified FGR rate expressions resolve a long standing ambiguity often encountered in using FGR and offer more reliable ways to model general rate processes. Finally, simple model calculations illustrate the utility and implications of new rate expressions.

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

Renormalization of states and quasiparticles in many-body downfolding

We explore the principles of many-body Hamiltonian complexity reduction via downfolding on an effective low-dimensional representation. We show that the renormalization factor provides a unique measure of the quality of the compression as it directly represents the projection between the approximate stationary state of the many-body Hamiltonian and the full many-body wavefunction. Hence, the renormalization factor is a measure of fidelity between the effective (reduced-rank) description and the full many-body treatment for arbitrary (i.e., ground and excited) states. When the entire problem is mapped on a system of interacting quasiparticles [Romanova et al., npj Comput. Mater. 9, 126 (2023)], the effective Hamiltonians can faithfully reproduce the physics only when a clear energy scale separation exists between the subsystems and their environment. We also demonstrate that it is necessary to include quasiparticle renormalization at distinct energy scales, capturing the distinct interaction between subsystems and their surrounding environments. Numerical results from simple, exactly solvable models highlight the limitations and strengths of this approach, particularly for ground and low-lying excited states. This work lays the groundwork for applying dynamical downfolding techniques to problems concerned with (quantum) interfaces.

Green-functions technique↗
DOE OSTI2020

First-Principles Quantum Transport Modeling of Spin-Transfer and Spin-Orbit Torques in Magnetic Multilayers

A unified approach for computing (i) spin-transfer torque in magnetic trilayers like spin valve and magnetic tunnel junction, where injected charge current flows perpendicularly to interfaces, and (ii) spin-orbit torque in magnetic bilayers of the type ferromagnet/spin-orbit-coupled material, where injected charge current flows parallel to the interface, was reviewed. The experimentally explored and technologically relevant spin-orbit-coupled materials include 5d heavy metals, topological insulators, Weyl semimetals, and transition metal dichalcogenides. This approach requires to construct the torque operator for a given Hamiltonian of the device and the steady-state nonequilibrium density matrix, where the latter is expressed in terms of the nonequilibrium Green’s functions and split into three contributions. Tracing these contributions with the torque operator automatically yields field-like and damping-like components of spin-transfer torque or spin-orbit torque vector, which is particularly advantageous for spin-orbit torque where the direction of these components depends on the unknown-in-advance orientation of the current-driven nonequilibrium spin density in the presence of spin-orbit coupling. Illustrative examples are provided by computing spin-transfer torque in a one-dimensional toy model of a magnetic tunnel junction and realistic Co/Cu/Co spin valve, both of which are described by first-principles Hamiltonians obtained from noncollinear density functional theory calculations, as well as by computing spin-orbit torque in a ferromagnetic layer described by a tight-binding Hamiltonian which includes spin-orbit proximity effect within ferromagnetic monolayers assumed to be generated by the adjacent monolayer transition metal dichalcogenide. In addition, it is shown here how spin-orbit proximity effect, quantified by computing (via first-principles retarded Green’s function) spectral functions and spin textures on monolayers of realistic ferromagnetic material like Co in contact with heavy metal or monolayer transition metal dichalcogenide, can be tailored to enhance the magnitude of spin-orbit torque. Furthermore, errors made in the calculation of spin-transfer torque are quantified when using Hamiltonian from collinear density functional theory, with rigidly rotated magnetic moments to create noncollinear magnetization configurations, instead of proper (but computationally more expensive) self-consistent Hamiltonian obtained from noncollinear density functional theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Anisotropic spin-wave excitations in multiferroic BiFeO 3

Polarized inelastic neutron-scattering experiments have been performed to elucidate the anisotropic behavior of the low-energy spin-wave excitations in a multiferroic BiFeO 3 , which shows a cycloidal spin structure below 640 K. Using neutron polarization analysis for single magnetic domain crystals, magnetic excitation modes in and out of the cycloidal plane below 6 meV were separated successfully. The magnetic excitation spectra were analyzed using linear spin-wave theory. The low-energy magnon density of states consist of several magnon modes, including the two anisotropic modes, Φ and Ψ modes, distributed in and out of the cycloidal plane, respectively, which were previously observed using optical spectroscopies. Furthermore, there are other magnon modes that are not active in optical measurements. Additionally, a model spin Hamiltonian, which reproduces the spin-wave frequencies observed using optical spectroscopies, explains the overall spectra reasonably well.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Quantum Time Dynamics Mediated by the Yang–Baxter Equation and Artificial Neural Networks

Quantum computing shows great potential, but errors pose a significant challenge. This study explores new strategies for mitigating quantum errors using artificial neural networks (ANNs) and the Yang–Baxter equation (YBE). Unlike traditional error mitigation methods, which are computationally intensive, we investigate artificial error mitigation. We developed a novel method that combines ANNs for noise mitigation combined with the YBE to generate noisy data. This approach effectively reduces noise in quantum simulations, enhancing the accuracy of the results. The YBE rigorously preserves quantum correlations and symmetries in spin chain simulations in certain classes of integrable lattice models, enabling effective compression of quantum circuits while retaining linear scalability with the number of qubits. This compression facilitates both full and partial implementations, allowing the generation of noisy quantum data on hardware alongside noiseless simulations using classical platforms. By introducing controlled noise through the YBE, we enhance the data set for error mitigation. We train an ANN model on partial data from quantum simulations, demonstrating its effectiveness in mitigating errors in time-evolving quantum states, providing a scalable framework to enhance quantum computation fidelity, particularly in noisy intermediate-scale quantum (NISQ) systems. We demonstrate the efficacy of this approach by performing quantum time dynamics simulations using the Heisenberg XY Hamiltonian on real quantum devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗