Many-body Research Articles

Wave-Function Network Description and Kolmogorov Complexity of Quantum Many-Body Systems

Programmable quantum devices are now able to probe wave functions at unprecedented levels. This is based on the ability to project the many-body state of atom and qubit arrays onto a measurement basis which produces snapshots of the system wave function. Extracting and processing information from such observations remains, however, an open quest. One often resorts to analyzing low-order correlation functions—that is, discarding most of the available information content. Here, we introduce wave-function networks—a mathematical framework to describe wave-function snapshots based on network theory. For many-body systems, these networks can become scale-free—a mathematical structure that has found tremendous success and applications in a broad set of fields, ranging from biology to epidemics to Internet science. We demonstrate the potential of applying these techniques to quantum science by introducing protocols to extract the Kolmogorov complexity corresponding to the output of a quantum simulator and implementing tools for fully scalable cross-platform certification based on similarity tests between networks. We demonstrate the emergence of scale-free networks analyzing experimental data obtained with a Rydberg quantum simulator manipulating up to 100 atoms. Our approach illustrates how, upon crossing a phase transition, the simulator complexity decreases while correlation length increases—a direct signature of buildup of universal behavior in data space. Comparing experiments with numerical simulations, we achieve cross-certification at the wave-function level up to timescales of 4 μs with a confidence level of 90% and determine experimental calibration intervals with unprecedented accuracy. Our framework is generically applicable to the output of quantum computers and simulators with access to the system wave function and requires probing accuracy and repetition rates accessible to most currently available platforms. Published by the American Physical Society 2024

Pseudo-Spin Symmetry and the Hints for Unstable and Superheavy Nuclei

The pseudo-spin symmetry (PSS) provides an important angle to understand nuclear microscopic structure and the novel phenomena found in unstable nuclei. The relativistic Hartree–Fock (RHF) theory, that takes the important degrees of freedom associated with the π-meson and ρ-tensor (ρ-T) couplings into account, provides an appropriate description of the PSS restoration in realistic nuclei, particularly for the pseudo-spin (PS) doublets with high angular momenta (l˜). The investigations of the PSS within the RHF theory are recalled in this paper by focusing on the effects of the Fock terms. Aiming at common artificial shell closures appearing in previous relativistic mean-field calculations, the mechanism responsible for the PSS restoration of high-l˜ orbits is stressed, revealing the manifestation of nuclear in-medium effects on the PSS, and thus, providing qualitative guidance on modeling the in-medium balance between nuclear attractions and repulsions. Moreover, the essential role played by the ρ-T coupling, that contributes mainly via the Fock terms, is introduced as combined with the relations between the PSS and various nuclear phenomena, including the shell structure and the evolution, novel halo and bubble-like phenomena, and the superheavy magicity. As the consequences of the nuclear force in complicated nuclear many-body systems, the PSS itself and the mechanism therein can not only deepen our understanding of nuclear microscopic structure and relevant phenomena, but also provide special insight into the nature of the nuclear force, which can further enrich our knowledge of nuclear physics.

Efficiency of neural-network state representations of one-dimensional quantum spin systems

Neural-network state representations of quantum many-body systems are attracting great attention and more rigorous quantitative analysis of their expressibility and complexity is warranted. Our analysis of the restricted Boltzmann machine (RBM) state representation of one-dimensional (1D) quantum spin systems provides new insight into their computational complexity. We define a class of long-range-fast-decay (LRFD) RBM states with quantifiable upper bounds on truncation errors and provide numerical evidence for a large class of 1D quantum systems that may be approximated by LRFD RBMs of at most polynomial complexities. These results lead us to conjecture that the ground states of a wide range of quantum systems may be exactly represented by LRFD RBMs or a variant of them, even in cases where other state representations become less efficient. At last, we provide the relations between multiple typical state manifolds. Our work proposes a paradigm for doing complexity analysis for generic long-range RBMs which naturally yields a further classification of this manifold. This paradigm and our characterization of their nonlocal structures may pave the way for understanding the natural measure of complexity for quantum many-body states described by RBMs and are generalizable for higher-dimensional systems and deep neural-network quantum states. Published by the American Physical Society 2024

