Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Effective 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 577 records · Page 32

Partial fillings of the bosonic $E$ 8 quantum Hall state

Here, we study bosonic topological phases constructed from electrons. In addition to a bulk excitation energy gap, these bosonic phases also have a fermion energy gap, below which all local excitations in the bulk and on the edge are even combinations of electrons. We focus on chiral phases, in which all low-energy edge excitations move in the same direction, that arise from the short-range entangled $E_8$ quantum Hall state, the bosonic analog of the filled lowest Landau level of electrons. The $E_8$ edge-state theory features an $E_8$ Kac-Moody symmetry that can be decomposed into ${\cal G}_A \times {\cal G}_B$ subalgebras, such as $SU(3) \times E_6$, $SO(M) \times SO(16-M)$, and $G_2 \times F_4$. (Here, $\{SO(M) \}$, $\{SU(N)\}$, and $\{E_8, G_2, F_4 \}$ denote orthogonal, unitary, and exceptional Lie algebras.) Using these symmetry decompositions, we construct exactly solvable coupled-wire model Hamiltonians for families of long-range entangled ${\cal G}_A$ or ${\cal G}_B$ bosonic fractional quantum Hall states that "partially fill" the $E_8$ state and are pairwise related by a generalized particle-hole symmetry. These long-range entangled states feature either Abelian or non-Abelian topological order. Some support the emergence of non-local Dirac and Majorana fermions, Ising anyons, metaplectic anyons, Fibonacci anyons, as well as deconfined $\mathbb{Z}_2$ gauge fluxes and charges.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The Singlet–Triplet Gap of Cyclobutadiene: The CIPSI-Driven CC( P ; Q ) Study

An accurate determination of singlet−triplet gaps in biradicals, including cyclobutadiene in the automerization barrier region where one has to balance the substantial nondynamical many-electron correlation effects characterizing the singlet ground state with the predominantly dynamical correlations of the lowest-energy triplet, remains a challenge for many quantum chemistry methods. High-level coupled-cluster (CC) approaches, such as the CC method with a full treatment of singly, doubly, and triply excited clusters (CCSDT), are often capable of providing reliable results, but routine application of such methods is hindered by their high computational costs. We have recently proposed a practical alternative to converging the CCSDT energetics at small fractions of the computational effort, even when electron correlations become stronger and connected triply excited clusters are larger and nonperturbative, by merging the CC(P;Q) moment expansions with the selected configuration interaction methodology abbreviated as CIPSI. We demonstrate that one can accurately approximate the highly accurate CCSDT potential surfaces characterizing the lowest singlet and triplet states of cyclobutadiene along the automerization coordinate and the gap between them using tiny fractions of triply excited cluster amplitudes identified with the help of relatively inexpensive CIPSI Hamiltonian diagonalizations.

Basis sets↗

Magnetic excitations of the hybrid multiferroic (ND 4 ) 2 FeCl 5 ·D 2 O

We report a comprehensive inelastic neutron scattering study of the hybrid molecule-based multiferroic compound ( ND 4 ) 2 FeCl 5 · D 2 O in the zero-field incommensurate cycloidal phase and the high-field quasicollinear phase. The spontaneous electric polarization changes its direction concurrently with the field-induced magnetic transition, from mostly aligned with the crystallographic a axis to the c axis. To account for such a change in polarization direction, the underlying multiferroic mechanism was proposed to switch from the spin-current model induced via the inverse Dzyaloshinskii-Moriya interaction to the p - d hybridization model. We perform a detailed analysis of the inelastic neutron data of ( ND 4 ) 2 FeCl 5 · D 2 O using linear spin-wave theory to quantify magnetic interaction strengths and investigate the possible impact of different multiferroic mechanisms on the magnetic couplings. Overall our result reveals that the spin dynamics of both multiferroic phases can be well described by a Heisenberg Hamiltonian with easy-plane anisotropy. We do not find notable differences between the optimal model parameters of the two phases. The hierarchy of exchange couplings and the balance among frustrated interactions remain the same between two phases, suggesting that magnetic interactions in ( ND 4 ) 2 FeCl 5 · D 2 O are much more robust than the electric polarization in response to delicate reorganizations of the electronic degrees of freedom in an applied magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Study on the magnetothermal properties of the D y 1 - x T b x A l 2 series of compounds

