Engineering Papers⌕ Search

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.

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

At least 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

QCD analysis of Δ S = 0 hadronic parity violation

We present a QCD analysis of the effective weak Hamiltonian at hadronic energy scales for strangeness-nonchanging (ΔS = 0) hadronic processes. Performing a leading-order renormalization group analysis in QCD from the W to the O(2 GeV) energy scale, we derive the pertinent effective Hamiltonian for hadronic parity violation, including the effects of both neutral and charged weak currents. We compute the complete renormalization group evolution of all isosectors and the evolution through heavy-flavor thresholds for the first time. We show that the additional four-quark operators that enter below the W mass scale from QCD operator mixing effects form a closed set, and they result in a 12 x 12 anomalous dimension matrix. Computing the resulting effective Hamiltonian and comparing to earlier results, we affirm the importance of operator mixing effects and find, as an example, that the parity-violating pion-nucleon coupling constant, using the factorization Ansatz and an assessment of the pertinent quark charge of the nucleon in lattice QCD at the 2 GeV scale, is in better agreement with recent experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Comparing light-front quantization with instant-time quantization

Here in this paper we compare light-front quantization and instant-time quantization both at the level of operators and at the level of their Feynman diagram matrix elements. At the level of operators light-front quantization and instant-time quantization lead to equal light-front time commutation (or anticommutation) relations that appear to be quite different from equal instant-time commutation (or anticommutation) relations. Despite this we show that at unequal times instant-time and light-front commutation (or anticommutation) relations actually can be transformed into each other, with it only being the restriction to equal times that makes the commutation (or anticommutation) relations appear to be so different. While our results are valid for both bosons and fermions, for fermions there are subtleties associated with tip of the light cone contributions that need to be taken care of. At the level of Feynman diagrams we show for non-vacuum Feynman diagrams that the pole terms in four-dimensional light-front Feynman diagrams reproduce the widely used three-dimensional light-front on-shell Hamiltonian Fock space formulation in which the light-front energy and light-front momentum are on shell. Moreover, we show that the contributions of pole terms in non-vacuum instant-time and non-vacuum light-front Feynman diagrams are equal. However, because of circle at infinity contributions we show that this equivalence of pole terms fails for four-dimensional light-front vacuum tadpole diagrams. Then, and precisely because of these circle at infinity contributions, we show that light-front vacuum tadpole diagrams are not only nonzero, they quite remarkably are actually equal to the pure pole term instant-time vacuum tadpole diagrams. Light-front vacuum diagrams are not correctly describable by the on-shell Hamiltonian formalism, and thus not by the closely related infinite momentum frame prescription either. Thus for the light-front vacuum sector we must use the off-shell Feynman formalism as it contains information that is not accessible in the on-shell Hamiltonian Fock space approach. We show that light-front quantization is intrinsically nonlocal, and that for fermions this nonlocality is present in Ward identities. One can project fermion spinors into so-called good and bad components, and both of these components contribute in Ward identities. Central to our analysis is that the transformation from instant-time coordinates and fields to light-front coordinates and fields is a unitary, spacetime-dependent translation. Consequently, not only are instant-time quantization and light-front quantization equivalent, because of general coordinate invariance they are unitarily equivalent.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Guiding centre motion for particles in a ponderomotive magnetostatic end plug

The Hamiltonian dynamics of a single particle in a rotating plasma column, interacting with an magnetic multipole is perturbatively solved for up to second order, using the method of Lie transformations. First, the exact Hamiltonian is expressed in terms of canonical action-angle variables, and then an approximate integrable Hamiltonian is introduced, using another set of actions and angles, which describe the centre of oscillation for the particle. The perturbation introduces an effective ponderomotive potential, which to leading order is positive. At the second order, the pseudopotential consists of a sum of terms of the Miller form, and can have either sign. Additionally, at second order, the ponderomotive interaction introduces a modification to the particle effective mass, when considering the motion along the column axis. It is found that particles can be axially confined by the ponderomotive potentials, but acquire radial excursions which scale as the confining potential. The radial excursions of the particle along its trajectory are investigated, and a condition for the minimal rotation frequency for which the particle remains radially confined is derived. Last, we comment on the changes to the aforementioned solution to the pseudopotentials and particle trajectory in the case of resonant motion, that is, a motion which has the same periodicity as the perturbation.