Exactly solvable Hamiltonian fragments obtained from a direct sum of Lie algebras.

Exactly solvable Hamiltonians are useful in the study of quantum many-body systems using quantum computers. In the variational quantum eigensolver, a decomposition of the target Hamiltonian into exactly solvable fragments can be used for the evaluation of the energies via repeated quantum measurements. In this work, we apply more general classes of exactly solvable qubit Hamiltonians than previously considered to address the Hamiltonian measurement problem. The most general exactly solvable Hamiltonians we use are defined by the condition that within each simultaneous eigenspace of a set of Pauli symmetries, the Hamiltonian acts effectively as an element of a direct sum of so(N) Lie algebras and can, therefore, be measured using a combination of unitaries in the associated Lie group, Clifford unitaries, and mid-circuit measurements. The application of such Hamiltonians to decomposing molecular electronic Hamiltonians via graph partitioning techniques shows a reduction in the total number of measurements required to estimate the expectation value compared to previously used exactly solvable qubit Hamiltonians.

Numerical investigation on the convergence of self-consistent Schrödinger-Poisson equations in semiconductor device transport simulation

Semiconductor devices at the nanoscale with low-dimensional materials as channels exhibit quantum transport characteristics, thereby their electrical simulation relies on the self-consistent solution of the Schrödinger-Poisson equations. While the non-equilibrium Green’s function (NEGF) method is widely used for solving this quantum many-body problem, its high computational cost and convergence challenges with the Poisson equation significantly limit its applicability. In this study, we investigate the stability of the NEGF method coupled with various forms of the Poisson equation, encompassing linear, analytical nonlinear, and numerical nonlinear forms Our focus lies on simulating carbon nanotube field-effect transistors (CNTFETs) under two distinct doping scenarios: electrostatic doping and ion implantation doping. The numerical experiments reveal that nonlinear formulas outperform linear counterpart. The numerical one demonstrates superior stability, particularly evident under high bias and ion implantation doping conditions. Additionally, we investigate different approaches for presolving potential, leveraging solutions from the Laplace equation and a piecewise guessing method tailored to each doping mode. These methods effectively reduce the number of iterations required for convergence.

Measuring the Boundary Gapless State and Criticality via Disorder Operator.

The disorder operator is often designed to reveal the conformal field theory (CFT) information in quantum many-body systems. By using large-scale quantum MonteCarlo simulation, we study the scaling behavior of disorder operators on the boundary in the two-dimensional Heisenberg model on the square-octagon lattice with gapless topological edge state. In the Affleck-Kennedy-Lieb-Tasaki phase, the disorder operator is shown to hold the perimeter scaling with a logarithmic term associated with the Luttinger liquid parameter K. This effective Luttinger liquid parameter K reflects the low-energy physics and CFT for (1+1)D boundary. At bulk critical point, the effective K is suppressed but it keeps finite value, indicating the coupling between the gapless edge state and bulk fluctuation. The logarithmic term numerically captures this coupling picture, which reveals the (1+1)D SU(2)_{1} CFT and (2+1)D O(3) CFT at boundary criticality. Our Letter paves a new way to study the exotic boundary state and boundary criticality.

A Dyson Brownian Motion Model for Weak Measurements in Chaotic Quantum Systems

We consider a toy model for the study of monitored dynamics in many-body quantum systems. We study the stochastic Schrödinger equation resulting from continuous monitoring with a rate Γ of a random Hermitian operator, drawn from the Gaussian unitary ensemble (GUE) at every time t. Due to invariance by unitary transformations, the dynamics of the eigenvalues {λα}α=1n of the density matrix decouples from that of the eigenvectors, and is exactly described by stochastic equations that we derive. We consider two regimes: in the presence of an extra dephasing term, which can be generated by imperfect quantum measurements, the density matrix has a stationary distribution, and we show that in the limit of large size n→∞ it matches with the inverse-Marchenko–Pastur distribution. In the case of perfect measurements, instead, purification eventually occurs and we focus on finite-time dynamics. In this case, remarkably, we find an exact solution for the joint probability distribution of λ’s at each time t and for each size n. Two relevant regimes emerge: at short times tΓ=O(1), the spectrum is in a Coulomb gas regime, with a well-defined continuous spectral distribution in the n→∞ limit. In that case, all moments of the density matrix become self-averaging and it is possible to exactly characterize the entanglement spectrum. In the limit of large times tΓ=O(n), one enters instead a regime in which the eigenvalues are exponentially separated log(λα/λβ)=O(Γt/n), but fluctuations ∼O(Γt/n) play an essential role. We are still able to characterize the asymptotic behaviors of the entanglement entropy in this regime.

