Search

Engineered dissipative reservoirs can steer many-body quantum systems toward correlated steady states useful for quantum simulation of high-temperature superconductivity or quantum magnetism. Using 49 superconducting qubits, we prepare low-energy states of the transverse-field Ising model via coupling to dissipative auxiliary qubits. In 1D, we observe long-range quantum correlations and a ground state fidelity that depends weakly on system sizes. In 2D, we find mutual information that extends beyond nearest neighbors. Lastly, by coupling the system to auxiliaries emulating reservoirs with different chemical potentials, we discover a new spin transport regime in the quantum Heisenberg model. Our results establish engineered dissipation as a scalable alternative to unitary evolution for preparing entangled many-body states on noisy quantum processors, and an essential tool for investigating nonequilibrium quantum phenomena.

Quantum computers hold the promise of executing tasks beyond the capability of classical computers. Noise competes with coherent evolution and destroys long-range correlations, making it an outstanding challenge to fully leverage the computation power of near-term quantum processors. We report Random Circuit Sampling (RCS) experiments where we identify distinct phases driven by the interplay between quantum dynamics and noise. Using cross-entropy benchmarking, we observe phase boundaries which can define the computational complexity of noisy quantum evolution. We conclude by presenting an RCS experiment with 70 qubits at 24 cycles. We estimate the computational cost against improved classical methods and demonstrate that our experiment is beyond the capabilities of existing classical supercomputers.

We examine the behavior of charged pions in neutron-rich matter using heavy-baryon chiral perturbation theory. This study is motivated by the prospect that pions, or pion-like, excitations, may be relevant in neutron-rich matter encountered in core-collapse supernovae and neutron star mergers. We find, as previously expected, that the $\pi^-$ mass increases with density and precludes s-wave condensation, and the mass of the $\pi^+$ mode decreases with increasing density. The energy difference between neutrons and protons in neutron-rich matter related to the nuclear symmetry energy alters the power counting. It enhances higher-order contributions to the pion self-energy. Previously unimportant but attractive diagrams are now enhanced. The net effect of this is the appearance of a new collective mode with the quantum numbers of the $\pi^+$, and a pronounced reduction of the $\pi^+$ mass., Comment: 13 pages, 4 figures, 1 appendix

The Monte Carlo evaluation of path integrals is one of a few general purpose methods to approach strongly coupled systems. It is used in all branches of Physics, from QCD/nuclear physics to the correlated electron systems. However, many systems of great importance (dense matter inside neutron stars, the repulsive Hubbard model away from half-filling, dynamical and non-equilibrium observables) are not amenable to the Monte Carlo method as it currently stands due to the so-called "sign-problem". We review a new set of ideas recently developed to tackle the sign problem based on the complexification of field space and the Picard-Lefshetz theory accompanying it. The mathematical ideas underpinning this approach, as well as the algorithms so far developed, are described together with non-trivial examples where the method has already been proved successful. Directions of future work, including the burgeoning use of machine learning techniques, are delineated., 35 pages, 21 figures

Monte Carlo simulations are useful tools for modeling quantum systems, but in some cases they suffer from a sign problem, leading to an exponential slow down in their convergence to a value. While solving the sign problem is generically NP-hard, many techniques exist for mitigating the sign problem in specific cases; in particular, the technique of deforming the Monte Carlo simulation's plane of integration onto Lefschetz thimbles (complex hypersurfaces of stationary phase) has seen significant success in the context of quantum field theories. We extend this methodology to spin systems by utilizing spin coherent state path integrals to re-express the spin system's partition function in terms of continuous variables. Using some toy systems, we demonstrate its effectiveness at lessening the sign problem in this setting, despite the fact that the initial mapping to spin coherent states introduces its own sign problem. The standard formulation of the spin coherent path integral is known to make use of uncontrolled approximations; despite this, for large spins they are typically considered to yield accurate results, so it is somewhat surprising that our results show significant systematic errors. Therefore, possibly of independent interest, our use of Lefschetz thimbles to overcome the intrinsic sign problem in spin coherent state path integral Monte Carlo enables a novel numerical demonstration of a breakdown in the spin coherent path integral., 11 pages, 5 figures