Rubin, T. (ORCID:0000000164850096)↗
DOE OSTI2024

Simulating nonlinear optical processes on a superconducting quantum device

Simulating plasma physics on quantum computers is difficult because most problems of interest are nonlinear, but quantum computers are not naturally suitable for nonlinear operations. In weakly nonlinear regimes, plasma problems can be modelled as wave–wave interactions. In this paper, we develop a quantization approach to convert nonlinear wave–wave interaction problems to Hamiltonian simulation problems. We demonstrate our approach using two qubits on a superconducting device. Unlike a photonic device, a superconducting device does not naturally have the desired interactions in its native Hamiltonian. Nevertheless, Hamiltonian simulations can still be performed by decomposing required unitary operations into native gates. To improve experimental results, we employ a range of error-mitigation techniques. Apart from readout error mitigation, we use randomized compilation to transform undiagnosed coherent errors into well-behaved stochastic Pauli channels. Moreover, to compensate for stochastic noise, we rescale exponentially decaying probability amplitudes using rates measured from cycle benchmarking. We carefully consider how different choices of product-formula algorithms affect the overall error and show how a trade-off can be made to best utilize limited quantum resources. This study provides an example of how plasma problems may be solved on near-term quantum computing platforms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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

Machine-learning-assisted insight into spin ice Dy 2 Ti 2 O 7

Complex behavior poses challenges in extracting models from experiment. An example is spin liquid formation in frustrated magnets like Dy 2 Ti 2 O 7 . Understanding has been hindered by issues including disorder, glass formation, and interpretation of scattering data. Here, we use an automated capability to extract model Hamiltonians from data, and to identify different magnetic regimes. This involves training an autoencoder to learn a compressed representation of three-dimensional diffuse scattering, over a wide range of spin Hamiltonians. The autoencoder finds optimal matches according to scattering and heat capacity data and provides confidence intervals. Validation tests indicate that our optimal Hamiltonian accurately predicts temperature and field dependence of both magnetic structure and magnetization, as well as glass formation and irreversibility in Dy 2 Ti 2 O 7 . The autoencoder can also categorize different magnetic behaviors and eliminate background noise and artifacts in raw data. Our methodology is readily applicable to other materials and types of scattering problems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Spin–spin interactions in defects in solids from mixed all-electron and pseudopotential first-principles calculations

Abstract Understanding the quantum dynamics of spin defects and their coherence properties requires an accurate modeling of spin-spin interaction in solids and molecules, for example by using spin Hamiltonians with parameters obtained from first principles calculations. We present a real-space approach based on density functional theory for the calculation of spin-Hamiltonian parameters, where only selected atoms are treated at the all-electron level, while the rest of the system is described with the pseudopotential approximation. Our approach permits calculations for systems containing more than 1000 atoms, as demonstrated for defects in diamond and silicon carbide. We show that only a small number of atoms surrounding the defect needs to be treated at the all-electron level, in order to obtain an overall all-electron accuracy for hyperfine and zero-field splitting tensors. We also present results for coherence times, computed with the cluster correlation expansion method, highlighting the importance of accurate spin-Hamiltonian parameters for quantitative predictions of spin dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Partial dynamical symmetry versus quasi dynamical symmetry examination within a quantum chaos analyses of spectral data for even–even nuclei