Selective Detection of Dynamics-Complete Set of Correlations via Quantum Channels.

Correlations of fluctuations are essential to understanding many-body systems and key information for advancing quantum technologies. To fully describe the dynamics of a physical system, all time-ordered correlations (TOCs), i.e., the dynamics-complete set of correlations are needed. The current measurement techniques can only access a limited set of TOCs, and there has been no systematic and feasible solution for extracting the dynamic-complete set of correlations hitherto. Here we propose a platform-universal protocol to selectively detect arbitrary types of TOCs via quantum channels. In our method, the quantum channels are synthesized with various controls, and engineer the evolution of a sensor-target system along a specific path that corresponds to a desired correlation. Using nuclear magnetic resonance, we experimentally demonstrate this protocol by detecting a specific type of fourth-order TOC that has never been accessed previously. We also show that the knowledge of the TOCs can be used to significantly improve the precision of quantum optimal control. Our method provides a new toolbox for characterizing the quantum many-body states and quantum noise, and hence for advancing the fields of quantum sensing and quantum computing.

Excited-state downfolding using ground-state formalisms

Downfolding coupled cluster (CC) techniques are powerful tools for reducing the dimensionality of many-body quantum problems. This work investigates how ground-state downfolding formalisms can target excited states using non-Aufbau reference determinants, paving the way for applications of quantum computing in excited-state chemistry. This study focuses on doubly excited states for which canonical equation-of-motion CC approaches struggle to describe unless one includes higher-than-double excitations. The downfolding technique results in state-specific effective Hamiltonians that, when diagonalized in their respective active spaces, provide ground- and excited-state total energies (and therefore excitation energies) comparable to high-level CC methods. The performance of this procedure is examined with doubly excited states of H2, Methylene, Formaldehyde, and Nitroxyl.

Wavefunction matching for solving quantum many-body problems

Ab initio calculations have an essential role in our fundamental understanding of quantum many-body systems across many subfields, from strongly correlated fermions1–3 to quantum chemistry4–6 and from atomic and molecular systems7–9 to nuclear physics10–14. One of the primary challenges is to perform accurate calculations for systems where the interactions may be complicated and difficult for the chosen computational method to handle. Here we address the problem by introducing an approach called wavefunction matching. Wavefunction matching transforms the interaction between particles so that the wavefunctions up to some finite range match that of an easily computable interaction. This allows for calculations of systems that would otherwise be impossible owing to problems such as Monte Carlo sign cancellations. We apply the method to lattice Monte Carlo simulations15,16 of light nuclei, medium-mass nuclei, neutron matter and nuclear matter. We use high-fidelity chiral effective field theory interactions17,18 and find good agreement with empirical data. These results are accompanied by insights on the nuclear interactions that may help to resolve long-standing challenges in accurately reproducing nuclear binding energies, charge radii and nuclear-matter saturation in ab initio calculations19,20.

Sensing Coherent Nuclear Spin Dynamics with an Ensemble of Paramagnetic Nitrogen Spins

The unpolarized spin environment surrounding a central spin qubit is typically considered as an incoherent source of dephasing, however, precise characterization and control of the spin bath can yield a resource for storing and sensing with quantum states. In this work, we use nitrogen-vacancy (NV) centers in diamond to measure the coherence of optically dark paramagnetic nitrogen defects (P1 centers) and detect coherent interactions between the P1 centers and a local bath of 13C nuclear spins. The dipolar coupling between the P1 centers and 13C nuclear spins is identified by signature periodic collapses and revivals in the P1 spin coherence signal. We then demonstrate, using a range of dynamical decoupling protocols, that the probing NV centers and the P1 spins are coupled to independent ensembles of 13C nuclear spins. Our work illustrates how the optically dark P1 spins can be used to extract information from their local environment and offers new insight into the interactions within a many-body system. Published by the American Physical Society 2024

