We construct a two-dimensional lattice model of fermions coupled to Ising ferromagnetic critical fluctuations. Using extensive sign-problem-free quantum Monte Carlo simulations, we show that the model realizes a continuous itinerant quantum phase transition. In comparison with other similar itinerant quantum critical points (QCPs), our QCP shows a much weaker superconductivity tendency with no superconducting state down to the lowest temperature investigated, hence making the system a good platform for the exploration of quantum critical fluctuations. Remarkably, clear signatures of non-Fermi liquid behavior in the fermion propagators are observed at the QCP. The critical fluctuations at the QCP partially resemble Hertz-Millis-Moriya behavior. However, careful scaling analysis reveals that the QCP belongs to a different universality class, deviating from both $(2+1)\mathrm{D}$ Ising and Hertz-Millis-Moriya predictions.

The Kondo insulator samarium hexaboride (${\mathrm{SmB}}_6$) has been intensely studied in recent years as a potential candidate of a strongly correlated topological insulator. One of the most exciting phenomena observed in ${\mathrm{SmB}}_6$ is the clear quantum oscillations appearing in magnetic torque at a low temperature despite the insulating behavior in resistance. These quantum oscillations show multiple frequencies and varied effective masses. The origin of quantum oscillation is, however, still under debate with evidence of both two-dimensional Fermi surfaces and three-dimensional Fermi surfaces. Here, we carry out angle-resolved torque magnetometry measurements in a magnetic field up to 45 T and a temperature range down to 40 mK. With the magnetic field rotated in the (010) plane, the quantum oscillation frequency of the strongest oscillation branch shows a fourfold rotational symmetry. However, in the angular dependence of the amplitude of the same branch, this fourfold symmetry is broken and, instead, a twofold symmetry shows up, which is consistent with the prediction of a two-dimensional Lifshitz-Kosevich model. No deviation of Lifshitz-Kosevich behavior is observed down to 40 mK. Our results suggest the existence of multiple light-mass surface states in ${\mathrm{SmB}}_6$, with their mobility significantly depending on the surface disorder level.

Geometric and electronic surface reconstructions determine the physical and chemical properties of surfaces and, consequently, their functionality in applications. The reconstruction of a surface minimizes its surface free energy in otherwise thermodynamically unstable situations, typically caused by dangling bonds, lattice stress, or a divergent surface potential, and it is achieved by a cooperative modification of the atomic and electronic structure. Here, we combined first-principles calculations and surface techniques (scanning tunneling microscopy, non-contact atomic force microscopy, scanning tunneling spectroscopy) to report that the repulsion between negatively charged polaronic quasiparticles, formed by the interaction between excess electrons and the lattice phonon field, plays a key role in surface reconstructions. As a paradigmatic example, we explain the ($1\times 1$) to ($1\times 2$) transition in rutile ${\mathrm{TiO}}_2(110)$.

Recently, significant progress has been made in ($2+1$)-dimensional conformal field theories without supersymmetry. In particular, it was realized that different Lagrangians may be related by hidden dualities; i.e., seemingly different field theories may actually be identical in the infrared limit. Among all the proposed dualities, one has attracted particular interest in the field of strongly correlated quantum-matter systems: the one relating the easy-plane noncompact ${\mathrm{CP}}^1$ model (${\mathrm{NCCP}}^1$) and noncompact quantum electrodynamics (QED) with two flavors ($N=2$) of massless two-component Dirac fermions. The easy-plane ${\mathrm{NCCP}}^1$ model is the field theory of the putative deconfined quantum-critical point separating a planar ($XY$) antiferromagnet and a dimerized (valence-bond solid) ground state, while $N=2$ noncompact QED is the theory for the transition between a bosonic symmetry-protected topological phase and a trivial Mott insulator. In this work, we present strong numerical support for the proposed duality. We realize the $N=2$ noncompact QED at a critical point of an interacting fermion model on the bilayer honeycomb lattice and study it using determinant quantum Monte Carlo (QMC) simulations. Using stochastic series expansion QMC simulations, we study a planar version of the $S=1/2$ $J\text{−}Q$ spin Hamiltonian (a quantum $XY$ model with additional multispin couplings) and show that it hosts a continuous transition between the $XY$ magnet and the valence-bond solid. The duality between the two systems, following from a mapping of their phase diagrams extending from their respective critical points, is supported by the good agreement between the critical exponents according to the proposed duality relationships. In the $J\text{−}Q$ model, we find both continuous and first-order transitions, depending on the degree of planar anisotropy, with deconfined quantum criticality surviving only up to moderate strengths of the anisotropy. This explains previous claims of no deconfined quantum criticality in planar two-component spin models, which were in the strong-anisotropy regime, and opens doors to further investigations of the global phase diagram of systems hosting deconfined quantum-critical points.