Path integral contour deformations have been shown to mitigate sign and signal-to-noise problems associated with phase fluctuations in lattice field theories. We define a family of contour deformations applicable to $SU(N)$ lattice gauge theory that can reduce sign and signal-to-noise problems associated with complex actions and complex observables. For observables, these contours can be used to define deformed observables with identical expectation value but different variance. As a proof-of-principle, we apply machine learning techniques to optimize the deformed observables associated with Wilson loops in two dimensional $SU(2)$ and $SU(3)$ gauge theory. We study loops consisting of up to 64 plaquettes and achieve variance reduction of up to 4 orders of magnitude., Comment: 25 pages, 12 figures

Monte Carlo studies of many quantum systems face exponentially severe signal-to-noise problems. We show that noise arising from complex phase fluctuations of observables can be reduced without introducing bias using path integral contour deformation techniques. A numerical study of contour deformations for correlation functions in Abelian gauge theory and complex scalar field theory demonstrates that variance can be reduced by orders of magnitude without modifying Monte Carlo sampling., 7 pages, 2 figures

We compute continuum and infinite volume limit extrapolations of the structure factors of neutron matter at finite temperature and density. Using a lattice formulation of leading-order pionless effective field theory, we compute the momentum dependence of the structure factors at finite temperature and at densities beyond the reach of the virial expansion. The Tan contact parameter is computed and the result agrees with the high momentum tail of the vector structure factor. All errors, statistical and systematic, are controlled for. This calculation is a first step towards a model-independent understanding of the linear response of neutron matter at finite temperature, a realm until now little explored., Comment: 5 pages, 3 figures, data included in submission; further experimental data included

We introduce a theoretical framework to analyze the thermal properties of a recently proposed "excluded volume" model of quarkyonic matter. This entails proposing finite temperature distribution functions and entropy functionals constrained by internal theoretical consistency. We use the formalism developed here to analyze the effect of temperature on the quark onset as a function of baryon density. We find that the baryon density at which the quarkyonic phase emerges decreases only slightly with temperature, an effect produced by large $N_c$. We furthermore find that a finite temperature increases the density of quarks produced. The framework presented here can be used to compute the equation of state of hot and dense quarkyonic matter which is of relevance to neutron star mergers., Comment: 14 pages, 6 figures

The interplay of interactions and strong disorder can lead to an exotic quantum many-body localized (MBL) phase. Beyond the absence of transport, the MBL phase has distinctive signatures, such as slow dephasing and logarithmic entanglement growth; they commonly result in slow and subtle modification of the dynamics, making their measurement challenging. Here, we experimentally characterize these properties of the MBL phase in a system of coupled superconducting qubits. By implementing phase sensitive techniques, we map out the structure of local integrals of motion in the MBL phase. Tomographic reconstruction of single and two qubit density matrices allowed us to determine the spatial and temporal entanglement growth between the localized sites. In addition, we study the preservation of entanglement in the MBL phase. The interferometric protocols implemented here measure affirmative correlations and allow us to exclude artifacts due to the imperfect isolation of the system. By measuring elusive MBL quantities, our work highlights the advantages of phase sensitive measurements in studying novel phases of matter., 5+28 pages, 5+22 figures, updated version

Simulations of gauge theories on quantum computers require the digitization of continuous field variables. Digitization schemes that uses the minimum amount of qubits are desirable. We present a practical scheme for digitizing $SU(3)$ gauge theories via its discrete subgroup $S(1080)$. The $S(1080)$ standard Wilson action cannot be used since a phase transition occurs as the coupling is decreased, well before the scaling regime. We proposed a modified action that allows simulations in the scaling window and carry out classical Monte Carlo calculations down to lattice spacings of order $a\approx 0.08$ fm. We compute a set of observables with sub-percent precision at multiple lattice spacings and show that the continuum extrapolated value agrees with the full $SU(3)$ results. This suggests that this digitization scheme provides sufficient precision for NISQ-era QCD simulations., Comment: 5 pages, 5 figures