Evidence for a finite-momentum Cooper pair in tricolor d-wave superconducting superlattices

Fermionic superfluidity with a nontrivial Cooper-pairing, beyond the conventional Bardeen-Cooper-Schrieffer state, is a captivating field of study in quantum many-body systems. In particular, the search for superconducting states with finite-momentum pairs has long been a challenge, but establishing its existence has long suffered from the lack of an appropriate probe to reveal its momentum. Recently, it has been proposed that the nonreciprocal electron transport is the most powerful probe for the finite-momentum pairs, because it directly couples to the supercurrents. Here we reveal such a pairing state by the non-reciprocal transport on tricolor superlattices with strong spin-orbit coupling combined with broken inversion-symmetry consisting of atomically thin d-wave superconductor CeCoIn5. We find that while the second-harmonic resistance exhibits a distinct dip anomaly at the low-temperature (T)/high-magnetic field (H) corner in the HT-plane for H applied to the antinodal direction of the d-wave gap, such an anomaly is absent for H along the nodal direction. By carefully isolating extrinsic effects due to vortex dynamics, we reveal the presence of a non-reciprocal response originating from intrinsic superconducting properties characterized by finite-momentum pairs. We attribute the high-field state to the helical superconducting state, wherein the phase of the order parameter is spontaneously spatially modulated.

Directly imaging spin polarons in a kinetically frustrated Hubbard system.

The emergence of quasiparticles in quantum many-body systems underlies the rich phenomenology in many strongly interacting materials. In the context of doped Mott insulators, magnetic polarons are quasiparticles that usually arise from an interplay between the kinetic energy of doped charge carriers and superexchange spin interactions1-8. However, in kinetically frustrated lattices, itinerant spin polarons-bound states of a dopant and a spin flip-have been theoretically predicted even in the absence of superexchange coupling9-14. Despite their important role in the theory of kinetic magnetism, a microscopic observation of these polarons is lacking. Here we directly image itinerant spin polarons in a triangular-lattice Hubbard system realized with ultracold atoms, revealing enhanced antiferromagnetic correlations in the local environment of a hole dopant. In contrast, around a charge dopant, we find ferromagnetic correlations, a manifestation of the elusive Nagaoka effect15,16. We study the evolution of these correlations with interactions and doping, and use higher-order correlation functions to further elucidate the relative contributions of superexchange and kinetic mechanisms. The robustness of itinerant spin polarons at high temperature paves the way for exploring potential mechanisms for hole pairing and superconductivity in frustrated systems10,11. Furthermore, our work provides microscopic insights into related phenomena in triangular-lattice moiré materials17-20.

High-Order Randomized Compiler for Hamiltonian Simulation

Hamiltonian simulation is known to be one of the fundamental building blocks of a variety of quantum algorithms such as its most immediate application, that of simulating many-body systems to extract their physical properties. In this work, we present qSWIFT, a high-order randomized algorithm for Hamiltonian simulation. In qSWIFT, the required number of gates for a given precision is independent of the number of terms in the Hamiltonian, while the systematic error is exponentially reduced with regard to the order parameter. In this respect, our qSWIFT is a higher-order counterpart of the previously proposed quantum stochastic drift protocol (qDRIFT), the number of gates in which scales linearly with the inverse of the precision required. We construct the qSWIFT channel and establish a rigorous bound for the systematic error quantified by the diamond norm. qSWIFT provides an algorithm to estimate given physical quantities by using a system with one ancilla qubit, which is as simple as other product-formula-based approaches such as regular Trotter-Suzuki decompositions and qDRIFT. Our numerical experiment reveals that the required number of gates in qSWIFT is significantly reduced compared to qDRIFT. In particular, the advantage is significant for problems where high precision is required; e.g., to achieve a systematic relative propagation error of 10−6, the required number of gates in third-order qSWIFT is 1000 times smaller than that of qDRIFT. Published by the American Physical Society 2024

