A long time ago, Bloch showed that in a system of interacting nonrelativistic particles the net particle-number current must vanish in any equilibrium state. Bloch’s argument does not generalize easily to the energy current. We devise an alternative argument which proves the vanishing of the net energy currents in equilibrium states of lattice systems as well as systems of nonrelativistic particles with finite-range potential interactions. We discuss some applications of these results. In particular, we show that neither a one-dimensional (1D) lattice system nor a 1D system of nonrelativistic particles with finite-range potential interactions can flow to a conformal field theory with unequal left-moving and right-moving central charges.

We consider simulating an $n$-qubit Hamiltonian with nearest-neighbor interactions evolving for time $t$ on a quantum computer. We show that this simulation has gate complexity $(nt{)}^{1+o(1)}$ using product formulas, a straightforward approach that has been demonstrated by several experimental groups. While it is reasonable to expect this complexity—in particular, this was claimed without rigorous justification by Jordan, Lee, and Preskill—we are not aware of a straightforward proof. Our approach is based on an analysis of the local error structure of product formulas, as introduced by Descombes and Thalhammer and significantly simplified here. We prove error bounds for canonical product formulas, which include well-known constructions such as the Lie-Trotter-Suzuki formulas. We also develop a local error representation for time-dependent Hamiltonian simulation, and we discuss generalizations to periodic boundary conditions, constant-range interactions, and higher dimensions. Combined with a previous lower bound, our result implies that product formulas can simulate lattice Hamiltonians with nearly optimal gate complexity.

Correlations of fluctuations are the driving forces behind the dynamics and thermodynamics in quantum many-body systems. For qubits embedded in a quantum bath, the correlations in the bath are key to understanding and combating decoherence—a critical issue in quantum information technology. However, there is no systematic method for characterizing the many-body correlations in quantum baths beyond the second order or the Gaussian approximation. Here we present a scheme to characterize the correlations in a quantum bath to arbitrary order. The scheme employs a weak measurement of the bath via the projective measurement of a central system. The bath correlations, including both the “classical” and the “quantum” parts, can be reconstructed from the correlations of the measurement outputs. The possibility of full characterization of many-body correlations in a quantum bath forms the basis for optimizing quantum control against decoherence in realistic environments, for studying the quantum characteristics of baths, and for the quantum sensing of correlated clusters in quantum baths.

Extendibility of bosonic Gaussian states is a key issue in continuous-variable quantum information. We show that a bosonic Gaussian state is $k$-extendible if and only if it has a Gaussian $k$-extension, and we derive a simple semidefinite program, whose size scales linearly with the number of local modes, to efficiently decide $k$-extendibility of any given bosonic Gaussian state. When the system to be extended comprises one mode only, we provide a closed-form solution. Implications of these results for the steerability of quantum states and for the extendibility of bosonic Gaussian channels are discussed. We then derive upper bounds on the distance of a $k$-extendible bosonic Gaussian state to the set of all separable states, in terms of trace norm and Rényi relative entropies. These bounds, which can be seen as “Gaussian de Finetti theorems,” exhibit a universal scaling in the total number of modes, independently of the mean energy of the state. Finally, we establish an upper bound on the entanglement of formation of Gaussian $k$-extendible states, which has no analogue in the finite-dimensional setting.

Quantum Brownian motion plays a fundamental role in many areas of modern physics. In the path-integral formulation, environmental fluctuations can be characterized by auxiliary stochastic fields. Intriguingly, for fermionic environments the stochastic fields must be Grassmann valued so as to memorize the order of the random forces exerted on the system. Such nonclassical fields cannot be represented by conventional means. We propose a strategy to map the Grassmann-number fields to conventional $c$-number noises and a set of quantized pseudolevels. The resulting stochastic equation of motion (SEOM) enables direct stochastic simulation of the fermionic dissipative dynamics. The SEOM gives exact physical observables of noninteracting systems, and yields accurate approximate results for interacting systems. The practicality and accuracy of the proposed strategy and the SEOM are exemplified by numerical studies conducted on a single-impurity Anderson model.

We find the asymptotic representation of the solution of the variable-order fractional diffusion equation, which remains unsolved since it was proposed by Chechkin, Gorenflo, and Sokolov [J. Phys. A, 38, L679 (2005)]. We identify a new advection term that causes ultraslow spatial aggregation of subdiffusive particles due to dominance over the standard advection and diffusion terms in the long-time limit. This uncovers the anomalous mechanism by which nonuniform distributions can occur. We perform Monte Carlo simulations of the underlying anomalous random walk and find excellent agreement with the asymptotic solution.

