We demonstrate laser wakefield acceleration of quasimonoenergetic electron bunches up to 15 MeV at 1-kHz repetition rate with 2.5-pC charge per bunch and a core with $<7\text{−}\mathrm{mrad}$ beam divergence. Acceleration is driven by 5-fs, $<2.7\text{−}\mathrm{mJ}$ laser incident on a thin, near-critical-density hydrogen gas jet. Low beam divergence is attributed to reduced sensitivity to laser carrier-envelope phase slip, achieved in two ways using gas jet positon control and laser polarization: (i) electron injection into the wake on the gas jet’s plasma density downramp and (ii) use of circularly polarized drive pulses. These results demonstrate the generation of high-quality electron beams from a few-cycle-pulse-driven laser plasma accelerator without the need for carrier-envelope phase stabilization.

In many unconventional superconductors, the presence of a pseudogap—a suppression in the electronic density of states extending above the critical temperature—has been a long-standing mystery. Here, we employ combined in situ electrical transport and angle-resolved photoemission spectroscopy measurements to reveal an unprecedentedly large pseudogap regime in single-layer $\mathrm{FeSe}/\mathrm{SrTi}{\mathrm{O}}_3$, an interfacial superconductor where incoherent Cooper pairs are initially formed above $T_{\mathrm{\Delta }}\approx 60\mathrm{K}$ but where a zero-resistance state is achieved only below $T_0<30\mathrm{K}$. We show that this behavior is accompanied by distinct transport signatures of two-dimensional phase fluctuating superconductivity, suggesting a mixed vortex state hosting incoherent Cooper pairs which persist well above the maximum clean limit $T_c$ of approximately 40 K. Our work establishes the critical role of reduced dimensionality in driving the complex interplay between Cooper pairing and phase coherence in two-dimensional high-$T_c$ superconductors, providing a paradigm for understanding and engineering higher-$T_c$ interfacial superconductors.

We report on gravitational-wave discoveries from compact binary coalescences detected by Advanced LIGO and Advanced Virgo in the first half of the third observing run (O3a) between 1 April 2019 $15:00$ UTC and 1 October 2019 $15:00$ UTC. By imposing a false-alarm-rate threshold of two per year in each of the four search pipelines that constitute our search, we present 39 candidate gravitational-wave events. At this threshold, we expect a contamination fraction of less than 10%. Of these, 26 candidate events were reported previously in near-real time through gamma-ray coordinates network notices and circulars; 13 are reported here for the first time. The catalog contains events whose sources are black hole binary mergers up to a redshift of approximately 0.8, as well as events whose components cannot be unambiguously identified as black holes or neutron stars. For the latter group, we are unable to determine the nature based on estimates of the component masses and spins from gravitational-wave data alone. The range of candidate event masses which are unambiguously identified as binary black holes (both objects $\ge 3M_{\odot }$) is increased compared to GWTC-1, with total masses from approximately $14M_{\odot }$ for GW190924_021846 to approximately $150M_{\odot }$ for GW190521. For the first time, this catalog includes binary systems with significantly asymmetric mass ratios, which had not been observed in data taken before April 2019. We also find that 11 of the 39 events detected since April 2019 have positive effective inspiral spins under our default prior (at 90% credibility), while none exhibit negative effective inspiral spin. Given the increased sensitivity of Advanced LIGO and Advanced Virgo, the detection of 39 candidate events in approximately 26 weeks of data (approximately 1.5 per week) is consistent with GWTC-1.

The design of periodic nanostructures allows to tailor the transport of photons, phonons, and matter waves for specific applications. Recent years have seen a further expansion of this field by engineering topological properties. However, what is missing currently are efficient ways to rapidly explore and optimize band structures and to classify their topological characteristics for arbitrary unit-cell geometries. In this work, we show how deep learning can address this challenge. We introduce an approach where a neural network first maps the geometry to a tight-binding model. The tight-binding model encodes not only the band structure but also the symmetry properties of the Bloch waves. This allows us to rapidly categorize a large set of geometries in terms of their band representations, identifying designs for fragile topologies. We demonstrate that our method is also suitable to calculate strong topological invariants, even when (like the Chern number) they are not symmetry indicated. Engineering of domain walls and optimization are accelerated by orders of magnitude. Our method directly applies to any passive linear material, irrespective of the symmetry class and space group. It is general enough to be extended to active and nonlinear metamaterials.