The deconfined quantum critical point (QCP), separating the Néel and valence bond solid phases in a 2D antiferromagnet, was proposed as an example of $(2+1)\mathrm{D}$ criticality fundamentally different from standard Landau-Ginzburg-Wilson-Fisher criticality. In this work, we present multiple equivalent descriptions of deconfined QCPs, and use these to address the possibility of enlarged emergent symmetries in the low-energy limit. The easy-plane deconfined QCP, besides its previously discussed self-duality, is dual to $N_f=2$ fermionic quantum electrodynamics, which has its own self-duality and hence may have an $\mathrm{O}(4)\times Z_2^T$ symmetry. We propose several dualities for the deconfined QCP with SU(2) spin symmetry which together make natural the emergence of a previously suggested SO(5) symmetry rotating the Néel and valence bond solid orders. These emergent symmetries are implemented anomalously. The associated infrared theories can also be viewed as surface descriptions of $(3+1)\mathrm{D}$ topological paramagnets, giving further insight into the dualities. We describe a number of numerical tests of these dualities. We also discuss the possibility of “pseudocritical” behavior for deconfined critical points, and the meaning of the dualities and emergent symmetries in such a scenario.

We present a scalable architecture for fault-tolerant topological quantum computation using networks of voltage-controlled Majorana Cooper pair boxes and topological color codes for error correction. Color codes have a set of transversal gates which coincides with the set of topologically protected gates in Majorana-based systems, namely, the Clifford gates. In this way, we establish color codes as providing a natural setting in which advantages offered by topological hardware can be combined with those arising from topological error-correcting software for full-fledged fault-tolerant quantum computing. We provide a complete description of our architecture, including the underlying physical ingredients. We start by showing that in topological superconductor networks, hexagonal cells can be employed to serve as physical qubits for universal quantum computation, and we present protocols for realizing topologically protected Clifford gates. These hexagonal-cell qubits allow for a direct implementation of open-boundary color codes with ancilla-free syndrome read-out and logical $T$ gates via magic-state distillation. For concreteness, we describe how the necessary operations can be implemented using networks of Majorana Cooper pair boxes, and we give a feasibility estimate for error correction in this architecture. Our approach is motivated by nanowire-based networks of topological superconductors, but it could also be realized in alternative settings such as quantum-Hall–superconductor hybrids.

We study scrambling, an avatar of chaos, in a weakly interacting metal in the presence of random potential disorder. It is well known that charge and heat spread via diffusion in such an interacting disordered metal. In contrast, we show within perturbation theory that chaos spreads in a ballistic fashion. The squared anticommutator of the electron-field operators inherits a light-cone-like growth, arising from an interplay of a growth (Lyapunov) exponent that scales as the inelastic electron scattering rate and a diffusive piece due to the presence of disorder. In two spatial dimensions, the Lyapunov exponent is universally related at weak coupling to the sheet resistivity. We are able to define an effective temperature-dependent butterfly velocity, a speed limit for the propagation of quantum information that is much slower than microscopic velocities such as the Fermi velocity and that is qualitatively similar to that of a quantum critical system with a dynamical critical exponent $z>1$.

Recent advances in engineered gradient metasurfaces have enabled unprecedented opportunities for light manipulation using optically thin sheets, such as anomalous refraction, reflection, or focusing of an incident beam. Here, we introduce a concept of multichannel functional metasurfaces, which are able to control incoming and outgoing waves in a number of propagation directions simultaneously. In particular, we reveal a possibility to engineer multichannel reflectors. Under the assumption of reciprocity and energy conservation, we find that there exist three basic functionalities of such reflectors: specular, anomalous, and retroreflections. Multichannel response of a general flat reflector can be described by a combination of these functionalities. To demonstrate the potential of the introduced concept, we design and experimentally test three different multichannel reflectors: three- and five-channel retroreflectors and a three-channel power splitter. Furthermore, by extending the concept to reflectors supporting higher-order Floquet harmonics, we forecast the emergence of other multichannel flat devices, such as isolating mirrors, complex splitters, and multi-functional gratings.

