Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2411.03872 · Physical Review Research

Electron correlation effects and spin-liquid state in the herbertsmithite kagome lattice

Sam Azadi, M. S. Bahramy, Thomas D. Kühne

Published Feb 14, 2025• 1 citation
arXivHTMLPDF

Abstract

We employ real-space variational quantum Monte Carlo methods with resonating valence bond many-body wave functions to investigate electron correlation effects in the kagome system Zn x Cu 4 − x O 6 . Using three trial wave functions of the Slater-Jastrow type, where (i) only the Jastrow correlation factor is optimized and the orbitals obtained by density functional theory, (ii) both the Jastrow factor and the Slater determinant are optimized, and (iii) additionally the Slater determinant is substituted by an antisymmetrized-geminal power-wave function, we analyze static and dynamic correlation energies across concentrations x = 0 , 1 3 , 2 3 , 1 . Optimizing the Slater determinant significantly enhances the correlation energy by approximately ∼ − 73 ( 4 ) mHa per electron. Eventually, the emergence of a quantum spin liquid state driven by a long-range correlation energy is discussed.

DOI
10.1103/physrevresearch.7.013165

Related papers

Correlation energy of the spin-polarized electron liquid studied using quantum Monte Carlo simulations

Sam Azadi, N. D. Drummond, Sam M. Vinko

Variational and diffusion quantum Monte Carlo (VMC and DMC) methods with Slater-Jastrow-backflow trial wave functions are used to study the spin-polarized three-dimensional uniform electron fluid. We report ground state VMC and DMC energies in the density range $0.5\ensuremath{\le}{r}_{\text{s}}\ensuremath{\le}20$. Finite-size errors are corrected using canonical-ensemble twist-averaged boundary conditions and extrapolation of the twist-averaged energy per particle calculated at three system sizes $(N=113, 259, \mathrm{and} 387)$ to the thermodynamic limit of infinite system size. The DMC energies in the thermodynamic limit are used to parametrize a local spin density approximation correlation function for inhomogeneous electron systems.

arXiv:2306.13426•Physical review. B./Physical review. B•Sep 15, 2023•2 citations

Novel quantum phases in correlated electron systems

Kyusung Hwang

The quantum nature of electrons emerges in various exotic ways when they are put together with quantum many-body interactions. In Mott insulating systems with strong electron correlation, frustration and strong quantum fluctuations promote localized electron spins to form nontrivial quantum paramagnetic phases such as quantum spin liquid and valence bond solid preventing the development of conventional magnetic orders. In this thesis, we study such exotic quantum phases in various different frustrated quantum magnets. First, encouraged by the spin liquid phenomenology observed in Ba3XSb2O9 (X=Cu, Ni), we investigate Z2 spin liquids on a three dimensional triangular lattice by employing the Schwinger boson mean-field theory with projective symmetry group analysis. It is found that a nematic spin liquid is energetically favorable on this lattice structure rather than symmetric spin liquids. Next, we study various valence bond solid phases that can appear in kagome lattice antiferromagnets by generalizing the recently proposed phenomenological Z2 gauge theory for the topological Z2 spin liquid state found in the nearest-neighbor kagome lattice Heisenberg model. This work provides a bridge between the topological spin liquid and valence bond solid phases previously proposed for kagome lattice antiferromagnets via vison condensation in the gauge theory. Then, we consider optical conductivity in Mott insulators as an experimental probe to detect quantum paramagnetic phases like valence bond solid. Applying this idea to deformed kagome lattice antiferromagnet Rb2Cu3SnF12 with a 12-site valence bond solid state, we estimate the size of optical conductivity in this quantum magnet. After this, another quantum paramagnet BiCu2PO6 is carefully investigated to clarify its puzzling ground state nature and quasiparticle excitations. Utilizing recent high resolution neutron scattering data, we construct a microscopic spin model and also an effective theory for the low energy spin excitations in the compound. We find significant influences of magnetic anisotropies on the ground state and quasiparticle dynamics. Lastly, we investigate emergent quantum phases in the thin film geometries of the pyrochlore iridates, where a number of exotic quantum ground states are proposed to occur in bulk materials as a result of the interplay between electron correlation and strong spin-orbit coupling. The fate of these bulk phases, as well as the novel quantum states that may arise only in a thin film setting, are studied via a theoretical model that allows layer-dependent magnetic structures.

•Nov 1, 2016•0 citations

Quantifying the phase diagram and Hamiltonian of S = 1/2 kagome antiferromagnets: bridging theory and experiment

Shengtao Jiang, Arthur C. Campello, Wei He +4 more

