Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MATRICES”

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 739 records · Page 41

Robust Millisecond Coherence Times of Erbium Electron Spins

Erbium doped solids are prime candidates for quantum memories in optical quantum networks given their telecom-compatible photon emission. An electron spin of erbium with millisecond coherence time is desirable for generating remote entanglement between adjacent quantum network nodes. Here, in this study, we report GHz-range electron spin transitions of 167 Er 3+ in yttrium oxide (Y 2 O 3 ) matrix with coherence times that are consistently longer than a millisecond. By polarizing para magnetic impurity spins we achieve a spin T 2 up to 1.46 ms, and up to 7.1 ms after dynamical decoupling. These coherence lifetimes are among the longest found for erbium in crystalline matrices despite the presence of host nuclear spins. We further enhance the coherence time beyond conventional dynamical decoupling, using customized sequences to simultaneously mitigate spectral diffusion and dipolar interactions. Our study not only establishes 167 Er 3+ : Y 2 O 3 as a significant quantum memory platform but also provides a guideline for engineering long-lived erbium spins in a variety of host materials for quantum technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simulation of Linear Non-Hermitian Boundary-Value Problems with Quantum Singular-Value Transformation

Herein we propose a quantum algorithm for simulating dissipative waves in inhomogeneous linear media as a boundary-value problem. Using the so-called quantum singular value transformation (QSVT), we construct a quantum circuit that models the propagation of electromagnetic waves in a one-dimensional system with outgoing boundary conditions. The corresponding measurement procedure is also discussed. Limitations of the QSVT algorithm are identified in connection with the large condition numbers that the dispersion matrices exhibit at weak dissipation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Determining topological order from infinite projected entangled pair states

In this work, we present a method of extracting information about the topological order from the ground state of a strongly correlated two-dimensional system computed with the infinite projected entangled pair state (iPEPS). For topologically ordered systems, the iPEPS wrapped on a torus becomes a superposition of degenerate, locally indistinguishable ground states. Projectors in the form of infinite matrix product operators onto states with well-defined anyon flux are used to compute topological S and T matrices (encoding mutual and self-statistics of emergent anyons). The algorithm is shown to be robust against a perturbation driving string-net toric code across a phase transition to a ferromagnetic phase. Our approach provides accurate results near quantum phase transition, where the correlation length is prohibitively large for other numerical methods. Moreover, we used numerically optimized iPEPS describing the ground state of the Kitaev honeycomb model in the toric code phase and obtained topological data in excellent agreement with theoretical prediction.

36 MATERIALS SCIENCE↗

Theory of the phase transition in random unitary circuits with measurements

Here, we present a theory of the entanglement transition tuned by measurement strength in qudit chains evolved by random unitary circuits and subject to either weak or random projective measurements. The transition can be understood as a nonanalytic change in the amount of information extracted by the measurements about the initial state of the system, quantified by the Fisher information. To compute the von Neumann entanglement entropy $\textit{S}$ and the Fisher information $\mathcal{F}$, we apply a replica method based on a sequence of quantities $\tilde{S}^{(n)}$ and $\mathcal{F}^{(n)}$ that depend on the $\textit{n}$th moments of density matrices and reduce to $\textit{S}$ and $\mathcal{F}$ in the limit $\textit{n}$ → 1. These quantities with $\textit{n}$ ≥ 2 are mapped to free energies of a classical spin model with $\textit{n}$! internal states in two dimensions with specific boundary conditions. In particular, $\tilde{S}^{(n)}$ is the excess free energy of a domain wall terminating on the top boundary, and $\mathcal{F}^{(n)}$ is related to the magnetization on the bottom boundary. Phase transitions occur as the spin models undergo ordering transitions in the bulk. Taking the limit of large local Hilbert space dimension $\textit{q}$ followed by the replica limit $\textit{n}$ → 1, we obtain the critical measurement probability $p_c$ = 1/2 and identify the transition as a bond percolation in the 2D square lattice in this limit. Finally, we show there is no phase transition if the measurements are allowed in an arbitrary nonlocal basis, thereby highlighting the relation between the phase transition and information scrambling. We establish an explicit connection between the entanglement phase transition and the purification dynamics of a mixed state evolution and discuss implications of our results to experimental observations of the transition and simulability of quantum dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Infrared fixed points of higher-spin fermions in topological semimetals