Statistical analyses of the spectral distributions of rotational bands in 51 deformed prolate even–even nuclei in the 152 ≤ A ≤ 250 mass region $R_{{4_{1}^{ + } /2_{1}^{ + } }} \ge 3.00$ are examined in terms of nearest neighbor spacing distributions. Specifically, the focus is on data for 0 + , 2 + , and 4 + energy levels of the ground, gamma, and beta bands. The chaotic behavior of the gamma band, especially the position of the $2_{\gamma }^{ + }$ band-head compared to other levels and bands, is clear. The levels are analyzed within the framework of two models, namely, a SU(3)-partial dynamical symmetry Hamiltonian and a SU(3) two-coupled quasi-dynamical symmetry Hamiltonian, with results that are further analyzed using random matrix theory. The partial and quasi dynamics both yield outcomes that are in reasonable agreement with the known experimental results. However, due to the degeneracy of the beta and gamma bands within the simplest SU(3) picture, the theory cannot be used to describe the fluctuation properties of excited bands. By changing relative weights of the different terms in the partial and quasi dynamical Hamiltonians, results are obtained that show more GOE-like statistics in the partial dynamical formalism as the strength of the pairing term is increased. Also, in the quasi-dynamical symmetry limit, more correlations are found because of the stronger couplings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Quantum algorithm for electronic band structures with local tight-binding orbitals

While the main thrust of quantum computing research in materials science is to accurately measure the classically intractable electron correlation effects due to Coulomb repulsion, designing optimal quantum algorithms for simpler problems with well-understood solutions is a useful tactic to advance our quantum “toolbox”. With this in mind, we consider the quantum calculation of a periodic system’s single-electron band structure over a path through reciprocal space. Previous efforts have used the Variational Quantum Eigensolver algorithm to solve the energy of each band, which involves numerically optimizing the parameters of a variational quantum circuit to minimize a cost function, constructed as the expectation value of a Hamiltonian operator. Traditionally, a unique Hamiltonian operator is constructed for each k-point, so that many cost functions, each with their own parameter space, must be optimized to generate a single band. Similarly, calculating higher bands than the first has traditionally involved modifying the cost function with additional overlap terms to ensure higher-energy eigenstates are orthogonal to those of lower bands. In this paper, we adopt a direct space approach, using a novel hybrid first/second-quantized qubit mapping which allows us to construct a single Hamiltonian, and a single cost-function, suitable for solving the entire electronic band structure. In contrast to previous approaches, the k-point and the band index are selected by additional parameters in our quantum circuit, rather than through modifications to the cost function. The result is a technically and conceptually simpler approach to band structure calculations on a quantum computer. Moreover, we expect that the tools developed herein will motivate new strategies for tackling highly-correlated materials beyond the grasp of classical computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Grand canonical inverse design of multicomponent colloidal crystals

Inverse design methods are powerful computational approaches for creating colloidal systems which self-assemble into a target morphology by reverse engineering the Hamiltonian of the system. Despite this, these optimization procedures tend to yield Hamiltonians which are too complex to be experimentally realized. An alternative route to complex structures involves the use of several different components, however, conventional inverse design methods do not explicitly account for the possibility of phase separation into compositionally distinct structures. We present an inverse design scheme for multicomponent colloidal systems by combining active learning with a method to directly compute their ground state phase diagrams. This explicitly accounts for phase separation and can locate stable regions of Hamiltonian parameter space which grid-based surveys are prone to miss. Using this we design low-density, binary structures with Lennard-Jones-like pairwise interactions that are simpler than in the single component case and potentially realizable in an experimental setting. This reinforces the concept that ground states of simple, multicomponent systems might be rich with previously unappreciated diversity, enabling the assembly of non-trivial structures with only few simple components instead of a single complex one.

36 MATERIALS SCIENCE↗
DOE OSTI2020

PT-symmetry entails pseudo-Hermiticity regardless of diagonalizability

We prove that in finite dimensions, a Parity-Time (PT)-symmetric Hamiltonian is necessarily pseudo-Hermitian regard-less of whether it is diagonalizable or not. This result is different from Mostafazadeh’s, which requires the Hamiltonian to be diagonalizable. PT-symmetry breaking often occurs at exceptional points where the Hamiltonian is not diagonalizable. Our result implies that PT-symmetry breaking is equivalent to the onset of instabilities of pseudo-Hermitian systems, which was systematically studied by Krein et al. in 1950s. In particular, we show that the mechanism of PT-symmetry breaking is the resonance between two eigen modes with opposite signs of actions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Continuum electromagnetic gyrokinetic simulations of turbulence in the tokamak scrape-off layer and laboratory devices