Abstract Spin-1/2 kagome antiferromagnets are leading candidates for realizing quantum spin liquid (QSL) ground states. While QSL ground states are predicted for the pure Heisenberg model, understanding the robustness of the QSL to additional interactions that may be present in real materials is a forefront question in the field. Here we employ large-scale density-matrix renormalization group simulations to investigate the effects of next-nearest neighbor exchange couplings J 2 and Dzyaloshinskii-Moriya interactions D , which are relevant to understanding the prototypical kagome materials herbertsmithite and Zn-barlowite. By utilizing clusters as large as XC12 and extrapolating the results to the thermodynamic limit, we precisely delineate the scope of the QSL phase, which remains robust across an expanded parameter range of J 2 and D . Direct comparison of the simulated static and dynamic spin structure factors with inelastic neutron scattering reveals the parameter space of the Hamiltonians for herbertsmithite and Zn-barlowite, and, importantly, provides compelling evidence that both materials exist within the QSL phase. These results establish a powerful convergence of theory and experiment in this most elusive state of matter.

arXiv:2504.07387•npj Computational Materials•Jan 15, 2026•1 citation

Improved wave function for strongly correlated electronic systems

Kenji Kobayashi, Kaoru Iguchi

A scheme of optimization of a trial wave function in the variational Monte Carlo calculation is developed. It provides a more precise description of the ground state for strongly correlated electronic systems with highly improved trial wave function. A general type of variational wave function constructed from a linear combination of Slater determinants multiplied by Jastrow correlation factors is proposed. They can explicitly contain low-energy fluctuations and long-range correlations. Using this optimization method, we examine the properties of these wave functions by applying them to a one-dimensional Hubbard model for which the exact solution is available for comparison. It is found that these wave functions prove to be substantial improvements compared to the conventional Jastrow-type wave functions.

Physical review. B, Condensed matter•Jan 15, 1993•14 citations

Variational Quantum Monte Carlo with Inclusion of Orbital Correlations

Shigenori Tanaka

Variational quantum Monte Carlo (VMC) method applied to interacting electron systems has played a pivotal role for the studies of electronic states of atoms, molecules and condensed matters. First, the many-body wave functions optimized by the VMC method provide good guiding functions in more accurate, diffusion quantum Monte Carlo (DMC) calculations for identical systems. Second, the VMC calculations themselves can give a good upper bound for the exact ground-state energy of many-electron systems when an appropriate form for the variational wave function is employed. In addition to a conventional form for the variational wave function that is expressed by the product of the Jastrow factor and the single Slater determinant, recent investigations have found that the expression in terms of many Slater determinants can provide a very accurate description of electronic states, especially for intrinsically multi-reference systems, which corresponds to the configuration interaction (CI) expansion method in conventional quantum chemistry. Although the CI method needs an extremely large number of expansion coefficients as variational parameters, the number of parameters can be reduced efficiently by using an idea of (molecular) orbital pair correlations. The purpose of the present paper is to examine the feasibility and usefulness of this idea within the framework of VMC method with the Jastrow–Slater wave function, where the ground state of a water molecule is investigated as a test case. In this study we rely on a variational wave function composed of a product of the Jastrow and Slater parts as 1⁄4 J S; ð1Þ

Journal of the Physical Society of Japan•Jun 10, 2013•3 citations

Static and dynamical spin correlations of the S = 1 2 random-bond antiferromagnetic Heisenberg model on the triangular and kagome lattices

Tokuro Shimokawa, Ken Watanabe, Hikaru Kawamura

Inspired by the recent theoretical suggestion that the random-bond $S=1/2$ antiferromagnetic Heisenberg model on the triangular and the kagome lattices might exhibit a randomness-induced quantum spin liquid (QSL) behavior when the strength of the randomness exceeds a critical value, and that this ``random-singlet state'' might be relevant to the QSL behaviors experimentally observed in triangular organic salts $\ensuremath{\kappa}\text{\ensuremath{-}}{(\mathrm{ET})}_{2}{\mathrm{Cu}}_{2}{(\mathrm{CN})}_{3}\phantom{\rule{4pt}{0ex}}\mathrm{and}\phantom{\rule{4pt}{0ex}}{\mathrm{EtMe}}_{3}\mathrm{Sb}{[\mathrm{Pd}{(\mathrm{dmit})}_{2}]}_{2}$ and in kagome herbertsmithite ${\mathrm{ZnCu}}_{3}{(\mathrm{OH})}_{6}{\mathrm{Cl}}_{2}$, we further investigate the nature of the static and the dynamical spin correlations of these models. We compute the static and the dynamical spin structure factors, $S(\mathbf{q})$ and $S(\mathbf{q},\ensuremath{\omega})$, by means of an exact diagonalization method. In both triangular and kagome models, the computed $S(\mathbf{q},\ensuremath{\omega})$ in the random-singlet state depends on the wave vector $\mathbf{q}$ only weakly, robustly exhibiting gapless behaviors accompanied by the broad distribution extending to higher energy $\ensuremath{\omega}$. Especially in the strongly random kagome model, $S(\mathbf{q},\ensuremath{\omega})$ hardly depends on $\mathbf{q}$, and exhibits an almost flat distribution for a wide range of $\ensuremath{\omega}$, together with a $\ensuremath{\omega}=0$ peak. These features agree semiquantitatively with the recent neutron-scattering data on a single-crystal herbertsmithite. Furthermore, the computed magnetization curve agrees almost quantitatively with the experimental one recently measured on a single-crystal herbertsmithite. These results suggest that the QSL state observed in herbertsmithite might indeed be the randomness-induced QSL state, i.e., the random-singlet state.