We determine the fate of interacting fermions described by the Hamiltonian H = p ∙ J in three-dimensional topological semimetals with linear band crossing, where p is momentum and J are the spin- j matrices for half-integer pseudospin j ≥ 3/2. While weak short-range interactions are irrelevant at the crossing point due to the vanishing density of states, weak long-range Coulomb interactions lead to a renormalization of the band structure. Using a self-consistent perturbative renormalization group approach, we show that band crossings of the type p ∙ J are unstable for j ≤ 7/2. Instead, through an intriguing interplay between cubic crystal symmetry, band topology, and interaction effects, the system is attracted to a variety of infrared fixed points. We also unravel several other properties of higher-spin fermions for general j , such as the relation between fermion self-energy and free energy, or the vanishing of the renormalized charge. An O(3) symmetric fixed point composed of equal chirality Weyl fermions is stable for j ≤ 7/2 and very likely so for all j. We then explore the rich fixed point structure for j = 5/2 in detail. We find additional attractive fixed points with enhanced O ( 3 ) symmetry that host both emergent Weyl or massless Dirac fermions, and identify a puzzling, infrared stable, anisotropic fixed point without enhanced symmetry in close analogy to the known case of j = 3/2.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Deconvolving the components of the sign problem

Auxiliary field quantum Monte Carlo simulations of interacting fermions require sampling over a Hubbard-Stratonovich field h introduced to decouple the interactions. The weight for a given configuration involves the products of the determinant of matrices $M_σ(h)$ where σ labels the species, and hence is typically not positive definite. Indeed, the average sign $\langle \mathscr {L} \rangle$ of the determinants goes to zero exponentially with increasing spatial size and decreasing temperature for most Hamiltonians of interest. This statement, however, does not explicitly separate two possible origins for the vanishing of $\langle \mathscr {L} \rangle$. Does $\langle \mathscr {L} \rangle$ → 0 because randomly chosen field configurations have det[M(h)] < 0, or does the sign problem arise because the specific subset of configurations chosen by the weighting function have a greater preponderance of negative values? In the latter case, the process of weighting the configurations with |det[M(h)]| might steer the simulation to a region of configuration space of h where positive and negative determinants are equally likely, even though randomly chosen h would preferentially have determinants with a single dominant sign. Here in this paper, we address the relative importance of these two mechanisms for the vanishing of $\langle \mathscr {L} \rangle$ in quantum simulations.

36 MATERIALS SCIENCE↗

Exact matrix product state representation and convergence of a fully correlated electronic wavefunction in the infinite-basis limit

Here In this paper we present the exact representation of a fully correlated electronic wavefunction as the single-particle basis approaches completeness. It consists of a half-infinite chain of matrices of exponentially increasing size. The complete basis limit is illustrated numerically using the density-matrix renormalization-group method by computing the core-valence entanglement in the C 2 ground state in increasing subsets of cc-pVTZ and pVQZ bases until convergence is reached.

36 MATERIALS SCIENCE↗

Integrable nonunitary quantum circuits

Here we show that the integrable Lindblad superoperators found recently can be used to build integrable nonunitary quantum circuits with two-site gates by demonstrating that the R(SIC) matrices are completely positive and trace preserving. Using the bond-site transformation, we obtain the corresponding integrable nonunitary quantum circuits with three-site gates. When restricted to the diagonals of the density matrix, these quantum models reduce to integrable classical cellular automata. This approach may pave the way for a systematic construction of circuits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantum phase transitions in a model Hamiltonian exhibiting entangled simultaneous fermion-pair and exciton condensations

Quantum states of a novel Bose-Einstein condensate, in which both fermion-pair and exciton condensations are simultaneously present, have recently been realized theoretically in a model Hamiltonian system. Here, in this study, we identify quantum phase transitions in that model between fermion-pair and exciton condensations based on a geometric analysis of the convex set of ground-state two-particle reduced density matrices (2-RDMs). The 2-RDM set provides a finite representation of the infinite parameter space of Hamiltonians that readily reveals a fermion-pair condensate phase and two distinct exciton condensate phases, as well as the emergence of first- and second-order phase transitions as the particle number of the system is increased. The set, furthermore, shows that the fermion-exciton condensate (FEC) lies along the second-order phase transition between the exciton and fermion-pair condensate phases. The detailed information about the exciton and fermion-pair phases, the forces behind these phases, as well as their associated transitions provides additional insight into the formation of the FEC condensate, which we anticipate will prove useful in its experimental realization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Windowed multipole representation of R -matrix cross sections

