Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “GROUND STATES”

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 1,063 records · Page 59

Error mitigation, optimization, and extrapolation on a trapped-ion testbed

Current noisy intermediate-scale quantum (NISQ) trapped-ion devices are subject to errors which can significantly impact the accuracy of calculations if left unchecked. A form of error mitigation called zero noise extrapolation (ZNE) can decrease an algorithm’s sensitivity to these errors without increasing the number of required qubits. Here we explore different methods for integrating this error mitigation technique into the Variational Quantum Eigensolver (VQE) algorithm for calculating the ground state of the HeH + molecule at 0.8 Å in the presence of experimental noise. Using the Quantum Scientific Computing Open User Testbed (QSCOUT) trapped-ion device, we test three methods of scaling noise for extrapolation: time stretching the two-qubit gates, scaling the sideband detuning parameter, and inserting two-qubit gate identity operations into the ansatz circuit. We find that time stretching and sideband detuning scaling fail to scale the noise on our particular hardware in a way that can be extrapolated to zero noise. Scaling our noise with global gate identity insertions and extrapolating after variational optimization, we achieve error suppression of 96.8%, resulting in an energy estimate within –0.004 ± 0.04 hartree of the ground state energy. This is an improvement, but still outside the chemical accuracy threshold of 0.0016 hartree. Furthermore, our results show that the efficacy of this error mitigation technique depends on choosing the correct implementation for a given device architecture.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Finite-temperature density matrix embedding theory

We describe a formulation of the density matrix embedding theory at finite temperature. We present a generalization of the ground-state bath orbital construction that embeds a mean-field finite-temperature density matrix up to a given order in the Hamiltonian, or the Hamiltonian up to a given order in the density matrix. In this work, we assess the performance of the finite-temperature density matrix embedding on the one-dimensional Hubbard model both at half-filling and away from it, and the two-dimensional Hubbard model at half-filling, comparing to exact data where available, as well as results from finite-temperature density matrix renormalization group, dynamical mean-field theory, and dynamical cluster approximations. The accuracy of finite-temperature density matrix embedding appears comparable to that of the ground-state theory, with, at most, a modest increase in bath size, and competitive with that of cluster dynamical mean-field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Towards Accurate Predictions of Martensitic Transition Temperatures for Shape Memory Alloys from Ab Initio Simulations

Experimentally, NiTi undergoes a single martensitic phase transition around 341 K from the lowtemperature (T) monoclinic B19’ phase (P21/m) to the high-temperature cubic B2 phase (Pm3m). Theoretically, an orthorhombic B33 (Cmcm) has also been proposed as the T=0 ground state structure, although this phase has never been observed. Accurate predictions of martensitic transition temperatures (MTT) have remained elusive in part due to several well-known theoretical complexities of these systems including low temperature instabilities of the B2 phase. Recently, we proposed a rigorous thermodynamic integration approach based on ab initio simulations to resolve many of these difficulties [1,2]. However, an unsatisfying overprediction of the MTT relative to experiment (by ~100 K) means a fully quantitative theory is still lacking. In this work, we report several new developments to our method that bring first principles theory and experiment much closer into agreement. Our calculations indicate that phonon free energies at low temperature stabilizes B19’ over B33, rationalizing B19’ as the ground state down to T=0. We also find that accurate computations of the electronic free energy, i.e. the change in energy and the appearance of electronic configurational entropy due to finite temperature, is crucial to obtain accurate MTT. Incorporating these corrections results in an MTT prediction of 365 K for binary NiTi, which is in very close agreement with experiment. Our theoretical approach is expected to be a broadly applicable and predictive theory for MTT of SMAs.

Zhigang Wu↗

Nuclear data for reactor production of 131 Ba and 133 Ba