Here, we present algorithms and results from Gkeyll, a full-f continuum, electromagnetic gyrokinetic code, designed to study turbulence in the edge region of fusion devices. The edge is computationally very challenging, requiring robust algorithms that can handle large-amplitude fluctuations and stable interactions with plasma sheaths. We present an energy-conserving high-order discontinuous Galerkin scheme that solves gyrokinetic equations in Hamiltonian form. Efficiency is improved by a careful choice of basis functions and automatically generated computation kernels. Previous verification tests were performed in the straight-field-line large plasma device [Shi et al., J. Plasma Phys. 83, 905830304 (2017)] and the Texas Helimak, a simple magnetized torus [Bernard et al., Phys. Plasmas 26, 042301 (2019)], including the effect of end-plate biasing on turbulence. Results for the scrape-off layer for NSTX parameters with a model helical magnetic geometry with bad curvature have been obtained [Shi et al., Phys. Plasmas 26, 012307 (2019)]. In this paper, we present algorithms for the two formulations of electromagnetic gyrokinetics: the Hamiltonian and the symplectic. We describe each formulation and show results of benchmark tests. Although our scheme works for the Hamiltonian formulation, the presence of spurious numerical modes for high-β and large $k^2_⊥ρ^2_s$ regimes shows that the symplectic formulation is more robust. We then review our recent algorithm for the symplectic formulation [Mandell et al., J. Plasma Phys. 86, 905860109 (2020)], along with example application of this new capability. Maintaining positivity of the distribution function can be challenging, and we describe a new and novel exponential recovery based algorithm to address this.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Multimode two-dimensional vibronic spectroscopy. II. Simulating and extracting vibronic coupling parameters from polarization-selective spectra

We report experimental demonstrations of polarization-selection two-dimensional Vibrational-Electronic (2D VE) and 2D Electronic-Vibrational (2D EV) spectroscopies aim to map the magnitudes and spatial orientations of coupled electronic and vibrational coordinates in complex systems. The realization of that goal depends on our ability to connect spectroscopic observables with molecular structural parameters. In this paper, we use a model Hamiltonian consisting of two anharmonically coupled vibrational modes in electronic ground and excited states with linear and bilinear vibronic coupling terms to simulate polarization-selective 2D EV and 2D VE spectra. We discuss the relationships between the linear vibronic coupling and two-dimensional Huang–Rhys parameters and between the bilinear vibronic coupling term and Duschinsky mixing. We develop a description of the vibronic transition dipoles and explore how the Hamiltonian parameters and non-Condon effects impact their amplitudes and orientations. Using simulated polarization-selective 2D EV and 2D VE spectra, we show how 2D peak positions, amplitudes, and anisotropy can be used to measure parameters of the vibronic Hamiltonian and non-Condon effects. This paper, along with the first in the series, provides the reader with a detailed description of reading, simulating, and analyzing multimode, polarization-selective 2D EV and 2D VE spectra with an emphasis on extracting vibronic coupling parameters from complex spectra.

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

Interpretation of multiple solutions in fully iterative GF2 and GW schemes using local analysis of two-particle density matrices

Due to the presence of non-linear equations, iterative Green’s function methods can result in multiple different solutions even for simple molecular systems. In contrast to the wave-function methods, a detailed and careful analysis of such molecular solutions was not performed before. In this work, we use two-particle density matrices to investigate local spin and charge correlators that quantify the charge resonance and covalent characters of these solutions. When applied within the unrestricted orbital set, spin correlators elucidate the broken symmetry of the solutions, containing necessary information for building effective magnetic Hamiltonians. Based on GW and GF2 calculations of simple molecules and transition metal complexes, we construct Heisenberg Hamiltonians, four-spin-four-center corrections, and biquadratic spin–spin interactions. These Hamiltonian parameterizations are compared to previous wave-function calculations.

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

Influence of non-adiabatic effects on linear absorption spectra in the condensed phase: Methylene blue