We study many-body quantum dynamics using Floquet quantum circuits in one space dimension as simple examples of systems with local interactions that support ergodic phases. Physical properties can be expressed in terms of multiple sums over Feynman histories, which for these models are paths or many-body orbits in Fock space. A natural simplification of such sums is the diagonal approximation, where the only terms that are retained are ones in which each path is paired with a partner that carries the complex conjugate weight. We identify the regime in which the diagonal approximation holds and the nature of the leading corrections to it. We focus on the behavior of the spectral form factor (SFF) and of matrix elements of local operators, averaged over an ensemble of random circuits, making comparisons with the predictions of random matrix theory (RMT) and the eigenstate thermalization hypothesis (ETH). We show that properties are dominated at long times by contributions to orbit sums in which each orbit is paired locally with a conjugate, as in the diagonal approximation, but that in large systems these contributions consist of many spatial domains, with distinct local pairings in neighboring domains. The existence of these domains is reflected in deviations of the SFF from RMT predictions, and of matrix element correlations from ETH predictions; deviations of both kinds diverge with system size. We demonstrate that our physical picture of orbit-pairing domains has a precise correspondence in the spectral properties of a transfer matrix that acts in the space direction to generate the ensemble-averaged SFF. In addition, we find that domains of a second type control non-Gaussian fluctuations of the SFF. These domains are separated by walls that are related to the entanglement membrane, known to characterize the scrambling of quantum information.

Metasurfaces have been enabling the miniaturization and integration of complex optical functionalities within an ultrathin platform by engineering the scattering features of localized modes. However, these efforts have mostly been limited to the manipulation of externally produced coherent light, e.g., from a laser. In parallel, the past two decades have seen the development of structured surfaces that emit partially coherent radiation via thermally populated, spatially extended (nonlocal) modes. However, the control over thermally emitted light is severely limited compared to optical metasurfaces, and even basic functionalities such as unidirectional emission to an arbitrary angle and polarization remain elusive. Here, we derive the necessary conditions to achieve full control over thermally emitted light, pointing to the need for simultaneously tailoring local and nonlocal scattering features across the structure. Based on these findings, we introduce a platform for thermal metasurfaces based on quasibound states in the continuum that satisfies these requirements and completes the program of compactification of optical systems by enabling a full degree of control of partially coherent light emission from structured thin films, including unidirectional emission of circularly polarized light, focusing, and control of spatial and temporal coherence, as well as wave-front control with designer spin and angular orbital momenta.

The acceleration of super-heavy ions (SHIs, mass number of about 200) from plasmas driven by ultrashort (tens of femtoseconds) laser pulses is a challenging topic awaiting a breakthrough. Detecting and controlling the ionization process and adopting the optimal acceleration scheme are crucial for the generation of highly energetic SHIs. Here, we report the experimental results on the generation of deeply ionized super-heavy ions (Au) with unprecedented energy of 1.2 GeV utilizing ultrathin targets and ultrashort laser pulses at an intensity of ${10}^{22}\mathrm{W}/{\mathrm{cm}}^2$. A novel self-calibrated diagnostic method was developed to acquire the absolute energy spectra and charge-state distributions of Au ions abundant at the charge state of $51+$ and extending to $61+$. The measured charge-state distributions supported by 2D particle-in-cell simulations serve as an additional tool to inspect the ionization dynamics associated with SHI acceleration, revealing that the laser intensity is the crucial parameter over the pulse duration for Au acceleration. Achieving a long acceleration time without sacrificing the strength of the acceleration field by utilizing composite targets can substantially increase the maximum energy of Au ions.

