Quantum computers can efficiently solve problems which are widely believed to lie beyond the reach of classical computers. In the near-term, hybrid quantum-classical algorithms, which efficiently embed quantum hardware in classical frameworks, are crucial in bridging the vast divide in the performance of the purely-quantum algorithms and their classical counterparts. Here, a hybrid quantum-classical algorithm is presented for the computation of non-perturbative characteristics of quantum field theories. The presented algorithm relies on a universal parametrized quantum circuit ansatz based on Euler and Cartan's decompositions of single and two-qubit operators. It is benchmarked by computing the energy spectra of lattice realizations of quantum field theories with both short and long range interactions. Low depth circuits are provided for false vacua as well as highly excited states corresponding to mesonic and baryonic excitations occurring in the analyzed models. The described algorithm opens a hitherto-unexplored avenue for the investigation of mass-ratios, scattering amplitudes and false-vacuum decays in quantum field theories.
We provide practical simulation methods for scalar field theories on a quantum computer that yield improved asymptotics as well as concrete gate estimates for the simulation and physical qubit estimates using the surface code. We achieve these improvements through two optimizations. First, we consider a finite volume approach for estimating the elements of the S-matrix. This approach is appropriate in general for 1+1D and for certain low-energy elastic collisions in higher dimensions. Second, we implement our approach using a series of different fault-tolerant simulation algorithms for Hamiltonians formulated both in the field occupation basis and field amplitude basis. Our algorithms are based on either second-order Trotterization or qubitization. The cost of Trotterization in occupation basis scales as O(lambda N7Q3/(M5/2F3/2)) where ) is the coupling strength, N is the occupation cutoff, Qis the volume of the spatial lattice, M is the mass of the particles and F is the uncertainty in the energy calculation used for the S-matrix determination. Qubitization in the field basis scales as O(Q2(k2A + kM2)/E), where k is the cutoff in the field and A is a scaled coupling constant. We find in both cases that the bounds suggest physically meaningful simulations can be performed using on the order of 4 & times; 106 physical qubits and 1012 T-gates which corresponds to roughly one day on a superconducting quantum computer with surface code and a cycle time of 100 ns. This places the simulation of scalar field theory within striking distance of the gate counts for the best available chemistry simulation results.
Using matrix product state techniques we study the nonequilibrium dynamical response of the half-filled Hubbard ladder when subject to an optical pump. Optical pumping offers a way of producing and manipulating new strongly correlated phenomena by suppressing existing magnetic correlations. The ladder allows the effects of pump directionality to be investigated, and compared to a single chain it has strong spin-charge coupling and a fully gapped excitation spectrum, promising different nonequilibrium physics. We compute time-dependent correlations, including the nonequilibrium dynamical structure factors for spin and charge. By deriving a combined spin-charge sum rule that applies both in and out-of-equilibrium, we show that spectral weight is pumped directly from the antiferromagnetic spin response into a low energy omega similar to 0 charge response below the Mott gap. The transfer of weight is pump direction dependent: pumping directed along the legs disrupts magnetic correlations more than pumping in the rung-direction, even if the post-pump energy density is similar. The charge correlation length is dramatically enhanced by the pump, whilst the spin correlations are most strongly suppressed at nearest- and next-nearest-neighbor spacings. After the pump the system is in a nonthermal correlated metallic state, with gapless charge excitations and approximately equal spin and charge correlation lengths, emphasizing the importance of treating these degrees of freedom on an equal footing in nonequilibrium systems.
Conventional wisdom dictates that quantum effects become unimportant at high temperatures. In magnets, when the thermal energy exceeds interactions between atomic magnetic moments, the moments are usually uncorrelated, and classical paramagnetic behavior is observed. This thermal decoherence of quantum spin behaviors is a major hindrance to quantum information applications of spin systems. Remarkably, our neutron scattering experiments on Yb chains in an insulating perovskite crystal defy these conventional expectations. We find a sharply defined spectrum of spinons, fractional quantum excitations of spin-1/2 chains, to persist to temperatures much higher than the scale of the interactions between Yb magnetic moments. The observed sharpness of the spinon continuum's dispersive upper boundary indicates a spinon mean free path exceeding ≈ 35 inter-atomic spacings at temperatures more than an order of magnitude above the interaction energy scale. We thus discover an important and highly unique quantum behavior, which expands the realm of quantumness to high temperatures where entropy-governed classical behaviors were previously believed to dominate. Our results have profound implications for spin systems in quantum information applications operating at finite temperatures and motivate new developments in quantum metrology.
Symmetry-resolved entanglement, capturing the refined structure of quantum entanglement in systems with global symmetries, has attracted a lot of attention recently. In this manuscript, introducing the notion of symmetry-resolved generalized entropies, we aim to develop a computational framework suitable for the study of excited state symmetry-resolved entanglement as well as the dynamical evolution of symmetry-resolved entanglement in symmetry-preserving out-of-equilibrium settings. We illustrate our framework using the example of (1+1)-d free massless compact boson theory, and benchmark our results using lattice computation in the XX chain. As a byproduct, our computational framework also provides access to the probability distribution of the symmetry charge contained within a subsystem and the corresponding full counting statistics.
Physical quantities with long lifetimes have both theoretical significance in the study of quantum many-body systems and practical implications for quantum technologies. In this manuscript, we investigate the roles played by topological defects in the construction of quasiconserved quantities, using as a prototypical example the Kramers-Wannier duality defect in a deformed one-dimensional quantum transverse-field Ising model. We construct the duality defect Hamiltonian in three different ways: half-chain Kramers-Wannier transformation, utilization of techniques in the Ising fusion category, and defect-modified weak integrability breaking deformation. The third method is also applicable for the study of generic integrable defects under weak integrability breaking deformations. We also work out the deformation of defect-modified higher charges in the model and study their slower decay behavior. Furthermore, we consider the corresponding duality defect twisted deformed Floquet transverse-field Ising model and investigate the stability of the isolated zero mode associated with the duality defect in the integrable Floquet Ising model, under such weak integrability breaking deformation.
Topological symmetries, invertible and otherwise, play a fundamental role in the investigation of quantum field theories. Despite their ubiquitous importance across a multitude of disciplines ranging from string theory to condensed matter physics, controlled realizations of models exhibiting these symmetries in physical systems are rare. Quantum simulators based on engineered solid-state devices provide a novel alternative to conventional condensed matter systems for realizing these models. In this work, eigenstates of impurity Hamiltonians and loop operators associated with the topological symmetries for the Ising conformal field theory in two space-time dimensions are realized on IBM's Kingston simulator. The relevant states are created on the quantum device using a hybrid quantum-classical algorithm. The latter is based on a variation of the quantum approximate optimization algorithm ansatz combined with the quantum natural gradient optimization method. Signatures of the topological symmetry are captured by measuring correlation functions of different qubit operators with results obtained from the quantum device in reasonable agreement with those obtained from classical computations. The current work demonstrates the viability of noisy quantum simulators as platforms for investigating low-dimensional quantum field theories with direct access to observables that are often difficult to probe in conventional condensed matter experiments.
Production of hard X-ray via inverse Compton scattering at photon energies below 100 keV range aimed at potential applications in medicine and material research is reported. Experiments have been performed at the Brookhaven National Laboratory, Accelerator Test Facility, employing the counter collision of a 70 MeV, 0.3 nC electron beam with a near infra-red Nd: YAG laser (1064 nm wavelength) pulse containing ~ 100 mJ in a single shot basis. The radiation distribution of the scattered photon beam is assessed to be sufficiently quasi monochromatic to produce clear contrast from the Au K- edge at 80.7 keV.
We propose a protocol for detecting multipartite entanglement in itinerant many-body electronic systems using the quantum Fisher information (QFI). We establish a connection between the QFI and single-particle Green's functions by identifying a set of non-trivial witness operators. To construct these operators, we employ a doubling of the system wherein we introduce two identical copies of the original model. While the witness operator hops electrons between copies, the copies do not interact with one another. We apply this methodology to a finite-sized fermionic system and showcase its effectiveness in detecting entanglement in spinless itinerant electron models. We show that the detected entanglement level is sensitive to the wave vector associated with the hopping process. We also demonstrate the important role that symmetry has in detecting levels of entanglement. Our protocol paves the way for detecting entanglement in many-body systems using scanning tunneling microscopy and angle-resolved photoemission spectroscopy, thus offering exciting prospects beyond the detection of entanglement via the dynamical spin response accessed in neutron scattering experiments.
Although entanglement is both a central ingredient in our understanding of quantum many-body systems and an essential resource for quantum technologies, we only have a limited ability to quantify entanglement in real quantum materials. Thus far, entanglement metrology in quantum materials has been limited to measurements involving Hermitian operators, such as the detection of spin entanglement using inelastic neutron scattering. Here, we devise a method to extract the quantum Fisher information (QFI) from non-Hermitian operators and formulate an entanglement witness for resonant inelastic x-ray scattering (RIXS). Our approach is then applied to the model iridate dimer system Ba$_3$CeIr$_2$O$_9$ and used to directly test for entanglement of the electronic orbitals between neighboring Ir sites. We find that entanglement is challenging to detect under standard conditions, but that it could be achieved by analyzing the outgoing x-ray polarization or via specific choices of momentum and energy. Our protocol provides a new handle for entanglement detection, which offers routes to related types of entanglement witness (such as orbitally-resolved measurements) and to the generalization to out-of-equilibrium settings accessed in ultrafast settings.
We consider the problem of the development of steady states in one-dimensional Bose gas tubes that are weakly coupled to one another through a density-density interaction. We analyze this development through a Boltzmann collision integral approach. We argue that when the leading order of the collision integral, where single particle-hole excitations are created in individual gases, is dominant, the state of the gas evolves first to a non-thermal fixed point, i.e. a prethermalization plateau. This order is dominant when a pair of tubes are inequivalent with, say, different temperatures or different effective interaction parameters, \gammaγ. When both tubes are in the strongly interacting regime we additionally characterize this non-thermal prethermalization plateau by constructing the quasi-conserved quantities that control the existence of this plateau as well as the associated generalized Gibbs ensemble.
Efficient preparation of arbitrary entangled quantum states is crucial for quantum computation. This is particularly important for noisy intermediate-scale quantum simulators relying on variational hybrid quantum-classical algorithms. To that end, we propose symmetry-conserving modified quantum approximate optimization algorithm (SCom-QAOA) circuits. The depths of these circuits depend not only on the desired fidelity to the target state but also on the amount of entanglement the state contains. The parameters of the SCom-QAOA circuits are optimized using the quantum natural gradient method based on the Fubini-Study metric. The SCom-QAOA circuit transforms an unentangled state into a ground state of a gapped one-dimensional Hamiltonian with a circuit depth that depends not on the system size but rather on the finite correlation length. In contrast, the circuit depth grows proportionally to the system size for preparing low-lying states of critical one-dimensional systems. Even in the latter case, SCom-QAOA circuits with depth less than the system size were sufficient to generate states with fidelity in excess of 99%, which is relevant for near-term applications. The proposed scheme enlarges the set of the initial states accessible for variational quantum algorithms and widens the scope of investigation of nonequilibrium phenomena in quantum simulators.
We propose herein an extension of truncated spectrum methodologies (TSMs), a non-perturbative numerical approach able to elucidate the low energy properties of quantum field theories. TSMs, in their various flavors, involve a division of a computational Hilbert space, $\mathcal{H}$, into two parts, one part, $\mathcal{H}_1$ that is `kept' for the numerical computations, and one part, $\mathcal{H}_2$, that is discarded or `truncated'. Even though $\mathcal{H}_2$ is discarded, TSMs will often try to incorporate the effects of $\mathcal{H}_2$ in some effective way. In these terms, we propose to keep the dimension of $\mathcal{H}_1$ small. We pair this choice of $\mathcal{H}_1$ with a Krylov subspace iterative approach able to take into account the effects of $\mathcal{H}_2$. This iterative approach can be taken to arbitrarily high order and so offers the ability to compute quantities to arbitrary precision. In many cases it also offers the advantage of not needing an explicit UV cutoff. To compute the matrix elements that arise in the Krylov iterations, we employ a Feynman diagrammatic representation that is then evaluated with Monte Carlo techniques. Each order of the Krylov iteration is variational and is guaranteed to improve upon the previous iteration. The first Krylov iteration is akin to the NLO approach of Elias-Mir\'o et al. \cite{Elias-Miro:2017tup}. To demonstrate this approach, we focus on the $1+1d$ dimensional $\phi^4$ model and compute the bulk energy and mass gaps in both the $\mathbb{Z}_2$-broken and unbroken sectors. We estimate the critical $\phi^4$ coupling in the broken phase to be $g_c=0.2645\pm0.002$.
In recent quantum algorithmic developments, a feedback-based approach has shown promise for preparing quantum many-body system ground states and solving combinatorial optimization problems. This method utilizes quantum Lyapunov control to iteratively construct quantum circuits. Here, we propose a substantial enhancement by implementing a protocol that uses ideas from quantum Lyapunov control and the counterdiabatic driving protocol, a key concept from quantum adiabaticity. Our approach introduces an additional control field inspired by counterdiabatic driving. We apply our algorithm to prepare ground states in one-dimensional quantum Ising spin chains. Comprehensive simulations demonstrate a remarkable acceleration in population transfer to low-energy states within a significantly reduced time frame compared to conventional feedback-based quantum algorithms. This acceleration translates to a reduced quantum circuit depth, a critical metric for potential quantum computer implementation. We validate our algorithm on the IBM cloud computer, highlighting its efficacy in expediting quantum computations for many-body systems and combinatorial optimization problems.
We study the time-dependent response of a Heisenberg spin ladder subjected to a time-dependent square form variation of its rung spin exchange coupling. To do so, we employ a field theoretic representation of the Heisenberg spin ladder consisting of a singlet and a triplet of Majorana fermions. Because this underlying description is free fermionic, we are able to develop closed form analytic expressions for dynamical quantities, both one-body measures of local spin correlations as well as two-time correlation functions. These expressions involve both the gaps to triplet and singlet excitations. We analyze these expressions obtaining both the time scales for their transients and long-time athermal steady state behaviors. We show that variations in the rung coupling are directly tied to changes in the local antiferromagnetic correlations. We further discuss the application of these results to pump-probe experiments on material realizations of low-dimensional magnetic systems.
Many experimentally relevant systems are quasi-one-dimensional, consisting of nearly decoupled chains. In these systems, there is a natural separation of scales between the strong intrachain interactions and the weak interchain coupling. When the intrachain interactions are integrable, weak interchain couplings play a crucial part in thermalizing the system. Here, we develop a Boltzmann-equation formalism involving a collision integral that is asymptotically exact for any interacting integrable system, and apply it to develop a quantitative theory of relaxation in coupled Bose gases in the experimentally relevant Newton's cradle setup. We find that relaxation involves a broad spectrum of timescales. We provide evidence that the Markov process governing relaxation at late times is gapless; thus, the approach to equilibrium is generally nonexponential, even for spatially uniform perturbations.
In this study, we further the thermodynamic bootstrap program which involves a set of recently developed ideas used to determine thermodynamic form factors of local operators in integrable quantum field theories. These form factors are essential building blocks for dynamic correlation functions at finite temperatures or non-equilibrium stationary states. In this work we extend this program in three ways. Firstly, we demonstrate that the conjectured annihilation pole axiom is valid in the low energy particle-hole excitations. Secondly, we introduce a crossing relation, which establishes a connection between form factors with different excitation content. Typically, the crossing relation is a consequence of Lorentz invariance, but due to the finite energy density of the considered states, Lorentz invariance is broken. Nonetheless a crossing relation involving excitations with both particles and holes can established using the finite volume representation of the thermodynamic form factors. Finally, we demonstrate that the thermodynamic form factors satisfy a reparameterization invariance, an invariance which encompasses crossing. Reparameterization invariance exploits the fact that the details of the representation of the thermodynamic state are unimportant. In the course of developing these results, we demonstrate the internal consistency of the thermodynamic form factor bootstrap program in a number of ways. Finally, we provide explicit computations of form factors of conserved charges and densities with crossed excitations and show our results can be used to infer information about thermodynamic form factors in the Lieb-Liniger model.
We analyze the behavior of quarks coupled to a $SU(N_c)$ gauge theory in 1+1 dimensions. In the limit of strong coupling, the model reduces to a Wess-Zumino-Novikov-Witten (WZNW) model. At nonzero density, excitations near the Fermi surface form a non-Fermi liquid. With $N_f$ flavors, the finite density of quarks reduce to a free $U(1)$ field, which governs fluctuations in baryon number, together with a WZNW $SU(N_f)$ nonlinear sigma model at level $N_c$, from the pion/kaon modes. We compute the singularity in the charge susceptibility at the Fermi surface and the attendant power law correlations. We suggest that this is relevant to the quarkyonic regime of cold, dense QCD in 3+1 dimensions, in the limit that the Fermi surface is covered by many small patches, and the theory is effectively one dimensional. In this regime the dominant excitations near the Fermi surface are not baryons, but gapless bosonic modes.
from the drift mode having wide measurement span covering energy jitter range from 10 −4 to 10 −1 . The key advantage of this method is that it allows us to extract the electron beam energy spread concurrently with the ongoing experiment and enables online optimization of the electron beam especially for future high charge single-shot ultrafast electron diffraction (UED) and ultrafast electron microscopy (UEM) experiments. Furthermore, real-time energy measurement enables the filtering process to remove off-energy shots, improving the resolution of time-resolved UED. As a result, this method can be applied to the entire UED user community, beyond the traditional electron beam diagnostics of accelerators used by accelerator physicists.
A bstract We introduce and study generalized Rényi entropies defined through the traces of products of Tr B (| Ψ i ⟩⟨Ψ j | ) where ∣Ψ i ⟩ are eigenstates of a two-dimensional conformal field theory (CFT). When ∣Ψ i ⟩ = ∣Ψ j ⟩ these objects reduce to the standard Rényi entropies of the eigenstates of the CFT. Exploiting the path integral formalism, we show that the second generalized Rényi entropies are equivalent to four point correlators. We then focus on a free bosonic theory for which the mode expansion of the fields allows us to develop an efficient strategy to compute the second generalized Rényi entropy for all eigenstates. As a byproduct, our approach also leads to new results for the standard Rényi and relative entropies involving arbitrary descendent states of the bosonic CFT.