Here, in this work, we developed a theoretical model Hamiltonian, in the mean field approximation, to describe the magnetic and magnetocaloric behavior of the series of compounds Dy 1-x Tb x Al 2 (x = 0.00, 0.15, 0.25, 0.30, 0.35, 0.40 and 0.75). We adjusted the exchange parameters λ DyDy , λ TbTb , and λ DyTb to obtain the spin reorientation temperatures (T SR ) and the critical temperature (T C ) for each compound in the series. The results obtained by the Hamiltonian model agree satisfactorily with the experimental results. The heat capacity curves with and without an applied magnetic field, adiabatic temperature variation and isothermal entropy variation were modeled and compared with the experimental data. As the experimental results show, our model was also able to reproduce the change in the spin reorientation process: a first order spin reorientation transition appears for concentrations x = 0.15, 0.25, 0.30, 0.35, and no spin reorientation transitions after x = 0.40.

36 MATERIALS SCIENCE↗

Magnetic properties of a staggered S = 1 chain with an alternating single-ion anisotropy direction

Materials composed of spin-1 antiferromagnetic (AFM) chains are known to adopt complex ground states that are sensitive to the single-ion-anisotropy (SIA) energy ( D ), and intrachain ( J 0 ) and interchain ( J 1 , 2 ′ ) exchange energy scales. While theoretical and experimental studies have extended this model to include various other energy scales, the effect of the lack of a common SIA axis is not well explored. Here we investigate the magnetic properties of Ni ( pyrimidine ) ( H 2 O ) 2 ( NO 3 ) 2 , a chain compound where the tilting of Ni octahedra leads to a twofold alternation of the easy-axis directions along the chain. Muon-spin relaxation measurements indicate a transition to long-range order at T N = 2.3 K and the magnetic structure is initially determined to be antiferromagnetic and collinear using elastic neutron diffraction experiments. Inelastic neutron scattering measurements were used to find J 0 = 5.107 ( 7 ) K , D = 2.79 ( 1 ) K , J 1 ′ = 0.00 ( 5 ) K , J 2 ′ = 0.18 ( 3 ) K , and a rhombic anisotropy energy E = 0.19 ( 9 ) K . Mean-field modeling reveals that the ground state structure hosts spin canting of ϕ ≈ 6 . 5 ∘ , which is not detectable above the noise floor of the elastic neutron diffraction data. Monte Carlo simulation of the powder-averaged magnetization, M ( H ) , is then used to confirm these Hamiltonian parameters, while single-crystal M ( H ) simulations provide insight into features observed in the data. Published by the American Physical Society 2025

Vaidya, S. (ORCID:000000033975497X)↗

Classical Benchmarks for Variational Quantum Eigensolver Simulations of the Hubbard Model

Simulating the Hubbard model is of great interest to a wide range of applications within condensed matter physics, however its solution on classical computers remains challenging in dimensions larger than one. The relative simplicity of this model, embodied by the sparseness of the Hamiltonian matrix, allows for its efficient implementation on quantum computers, and for its approximate solution using variational algorithms such as the variational quantum eigensolver. While these algorithms have been shown to reproduce the qualitative features of the Hubbard model, their quantitative accuracy in terms of producing true ground state energies and other properties, and the dependence of this accuracy on the system size and interaction strength, the choice of variational ansatz, and the degree of spatial inhomogeneity in the model, remains unknown. Here we present a rigorous classical benchmarking study, demonstrating the potential impact of these factors on the accuracy of the variational solution of the Hubbard model on quantum hardware, for systems with up to 32 qubits. We find that even when using the most accurate wavefunction ansätze for the Hubbard model, the error in its ground state energy and wavefunction plateaus for larger lattices, while stronger electronic correlations magnify this issue. Concurrently, spatially inhomogeneous parameters and the presence of off-site Coulomb interactions only have a small effect on the accuracy of the computed ground state energies. Our study highlights the capabilities and limitations of current approaches for solving the Hubbard model on quantum hardware, and we discuss potential future avenues of research.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