Successive quasienergy collapse and breakdown of photon blockade in the few-emitter limit

The emergent behavior that arises in many-body systems of increasing size follows universal laws that become apparent in order-to-disorder transitions. While this behavior has been traditionally studied for large numbers of emitters, recent progress allows for the exploration of the few-emitter limit, where correlations can be measured and connected to microscopic models to gain further insight into order-to-disorder transitions. We explore this few-body limit in the driven and damped Tavis–Cummings model, which describes a collection of atoms interacting with a driven and damped cavity mode. Our exploration revolves around the dressed states of the atomic ensemble and field, whose energies are shown to collapse as the driving field is increased to mark the onset of a dissipative quantum phase transition. The collapse occurs in stages and is an effect of light-matter correlations that are overlooked for single atoms and neglected in mean-field models. The implications of these correlations over the macroscopic observables of the system are presented. We encounter a shift in the expected transition point and an increased number of parity-broken states to choose from once the ordered phase is reached.

Simple and efficient methods for local structural analysis in polydisperse hard disk systems.

In nonequilibrium statistical physics, quantifying the nearest (and higher-order) neighbors and free volumes of particles in many-body systems is crucial to elucidating the origin of macroscopic collective phenomena, such as glass/granular jamming transitions and various aspects of the behavior of active matter. However, conventional techniques (based on a fixed-distance cutoff or the Voronoi construction) have mainly been applied to equilibrated, homogeneous, and monodisperse particle systems. In this paper, we implement simple and efficient methods for local structure analysis in nonequilibrium, inhomogeneous, and polydisperse hard disk systems. We show how these novel methods can overcome the difficulties encountered by conventional techniques as well as demonstrate some applications.

Tensor network representation and entanglement spreading in many-body localized systems: a novel approach

Tensor network representation and entanglement spreading in many-body localized systems: a novel approach

Insulating-to-conducting state transition in the bilayer Hubbard model induced by a perpendicular quench field

A many-body quantum system with varying parameters can exhibit two distinct quantum states within the same energy shell. This allows for a dynamic transition from the ground state of the pre-quench Hamiltonian to a steady state of the post-quench Hamiltonian. We investigate the dynamic response of the ground states in a two-layer half-filled Hubbard model to a perpendicular electric field. We demonstrate that the steady state exhibits conductivity when the field is in resonance with the on-site repulsion, while the initial state is a Mott-insulating state. Additionally, the two layers exhibit identical conducting behavior due to the formation of long-lived dopings, as evidenced by the charge fluctuation. The key factor in achieving this dynamic transition is the cooperative interplay between on-site interactions and the resonant field, rather than the individual roles they play. Our findings offer an alternative mechanism for field-induced conductivity in strongly correlated systems.

Spectral gaps of two- and three-dimensional many-body quantum systems in the thermodynamic limit

We present an expression for the spectral gap, opening up new possibilities for performing and accelerating spectral calculations of quantum many-body systems. We develop and demonstrate one such possibility in the context of tensor network simulations. Our approach requires only minor modifications of the widely used simple update method and is computationally lightweight relative to other approaches. We validate it by computing spectral gaps of the 2D and 3D transverse-field Ising models and find strong agreement with previously reported perturbation theory results. Published by the American Physical Society 2024

Separability criterion using one observable for special states: Entanglement detection via quantum quench

Detecting entanglement in many-body quantum systems is crucial but challenging, typically requiring multiple measurements. Here, we establish a class of states where measuring connected correlations in just basis gives a necessary condition to detect bipartite separability, thus offering a sufficient condition to detect entanglement, provided the appropriate basis and observables are chosen. This methodology leverages prior information about the state, which, although insufficient to reveal the complete state or its entanglement, enables our one basis approach to be effective. We discuss the possibility of one observable entanglement detection in a variety of systems, including those without conserved charges, such as the transverse Ising model, reaching the appropriate basis via quantum quench. This provides a much simpler pathway of detection than previous works. It also shows improved sensitivity from Pearson correlation detection techniques. Published by the American Physical Society 2024