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

Influence of hydrostatic pressure on hidden order, the Kondo lattice, and magnetism in URu 2 Si 2-x P x

Within the chemical substitution series URu 2 Si 2-x P x , there is an evolution in the ground-state behavior from hidden ordered (HO) for x ≲ 0.03, to Kondo lattice behavior with no ordering (NO) for 0.03 ≲ x ≲ 0.26, to antiferromagnetism (AFM-2) for 0.26 ≲ x ≲ 0.5 [A. Gallagher et al., Nat. Commun. 7, 10712 (2016); A. Gallagher et al., J. Phys.: Condens. Matter 29, 024004 (2016)]. To better understand what factors control this behavior, temperature-dependent electrical resistivity measurements are performed for this series under applied pressures P up to 20.5 kbar. Specimens in the HO x region show similarities to the parent compound, where HO transforms into antiferromagnetism (AFM-1) at a critical pressure (P c ). P c decreases with increasing x and collapses towards P = 0 near x ≈ 0.03, suggesting that AFM-1 occurs at ambient pressure for this concentration. No pressure-induced phase transitions are observed in the NO x region and the AFM-2 state is only weakly suppressed by P. Measurements further reveal that AFM-1 and AFM-2 are distinct from each other. Calculations of the wave functions using the tight-binding Hartree-Fock approximation are performed and show (i) that the radial probability distributions for the phosphorus ions are more tightly bound than those for the silicon and (ii) that the energy difference between the orbitals decreases with increasing x. The cumulative effect of these two factors is that Si → P substitution decreases the hybridization strength, which correlates with the weakening of HO. At large x, additional effects such as electrical charge tuning also play an important role in determining the ground-state

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Unconventional bipartite entanglement in the quantum dimer magnet Yb2Be 2 SiO 7

The quantum dimer magnet, with antiferromagnetic intradimer and interdimer Heisenberg exchange between spin-1/2 moments, is known to host an (|$\uparrow$$\downarrow$$\rangle$ - |$\downarrow$$\uparrow$$\rangle$)/$\sqrt{2}$ singlet ground state when the intradimer exchange is dominant. Rare-earth-based quantum dimer systems with strong spin-orbit coupling offer the opportunity for tuning their magnetic properties by using magnetic anisotropy as a control knob. Here, we present bulk characterization and neutron scattering measurements of the quantum dimer magnet Yb 2 Be 2 SiO 7 . We find that the Yb3+ ions can be described by an effective spin-1/2 model at low temperatures and the system does not show signs of magnetic order down to 50 mK. The magnetization, heat capacity, and neutron spectroscopy data can be well-described by an isolated dimer model with highly anisotropic exchange that stabilizes a singlet ground state with a wavefunction (|$\uparrow$$\uparrow$$\rangle$ - |$\downarrow$$\downarrow$$\rangle$)/$\sqrt{2}$ or (|$\uparrow$$\uparrow$$\rangle$ + |$\downarrow$$\downarrow$$\rangle$)/$\sqrt{2}$. Our results show that strong spin-orbit coupling can induce unusual entangled states of matter in quantum dimer magnets.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Magnetic field dependence of the spin fluctuations in CeCu 5.8 ⁢Ag 0.2

Quantum phase transitions are among the most intriguing phenomena that can occur when the electronic ground state of correlated metals are tuned by external parameters such as pressure, magnetic field, or chemical substitution. Such transitions between distinct states of matter are driven by quantum fluctuations, and can give rise to macroscopically coherent phases that are at the forefront of condensed matter research. However, the nature of the critical fluctuations, and thus the fundamental physics controlling many quantum phase transitions, remain poorly understood in numerous strongly correlated metals. Here we study the model material CeCu 5.8⁢ Ag 0.2 to gain insight into the implications of critical fluctuations originating from different regions in reciprocal space. By employing an external magnetic field along the crystallographic 𝑎 and 𝑐 axis as auxiliary tuning parameter, we observe a pronounced anisotropy in the suppression of the quantum critical fluctuations, reflecting the spin anisotropy of the long-range ordered ground state at larger silver concentration. Coupled with the temperature dependence of the quantum fluctuations, these results suggest that the quantum phase transition in CeCu 5.8⁢ Ag 0.2 is driven by three-dimensional spin-density wave fluctuations.

Boraley, Xavier [Paul Scherrer Inst. (PSI), Villig↗

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↗