Efficient simulations of epidemic models with tensor networks: Application to the one-dimensional susceptible-infected-susceptible model.

The contact process is an emblematic model of a nonequilibrium system, containing a phase transition between inactive and active dynamical regimes. In the epidemiological context, the model is known as the susceptible-infected-susceptible model, and it is widely used to describe contagious spreading. In this work, we demonstrate how accurate and efficient representations of the full probability distribution over all configurations of the contact process on a one-dimensional chain can be obtained by means of matrix product states (MPSs). We modify and adapt MPS methods from many-body quantum systems to study the classical distributions of the driven contact process at late times. We give accurate and efficient results for the distribution of large gaps, and illustrate the advantage of our methods over Monte Carlo simulations. Furthermore, we study the large deviation statistics of the dynamical activity, defined as the total number of configuration changes along a trajectory, and investigate quantum-inspired entropic measures, based on the second Rényi entropy.

Efficient Tensor Network Algorithm for Layered Systems.

Strongly correlated layered 2D systems are of central importance in condensed matter physics, but their numerical study is very challenging. Motivated by the enormous successes of tensor networks for 1D and 2D systems, we develop an efficient tensor network approach based on infinite projected entangled-pair states for layered 2D systems. Starting from an anisotropic 3D infinite projected entangled-pair state ansatz, we propose a contraction scheme in which the weakly interacting layers are effectively decoupled away from the center of the layers, such that they can be efficiently contracted using 2D contraction methods while keeping the center of the layers connected in order to capture the most relevant interlayer correlations. We present benchmark data for the anisotropic 3D Heisenberg model on a cubic lattice, which shows close agreement with quantum Monte Carlo and full 3D contraction results. Finally, we study the dimer to Néel phase transition in the Shastry-Sutherland model with interlayer coupling, a frustrated spin model that is out of reach of quantum Monte Carlo due to the negative sign problem.

Discovery of quantum phases in the Shastry-Sutherland compound SrCu2(BO3)2 under extreme conditions of field and pressure.

The 2-dimensional layered oxide material SrCu2(BO3)2, long studied as a realization of the Shastry-Sutherland spin topology, exhibits a range of intriguing physics as a function of both hydrostatic pressure and magnetic field, with a still debated intermediate plaquette phase appearing at approximately 20 kbar and a possible deconfined critical point at higher pressure. Here, we employ a tunnel diode oscillator (TDO) technique to probe the behavior in the combined extreme conditions of high pressure, high magnetic field, and low temperature. We reveal an extensive phase space consisting of multiple magnetic analogs of the elusive supersolid phase and a magnetization plateau. In particular, a 10 × 2 supersolid and a 1/5 plateau, identified by infinite Projected Entangled Pair States (iPEPS) calculations, are found to rely on the presence of both magnetic and non-magnetic particles in the sea of dimer singlets. These states are best understood as descendants of the full-plaquette phase, the leading candidate for the intermediate phase of SrCu2(BO3)2.

Emergent bound states and impurity pairs in chemically doped Shastry-Sutherland system.

Impurities often play a defining role in the ground states of frustrated quantum magnets. Studies of their effects are crucial in understanding of the phase diagram in these materials. SrCu2(BO3)2, an experimental realization of the Shastry-Sutherland (SS) lattice, provides a unique model system for such studies using both experimental and numerical approaches. Here we report effects of impurities on the crystals of bound states, and doping-induced emergent ground states in Mg-doped SrCu2(BO3)2, which remain stable in high magnetic fields. Using four complementary magnetometry techniques and theoretical simulations, a rich impurity-induced phenomenology at high fields is discovered. The results demonstrate a rare example in which even a small doping concentration interacts strongly with both triplets and bound states of triplets, and thus plays a significant role in the magnetization process even at high magnetic fields. Our findings provide insights into the study of impurity effects in geometrically frustrated quantum magnets.

Stripe order in the underdoped region of the two-dimensional Hubbard model.

Competing inhomogeneous orders are a central feature of correlated electron materials, including the high-temperature superconductors. The two-dimensional Hubbard model serves as the canonical microscopic physical model for such systems. Multiple orders have been proposed in the underdoped part of the phase diagram, which corresponds to a regime of maximum numerical difficulty. By combining the latest numerical methods in exhaustive simulations, we uncover the ordering in the underdoped ground state. We find a stripe order that has a highly compressible wavelength on an energy scale of a few kelvin, with wavelength fluctuations coupled to pairing order. The favored filled stripe order is different from that seen in real materials. Our results demonstrate the power of modern numerical methods to solve microscopic models, even in challenging settings.

Competing states in the t-J model: uniform D-wave state versus stripe state.

Variational studies of the t-J model on the square lattice based on infinite projected-entangled pair states confirm an extremely close competition between a uniform d-wave superconducting state and different stripe states. The site-centered stripe with an in-phase d-wave order has an equal or only slightly lower energy than the stripe with antiphase d-wave order. The optimal stripe filling is not constant but increases with J/t. A nematic anisotropy reduces the pairing amplitude and the energies of stripe phases are lowered relative to the uniform state with increasing nematicity.

Crystals of bound states in the magnetization plateaus of the Shastry-Sutherland model.

Using infinite projected entangled-pair states, we show that the Shastry-Sutherland model in an external magnetic field has low-magnetization plateaus which, in contrast to previous predictions, correspond to crystals of bound states of triplets, and not to crystals of triplets. The first sizable plateaus appear at magnetization 1/8, 2/15 and 1/6, in agreement with experiments on the orthogonal-dimer antiferromagnet SrCu2(BO3)2, and they can be naturally understood as regular patterns of bound states, including the intriguing 2/15 one. We also show that, even in a confined geometry, two triplets bind into a localized bound state with Sz=2. Finally, we discuss the role of competing domain-wall and supersolid phases, as well as that of additional anisotropic interactions.

Simultaneous dimerization and SU(4) symmetry breaking of 4-color fermions on the square lattice.

Using infinite projected entangled-pair states, exact diagonalization, and flavor-wave theory, we show that the SU(4) Heisenberg model undergoes a spontaneous dimerization on the square lattice, in contrast with its SU(2) and SU(3) counterparts, which develop Néel and three-sublattice stripelike long-range order. Since the ground state of a dimer is not a singlet for SU(4) but a 6-dimensional irreducible representation, this leaves the door open for further symmetry breaking. We provide evidence that, unlike in SU(4) ladders, where dimers pair up to form singlet plaquettes, here the SU(4) symmetry is additionally broken, leading to a gapless spectrum in spite of the broken translational symmetry.