The Seebeck coefficient $S$ of the cuprate ${\mathrm{YBa}}_2{\mathrm{Cu}}_3{\mathrm{O}}_y$ is measured in magnetic fields large enough to suppress superconductivity, at hole dopings $p=0.11$ and $p=0.12$, for heat currents along the $a$ and $b$ directions of the orthorhombic crystal structure. For both directions, $S/T$ decreases and becomes negative at low temperature, a signature that the Fermi surface undergoes a reconstruction due to broken translational symmetry. Above a clear threshold field, a strong new feature appears in $S_b$, for conduction along the $b$ axis only. We attribute this feature to the onset of 3D-coherent unidirectional charge-density-wave modulations seen by x-ray diffraction, also along the $b$ axis only. Because these modulations have a sharp onset temperature well below the temperature where $S/T$ starts to drop towards negative values, we infer that they are not the cause of Fermi-surface reconstruction. Instead, the reconstruction must be caused by the quasi-2D bidirectional modulations that develop at significantly higher temperature. The unidirectional order only confers an additional anisotropy to the already reconstructed Fermi surface, also manifest as an in-plane anisotropy of the resistivity.

The nitrogen-vacancy (NV) center in diamond has an optically addressable, highly coherent spin. However, a NV center even in high-quality single-crystalline material is a very poor source of single photons: Extraction out of the high-index diamond is inefficient, the emission of coherent photons represents just a few percent of the total emission, and the decay time is large. In principle, all three problems can be addressed with a resonant microcavity. In practice, it has proved difficult to implement this concept: Photonic engineering hinges on nanofabrication, yet it is notoriously difficult to process diamond without degrading the NV centers. Here, we present a microcavity scheme that uses minimally processed diamond, thereby preserving the high quality of the starting material and a tunable microcavity platform. We demonstrate a clear change in the lifetime for multiple individual NV centers on tuning both the cavity frequency and antinode position, a Purcell effect. The overall Purcell factor $F_P=2.0$ translates to a Purcell factor for the zero phonon line (ZPL) of $F_P^{\mathrm{ZPL}}\sim 30$ and an increase in the ZPL emission probability from about 3% to 46%. By making a step change in the NV’s optical properties in a deterministic way, these results pave the way for much enhanced spin-photon and spin-spin entanglement rates.

Machine learning offers an unprecedented perspective for the problem of classifying phases in condensed matter physics. We employ neural-network machine learning techniques to distinguish finite-temperature phases of the strongly correlated fermions on cubic lattices. We show that a three-dimensional convolutional network trained on auxiliary field configurations produced by quantum Monte Carlo simulations of the Hubbard model can correctly predict the magnetic phase diagram of the model at the average density of one (half filling). We then use the network, trained at half filling, to explore the trend in the transition temperature as the system is doped away from half filling. This transfer learning approach predicts that the instability to the magnetic phase extends to at least 5% doping in this region. Our results pave the way for other machine learning applications in correlated quantum many-body systems.

We report on the formation of a diverse family of transverse spatial polygon patterns in a microcavity polariton fluid under coherent driving by a blue-detuned pump. Patterns emerge spontaneously as a result of energy-degenerate polariton-polariton scattering from the pump state to interfering high-order vortex and antivortex modes, breaking azimuthal symmetry. The interplay between a multimode parametric instability and intrinsic optical bistability leads to a sharp spike in the value of second-order coherence $g^{(2)}(0)$ of the emitted light, which we attribute to the strongly superlinear kinetics of the underlying scattering processes driving the formation of patterns. We show numerically by means of a linear stability analysis how the growth of parametric instabilities in our system can lead to spontaneous symmetry breaking, predicting the formation and competition of different pattern states in good agreement with experimental observations.

Multidimensional coherent optical spectroscopy is one of the most powerful tools for investigating complex quantum mechanical systems. While it was conceived decades ago in magnetic resonance spectroscopy using microwaves and radio waves, it has recently been extended into the visible and UV spectral range. However, resolving MHz energy splittings with ultrashort laser pulses still remains a challenge. Here, we analyze two-dimensional Fourier spectra for resonant optical excitation of resident electrons to localized trions or donor-bound excitons in semiconductor nanostructures subject to a transverse magnetic field. Particular attention is devoted to Raman coherence spectra, which allow one to accurately evaluate tiny splittings of the electron ground state and to determine the relaxation times in the electron spin ensemble. A stimulated steplike Raman process induced by a sequence of two laser pulses creates a coherent superposition of the ground-state doublet which can be retrieved only optically because of selective excitation of the same subensemble with a third pulse. This provides the unique opportunity to distinguish between different complexes that are closely spaced in energy in an ensemble. The related experimental demonstration is based on photon-echo measurements in an $n$-type $\mathrm{CdTe}/(\mathrm{Cd},\mathrm{Mg})\mathrm{Te}$ quantum-well structure detected by a heterodyne technique. The difference in the sub-$\mu \mathrm{eV}$ range between the Zeeman splittings of donor-bound electrons and electrons localized at potential fluctuations can be resolved even though the homogeneous linewidth of the optical transitions is larger by 2 orders of magnitude.

It is well known that there is a particle-hole symmetry for spin-polarized electrons with two-body interactions in a partially filled Landau level, which becomes exact in the limit where the cyclotron energy is large compared to the interaction strength; thus, one can ignore mixing between Landau levels. This symmetry is explicit in the description of a half-filled Landau level recently introduced by Son, using Dirac fermions, but it was thought to be absent in the older fermion-Chern-Simons approach, developed by Halperin, Lee, and Read (HLR) and subsequent authors. We show here, however, that when properly evaluated, the HLR theory gives results for long-wavelength low-energy physical properties—including the Hall conductance in the presence of impurities and the positions of minima in the magnetoroton spectra for fractional quantized Hall states close to half-filling—that are identical to predictions of the Dirac formulation. In fact, the HLR theory predicts an emergent particle-hole symmetry near half-filling, even when the cyclotron energy is finite.

The antiferromagnetic spin-$1/2$ Heisenberg model on a kagome lattice is one of the most paradigmatic models in the context of spin liquids, yet the precise nature of its ground state is not understood. We use large-scale density matrix renormalization group simulations (DMRG) on infinitely long cylinders and find indications for the formation of a gapless Dirac spin liquid. First, we use adiabatic flux insertion to demonstrate that the spin gap is much smaller than estimated from previous DMRG simulation. Second, we find that the momentum-dependent excitation spectrum, as extracted from the DMRG transfer matrix, exhibits Dirac cones that match those of a $\pi$-flux free-fermion model [the parton mean-field ansatz of a $U(1)$ Dirac spin liquid].

Characterizing how entanglement grows with time in a many-body system, for example, after a quantum quench, is a key problem in nonequilibrium quantum physics. We study this problem for the case of random unitary dynamics, representing either Hamiltonian evolution with time-dependent noise or evolution by a random quantum circuit. Our results reveal a universal structure behind noisy entanglement growth, and also provide simple new heuristics for the “entanglement tsunami” in Hamiltonian systems without noise. In 1D, we show that noise causes the entanglement entropy across a cut to grow according to the celebrated Kardar-Parisi-Zhang (KPZ) equation. The mean entanglement grows linearly in time, while fluctuations grow like $(\text{time}{)}^{1/3}$ and are spatially correlated over a distance $\propto (\text{time}{)}^{2/3}$. We derive KPZ universal behavior in three complementary ways, by mapping random entanglement growth to (i) a stochastic model of a growing surface, (ii) a “minimal cut” picture, reminiscent of the Ryu-Takayanagi formula in holography, and (iii) a hydrodynamic problem involving the dynamical spreading of operators. We demonstrate KPZ universality in 1D numerically using simulations of random unitary circuits. Importantly, the leading-order time dependence of the entropy is deterministic even in the presence of noise, allowing us to propose a simple coarse grained minimal cut picture for the entanglement growth of generic Hamiltonians, even without noise, in arbitrary dimensionality. We clarify the meaning of the “velocity” of entanglement growth in the 1D entanglement tsunami. We show that in higher dimensions, noisy entanglement evolution maps to the well-studied problem of pinning of a membrane or domain wall by disorder.

We present a tensor network especially suited for multi-orbital Anderson impurity models and as an impurity solver for multi-orbital dynamical mean-field theory (DMFT). The solver works directly on the real-frequency axis and yields high spectral resolution at all frequencies. We use a large number $\mathbf{(}\mathcal{O}(100)\mathbf{)}$ of bath sites and therefore achieve an accurate representation of the bath. The solver can treat full rotationally invariant interactions with reasonable numerical effort. We show the efficiency and accuracy of the method by a benchmark for the three-orbital test-bed material ${\mathrm{SrVO}}_3$. There we observe multiplet structures in the high-energy spectrum, which are almost impossible to resolve by other multi-orbital methods. The resulting structure of the Hubbard bands can be described as a broadened atomic spectrum with rescaled interaction parameters. Additional features emerge when $U$ is increased. Finally, we show that our solver can be applied even to models with five orbitals. This impurity solver offers a new route to the calculation of precise real-frequency spectral functions of correlated materials.