Modeling linear absorption spectra of solvated chromophores is highly challenging as contributions are present both from coupling of the electronic states to nuclear vibrations and from solute–solvent interactions. In systems where excited states intersect in the Condon region, significant non-adiabatic contributions to absorption line shapes can also be observed. Here, we introduce a robust approach to model linear absorption spectra accounting for both environmental and non-adiabatic effects from first principles. This model parameterizes a linear vibronic coupling (LVC) Hamiltonian directly from energy gap fluctuations calculated along molecular dynamics (MD) trajectories of the chromophore in solution, accounting for both anharmonicity in the potential and direct solute–solvent interactions. The resulting system dynamics described by the LVC Hamiltonian are solved exactly using the thermalized time-evolving density operator with orthogonal polynomials algorithm (T-TEDOPA). The approach is applied to the linear absorption spectrum of methylene blue in water. We show that the strong shoulder in the experimental spectrum is caused by vibrationally driven population transfer between the bright S1 and the dark S2 states. The treatment of the solvent environment is one of many factors that strongly influence the population transfer and line shape; accurate modeling can only be achieved through the use of explicit quantum mechanical solvation. The efficiency of T-TEDOPA, combined with LVC Hamiltonian parameterizations from MD, leads to an attractive method for describing a large variety of systems in complex environments from first principles.

Dunnett, Angus J. (ORCID:0000000285797826)↗
DOE OSTI2023

Polarization consistent dielectric screening in polarizable continuum model calculations of solvation energies

A polarization consistent framework, where dielectric screening is affected consistently in polarizable continuum model (PCM) calculations, is employed for the study of solvation energies. The computational framework combines a screened range-separated-hybrid functional (SRSH) with PCM calculations, SRSH-PCM, where dielectric screening is imposed in both PCM self-consistent reaction field (SCRF) iterations and the electronic structure Hamiltonian. We begin by demonstrating the impact of modifying the Hamiltonian to include such dielectric screening in SCRF iterations by considering the solutions of electrostatically embedded Hartree–Fock (HF) exact exchange equations. Long-range screened HF-PCM calculations are shown to capture properly the linear dependence of gap energy of frontier orbitals on the inverse of the dielectric constant, whereas unscreened HF-PCM orbital energies are fallaciously semi-constant with respect to the dielectric constant and, therefore, inconsistent with the ionization energy gaps. Similar trends affect density functional theory (DFT) calculations that aim to achieve predictive quality. Importantly, the dielectric screened calculations are shown to significantly affect DFT- and HF PCM-based solvation energies, where screened solvation energies are smaller compared to the unscreened values. Importantly, SRSH-PCM, therefore, appears to reduce the tendency of DFT-PCM to overestimate solvation energies, where we find the effect to increase with the dielectric constant and the polarity of the molecular solute, trends that enhance the quality of DFT-PCM calculations of solvation energy. Understanding the relationship of dielectric screening in the Hamiltonian and DFT-PCM calculations can ultimately benefit on-going efforts for the design of predictive and parameter free descriptions of solvation energies.

Chemistry↗
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 OSTI2023

Restoring permutational invariance in the Jordan–Wigner transformation

The Jordan–Wigner transformation is a powerful tool for converting systems of spins into systems of fermions, or vice versa. While this mapping is exact, the transformation itself depends on the labelling of the spins. One consequence of this dependence is that approximate solutions of a Jordan–Wigner-transformed Hamiltonian may depend on the (physically inconsequential) labelling of the spins. In this work, we turn to an extended Jordan–Wigner transformation which remedies this problem and which may also introduce some correlation atop the Hartree–Fock solution of a transformed spin Hamiltonian. We demonstrate that this extended Jordan–Wigner transformation can be thought of as arising from a unitary version of the Lie algebraic similarity transformation (LAST) theory. Here, we show how these ideas, particularly in combination with the standard (non-unitary) version of LAST, can provide a potentially powerful tool for the treatment of the XXZ and J 1 –J 2 Heisenberg Hamiltonians.

74 ATOMIC AND MOLECULAR PHYSICS↗