Superconducting circuits have emerged as a powerful platform of quantum simulation, especially for emulating the dynamics of quantum many-body systems, because of their tunable interaction, long coherence time, and high-precision control. Here in experiments, we construct a Bose-Hubbard ladder with a ladder array of 20 qubits on a 24-qubit superconducting processor. We investigate theoretically and demonstrate experimentally the dynamics of single- and double-excitation states with distinct behaviors, indicating the uniqueness of the Bose-Hubbard ladder. We observe the linear propagation of photons in the single-excitation case, satisfying the Lieb-Robinson bounds. The double-excitation state, initially placed at the edge, localizes; while placed in the bulk, it splits into two single-excitation modes spreading linearly toward two boundaries, respectively. Remarkably, these phenomena, studied both theoretically and numerically as unique properties of the Bose-Hubbard ladder, are represented coherently by pairs of controllable qubits in experiments. Our results show that collective excitations, as a single mode, are not free. This work paves the way to simulation of exotic logic particles by subtly encoding physical qubits and exploration of rich physics by superconducting circuits.

We report on the realization and application of nondestructive three-qubit parity measurements on nuclear spin qubits in diamond. We use high-fidelity quantum logic to map the parity of the joint state of three nuclear spin qubits onto an electronic spin qubit that acts as an ancilla, followed by a single-shot nondestructive readout of the ancilla combined with an electron spin echo to ensure outcome-independent evolution of the nuclear spins. Through the sequential application of three such parity measurements, we demonstrate the generation of genuine multipartite entangled states out of the maximally mixed state. Furthermore, we implement a single-shot version of the Greenberger-Horne-Zeilinger experiment that can generate a quantum versus classical contradiction in each run. Finally, we test a state-independent noncontextuality inequality in eight dimensions. The techniques and insights developed are relevant for fundamental tests as well as for quantum information protocols such as quantum error correction.

The advantage of quantum metrology has been experimentally demonstrated for phase estimations where the dynamics are commuting. General noncommuting dynamics, however, can have distinct features. For example, the direct sequential scheme, which can achieve the Heisenberg scaling for the phase estimation under commuting dynamics, can have even worse performances than the classical scheme when the dynamics are noncommuting. Here we realize a scalable optimally controlled sequential scheme, which can achieve the Heisenberg precision under general noncommuting dynamics. We also present an intuitive geometrical framework for the controlled scheme and identify sweet spots in time at which the optimal controls used in the scheme can be prefixed without adaptation, which simplifies the experimental protocols significantly. We successfully implement the scheme up to eight controls in an optical platform and demonstrate a precision near the Heisenberg limit. Our work opens the avenue for harvesting the power of quantum control in quantum metrology, and provides a control-enhanced recipe to achieve the Heisenberg precision under general noncommuting dynamics.

Standard computation of size and credibility of a Bayesian credible region for certifying any point estimator of an unknown parameter (such as a quantum state, channel, phase, etc.) requires selecting points that are in the region from a finite parameter-space sample, which is infeasible for a large dataset or dimension as the region would then be extremely small. We solve this problem by introducing the in-region sampling theory to compute both region qualities just by sampling appropriate functions over the region itself using any Monte Carlo sampling method. We take in-region sampling to the next level by understanding the credible-region capacity (an alternative description for the region content to size) as the average $l_p$-norm distance $(p>0)$ between a random region point and the estimator, and present analytical formulas for $p=2$ to estimate both the capacity and credibility for any dimension and a sufficiently large dataset without Monte Carlo sampling, thereby providing a quick alternative to Bayesian certification. All results are discussed in the context of quantum-state tomography.

Originally defined for the optimal allocation of resources, optimal transport (OT) has found many theoretical and practical applications in multiple domains of science and physics. In this Letter we develop a new method for solving the discrete version of this problem using techniques derived from statistical physics. We derive a strongly concave free energy function that captures the constraints of the OT problem at a finite temperature. Its maximum defines an optimal transport plan, or registration between the two discrete probability measures that are compared, as well as a pseudodistance between those measures that satisfies the triangular inequalities. The computation of this pseudodistance is fast and numerically stable. The temperature dependent OT pseudodistance is shown to decrease monotonically with respect to the inverse of the temperature and to converge to the standard OT distance at zero temperature, providing a robust framework for temperature annealing. We illustrate applications of this framework to the problem of image comparison.

Quantum metrology employs quantum effects to attain a measurement precision surpassing the limit achievable in classical physics. However, it was previously found that the precision returns the shot-noise limit (SNL) from the ideal Zeno limit (ZL) due to the photon loss in quantum metrology based on Mech-Zehnder interferometry. Here, we find that not only can the SNL be beaten, but also the ZL can be asymptotically recovered in a long-encoding-time condition when the photon dissipation is exactly studied in its inherent non-Markovian manner. Our analysis reveals that it is due to the formation of a bound state of the photonic system and its dissipative noise. Highlighting the microscopic mechanism of the dissipative noise on the quantum optical metrology, our result supplies a guideline to realize the ultrasensitive measurement in practice by forming the bound state in the setting of reservoir engineering.

In the classical domain, it is well known that divisibility does not imply that a stochastic process is Markovian. However, for quantum processes, divisibility is often considered to be synonymous with Markovianity. We show that completely positive divisible quantum processes can still involve non-Markovian temporal correlations, that we then fully classify using the recently developed process tensor formalism, which generalizes the theory of stochastic processes to the quantum domain.

