-
Quantum Criticality in the infinite-range Transverse Field Ising Model
Authors:
Nicholas Curro,
Kaeshav Danesh,
Rajiv R. P. Singh
Abstract:
We study quantum criticality in the infinite range Transverse-Field Ising Model. We find subtle differences with respect to the well-known single-site mean-field theory, especially in terms of gap, entanglement and quantum criticality. The calculations are based on numerical diagonalization of Hamiltonians with up to a few thousand spins. This is made possible by the enhanced symmetries of the mod…
▽ More
We study quantum criticality in the infinite range Transverse-Field Ising Model. We find subtle differences with respect to the well-known single-site mean-field theory, especially in terms of gap, entanglement and quantum criticality. The calculations are based on numerical diagonalization of Hamiltonians with up to a few thousand spins. This is made possible by the enhanced symmetries of the model, which divide the Hamiltonian into many block-diagonal sectors. The finite temperature phase diagram and the characteristic jump in heat capacity closely resemble the behavior in mean-field theory. However, unlike mean-field theory where excitations are always gapped, the excitation gap in the infinite range model goes to zero from both the paramagnetic side and from the ferromagnetic side on approach to the quantum critical point. Also, contrary to mean-field theory, at the quantum critical point the Quantum Fisher Information becomes large, implying long-range multi-partite entanglement. We find that the main role of temperature is to shift statistical weights from one conserved sector to another. However, low energy excitations in each sector arise only near the quantum critical point implying that low energy quantum fluctuations can arise only in the vicinity of the quantum critical field where they can persist up to temperatures of order the exchange constant.
△ Less
Submitted 5 August, 2024;
originally announced August 2024.
-
Nuclear magnetic resonance studies in a model transverse field Ising system
Authors:
Y-H. Nian,
I. Vinograd,
C. Chaffey,
Y. Li,
M. P. Zic,
P. Massat,
R. R. P. Singh,
I. R. Fisher,
N. J. Curro
Abstract:
The suppression of ferroquadrupolar order in TmVO$_4$ in a magnetic field is well-described by the transverse field Ising model, enabling detailed studies of critical dynamics near the quantum phase transition. We describe nuclear magnetic resonance measurements in pure and Y-doped single crystals. The non-Kramers nature of the ground state doublet leads to a unique form of the hyperfine coupling…
▽ More
The suppression of ferroquadrupolar order in TmVO$_4$ in a magnetic field is well-described by the transverse field Ising model, enabling detailed studies of critical dynamics near the quantum phase transition. We describe nuclear magnetic resonance measurements in pure and Y-doped single crystals. The non-Kramers nature of the ground state doublet leads to a unique form of the hyperfine coupling that exclusively probes the transverse field susceptibility. Our results show that this quantity diverges at the critical field, in contrast to the mean-field prediction. Furthermore, we find evidence for quantum critical fluctuations present near Tm-rich regions in Y-doped crystals at levels beyond which long-range order is suppressed, suggesting the presence of quantum Griffiths phases.
△ Less
Submitted 28 February, 2024;
originally announced February 2024.
-
Competing Gauge Fields and Entropically-Driven Spin Liquid to Spin Liquid Transition in non-Kramers Pyrochlores
Authors:
Daniel Lozano-Gómez,
Vincent Noculak,
Jaan Oitmaa,
Rajiv R. P. Singh,
Yasir Iqbal,
Johannes Reuther,
Michel J. P. Gingras
Abstract:
Gauge theories are powerful tools in theoretical physics, allowing complex phenomena to be reduced to simple principles, and are used in both high-energy and condensed matter physics. In the latter context, gauge theories are becoming increasingly popular for capturing the intricate spin correlations in spin liquids, exotic states of matter in which the dynamics of quantum spins never ceases, even…
▽ More
Gauge theories are powerful tools in theoretical physics, allowing complex phenomena to be reduced to simple principles, and are used in both high-energy and condensed matter physics. In the latter context, gauge theories are becoming increasingly popular for capturing the intricate spin correlations in spin liquids, exotic states of matter in which the dynamics of quantum spins never ceases, even at absolute zero temperature. We consider a spin system on a three-dimensional pyrochlore lattice where emergent gauge fields not only describe the spin liquid behaviour at zero temperature but crucially determine the system's temperature evolution, with distinct gauge fields giving rise to different spin liquid phases in separate temperature regimes. Focusing first on classical spins, in an intermediate temperature regime, the system shows an unusual coexistence of emergent vector and matrix gauge fields where the former is known from classical spin ice systems while the latter has been associated with fractonic quasiparticles, a peculiar type of excitation with restricted mobility. Upon cooling, the system transitions into a low-temperature phase where an entropic selection mechanism depopulates the degrees of freedom associated with the matrix gauge field, rendering the system spin ice like. We further provide numerical evidence that in the corresponding quantum model, a spin liquid with coexisting vector and matrix gauge fields has a finite window of stability in the parameter space of spin interactions down to zero temperature. Finally, we discuss the relevance of our findings for non-Kramers pyrochlore materials.
△ Less
Submitted 3 September, 2024; v1 submitted 25 October, 2023;
originally announced October 2023.
-
Spin-echo and quantum versus classical critical fluctuations in TmVO$_4$
Authors:
Y-H. Nian,
I. Vinograd,
T. Green,
C. Chaffey,
P. Massat,
R. R. P. Singh,
M. P. Zic,
I. R. Fisher,
N. J. Curro
Abstract:
Using spin-echo Nuclear Magnetic Resonance in the model Transverse-Field Ising system TmVO$_4$, we show that low frequency quantum fluctuations at the quantum critical point have a very different effect on $^{51}$V nuclear-spins than classical low-frequency noise or fluctuations that arise at a finite temperature critical point. Spin-echos filter out the low frequency classical noise but not the q…
▽ More
Using spin-echo Nuclear Magnetic Resonance in the model Transverse-Field Ising system TmVO$_4$, we show that low frequency quantum fluctuations at the quantum critical point have a very different effect on $^{51}$V nuclear-spins than classical low-frequency noise or fluctuations that arise at a finite temperature critical point. Spin-echos filter out the low frequency classical noise but not the quantum fluctuations. This allows us to directly visualize the quantum critical fan and demonstrate the persistence of quantum fluctuations at the critical coupling strength in TmVO$_4$ to high temperatures in an experiment that remains transparent to finite temperature classical phase transitions. These results show that while dynamical decoupling schemes can be quite effective in eliminating classical noise in a qubit, a quantum critical environment may lead to rapid entanglement and decoherence.
△ Less
Submitted 22 June, 2023;
originally announced June 2023.
-
High-temperature magnetization and entropy of the triangular lattice Hubbard model in a Zeeman field
Authors:
Owen Bradley,
Yutan Zhang,
Jaan Oitmaa,
Rajiv R. P. Singh
Abstract:
We use strong coupling expansions to calculate the entropy function $S(T,h)$, the magnetization $M(T,h)$, and the double occupancy factor $D(T,h)$ for the half-filled triangular lattice Hubbard model as a function of temperature $T$ and Zeeman field $h$, for various values of the Hubbard parameter ratio $U/t$. These calculations converge well for temperatures larger than the exchange parameter…
▽ More
We use strong coupling expansions to calculate the entropy function $S(T,h)$, the magnetization $M(T,h)$, and the double occupancy factor $D(T,h)$ for the half-filled triangular lattice Hubbard model as a function of temperature $T$ and Zeeman field $h$, for various values of the Hubbard parameter ratio $U/t$. These calculations converge well for temperatures larger than the exchange parameter $J=\frac{4t^2}{U}$ for moderate to large $U/t$ values. Setting $μ=U/2$ suffices to obtain the density of half filling within a fraction of one percent at all temperatures studied for $U/t \geq 8$. We discuss the systematic variation of properties with $U/t$. The temperature dependence of entropy and the double occupancy parameter shows a mapping to an antiferromagnetic Mott insulating behavior at temperatures well above $T=J$ for $U/t\ge 10$. Convergence of the series is weaker at intermediate fields implying non-monotonic variation of spin-correlations with the Zeeman field. We discuss the relevance of the Hubbard model results to the triangular-lattice antiferromagnetic materials Lu$_3$Cu$_2$Sb$_3$O$_{14} $ (LCSO) studied recently by Yang et al [Yang et al arXiv:2102.09271 (2022)].
△ Less
Submitted 6 March, 2023;
originally announced March 2023.
-
Classical and quantum phases of the pyrochlore $S=1/2$ magnet with Heisenberg and Dzyaloshinskii-Moriya interactions
Authors:
Vincent Noculak,
Daniel Lozano-Gómez,
Jaan Oitmaa,
Rajiv R. P. Singh,
Yasir Iqbal,
Michel J. P. Gingras,
Johannes Reuther
Abstract:
We investigate the ground state and critical temperature phase diagrams of the classical and quantum $S=1/2$ pyrochlore lattice with nearest-neighbor Heisenberg and Dzyaloshinskii-Moriya interactions (DMI). We consider ferromagnetic and antiferromagnetic Heisenberg exchange as well as direct and indirect DMI. Classically, three ground states are found: all-in/all-out, ferromagnetic and a locally o…
▽ More
We investigate the ground state and critical temperature phase diagrams of the classical and quantum $S=1/2$ pyrochlore lattice with nearest-neighbor Heisenberg and Dzyaloshinskii-Moriya interactions (DMI). We consider ferromagnetic and antiferromagnetic Heisenberg exchange as well as direct and indirect DMI. Classically, 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. Quantum zero-point energy fluctuations are found to lift the classical ground state degeneracy and select the $ψ_3$ state in most parts of the $Γ_5$ regime. Likewise, thermal fluctuations treated classically, select the $ψ_3$ state at $T=0^+$. In contrast, classical Monte Carlo finds that the system orders at $T_c$ in the $ψ_2$ state of $Γ_5$ for antiferromagnetic Heisenberg exchange and indirect DMI with a transition from $ψ_2$ to $ψ_3$ at a temperature $T_{Γ_5} <T_c$. The same method finds that the system orders via a single transition at $T_c$ directly into the $ψ_3$ state for most of the region with ferromagnetic Heisenberg exchange and indirect DMI. Such ordering behavior at $T_c$ for the $S=1/2$ quantum model is corroborated by high-temperature series expansion. To investigate the $T=0$ quantum ground states, we apply the pseudo-fermion functional renormalization group (PFFRG). The quantum paramagnetic phase of the pure antiferromagnetic $S=1/2$ Heisenberg model is found to persist over a finite region in the phase diagram for both direct or indirect DMI. We find that near the boundary of ferromagnetism and $Γ_5$ antiferromagnetism the system may potentially realize a quantum ground state lacking conventional magnetic order. Otherwise, for the largest portion of the phase diagram, PFFRG finds the same ordered phases as in the classical model.
△ Less
Submitted 20 June, 2023; v1 submitted 16 November, 2022;
originally announced November 2022.
-
Onset of charge incompressibility and Mott gaps in the Honeycomb-Lattice SU(4) Hubbard Model: Lessons for Twisted Bilayer Graphene systems
Authors:
Rahul Hingorani,
Jaan Oitmaa,
Rajiv R. P. Singh
Abstract:
We use finite temperature strong coupling expansions to calculate thermodynamic properties of the Honeycomb-lattice SU(4) Hubbard model. We present numerical results for various properties including chemical potential, compressibility, entropy and specific heat as a function of temperature and density at several $U/t$ values. We study the onset of charge incompressibility and Mott gaps as the temp…
▽ More
We use finite temperature strong coupling expansions to calculate thermodynamic properties of the Honeycomb-lattice SU(4) Hubbard model. We present numerical results for various properties including chemical potential, compressibility, entropy and specific heat as a function of temperature and density at several $U/t$ values. We study the onset of charge incompressibility and Mott gaps as the temperature is lowered at integer densities. In the incompressible Mott regime, the expansions are recast into a high temperature expansion for a generalized spin model with SU(4) symmetry, which is then used to study the convergence of strong coupling expansions in t/U. We discuss lessons that can be drawn from high temperature properties of a simple Hubbard model regarding Twisted Bilayer Graphene (TBG) and other magic-angle flat-band systems.
△ Less
Submitted 27 March, 2022;
originally announced March 2022.
-
Order, Disorder and Monopole Confinement in the Spin-$1/2$ XXZ Model on a Pyrochlore Tube
Authors:
Chunhan Feng,
Alexander Wietek,
E. Miles Stoudenmire,
Rajiv R. P. Singh
Abstract:
We study the ground state and thermodynamic properties of the spin-half XXZ model, with an Ising interaction $J_z$ and a transverse exchange interaction $J_{x}$, on a pyrochlore tube obtained by joining together elementary cubes in a one-dimensional array. Periodic boundary conditions in the transverse directions ensure that the bulk of the system consists of corner-sharing tetrahedra, with the sa…
▽ More
We study the ground state and thermodynamic properties of the spin-half XXZ model, with an Ising interaction $J_z$ and a transverse exchange interaction $J_{x}$, on a pyrochlore tube obtained by joining together elementary cubes in a one-dimensional array. Periodic boundary conditions in the transverse directions ensure that the bulk of the system consists of corner-sharing tetrahedra, with the same local geometry as the pyrochlore lattice. We use exact diagonalization, the density matrix renormalization group (DMRG), and minimally entangled typical thermal states (METTS) methods to study the system. When $J_z$ is antiferromagnetic ($J_{z}>0$) and $J_x$ is ferromagnetic ($J_{x}<0$), we find a transition from a spin liquid to an XY ferromagnet, which has power-law correlations at $T=0$. For $J_{z}<0$ and $J_{x}>0$, spin-two excitations are found to have lower energy than spin-one at the transition away from the fully polarized state, showing evidence for incipient spin-nematic order. When both interactions are antiferromagnetic, we find a non-degenerate ground state with no broken symmetries and a robust energy gap. The low energy spectra evolve smoothly from predominantly Ising to predominantly XY interactions. In the spin-liquid regime of small $|J_{x}|$, we study the confinement of monopole-anti-monopole pairs and find that the confinement length scale is larger for $J_x<0$ than for $J_{x}>0$, although both length scales are very short. These results are consistent with a local spin-liquid phase for the Heisenberg antiferromagnet with no broken symmetries.
△ Less
Submitted 28 February, 2022;
originally announced March 2022.
-
On the Divergence of the Ferromagnetic Susceptibility in the SU(N) Nagaoka-Thouless Ferromagnet
Authors:
Rajiv R. P. Singh,
Jaan Oitmaa
Abstract:
Using finite temperature strong coupling expansions for the SU(N) Hubbard Model, we calculate the thermodynamic properties of the model in the infinite-$U$ limit for arbitrary density $0\leq ρ\leq 1$ and all $N$. We express the ferromagnetic susceptibility of the model as a Curie term plus a $Δχ$, an excess susceptibility above the Curie-behavior. We show that, on a bipartite lattice, graph by gra…
▽ More
Using finite temperature strong coupling expansions for the SU(N) Hubbard Model, we calculate the thermodynamic properties of the model in the infinite-$U$ limit for arbitrary density $0\leq ρ\leq 1$ and all $N$. We express the ferromagnetic susceptibility of the model as a Curie term plus a $Δχ$, an excess susceptibility above the Curie-behavior. We show that, on a bipartite lattice, graph by graph the contributions to $Δχ$ are non-negative in the limit that the hole density $δ=1-ρ$ goes to zero. By summing the contributions from all graphs consisting of closed loops we find that the low hole-density ferromagnetic susceptibility diverges exponentially as $\exp{Δ/T}$ as $T \to 0$ in two and higher dimensions. This demonstrates that Nagaoka-Thouless ferromagnetic state exists as a thermodynamic state of matter at low enough density of holes and sufficiently low temperatures. The constant $Δ$ scales with the SU(N) parameter $N$ as $1/N$ implying that ferromagnetism is gradually weakened with increasing $N$ as the characteristic temperature scale for ferromagnetic order goes down.
△ Less
Submitted 3 February, 2022;
originally announced February 2022.
-
Finite Temperature Strong Coupling Expansions for the SU(N) Hubbard Model
Authors:
Rajiv R. P. Singh,
Jaan Oitmaa
Abstract:
We develop finite temperature strong coupling expansions for the SU(N) Hubbard Model in powers of $βt$, $w=\exp{(-βU)}$ and ${1\over βU}$ for arbitrary filling. The expansions are done in the grand canonical ensemble and are most useful at a density of one particle per site, where for $U$ larger than or of order the Bandwidth, the expansions converge over a wide temperature range…
▽ More
We develop finite temperature strong coupling expansions for the SU(N) Hubbard Model in powers of $βt$, $w=\exp{(-βU)}$ and ${1\over βU}$ for arbitrary filling. The expansions are done in the grand canonical ensemble and are most useful at a density of one particle per site, where for $U$ larger than or of order the Bandwidth, the expansions converge over a wide temperature range $t^2/U \ \lesssim \ T \ \lesssim \ 10 U$. By taking the limit $w\to 0$, valid at temperatures much less than $U$, the expansions turn into a high temperature expansion for a dressed SU(N) Heisenberg model that includes nearest-neighbor exchange, further neighbor exchanges and ring exchanges known from the $T=0$ perturbation theory of the SU(2) Hubbard model. Below a filling of one particle per site, the $w\to 0$ limit corresponds to an effective $t-J$ model. The onset of strong correlations can be identified by a plateau-like behavior in the entropy as a function of temperature. At small deviations from one particle per site, the expansions can be arranged in powers of a small parameter $δ=1-n$, the deviation from one particle per site, where the leading $βt$ dependent terms correspond to holes sloshing around in a disordered SU(N) background. We use these expansions to calculate the thermodynamic properties of the model at moderate and high temperatures over a wide parameter range.
△ Less
Submitted 20 January, 2022; v1 submitted 17 January, 2022;
originally announced January 2022.
-
Magnetism of competing high-spin/low-spin states in Ba$_2$NiO$_2$(AgSe)$_2$ and related two-orbital two-electron systems
Authors:
Rajiv R. P. Singh
Abstract:
We discuss the thermodynamic and magnetic properties of a competing high-spin/low-spin two-orbital two-electron model on a square-lattice possibly relevant to the nickelates such as Ba$_2$NiO$_2$(AgSe)$_2$ (BNOAS). We focus on parameter regimes where a high-spin ($S=1$) and a low-spin ($S=0$) state are energetically close to each other and discuss various exchange processes in such a system. The m…
▽ More
We discuss the thermodynamic and magnetic properties of a competing high-spin/low-spin two-orbital two-electron model on a square-lattice possibly relevant to the nickelates such as Ba$_2$NiO$_2$(AgSe)$_2$ (BNOAS). We focus on parameter regimes where a high-spin ($S=1$) and a low-spin ($S=0$) state are energetically close to each other and discuss various exchange processes in such a system. The model we study is a variant of, but different from the Kondo-necklace model proposed for the system by Jin et al (H. S. Jin, W. E. Pickett and K. W. Lee, Phys. Rev. Research 2, 033197 (2020)). Although there are similarities between the two models in terms of the ground-state phases and their symmetries, the detailed properties of the different phases and the phase transitions are entirely different and should be easy to distinguish from an experimental point of view.
△ Less
Submitted 30 August, 2021; v1 submitted 22 August, 2021;
originally announced August 2021.
-
Instabilities of spin-1 Kitaev spin liquid phase in presence of single-ion anisotropies
Authors:
Owen Bradley,
Rajiv R. P. Singh
Abstract:
We study the spin-one Kitaev model on the honeycomb lattice in the presence of single-ion anisotropies. We consider two types of single ion anisotropies: A $D_{111}$ anisotropy which preserves the symmetry between $X$, $Y$, and $Z$ bonds but violates flux conservation and a $D_{100}$ anisotropy that breaks the symmetry between $X$, $Y$, and $Z$ bonds but preserves flux conservation. We use series…
▽ More
We study the spin-one Kitaev model on the honeycomb lattice in the presence of single-ion anisotropies. We consider two types of single ion anisotropies: A $D_{111}$ anisotropy which preserves the symmetry between $X$, $Y$, and $Z$ bonds but violates flux conservation and a $D_{100}$ anisotropy that breaks the symmetry between $X$, $Y$, and $Z$ bonds but preserves flux conservation. We use series expansion methods, degenerate perturbation theory, and exact diagonalization to study these systems. Large positive $D_{111}$ anisotropy leads to a simple product ground state with conventional magnon-like excitations, while large negative $D_{111}$ leads to a broken symmetry and degenerate ground states. For both signs there is a phase transition at a small $|D_{111}| \approx 0.12$ separating the more conventional phases from the Kitaev spin liquid phase. With large $D_{100}$ anisotropy, the ground state is a simple product state, but the model lacks conventional dispersive excitations due to the large number of conservation laws. Large negative $D_{100}$ leads to decoupled one-dimensional systems and many degenerate ground states. No evidence of a phase transition is seen in our numerical studies at any finite $D_{100}$. Convergence of the series expansion extrapolations all the way to $D_{100}=0$ suggests that the nontrivial Kitaev spin liquid is a singular limit of this type of single-ion anisotropy going to zero, which also restores symmetry between the $X$, $Y$, and $Z$ bonds.
△ Less
Submitted 20 April, 2022; v1 submitted 11 August, 2021;
originally announced August 2021.
-
Thermodynamic behavior of modified integer-spin Kitaev models on the honeycomb lattice
Authors:
Owen Bradley,
Jaan Oitmaa,
Diptiman Sen,
Rajiv R. P. Singh
Abstract:
We study the thermodynamic properties of modified spin-$S$ Kitaev models introduced by Baskaran, Sen and Shankar (Phys. Rev. B 78, 115116 (2008)). These models have the property that for half-odd-integer spins their eigenstates map on to those of spin-1/2 Kitaev models, with well-known highly entangled quantum spin-liquid states and Majorana fermions. For integer spins, the Hamiltonian is made out…
▽ More
We study the thermodynamic properties of modified spin-$S$ Kitaev models introduced by Baskaran, Sen and Shankar (Phys. Rev. B 78, 115116 (2008)). These models have the property that for half-odd-integer spins their eigenstates map on to those of spin-1/2 Kitaev models, with well-known highly entangled quantum spin-liquid states and Majorana fermions. For integer spins, the Hamiltonian is made out of commuting local operators. Thus, the eigenstates can be chosen to be completely unentangled between different sites, though with a significant degeneracy for each eigenstate. For half-odd-integer spins, the thermodynamic properties can be related to the spin-1/2 Kitaev models apart from an additional degeneracy. Hence we focus here on the case of integer spins. We use transfer matrix methods, high temperature expansions and Monte Carlo simulations to study the thermodynamic properties of ferromagnetic and antiferromagnetic models with spin $S=1$ and $S=2$. Apart from large residual entropies, which all the models have, we find that they can have a variety of different behaviors. Transfer matrix calculations show that for the different models, the correlation lengths can be finite as $T\to 0$, become critical as $T\to 0$ or diverge exponentially as $T\to 0$. There is a conserved $Z_2$ flux variable associated with each hexagonal plaquette which saturates at the value $+1$ as $T\rightarrow0$ in all models except the $S=1$ antiferromagnet where the mean flux remains zero as $T\to 0$. We provide qualitative explanations for these results.
△ Less
Submitted 8 February, 2021; v1 submitted 3 May, 2020;
originally announced May 2020.
-
On two phases inside the Bose condensation dome of Yb$_2$Si$_2$O$_7$
Authors:
Michael O. Flynn,
Thomas E. Baker,
Siddharth Jindal,
Rajiv R. P. Singh
Abstract:
Recent experimental data on Bose-Einstein Condensation (BEC) of magnons in the spin-gap compound Yb$_2$Si$_2$O$_7$ revealed an asymmetric BEC dome arXiv:1810.13096v2. We examine modifications to the Heisenberg model on a breathing honeycomb lattice, showing that this physics can be explained by competing forms of weak anisotropy. We employ a gamut of analytical and numerical techniques to show tha…
▽ More
Recent experimental data on Bose-Einstein Condensation (BEC) of magnons in the spin-gap compound Yb$_2$Si$_2$O$_7$ revealed an asymmetric BEC dome arXiv:1810.13096v2. We examine modifications to the Heisenberg model on a breathing honeycomb lattice, showing that this physics can be explained by competing forms of weak anisotropy. We employ a gamut of analytical and numerical techniques to show that the anisotropy yields a field driven phase transition from a state with broken Ising symmetry to a phase which breaks no symmetries and crosses over to the polarized limit.
△ Less
Submitted 3 October, 2020; v1 submitted 22 January, 2020;
originally announced January 2020.
-
Confinement, reduced entanglement, and spin-glass order in a random quantum spin-ice model
Authors:
Anirudha Menon,
Tom Pardini,
Rajiv R. P. Singh
Abstract:
We study an effective spin model derived perturbatively from random transverse-field Ising model on the pyrochlore lattice. The model consists of spin-configurations on the pyrochlore lattice, restricted to the spin-ice subspace, with spins interacting with random Ising exchange couplings as well as ring exchanges along the hexagons of the lattice. This model is studied by exact diagonalization up…
▽ More
We study an effective spin model derived perturbatively from random transverse-field Ising model on the pyrochlore lattice. The model consists of spin-configurations on the pyrochlore lattice, restricted to the spin-ice subspace, with spins interacting with random Ising exchange couplings as well as ring exchanges along the hexagons of the lattice. This model is studied by exact diagonalization upto N=64 site systems. We calculate spin-glass correlation functions and local entanglement entropy $S_T$ between spins in a single tetrahedron and the rest of the system. We find that the model undergoes two phase transitions. At weak randomness the model is in a quantum spin-ice phase where $S_T=\ln{6}$. Increasing randomness first leads to a frozen phase, with long-range spin-glass order and $S_T=\ln{2}$ corresponding to the Cat states associated with Ising order. Further increase in randomness leads to a random resonating-hexagon phase with a frozen backbone of spins and a broad distribution of entanglement entropies. The implications of these studies for non-Kramers rare-earth pyrochlores are discussed.
△ Less
Submitted 22 November, 2019;
originally announced November 2019.
-
High-Field Expansion Approach to Kagome Antiferromagnets with Dzyaloshinskii-Moriya Interactions
Authors:
Michael O. Flynn,
Rajiv R. P. Singh
Abstract:
We apply linked cluster expansion techniques to study the polarized high-field phase of a spin-half antiferromagnet on the Kagome lattice with Heisenberg and Dzyaloshinskii-Moriya interactions (DMI). We find that the Dirac points of the single-magnon spectrum without DMI are robust against arbitrary DMI when the magnetic field lies in the Kagome plane. Unlike the typical case where DMI gaps the sp…
▽ More
We apply linked cluster expansion techniques to study the polarized high-field phase of a spin-half antiferromagnet on the Kagome lattice with Heisenberg and Dzyaloshinskii-Moriya interactions (DMI). We find that the Dirac points of the single-magnon spectrum without DMI are robust against arbitrary DMI when the magnetic field lies in the Kagome plane. Unlike the typical case where DMI gaps the spectrum, here we find that varying the DMI merely shifts the location of the Dirac points. In contrast, a magnetic field with a component out of the Kagome plane gaps the spectrum, leading to topological magnon bands. We map out a topological phase diagram as the couplings are varied by computing the band Chern numbers. A pair of phase transitions are observed and we find an enhanced thermal Hall conductivity near the phase boundary.
△ Less
Submitted 4 February, 2020; v1 submitted 21 May, 2019;
originally announced May 2019.
-
Robustness of entropy plateaus: A case study of triangular Ising antiferromagnets
Authors:
Owen Bradley,
Chunhan Feng,
Richard T. Scalettar,
Rajiv R. P. Singh
Abstract:
Residual entropy is a key feature associated with emergence in many-body systems. From a variety of frustrated magnets to the onset of spin-charge separation in Hubbard models and fermion-$Z_2$-flux variables in Kitaev models, the freezing of one set of degrees of freedom and establishment of local constraints are marked by a plateau in entropy as a function of temperature. Yet, with the exception…
▽ More
Residual entropy is a key feature associated with emergence in many-body systems. From a variety of frustrated magnets to the onset of spin-charge separation in Hubbard models and fermion-$Z_2$-flux variables in Kitaev models, the freezing of one set of degrees of freedom and establishment of local constraints are marked by a plateau in entropy as a function of temperature. Yet, with the exception of rare-earth pyrochlore family of spin-ice materials, evidence for such plateaus is rarely seen in real materials, raising questions about their robustness. Following recent experimental findings of the absence of such plateaus in triangular-lattice Ising antiferromagnet (TIAF) TmMgGaO$_4$ by Li et al, we explore in detail the existence and rounding of entropy plateaus in TIAF. We use a transfer matrix method to numerically calculate the properties of the system at different temperatures and magnetic fields, with further neighbor interactions and disorder. We find that temperature windows of entropy plateaus exist only when second-neighbor interactions are no more than a couple of percent of the nearest-neighbor ones, and they are also easily destroyed by disorder in the nearest-neighbor exchange variable, thereby explaining the challenge in observing such effects.
△ Less
Submitted 20 August, 2019; v1 submitted 23 January, 2019;
originally announced January 2019.
-
Griffiths-McCoy Singularities in the Dilute Transverse-Field Ising Model: A Numerical Linked Cluster Expansion Study
Authors:
Foster Thompson,
Rajiv R. P. Singh
Abstract:
We use Numerical Linked Cluster Expansions (NLCE) to study the site diluted transverse-field Ising model on the square-lattice at $T=0$. NLCE with a self-consistent mean-field on the boundary of the clusters is used to obtain the ground state magnetization, susceptibility, and structure factor as a function of transverse field $h$ and exchange constant $J$. Adding site-dilution to the model turns…
▽ More
We use Numerical Linked Cluster Expansions (NLCE) to study the site diluted transverse-field Ising model on the square-lattice at $T=0$. NLCE with a self-consistent mean-field on the boundary of the clusters is used to obtain the ground state magnetization, susceptibility, and structure factor as a function of transverse field $h$ and exchange constant $J$. Adding site-dilution to the model turns NLCE into a series expansion in the dilution parameter $p$. Studying the divergence of the structure factor allows us to establish the phase-diagram in the $h/J$ and $p$ plane. By studying the magnetization of the system in a longitudinal field, we investigate the Griffiths-McCoy (GM) singularities. We find that the magnetization develops non-linearities in the Griffiths phase with exponents that vary continuously with $h$. Additionally, the probability distribution of the local susceptibility develops long tails in the Griffith's phase, which is studied in terms of its moments.
△ Less
Submitted 13 November, 2018;
originally announced November 2018.
-
Incipient and well-developed entropy plateaus in spin-S Kitaev models
Authors:
J. Oitmaa,
A. Koga,
R. R. P. Singh
Abstract:
We present results on entropy and heat-capacity of the spin-S honeycomb-lattice Kitaev models using high-temperature series expansions and thermal pure quantum (TPQ) state methods. We study models with anisotropic couplings $J_z=1\ge J_x=J_y$ for spin values 1/2, 1, 3/2, and 2. We show that for $S>1/2$, any anisotropy leads to well developed plateaus in the entropy function at an entropy value of…
▽ More
We present results on entropy and heat-capacity of the spin-S honeycomb-lattice Kitaev models using high-temperature series expansions and thermal pure quantum (TPQ) state methods. We study models with anisotropic couplings $J_z=1\ge J_x=J_y$ for spin values 1/2, 1, 3/2, and 2. We show that for $S>1/2$, any anisotropy leads to well developed plateaus in the entropy function at an entropy value of $\frac{1}{2}\ln{2}$, independent of $S$. However, in the absence of anisotropy, there is an incipient entropy plateau at $S_{max}/2$, where $S_{max}$ is the infinite temperature entropy of the system. We discuss possible underlying microscopic reasons for the origin and implications of these entropy plateaus.
△ Less
Submitted 10 December, 2018; v1 submitted 26 September, 2018;
originally announced September 2018.
-
Rapid filling of the spin gap with temperature in the Schwinger-boson mean-field theory of the antiferromagnetic Heisenberg kagome model
Authors:
Jad C. Halimeh,
Rajiv R. P. Singh
Abstract:
Using Schwinger-boson mean-field theory, we calculate the dynamic spin structure factor at low temperatures $0<T\ll J$ for the spin-$1/2$ antiferromagnetic Heisenberg kagome model, within the gapped $\mathbb{Z}_2$ spin liquid phase Ansatz. We find that the spectral gap rapidly fills with temperature, with robust low-energy spectral weight developing by a temperature of $Δ/3$, where the spin gap is…
▽ More
Using Schwinger-boson mean-field theory, we calculate the dynamic spin structure factor at low temperatures $0<T\ll J$ for the spin-$1/2$ antiferromagnetic Heisenberg kagome model, within the gapped $\mathbb{Z}_2$ spin liquid phase Ansatz. We find that the spectral gap rapidly fills with temperature, with robust low-energy spectral weight developing by a temperature of $Δ/3$, where the spin gap is $2Δ$ (i.e., $Δ$ is the spinon gap), before any appreciable rise in spinon density or change in zero-temperature mean-field parameters. This is due to deconfinement of spinons which leads to terms suppressed only by $\exp(-Δ/T)$. At still higher temperatures, the spinon density increases rapidly leading to a breakdown of the Schwinger-boson mean-field approach. We suggest that if the impurity-free spectral functions can be obtained through neutron scattering experiments on kagome herbertsmithites, temperature dependence of the subgap weight can provide distinct signatures of a $\mathbb{Z}_2$ quantum spin liquid.
△ Less
Submitted 17 April, 2019; v1 submitted 18 July, 2018;
originally announced July 2018.
-
Structure Factors of the Kagome-Lattice Heisenberg antiferromagnets at finite temperatures
Authors:
Nicholas E. Sherman,
Rajiv R. P. Singh
Abstract:
We compute the real-space spin correlations and frquency and wave-vector resolved dynamic structure factors $S(\vec q,ω)$ for the nearest-neighbor Kagome-Lattice Heisenberg Antiferromagnet (KLHM) at finite temperatures using Numerical Linked Cluster Expansion (NLCE) method. A triangle-based NLCE is used to calculate frequency moments of the dynamic structure factors in the thermodynamic limit, whi…
▽ More
We compute the real-space spin correlations and frquency and wave-vector resolved dynamic structure factors $S(\vec q,ω)$ for the nearest-neighbor Kagome-Lattice Heisenberg Antiferromagnet (KLHM) at finite temperatures using Numerical Linked Cluster Expansion (NLCE) method. A triangle-based NLCE is used to calculate frequency moments of the dynamic structure factors in the thermodynamic limit, which show excellent convergence for $T> J/4$. A Gaussian approximation and the fluctuation-dissipation relation are used to to reconstruct the frequency dependence. We find that some features of the low temperature KLHM structure-factors begin to set in at temperatures of order $J$. Our results are in very good agreement with powder diffraction measurements reported earlier on the Herbertsmithite materials ZnCu$_3$(OH)$_6$Cl$_2$. However, the calculated properties differ from the low temperature ($T\approx J/100$) experimental measurements in one important regard. In line with the experimental observations, the spectral weight has a diffuse nature, which is predominantly spread along the extended Brillouin-Zone boundary. However, the maximum intensity is found in our calculations to be at the $K$ point of the extended Brillouin Zone in contrast to the low temperature experiments, where it is at the $M$ point. We suggest that experiments should be done at various temperatures to look for such a crossover of the maximum from the $K$ point to the $M$ point. In the absence of such a crossover, the Herbertsmithite materials must differ from KLHM in a significant manner.
△ Less
Submitted 14 November, 2017;
originally announced November 2017.
-
Efficient generation of series expansions for $\pm J$ Ising spin-glasses in a classical or a quantum (transverse) field
Authors:
R. R. P. Singh,
A. P. Young
Abstract:
We discuss generation of series expansions for Ising spin-glasses with a symmetric $\pm J$ (i.e. bimodal) distribution on d-dimensional hypercubic lattices using linked-cluster methods. Simplifications for the bimodal distribution allow us to go to higher order than for a general distribution. We discuss two types of problem, one classical and one quantum. The classical problem is that of the Isin…
▽ More
We discuss generation of series expansions for Ising spin-glasses with a symmetric $\pm J$ (i.e. bimodal) distribution on d-dimensional hypercubic lattices using linked-cluster methods. Simplifications for the bimodal distribution allow us to go to higher order than for a general distribution. We discuss two types of problem, one classical and one quantum. The classical problem is that of the Ising spin glass in a longitudinal magnetic field, $h$, for which we obtain high temperature series expansions in variables $\tanh(J/T)$ and $\tanh(h/T)$. The quantum problem is a $T=0$ study of the Ising spin glass in a transverse magnetic field $h_T$ for which we obtain a perturbation theory in powers of $J/h_T$. These methods require (i) enumeration and counting of \textit{all} connected clusters that can be embedded in the lattice up to some order $n$, and (ii) an evaluation of the contribution of each cluster for the quantity being calculated, known as the weight. We discuss a general method that takes the much smaller list (and count) of all no free-end (NFE) clusters on a lattice up to some order $n$, and automatically generates all other clusters and their counts up to the same order. The weights for finite clusters in both cases have a simple graphical interpretation that allows us to proceed efficiently for a general configuration of the $\pm J$ bonds, and at the end perform suitable disorder averaging. The order of our computations is limited by the weight calculations for the high-temperature expansions of the classical model, while they are limited by graph counting for the $T=0$ quantum system. Details of the calculational methods are presented.
△ Less
Submitted 11 September, 2017;
originally announced September 2017.
-
Exotic pairing symmetry of interacting Dirac fermions on a $π$ flux lattice
Authors:
Huaiming Guo,
Ehsan Khatami,
Yao Wang,
Thomas P. Devereaux,
Rajiv R. P. Singh,
Richard T. Scalettar
Abstract:
The pairing symmetry of interacting Dirac fermions on the $π$-flux lattice is studied with the determinant quantum Monte Carlo and numerical linked cluster expansion methods. The extended $s^*$- (i.e. extended $s$-) and d-wave pairing symmetries, which are distinct in the conventional square lattice, are degenerate under the Landau gauge. We demonstrate that the dominant pairing channel at strong…
▽ More
The pairing symmetry of interacting Dirac fermions on the $π$-flux lattice is studied with the determinant quantum Monte Carlo and numerical linked cluster expansion methods. The extended $s^*$- (i.e. extended $s$-) and d-wave pairing symmetries, which are distinct in the conventional square lattice, are degenerate under the Landau gauge. We demonstrate that the dominant pairing channel at strong interactions is an exotic $ds^*$-wave phase consisting of alternating stripes of $s^*$- and d-wave phases. A complementary mean-field analysis shows that while the $s^*$- and d-wave symmetries individually have nodes in the energy spectrum, the $ds^*$ channel is fully gapped. The results represent a new realization of pairing in Dirac systems, connected to the problem of chiral d-wave pairing on the honeycomb lattice, which might be more readily accessed by cold-atom experiments.
△ Less
Submitted 10 April, 2018; v1 submitted 23 August, 2017;
originally announced August 2017.
-
Principal Component Analysis for Fermionic Critical Points
Authors:
Natanael C. Costa,
Wenjian Hu,
Z. J. Bai,
Richard T. Scalettar,
Rajiv R. P. Singh
Abstract:
We use determinant Quantum Monte Carlo (DQMC), in combination with the principal component analysis (PCA) approach to unsupervised learning, to extract information about phase transitions in several of the most fundamental Hamiltonians describing strongly correlated materials. We first explore the zero temperature antiferromagnet to singlet transition in the Periodic Anderson Model, the Mott insul…
▽ More
We use determinant Quantum Monte Carlo (DQMC), in combination with the principal component analysis (PCA) approach to unsupervised learning, to extract information about phase transitions in several of the most fundamental Hamiltonians describing strongly correlated materials. We first explore the zero temperature antiferromagnet to singlet transition in the Periodic Anderson Model, the Mott insulating transition in the Hubbard model on a honeycomb lattice, and the magnetic transition in the 1/6-filled Lieb lattice. We then discuss the prospects for learning finite temperature superconducting transitions in the attractive Hubbard model, for which there is no sign problem. Finally, we investigate finite temperature charge density wave (CDW) transitions in the Holstein model, where the electrons are coupled to phonon degrees of freedom, and carry out a finite size scaling analysis to determine $T_c$. We examine the different behaviors associated with Hubbard-Stratonovich auxiliary field configurations on both the entire space-time lattice and on a single imaginary time slice, or other quantities, such as equal-time Green's and pair-pair correlation functions.
△ Less
Submitted 17 December, 2018; v1 submitted 15 August, 2017;
originally announced August 2017.
-
Critical and Griffiths-McCoy singularities in quantum Ising spin-glasses on d-dimensional hypercubic lattices: A series expansion study
Authors:
R. R. P. Singh,
A. P. Young
Abstract:
We study the $\pm J$ transverse-field Ising spin glass model at zero temperature on d-dimensional hypercubic lattices and in the Sherrington-Kirkpatrick (SK) model, by series expansions around the strong field limit. In the SK model and in high-dimensions our calculated critical properties are in excellent agreement with the exact mean-field results, surprisingly even down to dimension $d = 6$ whi…
▽ More
We study the $\pm J$ transverse-field Ising spin glass model at zero temperature on d-dimensional hypercubic lattices and in the Sherrington-Kirkpatrick (SK) model, by series expansions around the strong field limit. In the SK model and in high-dimensions our calculated critical properties are in excellent agreement with the exact mean-field results, surprisingly even down to dimension $d = 6$ which is below the upper critical dimension of $d=8$. In contrast, in lower dimensions we find a rich singular behavior consisting of critical and Griffiths-McCoy singularities. The divergence of the equal-time structure factor allows us to locate the critical coupling where the correlation length diverges, implying the onset of a thermodynamic phase transition. We find that the spin-glass susceptibility as well as various power-moments of the local susceptibility become singular in the paramagnetic phase $\textit{before}$ the critical point. Griffiths-McCoy singularities are very strong in two-dimensions but decrease rapidly as the dimension increases. We present evidence that high enough powers of the local susceptibility may become singular at the pure-system critical point.
△ Less
Submitted 31 October, 2017; v1 submitted 5 July, 2017;
originally announced July 2017.
-
High temperature thermodynamics of the honeycomb-lattice Kitaev-Heisenberg model: A high temperature series expansion study
Authors:
Rajiv R. P. Singh,
Jaan Oitmaa
Abstract:
We develop high temperature series expansions for the thermodynamic properties of the honeycomb-lattice Kitaev-Heisenberg model. Numerical results for uniform susceptibility, heat capacity and entropy as a function of temperature for different values of the Kitaev coupling $K$ and Heisenberg exachange coupling $J$ (with $|J|\le |K|$) are presented. These expansions show good convergence down to a…
▽ More
We develop high temperature series expansions for the thermodynamic properties of the honeycomb-lattice Kitaev-Heisenberg model. Numerical results for uniform susceptibility, heat capacity and entropy as a function of temperature for different values of the Kitaev coupling $K$ and Heisenberg exachange coupling $J$ (with $|J|\le |K|$) are presented. These expansions show good convergence down to a temperature of a fraction of $K$ and in some cases down to $T=K/10$. In the Kitaev exchange dominated regime, the inverse susceptibility has a nearly linear temperature dependence over a wide temperature range. However, we show that already at temperatures $10$-times the Curie-Weiss temperature, the effective Curie-Weiss constant estimated from the data can be off by a factor of 2. We find that the magnitude of the heat capacity maximum at the short-range order peak, is substantially smaller for small $J/K$ than for $J$ of order or larger than $K$. We suggest that this itself represents a simple marker for the relative importance of the Kitaev terms in these systems. Somewhat surprisingly, both heat capacity and susceptibility data on Na$_2$IrO$_3$ are consistent with a dominant {\it antiferromagnetic} Kitaev exchange constant of about $300-400$ $K$.
△ Less
Submitted 4 July, 2017;
originally announced July 2017.
-
de Almeida-Thouless instability in short-range Ising spin-glasses
Authors:
R. R. P. Singh,
A. P. Young
Abstract:
We use high temperature series expansions to study the $\pm J$ Ising spin-glass in a magnetic field in $d$-dimensional hypercubic lattices for $d=5, 6, 7$ and $8$, and in the infinite-range Sherrington-Kirkpatrick (SK) model. The expansions are obtained in the variable $w=\tanh^2{J/T}$ for arbitrary values of $u=\tanh^2{h/T}$ complete to order $w^{10}$. We find that the scaling dimension $Δ$ assoc…
▽ More
We use high temperature series expansions to study the $\pm J$ Ising spin-glass in a magnetic field in $d$-dimensional hypercubic lattices for $d=5, 6, 7$ and $8$, and in the infinite-range Sherrington-Kirkpatrick (SK) model. The expansions are obtained in the variable $w=\tanh^2{J/T}$ for arbitrary values of $u=\tanh^2{h/T}$ complete to order $w^{10}$. We find that the scaling dimension $Δ$ associated with the ordering-field $h^2$ equals $2$ in the SK model and for $d\ge 6$. However, in agreement with the work of Fisher and Sompolinsky, there is a violation of scaling in a finite field, leading to an anomalous $h$-$T$ dependence of the Almeida-Thouless (AT) line in high dimensions, while scaling is restored as $d \to 6$. Within the convergence of our series analysis, we present evidence supporting an AT line in $d\ge 6$. In $d=5$, the exponents $γ$ and $Δ$ are substantially larger than mean-field values, but we do not see clear evidence for the AT line in $d=5$.
△ Less
Submitted 17 July, 2017; v1 submitted 2 May, 2017;
originally announced May 2017.
-
Discovering Phases, Phase Transitions and Crossovers through Unsupervised Machine Learning: A critical examination
Authors:
Wenjian Hu,
Rajiv R. P. Singh,
Richard T. Scalettar
Abstract:
We apply unsupervised machine learning techniques, mainly principal component analysis (PCA), to compare and contrast the phase behavior and phase transitions in several classical spin models - the square and triangular-lattice Ising models, the Blume-Capel model, a highly degenerate biquadratic-exchange spin-one Ising (BSI) model, and the 2D XY model, and examine critically what machine learning…
▽ More
We apply unsupervised machine learning techniques, mainly principal component analysis (PCA), to compare and contrast the phase behavior and phase transitions in several classical spin models - the square and triangular-lattice Ising models, the Blume-Capel model, a highly degenerate biquadratic-exchange spin-one Ising (BSI) model, and the 2D XY model, and examine critically what machine learning is teaching us. We find that quantified principal components from PCA not only allow exploration of different phases and symmetry-breaking, but can distinguish phase transition types and locate critical points. We show that the corresponding weight vectors have a clear physical interpretation, which is particularly interesting in the frustrated models such as the triangular antiferromagnet, where they can point to incipient orders. Unlike the other well-studied models, the properties of the BSI model are less well known. Using both PCA and conventional Monte Carlo analysis, we demonstrate that the BSI model shows an absence of phase transition and macroscopic ground-state degeneracy. The failure to capture the `charge' correlations (vorticity) in the BSI model (XY model) from raw spin configurations points to some of the limitations of PCA. Finally, we employ a nonlinear unsupervised machine learning procedure, the `antoencoder method', and demonstrate that it too can be trained to capture phase transitions and critical points.
△ Less
Submitted 23 May, 2017; v1 submitted 31 March, 2017;
originally announced April 2017.
-
Cubic trihedral corner entanglement for a free scalar
Authors:
Lauren E. Hayward Sierens,
Pablo Bueno,
Rajiv R. P. Singh,
Robert C. Myers,
Roger G. Melko
Abstract:
We calculate the universal contribution to the $α$-Renyi entropy from a cubic trihedral corner in the boundary of the entangling region in 3+1 dimensions for a massless free scalar. The universal number, $v_α$, is manifest as the coefficient of a scaling term that is logarithmic in the size of the entangling region. Our numerical calculations find that this universal coefficient has both larger ma…
▽ More
We calculate the universal contribution to the $α$-Renyi entropy from a cubic trihedral corner in the boundary of the entangling region in 3+1 dimensions for a massless free scalar. The universal number, $v_α$, is manifest as the coefficient of a scaling term that is logarithmic in the size of the entangling region. Our numerical calculations find that this universal coefficient has both larger magnitude and the opposite sign to that induced by a smooth spherical entangling boundary in 3+1 dimensions, for which there is a well-known subleading logarithmic scaling. Despite these differences, up to the uncertainty of our finite-size lattice calculations, the functional dependence of the trihedral coefficient $v_α$ on the Rényi index $α$ is indistinguishable from that for a sphere, which is known analytically for a massless free scalar. We comment on the possible source of this $α$-dependence arising from the general structure of (3+1)-dimensional conformal field theories, and suggest calculations past the free scalar which could further illuminate the general structure of the trihedral divergence in the Rényi entropy.
△ Less
Submitted 9 March, 2017;
originally announced March 2017.
-
All-atom Molecular Dynamics Simulations of the Projection Domain of the Intrinsically Disordered htau40 Protein
Authors:
Natalie E. Sanborn,
N. Robert Hayre,
Rajiv R. P Singh,
Daniel L. Cox
Abstract:
We have performed all atom molecular dynamics simulations on the projection domain of the intrinsically disordered htau40 protein. After generating a suitable ensemble of starting conformations at high temperatures, at room temperature in an adaptive box algorithm we have generated histograms for the radius of gyration, secondary structure time series, generated model small angle x-ray scattering…
▽ More
We have performed all atom molecular dynamics simulations on the projection domain of the intrinsically disordered htau40 protein. After generating a suitable ensemble of starting conformations at high temperatures, at room temperature in an adaptive box algorithm we have generated histograms for the radius of gyration, secondary structure time series, generated model small angle x-ray scattering intensities, and model chemical shift plots for comparison to nuclear magnetic resonance data for solvated and filamentous tau. Significantly, we find that the chemical shift spectrum is more consistent with filamentous tau than full length solution based tau. We have also carried out principle component analysis and find three basics groups: compact globules, tadpoles, and extended hinging structures. To validate the adaptive box and our force field choice, we have run limited simulations in a large conventional box with varying force fields and find that our essential results are unchanged. We also performed two simulations with the TIP4P-D water model, the effects of which depended on whether the initial configuration was compact or extended.
△ Less
Submitted 30 October, 2016;
originally announced October 2016.
-
Nuclear relaxation rates in the Herbertsmithite Kagome antiferromagnets ZnCu3(OH)6Cl2
Authors:
Nicholas E. Sherman,
Takashi Imai,
Rajiv R. P. Singh
Abstract:
Local spectral functions and Nuclear Magnetic Relaxation (NMR) rates, 1/T1, for the spin-half Heisenberg antiferromagnet on the Kagome Lattice are calculated using Moriyas Gaussian approximation, as well as through an extrapolation of multiple frequency moments. The temperature dependence of the calculated rates is compared with the oxygen 1/T1 NMR data in Herbertsmithite. We find that the Gaussia…
▽ More
Local spectral functions and Nuclear Magnetic Relaxation (NMR) rates, 1/T1, for the spin-half Heisenberg antiferromagnet on the Kagome Lattice are calculated using Moriyas Gaussian approximation, as well as through an extrapolation of multiple frequency moments. The temperature dependence of the calculated rates is compared with the oxygen 1/T1 NMR data in Herbertsmithite. We find that the Gaussian approximation for 1/T1 shows behavior qualitatively similar to experiments with a sharp drop in rates at low temperatures, consistent with a spin-gapped behavior. However, this approximation significantly underestimates the magnitude of 1/T1 even at room temperature. Rates obtained from extrapolation of multiple frequency moments give very good agreement with the room temperature NMR data with J = 195 +/- 20 K and hyperfine couplings determined independently from other measurements. The use of multiple frequency moments also leads to additional low frequency weight in the local structure factors. The convergence of our calculations with higher frequency moments breaks down at low temperatures suggesting the existence of longer range dynamic correlations in the system despite the very short-range static correlations.
△ Less
Submitted 1 September, 2016;
originally announced September 2016.
-
Tricriticality in Crossed Ising Chains
Authors:
T. Cary,
R. R. P. Singh,
R. T. Scalettar
Abstract:
We explore the phase diagram of Ising spins on one-dimensional chains which criss-cross in two perpendicular directions and which are connected by interchain couplings. This system is of interest as a simpler, classical analog of a quantum Hamiltonian which has been proposed as a model of magnetic behavior in Nb$_{12}$O$_{29}$ and also, conceptually, as a geometry which is intermediate between one…
▽ More
We explore the phase diagram of Ising spins on one-dimensional chains which criss-cross in two perpendicular directions and which are connected by interchain couplings. This system is of interest as a simpler, classical analog of a quantum Hamiltonian which has been proposed as a model of magnetic behavior in Nb$_{12}$O$_{29}$ and also, conceptually, as a geometry which is intermediate between one and two dimensions. Using mean field theory as well as Metropolis Monte Carlo and Wang-Landau simulations, we locate quantitatively the boundaries of four ordered phases. Each becomes an effective Ising model with unique effective couplings at large interchain coupling. Away from this limit we demonstrate non-trivial critical behavior, including tricritical points which separate first and second order phase transitions. Finally, we present evidence that this model belongs to the 2D Ising universality class.
△ Less
Submitted 23 September, 2017; v1 submitted 9 August, 2016;
originally announced August 2016.
-
Spiral versus modulated collinear phases in the quantum ANNNH model
Authors:
Jaan Oitmaa,
Rajiv R. P. Singh
Abstract:
Motivated by the discovery of spiral and modulated collinear phases in several magnetic materials, we investigate the magnetic properties of Heisenberg spin S=1/2 antiferromagnets in 2 and 3 dimensions, with frustration arising from 2nd-neighbor couplings in one axial direction (the ANNNH model). Our results clearly demonstrate the presence of an incommensurate spiral phase at T=0 in 2 dimensions,…
▽ More
Motivated by the discovery of spiral and modulated collinear phases in several magnetic materials, we investigate the magnetic properties of Heisenberg spin S=1/2 antiferromagnets in 2 and 3 dimensions, with frustration arising from 2nd-neighbor couplings in one axial direction (the ANNNH model). Our results clearly demonstrate the presence of an incommensurate spiral phase at T=0 in 2 dimensions, extending to finite temperatures in 3 dimensions. The crossover between Néel and spiral order occurs at a value of the frustration parameter considerably above the classical value 0.25, a sign of substantial quantum fluctuations. We also investigate a possible modulated collinear phase with a wavelength of 4 lattice spacings, and find that it has substantially higher energy and hence is not realized in the model.
△ Less
Submitted 30 June, 2016;
originally announced June 2016.
-
Magnetic excitations and electronic interactions in Sr$_2$CuTeO$_6$: a spin-1/2 square lattice Heisenberg antiferromagnet
Authors:
P. Babkevich,
Vamshi M. Katukuri,
B. Fåk,
S. Rols,
T. Fennell,
D. Pajić,
H. Tanaka,
T. Pardini,
R. R. P. Singh,
A. Mitrushchenkov,
O. V. Yazyev,
H. M. Rønnow
Abstract:
Sr$_2$CuTeO$_6$ presents an opportunity for exploring low-dimensional magnetism on a square lattice of $S=1/2$ Cu$^{2+}$ ions. We employ ab initio multi-reference configuration interaction calculations to unravel the Cu$^{2+}$ electronic structure and to evaluate exchange interactions in Sr$_2$CuTeO$_6$. The latter results are validated by inelastic neutron scattering using linear spin-wave theory…
▽ More
Sr$_2$CuTeO$_6$ presents an opportunity for exploring low-dimensional magnetism on a square lattice of $S=1/2$ Cu$^{2+}$ ions. We employ ab initio multi-reference configuration interaction calculations to unravel the Cu$^{2+}$ electronic structure and to evaluate exchange interactions in Sr$_2$CuTeO$_6$. The latter results are validated by inelastic neutron scattering using linear spin-wave theory and series-expansion corrections for quantum effects to extract true coupling parameters. Using this methodology, which is quite general, we demonstrate that Sr$_2$CuTeO$_6$ is an almost realization of a nearest-neighbor Heisenberg antiferromagnet but with relatively weak coupling of 7.18(5) meV.
△ Less
Submitted 6 December, 2016; v1 submitted 31 May, 2016;
originally announced May 2016.
-
Competeing orders in spin-1 and spin-3/2 XXZ Kagome antiferromagnets: A series expansion study
Authors:
Jaan Oitmaa,
Rajiv R. P. Singh
Abstract:
We study the competition between $\sqrt{3} \times \sqrt{3}$ (RT3) and $q=0$ (Q0) magnetic orders in spin-one and spin-$3/2$ Kagome-lattice XXZ antiferromagnets with varying XY anisotropy parameter $Δ$, using series expansion methods. The Hamiltonian is split into two parts: an $H_0$ which favors the classical order in the desired pattern and an $H_1$, which is treated in perturbation theory by a s…
▽ More
We study the competition between $\sqrt{3} \times \sqrt{3}$ (RT3) and $q=0$ (Q0) magnetic orders in spin-one and spin-$3/2$ Kagome-lattice XXZ antiferromagnets with varying XY anisotropy parameter $Δ$, using series expansion methods. The Hamiltonian is split into two parts: an $H_0$ which favors the classical order in the desired pattern and an $H_1$, which is treated in perturbation theory by a series expansion. We find that the ground state energy series for the RT3 and Q0 phases are identical up to sixth order in the expansion, but ultimately a selection occurs, which depends on spin and the anisotropy $Δ$. Results for ground state energy and the magnetization are presented. These results are compared with recent spin-wave theory and coupled-cluster calculations. The series results for the phase diagram are close to the predictions of spin-wave theory. For the spin-one model at the Heisenberg point ($Δ=1$), our results are consistent with a vanishing order parameter, that is an absence of a magnetically ordered phase. We also develop series expansions for the ground state energy of the spin-one Heisenberg model in the trimerized phase. We find that the ground state energy in this phase is lower than those of magnetically ordered ones, supporting the existence of a spontaneously trimerized phase in this model.
△ Less
Submitted 28 October, 2015;
originally announced October 2015.
-
Nonzero temperature Entanglement Negativity of quantum spin models: Area law, Linked Cluster Expansions and Sudden Death
Authors:
Nicholas E. Sherman,
Trithep Devakul,
Matthew B. Hastings,
Rajiv R. P. Singh
Abstract:
We show that the bipartite logarithmic entanglement negativity (EN) of quantum spin models obeys an area law at all nonzero temperatures. We develop numerical linked cluster (NLC) expansions for the `area-law' logarithmic entanglement negativity as a function of temperature and other parameters. For one-dimensional models the results of NLC are compared with exact diagonalization on finite systems…
▽ More
We show that the bipartite logarithmic entanglement negativity (EN) of quantum spin models obeys an area law at all nonzero temperatures. We develop numerical linked cluster (NLC) expansions for the `area-law' logarithmic entanglement negativity as a function of temperature and other parameters. For one-dimensional models the results of NLC are compared with exact diagonalization on finite systems and are found to agree very well. The NLC results are also obtained for two dimensional XXZ and transverse-field Ising models. In all cases, we find a sudden onset (or sudden death) of negativity at a finite temperature above which the negativity is zero. We use perturbation theory to develop a physical picture for this sudden onset (or sudden death). The onset of EN or its magnitude are insensitive to classical finite-temperature phase transitions, supporting the argument for absence of any role of quantum mechanics at such transitions. On approach to a quantum critical point at $T=0$, negativity shows critical scaling in size and temperature.
△ Less
Submitted 27 October, 2015;
originally announced October 2015.
-
Unusual Corrections to Scaling and Convergence of Universal Renyi Properties at Quantum Critical Points
Authors:
Sharmistha Sahoo,
E. Miles Stoudenmire,
Jean-Marie Stéphan,
Trithep Devakul,
Rajiv R. P. Singh,
Roger G. Melko
Abstract:
At a quantum critical point, bipartite entanglement entropies have universal quantities which are subleading to the ubiquitous area law. For Renyi entropies, these terms are known to be similar to the von Neumann entropy, while being much more amenable to numerical and even experimental measurement. We show here that when calculating universal properties of Renyi entropies, it is important to acco…
▽ More
At a quantum critical point, bipartite entanglement entropies have universal quantities which are subleading to the ubiquitous area law. For Renyi entropies, these terms are known to be similar to the von Neumann entropy, while being much more amenable to numerical and even experimental measurement. We show here that when calculating universal properties of Renyi entropies, it is important to account for unusual corrections to scaling that arise from relevant local operators present at the conical singularity in the multi-sheeted Riemann surface. These corrections grow in importance with increasing Renyi index. We present studies of Renyi correlation functions in the 1+1 transverse-field Ising model (TFIM) using conformal field theory, mapping to free fermions, and series expansions, and the logarithmic entropy singularity at a corner in 2+1 for both free bosonic field theory and the TFIM, using numerical linked cluster expansions. In all numerical studies, accurate results are only obtained when unusual corrections to scaling are taken into account. In the worst case, an analysis ignoring these corrections can get qualitatively incorrect answers, such as predicting a decrease in critical exponents with the Renyi index, when they are actually increasing. We discuss a two-step extrapolation procedure that can be used to account for the unusual corrections to scaling.
△ Less
Submitted 1 September, 2015;
originally announced September 2015.
-
Early breakdown of area-law entanglement at many-body delocalization transition
Authors:
Trithep Devakul,
Rajiv R. P. Singh
Abstract:
We introduce the numerical linked cluster (NLC) expansion as a controlled numerical tool for the study of the many-body localization (MBL) transition in a disordered system with continuous non-perturbative disorder. Our approach works directly in the thermodynamic limit, in any spatial dimension, and does not rely on any finite size scaling procedure. We study the onset of many-body delocalization…
▽ More
We introduce the numerical linked cluster (NLC) expansion as a controlled numerical tool for the study of the many-body localization (MBL) transition in a disordered system with continuous non-perturbative disorder. Our approach works directly in the thermodynamic limit, in any spatial dimension, and does not rely on any finite size scaling procedure. We study the onset of many-body delocalization through the breakdown of area-law entanglement in a generic many-body eigenstate. By looking for initial signs of an instability of the localized phase, we obtain a value for the critical disorder, which we believe should be a lower bound for the true value, that is higher than current best estimates from finite size studies. This implies that most current methods tend to overestimate the extent of the localized phase due to finite size effects making the localized phase appear stable at small length scales. We also study the mobility edge in these systems as a function of energy density, and find that our conclusion is the same at all examined energies.
△ Less
Submitted 28 October, 2015; v1 submitted 19 August, 2015;
originally announced August 2015.
-
The many faces of quantum kagome materials: Interplay of further-neighbour exchange and Dzyaloshinskii-Moriya interaction
Authors:
Tsezar F. Seman,
Cheng-Chien Chen,
Rajiv R. P. Singh,
Michel van Veenendaal
Abstract:
The field of frustrated magnetism has been enriched significantly by the discovery of various kagome lattice compounds. These materials exhibit a great variety of macroscopic behaviours ranging from magnetic orders to quantum spin liquids. Using large-scale exact diagonalization, we construct the phase diagram of the $S=1/2$ $J_1$-$J_2$ kagome Heisenberg model with $z$-axis Dzyaloshinskii-Moriya i…
▽ More
The field of frustrated magnetism has been enriched significantly by the discovery of various kagome lattice compounds. These materials exhibit a great variety of macroscopic behaviours ranging from magnetic orders to quantum spin liquids. Using large-scale exact diagonalization, we construct the phase diagram of the $S=1/2$ $J_1$-$J_2$ kagome Heisenberg model with $z$-axis Dzyaloshinskii-Moriya interaction $D_z$. We show that this model can systematically account for many of the experimentally observed phases. Small $J_2$ and $D_z$ can stabilize respectively a gapped and a gapless spin liquid. When $J_2$ or $D_z$ is substantial, the ground state develops a $\mathbf{Q}=0$, $120^\textrm{o}$ antiferromagnetic order. The critical strengths for inducing magnetic transition are $D^c_z \sim 0.1\, J_1$ at $J_2=0$, and $J^c_2\sim 0.4\, J_1$ at $D_z=0$. The previously reported values of $D_z$ and $J_2$ for herbertsmithite [ZnCu$_3$(OH)$_6$Cl$_2$] place the compound in close proximity to a quantum critical point.
△ Less
Submitted 6 August, 2015;
originally announced August 2015.
-
Interplay of Magnetic Order, Pairing and Phase Separation in a One Dimensional Spin Fermion Model
Authors:
Wenjian Hu,
Richard T. Scalettar,
Rajiv R. P. Singh
Abstract:
We consider a lattice model of itinerant electrons coupled to an array of localized classical Heisenberg spins. The nature of the ground state ordered magnetic phases that result from the indirect spin-spin coupling mediated by the electrons is determined as a function of density and the spin-fermion coupling $J$. At a fixed chemical potential, spiral phases exist only up to values of $J$ which ar…
▽ More
We consider a lattice model of itinerant electrons coupled to an array of localized classical Heisenberg spins. The nature of the ground state ordered magnetic phases that result from the indirect spin-spin coupling mediated by the electrons is determined as a function of density and the spin-fermion coupling $J$. At a fixed chemical potential, spiral phases exist only up to values of $J$ which are less than roughly half the electronic bandwidth. At a fixed electron density and near half filling, the system phase separates into a half-filled antiferromagnetic phase and a spiral phase. The ferromagnetic phases are shown to be fully polarized, while the spiral phases have equal admixture of up and down spins. Phase separation survives in the presence of weak pairing field $Δ$ but disappears when $Δ$ exceeds a critical value $Δ_c$. If pairing fields are large enough, an additional spiral state arises at strong coupling $J$. The relevance of this study, especially the phase separation, to artificially engineered systems of adjacent itinerant electrons and localized spins is discussed. In particular, we propose a method which might allow for the braiding of Majorana fermions by changing the density and moving their location as they are pulled along by a phase separation boundary.
△ Less
Submitted 23 May, 2017; v1 submitted 15 June, 2015;
originally announced June 2015.
-
Are multiphase competition & order-by-disorder the keys to understanding Yb2Ti2O7?
Authors:
L. D. C. Jaubert,
Owen Benton,
Jeffrey G. Rau,
J. Oitmaa,
R. R. P. Singh,
Nic Shannon,
Michel J. P. Gingras
Abstract:
If magnetic frustration is most commonly known for undermining long-range order, as famously illustrated by spin liquids, the ability of matter to develop new collective mechanisms in order to fight frustration is no less fascinating, providing an avenue for the exploration and discovery of unconventional properties of matter. Here we study an ideal minimal model of such mechanisms which, incident…
▽ More
If magnetic frustration is most commonly known for undermining long-range order, as famously illustrated by spin liquids, the ability of matter to develop new collective mechanisms in order to fight frustration is no less fascinating, providing an avenue for the exploration and discovery of unconventional properties of matter. Here we study an ideal minimal model of such mechanisms which, incidentally, pertains to the perplexing quantum spin ice candidate Yb2Ti2O7. Specifically, we explain how thermal and quantum fluctuations, optimized by order-by-disorder selection, conspire to expand the stability region of an accidentally degenerate continuous symmetry U(1) manifold against the classical splayed ferromagnetic ground state that is displayed by the sister compound Yb2Sn2O7. The resulting competition gives rise to multiple phase transitions, in striking similitude with recent experiments on Yb2Ti2O7 [Lhotel et al., Phys. Rev. B 89 224419 (2014)]. Considering the effective Hamiltonian determined for Yb2Ti2O7, we provide, by combining a gamut of numerical techniques, compelling evidence that such multiphase competition is the long-sought missing key to understanding the intrinsic properties of this material. As a corollary, our work offers a pertinent illustration of the influence of chemical pressure in rare-earth pyrochlores.
△ Less
Submitted 20 May, 2015;
originally announced May 2015.
-
Finite-temperature superconducting correlations of the Hubbard model
Authors:
Ehsan Khatami,
Richard T. Scalettar,
Rajiv R. P. Singh
Abstract:
We utilize numerical linked-cluster expansions (NLCEs) and the determinantal quantum Monte Carlo algorithm to study pairing correlations in the square lattice Hubbard model. To benchmark the NLCE, we first locate the finite-temperature phase transition of the attractive model to a superconducting state away from half filling. We then explore the superconducting properties of the repulsive model fo…
▽ More
We utilize numerical linked-cluster expansions (NLCEs) and the determinantal quantum Monte Carlo algorithm to study pairing correlations in the square lattice Hubbard model. To benchmark the NLCE, we first locate the finite-temperature phase transition of the attractive model to a superconducting state away from half filling. We then explore the superconducting properties of the repulsive model for the d-wave and extended s-wave pairing symmetries. The pairing structure factor shows a strong tendency to d-wave pairing and peaks at an interaction strength comparable to the bandwidth. The extended s-wave structure factor and correlation length are larger at higher temperatures but clearly saturate as temperature is lowered, whereas the d-wave counterparts, which start off lower at high temperatures, continue to rise near half filling. This rise is even more dramatic in the d-wave susceptibility. The convergence of NLCEs breaks down as the susceptibilities and correlation lengths become large, so we are unable to determine the onset of long-range order. However, our results extend the conclusion, previously restricted to only magnetic and charge correlations, that NLCEs offer unique window into pairing in the Hubbard model at strong coupling.
△ Less
Submitted 12 July, 2015; v1 submitted 20 March, 2015;
originally announced March 2015.
-
Magnons and Excitation Continuum in XXZ triangular antiferromagnetic model: Application to $Ba_3CoSb_2O_9$
Authors:
E. A. Ghioldi,
A. Mezio,
L. O. Manuel,
R. R. P. Singh,
J. Oitmaa,
A. E. Trumper
Abstract:
We investigate the excitation spectrum of the triangular-lattice antiferromagnetic $XXZ$ model using series expansions and mean field Schwinger bosons approaches. The single-magnon spectrum computed with series expansions exhibits rotonic minima at the middle points of the edges of the Brillouin zone, for all values of the anisotropy parameter in the range $0\leq J^z/J\leq1$. Based on the good agr…
▽ More
We investigate the excitation spectrum of the triangular-lattice antiferromagnetic $XXZ$ model using series expansions and mean field Schwinger bosons approaches. The single-magnon spectrum computed with series expansions exhibits rotonic minima at the middle points of the edges of the Brillouin zone, for all values of the anisotropy parameter in the range $0\leq J^z/J\leq1$. Based on the good agreement with series expansions for the single-magnon spectrum, we compute the full dynamical magnetic structure factor within the mean field Schwinger boson approach to investigate the relevance of the $XXZ$ model for the description of the unusual spectrum found recently in $Ba_3CoSb_2O_9$. In particular, we obtain an extended continuum above the spin wave excitations, which is further enhanced and brought closer to those observed in $Ba_3CoSb_2O_9$ with the addition of a second neighbor exchange interaction approximately 15% of the nearest-neighbor value. Our results support the idea that excitation continuum with substantial spectral-weight are generically present in two-dimensional frustrated spin systems and fractionalization in terms of {\it bosonic} spinons presents an efficient way to describe them.
△ Less
Submitted 23 April, 2015; v1 submitted 19 November, 2014;
originally announced November 2014.
-
Spin Generation Via Bulk Spin Current in Three Dimensional Topological Insulators
Authors:
Xingyue Peng,
Yiming Yang,
Rajiv R. P. Singh,
Sergey Y. Savrasov,
Dong Yu
Abstract:
To date, spin generation in three-dimensional topological insulators is primarily modeled as a single-surface phenomenon, attributed to the momentum-spin locking on each individual surface. In this article we propose a mechanism of spin generation where the role of the insulating yet topologically non-trivial bulk becomes explicit: an external electric field creates a transverse pure spin current…
▽ More
To date, spin generation in three-dimensional topological insulators is primarily modeled as a single-surface phenomenon, attributed to the momentum-spin locking on each individual surface. In this article we propose a mechanism of spin generation where the role of the insulating yet topologically non-trivial bulk becomes explicit: an external electric field creates a transverse pure spin current through the bulk of a three-dimensional topological insulator, which transports spins between the top and bottom surfaces. Under sufficiently high surface disorder, the spin relaxation time can be extended via the Dyakonov-Perel mechanism. Consequently both the spin generation efficiency and surface conductivity are largely enhanced. Numerical simulation confirms that this spin generation mechanism originates from the unique topological connection of the top and bottom surfaces and is absent in other two dimensional systems such as graphene, even though they possess a similar Dirac cone-type dispersion.
△ Less
Submitted 6 January, 2016; v1 submitted 20 October, 2014;
originally announced October 2014.
-
Simulated Cytoskeletal Collapse via Tau Degradation
Authors:
Austin Sendek,
Henry R. Fuller,
N. Robert Hayre,
Rajiv R. P. Singh,
Daniel L. Cox
Abstract:
We present a coarse-grained two dimensional mechanical model for the microtubule-tau bundles in neuronal axons in which we remove taus, as can happen in various neurodegenerative conditions such as Alzheimer's disease, tauopathies, and chronic traumatic encephalopathy. Our simplified model includes (i) taus modeled as entropic springs between microtubules, (ii) removal of taus from the bundles due…
▽ More
We present a coarse-grained two dimensional mechanical model for the microtubule-tau bundles in neuronal axons in which we remove taus, as can happen in various neurodegenerative conditions such as Alzheimer's disease, tauopathies, and chronic traumatic encephalopathy. Our simplified model includes (i) taus modeled as entropic springs between microtubules, (ii) removal of taus from the bundles due to phosphorylation, and (iii) a possible depletion force between microtubules due to these dissociated phosphorylated taus. We equilibrate upon tau removal using steepest descent relaxation. In the absence of the depletion force, the transverse rigidity to radial compression of the bundle falls to zero at about 60% tau occupancy, in agreement with standard percolation theory results. However, with the attractive depletion force, spring removal leads to a first order collapse of the bundles over a wide range of tau occupancies for physiologically realizable conditions. While our simplest calculations assume a constant concentration of microtubule intercalants to mediate the depletion force, including a dependence that is linear in the detached taus yields the same collapse. Applying percolation theory to removal of taus at microtubule tips, which are likely to be the protective sites against dynamic instability, we argue that the microtubule instability can only obtain at low tau occupancy, from 0.06-0.30 depending upon the tau coordination at the microtubule tips. Hence, the collapse we discover is likely to be more robust over a wide range of tau occupancies than the dynamic instability. We suggest in vitro tests of our predicted collapse.
△ Less
Submitted 29 September, 2014;
originally announced September 2014.
-
Entanglement properties of the antiferromagnetic-singlet transition in the Hubbard model on bilayer square lattices
Authors:
Chia-Chen Chang,
Rajiv R. P. Singh,
Richard T. Scalettar
Abstract:
We calculate the bipartite \Renyi entanglement entropy of an $L\times L\times 2$ bilayer Hubbard model using a determinantal quantum Monte Carlo method recently proposed by Grover [Phys. Rev. Lett. {\bf 111}, 130402 (2013)]. Two types of bipartition are studied: (i) One that divides the lattice into two $L \times L$ planes, and (ii) One that divides the lattice into two equal-size (…
▽ More
We calculate the bipartite \Renyi entanglement entropy of an $L\times L\times 2$ bilayer Hubbard model using a determinantal quantum Monte Carlo method recently proposed by Grover [Phys. Rev. Lett. {\bf 111}, 130402 (2013)]. Two types of bipartition are studied: (i) One that divides the lattice into two $L \times L$ planes, and (ii) One that divides the lattice into two equal-size ($L\times L/2\times 2$) bilayers. We compare our calculations with those for the tight-binding model studied by the correlation matrix method. As expected, the entropy for bipartition (i) scales as $L^2$, while the latter scales with $L$ with possible logarithmic corrections. The onset of the antiferromagnet to singlet transition shows up by a saturation of the former to a maximal value and the latter to a small value in the singlet phase. We comment on the large uncertainties in the numerical results with increasing $U$, which would have to be overcome before the critical behavior and logarithmic corrections can be quantified.
△ Less
Submitted 17 September, 2014;
originally announced September 2014.
-
Ground state and spectral properties across a Charge Density Wave transition in a triangular-lattice spinless Fermion model
Authors:
Trithep Devakul,
Rajiv R. P. Singh
Abstract:
We study ground state properties and particle excitation spectra across a commensurate charge density wave transition in a system of strongly interacting fermions, using series expansion methods and mean-field theory. We consider a 1/3-filled system of spinless fermions on a triangular-lattice, with hopping parameter $t$, nearest-neighbor repulsion $V$, and a sublattice dependent chemical potentia…
▽ More
We study ground state properties and particle excitation spectra across a commensurate charge density wave transition in a system of strongly interacting fermions, using series expansion methods and mean-field theory. We consider a 1/3-filled system of spinless fermions on a triangular-lattice, with hopping parameter $t$, nearest-neighbor repulsion $V$, and a sublattice dependent chemical potential $μ_s$. The phase transition is found to be first order for $μ_s=0$, but becomes continuous with increasing $μ_s$. The particle and hole excitation spectra exhibit dramatic changes in the vicinity of the phase transitions and in the charge-density wave ordered phase. We discuss the relevance of this study to the Pinball Fermi liquid phase postulated theoretically in earlier studies as well as to various strongly correlated triangular-lattice materials.
△ Less
Submitted 17 September, 2014;
originally announced September 2014.
-
On the entanglement across a cubic interface in 3+1 dimensions
Authors:
Trithep Devakul,
Rajiv R. P. Singh
Abstract:
We calculate the area, edge and corner Renyi entanglement entropies in the ground state of the transverse-field Ising model, on a simple-cubic lattice, by high-field and low-field series expansions. We find that while the area term is positive and the line term is negative as required by strong subadditivity, the corner contributions are positive in 3-dimensions. Analysis of the series suggests th…
▽ More
We calculate the area, edge and corner Renyi entanglement entropies in the ground state of the transverse-field Ising model, on a simple-cubic lattice, by high-field and low-field series expansions. We find that while the area term is positive and the line term is negative as required by strong subadditivity, the corner contributions are positive in 3-dimensions. Analysis of the series suggests that the expansions converge up to the physical critical point from both sides. The leading area-law Renyi entropies match nicely from the high and low field expansions at the critical point, forming a sharp cusp there. We calculate the coefficients of the logarithmic divergence associated with the corner entropy and compare them with conformal field theory results with smooth interfaces and find a striking correspondence.
△ Less
Submitted 17 September, 2014; v1 submitted 30 June, 2014;
originally announced July 2014.
-
Quantum Critical Universality and Singular Corner Entanglement Entropy of Bilayer Heisenberg-Ising model
Authors:
Trithep Devakul,
Rajiv R. P. Singh
Abstract:
We consider a bilayer quantum spin model with anisotropic intra-layer exchange couplings. By varying the anisotropy, the quantum critical phenomena changes from XY to Heisenberg to Ising universality class, with two, three and one modes respectively becoming gapless simultaneously. We use series expansion methods to calculate the second and third Renyi entanglement entropies when the system is bip…
▽ More
We consider a bilayer quantum spin model with anisotropic intra-layer exchange couplings. By varying the anisotropy, the quantum critical phenomena changes from XY to Heisenberg to Ising universality class, with two, three and one modes respectively becoming gapless simultaneously. We use series expansion methods to calculate the second and third Renyi entanglement entropies when the system is bipartitioned into two parts. Leading area-law terms and subleading entropies associated with corners are separately calculated. We find clear evidence that the logarithmic singularity associated with the corners is universal in each class. Its coefficient along the Ising critical line is in excellent agreement with those obtained previously for the transverse-field Ising model. Our results provide strong evidence for the idea that the universal terms in the entanglement entropy provide an approximate measure of the low energy degrees of freedom in the system.
△ Less
Submitted 17 September, 2014; v1 submitted 1 June, 2014;
originally announced June 2014.
-
Magnetic correlations and pairing in the 1/5-depleted square lattice Hubbard model
Authors:
Ehsan Khatami,
Rajiv R. P. Singh,
Warren E. Pickett,
Richard T. Scalettar
Abstract:
We study the single-orbital Hubbard model on the 1/5-depleted square lattice geometry, which arises in such diverse systems as the spin-gap magnetic insulator CaV4O9 and ordered-vacancy iron selenides, presenting new issues regarding the origin of both magnetic ordering and superconductivity in these materials. We find a rich phase diagram that includes a plaquette singlet phase, a dimer singlet p…
▽ More
We study the single-orbital Hubbard model on the 1/5-depleted square lattice geometry, which arises in such diverse systems as the spin-gap magnetic insulator CaV4O9 and ordered-vacancy iron selenides, presenting new issues regarding the origin of both magnetic ordering and superconductivity in these materials. We find a rich phase diagram that includes a plaquette singlet phase, a dimer singlet phase, a Neel and a block-spin antiferromagnetic phase, and stripe phases. Quantum Monte Carlo simulations show that the dominant pairing correlations at half filling change character from d-wave in the plaquette phase to extended s-wave upon transition to the Neel phase. These findings have intriguing connections to iron-based superconductors, and suggest that some physics of multiorbital systems can be captured by a single-orbital model at different dopings.
△ Less
Submitted 25 September, 2014; v1 submitted 14 April, 2014;
originally announced April 2014.