Source author record

Wen-Jun Hu

Wen-Jun Hu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2020arXiv

Density Matrix Renormalization Group Study of Nematicity in Two Dimensions: Application to a Spin-$1$ Bilinear-Biquadratic Model on the Square Lattice

Nematic order is an exotic property observed in several strongly correlated systems, such as the iron-based superconductors. Using large-scale density matrix renormalization group (DMRG) techniques, we study at zero-temperature the nematic spin liquid that competes with spin dipolar and quadrupolar orders. We use these nematic orders to characterize different quantum phases and quantum phase transitions. More specifically, we study a spin-$1$ bilinear-biquadratic Heisenberg model on the square lattice with couplings beyond nearest neighbors. We focus on parameter regions around the highly symmetric $SU(3)$ point where the bilinear and biquadratic interactions are equal. With growing further-neighbor biquadratic interactions, we identify different spin dipolar and quadrupolar orders. We find that the DMRG results on cylindrical geometries correctly detect nematicity in different quantum states and accurately characterize the quantum phase transitions among them. Therefore, spin-driven nematicity -- here defined as the spontaneous breaking of the lattice invariance under a 90$^o$ rotation -- is an order parameter which can be studied directly in DMRG calculations in two dimensions in different quantum states.

preprint2020arXiv

Fractionalized Excitations Revealed by Entanglement Entropy

Fractionalized excitations develop in many unusual many-body states such as quantum spin liquids, disordered phases that cannot be described using any local order parameter. Because these exotic excitations correspond to emergent degrees of freedom, how to probe them and establish their existence is a long-standing challenge. We present a general procedure to reveal the fractionalized excitations using real-space entanglement entropy in critical spin liquids that are particularly relevant to experiments. Moreover, we show how to use the entanglement entropy to construct the corresponding spinon Fermi surface. Our work defines a new pathway to establish and characterize exotic excitations in novel quantum phases of matter.

preprint2020arXiv

Origin of the Magnetic and Orbital ordering in $α$-Sr$_2$CrO$_4$

Motivated by recent experimental progress in transition metal oxides with the K$_2$NiF$_4$ structure, we investigate the magnetic and orbital ordering in $α$-Sr$_2$CrO$_4$. Using first principles calculations, first we derive a three-orbital Hubbard model, which reproduces the {\it ab initio} band structure near the Fermi level. The unique reverse splitting of $t_{2g}$ orbitals in $α$-Sr$_2$CrO$_4$, with the $3d^2$ electronic configuration for the Cr$^{4+}$ oxidation state, opens up the possibility of orbital ordering in this material. Using real-space Hartree-Fock for multi-orbital systems, we constructed the ground-state phase diagram for the two-dimensional compound $α$-Sr$_2$CrO$_4$. We found stable ferromagnetic, antiferromagnetic, antiferro-orbital, and staggered orbital stripe ordering in robust regions of the phase diagram. Furthermore, using the density matrix renormalization group method for two-leg ladders with the realistic hopping parameters of $α$-Sr$_2$CrO$_4$, we explore magnetic and orbital ordering for experimentally relevant interaction parameters. Again, we find a clear signature of antiferromagnetic spin ordering along with antiferro-orbital ordering at moderate to large Hubbard interaction strength. We also explore the orbital-resolved density of states with Lanczos, predicting insulating behavior for the compound $α$-Sr$_2$CrO$_4$, in agreement with experiments. Finally, an intuitive understanding of the results is provided based on a hierarchy between orbitals, with $d_{xy}$ driving the spin order, while electronic repulsion and the effective one dimensionality of the movement within the $d_{xz}$ and $d_{yz}$ orbitals driving the orbital order.

preprint2020arXiv

Quantum Transitions of Nematic Phases in a Spin-$1$ Bilinear-Biquadratic Model and Their Implications for FeSe