We introduce a near-term experimental platform for realizing an associative memory. It can simultaneously store many memories by using spinful bosons coupled to a degenerate multimode optical cavity. The associative memory is realized by a confocal cavity QED neural network, with the modes serving as the synapses, connecting a network of superradiant atomic spin ensembles,which serve as the neurons. Memories are encoded in the connectivity matrix between the spins and can be accessed through the input and output of patterns of light. Each aspect of the scheme is based on recently demonstrated technology using a confocal cavity and Bose-condensed atoms. Our scheme has two conceptually novel elements. First, it introduces a new form of random spin system that interpolates between a ferromagnetic and a spin glass regime as a physical parameter is tuned—the positions of ensembles within the cavity. Second, and more importantly, the spins relax via deterministic steepest-descent dynamics rather than Glauber dynamics. We show that this nonequilibrium quantum-optical scheme has significant advantages for associative memory over Glauber dynamics: These dynamics can enhance the network’s ability to store and recall memories beyond that of the standard Hopfield model. Surprisingly, the cavity QED dynamics can retrieve memories even when the system is in the spin glass phase. Thus, the experimental platform provides a novel physical instantiation of associative memories and spin glasses as well as provides an unusual form of relaxational dynamics that is conducive to memory recall even in regimes where it was thought to be impossible.

Variational quantum circuits (VQCs) built upon noisy intermediate-scale quantum (NISQ) hardware, in conjunction with classical processing, constitute a promising architecture for quantum simulations, classical optimization, and machine learning. However, the required VQC depth to demonstrate a quantum advantage over classical schemes is beyond the reach of available NISQ devices. Supervised learning assisted by an entangled sensor network (SLAEN) is a distinct paradigm that harnesses VQCs trained by classical machine-learning algorithms to tailor multipartite entanglement shared by sensors for solving practically useful data-processing problems. Here, we report the first experimental demonstration of SLAEN and show an entanglement-enabled reduction in the error probability for classification of multidimensional radio-frequency signals. Our work paves a new route for quantum-enhanced data processing and its applications in the NISQ era.

Materials and machines are often designed with particular goals in mind, so that they exhibit desired responses to given forces or constraints. Here we explore an alternative approach, namely physical coupled learning. In this paradigm, the system is not initially designed to accomplish a task, but physically adapts to applied forces to develop the ability to perform the task. Crucially, we require coupled learning to be facilitated by physically plausible learning rules, meaning that learning requires only local responses and no explicit information about the desired functionality. We show that such local learning rules can be derived for any physical network, whether in equilibrium or in steady state, with specific focus on two particular systems, namely disordered flow networks and elastic networks. By applying and adapting advances of statistical learning theory to the physical world, we demonstrate the plausibility of new classes of smart metamaterials capable of adapting to users’ needs in situ.

Triangular lattice of rare-earth ions with interacting effective spin-$1/2$ local moments is an ideal platform to explore the physics of quantum spin liquids (QSLs) in the presence of strong spin-orbit coupling, crystal electric fields, and geometrical frustration. The Yb delafossites, $\mathrm{NaYb}C{h}_2$ ($Ch=\mathrm{O}$, S, Se) with Yb ions forming a perfect triangular lattice, have been suggested to be candidates for QSLs. Previous thermodynamics, nuclear magnetic resonance, and powder-sample neutron scattering measurements on ${\mathrm{NaYbCh}}_2$ have supported the suggestion of the QSL ground states. The key signature of a QSL, the spin excitation continuum, arising from the spin quantum number fractionalization, has not been observed. Here we perform both elastic and inelastic neutron scattering measurements as well as detailed thermodynamic measurements on high-quality single-crystal ${\mathrm{NaYbSe}}_2$ samples to confirm the absence of long-range magnetic order down to 40 mK, and further reveal a clear signature of magnetic excitation continuum extending from 0.1 to 2.5 meV. The comparison between the structure of the magnetic excitation spectra and the theoretical expectation from the spinon continuum suggests that the ground state of ${\mathrm{NaYbSe}}_2$ is a QSL with a spinon Fermi surface.