We compute with lattice field theory the vector and axial static structure factors of the unitary gas for arbitrary temperature above the superfluid transition and for fugacities 0.1 < z < 1.0. Using the lattice formulation, we calculate beyond the validity of the virial expansion, a commonly used technique in many-body physics. We find qualitative differences in the behavior of the structure factors at high fugacity compared to the predictions of the virial expansion. Due to the large scattering length of neutrons, we expect the unitary gas structure factors to approximate the structure factors of hot neutron gases, and we therefore expect our calculations to be useful in supernova simulations, where neutron gas structure factors are needed to compute in-medium neutrino-neutron scattering rates., 10 pages, 6 figures, 1 table

We present Monte Carlo calculations of the thermodynamics of the (2+1) dimensional Thirring model at finite density. We bypass the sign problem by deforming the domain of integration of the path integral into complex space in such a way as to maximize the average sign within a parameterized family of manifolds. We present results for lattice sizes up to $10^3$ and we find that at high densities and/or temperatures the chiral condensate is abruptly reduced., 5 pages, 4 figures

We study by Monte Carlo methods the thermodynamics of a spin polarized gas of non-relativistic fermions in 1+1 dimensions. The main result of this work is that our action suffers no significant sign problem for any spin polarization in the region relevant for dilute degenerate fermi gases. This lack of sign problem allows us to study attractive spin polarized fermions non-perturbatively at spin polarizations not previously explored. For some parameters values we verify results previously obtained by methods which include an uncontrolled step like complex Langevin and/or analytical continuation from imaginary chemical potential. For others, larger values of the polarization, we deviate from these previous results., 9 pages, 5 figures

We calculate the differential scattering rate for thermal neutrinos in a hot and dilute gas of interacting neutrons using linear response theory. The dynamical structure factors for density and spin fluctuations of the strongly interacting neutron matter, expected in the neutrino decoupling regions of supernovae and neutron star mergers, are calculated in the virial expansion for the first time. Correlations due to nucleon-nucleon interactions are taken into account using a pseudo-potential that reproduces measured nucleon-nucleon phase shifts, and we find that attractive s-wave interactions enhance the density response and suppress the spin response of neutron matter. The net effect of neutron correlations is to strongly suppress backscattering. Moreover, we find nearly exact scaling laws for the response functions, valid for the range $T = 5 - 10$ MeV and q < 30 MeV, allowing us to obtain analytic results for the dynamic structure factors at second-order in the fugacity of the neutron gas. We find that the modification of scattering rates depends on the energy and momentum exchanged, implying that dynamical structure factors are essential to describe neutrino decoupling in supernovae and neutron star mergers.

