Yasir Iqbal

Fracton phases are a particularly exotic type of quantum spin liquids where the elementary quasiparticles are intrinsically immobile. These phases may be described by unconventional gauge theories known as tensor or multipolar gauge theories, characteristic for so-called type-I or type-II fracton phases, respectively. Both variants have been associated with distinctive singular patterns in the spin structure factor, such as multifold pinch points for type-I and quadratic pinch points for type-II fracton phases. Here, we assess the impact of quantum fluctuations on these patterns by numerically investigating the spin S=1/2 quantum version of a classical spin model on the octahedral lattice featuring exact realizations of multifold and quadratic pinch points, as well as an unusual pinch line singularity. Based on large scale pseudofermion and pseudo-Majorana functional renormalization group calculations, we take the intactness of these spectroscopic signatures as a measure for the stability of the corresponding fracton phases. We find that in all three cases, quantum fluctuations significantly modify the shape of pinch points or lines by smearing them out and shifting signal away from the singularities in contrast to effects of pure thermal fluctuations. This indicates possible fragility of these phases and allows us to identify characteristic fingerprints of their remnants.

We theoretically investigate the lowerature phase of the recently synthesized Lu2Mo2O5N2 material, an extraordinarily rare realization of a S=1/2 three-dimensional pyrochlore Heisenberg antiferromagnet in which Mo5+ are the S=1/2 magnetic species. Despite a Curie-Weiss temperature (ΘCW) of -121(1) K, experiments have found no signature of magnetic ordering or spin freezing down to T∗≈0.5 K. Using density functional theory, we find that the compound is well described by a Heisenberg model with exchange parameters up to third nearest neighbors. The analysis of this model via the pseudofermion functional renormalization group method reveals paramagnetic behavior down to a temperature of at least T=|ΘCW|/100, in agreement with the experimental findings hinting at a possible three-dimensional quantum spin liquid. The spin susceptibility profile in reciprocal space shows momentum-dependent features forming a "gearwheel" pattern, characterizing what may be viewed as a molten version of a chiral noncoplanar incommensurate spiral order under the action of quantum fluctuations. Our calculated reciprocal space susceptibility maps provide benchmarks for future neutron scattering experiments on single crystals of Lu2Mo2O5N2.

Pinch point singularities, associated with flat band magnetic excitations, are tell-tale signatures of Coulomb spin liquids. While their properties in the presence of quantum fluctuations have been widely studied, the fate of the complementary nonanalytic features - shaped as half moons and stars - arising from adjacent shallow dispersive bands has remained unexplored. Here, we address this question for the spin S=1/2 Heisenberg antiferromagnet on the kagome lattice with second and third neighbor couplings, which allows one to tune the classical ground state characterized by flat bands to one that is governed by shallow dispersive bands for intermediate coupling strengths. Employing the complementary strengths of variational Monte Carlo, pseudofermion functional renormalization group, and density-matrix renormalization group, we establish the quantum phase diagram of the model. The U(1) Dirac spin liquid ground state of the nearest-neighbor antiferromagnet remains remarkably robust till intermediate coupling strengths when it transitions into a pinwheel valence bond crystal displaying signatures of half moons in its structure factor. Our Letter thus identifies a microscopic setting that realizes one of the proximate orders of the Dirac spin liquid identified in a recent work [Song, Wang, Vishwanath, and He, Nat. Commun. 10, 4254 (2019)2041-172310.1038/s41467-019-11727-3]. For larger couplings, we obtain a collinear magnetically ordered ground state characterized by starlike patterns.

We investigate the nature of the ground state of the spin-12 Heisenberg antiferromagnet on the shuriken lattice by complementary state-of-the-art numerical techniques, such as variational Monte Carlo (VMC) with versatile Gutzwiller-projected Jastrow wave functions, unconstrained multivariable variational Monte Carlo (mVMC), and pseudofermion/pseudo-Majorana functional renormalization group (PFFRG/PMFRG) methods. We establish the presence of a quantum paramagnetic ground state and investigate its nature, by classifying symmetric and chiral quantum spin liquids, and inspecting their instabilities towards competing valence bond crystal (VBC) orders. Our VMC analysis reveals that a VBC with a pinwheel structure emerges as the lowest-energy variational ground state, and it is obtained as an instability of the U(1) Dirac spin liquid. Analogous conclusions are drawn from mVMC calculations employing accurate BCS pairing states supplemented by symmetry projectors, which confirm the presence of pinwheel VBC order by a thorough analysis of dimer-dimer correlation functions. Our work highlights the nontrivial role of quantum fluctuations via the Gutzwiller projector in resolving the subtle interplay between competing orders.

We investigate the ground-state properties of the spin-1/2 pyrochlore Heisenberg antiferromagnet using pseudofermion functional renormalization group techniques. The first part of our analysis is based on an enhanced parton mean-field approach, which takes into account fluctuation effects from renormalized vertex functions. Our implementation of this technique extends earlier approaches and resolves technical difficulties associated with a diagrammatic overcounting. Using various parton ansätze for quantum spin liquids, dimerized and nematic states our results indicate a tendency for lattice symmetry breaking in the ground state. While overall quantum spin liquids seem unfavorable in this system, the recently proposed monopole state still shows the strongest support among all spin liquid ansätze that we have tested, which is further confirmed by our complementary variational Monte Carlo calculations. In the second part of our investigation, we probe lattice symmetry breaking more directly by applying the pseudofermion functional renormalization group to perturbed systems. Our results from this technique confirm that the system's ground state either exhibits broken C3 rotation symmetry, or a combination of inversion and C3 symmetry breaking.