High field magnetization of FePS 3

High field magnetization measurements in pulsed fields up to 65 T have been performed on FePS 3 , which is nominally a good example of a two-dimensional Ising-like antiferromagnet on a honeycomb lattice. Measurements with the field parallel to the moment direction confirm the presence of two first-order transitions above 35 T, to M/M sat = 1/2 and M/M sat = 1, respectively, at 4 K. The measurements are in contradiction with published estimates for the magnetic exchange parameters, but the contradiction can be resolved by allowing for anisotropic exchange parameters in the Hamiltonian. Furthermore, the magnetization with the field perpendicular to the moment direction is anisotropic, with no transitions observed for fields along the a axis while a cascade of first-order transitions is observed for fields above 50 T along the b axis, the latter case also showing a strong degradation of the sample after repeated pulses. The results indicate a strong magnetolattice coupling in FePS 3 . Temperature-dependent measurements hint at a possible tricritical point.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Characterization and thermometry of dissipatively stabilized steady states

In this work we study the properties of dissipatively stabilized steady states of noisy quantum algorithms, exploring the extent to which they can be well approximated as thermal distributions, and proposing methods to extract the effective temperature T. We study an algorithm called the relaxational quantum eigensolver (RQE), which is one of a family of algorithms that attempt to find ground states and balance error in noisy quantum devices. In RQE, we weakly couple a second register of auxiliary ‘shadow’ qubits to the primary system in Trotterized evolution, thus engineering an approximate zero-temperature bath by periodically resetting the auxiliary qubits during the algorithm’s runtime. Balancing the infinite temperature bath of random gate error, RQE returns states with an average energy equal to a constant fraction of the ground state. We probe the steady states of this algorithm for a range of base error rates, using several methods for estimating both T and deviations from thermal behavior. In particular, we both confirm that the steady states of these systems are often well-approximated by thermal distributions, and show that the same resources used for cooling can be adopted for thermometry, yielding a fairly reliable measure of the temperature. These methods could be readily implemented in near-term quantum hardware, and for stabilizing and probing Hamiltonians where simulating approximate thermal states is hard for classical computers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Spin pumping from antiferromagnetic insulator spin-orbit-proximitized by adjacent heavy metal: a first-principles Floquet-nonequilibrium Green function study

Motivated by the recent experiment on spin pumping from sub-THz radiation-driven uniaxial antiferromagnetic insulator (AFI) MnF 2 into heavy metal (HM) Pt hosting strong spin-orbit (SO) coupling, we compute and compare pumped spin currents in Cu/MnF 2 /Cu and Pt/MnF 2 /Cu heterostructures. Recent theories of spin pumping by AFI have relied on simplistic Hamiltonians (such as tight-binding) and the scattering approach to quantum transport yielding the so-called interfacial spin mixing conductance (SMC), but the concept of SMC ceases to be applicable when SO coupling is present directly at the interface. In contrast, we use a more general first-principles quantum transport approach which combines noncollinear density functional theory with Floquet-nonequilibrium Green's functions in order to take into account: SO-proximitized AFI as a new type of quantum material, different from isolated AFI and brought about by AFI hybridization with adjacent HM layer; strong SO coupling at the interface; and evanescent wavefunctions penetrating from Pt or Cu into AFI layer to make its interfacial region conducting rather than insulating as in the isolated AFI. The DC component of pumped spin current $I_\mathrm{DC}^{S_z}$ vs. precession cone angle $\theta_{\boldsymbol{l}}$ of the Néel vector l of AFI does not follow putative $I^{S_z}_\mathrm{DC} \propto \sin^2 \theta_{\boldsymbol{l}}$, except for very small angles $\theta_{\boldsymbol{l}} \lesssim 10^\circ$ for which we define an effective SMC from the prefactor and find that it doubles from MnF2/Cu to MnF2/Pt interface. In addition, the angular dependence $I^{S_z}_\mathrm{DC}(\theta_{\boldsymbol{l}})$ differs for opposite directions of precession of the Néel vector, leading to twice as large SMC for the right-handed than for the left-handed chirality of the precession mode.