Since its discovery, iron-based superconductivity has been known to develop near an antiferromagnetic order, but this paradigm fails in the iron chalcogenide FeSe, whose single-layer version holds the record for the highest superconducting transition temperature in the iron-based superconductors. The striking puzzle that FeSe displays nematic order (spontaneously broken lattice rotational symmetry) while being non-magnetic, has led to several competing proposals for its origin in terms of either the $3d$-electron's orbital degrees of freedom or spin physics in the form of frustrated magnetism. Here we argue that the phase diagram of FeSe under pressure could be qualitatively described by a quantum spin model with highly frustrated interactions. We implement both the site-factorized wave-function analysis and the large-scale density matrix renormalization group (DMRG) in cylinders to study the spin-$1$ bilinear-biquadratic model on the square lattice, and identify quantum transitions from the well-known $(π,0)$ antiferromagnetic state to an exotic $(π,0)$ antiferroquadrupolar order, either directly or through a $(π/2,π)$ antiferromagnetic state. These many phases, while distinct, are all nematic. We also discuss our theoretical ground-state phase diagram for the understanding of the experimental low-temperature phase diagram obtained by the NMR [P. S. Wang {\it et al.}, Phys. Rev. Lett. 117, 237001 (2016)] and X-ray scattering [K. Kothapalli {\it et al.}, Nature Communications 7, 12728 (2016)] measurements in pressurized FeSe. Our results suggest that superconductivity in a wide range of iron-based materials has a common origin in the antiferromagnetic correlations of strongly correlated electrons.

preprint2019arXiv

Deep Learning-Enhanced Variational Monte Carlo Method for Quantum Many-Body Physics

Artificial neural networks have been successfully incorporated into variational Monte Carlo method (VMC) to study quantum many-body systems. However, there have been few systematic studies of exploring quantum many-body physics using deep neural networks (DNNs), despite of the tremendous success enjoyed by DNNs in many other areas in recent years. One main challenge of implementing DNN in VMC is the inefficiency of optimizing such networks with large number of parameters. We introduce an importance sampling gradient optimization (ISGO) algorithm, which significantly improves the computational speed of training DNN in VMC. We design an efficient convolutional DNN architecture to compute the ground state of a one-dimensional (1D) SU($N$) spin chain. Our numerical results of the ground-state energies with up to 16 layers of DNN show excellent agreement with the Bethe-Ansatz exact solution. Furthermore, we also calculate the loop correlation function using the wave function obtained. Our work demonstrates the feasibility and advantages of applying DNNs to numerical quantum many-body calculations.

preprint2016arXiv

Frustrated magnetism and bicollinear antiferromagnetic order in FeTe

Iron chalcogenides display a rich variety of electronic orders in their phase diagram. A particularly enigmatic case is FeTe, a metal which possesses co-existing hole and electron Fermi surfaces as in the iron pnictides but has a distinct ($π$/2,$π/2$) bicollinear antiferromagnetic order in the Fe square lattice. While local-moment physics has been recognized as essential for understanding the electronic order, it has been a long-standing challenge to understand how the bicollinear antiferromagnetic ground state emerges in a proper quantum spin model. We show here that a bilinear-biquadratic spin-$1$ model on a square lattice with nonzero ring-exchange interactions exhibits the bicollinear antiferromagnetic order over an extended parameter space in its phase diagram. Our work shows that frustrated magnetism in the quantum spin model provides a unified description of the electronic orders in the iron chalcogenides and iron pnictides.

preprint2016arXiv

Spin liquid nature in the Heisenberg $J_{1}$-$J_{2}$ triangular antiferromagnet

We investigate the spin-$\frac{1}{2}$ Heisenberg model on the triangular lattice in the presence of nearest-neighbor $J_1$ and next-nearest-neighbor $J_2$ antiferromagnetic couplings. Motivated by recent findings from density-matrix renormalization group (DMRG) claiming the existence of a gapped spin liquid with signatures of spontaneously broken lattice point group symmetry [Zhu and White, Phys. Rev. B 92, 041105 (2015); Hu, Gong, Zhu, and Sheng, Phys. Rev. B 92, 140403 (2015)], we employ the variational Monte Carlo (VMC) approach to analyze the model from an alternative perspective that considers both magnetically ordered and paramagnetic trial states. We find a quantum paramagnet in the regime $0.08\lesssim J_2/J_1\lesssim 0.16$, framed by $120^{\circ}$ coplanar (stripe collinear) antiferromagnetic order for smaller (larger) $J_2/J_1$. By considering the optimization of spin-liquid wave functions of a different gauge group and lattice point group content as derived from Abrikosov mean-field theory, we obtain the gapless $U(1)$ Dirac spin liquid as the energetically most preferable state in comparison to all symmetric or nematic gapped $\mathbb{Z}_{2}$ spin liquids so far advocated by DMRG. Moreover, by the application of few Lanczos iterations, we find the energy to be the same as the DMRG result within error-bars. To further resolve the intriguing disagreement between VMC and DMRG, we complement our methodological approach by the pseudofermion functional renormalization group (PFFRG) to compare the spin structure factors for the paramagnetic regime calculated by VMC, DMRG, and PFFRG. This model promises to be an ideal test-bed for future numerical refinements in tracking the long-range correlations in frustrated magnets.

