No Arabic abstract
The exotic normal state of iron chalcogenide superconductor FeSe, which exhibits vanishing magnetic order and possesses an electronic nematic order, triggered extensive explorations of its magnetic ground state. To understand its novel properties, we study the ground state of a highly frustrated spin-$1$ system with bilinear-biquadratic interactions using unbiased large-scale density matrix renormalization group. Remarkably, with increasing biquadratic interactions, we find a paramagnetic phase between Neel and stripe magnetic ordered phases. We identify this phase as a candidate of nematic quantum spin liquid by the compelling evidences, including vanished spin and quadrupolar orders, absence of lattice translational symmetry breaking, and a persistent non-zero lattice nematic order in the thermodynamic limit. The established quantum phase diagram natually explains the observations of enhanced spin fluctuations of FeSe in neutron scattering measurement and the phase transition with increasing pressure. This identified paramagnetic phase provides a new possibility to understand the novel properties of FeSe.
Frustration in quantum spin systems promote a variety of novel quantum phases. An important example is the frustrated spin-$1$ model on the square lattice with the nearest-neighbor bilinear ($J_1$) and biquadratic ($K_1$) interactions. We provide strong evidence for a nematic spin liquid phase in a range of $K_1/J_1$ near the SU(3)-symmetric point ($J_1 = K_1$), based on the linear flavor-wave theory and extensive density matrix renormalization group calculation. This phase displays no spin dipolar or quadrupolar order, preserves translational symmetry but spontaneously breaks $C_4$ lattice rotational symmetry, and possesses fluctuations peaked at the wavevector $(pi, 2pi/3)$. The spin excitation gap drops rapidly with system size and appears to be gapless, and the nematic order is attributed to the dominant $(pi, 2pi/3)$ fluctuations. Our results provide a novel mechanism for electronic nematic order and, more generally, open up a new avenue to explore frustration-induced exotic ground states.
The interplay between the nematic order and magnetism in FeSe is not yet well understood. There is a controversy concerning the relationship between orbital and spin degrees of freedom in FeSe and their relevance for superconductivity. Here we investigate the effect of S substitution on the nematic transition temperature ($T_{rm n}$) and the low-energy spin fluctuations (SF) in FeSe single crystals. We show that the low-energy SF emerge below the nematic transition. The difference between the onset temperature for the critical SF ($T_{rm SF}$) and $T_{rm n}$ is small for FeSe but significantly increases with S substitution. Below $T_{rm SF}$ the Korringa relation is violated and the effective muon hyperfine coupling constant changes a sign. Our results exclude a direct coupling of the low-energy SF to the electronic nematic order indicating a presence of multiple spin degrees of freedom in FeSe$_{rm 1-x}$S$_{rm x}$.
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$-electrons 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 $(pi,0)$ antiferromagnetic state to an exotic $(pi,0)$ antiferroquadrupolar order, either directly or through a $(pi/2,pi)$ 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.
We use inelastic neutron scattering to study acoustic phonons and spin excitations in single crystals of NaFeAs, a parent compound of iron pnictide superconductors. NaFeAs exhibits a tetragonal-to-orthorhombic structural transition at $T_sapprox 58$ K and a collinear antiferromagnetic (AF) order at $T_Napprox 45$ K. While longitudinal and out-of-plane transverse acoustic phonons behave as expected, the in-plane transverse acoustic phonons reveal considerable softening on cooling to $T_s$, and then harden on approaching $T_N$ before saturating below $T_N$. In addition, we find that spin-spin correlation lengths of low-energy magnetic excitations within the FeAs layer and along the $c$-axis increase dramatically below $T_s$, and show weak anomaly across $T_N$. These results suggest that the electronic nematic phase present in the paramagnetic tetragonal phase is closely associated with dynamic spin-lattice coupling, possibly arising from the one-phonon-two-magnon mechanism.
We use the state-of-the-art tensor network state method, specifically, the finite projected entangled pair state (PEPS) algorithm, to simulate the global phase diagram of spin-$1/2$ $J_1$-$J_2$ Heisenberg model on square lattices up to $24times 24$. We provide very solid evidences to show that the nature of the intermediate nonmagnetic phase is a gapless quantum spin liquid (QSL), whose spin-spin and dimer-dimer correlations both decay with a power law behavior. There also exists a valence-bond solid (VBS) phase in a very narrow region $0.56lesssim J_2/J_1leq0.61$ before the system enters the well known collinear antiferromagnetic phase. We stress that our work gives rise to the first solid PEPS results beyond the well established density matrix renormalization group (DMRG) through one-to-one direct benchmark for small system sizes. Thus our numerical evidences explicitly demonstrate the huge power of PEPS for solving long-standing 2D quantum many-body problems. The physical nature of the discovered gapless QSL and potential experimental implications are also addressed.