The newest radioisotope for brachytherapy treatment of prostate cancer is 131 Cs (t 1/2 = 9.69 d, 100% EC). Generated via electron capture decay of 131 Ba (t 1/2 = 11.6 d, 100% EC), 131 Cs has been used in brachytherapy for prostate cancer since 2004. The 131 Ba parent is produced through neutron capture of enriched 130 Ba in a nuclear reactor. For large-scale production of 131 Ba, an accurate knowledge of production and burnup cross sections of 131 Ba are essential. Here, we report two group cross sections (thermal and resonance integrals) for 130 Ba and 131 Ba and a new measure of the half-life of 131 Ba. Targets consisting of milligram quantities of enriched 130Ba (~35%) were irradiated in Oak Ridge National Laboratory's High Flux Isotope Reactor at thermal and resonance neutron fluxes of (1.9–2.1) × 10 15 and (5.8–7.0) × 10 13 neutrons·cm -2 s -1 , respectively, for durations ranging from 3 to 26 days. In addition, cadmium covered samples of 130 Ba were irradiated for 1 hour at 12.6% full reactor power (10.7 MW). The yield of 131 Ba approaches a saturation value of ~60 GBq (~1.6 Ci) per mg of 130 Ba for 20 days irradiation at a thermal neutron flux of 1.8 × 10 15 n·s -1 ·cm -2 , with a thermal/epithermal ratio of ~30. Under the above experimental conditions, the two group cross sections of 130 Ba are 6.9 ± 0.5 b (thermal, σ 0 ) and 173 ± 7 b (resonance, I 0 ). These values represent the sum of cross sections to metastable and ground states of 131 Ba. For 131 Ba, the empirically measured thermal cross section is 200 ± 50 b assuming an I 0 /σ 0 of 10. This cross section is reported for the first time. Further, the half-life of 131 Ba was remeasured to be 11.657 ± 0.008 d. Lastly, this study also resulted in the co-production of 133 Ba (t 1/2 = 10.52 y, 100% EC). The experimental yield of 133 Ba is ~370 MBq (~10 mCi) per mg of 132 Ba (thin target) for one cycle irradiation in the High Flux Isotope Reactor, and measured two-group 132 Ba cross sections are 7.2 ± 0.2 b and 39.9 ± 1.3 b. These values also represent the sum of cross sections to metastable and ground states of 133 Ba.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Imaginary Time Propagation on a Quantum Chip

We report evolution in imaginary time is a prominent technique for finding the ground state of quantum many-body systems, and the heart of a number of numerical methods that have been used with great success in quantum chemistry, condensed matter, and nuclear physics. We propose an algorithm to implement imaginary time propagation on a quantum computer. Our algorithm is devised in the context of an efficient encoding into an optimized gate, drawing on the underlying characteristics of the quantum device of a unitary operation in an extended Hilbert space. However, we prove that for simple problems it can also be successfully applied to standard digital quantum machines. This work paves the way for porting quantum many-body methods based on imaginary-time propagation to near-term quantum devices, enabling the future quantum simulation of the ground states of a broad class of microscopic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ferrimagnetic spin waves in honeycomb and triangular layers of Mn{sub 3}Si{sub 2}Te{sub 6}.

A detailed analysis of the ferrimagnetic ground state of Mn 3 Si 2 Te 6 has been performed using inelastic neutron scattering. Although the proposed valence of the nominal Mn 2+ ions would have quenched orbital angular momentum, a significant exchange anisotropy exists in Mn 3 Si 2 Te 6 . This apparent exchange anisotropy is a manifestation of a weak spin-orbit coupling in the layered material. We employ a detailed simulation of the spin-wave spectrum coupling traditional refinement of dispersion parameters to image analysis techniques, while including Monte Carlo simulations of the instrumental resolution to accurately identify the exchange couplings to the third nearest neighbor. An independent validation of our results is made by comparing our final Hamiltonian to heat capacity measurements.

Sala, G.↗

Anisotropic field-induced changes in the superconducting order parameter of UTe 2

UTe 2 is a newly discovered unconventional superconductor, where electron Cooper pairs combine into a spin-triplet ground state. We study the specific heat C( H,T ) of a high-quality UTe 2 single crystal, with a single anomaly at a high superconducting transition temperature T c = 2.1 K and a low residual Sommerfeld coefficient γ 0 . Applying magnetic fields up to 12 T along the crystallographic axes a, b , and c , we observe highly anisotropic behavior in γ 0 ( H ). Below 4 T, γ 0 ( H ) ≈ γ 0 + λi √H (i = a, b, c) , which is expected for an unconventional superconductor with nodes on the Fermi surface away from the direction of the magnetic field. At 4 T ( a and b axes) and 5.5 T ( c axis), distinctive changes in the behavior of γ 0 (H) suggest that a stable low-field ground state begins to evolve above 4 T. The data analysis indicates that the nodes are located near the a–b plane, with possible order parameters d B2u + iεd B1u or d B2u + iεd Au . Our findings support B 2u as the primary superconducting order parameter in UTe 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Grid-based diffusion Monte Carlo for fermions without the fixed-node approximation