preprint2016arXiv

Variational Monte Carlo study of chiral spin liquid in quantum antiferromagnet on the triangular lattice

By using Gutzwiller projected fermionic wave functions and variational Monte Carlo technique, we study the spin-$1/2$ Heisenberg model with the first-neighbor ($J_1$), second-neighbor ($J_2$), and additional scalar chiral interaction $J_χ{\bf S}_i \cdot ({\bf S}_j \times {\bf S}_k)$ on the triangular lattice. In the non-magnetic phase of the $J_1-J_2$ triangular model with $0.08 \lesssim J_2/J_1 \lesssim 0.16$, recent density-matrix renormalization group (DMRG) studies [Zhu and White, Phys. Rev. B {\bf 92}, 041105 (2015); Hu, Gong, Zhu, and Sheng, Phys. Rev. B {\bf 92}, 140403 (2015)] find a possible gapped spin liquid with the signal of a competition between a chiral and a $Z_2$ spin liquid. Motivated by the DMRG results, we consider the chiral interaction $J_χ{\bf S}_i \cdot ({\bf S}_j \times {\bf S}_k)$ as a pertubation for this non-magnetic phase. We find that with growing $J_χ$, the gapless U(1) Dirac spin liquid, which has the best variational energy for $J_χ=0$, exhibits the energy instability towards a gapped spin liquid with non-trivial magnetic fluxes and nonzero chiral order. We calculate topological Chern number and ground-state degeneracy, both of which identify this flux state as the chiral spin liquid with fractionalized Chern number $C=1/2$ and two-fold topological degeneracy. Our results indicate a positive direction to stabilize a chiral spin liquid near the non-magnetic phase of the $J_1-J_2$ triangular model.

preprint2015arXiv

Competing Spin Liquid States in the Spin-$1/2$ Heisenberg Model On Triangular Lattice

We study the spin-$1/2$ Heisenberg model on the triangular lattice with the antiferromagnetic first ($J_1$) and second ($J_2$) nearest-neighbor interactions using density matrix renormalization group. By studying the spin correlation function, we find a $120^{\circ}$ magnetic order phase for $J_2 \lesssim 0.07 J_1$ and a stripe antiferromagnetic phase for $J_2 \gtrsim 0.15 J_1$. Between these two phases, we identify a spin liquid region characterized by the exponential decaying spin and dimer correlations, as well as the large spin singlet and triplet excitation gaps on finite-size systems. We find two near degenerating ground states with distinct properties in two sectors, which indicates more than one spin liquid candidates in this region. While the sector with spinon is found to respect the time reversal symmetry, the even sector without a spinon breaks such a symmetry for finite-size systems. Furthermore, we detect the signature of the fractionalization by following the evolution of different ground states with inserting spin flux into the cylinder system. Moreover, by tuning the anisotropic bond coupling, we explore the nature of the spin liquid phase and find the optimal parameter region for the gapped $Z_2$ spin liquid.

preprint2015arXiv

Variational Monte Carlo study of chiral spin liquid in the extended Heisenberg model on the Kagome lattice

We investigate the extended Heisenberg model on the Kagome lattice by using Gutzwiller projected fermionic states and the variational Monte Carlo technique. In particular, when both second- and third-neighbor super-exchanges are considered, we find that a gapped spin liquid described by non-trivial magnetic fluxes and long-range chiral-chiral correlations is energetically favored compared to the gapless U(1) Dirac state. Furthermore, the topological Chern number, obtained by integrating the Berry curvature, and the degeneracy of the ground state, by constructing linearly independent states, lead us to identify this flux state as the chiral spin liquid with $C=1/2$ fractionalized Chern number.

preprint2015arXiv

Variational Monte Carlo study of gapless spin liquid in the spin-$1/2$ XXZ antiferromagnetic model on the kagome lattice