We consider the quench dynamics of a two-dimensional quantum dimer model and determine the role of its kinematic constraints. We interpret the nonequilibrium dynamics in terms of the underlying equilibrium phase transitions consisting of a Berezinskii-Kosterlitz-Thouless (BKT) transition between a columnar ordered valence bond solid (VBS) and a valence bond liquid (VBL), as well as a first-order transition between a staggered VBS and the VBL. We find that quenches from a columnar VBS are ergodic and both order parameters and spatial correlations quickly relax to their thermal equilibrium. By contrast, the staggered side of the first-order transition does not display thermalization on numerically accessible timescales. Based on the model’s kinematic constraints, we uncover a mechanism of relaxation that rests on emergent, highly detuned multidefect processes in a staggered background, which gives rise to slow, glassy dynamics at low temperatures even in the thermodynamic limit.

Benchmarking methods that can be adapted to multiqubit systems are essential for assessing the overall or “holistic” performance of nascent quantum processors. The current industry standard is Clifford randomized benchmarking (RB), which measures a single error rate that quantifies overall performance. But, scaling Clifford RB to many qubits is surprisingly hard. It has only been performed on one, two, and three qubits as of this writing. This reflects a fundamental inefficiency in Clifford RB: the $n$-qubit Clifford gates at its core have to be compiled into large circuits over the one- and two-qubit gates native to a device. As $n$ grows, the quality of these Clifford gates quickly degrades, making Clifford RB impractical at relatively low $n$. In this Letter, we propose a direct RB protocol that mostly avoids compiling. Instead, it uses random circuits over the native gates in a device, which are seeded by an initial layer of Clifford-like randomization. We demonstrate this protocol experimentally on two to five qubits using the publicly available ibmqx5. We believe this to be the greatest number of qubits holistically benchmarked, and this was achieved on a freely available device without any special tuning up. Our protocol retains the simplicity and convenient properties of Clifford RB: it estimates an error rate from an exponential decay. But, it can be extended to processors with more qubits—we present simulations on 10+ qubits—and it reports a more directly informative and flexible error rate than the one reported by Clifford RB. We show how to use this flexibility to measure separate error rates for distinct sets of gates, and we use this method to estimate the average error rate of a set of cnot gates.

We investigate the relaxation dynamics of an interacting Stark-localized system coupled to a dephasing bath, and compare its behavior to the conventional disorder-induced many body localized system. Specifically, we study the dynamics of population imbalance between even and odd sites, and the growth of the von Neumann entropy. For a large potential gradient, the imbalance is found to decay on a timescale $\tau$ that grows quadratically with the Wannier-Stark tilt. For the noninteracting system, it shows an exponential decay, which becomes a stretched exponential decay in the presence of finite interactions. This is different from a system with disorder-induced localization, where the imbalance exhibits a stretched exponential decay also for vanishing interactions. As another clear qualitative difference, we do not find a logarithmically slow growth of the von Neumann entropy as it is found for the disordered system. Our findings can immediately be tested experimentally with ultracold atoms in optical lattices.

We show how, upon heating the spin degrees of freedom of the Hubbard model to infinite temperature, the symmetries of the system allow the creation of steady states with long-range correlations between $\eta$ pairs. We induce this heating with either dissipation or periodic driving and evolve the system towards a nonequilibrium steady state, a process which melts all spin order in the system. The steady state is identical in both cases and displays distance-invariant off-diagonal $\eta$ correlations. These correlations were first recognized in the superconducting eigenstates described in Yang’s seminal Letter [Phys. Rev. Lett. 63, 2144 (1989)], which are a subset of our steady states. We show that our results are a consequence of symmetry properties and entirely independent of the microscopic details of the model and the heating mechanism.

The generation and distribution of entanglement are key resources in quantum repeater schemes. Temporally multiplexed systems offer time-bin encoding of quantum information which provides robustness against decoherence in fibers, crucial in long distance communication. Here, we demonstrate the direct generation of entanglement in time between a photon and a collective spin excitation in a rare earth ion doped ensemble. We analyze the entanglement by mapping the atomic excitation onto a photonic qubit and by using time-bin qubit analyzers implemented with another doped crystal using the atomic frequency comb technique. Our results provide a solid-state source of entangled photons with embedded quantum memory. Moreover, the quality of the entanglement is high enough to enable a violation of a Bell inequality by more than two standard deviations.

Quantum pseudorandomness, also known as unitary designs, comprises a powerful resource for emergent quantum technologies. Although in theory pseudorandom unitary operators can be constructed efficiently, realizing these objects in realistic physical systems is a challenging task. Here, we demonstrate experimental generation and detection of quantum pseudorandomness on a 12-qubit nuclear magnetic resonance system. We first apply random sequences to the interacting nuclear spins, leading to random quantum evolutions that can quickly form unitary designs. Then, in order to probe the growth of quantum pseudorandomness during the time evolutions, we propose the idea of using the system’s multiple-quantum coherence distribution as an indicator. Based on this indicator, we measure the spreading of quantum coherences and find that substantial quantum pseudorandomness has been achieved at the 12-qubit scale. This may open up a path to experimentally explore quantum randomness on forthcoming large-scale quantum processors.