ترغب بنشر مسار تعليمي؟ اضغط هنا

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

201   0   0.0 ( 0 )
 نشر من قبل Wen-Jun Hu
 تاريخ النشر 2020
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

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.



قيم البحث

اقرأ أيضاً

320 - Annika Voll , Stefan Wessel 2014
We study thermodynamic properties as well as the dynamical spin and quadrupolar structure factors of the O(3)-symmetric spin-1 Heisenberg model with bilinear-biquadratic exchange interactions on the triangular lattice. Based on a sign-problem-free qu antum Monte Carlo approach, we access both the ferromagnetic and the ferroquadrupolar ordered, spin nematic phase as well as the SU(3)-symmetric point which separates these phases. Signatures of Goldstone soft-modes in the dynamical spin and the quadrupolar structure factors are identified, and the properties of the low-energy excitations are compared to the thermodynamic behavior observed at finite temperatures as well as to Schwinger-boson flavor-wave theory.
The ground state of a hole-doped t-t-J ladder with four legs favors a striped charge distribution. Spin excitation from the striped ground state is known to exhibit incommensurate spin excitation near q=(pi,pi) along the leg direction (qx direction). However, an outward dispersion from the incommensurate position toward q=(0,pi) is strong in intensity, inconsistent with inelastic neutron scattering (INS) experiment in hole-doped cuprates. Motivated by this inconsistency, we use the t-t-J model with m x n=96 lattice sites by changing lattice geometry from four-leg (24x4) to rectangle (12x8) shape and investigate the dynamical spin structure factor by using the dynamical density matrix renormalization group. We find that the outward dispersion has weak spectral weights in the 12x8 lattice, accompanied with the decrease of excitation energy close to q=(pi,pi), being consistent with the INS data. In the 12x8 lattice, weakening of incommensurate spin correlation is realized even in the presence of the striped charge distribution. For further investigation of geometry related spin dynamics, we focus on direction dependent spin excitation reported by recent resonant inelastic x-ray scattering (RIXS) for cuprate superconductors and obtain a consistent result with RIXS by examining an 8x8 t-t-J square lattice.
The channel-decomposed functional renormalization group (FRG) approach, most recently in the variant of truncated-unity-(TU-)FRG, has so far been used for various two-dimensional model systems. Yet, for many interesting material systems the third spa tial dimension is of clear relevance. Therefore FRG schemes working in three spatial dimensions (3D) are definitely on the wishlist. Here we demonstrate that a 3D TUFRG scheme can be set up in straightforward extension of previous 2D codes and gives physically sensible results with affordable numerical effort, both regarding the qualitative as well as the quantitative description. The computed phase diagram of the three-dimensional Hubbard model at half filling or perfect nesting shows a phase transition to a ((pi,pi,pi))-ordered antiferromagnetic ground state for repulsive interactions at an energy scale that compares well with other numerical approaches in the literature. Furthermore, the method allowed us to detect a (d)-wave pairing and a concurring ((pi,pi,0)) antiferromagnetic ground state in the hole doped Hubbard model.
We investigate the application of the Density Matrix Renormalization Group (DMRG) to the Hubbard model in momentum-space. We treat the one-dimensional models with dispersion relations corresponding to nearest-neighbor hopping and $1/r$ hopping and th e two-dimensional model with isotropic nearest-neighbor hopping. By comparing with the exact solutions for both one-dimensional models and with exact diagonalization in two dimensions, we first investigate the convergence of the ground-state energy. We find variational convergence of the energy with the number of states kept for all models and parameter sets. In contrast to the real-space algorithm, the accuracy becomes rapidly worse with increasing interaction and is not significantly better at half filling. We compare the results for different dispersion relations at fixed interaction strength over bandwidth and find that extending the range of the hopping in one dimension has little effect, but that changing the dimensionality from one to two leads to lower accuracy at weak to moderate interaction strength. In the one-dimensional models at half-filling, we also investigate the behavior of the single-particle gap, the dispersion of spinon excitations, and the momentum distribution function. For the single-particle gap, we find that proper extrapolation in the number of states kept is important. For the spinon dispersion, we find that good agreement with the exact forms can be achieved at weak coupling if the large momentum-dependent finite-size effects are taken into account for nearest-neighbor hopping. For the momentum distribution, we compare with various weak-coupling and strong-coupling approximations and discuss the importance of finite-size effects as well as the accuracy of the DMRG.
The real part of optical conductivity, $text{Re}sigma(omega)$, of the Mott insulators has a large amount of information on how spin and charge degrees of freedom interact with each other. By using the time-dependent density-matrix renormalization gro up, we study $text{Re}sigma(omega)$ of the two-dimensional Hubbard model on a square lattice at half filling. We find an excitonic peak at the Mott-gap edge of $text{Re}sigma(omega)$ not only for the two-dimensional square lattice but also for two- and four-leg ladders. For the square lattice, however, we do not clearly find a gap between an excitonic peak and continuum band, which indicates that a bound state is not well-defined. The emergence of an excitonic peak in $text{Re}sigma(omega)$ implies the formation of a spin polaron. Examining the dependence of $text{Re}sigma(omega)$ on the on-site Coulomb interaction and next-nearest neighbor hoppings, we confirm that an excitonic peak is generated from a magnetic effect. Electron scattering due to an electron-phonon interaction is expected to easily suppress an excitonic peak since spectral width of an excitonic peak is very narrow. Introducing a large broadening in $text{Re}sigma(omega)$ by modeling the electron-phonon coupling present in La$_{2}$CuO$_{4}$ and Nd$_{2}$CuO$_{4}$, we obtain $text{Re}sigma(omega)$ comparable with experiments.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا