We examine the magnetic correlations in quantum spin models that were derived recently as effective low-energy theories for electronic correlation effects on the edge states of graphene nanoribbons. For this purpose, we employ quantum Monte Carlo simulations to access the large-distance properties, accounting for quantum fluctuations beyond mean-field-theory approaches to edge magnetism. For certain chiral nanoribbons, antiferromagnetic interedge couplings were previously found to induce a gapped quantum disordered ground state of the effective spin model. We find that the extended nature of the intraedge couplings in the effective spin model for zigzag nanoribbons leads to a quantum phase transition at a large, finite value of the interedge coupling. This quantum critical point separates the quantum disordered region from a gapless phase of stable edge magnetism at weak intraedge coupling, which includes the ground states of spin-ladder models for wide zigzag nanoribbons. To study the quantum critical behavior, the effective spin model can be related to a model of two antiferromagnetically coupled Haldane-Shastry spin-half chains with long-ranged ferromagnetic intrachain couplings. The results for the critical exponents are compared also to several recent renormalization-group calculations for related long-ranged interacting quantum systems.
Three models of strongly correlated electron systems have been studied in this thesis: the Shastry-Sutherland quantum antiferromagnet, the quantum spin liquid emerging out of correlated Dirac fermions on the honeycomb lattice, and the edge-state magnetism in the graphene nanoribbons. The methods applied are quantum Monte Carlo simulations, namely, stochastic series expansion quantum Monte Carlo for bosonic (spin) systems and projector auxiliary field quantum Monte Carlo for fermionic systems. The stochastic series expansion quantum Monte Carlo method is applied to study the ground state properties as well as the magnetization process of the spin-1/2 easy-axis Heisenberg model on the Shastry-Sutherland lattice. The quantum phase transitions between a N\'eel state, a dimer triplet state and a ferromagnetic XY phase have been characterized by finite size scaling analysis. The crystalline structures of the triplon magnetic excitations in the magnetic plateaus have been identified. The projector auxiliary field quantum Monte Carlo is applied to explore the ground state phase diagram of the SU(2) Hubbard model on the honeycomb lattice. The interplay between the relativistic Dirac dispersion, the vanishing density of states at half-filling and the strong quantum fluctuations, drives the system from a semimetal at weak coupling to an antiferromagnetic Mott insulator in the strong coupling regime. In between these phases there emerges an unexpected quantum spin liquid with no symmetry breaking in any channel. After careful finite size scaling analysis of various correlation functions, the quantum spin liquid phase is characterized as a short-range resonating valence bond state with a finite spin gap to all magnetic excitations. The study of edge-state magnetism in graphene nanoribbons is more experimentally motivated. The projector auxiliary field quantum Monte Carlo method is used to measure the static spin-spin correlation function, the single-particle spectral function, as well as the dynamic spin structure factor along the zigzag edge of graphene nanoribbons. A dominant low energy peak generated by the Coulomb interaction in the single-particle spectral function is identified as the dynamical signature of edge-state magnetism. The magnetic excitations along the zigzag edge show a linear dispersion at low energies due to the effective antiferromagnetic coupling of the edges across the ribbon. The magnetic excitations are damped by scattering processes into the particle-hole continuum at high energies.
We investigate the low-lying excitation spectrum and ground-state properties of narrow graphene nanoribbons with zigzag edge configurations. Nanoribbons of comparable widths have been synthesized very recently [P. Ruffieux et al., Nature (London) 531, 489 (2016)], and their descriptions require more sophisticated methods since in this regime conventional methods, like mean-field or density-functional theory with local-density approximation, fail to capture the enhanced quantum fluctuations. Using the unbiased density-matrix renormalization-group algorithm, we calculate the charge gaps with high accuracy for different widths and interaction strengths and compare them with mean-field results. It turns out that the gaps are much smaller in the former case due to the proper treatment of quantum fluctuations. Applying the elements of quantum information theory, we also reveal the entanglement structure inside a ribbon and examine the spectrum of subsystem density matrices to understand the origin of entanglement. We examine the possibility of magnetic ordering and the effect of magnetic field. Our findings are relevant for understanding the gap values in different recent experiments and the deviations between them.
The interplay between the Kondo effect, indirect magnetic interaction, and geometrical frustration is studied in the Kondo lattice on the one-dimensional zigzag ladder. Using the density-matrix renormalization group, the ground-state and various short- and long-range spin- and density-correlation functions are calculated for the model at half filling as a function of the antiferromagnetic Kondo interaction down to
J=0.3t, where
t is the nearest-neighbor hopping on the zigzag ladder. Geometrical frustration is shown to lead to at least two critical points: Starting from the strong-
J limit, where almost local Kondo screening dominates and where the system is a nonmagnetic Kondo insulator, antiferromagnetic correlations between nearest-neighbor and next-nearest-neighbor local spins become stronger and stronger, until at ${J}_{\mathrm{c}}^{\mathrm{dim}}\ensuremath{\approx}0.89t$ frustration is alleviated by a spontaneous breaking of translational symmetry and a corresponding transition to a dimerized state. This is characterized by antiferromagnetic correlations along the legs and by alternating antiferro- and ferromagnetic correlations on the rungs of the ladder. A mechanism of partial Kondo screening that has been suggested for the Kondo lattice on the two-dimensional triangular lattice is not realized in the one-dimensional case. Furthermore, within the symmetry-broken dimerized state, there is a magnetic transition to a ${90}^{\ensuremath{\circ}}$ quantum spin spiral with quasi-long-range order at ${J}_{\mathrm{c}}^{\mathrm{mag}}\ensuremath{\approx}0.84t$. The quantum-critical point is characterized by a closure of the spin gap (with decreasing
J) and a divergence of the spin-correlation length and of the spin-structure factor
S(q) at wave vector $q=\ensuremath{\pi}/2$. This is opposed to the model on the one-dimensional bipartite chain, which is known to have a finite spin gap for all
J>0 at half filling.
We study the nature of edge states in extrinsically and spontaneously dimerized states of two-dimensional spin-
21 antiferromagnets, by performing quantum Monte Carlo simulation. We show that a gapless edge mode emerges in the wide region of the dimerized phases, and the critical exponent of spin correlators along the edge deviates from the value of Tomonaga-Luttinger liquid (TLL) universality in large but finite systems at low temperatures. We also demonstrate that the gapless nature at edges is stable against several perturbations such as external magnetic field, easy-plane XXZ anisotropy, Dzyaloshinskii-Moriya interaction, and further-neighbor exchange interactions. The edge states exhibit non-TLL behavior, depending strongly on model parameters and kinds of perturbations. Possible ways of detecting these edge states are discussed. Properties of edge states we show in this paper could also be used as reference points to study other edge states of more exotic gapped magnetic phases such as spin liquids.
Based on the effective spin-ladder model, we compute the negativity and geometric quantum discord (GQD) between two edge spins in zigzag graphene nanoribbons (ZGNRs) thermalized with a reservoir at temperature T (canonical ensemble). The results show that the entanglement only exists in antiferromagnetically coupled spin pairs, and completely vanishes between any two spins at the same ribbon edge. Intriguingly, the dominant entanglement always occurs in the nearest inter-edge coupled spin pairs, which is very robust even at room temperature. By contrast, almost all edge spin pairs can hold non-zero GQD, which strongly depends on the ribbon width. To get further insight into the properties of entanglement, the multiparticle entanglement (ME) shared between three spins is studied. The results show that ME always exists among spins located at different edges, and it is more robust than the corresponding two-spin entanglement in against the temperature. These remarkable properties make ZGNRs very promising for possible applications in quantum information processing.
The effect of next-nearest-neighbor hopping
t2 on the ground-state phase diagram of the one-dimensional Kondo lattice is studied with density-matrix renormalization-group techniques and by comparing with the phase diagram of the classical-spin variant of the same model. For a finite
t2, i.e., for a zigzag-ladder geometry, indirect antiferromagnetic interactions between the localized spins are geometrically frustrated. We demonstrate that
t2 at the same time triggers several magnetic phases which are absent in the model with nearest-neighbor hopping only. For strong
J, we find a transition from antiferromagnetic to incommensurate magnetic short-range order, which can be understood entirely in the classical-spin picture. For weaker
J, a spin-dimerized phase emerges, which spontaneously breaks the discrete translation symmetry. The phase is not accessible to perturbative means but is explained, on a qualitative level, by the classical-spin model as well. Spin dimerization alleviates magnetic frustration and is interpreted as a key to understand the emergence of quasi-long-range spiral magnetic order, which is found at weaker couplings. The phase diagram at weak
J, with gapless quasi-long-range order on top of the twofold degenerate spin-dimerized ground state, competing with a nondegenerate phase with gapped spin (and charge) excitations, is unconventional and eludes an effective low-energy spin-only theory.