By using the variational Monte Carlo technique, we study the spin-$1/2$ XXZ antiferromagnetic model (with easy-plane anisotropy) on the kagome lattice. A class of Gutzwiller projected fermionic states with a spin Jastrow factor is considered to describe either spin liquids (with $U(1)$ or $Z_2$ symmetry) or magnetically ordered phases (with ${\bf q}=(0,0)$ or ${\bf q}=(4π/3,0)$). We find that the magnetic states are not stable in the thermodynamic limit. Moreover, there is no energy gain to break the gauge symmetry from $U(1)$ to $Z_2$ within the spin-liquid states, as previously found in the Heisenberg model. The best variational wave function is therefore the $U(1)$ Dirac state, supplemented by the spin Jastrow factor. Furthermore, a vanishing $S=2$ spin gap is obtained at the variational level, in the whole regime from the $XY$ to the Heisenberg model.

preprint2014arXiv

Lanczos steps to improve variational wave functions

Gutzwiller-projected fermionic states can be efficiently implemented within quantum Monte Carlo calculations to define extremely accurate variational wave functions for Heisenberg models on frustrated two-dimensional lattices, not only for the ground state but also for low-energy excitations. The application of few Lanczos steps on top of these states further improves their accuracy, allowing calculations on large clusters. In addition, by computing both the energy and its variance, it is possible to obtain reliable estimations of exact results. Here, we report the cases of the frustrated Heisenberg models on square and Kagome lattices.

preprint2013arXiv

Direct evidence for a gapless $Z_2$ spin liquid by frustrating Néel antiferromagnetism

By direct calculations of the spin gap in the frustrated Heisenberg model on the square lattice, with nearest- ($J_1$) and next-nearest-neighbor ($J_2$) super-exchange couplings, we provide a solid evidence that the spin-liquid phase in the frustrated regime $0.45 \lesssim J_2/J_1 \lesssim 0.6$ is gapless. Our numerical method is based on a variational wave function that is {\it systematically} improved by the application of few Lanczos steps and allows us to obtain reliable extrapolations in the thermodynamic limit. The peculiar nature of the non-magnetic state is unveiled by the existence of S=1 gapless excitations at $k=(π,0)$ and $(0,π)$. The magnetic transition can be described and interpreted by a variational state that is built from Abrikosov fermions having a $Z_2$ gauge structure and four Dirac points in the spinon spectrum.

preprint2012arXiv

Absence of static stripes in the two-dimensional $t{-}J$ model by an accurate and systematic quantum Monte Carlo approach

We examine the two-dimensional $t{-}J$ model by using variational approach combined with well established quantum Monte Carlo techniques [S. Sorella {\it et al.}, \prl {\bf 88}, 117002 (2002)] that are used to improve systematically the accuracy of the variational ansatz. Contrary to recent density-matrix renormalization group and projected entangled-pair state calculations [P. Corboz {\it et al.}, \prb {\bf 84}, 041108(R) (2011)], a uniform phase is found for $J/t=0.4$, even when the calculation is biased with an ansatz that explicitly contains stripe order. Moreover, in the small hole doping regime, i.e., $δ\lesssim 0.1$, our results support the coexistence of antiferromagnetism and superconductivity.

preprint2009arXiv

Dynamical Mean Field Theory for the Bose-Hubbard Model

The dynamical mean field theory (DMFT), which is successful in the study of strongly correlated fermions, was recently extended to boson systems [Phys. Rev. B {\textbf 77}, 235106 (2008)]. In this paper, we employ the bosonic DMFT to study the Bose-Hubbard model which describes on-site interacting bosons in a lattice. Using exact diagonalization as the impurity solver, we get the DMFT solutions for the Green's function, the occupation density, as well as the condensate fraction on a Bethe lattice. Various phases are identified: the Mott insulator, the Bose-Einstein condensed (BEC) phase, and the normal phase. At finite temperatures, we obtain the crossover between the Mott-like regime and the normal phase, as well as the BEC-to-normal phase transition. Phase diagrams on the $μ/U-\tilde{t}/U$ plane and on the $T/U-\tilde{t}/U$ plane are produced ($\tilde{t}$ is the scaled hopping amplitude). We compare our results with the previous ones, and discuss the implication of these results to experiments.