arXiv:1506.03576•Physical Review B•Oct 9, 2015•102 citations
We employ real-space variational quantum Monte Carlo methods with resonating valence bond many-body wave functions to investigate electron correlation effects in the kagome system Zn x Cu 4 − x O 6 . Using three trial wave functions of the Slater-Jastrow type, where (i) only the Jastrow correlation factor is optimized and the orbitals obtained by density functional theory, (ii) both the Jastrow factor and the Slater determinant are optimized, and (iii) additionally the Slater determinant is substituted by an antisymmetrized-geminal power-wave function, we analyze static and dynamic correlation energies across concentrations x = 0 , 1 3 , 2 3 , 1 . Optimizing the Slater determinant significantly enhances the correlation energy by approximately ∼ − 73 ( 4 ) mHa per electron. Eventually, the emergence of a quantum spin liquid state driven by a long-range correlation energy is discussed.
Variational and diffusion quantum Monte Carlo (VMC and DMC) methods with Slater-Jastrow-backflow trial wave functions are used to study the spin-polarized three-dimensional uniform electron fluid. We report ground state VMC and DMC energies in the density range $0.5\ensuremath{\le}{r}_{\text{s}}\ensuremath{\le}20$. Finite-size errors are corrected using canonical-ensemble twist-averaged boundary conditions and extrapolation of the twist-averaged energy per particle calculated at three system sizes (N=113,259,and387)(N=113, 259, \mathrm{and} 387)(N=113,259,and387) to the thermodynamic limit of infinite system size. The DMC energies in the thermodynamic limit are used to parametrize a local spin density approximation correlation function for inhomogeneous electron systems.
The quantum nature of electrons emerges in various exotic ways when they are put together with quantum many-body interactions. In Mott insulating systems with strong electron correlation, frustration and strong quantum fluctuations promote localized electron spins to form nontrivial quantum paramagnetic phases such as quantum spin liquid and valence bond solid preventing the development of conventional magnetic orders. In this thesis, we study such exotic quantum phases in various different frustrated quantum magnets. First, encouraged by the spin liquid phenomenology observed in Ba3XSb2O9 (X=Cu, Ni), we investigate Z2 spin liquids on a three dimensional triangular lattice by employing the Schwinger boson mean-field theory with projective symmetry group analysis. It is found that a nematic spin liquid is energetically favorable on this lattice structure rather than symmetric spin liquids. Next, we study various valence bond solid phases that can appear in kagome lattice antiferromagnets by generalizing the recently proposed phenomenological Z2 gauge theory for the topological Z2 spin liquid state found in the nearest-neighbor kagome lattice Heisenberg model. This work provides a bridge between the topological spin liquid and valence bond solid phases previously proposed for kagome lattice antiferromagnets via vison condensation in the gauge theory. Then, we consider optical conductivity in Mott insulators as an experimental probe to detect quantum paramagnetic phases like valence bond solid. Applying this idea to deformed kagome lattice antiferromagnet Rb2Cu3SnF12 with a 12-site valence bond solid state, we estimate the size of optical conductivity in this quantum magnet. After this, another quantum paramagnet BiCu2PO6 is carefully investigated to clarify its puzzling ground state nature and quasiparticle excitations. Utilizing recent high resolution neutron scattering data, we construct a microscopic spin model and also an effective theory for the low energy spin excitations in the compound. We find significant influences of magnetic anisotropies on the ground state and quasiparticle dynamics. Lastly, we investigate emergent quantum phases in the thin film geometries of the pyrochlore iridates, where a number of exotic quantum ground states are proposed to occur in bulk materials as a result of the interplay between electron correlation and strong spin-orbit coupling. The fate of these bulk phases, as well as the novel quantum states that may arise only in a thin film setting, are studied via a theoretical model that allows layer-dependent magnetic structures.
Abstract Spin-1/2 kagome antiferromagnets are leading candidates for realizing quantum spin liquid (QSL) ground states. While QSL ground states are predicted for the pure Heisenberg model, understanding the robustness of the QSL to additional interactions that may be present in real materials is a forefront question in the field. Here we employ large-scale density-matrix renormalization group simulations to investigate the effects of next-nearest neighbor exchange couplings J 2 and Dzyaloshinskii-Moriya interactions D , which are relevant to understanding the prototypical kagome materials herbertsmithite and Zn-barlowite. By utilizing clusters as large as XC12 and extrapolating the results to the thermodynamic limit, we precisely delineate the scope of the QSL phase, which remains robust across an expanded parameter range of J 2 and D . Direct comparison of the simulated static and dynamic spin structure factors with inelastic neutron scattering reveals the parameter space of the Hamiltonians for herbertsmithite and Zn-barlowite, and, importantly, provides compelling evidence that both materials exist within the QSL phase. These results establish a powerful convergence of theory and experiment in this most elusive state of matter.
A scheme of optimization of a trial wave function in the variational Monte Carlo calculation is developed. It provides a more precise description of the ground state for strongly correlated electronic systems with highly improved trial wave function. A general type of variational wave function constructed from a linear combination of Slater determinants multiplied by Jastrow correlation factors is proposed. They can explicitly contain low-energy fluctuations and long-range correlations. Using this optimization method, we examine the properties of these wave functions by applying them to a one-dimensional Hubbard model for which the exact solution is available for comparison. It is found that these wave functions prove to be substantial improvements compared to the conventional Jastrow-type wave functions.
Variational quantum Monte Carlo (VMC) method applied to interacting electron systems has played a pivotal role for the studies of electronic states of atoms, molecules and condensed matters. First, the many-body wave functions optimized by the VMC method provide good guiding functions in more accurate, diffusion quantum Monte Carlo (DMC) calculations for identical systems. Second, the VMC calculations themselves can give a good upper bound for the exact ground-state energy of many-electron systems when an appropriate form for the variational wave function is employed. In addition to a conventional form for the variational wave function that is expressed by the product of the Jastrow factor and the single Slater determinant, recent investigations have found that the expression in terms of many Slater determinants can provide a very accurate description of electronic states, especially for intrinsically multi-reference systems, which corresponds to the configuration interaction (CI) expansion method in conventional quantum chemistry. Although the CI method needs an extremely large number of expansion coefficients as variational parameters, the number of parameters can be reduced efficiently by using an idea of (molecular) orbital pair correlations. The purpose of the present paper is to examine the feasibility and usefulness of this idea within the framework of VMC method with the Jastrow–Slater wave function, where the ground state of a water molecule is investigated as a test case. In this study we rely on a variational wave function composed of a product of the Jastrow and Slater parts as 1⁄4 J S; ð1Þ
Inspired by the recent theoretical suggestion that the random-bond S=1/2S=1/2S=1/2 antiferromagnetic Heisenberg model on the triangular and the kagome lattices might exhibit a randomness-induced quantum spin liquid (QSL) behavior when the strength of the randomness exceeds a critical value, and that this ``random-singlet state'' might be relevant to the QSL behaviors experimentally observed in triangular organic salts $\ensuremath{\kappa}\text{\ensuremath{-}}{(\mathrm{ET})}_{2}{\mathrm{Cu}}_{2}{(\mathrm{CN})}_{3}\phantom{\rule{4pt}{0ex}}\mathrm{and}\phantom{\rule{4pt}{0ex}}{\mathrm{EtMe}}_{3}\mathrm{Sb}{[\mathrm{Pd}{(\mathrm{dmit})}_{2}]}_{2}$ and in kagome herbertsmithite ZnCu3(OH)6Cl2{\mathrm{ZnCu}}_{3}{(\mathrm{OH})}_{6}{\mathrm{Cl}}_{2}ZnCu3​(OH)6​Cl2​, we further investigate the nature of the static and the dynamical spin correlations of these models. We compute the static and the dynamical spin structure factors, S(q)S(\mathbf{q})S(q) and $S(\mathbf{q},\ensuremath{\omega})$, by means of an exact diagonalization method. In both triangular and kagome models, the computed $S(\mathbf{q},\ensuremath{\omega})$ in the random-singlet state depends on the wave vector q\mathbf{q}q only weakly, robustly exhibiting gapless behaviors accompanied by the broad distribution extending to higher energy $\ensuremath{\omega}$. Especially in the strongly random kagome model, $S(\mathbf{q},\ensuremath{\omega})$ hardly depends on q\mathbf{q}q, and exhibits an almost flat distribution for a wide range of $\ensuremath{\omega}$, together with a $\ensuremath{\omega}=0$ peak. These features agree semiquantitatively with the recent neutron-scattering data on a single-crystal herbertsmithite. Furthermore, the computed magnetization curve agrees almost quantitatively with the experimental one recently measured on a single-crystal herbertsmithite. These results suggest that the QSL state observed in herbertsmithite might indeed be the randomness-induced QSL state, i.e., the random-singlet state.