-
Origin and stability of generalized Wigner crystallinity in triangular moiré systems
Authors:
Aman Kumar,
Cyprian Lewandowski,
Hitesh J. Changlani
Abstract:
Generalized Wigner crystals (GWC) on triangular moiré superlattices, formed from stacking two layers of transition metal chalcogenides, have been recently observed at multiple fractional fillings of the moiré unit cell [Nature 587, 214 - 218 (2020), Nat. Phys. 17, 715 - 719 (2021), Nature 597, 650 - 654 (2021)]. Motivated by these experiments, tied with the need for an accurate microscopic descrip…
▽ More
Generalized Wigner crystals (GWC) on triangular moiré superlattices, formed from stacking two layers of transition metal chalcogenides, have been recently observed at multiple fractional fillings of the moiré unit cell [Nature 587, 214 - 218 (2020), Nat. Phys. 17, 715 - 719 (2021), Nature 597, 650 - 654 (2021)]. Motivated by these experiments, tied with the need for an accurate microscopic description of these materials, we explore the theoretical origins of GWC at $n=1/3$ and $2/3$ filling. We demonstrate the limitations of theoretical descriptions of these moiré GWCs relying on truncated (finite-range) electron-electron interactions instead of a long-range Coulomb interaction. We validate our findings by studying both classical and quantum mechanical effects at zero and finite temperatures. More generally, we discuss the role of charge frustration in the theoretical extended Hubbard-model phase diagram, identifying a "pinball" phase, a partially quantum melted generalized Wigner crystal with coexisting solid and liquid-like features, with no classical analog. Quantum effects also explain the small, but experimentally detectable, asymmetry in the transition temperatures of the $n=1/3$ and $2/3$ crystals. Our calculations reveal that the charge ordering temperature and the magnetic crossover temperature predicted using the long-range Coulomb interaction are correctly captured with appropriately renormalized nearest-neighbor interactions. The effective nearest-neighbor interaction strength is significantly weaker than previously reported, placing it close to a metal-insulator phase boundary. This observation has implications for future experiments that we discuss. We conclude by studying the dependence of melting temperatures on gate-to-sample separation and we also predict temperature scales at which magnetic crossovers should be observed.
△ Less
Submitted 20 September, 2024;
originally announced September 2024.
-
Stable infinite-temperature eigenstates in SU(2)-symmetric nonintegrable models
Authors:
Christopher J. Turner,
Marcin Szyniszewski,
Bhaskar Mukherjee,
Ronald Melendrez,
Hitesh J. Changlani,
Arijeet Pal
Abstract:
Nonintegrable many-body quantum systems typically thermalize at long times through the mechanism of quantum chaos. However, some exceptional systems, such as those harboring quantum scars, break thermalization, serving as testbeds for foundational problems of quantum statistical physics. Here, we investigate a class of nonintegrable bond-staggered models that is endowed with a large number of zero…
▽ More
Nonintegrable many-body quantum systems typically thermalize at long times through the mechanism of quantum chaos. However, some exceptional systems, such as those harboring quantum scars, break thermalization, serving as testbeds for foundational problems of quantum statistical physics. Here, we investigate a class of nonintegrable bond-staggered models that is endowed with a large number of zero-energy eigenstates and possesses a non-Abelian internal symmetry. We use character theory to give a lower bound on the zero-energy degeneracy, which matches exact diagonalization results, and is found to grow exponentially with the system size. We also show that few-magnon zero-energy states have an exact analytical description, allowing us to build a basis of low-entangled fixed-separation states, which is stable to most perturbations found in experiments. This remarkable dynamical stability of special states elucidates our understanding of nonequilibrium processes in non-Abelian chaotic quantum models.
△ Less
Submitted 16 July, 2024;
originally announced July 2024.
-
Thermodynamics of the dipole-octupole pyrochlore magnet Ce$_2$Hf$_2$O$_{7}$ in applied magnetic fields
Authors:
Anish Bhardwaj,
Victor Porée,
Han Yan,
Nicolas Gauthier,
Elsa Lhotel,
Sylvain Petit,
Jeffrey A. Quilliam,
Andriy H. Nevidomskyy,
Romain Sibille,
Hitesh J. Changlani
Abstract:
The recently discovered dipole-octupole pyrochlore magnet Ce$_2$Hf$_2$O$_7$ is a promising three-dimensional quantum spin liquid candidate which shows no signs of ordering at low temperature. The low energy effective pseudospin-1/2 description in a magnetic field is characterized by the XYZ Hamiltonian and a Zeeman term where the dipolar local z-component of the pseudospin couples to the local z-c…
▽ More
The recently discovered dipole-octupole pyrochlore magnet Ce$_2$Hf$_2$O$_7$ is a promising three-dimensional quantum spin liquid candidate which shows no signs of ordering at low temperature. The low energy effective pseudospin-1/2 description in a magnetic field is characterized by the XYZ Hamiltonian and a Zeeman term where the dipolar local z-component of the pseudospin couples to the local z-component of the applied magnetic field, while the local x- and y-components of the pseudospin remain decoupled as a consequence of their octupolar character. Using effective parameters determined in V. Poree et al., arXiv:2305.08261 (2023), remarkable experimental features can be reproduced, as for instance the specific heat and magnetization data as well as the continuum of states seen in neutron scattering. Here we investigate the thermodynamic response to magnetic fields applied along the global [110] direction using specific heat measurements and fits using numerical methods, and solve the corresponding magnetic structure using neutron diffraction. Specific heat data in moderate fields are reproduced well, however, at high fields the agreement is not satisfactory. We especially observe a two-step release of entropy, a finding that demands a review of both theory and experiment. We address it within the framework of three possible scenarios, including an analysis of the crystal field Hamiltonian not restricted to the two-dimensional single-ion doublet subspace. We conclusively rule out two of these scenarios and find qualitative agreement with a simple model of field misalignment with respect to the crystalline direction. We discuss the implications of our findings for [111] applied fields and for future experiments on Ce$_2$Hf$_2$O$_7$ and its sister compounds.
△ Less
Submitted 13 February, 2024;
originally announced February 2024.
-
Locally controlled arrested thermalization
Authors:
Ken K. W. Ma,
Hitesh J. Changlani
Abstract:
The long-time dynamics of quantum systems, typically, but not always, results in a thermal steady state. The microscopic processes that lead to or circumvent this fate are of interest, since everyday experience tells us that not all spatial regions of a system heat up or cool down uniformly. This motivates the question: under what conditions can one slow down or completely arrest thermalization lo…
▽ More
The long-time dynamics of quantum systems, typically, but not always, results in a thermal steady state. The microscopic processes that lead to or circumvent this fate are of interest, since everyday experience tells us that not all spatial regions of a system heat up or cool down uniformly. This motivates the question: under what conditions can one slow down or completely arrest thermalization locally? Is it possible to construct realistic Hamiltonians and initial states such that a local region is effectively insulated from the rest, or acts like a barrier between two or more regions? We answer this in the affirmative by outlining the conditions that govern the flow of energy and entropy between subsystems. Using these ideas we provide a representative example for how simple few-body states can be used to engineer a ``thermal switch" between interacting regions.
△ Less
Submitted 9 May, 2024; v1 submitted 12 June, 2023;
originally announced June 2023.
-
Emergent strong zero mode through local Floquet engineering
Authors:
Bhaskar Mukherjee,
Ronald Melendrez,
Marcin Szyniszewski,
Hitesh J. Changlani,
Arijeet Pal
Abstract:
Periodically driven quantum systems host exotic phenomena which often do not have any analog in undriven systems. Floquet prethermalization and dynamical freezing of certain observables, via the emergence of conservation laws, are realized by controlling the drive frequency. These dynamical regimes can be leveraged to construct quantum memories and have potential applications in quantum informatio…
▽ More
Periodically driven quantum systems host exotic phenomena which often do not have any analog in undriven systems. Floquet prethermalization and dynamical freezing of certain observables, via the emergence of conservation laws, are realized by controlling the drive frequency. These dynamical regimes can be leveraged to construct quantum memories and have potential applications in quantum information processing. Solid state and cold atom experimental architectures have opened avenues for implementing local Floquet engineering which can achieve spatially modulated quantum control of states. Here, we uncover the novel memory effects of local periodic driving in a nonintegrable spin-half staggered Heisenberg chain. For a boundary-driven protocol at the dynamical freezing frequency, we show the formation of an approximate strong zero mode, a prethermal quasi-local operator, due to the emergence of a discrete global $\mathbb{Z}_2$ symmetry. This is captured by constructing an accurate effective Floquet Hamiltonian using a higher-order partially resummed Floquet-Magnus expansion. The lifetime of the boundary spin can be exponentially enhanced by enlarging the set of suitably chosen driven sites. We demonstrate that in the asymptotic limit, achieved by increasing the number of driven sites, a strong zero mode emerges, where the lifetime of the boundary spin grows exponentially with system size. The non-local processes in the Floquet Hamiltonian play a pivotal role in the total freezing of the boundary spin in the thermodynamic limit. The novel dynamics of the boundary spin is accompanied by a rich structure of entanglement in the Floquet eigenstates where specific bipartitions yield an area-law scaling while the entanglement for random bipartitions scales as a volume-law.
△ Less
Submitted 2 June, 2023;
originally announced June 2023.
-
Dipolar-octupolar correlations and hierarchy of exchange interactions in Ce$_2$Hf$_2$O$_7$
Authors:
Victor Porée,
Anish Bhardwaj,
Elsa Lhotel,
Sylvain Petit,
Nicolas Gauthier,
Han Yan,
Vladimir Pomjakushin,
Jacques Ollivier,
Jeffrey A. Quilliam,
Andriy H. Nevidomskyy,
Hitesh J. Changlani,
Romain Sibille
Abstract:
We investigate the correlated state of Ce$_2$Hf$_2$O$_7$ using neutron scattering, finding signatures of correlations of both dipolar and octupolar character. A dipolar inelastic signal is also observed, as expected for spinons in a quantum spin ice (QSI). Fits of thermodynamic data using exact diagonalization methods indicate that the largest interaction is an octupolar exchange, with a strength…
▽ More
We investigate the correlated state of Ce$_2$Hf$_2$O$_7$ using neutron scattering, finding signatures of correlations of both dipolar and octupolar character. A dipolar inelastic signal is also observed, as expected for spinons in a quantum spin ice (QSI). Fits of thermodynamic data using exact diagonalization methods indicate that the largest interaction is an octupolar exchange, with a strength roughly twice as large as other terms. A hierarchy of exchange interactions with dominant octupolar and significant dipolar exchange, still in the octupolar QSI phase, rationalises neutron scattering observations. Our results reveal a `quantum multipolar liquid' where correlations involve multiple terms in moment series expansion, opening questions about their intertwining and possible hierarchy.
△ Less
Submitted 20 February, 2024; v1 submitted 14 May, 2023;
originally announced May 2023.
-
Real space thermalization of locally driven quantum magnets
Authors:
Ronald Melendrez,
Bhaskar Mukherjee,
Prakash Sharma,
Arijeet Pal,
Hitesh J. Changlani
Abstract:
The study of thermalization and its breakdown in isolated systems has led to a deeper understanding of non-equilibrium quantum states and their dependence on initial conditions. The role of initial conditions is prominently highlighted by the existence of quantum many-body scars, special athermal states with an underlying effective superspin structure, embedded in an otherwise chaotic many-body sp…
▽ More
The study of thermalization and its breakdown in isolated systems has led to a deeper understanding of non-equilibrium quantum states and their dependence on initial conditions. The role of initial conditions is prominently highlighted by the existence of quantum many-body scars, special athermal states with an underlying effective superspin structure, embedded in an otherwise chaotic many-body spectrum. Spin Heisenberg and $XXZ$ models and their variants in one and higher dimension have been shown to host exact quantum many-body scars, exhibiting perfect revivals of spin helix states that are realizable in synthetic and condensed matter systems. Motivated by these advances, we propose experimentally accessible, local, time-dependent protocols to explore the spatial thermalization profile and highlight how different parts of the system thermalize and affect the fate of the superspin. We identify distinct parametric regimes for the ferromagnetic ($X$-polarized) initial state based on the interplay between the driven spin and the rest, including local athermal behavior where the driven spin effectively decouples, acting like a ``cold" spot while being instrumental in heating up the other spins. We also identify parameter regimes where the superspin remains resilient to local driving for long time scales. We develop a real and Floquet space picture that explains our numerical observations, and make predictions that can be tested in various experimental setups.
△ Less
Submitted 26 September, 2024; v1 submitted 28 December, 2022;
originally announced December 2022.
-
Triangular lattice Hubbard model physics at intermediate temperatures
Authors:
Kyungmin Lee,
Prakash Sharma,
Oskar Vafek,
Hitesh J. Changlani
Abstract:
Moire systems offer an exciting playground to study many-body effects of strongly correlated electrons in regimes that are not easily accessible in conventional material settings. Motivated by a recent experiment on $\text{WSe}_2/\text{WS}_2$ moire bilayers [Y. Tang et al., Nature 579, 353-358 (2020)], which realizes a triangular superlattice with a small hopping (of approximately 10 Kelvin), with…
▽ More
Moire systems offer an exciting playground to study many-body effects of strongly correlated electrons in regimes that are not easily accessible in conventional material settings. Motivated by a recent experiment on $\text{WSe}_2/\text{WS}_2$ moire bilayers [Y. Tang et al., Nature 579, 353-358 (2020)], which realizes a triangular superlattice with a small hopping (of approximately 10 Kelvin), with tunable density of holes, we explore the Hubbard model on the triangular lattice for intermediate temperatures $t \lesssim T < U$. Employing finite temperature Lanczos calculations, and closely following the fitting protocols used in the experiment, we recover the observed trends in the reported Curie-Weiss temperature $Θ$ with filling, using the reported interaction strength $U/t=20$. We focus on the large increase of $|Θ|$ on decreasing the density below half filling and the sign change of $Θ$ at higher fillings, which signals the onset of ferromagnetism. The increase in $|Θ|$ is also seen in the $t$-$J$ model (the low energy limit of the Hubbard model) in the intermediate temperature range, which we clarify is opposite to the trend in its high temperature limit. Differences between the low, intermediate and high temperature regimes are discussed. Our numerical calculations also capture the crossover between short-range antiferromagnetic to ferromagnetic order in the intermediate temperature regime, a result broadly consistent with the experimental findings. We find that this behavior is a finite-temperature remnant of the underlying zero temperature phase transition, which we explore with ground state density matrix renormalization group calculations. We provide evidence of ferromagnetism characterized by weak (but robust) correlations that explain the small $Θ$ seen in the experiment.
△ Less
Submitted 5 June, 2023; v1 submitted 1 September, 2022;
originally announced September 2022.
-
Enigma of the vortex state in a strongly correlated d-wave superconductor
Authors:
Anushree Datta,
Hitesh J. Changlani,
Kun Yang,
Amit Ghosal
Abstract:
We show that strong electronic repulsion transforms a vortex core from a metallic-type in overdoped regime to a Mott-insulator at underdoping of a strongly correlated d-wave superconductor. This changeover is accompanied by an accumulation of electron density at the vortex core towards local half-filling in the underdoped region, which in turn facilitates the formation of the Mott insulating core.…
▽ More
We show that strong electronic repulsion transforms a vortex core from a metallic-type in overdoped regime to a Mott-insulator at underdoping of a strongly correlated d-wave superconductor. This changeover is accompanied by an accumulation of electron density at the vortex core towards local half-filling in the underdoped region, which in turn facilitates the formation of the Mott insulating core. We find that the size of vortices evolves non-monotonically with doping. A similar non-monotonicity of critical field $H_{\rm c2}$, as extracted from superfluid stiffness, is also found. Our results explain some recent experimental puzzles of cuprate superconductors.
△ Less
Submitted 2 May, 2023; v1 submitted 15 June, 2022;
originally announced June 2022.
-
Response to comment on "Spin-lattice coupling and the emergence of the trimerized phase in the $S=1$ kagome antiferromagnet Na$_2$Ti$_3$Cl$_8$", Phys. Rev. Lett. 124, 167203 (2020)
Authors:
Arpita Paul,
Chia-Min Chung,
Turan Birol,
Hitesh J. Changlani
Abstract:
In the context of the $S=1$ kagome antiferromagnet Na$_2$Ti$_3$Cl$_8$, we respond to the comment by Khomskii et al. [D.I. Khomskii, T. Mizokawa and S.V. Streltsov, Phys. Rev. Lett. 127, 049701 (2021)] on previous work by Paul et al. [A. Paul, C.-M. Chung, T. Birol, and H. J. Changlani, Phys. Rev. Lett. 124, 167203 (2020)].
In the context of the $S=1$ kagome antiferromagnet Na$_2$Ti$_3$Cl$_8$, we respond to the comment by Khomskii et al. [D.I. Khomskii, T. Mizokawa and S.V. Streltsov, Phys. Rev. Lett. 127, 049701 (2021)] on previous work by Paul et al. [A. Paul, C.-M. Chung, T. Birol, and H. J. Changlani, Phys. Rev. Lett. 124, 167203 (2020)].
△ Less
Submitted 5 August, 2021;
originally announced August 2021.
-
Sleuthing out exotic quantum spin liquidity in the pyrochlore magnet Ce$_2$Zr$_2$O$_7$
Authors:
Anish Bhardwaj,
Shu Zhang,
Han Yan,
Roderich Moessner,
Andriy H. Nevidomskyy,
Hitesh J. Changlani
Abstract:
The search for quantum spin liquids (QSL) -- topological magnets with fractionalized excitations -- has been a central theme in condensed matter and materials physics. While theories are no longer in short supply, tracking down materials has turned out to be remarkably tricky, in large part because of the difficulty to diagnose experimentally a state with only topological, rather than conventional…
▽ More
The search for quantum spin liquids (QSL) -- topological magnets with fractionalized excitations -- has been a central theme in condensed matter and materials physics. While theories are no longer in short supply, tracking down materials has turned out to be remarkably tricky, in large part because of the difficulty to diagnose experimentally a state with only topological, rather than conventional, forms of order. Pyrochlore systems have proven particularly promising, hosting a classical Coulomb phase in the spin ices Dy/Ho$_2$Ti$_2$O$_7$, with subsequent proposals of candidate QSLs in other pyrochlores. Connecting experiment with detailed theory exhibiting a robust QSL has remained a central challenge. Here, focusing on the strongly spin-orbit coupled effective $S=1/2$ pyrochlore Ce$_2$Zr$_2$O$_7$, we analyse recent thermodynamic and neutron scattering experiments, to identify a microscopic effective Hamiltonian through a combination of finite temperature Lanczos, Monte Carlo and analytical spin dynamics calculations. Its parameter values suggest a previously unobserved exotic phase, a $π$-flux U(1) QSL. Intriguingly, the octupolar nature of the moments makes them less prone to be affected by crystal imperfections or magnetic impurities, while also hiding some otherwise characteristic signatures from neutrons, making this QSL arguably more stable than its more conventional counterparts.
△ Less
Submitted 2 August, 2021;
originally announced August 2021.
-
Multimagnon dynamics and thermalization in the $S=1$ easy-axis ferromagnetic chain
Authors:
Prakash Sharma,
Kyungmin Lee,
Hitesh J. Changlani
Abstract:
Quasiparticles are physically motivated mathematical constructs for simplifying the seemingly complicated many-body description of solids. A complete understanding of their dynamics and the nature of the effective interactions between them provides rich information on real material properties at the microscopic level. In this work, we explore the dynamics and interactions of magnon quasiparticles…
▽ More
Quasiparticles are physically motivated mathematical constructs for simplifying the seemingly complicated many-body description of solids. A complete understanding of their dynamics and the nature of the effective interactions between them provides rich information on real material properties at the microscopic level. In this work, we explore the dynamics and interactions of magnon quasiparticles in a ferromagnetic spin-1 Heisenberg chain with easy-axis onsite anisotropy, a model relevant for the explanation of recent terahertz optics experiments on NiNb$_2$O$_6$ [P. Chauhan et al., Phys. Rev. Lett. 124, 037203 (2020)],and nonequilibrium dynamics in ultracold atomic settings [W.C. Chung et al., Phys. Rev. Lett. 126, 163203 (2021)]. We build a picture for the properties of clouds of a few magnons with the help of exact diagonalization and density matrix renormalization group calculations supported by physically motivated Jastrow wavefunctions. We show how the binding energy of magnons effectively reduces with their number and explain how this energy scale is of direct relevance for dynamical magnetic susceptibility measurements. This understanding is used to make predictions for ultracold-atomic platforms which are ideally suited to study the thermalization of multimagnon states. We simulate the nonequilibrium dynamics of these chains using the matrix product state based time-evolution block decimation algorithm and explore the dependence of revivals and thermalization on magnon density and easy-axis onsite anisotropy (which controls the strength of effective magnon interactions). We observe behaviors akin to those reported for many-body quantum scars which we explain with an analytic approximation that is accurate in the limit of small anisotropy.
△ Less
Submitted 1 March, 2022; v1 submitted 19 July, 2021;
originally announced July 2021.
-
Many-body energy invariant for $T$-linear resistivity
Authors:
Aavishkar A. Patel,
Hitesh J. Changlani
Abstract:
The description of dynamics of strongly correlated quantum matter is a challenge, particularly in physical situations where a quasiparticle description is absent. In such situations, however, the many-body Kubo formula from linear response theory, involving matrix elements of the current operator computed with many-body wavefunctions, remains valid. Working directly in the many-body Hilbert space…
▽ More
The description of dynamics of strongly correlated quantum matter is a challenge, particularly in physical situations where a quasiparticle description is absent. In such situations, however, the many-body Kubo formula from linear response theory, involving matrix elements of the current operator computed with many-body wavefunctions, remains valid. Working directly in the many-body Hilbert space and not making any reference to quasiparticles (or lack thereof), we address the puzzle of linear in temperature ($T$-linear) resistivity seen in non-Fermi liquid phases that occur in several strongly correlated condensed matter systems. We derive a simple criterion for the occurrence of $T$-linear resistivity based on an analysis of the contributions to the many-body Kubo formula, determined by an energy invariant "$f$-function" involving current matrix elements and energy eigenvalues that describes the DC conductivity of the system in the microcanonical ensemble. Using full diagonalization, we test this criterion for the $f$-function in the spinless nearest neighbor Hubbard model and in a system of Sachdev-Ye-Kitaev dots coupled by weak single particle hopping. We also study the $f$-function for the spin conductivity in the 2D Heisenberg model and arrive at similar conclusions. Our work suggests that a general principle, formulated in terms of many-body Hilbert space concepts, is at the core of the occurrence of $T$-linear resistivity in a wide range of systems, and precisely translates $T$-linear resistivity into a notion of energy scale invariance far beyond what is typically associated with quantum critical points.
△ Less
Submitted 1 March, 2022; v1 submitted 2 June, 2021;
originally announced June 2021.
-
Investigation of the monopole magneto-chemical potential in spin ices using capacitive torque magnetometry
Authors:
Naween Anand,
Kevin Barry,
Jennifer N. Neu,
David E. Graf,
Qing Huang,
Haidong Zhou,
Theo Siegrist,
Hitesh J. Changlani,
Christianne Beekman
Abstract:
The single-ion anisotropy and magnetic interactions in spin-ice systems give rise to unusual non-collinear spin textures, such as Pauling states and magnetic monopoles. The effective spin correlation strength ($J_{eff}$) determines the relative energies of the different spin-ice states. With this work, we display the capability of capacitive torque magnetometry in characterizing the magneto-chemic…
▽ More
The single-ion anisotropy and magnetic interactions in spin-ice systems give rise to unusual non-collinear spin textures, such as Pauling states and magnetic monopoles. The effective spin correlation strength ($J_{eff}$) determines the relative energies of the different spin-ice states. With this work, we display the capability of capacitive torque magnetometry in characterizing the magneto-chemical potential associated with monopole formation. We build a magnetic phase diagram of Ho$_2$Ti$_2$O$_7$, and show that the magneto-chemical potential depends on the spin-sublattice ($α$ or $β$), i.e., the Pauling state, involved in the transition. Monte-Carlo simulations using the dipolar-spin-ice Hamiltonian support our findings of a sublattice-dependent magneto-chemical potential, but the model underestimates the $J_{eff}$ for the $β$-sublattice. Additional simulations, including next-nearest neighbor interactions ($J_2$), show that long-range exchange terms in the Hamiltonian are needed to describe the measurements. This demonstrates that torque magnetometry provides a sensitive test for $J_{eff}$ and the spin-spin interactions that contribute to it.
△ Less
Submitted 5 July, 2022; v1 submitted 12 May, 2021;
originally announced May 2021.
-
Colorful points in the XY regime of XXZ quantum magnets
Authors:
Santanu Pal,
Prakash Sharma,
Hitesh J. Changlani,
Sumiran Pujari
Abstract:
In the XY regime of the XXZ Heisenberg model phase diagram, we demonstrate that the origin of magnetically ordered phases is influenced by the presence of solvable points with exact quantum coloring ground states featuring a quantum-classical correspondence. Using exact diagonalization and density matrix renormalization group calculations, for both the square and the triangular lattice magnets, we…
▽ More
In the XY regime of the XXZ Heisenberg model phase diagram, we demonstrate that the origin of magnetically ordered phases is influenced by the presence of solvable points with exact quantum coloring ground states featuring a quantum-classical correspondence. Using exact diagonalization and density matrix renormalization group calculations, for both the square and the triangular lattice magnets, we show that the ordered physics of the solvable points in the extreme XY regime, at $\frac{J_z}{J_\perp}=-1$ and $\frac{J_z}{J_\perp}=-\frac{1}{2}$ respectively with $J_\perp > 0$, adiabatically extends to the more isotropic regime $\frac{J_z}{J_\perp} \sim 1$. We highlight the projective structure of the coloring ground states to compute the correlators in fixed magnetization sectors which enables an understanding of the features in the static spin structure factors and correlation ratios. These findings are contrasted with an anisotropic generalization of the celebrated one-dimensional Majumdar-Ghosh model, which is also found to be (ground state) solvable. For this model, both exact dimer and three-coloring ground states exist at $\frac{J_z}{J_\perp}=-\frac{1}{2}$ but only the two dimer ground states survive for any $\frac{J_z}{J_\perp} > -\frac{1}{2}$.
△ Less
Submitted 12 April, 2021; v1 submitted 10 December, 2020;
originally announced December 2020.
-
Frustration-induced emergent Hilbert space fragmentation
Authors:
Kyungmin Lee,
Arijeet Pal,
Hitesh J. Changlani
Abstract:
Although most quantum systems thermalize locally on short time scales independent of initial conditions, recent developments have shown this is not always the case. Lattice geometry and quantum mechanics can conspire to produce constrained quantum dynamics and associated glassy behavior, a phenomenon that falls outside the rubric of the eigenstate thermalization hypothesis. Constraints "fragment"…
▽ More
Although most quantum systems thermalize locally on short time scales independent of initial conditions, recent developments have shown this is not always the case. Lattice geometry and quantum mechanics can conspire to produce constrained quantum dynamics and associated glassy behavior, a phenomenon that falls outside the rubric of the eigenstate thermalization hypothesis. Constraints "fragment" the many-body Hilbert space due to which some states remain insulated from others and the system fails to attain thermal equilibrium. Such fragmentation is a hallmark of geometrically frustrated magnets with low-energy "icelike manifolds" exhibiting a broad range of relaxation times for different initial states. Focusing on the highly frustrated kagome lattice, we demonstrate these phenomena in the Balents-Fisher-Girvin Hamiltonian (easy-axis regime), and a three-coloring model (easy-plane regime), both with constrained Hilbert spaces. We study their level statistics and relaxation dynamics to develop a coherent picture of fragmentation in various limits of the XXZ model on the kagome lattice.
△ Less
Submitted 16 June, 2021; v1 submitted 3 November, 2020;
originally announced November 2020.
-
The Future of the Correlated Electron Problem
Authors:
A. Alexandradinata,
N. P. Armitage,
Andrey Baydin,
Wenli Bi,
Yue Cao,
Hitesh J. Changlani,
Eli Chertkov,
Eduardo H. da Silva Neto,
Luca Delacretaz,
Ismail El Baggari,
G. M. Ferguson,
William J. Gannon,
Sayed Ali Akbar Ghorashi,
Berit H. Goodge,
Olga Goulko,
G. Grissonnanche,
Alannah Hallas,
Ian M. Hayes,
Yu He,
Edwin W. Huang,
Anshul Kogar,
Divine Kumah,
Jong Yeon Lee,
A. Legros,
Fahad Mahmood
, et al. (22 additional authors not shown)
Abstract:
A central problem in modern condensed matter physics is the understanding of materials with strong electron correlations. Despite extensive work, the essential physics of many of these systems is not understood and there is very little ability to make predictions in this class of materials. In this manuscript we share our personal views on the major open problems in the field of correlated electro…
▽ More
A central problem in modern condensed matter physics is the understanding of materials with strong electron correlations. Despite extensive work, the essential physics of many of these systems is not understood and there is very little ability to make predictions in this class of materials. In this manuscript we share our personal views on the major open problems in the field of correlated electron systems. We discuss some possible routes to make progress in this rich and fascinating field. This manuscript is the result of the vigorous discussions and deliberations that took place at Johns Hopkins University during a three-day workshop January 27, 28, and 29, 2020 that brought together six senior scientists and 46 more junior scientists. Our hope, is that the topics we have presented will provide inspiration for others working in this field and motivation for the idea that significant progress can be made on very hard problems if we focus our collective energies.
△ Less
Submitted 13 July, 2022; v1 submitted 1 October, 2020;
originally announced October 2020.
-
Exact three-colored quantum scars from geometric frustration
Authors:
Kyungmin Lee,
Ronald Melendrez,
Arijeet Pal,
Hitesh J. Changlani
Abstract:
Non-equilibrium properties of quantum materials present many intriguing properties, among them athermal behavior, which violates the eigenstate thermalization hypothesis. Such behavior has primarily been observed in disordered systems. More recently, experimental and theoretical evidence for athermal eigenstates, known as "quantum scars" has emerged in non-integrable disorder-free models in one di…
▽ More
Non-equilibrium properties of quantum materials present many intriguing properties, among them athermal behavior, which violates the eigenstate thermalization hypothesis. Such behavior has primarily been observed in disordered systems. More recently, experimental and theoretical evidence for athermal eigenstates, known as "quantum scars" has emerged in non-integrable disorder-free models in one dimension with constrained dynamics. In this work, we show the existence of quantum scar eigenstates and investigate their dynamical properties in many simple two-body Hamiltonians with "staggered" interactions, involving ferromagnetic and antiferromagnetic motifs, in arbitrary dimensions. These magnetic models include simple modifications of widely studied ones (e.g., the XXZ model) on a variety of frustrated and unfrustrated lattices. We demonstrate our ideas by focusing on the two dimensional frustrated spin-1/2 kagome antiferromagnet, which was previously shown to harbor a special exactly solvable point with "three-coloring" ground states in its phase diagram. For appropriately chosen initial product states -- for example, those which correspond to any state of valid three-colors -- we show the presence of robust quantum revivals, which survive the addition of anisotropic terms. We also suggest avenues for future experiments which may see this effect in real materials.
△ Less
Submitted 20 February, 2020;
originally announced February 2020.
-
Orientation Dependence of the Magnetic Phase Diagram of Yb$_2$Ti$_2$O$_7$
Authors:
S. Säubert,
A. Scheie,
C. Duvinage,
J. Kindervater,
S. Zhang,
H. J. Changlani,
Guangyong Xu,
S. M. Koohpayeh,
O. Tchernyshyov,
C. L. Broholm,
C. Pfleiderer
Abstract:
In the quest to realize a quantum spin liquid (QSL), magnetic long-range order is hardly welcome. Yet it can offer deep insights into a complex world of strong correlations and fluctuations. Much hope was placed in the cubic pyrochlore Yb$_2$Ti$_2$O$_7$ as a putative U(1) QSL but a new class of ultra-pure single crystals make it abundantly clear the stoichiometric compound is a ferromagnet. Here w…
▽ More
In the quest to realize a quantum spin liquid (QSL), magnetic long-range order is hardly welcome. Yet it can offer deep insights into a complex world of strong correlations and fluctuations. Much hope was placed in the cubic pyrochlore Yb$_2$Ti$_2$O$_7$ as a putative U(1) QSL but a new class of ultra-pure single crystals make it abundantly clear the stoichiometric compound is a ferromagnet. Here we present a detailed experimental and theoretical study of the corresponding field-temperature phase diagram. We find it to be richly anisotropic with a critical endpoint for $\vec{B}\,\parallel\,\langle 100\rangle$, while field parallel to $\langle 110 \rangle$ and $\langle 111 \rangle$ enhances the critical temperature by up to a factor of two and shifts the onset of the field-polarized state to finite fields. Landau theory shows that Yb$_2$Ti$_2$O$_7$ in some ways is remarkably similar to pure iron. However, it also pinpoints anomalies that cannot be accounted for at the classical mean-field level including a dramatic enhancement of $T_{\mathrm{C}}$ and reentrant phase boundary by fields with a component transverse to the easy axes, as well as the anisotropy of the upper critical field in the quantum limit.
△ Less
Submitted 23 December, 2019;
originally announced December 2019.
-
Multiphase Magnetism in Yb2Ti2O7
Authors:
A. Scheie,
J. Kindervater,
S. Zhang,
H. J. Changlani,
G. Sala,
G. Ehlers,
A. Heinemann,
G. S. Tucker,
S. M. Koohpayeh,
C. Broholm
Abstract:
We document the coexistence of ferro- and anti-ferromagnetism in pyrochlore $\rm Yb_2Ti_2O_7$ using three neutron scattering techniques on stoichiometric crystals: elastic neutron scattering shows a canted ferromagnetic ground state, neutron scattering shows spin wave excitations from both a ferro-and an antiferro-magnetic state, and field and temperature dependent small angle neutron scattering r…
▽ More
We document the coexistence of ferro- and anti-ferromagnetism in pyrochlore $\rm Yb_2Ti_2O_7$ using three neutron scattering techniques on stoichiometric crystals: elastic neutron scattering shows a canted ferromagnetic ground state, neutron scattering shows spin wave excitations from both a ferro-and an antiferro-magnetic state, and field and temperature dependent small angle neutron scattering reveals the corresponding anisotropic magnetic domain structure. High-field $\langle 111 \rangle$ spin wave fits show that $\rm Yb_2Ti_2O_7$ is extremely close to an antiferromagnetic phase boundary. Classical Monte Carlo simulations based on the interactions inferrred from high field spin wave measurements confirm $ψ_2$ antiferromagnetism is metastable within the FM ground state.
△ Less
Submitted 10 December, 2019;
originally announced December 2019.
-
Spin--lattice coupling and the emergence of the trimerized phase in the $S=1$ Kagome antiferromagnet Na$_2$Ti$_3$Cl$_8$
Authors:
Arpita Paul,
Chia-Min Chung,
Turan Birol,
Hitesh J. Changlani
Abstract:
Spin-1 antiferromagnets are abundant in nature, but few theories or results exist to understand their general properties and behavior, particularly in situations when geometric frustration is present. Here we study the $S=1$ Kagome compound Na$_2$Ti$_3$Cl$_8$ using a combination of Density Functional Theory, Exact Diagonalization, and Density Matrix Renormalization Group methods to achieve a first…
▽ More
Spin-1 antiferromagnets are abundant in nature, but few theories or results exist to understand their general properties and behavior, particularly in situations when geometric frustration is present. Here we study the $S=1$ Kagome compound Na$_2$Ti$_3$Cl$_8$ using a combination of Density Functional Theory, Exact Diagonalization, and Density Matrix Renormalization Group methods to achieve a first principles supported explanation of exotic magnetic phases in this compound. We find that the effective magnetic Hamiltonian includes essential non-Heisenberg terms that do not stem from spin-orbit coupling, and both trimerized and spin-nematic magnetic phases are relevant. The experimentally observed structural transition to a breathing Kagome phase is driven by spin--lattice coupling, which favors the trimerized magnetic phase against the quadrupolar one. We thus show that lattice effects can be necessary to understand the magnetism in frustrated magnetic compounds, and surmise that Na$_2$Ti$_3$Cl$_8$ is a compound which cannot be understood from only electronic or only lattice Hamiltonians, very much like VO$_2$.
△ Less
Submitted 4 September, 2019;
originally announced September 2019.
-
Probing magnon dynamics and interactions in a ferromagnetic spin-1 chain
Authors:
Prashant Chauhan,
Fahad Mahmood,
Hitesh J. Changlani,
S. M. Koohpayeh,
N. P. Armitage
Abstract:
NiNb$_{2}$O$_{6}$ is an almost ideal realization of a 1D spin-1 ferromagnetic Heisenberg chain compound with weak unidirectional anisotropy. Using time-domain THz spectroscopy, we measure the low-energy electrodynamic response of NiNb$_{2}$O$_{6}$ as a function of temperature and external magnetic field. At low temperatures, we find a magnon-like spin-excitation, which corresponds to the lowest en…
▽ More
NiNb$_{2}$O$_{6}$ is an almost ideal realization of a 1D spin-1 ferromagnetic Heisenberg chain compound with weak unidirectional anisotropy. Using time-domain THz spectroscopy, we measure the low-energy electrodynamic response of NiNb$_{2}$O$_{6}$ as a function of temperature and external magnetic field. At low temperatures, we find a magnon-like spin-excitation, which corresponds to the lowest energy excitation at $q\sim0$. At higher temperatures, we unexpectedly observe a temperature-dependent renormalization of the spin-excitation energy, which has a strong dependence on field direction. Using theoretical arguments, exact diagonalizations and finite temperature dynamical Lanczos calculations, we construct a picture of magnon-magnon interactions that naturally explains the observed renormalization. This unique scenario is a consequence of the spin-1 nature and has no analog in the more widely studied spin-1/2 systems.
△ Less
Submitted 15 July, 2019;
originally announced July 2019.
-
Dynamical structure factor of the three-dimensional quantum spin liquid candidate NaCaNi$_2$F$_7$
Authors:
Shu Zhang,
Hitesh J. Changlani,
Kemp W. Plumb,
Oleg Tchernyshyov,
Roderich Moessner
Abstract:
We study the spin-1 pyrochlore material NaCaNi$_2$F$_7$ with a combination of molecular dynamics simulations, stochastic dynamical theory and linear spin wave theory. The dynamical structure factor from inelastic neutron scattering is well described with a near-ideal Heisenberg Hamiltonian incorporating small anisotropic terms {and weak second-neighbor interactions}. We find that all three approac…
▽ More
We study the spin-1 pyrochlore material NaCaNi$_2$F$_7$ with a combination of molecular dynamics simulations, stochastic dynamical theory and linear spin wave theory. The dynamical structure factor from inelastic neutron scattering is well described with a near-ideal Heisenberg Hamiltonian incorporating small anisotropic terms {and weak second-neighbor interactions}. We find that all three approaches reproduce remarkably well the momentum dependence of the scattering intensity as well as its energy dependence with the exception of the lowest energies. These results are notable in that (i) the data show a complete lack of sharp quasiparticle excitations in momentum space over much, if not all, of the energy range; (ii) linear spin-wave theory appears to apply in a regime where it would be expected to fail for a number of reasons. We elucidate what underpins these surprises, and note that basic questions about the nature of quantum spin liquidity in such systems pose themselves as a result.
△ Less
Submitted 22 October, 2018;
originally announced October 2018.
-
Resonating quantum three-coloring wavefunctions for the kagome quantum antiferromagnet
Authors:
Hitesh J. Changlani,
Sumiran Pujari,
Chia-Min Chung,
Bryan K. Clark
Abstract:
Motivated by the recent discovery of a macroscopically degenerate exactly solvable point of the spin-$1/2$ $XXZ$ model for $J_z/J=-1/2$ on the kagome lattice [H. J. Changlani et al. Phys. Rev. Lett 120, 117202 (2018)] -- a result that holds for arbitrary magnetization -- we develop an exact mapping between its exact quantum three-coloring wavefunctions and the characteristic localized and topologi…
▽ More
Motivated by the recent discovery of a macroscopically degenerate exactly solvable point of the spin-$1/2$ $XXZ$ model for $J_z/J=-1/2$ on the kagome lattice [H. J. Changlani et al. Phys. Rev. Lett 120, 117202 (2018)] -- a result that holds for arbitrary magnetization -- we develop an exact mapping between its exact quantum three-coloring wavefunctions and the characteristic localized and topological magnons. This map, involving resonating two-color loops, is developed to represent exact many-body ground state wavefunctions for special high magnetizations. Using this map we show that these exact ground state solutions are valid for any $J_z/J \geq -1/2$. This demonstrates the equivalence of the ground-state wavefunction of the Ising, Heisenberg and $XY$ regimes all the way to the $J_z/J=-1/2$ point for these high magnetization sectors. In the hardcore bosonic language, this means that a certain class of exact many-body solutions, previously argued to hold for purely repulsive interactions ($J_z \geq 0$), actually hold for attractive interactions as well, up to a critical interaction strength. For the case of zero magnetization, where the ground state is not exactly known, we perform density matrix renormalization group calculations. Based on the calculation of the ground state energy and measurement of order parameters, we provide evidence for a lack of any qualitative change in the ground state on finite clusters in the Ising ($J_z \gg J$), Heisenberg ($J_z=J$) and $XY$ ($J_z=0$) regimes, continuing adiabatically to the vicinity of the macroscopically degenerate $J_z/J=-1/2$ point. These findings offer a framework for recent results in the literature, and also suggest that the $J_z/J=-1/2$ point is an unconventional quantum critical point whose vicinity may contain the key to resolving the spin-$1/2$ kagome problem.
△ Less
Submitted 26 August, 2018;
originally announced August 2018.
-
From real materials to model Hamiltonians with density matrix downfolding
Authors:
Huihuo Zheng,
Hitesh J. Changlani,
Kiel T. Williams,
Brian Busemeyer,
Lucas K. Wagner
Abstract:
Due to advances in computer hardware and new algorithms, it is now possible to perform highly accurate many-body simulations of realistic materials with all their intrinsic complications. The success of these simulations leaves us with a conundrum: how do we extract useful physical models and insight from these simulations? In this article, we present a formal theory of downfolding--extracting an…
▽ More
Due to advances in computer hardware and new algorithms, it is now possible to perform highly accurate many-body simulations of realistic materials with all their intrinsic complications. The success of these simulations leaves us with a conundrum: how do we extract useful physical models and insight from these simulations? In this article, we present a formal theory of downfolding--extracting an effective Hamiltonian from first-principles calculations. The theory maps the downfolding problem into fitting information derived from wave functions sampled from a low-energy subspace of the full Hilbert space. Since this fitting process most commonly uses reduced density matrices, we term it density matrix downfolding (DMD).
△ Less
Submitted 1 December, 2017;
originally announced December 2017.
-
Continuum of quantum fluctuations in a three-dimensional $S\!=\!1$ Heisenberg magnet
Authors:
K. W. Plumb,
Hitesh J. Changlani,
A. Scheie,
Shu Zhang,
J. W. Kriza,
J. A. Rodriguez-Rivera,
Yiming Qiu,
B. Winn,
R. J. Cava,
C. L. Broholm
Abstract:
Conventional crystalline magnets are characterized by symmetry breaking and normal modes of excitation called magnons with quantized angular momentum $\hbar$. Neutron scattering correspondingly features extra magnetic Bragg diffraction at low temperatures and dispersive inelastic scattering associated with single magnon creation and annihilation. Exceptions are anticipated in so-called quantum spi…
▽ More
Conventional crystalline magnets are characterized by symmetry breaking and normal modes of excitation called magnons with quantized angular momentum $\hbar$. Neutron scattering correspondingly features extra magnetic Bragg diffraction at low temperatures and dispersive inelastic scattering associated with single magnon creation and annihilation. Exceptions are anticipated in so-called quantum spin liquids as exemplified by the one-dimensional spin-1/2 chain which has no magnetic order and where magnons accordingly fractionalize into spinons with angular momentum $\hbar/2$. This is spectacularly revealed by a continuum of inelastic neutron scattering associated with two-spinon processes and the absence of magnetic Bragg diffraction. Here, we report evidence for these same key features of a quantum spin liquid in the three-dimensional Heisenberg antiferromagnet NaCaNi$_2$F$_7$. Through specific heat and neutron scattering measurements, Monte Carlo simulations, and analytic approximations to the equal time correlations, we show that NaCaNi$_2$F$_7$ is an almost ideal realization of the spin-1 antiferromagnetic Heisenberg model on a pyrochlore lattice with weak connectivity and frustrated interactions. Magnetic Bragg diffraction is absent and 90\% of the spectral weight forms a continuum of magnetic scattering not dissimilar to that of the spin-1/2 chain but with low energy pinch points indicating NaCaNi$_2$F$_7$ is in a Coulomb phase. The residual entropy and diffuse elastic scattering points to an exotic state of matter driven by frustration, quantum fluctuations and weak exchange disorder.
△ Less
Submitted 20 November, 2017;
originally announced November 2017.
-
Quantum versus classical effects at zero and finite temperature in the quantum pyrochlore Yb$_2$Ti$_2$O$_7$
Authors:
Hitesh J. Changlani
Abstract:
We study the finite temperature properties of the candidate quantum spin ice material Yb$_2$Ti$_2$O$_7$ within the framework of an anisotropic nearest-neighbor spin $1/2$ model on the pyrochlore lattice. Using a combination of finite temperature Lanczos and classical Monte Carlo methods, we highlight the importance of quantum mechanical effects for establishing the existence and location of the lo…
▽ More
We study the finite temperature properties of the candidate quantum spin ice material Yb$_2$Ti$_2$O$_7$ within the framework of an anisotropic nearest-neighbor spin $1/2$ model on the pyrochlore lattice. Using a combination of finite temperature Lanczos and classical Monte Carlo methods, we highlight the importance of quantum mechanical effects for establishing the existence and location of the low-temperature ordering transition. We perform simulations of the 32 site cluster, which capture the essential features of the specific heat curve seen in the cleanest known samples of this material. Focusing on recent experimental findings on Yb$_2$Ti$_2$O$_7$ [A. Scheie et al., Phys. Rev. Lett. 119, 127201 (2017) and J. D. Thompson et al., Phys. Rev. Lett. 119, 057203 (2017)], we then address the question of how the phase boundary between the ferromagnetic and paramagnetic phases changes when subjected to a magnetic field. We find that the quantum calculations explain discrepancies observed with a completely classical treatment and show that Yb$_2$Ti$_2$O$_7$ displays significant renormalization effects, which are at the heart of its reentrant lobed phase diagram. Finally, we develop a qualitative understanding of the existence of a ferromagnet by relating it to its counterpart that exists in the vicinity of the classical ice manifold.
△ Less
Submitted 18 May, 2018; v1 submitted 5 October, 2017;
originally announced October 2017.
-
Reentrant Phase Diagram of $\rm{Yb_2Ti_2O_7}$ in $\langle 111 \rangle$ Magnetic Field
Authors:
A. Scheie,
J. Kindervater,
S. Säubert,
C. Duvinage,
C. Pfleiderer,
H. J. Changlani,
S. Zhang,
L. Harriger,
K. Arpino,
S. M. Koohpayeh,
O. Tchernyshyov,
C. Broholm
Abstract:
We present a magnetic phase diagram of rare-earth pyrochlore $\rm{Yb_2Ti_2O_7}$ in a $\langle 111 \rangle$ magnetic field. Using heat capacity, magnetization, and neutron scattering data, we show an unusual field-dependence of a first-order phase boundary, wherein a small applied field increases the ordering temperature. The zero-field ground state has ferromagnetic domains, while the spins polari…
▽ More
We present a magnetic phase diagram of rare-earth pyrochlore $\rm{Yb_2Ti_2O_7}$ in a $\langle 111 \rangle$ magnetic field. Using heat capacity, magnetization, and neutron scattering data, we show an unusual field-dependence of a first-order phase boundary, wherein a small applied field increases the ordering temperature. The zero-field ground state has ferromagnetic domains, while the spins polarize along $\langle 111 \rangle$ above 0.65T. A classical Monte Carlo analysis of published Hamiltonians does account for the critical field in the low T limit. However, this analysis fails to account for the large bulge in the reentrant phase diagram, suggesting that either long-range interactions or quantum fluctuations govern low field properties.
△ Less
Submitted 22 June, 2017; v1 submitted 20 March, 2017;
originally announced March 2017.
-
The mother of all states of the kagome quantum antiferromagnet
Authors:
Hitesh J. Changlani,
Dmitrii Kochkov,
Krishna Kumar,
Bryan K. Clark,
Eduardo Fradkin
Abstract:
Frustrated quantum magnets are a central theme in condensed matter physics due to the richness of their phase diagrams. They support a panoply of phases including various ordered states and topological phases. Yet, this problem has defied a solution for a long time due to the lack of controlled approximations which make it difficult to distinguish between competing phases. Here we report the disco…
▽ More
Frustrated quantum magnets are a central theme in condensed matter physics due to the richness of their phase diagrams. They support a panoply of phases including various ordered states and topological phases. Yet, this problem has defied a solution for a long time due to the lack of controlled approximations which make it difficult to distinguish between competing phases. Here we report the discovery of a special quantum macroscopically degenerate point in the XXZ model on the spin-1/2 kagome quantum antiferromagnet for the ratio of Ising to antiferromagnetic transverse coupling Jz/J=-1/2. This point is proximate to many competing phases explaining the source of the complexity of the phase diagram. We identify five phases near this point including both spin-liquid and broken-symmetry phases and give evidence that the kagome Heisenberg antiferromagnet is close to a transition between two phases.
△ Less
Submitted 15 March, 2018; v1 submitted 14 March, 2017;
originally announced March 2017.
-
Quantum spin liquid with 7 elementary particles
Authors:
Haoyu Wang,
Hitesh J. Changlani,
Yuan Wan,
Oleg Tchernyshyov
Abstract:
We present an exactly solvable model of a quantum spin liquid with Abelian anyons in d=2 spatial dimensions. With spins 1/2 on a triangular lattice and six-body interactions, our model has zero spin correlation length and localized elementary excitations like the toric codes of Kitaev and Wen. In contrast to those earlier models, it has more elementary particles---4 bosons and 3 fermions---and hig…
▽ More
We present an exactly solvable model of a quantum spin liquid with Abelian anyons in d=2 spatial dimensions. With spins 1/2 on a triangular lattice and six-body interactions, our model has zero spin correlation length and localized elementary excitations like the toric codes of Kitaev and Wen. In contrast to those earlier models, it has more elementary particles---4 bosons and 3 fermions---and higher topological degeneracy of 64 on a torus. Elementary excitations are boson-fermion pairs that come in 12 distinct flavors. We use string operators to expose the topological nature of the model.
△ Less
Submitted 29 March, 2017; v1 submitted 4 February, 2017;
originally announced February 2017.
-
Numerical evidence for a chiral spin liquid in the XXZ antiferromagnetic Heisenberg model on the kagome lattice at $m=\frac{2}{3}$ magnetization
Authors:
Krishna Kumar,
Hitesh J. Changlani,
Bryan K. Clark,
Eduardo Fradkin
Abstract:
We perform an exact diagonalization study of the spin-$1/2$ XXZ Heisenberg antiferromagnet on the kagome lattice at finite magnetization $m = \frac{2}{3}$ with an emphasis on the XY point ($J_z = 0$), and in the presence of a small chiral term. Recent analytic work by Kumar, Sun and Fradkin [Phys. Rev. B 90, 174409 (2014)] on the same model, using a newly developed flux attachment transformation,…
▽ More
We perform an exact diagonalization study of the spin-$1/2$ XXZ Heisenberg antiferromagnet on the kagome lattice at finite magnetization $m = \frac{2}{3}$ with an emphasis on the XY point ($J_z = 0$), and in the presence of a small chiral term. Recent analytic work by Kumar, Sun and Fradkin [Phys. Rev. B 90, 174409 (2014)] on the same model, using a newly developed flux attachment transformation, predicts a plateau at this value of the magnetization described by a chiral spin liquid (CSL) with a spin Hall conductance of $σ_{xy} = \frac{1}{2}$. Such a state is topological in nature, has a ground state degeneracy and exhibits fractional excitations. We analyze the degeneracy structure in the low energy manifold, identify the candidate topological states and use them to compute the modular matrices and Chern numbers all of which strongly agree with expected theoretical behavior for the $σ_{xy} = \frac{1}{2}$ CSL. We argue that the evidence suggests the CSL is robust even in the limit of zero external chirality.
△ Less
Submitted 13 June, 2016;
originally announced June 2016.
-
Charge density waves in disordered media circumventing the Imry-Ma argument
Authors:
Hitesh J. Changlani,
Norm M. Tubman,
Taylor L. Hughes
Abstract:
Two powerful theoretical predictions, Anderson localization and the Imry-Ma argument, impose significant restrictions on the phases of matter that can exist in the presence of even the smallest amount of disorder in one-dimensional systems. These predictions forbid conducting states and ordered states respectively. It was thus remarkable that a mechanism to circumvent Anderson localization relying…
▽ More
Two powerful theoretical predictions, Anderson localization and the Imry-Ma argument, impose significant restrictions on the phases of matter that can exist in the presence of even the smallest amount of disorder in one-dimensional systems. These predictions forbid conducting states and ordered states respectively. It was thus remarkable that a mechanism to circumvent Anderson localization relying on the presence of correlated disorder was found, that is also realized in certain biomolecular systems. In a similar manner, we show that the Imry-Ma argument can be circumvented resulting in the formation of stable ordered states with discrete broken symmetries in disordered one dimensional systems. Specifically, we simulate a family of Hamiltonians of spinless fermions with correlated disorder and interactions, where we find that a charge density wave is stable up to a finite critical disorder strength. Having circumvented the Imry-Ma mechanism, we then investigate other mechanisms by which disorder can destroy an ordered state.
△ Less
Submitted 12 January, 2016;
originally announced January 2016.
-
Efficient heat-bath sampling in Fock space
Authors:
Adam Holmes,
Hitesh J. Changlani,
C. J. Umrigar
Abstract:
We introduce an algorithm for sampling many-body quantum states in Fock space. The algorithm efficiently samples states with probability approximately proportional to an arbitrary function of the second-quantized Hamiltonian matrix element connecting the sampled state to the current state. We apply the new sampling algorithm to the recently-developed Semistochastic Full Configuration Interaction Q…
▽ More
We introduce an algorithm for sampling many-body quantum states in Fock space. The algorithm efficiently samples states with probability approximately proportional to an arbitrary function of the second-quantized Hamiltonian matrix element connecting the sampled state to the current state. We apply the new sampling algorithm to the recently-developed Semistochastic Full Configuration Interaction Quantum Monte Carlo method (S-FCIQMC), a semistochastic implementation of the power method for projecting out the ground state energy in a basis of Slater determinants. The heat-bath sampling requires modest additional computational time and memory compared to uniform sampling but results in newly-spawned weights that are approximately of the same magnitude, thereby greatly improving the efficiency of projection. A comparison in efficiency between uniform and approximate heat-bath sampling is performed on the all-electron nitrogen dimer at equilibrium in Dunning's cc-pVXZ basis sets with X in {D, T, Q, 5}, demonstrating a large gain in efficiency that increases with basis set size. In addition, a comparison in efficiency is performed on three all-electron first-row dimers, B_2, N_2, and F_2, in a cc-pVQZ basis, demonstrating that the gain in efficiency compared to uniform sampling also increases dramatically with the number of electrons.
△ Less
Submitted 11 December, 2015;
originally announced December 2015.
-
Schwinger boson mean field perspective on emergent spins in diluted Heisenberg antiferromagnets
Authors:
Shivam Ghosh,
Hitesh J. Changlani,
Christopher L. Henley
Abstract:
Using an adaptation of Schwinger boson mean field theory (SBMFT) for nonuniform systems, we study the nature of low-energy spin excitations on the square and Bethe lattice at their percolation threshold. The optimal SBMFT parameters are interpreted as on-site potentials and pairing amplitudes, which enables an explanation of why emergent local moments develop in this system on dilution [L. Wang an…
▽ More
Using an adaptation of Schwinger boson mean field theory (SBMFT) for nonuniform systems, we study the nature of low-energy spin excitations on the square and Bethe lattice at their percolation threshold. The optimal SBMFT parameters are interpreted as on-site potentials and pairing amplitudes, which enables an explanation of why emergent local moments develop in this system on dilution [L. Wang and A. W. Sandvik, Phys. Rev. Lett. 97, 117204 (2006); H.J. Changlani et al., ibid 111, 157201 (2013)] and why the corresponding single particle frequencies are driven to anomalously low values. We discuss how our mean field calculations suggest the strong link between the presence of sublattice imbalance and long range antiferromagnetic order, and why linear spin wave theory is inadequate for capturing this relation. Within the SBMFT framework, we also extract an energy scale for the interaction between emergent moments, which show qualitative agreement with many-body calculations.
△ Less
Submitted 3 August, 2015; v1 submitted 21 May, 2015;
originally announced May 2015.
-
Density-matrix based determination of low-energy model Hamiltonians from ab initio wavefunctions
Authors:
Hitesh J. Changlani,
Huihuo Zheng,
Lucas K. Wagner
Abstract:
We propose a way of obtaining effective low energy Hubbard-like model Hamiltonians from ab initio Quantum Monte Carlo calculations for molecular and extended systems. The Hamiltonian parameters are fit to best match the ab initio two-body density matrices and energies of the ground and excited states, and thus we refer to the method as ab initio density matrix based downfolding (AIDMD). For benzen…
▽ More
We propose a way of obtaining effective low energy Hubbard-like model Hamiltonians from ab initio Quantum Monte Carlo calculations for molecular and extended systems. The Hamiltonian parameters are fit to best match the ab initio two-body density matrices and energies of the ground and excited states, and thus we refer to the method as ab initio density matrix based downfolding (AIDMD). For benzene (a finite system), we find good agreement with experimentally available energy gaps without using any experimental inputs. For graphene, a two dimensional solid (extended system) with periodic boundary conditions, we find the effective on-site Hubbard $U^{*}/t$ to be $1.3 \pm 0.2$, comparable to a recent estimate based on the constrained random phase approximation. For molecules, such parameterizations enable calculation of excited states that are usually not accessible within ground state approaches. For solids, the effective Hamiltonian enables large scale calculations using techniques designed for lattice models.
△ Less
Submitted 4 August, 2015; v1 submitted 14 April, 2015;
originally announced April 2015.
-
Phase diagram of the $Z_3$ Parafermionic Chain with Chiral Interactions
Authors:
Ye Zhuang,
Hitesh J. Changlani,
Norm M. Tubman,
Taylor L. Hughes
Abstract:
Parafermions are exotic quasiparticles with non-Abelian fractional statistics that can be realized and stabilized in 1-dimensional models that are generalizations of the Kitaev p-wave wire. We study the simplest generalization, i.e. the $Z_3$ parafermionic chain. Using a Jordan-Wigner transform we focus on the equivalent three-state chiral clock model, and study its rich phase diagram using the de…
▽ More
Parafermions are exotic quasiparticles with non-Abelian fractional statistics that can be realized and stabilized in 1-dimensional models that are generalizations of the Kitaev p-wave wire. We study the simplest generalization, i.e. the $Z_3$ parafermionic chain. Using a Jordan-Wigner transform we focus on the equivalent three-state chiral clock model, and study its rich phase diagram using the density matrix renormalization group technique. We perform our analyses using quantum entanglement diagnostics which allow us to determine phase boundaries, and the nature of the phase transitions. In particular, we study the transition between the topological and trivial phases, as well as to an intervening incommensurate phase which appears in a wide region of the phase diagram. The phase diagram is predicted to contain a Lifshitz type transition which we confirm using entanglement measures. We also attempt to locate and characterize a putative tricritical point in the phase diagram where the three above mentioned phases meet at a single point.
△ Less
Submitted 31 July, 2015; v1 submitted 17 February, 2015;
originally announced February 2015.
-
Determination of Tomonaga-Luttinger parameters for a two-component liquid
Authors:
Olabode M. Sule,
Hitesh J. Changlani,
Isao Maruyama,
Shinsei Ryu
Abstract:
We provide evidence for the mapping of critical spin-1 chains, in particular the SU(3) symmetric bilinear-biquadratic model with additional interactions, to free boson theories using exact diagonalization and the density matrix renormalization group algorithm. Using the correspondence with a conformal field theory with central charge c=2, we determine the analytic formulae for the scaling dimensio…
▽ More
We provide evidence for the mapping of critical spin-1 chains, in particular the SU(3) symmetric bilinear-biquadratic model with additional interactions, to free boson theories using exact diagonalization and the density matrix renormalization group algorithm. Using the correspondence with a conformal field theory with central charge c=2, we determine the analytic formulae for the scaling dimensions in terms of four Tomonaga-Luttinger liquid parameters. By matching the lowest scaling dimensions, we numerically calculate these field-theoretic parameters and track their evolution as a function of the parameters of the lattice model.
△ Less
Submitted 28 January, 2016; v1 submitted 8 February, 2015;
originally announced February 2015.
-
Density-matrix based numerical methods for discovering order and correlations in interacting systems
Authors:
Christopher L. Henley,
Hitesh J. Changlani
Abstract:
We review recently introduced numerical methods for the unbiased detection of the order parameter and/or dominant correlations, in many-body interacting systems, by using reduced density matrices. Most of the paper is devoted to the "quasi-degenerate density matrix" (QDDM) which is rooted in Anderson's observation that the degenerate symmetry-broken states valid in the thermodynamic limit, are man…
▽ More
We review recently introduced numerical methods for the unbiased detection of the order parameter and/or dominant correlations, in many-body interacting systems, by using reduced density matrices. Most of the paper is devoted to the "quasi-degenerate density matrix" (QDDM) which is rooted in Anderson's observation that the degenerate symmetry-broken states valid in the thermodynamic limit, are manifested in finite systems as a set of low-energy "quasi-degenerate" states (in addition to the ground state). This method, its original form due to Furukawa et al.[Phys. Rev. Lett. 96, 047211 (2006)], is given a number of improvements here, above all the extension from two-fold symmetry breaking to arbitrary cases. This is applied to two test cases (1) interacting spinless hardcore bosons on the triangular lattice and (2) a spin-1/2 antiferromagnetic system at the percolation threshold. In addition, we survey a different method called the "correlation density matrix", which detects (possibly long-range) correlations only from the ground state, but using the reduced density matrix from a cluster consisting of two spatially separated regions.
△ Less
Submitted 15 July, 2014;
originally announced July 2014.
-
Trimerized ground state of the spin-1 Heisenberg antiferromagnet on the kagome lattice
Authors:
Hitesh J. Changlani,
Andreas M. Läuchli
Abstract:
We study the phase diagram of the spin-1 quantum bilinear-biquadratic antiferromagnet on the kagome lattice using exact diagonalization and the density matrix renormalization group algorithm. The SU(3)-symmetric point of this model Hamiltonian is a spontaneously trimerized state whose qualitative nature persists even at the Heisenberg point, a finding that contrasts previous proposals. We report t…
▽ More
We study the phase diagram of the spin-1 quantum bilinear-biquadratic antiferromagnet on the kagome lattice using exact diagonalization and the density matrix renormalization group algorithm. The SU(3)-symmetric point of this model Hamiltonian is a spontaneously trimerized state whose qualitative nature persists even at the Heisenberg point, a finding that contrasts previous proposals. We report the ground state energy per site of the Heisenberg model to be -1.410(2) and establish the presence of a spin gap.
△ Less
Submitted 28 January, 2016; v1 submitted 18 June, 2014;
originally announced June 2014.
-
Stochastically Projecting Tensor Networks
Authors:
Bryan K. Clark,
Hitesh J. Changlani
Abstract:
We apply a series of projection techniques on top of tensor networks to compute energies of ground state wave functions with higher accuracy than tensor networks alone with minimal additional cost. We consider both matrix product states as well as tree tensor networks in this work. Building on top of these approaches, we apply fixed-node quantum Monte Carlo, Lanczos steps, and exact projection. We…
▽ More
We apply a series of projection techniques on top of tensor networks to compute energies of ground state wave functions with higher accuracy than tensor networks alone with minimal additional cost. We consider both matrix product states as well as tree tensor networks in this work. Building on top of these approaches, we apply fixed-node quantum Monte Carlo, Lanczos steps, and exact projection. We demonstrate these improvements for the triangular lattice Heisenberg model, where we capture up to 57 percent of the remaining energy not captured by the tensor network alone. We conclude by discussing further ways to improve our approach.
△ Less
Submitted 8 April, 2014;
originally announced April 2014.
-
Emergent Spin Excitations in a Bethe Lattice at Percolation
Authors:
Hitesh J. Changlani,
Shivam Ghosh,
Sumiran Pujari,
Christopher L. Henley
Abstract:
We study the spin 1/2 quantum Heisenberg antiferromagnet on a Bethe lattice diluted to the percolation threshold. Dilution creates areas of even/odd sublattice imbalance resulting in "dangling spins" (L. Wang and A. W. Sandvik, Phys. Rev. Lett. 97, 117204 (2006); Phys. Rev. B 81, 054417 (2010)). These collectively act as "emergent" spin 1/2 degrees of freedom and are responsible for the creation o…
▽ More
We study the spin 1/2 quantum Heisenberg antiferromagnet on a Bethe lattice diluted to the percolation threshold. Dilution creates areas of even/odd sublattice imbalance resulting in "dangling spins" (L. Wang and A. W. Sandvik, Phys. Rev. Lett. 97, 117204 (2006); Phys. Rev. B 81, 054417 (2010)). These collectively act as "emergent" spin 1/2 degrees of freedom and are responsible for the creation of a set of low lying "quasidegenerate states". Using Density Matrix Renormalization Group (DMRG) calculations, we detect the presence and location of these emergent spins. We find an effective Hamiltonian of these emergent spins, with Heisenberg interactions that decay exponentially with the distance between them.
△ Less
Submitted 8 October, 2013; v1 submitted 16 October, 2012;
originally announced October 2012.
-
Heisenberg antiferromagnet on Cayley trees: low-energy spectrum and even/odd site imbalance
Authors:
Hitesh J. Changlani,
Shivam Ghosh,
Christopher L. Henley,
Andreas M. Läuchli
Abstract:
To understand the role of local sublattice imbalance in low-energy spectra of s=1/2 quantum antiferromagnets, we study the s=1/2 quantum nearest neighbor Heisenberg antiferromagnet on the coordination 3 Cayley tree. We perform many-body calculations using an implementation of the density matrix renormalization group (DMRG) technique for generic tree graphs. We discover that the bond-centered Cayle…
▽ More
To understand the role of local sublattice imbalance in low-energy spectra of s=1/2 quantum antiferromagnets, we study the s=1/2 quantum nearest neighbor Heisenberg antiferromagnet on the coordination 3 Cayley tree. We perform many-body calculations using an implementation of the density matrix renormalization group (DMRG) technique for generic tree graphs. We discover that the bond-centered Cayley tree has a quasidegenerate set of a low-lying tower of states and an "anomalous" singlet-triplet finite-size gap scaling. For understanding the construction of the first excited state from the many-body ground state, we consider a wave function ansatz given by the single-mode approximation, which yields a high overlap with the DMRG wave function. Observing the ground-state entanglement spectrum leads us to a picture of the low-energy degrees of freedom being "giant spins" arising out of sublattice imbalance, which helps us analytically understand the scaling of the finite-size spin gap. The Schwinger-boson mean-field theory has been generalized to nonuniform lattices, and ground states have been found which are spatially inhomogeneous in the mean-field parameters.
△ Less
Submitted 11 March, 2013; v1 submitted 8 August, 2012;
originally announced August 2012.
-
Semistochastic Projector Monte Carlo Method
Authors:
F. R. Petruzielo,
A. A. Holmes,
Hitesh J. Changlani,
M. P. Nightingale,
C. J. Umrigar
Abstract:
We introduce a semistochastic implementation of the power method to compute, for very large matrices, the dominant eigenvalue and expectation values involving the corresponding eigenvector. The method is semistochastic in that the matrix multiplication is partially implemented numerically exactly and partially with respect to expectation values only. Compared to a fully stochastic method, the semi…
▽ More
We introduce a semistochastic implementation of the power method to compute, for very large matrices, the dominant eigenvalue and expectation values involving the corresponding eigenvector. The method is semistochastic in that the matrix multiplication is partially implemented numerically exactly and partially with respect to expectation values only. Compared to a fully stochastic method, the semistochastic approach significantly reduces the computational time required to obtain the eigenvalue to a specified statistical uncertainty. This is demonstrated by the application of the semistochastic quantum Monte Carlo method to systems with a sign problem: the fermion Hubbard model and the carbon dimer.
△ Less
Submitted 22 October, 2012; v1 submitted 25 July, 2012;
originally announced July 2012.
-
Approximating strongly correlated spin and fermion wavefunctions with correlator product states
Authors:
Hitesh J. Changlani,
Jesse M. Kinder,
Cyrus J. Umrigar,
Garnet Kin-Lic Chan
Abstract:
We explore correlator product states for the approximation of correlated wavefunctions in arbitrary dimensions. We show that they encompass many interesting states including Laughlin's quantum Hall wavefunction, Huse and Elser's frustrated spin states, and Kitaev's toric code. We further establish their relation to common families of variational wavefunctions, such as matrix and tensor product s…
▽ More
We explore correlator product states for the approximation of correlated wavefunctions in arbitrary dimensions. We show that they encompass many interesting states including Laughlin's quantum Hall wavefunction, Huse and Elser's frustrated spin states, and Kitaev's toric code. We further establish their relation to common families of variational wavefunctions, such as matrix and tensor product states and resonating valence bond states. Calculations on the Heisenberg and spinless Hubbard models show that correlator product states capture both two-dimensional correlations (independent of system width) as well as non-trivial fermionic correlations (without sign problems). In one-dimensional simulations, correlator product states appear competitive with matrix product states with a comparable number of variational parameters, suggesting they may eventually provide a route to practically generalise the density matrix renormalisation group to higher dimensions.
△ Less
Submitted 29 October, 2009; v1 submitted 27 July, 2009;
originally announced July 2009.