Causality is a seminal concept in science: Any research discipline, from sociology and medicine to physics and chemistry, aims at understanding the causes that could explain the correlations observed among some measured variables. While several methods exist to characterize classical causal models, no general construction is known for the quantum case. In this work, we present quantum inflation, a systematic technique to falsify if a given quantum causal model is compatible with some observed correlations. We demonstrate the power of the technique by reproducing known results and solving open problems for some paradigmatic examples of causal networks. Our results may find applications in many fields: from the characterization of correlations in quantum networks to the study of quantum effects in thermodynamic and biological processes.

In twisted homobilayer transition metal dichalcogenides, intra- and interlayer valley excitons hybridize with the layer configurations spatially varying in the moiré. The ground state valley excitons are trapped at two high-symmetry points with opposite electric dipoles in a moiré supercell, forming a honeycomb superlattice of nearest-neighbor dipolar attraction. We find that the spatial texture of layer configuration results in a luminescence anomaly of the moiré trapped excitons, where a tiny displacement by interactions dramatically increases the brightness and changes polarization from circular to linear. At full filling, radiative recombination predominantly occurs at edges and vacancies of the exciton superlattice. The anomaly also manifests in the cascaded emission of small clusters, producing chains of polarization entangled photons. An interlayer bias can switch the superlattice into a single-orbital triangular lattice with repulsive interactions only, where the luminescence anomaly can be exploited to distinguish ordered states and domain boundaries at fractional filling.

Two widely used but distinct approaches to the dynamics of open quantum systems are the Nakajima-Zwanzig and time-convolutionless quantum master equation, respectively. Although both describe identical quantum evolutions with strong memory effects, the first uses a time-nonlocal memory kernel $\mathcal{K}$, whereas the second achieves the same using a time-local generator $\mathcal{G}$. Here we show that the two are connected by a simple yet general fixed-point relation: $\mathcal{G}=\widehat{\mathcal{K}}[\mathcal{G}]$. This allows one to extract nontrivial relations between the two completely different ways of computing the time evolution and combine their strengths. We first discuss the stationary generator, which enables a Markov approximation that is both nonperturbative and completely positive for a large class of evolutions. We show that this generator is not equal to the low-frequency limit of the memory kernel, but additionally “samples” it at nonzero characteristic frequencies. This clarifies the subtle roles of frequency dependence and semigroup factorization in existing Markov approximation strategies. Second, we prove that the fixed-point equation sums up the time-domain gradient or Moyal expansion for the time-nonlocal quantum master equation, providing nonperturbative insight into the generation of memory effects. Finally, we show that the fixed-point relation enables a direct iterative numerical computation of both the stationary and the transient generator from a given memory kernel. For the transient generator this produces nonsemigroup approximations which are constrained to be both initially and asymptotically accurate at each iteration step.

Recent experimental and theoretical works have made much progress toward understanding nonequilibrium phenomena in thermalizing systems, which act as thermal baths for their small subsystems, and many-body localized systems, which fail to do so. The description of time evolution in many-body systems is generally challenging due to the dynamical generation of quantum entanglement. In this work, we introduce an approach to study quantum many-body dynamics, inspired by the Feynman-Vernon influence functional. Focusing on a family of interacting, Floquet spin chains, we consider a Keldysh path-integral description of the dynamics. The central object in our approach is the influence matrix, which describes the effect of the system on the dynamics of a local subsystem. For translationally invariant models, we formulate a self-consistency equation for the influence matrix. For certain special values of the model parameters, we obtain an exact solution which represents a perfect dephaser (PD). Physically, a PD corresponds to a many-body system that acts as a perfectly Markovian bath on itself: at each period, it measures every spin. For the models considered here, we establish that PD points include dual-unitary circuits investigated in recent works. In the vicinity of PD points, the system is not perfectly Markovian, but rather acts as a bath with a short memory time. In this case, we demonstrate that the self-consistency equation can be solved using matrix-product states (MPS) methods, as the influence matrix temporal entanglement is low. A combination of analytical insights and MPS computations allows us to characterize the structure of the influence matrix in terms of an effective “statistical-mechanics” description. We finally illustrate the predictive power of this description by analytically computing how quickly an embedded impurity spin thermalizes. The influence matrix approach formulated here provides an intuitive view of the quantum many-body dynamics problem, opening a path to constructing models of thermalizing dynamics that are solvable or can be efficiently treated by MPS-based methods and to further characterizing quantum ergodicity or lack thereof.