The Rashba and Dresselhaus spin-orbit (SO) interactions in 2D electron gases act as effective magnetic fields with momentum-dependent directions, which cause spin decay as the spins undergo arbitrary precessions about these randomly oriented SO fields due to momentum scattering. Theoretically and experimentally, it has been established that by fine-tuning the Rashba $\alpha$ and renormalized Dresselhaus $\beta$ couplings to equal fixed strengths $\alpha =\beta$, the total SO field becomes unidirectional, thus rendering the electron spins immune to decay due to momentum scattering. A robust persistent spin helix (PSH), i.e., a helical spin-density wave excitation with constant pitch $P=2\pi /Q$, $Q=4m\alpha /{\hslash }^2$, has already been experimentally realized at this singular point $\alpha =\beta$, enhancing the spin lifetime by up to 2 orders of magnitude. Here, we employ the suppression of weak antilocalization as a sensitive detector for matched SO fields together with independent electrical control over the SO couplings via top gate voltage $V_T$ and back gate voltage $V_B$ to extract all SO couplings when combined with detailed numerical simulations. We demonstrate for the first time the gate control of the renormalized $\beta$ and the continuous locking of the SO fields at $\alpha =\beta$; i.e., we are able to vary both $\alpha$ and $\beta$ controllably and continuously with $V_T$ and $V_B$, while keeping them locked at equal strengths. This makes possible a new concept: “stretchable PSHs,” i.e., helical spin patterns with continuously variable pitches $P$ over a wide parameter range. Stretching the PSH, i.e., gate controlling $P$ while staying locked in the PSH regime, provides protection from spin decay at the symmetry point $\alpha =\beta$, thus offering an important advantage over other methods. This protection is limited mainly by the cubic Dresselhaus term, which breaks the unidirectionality of the total SO field and causes spin decay at higher electron densities. We quantify the cubic term, and find it to be sufficiently weak so that the extracted spin-diffusion lengths and decay times show a significant enhancement near $\alpha =\beta$. Since within the continuous-locking regime quantum transport is diffusive (2D) for charge while ballistic (1D) for spin and thus amenable to coherent spin control, stretchable PSHs could provide the platform for the much heralded long-distance communication $\sim 8–25\mu \mathrm{m}$ between solid-state spin qubits, where the spin diffusion length for $\alpha \ne \beta$ is an order of magnitude smaller.

We study spinless electrons in a single-channel quantum wire interacting through attractive interaction, and the quantum Hall states that may be constructed by an array of such wires. For a single wire, the electrons may form two phases, the Luttinger liquid and the strongly paired phase. The Luttinger liquid is gapless to one- and two-electron excitations, while the strongly paired state is gapped to the former and gapless to the latter. In contrast to the case in which the wire is proximity coupled to an external superconductor, for an isolated wire there is no separate phase of a topological, weakly paired superconductor. Rather, this phase is adiabatically connected to the Luttinger liquid phase. The properties of the one-dimensional topological superconductor emerge when the number of channels in the wire becomes large. The quantum Hall states that may be formed by an array of single-channel wires depend on the Landau-level filling factors. For odd-denominator fillings $\nu =1/(2n+1)$, wires at the Luttinger phase form Laughlin states, while wires in the strongly paired phase form a bosonic fractional quantum Hall state of strongly bound pairs at a filling of $1/(8n+4)$. The transition between the two is of the universality class of Ising transitions in three dimensions. For even-denominator fractions $\nu =1/2n$, the two single-wire phases translate into four quantum Hall states. Two of those states are bosonic fractional quantum Hall states of weakly and strongly bound pairs of electrons. The other two are non-Abelian quantum Hall states, which originate from coupling wires close to their critical point. One of these non-Abelian states is the Moore-Read state. The transitions between all of these states are of the universality class of Majorana transitions. We point out some of the properties that characterize the different phases and the phase transitions.

A system of Majorana zero modes with random infinite-range interactions—the Sachdev-Ye-Kitaev (SYK) model—is thought to exhibit an intriguing relation to the horizons of extremal black holes in two-dimensional anti–de Sitter space. This connection provides a rare example of holographic duality between a solvable quantum-mechanical model and dilaton gravity. Here, we propose a physical realization of the SYK model in a solid-state system. The proposed setup employs the Fu-Kane superconductor realized at the interface between a three-dimensional topological insulator and an ordinary superconductor. The requisite $N$ Majorana zero modes are bound to a nanoscale hole fabricated in the superconductor that is threaded by $N$ quanta of magnetic flux. We show that when the system is tuned to the surface neutrality point (i.e., chemical potential coincident with the Dirac point of the topological insulator surface state) and the hole has sufficiently irregular shape, the Majorana zero modes are described by the SYK Hamiltonian. We perform extensive numerical simulations to demonstrate that the system indeed exhibits physical properties expected of the SYK model, including thermodynamic quantities and two-point as well as four-point correlators, and discuss ways in which these can be observed experimentally.