Nuclear cross sections are basic inputs to any nuclear computation. Campaigns of experiments are fitted with the parametric R-matrix model of quantum nuclear interactions, and the resulting cross sections are documented—both pointwise and as resonance parameters (with uncertainties)—in standard evaluated nuclear data libraries (ENDF, JEFF, BROND, JENDL, CENDL, TENDL): these constitute our common knowledge of fundamental low-energy nuclear cross sections. In the past decade, a collaborative effort has been deployed to establish a new nuclear cross-section library format—the Windowed Multipole Library—with the goal of considerably reducing the computational cost of cross-section calculations in nuclear transport simulations. This work lays the theoretical foundations underpinning these efforts. From general R-matrix scattering theory, we derive the windowed multipole representation of nuclear cross sections. Though physically and mathematically equivalent to R-matrix cross sections, the windowed multipole representation is particularly well suited for subsequent temperature treatment of angle-integrated cross sections, in particular Doppler broadening, which is the averaging of cross sections over the thermal motion of the target atoms. Doppler broadening is of critical importance in neutron transport applications, as it ensures the stability of many nuclear reactors (negative thermal reactivity). Yet, Doppler broadening of nuclear cross sections has been a considerable bottleneck for nuclear transport computations, often requiring memory-costly pretabulations. We show that the windowed multipole representation can perform accurate Doppler broadening analytically (up to the first reaction threshold), from which we derive cross-section temperature derivatives to any order—all computable on the fly (without precalculations stored in memory). Furthermore, we here establish a way of converting the R-matrix resonance parameters uncertainty (covariance matrices) into windowed multipole parameters uncertainty. We show that generating stochastic nuclear cross sections by sampling from the resulting windowed multipole covariance matrix can reproduce the cross-section uncertainty in the original nuclear data file. The windowed multipole representation is therefore a novel nuclear physics formalism able to generate Doppler broadened stochastic nuclear cross sections on the fly, unlocking breakthrough computational gains for nuclear computations. Through this foundational paper, we hope to make the windowed multipole representation accessible, reproducible, and usable for the nuclear physics community, as well as provide the theoretical basis for future research on expanding its capabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Noniterative finite amplitude methods for E 1 and M 1 giant resonances

The finite amplitude method (FAM) is a very efficient approach for solving the fully self-consistent random-phase approximation (RPA) equations. We use FAM to rederive the RPA matrices for Skyrme functionals, calculate the electric dipole ( E 1) and the magnetic dipole ( M 1) giant resonances, and compare the results with available experimental and evaluated data. For the E 1 transitions in heavy nuclei, the calculations reproduce well the resonance energy of the photoabsorption cross sections. In the case of M 1 transitions, we show that the residual interaction does not affect the transition strength of double-magic nuclei, which suggests that the spin terms in the Skyrme force currently neglected in the present computation could improve the agreement between FAM and experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Examining the justification for the introduction of a fermion localization function

Becke and Edgecombe suggested in 1990 a theoretical tool to describe electron localization in atoms and molecules, an idea which was borrowed by a large number of nuclear theorists since 2011 to describe nucleon localization in nuclear systems. I argue here that these arguments are highly questionable and cannot be used in interacting systems where effects beyond the naive mean field or the simple Hartree-Fock framework are important and the inclusion of correlations induced by particle interactions is necessary in order to introduce such a localization function. Furthermore, I also describe several aspects of the exchange and irreducible two-body density matrices, which depend on the character and strength of the two-particle interaction and which can be useful in justifying the derivation of an appropriate energy density functional.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Using DAPPER to extract the photon strength function of 58 Fe using the inverse Oslo and shape methods

The photon strength function of 58 Fe has been extracted using both the Oslo and Shape methods from particle–γ coincidence data measured using the Detector Array for Photons, Protons, and Exotic Residues, which probes nuclei utilizing (d,p) reactions in inverse kinematics. Four particle–γ coincidence matrices, each constructed with different treatments of the γ–ray energies, are explored in order to observe the impact on the resulting nuclear level density and photon strength. The final photon strength function reported is found to agree well with previous Oslo measurements of other iron isotopes. Systematic uncertainties are included, using different model parameters and their reported errors to perform the Oslo method normalization. The model-independent Shape method is explored and the functional form of the photon strength function obtained is in agreement with the Oslo method results. A low-energy enhancement is not reported for 58 Fe in this work given possible subtraction issues originating from strongly populated states.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Khuri-Treiman equations for 3$π$ decays of particles with spin