Quantum many-body states that frequently appear in physics often obey an entropy scaling law, meaning that an entanglement entropy of a subsystem can be expressed as a sum of terms that scale linearly with its volume and area, plus a correction term that is independent of its size. We conjecture that these states have an efficient dual description in terms of a set of marginal density matrices on bounded regions, obeying the same entropy scaling law locally. We prove a restricted version of this conjecture for translationally invariant systems in two spatial dimensions. Specifically, we prove that a translationally invariant marginal obeying three nonlinear constraints—all of which follow from the entropy scaling law straightforwardly—must be consistent with some global state on an infinite lattice. Moreover, we derive a closed-form expression for the maximum-entropy density compatible with those marginals, deriving a variational upper bound on the thermodynamic free energy. The main assumptions of our construction are satisfied exactly by solvable models of topological order and approximately by finite-temperature Gibbs states of certain quantum spin Hamiltonians.

The mechanism behind the nematicity of FeSe is not known. Through elastoresitivity measurements it has been shown to be an electronic instability. However, thus far measurements have extended only to small strains, where the response is linear. Here, we apply large elastic strains to FeSe and perform two types of measurement. (1) Using applied strain to control twinning, the nematic resistive anisotropy at temperatures below the nematic transition temperature $T_s$ is determined. (2) Resistive anisotropy is measured as nematicity is induced through applied strain at fixed temperature above $T_s$. In both cases, as nematicity strengthens, the resistive anisotropy peaks at about 7%, then decreases. Below $\approx 40\mathrm{K}$, the nematic resistive anisotropy changes sign. We discuss possible implications of this behavior for theories of nematicity. In addition, we report the following. (1) Under experimentally accessible conditions with bulk crystals, stress, rather than strain, is the conjugate field to the nematicity of FeSe. (2) At low temperatures the twin boundary resistance is $\sim 10\%$ of the sample resistance, and must be properly subtracted to extract intrinsic resistivities. (3) Biaxial in-plane compression increases both in-plane resistivity and the superconducting critical temperature $T_c$, consistent with a strong role of the $yz$ orbital in the electronic correlations.

We present a full symmetry classification of fermion matter in and out of thermal equilibrium. Our approach starts from first principles, the ten different classes of linear and antilinear state transformations in fermionic Fock spaces, and symmetries defined via invariance properties of the dynamical equation for the density matrix. The object of classification is then the generators of reversible dynamics, dissipation and fluctuations, featuring’ in the generally irreversible and interacting dynamical equations. A sharp distinction between the symmetries of equilibrium and out-of-equilibrium dynamics, respectively, arises from the different role played by “time” in these two cases: In unitary quantum mechanics as well as in “microreversible” thermal equilibrium, antilinear transformations combined with an inversion of time define time-reversal symmetry. However, out of equilibrium an inversion of time becomes meaningless, while antilinear transformations in Fock space remain physically significant, and hence must be considered in autonomy. The practical consequence of this dichotomy is a novel realization of antilinear symmetries (six out of the ten fundamental classes) in nonequilibrium quantum dynamics that is fundamentally different from the established rules of thermal equilibrium. At large times, the dynamical generators thus symmetry classified determine the steady-state nonequilibrium distributions for arbitrary interacting systems. To illustrate this principle, we consider the fixation of a symmetry protected topological phase in a system of interacting lattice fermions. More generally, we consider the practically important class of mean field interacting systems, represented by Gaussian states. This class is naturally described in the language of non-Hermitian matrices, which allows us to compare to previous classification schemes in the literature.