A diffusion Monte Carlo algorithm is introduced that can determine the correct nodal structure of the wave function of a few-fermion system and its ground-state energy without an uncontrolled bias. This is achieved by confining signed random walkers to the points of a uniform infinite spatial grid, allowing them to meet and annihilate one another to establish the nodal structure without the fixed-node approximation. An imaginary-time propagator is derived rigorously from a discretized Hamiltonian, governing a non-Gaussian, sign-flipping, branching, and mutually annihilating random walk of particles. The accuracy of the resulting stochastic representations of a fermion wave function is limited only by the grid and imaginary-time resolutions and can be improved in a controlled manner. Here, the method is tested for a series of model problems including fermions in a harmonic trap as well as the He atom in its singlet or triplet ground state. For the latter case, the energies approach from above with increasing grid resolution and converge within 0.015 E h of the exact basis-set-limit value for the grid spacing of 0.08 a.u. with a statistical uncertainty of 10 –5 E h without an importance sampling or Jastrow factor.

74 ATOMIC AND MOLECULAR PHYSICS↗

Static and dynamic properties of atomic nuclei with high-resolution potentials

Here, we compute ground-state and dynamical properties of 4 He and 16 O nuclei using as input high-resolution, phenomenological nucleon-nucleon and three-nucleon forces that are local in coordinate space. The nuclear Schrodinger equation for both nuclei is accurately solved employing the auxiliary-field diffusion Monte Carlo approach. For the 4 He nucleus, detailed benchmarks are carried out with the hyperspherical harmonics method. In addition to presenting results for the binding energies and radii, we also analyze the momentum distributions of these nuclei and their Euclidean response function corresponding to the isoscalar density transition. The latter quantity is particularly relevant for lepton-nucleus scattering experiments, as it paves the way to quantum Monte Carlo calculations of electroweak response functions of 16 O.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Interfaces and the extended Hilbert space of Chern-Simons theory

The low energy effective field theories of (2 + 1) dimensional topological phases of matter provide powerful avenues for investigating entanglement in their ground states. In [1] the entanglement between distinct Abelian topological phases was investigated through Abelian Chern-Simons theories equipped with a set of topological boundary conditions (TBCs). In the present paper we extend the notion of a TBC to non-Abelian Chern-Simons theories, providing an effective description for a class of gapped interfaces across non-Abelian topological phases. These boundary conditions furnish a defining relation for the extended Hilbert space of the quantum theory and allow the calculation of entanglement directly in the gauge theory. Because we allow for trivial interfaces, this includes a generic construction of the extended Hilbert space in any (compact) Chern-Simons theory quantized on a Riemann surface. Additionally, this provides a constructive and principled definition for the Hilbert space of effective ground states of gapped phases of matter glued along gapped interfaces. Lastly, we describe a generalized notion of surgery, adding a powerful tool from topological field theory to the gapped interface toolbox.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantum Many-Body Theory from a Solution of the N -Representability Problem

Here, in this study, we present a many-body theory based on a solution of the N-representability problem in which the ground-state two-particle reduced density matrix (2-RDM) is determined directly without the many-particle wave function. We derive an equation that re-expresses physical constraints on higher-order RDMs to generate direct constraints on the 2-RDM, which are required for its derivation from an N-particle density matrix, known as N-representability conditions. The approach produces a complete hierarchy of 2-RDM constraints that do not depend explicitly upon the higher RDMs or the wave function. By using the two-particle part of a unitary decomposition of higher order constraint matrices, we can solve the energy minimization by semidefinite programming in a form where the low-rank structure of these matrices can be potentially exploited. We illustrate by computing the ground-state electronic energy and properties of the H 8 ring.

74 ATOMIC AND MOLECULAR PHYSICS↗

Towards utility-scale electronic structure with sample-based quantum bootstrap embedding

One of the main applications for which quantum computers are hoped to find utility is in simulating ground state energies and other observables of molecular chemical systems. The recently proposed sample-based diagonalization method is a readily implementable method for this task on current-day hardware using short circuit depths and has been demonstrated on as many as 85 qubits in recent studies. In this work, we combine the recently proposed quantum bootstrap embedding (QBE) method with sampled-based diagonalization (QBE-SQD) and present the first benchmarking study of the QBE method on real quantum hardware, ibm_pittsburgh, a Heron r3 processor with 156 qubits. Our test system is a hydrogen ring with 8 hydrogen atoms in the cc-pVDZ basis. We show that for this system, QBE-SQD using an active space of (8e, 19o) per fragment with a 43 qubit footprint produces a ground state energy accuracy which exceeds that of an SQD calculation with an (8e, 30o) active space with a 67 qubit footprint when using a comparable number of Slater determinants. This demonstrates that the use of quantum bootstrap embedding techniques is a promising path towards extending the capabilities of state-of-the-art quantum eigensolvers on near-term devices.