Khuri-Treiman equations have proven to be a useful theoretical tool in the analysis of three-body decays, especially into the 3$π$ final state. In this work we present in full detail the necessary generalization of the formalism to study the decays of particles with arbitrary spin, parity, and charge conjugation. To this extent, we find it most convenient to work with helicity amplitudes instead of the so-called invariant amplitudes, especially when dealing with the unitarity relations. The isobar expansions in the three possible ( s-, t-, and u- ) final channels are related with the appropriate crossing matrices. We pay special attention to the kinematical singularities and constraints of the helicity amplitudes, showing that these can be derived by means of the crossing matrix.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Data compression and covariance matrix inspection: Cosmic shear

Covariance matrices are among the most difficult pieces of end-to-end cosmological analyses. In principle, for two-point functions, each component involves a four-point function, and the resulting covariance often has hundreds of thousands of elements. We investigate various compression mechanisms capable of vastly reducing the size of the covariance matrix in the context of cosmic shear statistics. This helps identify which of its parts are most crucial to parameter estimation. We start with simple compression methods, by isolating and “removing” 200 modes associated with the lowest eigenvalues, then those with the lowest signal-to-noise ratio, before moving on to more sophisticated schemes like compression at the tomographic level and, finally, with the massively optimized parameter estimation and data compression (MOPED). We find that, while most of these approaches prove useful for a few parameters of interest, like Ω m , the simplest yield a loss of constraining power on the intrinsic alignment (IA) parameters as well as S 8 . For the case considered—cosmic shear from the first year of data from the Dark Energy Survey—only MOPED was able to replicate the original constraints in the 16-parameter space. Finally, we apply a tolerance test to the elements of the compressed covariance matrix obtained with MOPED and confirm that the IA parameter A IA is the most susceptible to inaccuracies in the covariance matrix.

79 ASTRONOMY AND ASTROPHYSICS↗

Simple Hamiltonian for quantum simulation of strongly coupled $(2+1)D$ SU(2) lattice gauge theory on a honeycomb lattice

Here, we find a simple spin Hamiltonian to describe physical states of $(2+1)$-dimensional SU(2) lattice gauge theory on a honeycomb lattice with a truncation of the electric field representation at $j_{max}=\frac{1}{2}$. The simple spin Hamiltonian contains only local products of Pauli matrices, even though Gauss’s law has been completely integrated out.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Lattice QCD study of π Σ − K ¯ N scattering and the Λ ( 1405 ) resonance

A lattice QCD computation of the coupled channel π Σ – K ¯ N scattering amplitudes in the Λ ( 1405 ) region is detailed. Results are obtained using a single ensemble of gauge field configurations with N f = 2 + 1 dynamical quark flavors and m π ≈ 200 MeV and m K ≈ 487 MeV . Hermitian correlation matrices using both single baryon and meson-baryon interpolating operators for a variety of different total momenta and irreducible representations are used. Several parametrizations of the two-channel scattering K -matrix are utilized to obtain the scattering amplitudes from the finite-volume spectrum. The amplitudes, continued to the complex energy plane, exhibit a virtual bound state below the π Σ threshold and a resonance pole just below the K ¯ N threshold. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Relativistic spin dynamics for vector mesons

We propose a relativistic theory for spin density matrices of vector mesons based on Kadanoff-Baym equations in the closed-time-path formalism. The theory puts the calculation of spin observables such as the spin density matrix element ρ 00 for vector mesons on a solid ground. Within the theory we formulate ρ 00 for ϕ mesons into a factorization form in separation of momentum and spacetime variables. We argue that the main contribution to ρ 00 at lower energies should be from the ϕ fields that can polarize the strange quark and antiquark in the same way as electromagnetic fields. The key observation is that there is correlation inside the ϕ meson wave function between the ϕ field that polarizes the strange quark and that polarizes the strange antiquark. This is reflected by the fact that the contributions to ρ 00 are all in squares of fields that are nonvanishing even if the fields may strongly fluctuate in spacetime. The fluctuation of strong force fields can be extracted from ρ 00 of unflavored vector mesons as links to fundamental properties of quantum chromodynamics. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