36 MATERIALS SCIENCE↗

Variational Self-Consistent Theory for Beam-Loaded Cavities

A new variational theory is presented for beam loading in microwave cavities. The beam{field interaction is formulated as a dynamical interaction whose stationarity according to Hamilton’s principle will naturally lead to steady-state solutions that indicate how a cavity’s resonant frequency, Q and optimal coupling coefficient will detune as a result of the beam loading. A driven cavity Lagrangian is derived from first principles, including the effects of cavity wall losses, input power and beam interaction. The general formulation is applied to a typical klystron input cavity to predict the appropriate detuning parameters required to maximize the gain (or modulation depth) in the average Lorentz factor boost, (Δ γ). Numerical examples are presented, showing agreement with the general detuning trends previously observed in the literature. The developed formulation carries several advantages for beam-loaded cavity structures. It provides a self-consistent model for the dynamical (nonlinear) beam{field interaction, a procedure for maximizing gain under beam-loading conditions, and a useful set of parameters to guide cavity-shape optimization during the design of beam-loaded systems. Enhanced clarity of the physical picture underlying the problem seems to be gained using this approach, allowing straightforward inclusion or exclusion of different field configurations in the calculation and expressing the final results in terms of measurable quantities. Two field configurations are discussed for the klystron input cavity, using finite magnetic confinement or no confinement at all. Formulating the problem in a language that is directly accessible to the powerful techniques found in Hamiltonian dynamics and canonical transformations may potentially carry an additional advantage in terms of analytical computational gains, under suitable conditions

43 PARTICLE ACCELERATORS↗

Massless and massive photons within basis light-front quantization

We compute the structure function, transverse momentum dependent parton distributions (TMDs), and generalized parton distributions (GPDs) for the photon from the light-front quantum electrodynamics (QED) Hamiltonian, determined for its constituent bare photon and electron-positron Fock components within the basis light-front quantization framework. While the structure function and GPDs align well with expectations, some oscillatory behavior in the TMDs emerges as a finite-basis artifact. We also compute the GPD and the corresponding structure function for the longitudinal polarization mode of the massive photon. We further study the photon mass dependence of the calculated observables. After performing nonperturbative renormalization, we obtain a good quality description of the perturbative QED properties of the photon, as well as good agreement with the experimental data for the massless photon’s structure function.

Ab initio calculations↗

Molecular structure determination: Equilibrium structure of pyrimidine (m-C 4 H 4 N 2 ) from rotational spectroscopy (reSE) and high-level ab initio calculation (re) agree within the uncertainty of experimental measurement

The pure rotational spectrum of pyrimidine (m-C 4 H 4 N 2 ), the meta-substituted dinitrogen analog of benzene, has been studied in the millimeter-wave region from 235 GHz to 360 GHz. The rotational spectrum of the ground vibrational state has been assigned and fit to yield accurate rotational and distortion constants. Over 1700 distinct transitions were identified for the normal isotopologue in its ground vibrational state and least-squares fit to a partial sextic S-reduced Hamiltonian. Transitions for all four singly substituted 13 C and 15 N isotopologues were observed at natural abundance and were likewise fit. Deuterium-enriched samples of pyrimidine were synthesized, giving access to all eleven possible deuterium-substituted isotopologues, ten of which were previously unreported. Experimental values of rotational constants and computed values of vibration–rotation interaction constants and electron-mass corrections were used to determine semi-experimental equilibrium structures (r e SE ) of pyrimidine. The r e SE structure obtained using coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] corrections shows exceptional agreement with the re structure computed at the CCSD(T)/cc-pCV5Z level (≤0.0002 Å in bond distance and ≤0.03° in bond angle). Of the various computational methods examined, CCSD(T)/cc-pCV5Z is the only method for which the computed value of each geometric parameter lies within the statistical experimental uncertainty (2σ) of the corresponding semi-experimental coordinate. The exceptionally high accuracy and precision of the structure determination is a consequence of the large number of isotopologues measured, the precision and extent of the experimental frequency measurements, and the sophisticated theoretical treatment of the effects of vibration–rotation coupling and electron mass. Taken together, these demanding experimental and computational studies establish the capabilities of modern structural analysis for a prototypical monocyclic aromatic compound.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Conservative Spin-Magnitude Change in Orbital Evolution in General Relativity

We show that physical scattering observables for compact spinning objects in general relativity can depend on additional degrees of freedom in the spin tensor beyond those described by the spin vector alone. The impulse, spin kick, and leading-order waveforms exhibit such a nontrivial dependence. A signal of this additional structure is the change in the magnitude of the spin vector under conservative Hamiltonian evolution, similar to our previous studies in electrodynamics. These additional degrees of freedom describe dynamical mass multipoles of compact objects and decouple for black holes. We also show that the conservative impulse, spin kick, and change of the additional degrees of freedom are encoded in the eikonal phase.

Classical black holes↗

On the S-matrix of Ising field theory in two dimensions

We explore the analytic structure of the non-perturbative S-matrix in arguably the simplest family of massive non-integrable quantum field theories: the Ising field theory (IFT) in two dimensions, which may be viewed as the Ising CFT deformed by its two relevant operators, or equivalently, the scaling limit of the Ising model in a magnetic field. Our strategy is that of collider physics: we employ Hamiltonian truncation method (TFFSA) to extract the scattering phase of the lightest particles in the elastic regime, and combine it with S-matrix bootstrap methods based on unitarity and analyticity assumptions to determine the analytic continuation of the 2 → 2 S-matrix element to the complex s-plane. Focusing primarily on the “high temperature” regime in which the IFT interpolates between that of a weakly coupled massive fermion and the E 8 affine Toda theory, we will numerically determine 3-particle amplitudes, follow the evolution of poles and certain resonances of the S-matrix, and exclude the possibility of unknown wide resonances up to reasonably high energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Approaching the self-dual point of the sinh-Gordon model

One of the most striking but mysterious properties of the sinh-Gordon model (ShG) is the b → 1/b self-duality of its S-matrix, of which there is no trace in its Lagrangian formulation. Here b is the coupling appearing in the model’s eponymous hyperbolic cosine present in its Lagrangian, cosh(bΦ). In this paper we develop truncated spectrum methods (TSMs) for studying the sinh-Gordon model at a finite volume as we vary the coupling constant. We obtain the expected results for b << 1 and intermediate values of b, but as the self-dual point b = 1 is approached, the basic application of the TSM to the ShG breaks down. We find that the TSM gives results with a strong cutoff Ec dependence, which disappears according only to a very slow power law in Ec. Standard renormalization group strategies — whether they be numerical or analytic — also fail to improve upon matters here. We thus explore three strategies to address the basic limitations of the TSM in the vicinity of b = 1. In the first, we focus on the small-volume spectrum. We attempt to understand how much of the physics of the ShG is encoded in the zero mode part of its Hamiltonian, in essence how ‘quantum mechanical’ vs ‘quantum field theoretic’ the problem is. In the second, we identify the divergencies present in perturbation theory and perform their resummation using a supra-Borel approximate. In the third approach, we use the exact form factors of the model to treat the ShG at one value of b as a perturbation of a ShG at a different coupling. In the light of this work, we argue that the strong coupling phase b > 1 of the Lagrangian formulation of model may be different from what is naïvely inferred from its S-matrix. In particular, we present an argument that the theory is massless for b > 1.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Simple relativistic quark models

A class of phenomenological relativistic models of hadronic systems motivated by quantum chromodynamics that have dual representations as models of mesons and nucleons or quarks and gluons is investigated. These models are designed to provided qualitative insight into the role of sea quarks in hadronic structure and reactions. The model assumption is that the Hamiltonian can be divided into two parts: one that involves degrees of freedom in the same connected local and global color singlet and the remaining interactions that allow the connected local and global color singlets to interact. The first class of interactions results in infinite towers of bare “particles” with hadronic quantum numbers. All but a finite number of these remain stable when the second class of interactions is included. The model interactions are expressed in terms of subhadronic degrees of freedom, which determine the bare hadronic spectrum and the interactions involving the bare hadrons in terms of a small number of subhadronic model parameters. As a first test, this paper considers the simplest case of mesons that interact via a string-breaking interaction. One virtue of this model is that all of the bare meson masses and eigenfunctions can be computed analytically. In addition, the string-breaking interaction leads to production vertices that can also be computed analytically. The relativistic wave functions have a light-front kinematic symmetry. The goal is to find a simple relativistic quantum mechanical model based on subhadronic degrees of freedom that can provide an efficient, qualitatively consistent description of hadronic masses, lifetimes, cross sections, sea quark effects, and electromagnetic properties. The simplicity of the model makes it a potentially useful tool to study the impact of sea quarks on hadronic structure and reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Enhancing scalability of a matrix-free eigensolver for studying many-body localization

We propose several techniques to enhance the parallel scalability of a matrix-free eigensolver designed for studying many-body localization (MBL) of quantum spin chain models with nearest-neighbor interactions and on-site disorder. This type of problem is computationally challenging because the dimension of the associated Hamiltonian matrix grows exponentially with respect to the number of spins L, and we need to average over different realizations of the random disorder to obtain relevant statistical behavior. For each disorder realization, we need to compute eigenvalues from different regions of the spectrum and their corresponding eigenvectors. In previous work, the interior eigenstates for a single eigenvalue problem are computed via the shift-and-invert Lanczos algorithm. Due to the extremely high memory footprint of the LU factorizations, this technique is not well suited for large L’s. For example, we need thousands of compute nodes on modern high performance computing infrastructures to go beyond L = 24. The matrix-free approach does not suffer from this memory bottleneck, however, its scalability is limited by a computation and communication load imbalance. To reduce this imbalance and to significantly enhance the scalability of the matrix-free eigensolver, we reorder the matrix and leverage the consistent space runtime, CSPACER. We also show its efficiency in managing irregular communication patterns at scale compared to optimized MPI non-blocking two-sided and one-sided RMA implementation variants. This effort enables us to study MBL for spin chains with a larger number of spins. The efficiency and effectiveness of the proposed algorithm is demonstrated by computing eigenstates on a massively parallel many-core high performance computer.

METIS↗

Improved Fast Randomized Iteration Approach to Full Configuration Interaction

In this work, we present three modifications to our recently introduced fast randomized iteration method for full configuration interaction (FCI-FRI) and investigate their effects on the method’s performance for Ne, H 2 O, and N 2 . The initiator approximation, originally developed for full configuration interaction quantum Monte Carlo, significantly reduces statistical error in FCI-FRI when few samples are used in compression operations, enabling its application to larger chemical systems. The semistochastic extension, which involves exactly preserving a fixed subset of elements in each compression, improves statistical efficiency in some cases but reduces it in others. We also developed a new approach to sampling excitations that yields consistent improvements in statistical efficiency and reductions in computational cost. We discuss possible strategies based on our findings for improving the performance of stochastic quantum chemistry methods more generally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