Complex integrated circuits require multiple wiring layers. In complementary metal-oxide-semiconductor (CMOS) processing, these layers are robustly separated by amorphous dielectrics. These dielectrics would dominate energy loss in superconducting integrated circuits. Here we demonstrate a procedure that capitalizes on the structural benefits of inter-layer dielectrics during fabrication and mitigates the added loss. We separate and support multiple wiring layers throughout fabrication using SiO$_2$ scaffolding, then remove it post-fabrication. This technique is compatible with foundry level processing and the can be generalized to make many different forms of low-loss multi-layer wiring. We use this technique to create freestanding aluminum vacuum gap crossovers (airbridges). We characterize the added capacitive loss of these airbridges by connecting ground planes over microwave frequency $\lambda/4$ coplanar waveguide resonators and measuring resonator loss. We measure a low power resonator loss of $\sim 3.9 \times 10^{-8}$ per bridge, which is 100 times lower than dielectric supported bridges. We further characterize these airbridges as crossovers, control line jumpers, and as part of a coupling network in gmon and fuxmon qubits. We measure qubit characteristic lifetimes ($T_1$'s) in excess of 30 $\mu$s in gmon devices.

Quantum field theories with complex actions cannot be investigated using importance sampling due to the sign problem. One possible solution is to use the holomorphic gradient flow, a method we introduced related to the Lefschetz thimbles idea. In many cases the probability distribution generated by this method is multi-modal and standard Monte-Carlo sampling fails. We propose an algorithm that incorporates tempered proposals to solve this problem. We apply this algorithm to the 0+1 dimensional Thirring model at finite density for a parameter set where standard sampling fails and show that tempered proposals cure this problem., 10 pages, 6 figures

Statistical mechanics is founded on the assumption that a system can reach thermal equilibrium, regardless of the starting state. Interactions between particles facilitate thermalization, but, can interacting systems always equilibrate regardless of parameter values\,? The energy spectrum of a system can answer this question and reveal the nature of the underlying phases. However, most experimental techniques only indirectly probe the many-body energy spectrum. Using a chain of nine superconducting qubits, we implement a novel technique for directly resolving the energy levels of interacting photons. We benchmark this method by capturing the intricate energy spectrum predicted for 2D electrons in a magnetic field, the Hofstadter butterfly. By increasing disorder, the spatial extent of energy eigenstates at the edge of the energy band shrink, suggesting the formation of a mobility edge. At strong disorder, the energy levels cease to repel one another and their statistics approaches a Poisson distribution - the hallmark of transition from the thermalized to the many-body localized phase. Our work introduces a new many-body spectroscopy technique to study quantum phases of matter.

A relativistic Bose gas at finite density suffers from a sign problem that makes direct numerical simulations not feasible. One possible solution to the sign problem is to re-express the path integral in terms of Lefschetz thimbles. Using this approach we study the relativistic Bose gas both in the symmetric phase (low-density) and the spontaneously broken phase (high-density). In the high-density phase we break explicitly the symmetry and determine the dependence of the order parameter on the breaking. We study the relative contributions of the dominant and sub-dominant thimbles in this phase. We find that the sub-dominant thimble only contributes substantially when the explicit symmetry breaking is small, a regime that is dominated by finite volume effects. In the regime relevant for the thermodynamic limit, this contribution is negligible., 12 pages, 6 figures, 1 table

A solution to the sign problem is the so-called "Lefschetz thimble approach" where the domain of integration for field variables in the path integral is deformed from the real axis to a sub-manifold in the complex space. For properly chosen sub-manifolds ("thimbles") the sign problem disappears or is drastically alleviated. The parametrization of the thimble by real coordinates require the calculation of a jacobian with a computational cost of order O(V^3), where V is proportional to the spacetime volume. In this note we propose two estimators for this jacobian with a computational cost of order O(V). We discuss analytically the regimes where we expect the estimator to work and show numerical examples in two different models., 10 pages, 3 figures

We present a study of quasar selection using the DES supernova fields. We used a quasar catalog from an overlapping portion of the SDSS Stripe 82 region to quantify the completeness and efficiency of selection methods involving color, probabilistic modeling, variability, and combinations of color/probabilistic modeling with variability. We only considered objects that appear as point sources in the DES images. We examine color selection methods based on the WISE mid-IR W1-W2 color, a mixture of WISE and DES colors (g-i and i-W1) and a mixture of VHS and DES colors (g-i and i-K). For probabilistic quasar selection, we used XDQSOz, an algorithm that employs an empirical multi-wavelength flux model of quasars to assign quasar probabilities. Our variability selection uses the multi-band chi2-probability that sources are constant in the DES Year 1 griz-band light curves. The completeness and efficiency are calculated relative to an underlying sample of point sources that are detected in the required selection bands and pass our data quality and photometric error cuts. We conduct our analyses at two magnitude limits, i85% for both i-band magnitude limits and efficiencies of >80% to the bright limit and >60% to the faint limit; however, the giW1 and giW1+variability methods give the highest quasar surface densities. The XDQSOz method and combinations of W1W2/giW1/XDQSOz with variability are among the better selection methods when both high completeness and high efficiency are desired. We also present the OzDES Quasar Catalog of 1,263 spectroscopically-confirmed quasars taken by the OzDES survey. The catalog includes quasars with redshifts up to z~4 and brighter than i=22 mag, although the catalog is not complete up this magnitude limit., 16 pages, 12 figures, 3 tables; Submitted to ApJ. The OzDES Quasar Catalog can be downloaded at http://bit.ly/2fWsK9t . For a video summary of the paper, please see https://www.youtube.com/watch?v=JeQNdUATze0

We present results of the numerical simulation of the two-dimensional Thirring model at finite density and temperature. The severe sign problem is dealt with by deforming the domain of integration into complex field space. This is the first example where a fermionic sign problem is solved in a quantum field theory by using the holomorphic gradient flow approach, a generalization of the Lefschetz thimble method., Comment: 6 pages, 5 figures

Monte Carlo studies involving real time dynamics are severely restricted by the sign problem that emerges from highly oscillatory phase of the path integral. In this letter, we present a new method to compute real time quantities on the lattice using the Schwinger-Keldysh formalism via Monte Carlo simulations. The key idea is to deform the path integration domain to a complex manifold where the phase oscillations are mild and the sign problem is manageable. We use the previously introduced "contraction algorithm" to create a Markov chain on this alternative manifold. We substantiate our approach by analyzing the quantum mechanical anharmonic oscillator. Our results are in agreement with the exact ones obtained by diagonalization of the Hamiltonian. The method we introduce is generic and in principle applicable to quantum field theory albeit very slow. We discuss some possible improvements that should speed up the algorithm., Comment: 5 pages, 3 figures

We point out that Monte Carlo simulations of theories with severe sign problems can be profitably performed over manifolds in complex space different from the one with fixed imaginary part of the action. We describe a family of such manifolds that interpolate between the tangent space at one critical point, where the sign problem is milder compared to the real plane but in some cases still severe, and the union of relevant thimbles, where the sign problem is mild but a multimodal distribution function complicates the Monte Carlo sampling. We exemplify this approach using a simple 0 + 1 dimensional fermion model previously used on sign problem studies and show that it can solve the model for some parameter values where a solution using Lefshetz thimbles was elusive., Comment: 12 pages, 7 figures, figure added, corrected typos

Weak measurement has provided new insight into the nature of quantum measurement, by demonstrating the ability to extract average state information without fully projecting the system. For single qubit measurements, this partial projection has been demonstrated with violations of the Leggett-Garg inequality. Here we investigate the effects of weak measurement on a maximally entangled Bell state through application of the Hybrid Bell-Leggett-Garg inequality (BLGI) on a linear chain of four transmon qubits. By correlating the results of weak ancilla measurements with subsequent projective readout, we achieve a violation of the BLGI with 27 standard deviations of certainty.

A quantum computer can solve hard problems - such as prime factoring, database searching, and quantum simulation - at the cost of needing to protect fragile quantum states from error. Quantum error correction provides this protection, by distributing a logical state among many physical qubits via quantum entanglement. Superconductivity is an appealing platform, as it allows for constructing large quantum circuits, and is compatible with microfabrication. For superconducting qubits the surface code is a natural choice for error correction, as it uses only nearest-neighbour coupling and rapidly-cycled entangling gates. The gate fidelity requirements are modest: The per-step fidelity threshold is only about 99%. Here, we demonstrate a universal set of logic gates in a superconducting multi-qubit processor, achieving an average single-qubit gate fidelity of 99.92% and a two-qubit gate fidelity up to 99.4%. This places Josephson quantum computing at the fault-tolerant threshold for surface code error correction. Our quantum processor is a first step towards the surface code, using five qubits arranged in a linear array with nearest-neighbour coupling. As a further demonstration, we construct a five-qubit Greenberger-Horne-Zeilinger (GHZ) state using the complete circuit and full set of gates. The results demonstrate that Josephson quantum computing is a high-fidelity technology, with a clear path to scaling up to large-scale, fault-tolerant quantum circuits., 15 pages, 13 figures, including supplementary material

Understanding complex quantum matter presents a central challenge in condensed matter physics. The difficulty lies in the exponential scaling of the Hilbert space with the system size, making solutions intractable for both analytical and conventional numerical methods. As originally envisioned by Richard Feynman, this class of problems can be tackled using controllable quantum simulators. Despite many efforts, building an quantum emulator capable of solving generic quantum problems remains an outstanding challenge, as this involves controlling a large number of quantum elements. Here, employing a multi-element superconducting quantum circuit and manipulating a single microwave photon, we demonstrate that we can simulate the weak localization phenomenon observed in mesoscopic systems. By engineering the control sequence in our emulator circuit, we are also able to reproduce the well-known temperature dependence of weak localization. Furthermore, we can use our circuit to continuously tune the level of disorder, a parameter that is not readily accessible in mesoscopic systems. By demonstrating a high level of control and complexity, our experiment shows the potential for superconducting quantum circuits to realize scalable quantum simulators., Comment: 9 pages, including supplement

Progress in superconducting qubit experiments with greater numbers of qubits or advanced techniques such as feedback requires faster and more accurate state measurement. We have designed a multiplexed measurement system with a bandpass filter that allows fast measurement without increasing environmental damping of the qubits. We use this to demonstrate simultaneous measurement of four qubits on a single superconducting integrated circuit, the fastest of which can be measured to 99.8% accuracy in 140ns. This accuracy and speed is suitable for advanced multi-qubit experiments including surface code error correction., Five figures

We demonstrate a high efficiency deterministic quantum receiver to convert flying qubits to logic qubits. We employ a superconducting resonator, which is driven with a shaped pulse through an adjustable coupler. For the ideal "time reversed" shape, we measure absorption and receiver fidelities at the single microwave photon level of, respectively, 99.41% and 97.4%. These fidelities are comparable with gates and measurement and exceed the deterministic quantum communication and computation fault tolerant thresholds., Main paper: 5 pages, 4 figures. Supplement: 11 pages, 12 figures. Revised abstract and introduction. Minor changes to Figure 1 and figure captions

The quantum behavior of superconducting qubits coupled to resonators is very similar to that of atoms in optical cavities [1, 2], in which the resonant cavity confines photons and promotes strong light-matter interactions. The cavity end-mirrors determine the performance of the coupled system, with higher mirror reflectivity yielding better quantum coherence, but higher mirror transparency giving improved measurement and control, forcing a compromise. An alternative is to control the mirror transparency, enabling switching between long photon lifetime during quantum interactions and large signal strength when performing measurements. Here we demonstrate the superconducting analogue, using a quantum system comprising a resonator and a qubit, with variable coupling to a measurement transmission line. The coupling can be adjusted through zero to a photon emission rate 1,000 times the intrinsic photon decay rate. We use this system to control photons in coherent states as well as in non-classical Fock states, and dynamically shape the waveform of released photons. This has direct applications to circuit quantum electrodynamics [3], and may enable high-fidelity quantum state transfer between distant qubits, for which precisely-controlled waveform shaping is a critical and non-trivial requirement [4, 5].

We present photometric redshift estimates for galaxies used in the weak lensing analysis of the Dark Energy Survey Science Verification (DES SV) data. Four model- or machine learning-based photometric redshift methods -- ANNZ2, BPZ calibrated against BCC-Ufig simulations, SkyNet, and TPZ -- are analysed. For training, calibration, and testing of these methods, we construct a catalogue of spectroscopically confirmed galaxies matched against DES SV data. The performance of the methods is evaluated against the matched spectroscopic catalogue, focusing on metrics relevant for weak lensing analyses, with additional validation against COSMOS photo-zs. From the galaxies in the DES SV shear catalogue, which have mean redshift $0.72\pm0.01$ over the range $0.3, Comment: high-resolution versions of figures can be downloaded from http://deswl.github.io