We investigate the quantum phases of the frustrated spin-12J1-J2-J3 Heisenberg model on the square lattice with ferromagnetic J1 and antiferromagnetic J2 and J3 interactions. Using the pseudofermion functional renormalization group technique, we find an intermediate paramagnetic phase located between classically ordered ferromagnetic, stripy antiferromagnetic, and incommensurate spiral phases. We observe that quantum fluctuations lead to significant shifts of the spiral pitch angles compared to the classical limit. By computing the response of the system with respect to various spin rotation and lattice symmetry-breaking perturbations, we identify a complex interplay between different nematic spin states in the paramagnetic phase. While retaining time-reversal invariance, these phases either break spin-rotation symmetry, lattice-rotation symmetry, or a combination of both. We therefore propose the J1-J2-J3 Heisenberg model on the square lattice as a paradigmatic example where different intimately connected types of nematic orders emerge in the same model.

We investigate the spin S=1/2 Heisenberg model on the body centered cubic lattice in the presence of ferromagnetic and antiferromagnetic nearest-neighbor J1, second-neighbor J2, and third-neighbor J3 exchange interactions. The classical ground state phase diagram obtained by a Luttinger-Tisza analysis is shown to host six different (noncollinear) helimagnetic orders in addition to ferromagnetic, Néel, stripe, and planar antiferromagnetic orders. Employing the pseudofermion functional renormalization group (PFFRG) method for quantum spins (S=1/2) we find an extended nonmagnetic region, and significant shifts to the classical phase boundaries and helimagnetic pitch vectors caused by quantum fluctuations, while no new long-range dipolar magnetic orders are stabilized. The nonmagnetic phase is found to disappear for S=1. We calculate the magnetic ordering temperatures from PFFRG and quantum Monte Carlo methods, and make comparisons to available data.

The quantum spin liquid is a highly entangled magnetic state characterized by the absence of static magnetism in its ground state. Instead, the spins fluctuate in a highly correlated way down to the lowest temperatures. Quantum spin liquids are very rare and are confined to a few specific cases where the interactions between the magnetic ions cannot be simultaneously satisfied (known as frustration). Lattices with magnetic ions in triangular or tetrahedral arrangements, which interact via isotropic antiferromagnetic interactions, can generate such a frustration. Three-dimensional isotropic spin liquids have mostly been sought in materials where the magnetic ions form pyrochlore or hyperkagome lattices. Here we present a three-dimensional lattice called the hyper-hyperkagome that enables spin liquid behaviour and manifests in the compound PbCuTe2O6. Using a combination of experiment and theory, we show that this system exhibits signs of being a quantum spin liquid with no detectable static magnetism together with the presence of diffuse continua in the magnetic spectrum suggestive of fractional spinon excitations.

We investigate the ground state and critical temperature (Tc) phase diagrams of the classical and quantum S=12 pyrochlore lattice with nearest-neighbor Heisenberg and Dzyaloshinskii-Moriya interactions (DMI). We consider ferromagnetic and antiferromagnetic Heisenberg exchange interaction as well as direct and indirect DMI. At the classical level, three ground states are found: all-in/all-out, ferromagnetic, and a locally ordered XY phase, known as Γ5, which displays an accidental classical U(1) degeneracy at the mean-field level. Quantum zero-point energy fluctuations computed to order 1/S are found to lift the classical ground-state degeneracy and select the so-called ψ3 state out of the degenerate manifold in most parts of the Γ5 regime. Likewise, thermal fluctuations treated classically at the Gaussian level entropically select the ψ3 state at T=0+. In contrast to this low-temperature state-selection behavior, classical Monte Carlo simulations find that the system orders at Tc in the noncoplanar ψ2 state of Γ5 for antiferromagnetic Heisenberg exchange and indirect DMI with a transition from ψ2 to ψ3 at a temperature TΓ5

We address the long-standing problem of the microscopic origin of the richly diverse phenomena in the chromium breathing pyrochlore material family. Combining electronic structure and renormalization group techniques we resolve the magnetic interactions and analyze their reciprocal-space susceptibility. We show that the physics of these materials is principally governed by long-range Heisenberg Hamiltonian interactions, a hitherto unappreciated fact. Our calculations uncover that in these isostructural compounds, the choice of chalcogen triggers a proximity of the materials to classical spin liquids featuring degenerate manifolds of wave-vectors of different dimensions: A Coulomb phase with three-dimensional degeneracy for LiInCr4O8 and LiGaCr4O8, a spiral spin liquid with two-dimensional degeneracy for CuInCr4Se8 and one-dimensional line degeneracies characteristic of the face-centered cubic antiferromagnet for LiInCr4S8, LiGaCr4S8, and CuInCr4S8. The surprisingly complex array of prototypical pyrochlore behaviors we discovered in chromium spinels may inspire studies of transition paths between different semi-classical spin liquids by doping or pressure.