Bierman, Joel [North Carolina State University, Ra↗

Influence of 73 Rb on the ashes of accreting neutron stars

We find that the proton separation energy, S(p), of 73 Rb is –640(40) keV, deduced from the observation of β-delayed ground-state protons following the decay of 73 Sr. This lower-limit determination of the proton separation energy of 73 Rb coupled with previous upper limits from nonobservation, provides a full constraint on the mass excess with ΔM ( 73 Rb) = –46.01 ± 0.04 MeV. With this new mass excess and the excitation energy of the J π = 5/2 – isobaric-analog state (T = 3/2) in 73 Rb, an improved constraint can be put on the mass excess of 73Sr using the isobaric-multiplet mass equation (IMME), and we find ΔM( 73 Sr) = –31.98 ± 0.37 MeV. Furthermore, these new data were then used to study the composition of ashes on accreting neutron stars following Type I x-ray bursts. Counterintuitively, we find that there should be an enhanced fraction of A > 102 nuclei with more negative proton separation energies at the 72 Kr rp-process waiting point. Larger impurities of heavier nuclei in the ashes of accreting neutron stars will impact the cooling models for such astrophysical scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗

Spin-orbit-controlled metal–insulator transition in Sr 2 IrO 4

In the context of correlated insulators, where electron–electron interactions (U) drive the localization of charge carriers, the metal–insulator transition is described as either bandwidth- or filling-controlled. Motivated by the challenge of the insulating phase in Sr 2 IrO 4 , a new class of correlated insulators has been proposed, in which spin–orbit coupling (SOC) is believed to renormalize the bandwidth of the half-filled j eff = 1/2 doublet, allowing a modest U to induce a charge-localized phase. Although this framework has been tacitly assumed, a thorough characterization of the ground state has been elusive. Furthermore, direct evidence for the role of SOC in stabilizing the insulating state has not been established, because previous attempts at revealing the role of SOC have been hindered by concurrently occurring changes to the filling. Here, we overcome this challenge by employing multiple substituents that introduce well-defined changes to the signatures of SOC and carrier concentration in the electronic structure, as well as a new methodology that allows us to monitor SOC directly. Specifically, we study Sr 2 Ir 1-x T x O 4 (T = Ru, Rh) by angle-resolved photoemission spectroscopy, combined with ab initio and supercell tight-binding calculations. This allows us to distinguish relativistic and filling effects, thereby establishing conclusively the central role of SOC in stabilizing the insulating state of Sr 2 IrO 4 . Most importantly, we estimate the critical value for SOC in this system to be λ c = 0.42 ± 0.01 eV, and provide the first demonstration of a spin–orbit-controlled metal–insulator transition.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Rotational excitations in rare-earth nuclei: A comparative study within three cranking models with different mean fields and treatments of pairing correlations

High-spin rotational bands in rare-earth Er (Z = 68), Tm (Z = 69), and Yb (Z = 70) isotopes are investigated by three different nuclear models. These are (i) the cranked relativistic Hartree-Bogoliubov (CRHB) approach with approximate particle number projection by means of the Lipkin-Nogami (LN) method, (ii) the cranking covariant density functional theory (CDFT) with pairing correlations treated by a shell-model-like approach (SLAP) or the so-called particle-number conserving (PNC) method, and (iii) cranked shell model (CSM) based on the Nilsson potential with pairing correlations treated by the PNC method. A detailed comparison between these three models in the description of the ground-state rotational bands of even-even Er and Yb isotopes is performed. The similarities and differences between these models in the description of the moments of inertia, the features of band crossings, and equilibrium deformations and pairing energies of the even-even nuclei under study are discussed. These quantities are considered as a function of rotational frequency and proton and neutron numbers. The changes in the properties of the first band crossings with increasing neutron number in this mass region are investigated. On average, a comparable accuracy of the description of available experimental data is achieved in these models. However, the differences between model predictions become larger above the first band crossings. Because of the time-consuming nature of numerical calculations in the CDFT-based models, a systematic study of the rotational properties of both ground-state and excited-state bands in odd-mass Tm nuclei is carried out only by the PNC-SCM. With few exceptions, the rotational properties of experimental 1-quasiparticle and 3-quasiparticle bands in 165,167,169,171 Tm are reproduced reasonably well. As a result, the appearance of backbendings or upbendings in these nuclei is well understood from the analysis of the variations of the occupation probabilities of the single-particle states and their contributions to total angular momentum alignment with rotational frequency.

150 ≤ A ≤ 189↗

Uniaxial Strain Control of Bulk Ferromagnetism in Rare-Earth Titanates

The perovskite rare-earth titanates are model Mott insulators with magnetic ground states that are very sensitive to structural distortions. These distortions couple strongly to the orbital degrees of freedom and, in principle, it should be possible to tune the superexchange and the magnetic transition with strain. We investigate the representative system (Y,La,Ca)TiO 3 , which exhibits low crystallographic symmetry and no structural instabilities. From magnetic susceptibility measurements of the Curie temperature, we demonstrate direct, reversible, and continuous control of ferromagnetism by influencing the TiO 6 octahedral tilts and rotations with uniaxial strain. The relative change in T C as a function of strain is well described by ab initio calculations, which provides detailed understanding of the complex interactions among structural, orbital, and magnetic properties in rare-earth titanates. The demonstrated manipulation of octahedral distortions opens up far-reaching possibilities for investigations of electron-lattice coupling, competing ground states, and magnetic quantum phase transitions in a wide range of quantum materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Spin fluctuations, absence of magnetic order, and crystal electric field studies in the Yb 3+ -based triangular lattice antiferromagnet Rb 3 ⁢Yb⁢(VO 4 ) 2

Here, we report a comprehensive experimental investigation of the structural, thermodynamic, static, and dynamic properties of a triangular lattice antiferromagnet Rb 3 ⁢Yb(VO 4 ) 2 . Through the analysis of magnetic susceptibility, magnetization, and specific heat, complemented by crystal electric field (CEF) calculations, we confirm the Kramers' doublet with effective spin 𝐽 eff = 1/2 ground state. Magnetic susceptibility and isothermal magnetization analysis reveal a weak antiferromagnetic interaction among the 𝐽 eff = 1/2 spins, characterized by a small Curie-Weiss temperature (𝜃$^{\textrm{LT}}_{\textrm{CW}}$ ≃−0.26 K) or a reduced exchange coupling (𝐽/𝑘 B ≃ 0.18 K). The 51 V NMR spectra and spin-lattice relaxation rate (1/𝑇 1 ) show no evidence of magnetic long-range-order down to 1.6 K but reflect strong influence of CEF excitations in the intermediate temperatures. At low temperatures, 1/𝑇 1 ⁡(𝑇) shows pronounced frequency dependence and 1/𝑇 1 vs field at different temperatures follows the scaling behavior, highlighting the role of paramagnetic fluctuations. The CEF calculations using the point charge approximation divulge a large energy gap ( ∼18.61 meV) between the lowest and second lowest energy doublets, further establishing Kramers' doublet as the ground state. Our calculations also reproduce the experimental magnetization and specific heat data and indicate an in-plane magnetic anisotropy. These findings position Rb 3⁢ Yb(VO 4 ) 2 as an ideal candidate to explore intrinsic quantum fluctuations and possible quantum spin-liquid physics in a Yb 3+ -based triangular lattice antiferromagnet.

Sebastian, Sebin J. [Indian Institute of Science E↗

Collective excitations and quantum incompressibility in electron-hole bilayers

In this study, we apply quantum continuum mechanics to the calculation of the excitation spectrum of a coupled electron-hole bilayer. The theory expresses excitation energies in terms of ground-state intra- and interlayer pair correlation functions, which are available from quantum Monte Carlo calculations. The final formulas for the collective modes deduced from this approach coincide with the formulas obtained in the “quasilocalized particle approximation” by Kalman et al. [G. Kalman, V. Valtchinov, and K. I. Golden, Phys. Rev. Lett. 82 , 3124 (1999) ], and likewise the theory predicts the existence of gapped excitations in the charged channel, with the gap arising from electron-hole correlation. An immediate consequence of the gap is that the static density-density response function of the charged channel vanishes as q 2 for wave vector q → 0 , rather than linearly in q , as commonly expected. In this sense, the system is incompressible . This feature, which has no analog in the classical electron-hole plasma, is consistent with the existence of an excitonic ground state and implies the existence of a discontinuity in the chemical potential of electrons and holes when the numbers of electrons and holes are equal. It should be experimentally observable by monitoring the densities of electrons and holes in response to potentials that attempt to change these densities in opposite directions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