-
Phase diagram of the antiferromagnetic $J_1$-$J_2$ spin-$1$ pyrochlore Heisenberg model
Authors:
Imre Hagymási,
Nils Niggemann,
Johannes Reuther
Abstract:
We study the phase diagram of the antiferromagnetic $J_1$-$J_2$ Heisenberg model on the pyrochlore lattice with $S=1$ spins at zero and finite temperatures. We use a combination of complementary state-of-the-art quantum many-body approaches such as density matrix renormalization group (DMRG), density-matrix purification and pseudo-Majorana functional renormalization group (PMFRG). We present an ef…
▽ More
We study the phase diagram of the antiferromagnetic $J_1$-$J_2$ Heisenberg model on the pyrochlore lattice with $S=1$ spins at zero and finite temperatures. We use a combination of complementary state-of-the-art quantum many-body approaches such as density matrix renormalization group (DMRG), density-matrix purification and pseudo-Majorana functional renormalization group (PMFRG). We present an efficient approach to preserve the applicability of the PMFRG for spin-1 systems at finite temperatures despite the inevitable presence of unphysical spin states. The good performance of our methods is first demonstrated for the nearest-neighbor pyrochlore Heisenberg model where the finite temperature behavior of the specific heat and uniform susceptibility show excellent agreement within PMFRG and density-matrix purification. Including an antiferromagnetic second neighbor coupling we find that the non-magnetic ground-state phase of the nearest neighbor model extents up to $J_2/J_1 \sim 0.02 $ within DMRG, beyond which magnetic $\boldsymbol{k}=0$ long-range order sets in. Our PMFRG calculations find the phase transition in a similar regime $J_2/J_1\sim 0.035(8)$ which, together with the DMRG result, provides a strong argument for the existence of a small but finite non-magnetic ground-state phase in the spin-1 pyrochlore Heisenberg model. We also discuss the origin of discrepancies between different versions of the functional renormalization group concerning the location of this phase transition.
△ Less
Submitted 21 May, 2024;
originally announced May 2024.
-
Enhanced symmetry-breaking tendencies in the $S=1$ pyrochlore antiferromagnet
Authors:
Imre Hagymási,
Vincent Noculak,
Johannes Reuther
Abstract:
We investigate the ground-state properties of the nearest-neighbor $S=1$ pyrochlore Heisenberg antiferromagnet using two complementary numerical methods, density-matrix renormalization group (DMRG) and pseudofermion functional renormalization group (PFFRG). Within DMRG, we are able to reliably study clusters with up to 48 spins by keeping 20 000 SU(2) states. The investigated 32-site and 48-site c…
▽ More
We investigate the ground-state properties of the nearest-neighbor $S=1$ pyrochlore Heisenberg antiferromagnet using two complementary numerical methods, density-matrix renormalization group (DMRG) and pseudofermion functional renormalization group (PFFRG). Within DMRG, we are able to reliably study clusters with up to 48 spins by keeping 20 000 SU(2) states. The investigated 32-site and 48-site clusters both show indications of a robust $C_3$ rotation symmetry breaking of the ground-state spin correlations and the 48-site cluster additionally features inversion symmetry breaking. Our PFFRG analysis of various symmetry-breaking perturbations corroborates the findings of either $C_3$ or a combined $C_3$/inversion symmetry breaking. Moreover, in both methods the symmetry-breaking tendencies appear to be more pronounced than in the $S=1/2$ system.
△ Less
Submitted 4 July, 2022;
originally announced July 2022.
-
Observation of competing, correlated ground states in the flat band of rhombohedral graphite
Authors:
Imre Hagymási,
Mohammad Syahid Mohd Isa,
Zoltán Tajkov,
Krisztián Márity,
Oroszlány László,
János Koltai,
Assem Alassaf,
Péter Kun,
Konrád Kandrai,
András Pálinkás,
Péter Vancsó,
Levente Tapasztó,
Péter Nemes-Incze
Abstract:
In crystalline solids the interactions of charge and spin can result in a variety of emergent quantum ground states, especially in partially filled, topological flat bands such as Landau levels or in 'magic-angle' bilayer graphene. Much less explored is rhombohedral graphite (RG), perhaps the simplest and structurally most perfect condensed matter system to host a flat band protected by symmetry.…
▽ More
In crystalline solids the interactions of charge and spin can result in a variety of emergent quantum ground states, especially in partially filled, topological flat bands such as Landau levels or in 'magic-angle' bilayer graphene. Much less explored is rhombohedral graphite (RG), perhaps the simplest and structurally most perfect condensed matter system to host a flat band protected by symmetry. By scanning tunneling microscopy we map the flat band charge density of 8, 10 and 17 layers and identify a domain structure emerging from a competition between a sublattice antiferromagnetic insulator and a gapless correlated paramagnet. Our density-matrix renormalization group calculations explain the observed features and demonstrate that the correlations are fundamentally different from graphene based magnetism identified until now, forming the ground state of a quantum magnet. Our work establishes RG as a new platform to study many-body interactions beyond the mean-field approach, where quantum fluctuations and entanglement dominate.
△ Less
Submitted 15 July, 2022; v1 submitted 26 January, 2022;
originally announced January 2022.
-
Magnetization process and ordering of the $S=1/2$ pyrochlore magnet in a field
Authors:
Imre Hagymási,
Robin Schäfer,
Roderich Moessner,
David J. Luitz
Abstract:
We study the $S=1/2$ pyrochlore Heisenberg antiferromagnet in a magnetic field. Using large scale density-matrix renormalization group (DMRG) calculations for clusters up to $128$ spins, we find indications for a finite triplet gap, causing a threshold field to nonzero magnetization in the magnetization curve. We obtain a robust saturation field consistent with a magnon crystal, although the corre…
▽ More
We study the $S=1/2$ pyrochlore Heisenberg antiferromagnet in a magnetic field. Using large scale density-matrix renormalization group (DMRG) calculations for clusters up to $128$ spins, we find indications for a finite triplet gap, causing a threshold field to nonzero magnetization in the magnetization curve. We obtain a robust saturation field consistent with a magnon crystal, although the corresponding $5/6$ magnetization plateau is very slim and possibly unstable. Most remarkably, there is a pronounced and apparently robust $\frac 1 2$ magnetization plateau where the groundstate breaks (real-space) rotational symmetry, exhibiting oppositely polarised spins on alternating kagome and triangular planes. Reminiscent of the kagome ice plateau of the pyrochlore Ising magnet known as spin ice, it arises via a much more subtle `quantum order by disorder' mechanism.
△ Less
Submitted 13 September, 2022; v1 submitted 17 June, 2021;
originally announced June 2021.
-
Possible inversion symmetry breaking in the $S=1/2$ pyrochlore Heisenberg magnet
Authors:
Imre Hagymási,
Robin Schäfer,
Roderich Moessner,
David J. Luitz
Abstract:
We address the ground-state properties of the long-standing and much-studied three-dimensional quantum spin liquid candidate, the $S=\frac 1 2$ pyrochlore Heisenberg antiferromagnet. By using $SU(2)$ density-matrix renormalization group (DMRG), we are able to access cluster sizes of up to 128 spins. Our most striking finding is a robust spontaneous inversion symmetry breaking, reflected in an ener…
▽ More
We address the ground-state properties of the long-standing and much-studied three-dimensional quantum spin liquid candidate, the $S=\frac 1 2$ pyrochlore Heisenberg antiferromagnet. By using $SU(2)$ density-matrix renormalization group (DMRG), we are able to access cluster sizes of up to 128 spins. Our most striking finding is a robust spontaneous inversion symmetry breaking, reflected in an energy density difference between the two sublattices of tetrahedra, familiar as a starting point of earlier perturbative treatments. We also determine the ground-state energy, $E_0/N_\text{sites} = -0.490(6) J$, by combining extrapolations of DMRG with those of a numerical linked cluster expansion. These findings suggest a scenario in which a finite-temperature spin liquid regime gives way to a symmetry-broken state at low temperatures.
△ Less
Submitted 19 March, 2021; v1 submitted 7 October, 2020;
originally announced October 2020.
-
The pyrochlore S=1/2 Heisenberg antiferromagnet at finite temperature
Authors:
Robin Schäfer,
Imre Hagymási,
Roderich Moessner,
David J. Luitz
Abstract:
Frustrated three dimensional quantum magnets are notoriously impervious to theoretical analysis. Here we use a combination of three computational methods to investigate the three dimensional pyrochlore $S=1/2$ quantum antiferromagnet, an archetypical frustrated magnet, at finite temperature, $T$: canonical typicality for a finite cluster of $2\times 2 \times 2$ unit cells (i.e. $32$ sites), a fini…
▽ More
Frustrated three dimensional quantum magnets are notoriously impervious to theoretical analysis. Here we use a combination of three computational methods to investigate the three dimensional pyrochlore $S=1/2$ quantum antiferromagnet, an archetypical frustrated magnet, at finite temperature, $T$: canonical typicality for a finite cluster of $2\times 2 \times 2$ unit cells (i.e. $32$ sites), a finite-$T$ matrix product state method on a larger cluster with $48$ sites, and the numerical linked cluster expansion (NLCE) using clusters up to $25$ lattice sites, which include non-trivial hexagonal and octagonal loops. We focus on thermodynamic properties (energy, specific heat capacity, entropy, susceptibility, magnetisation) next to the static structure factor. We find a pronounced maximum in the specific heat at $T = 0.57 J$, which is stable across finite size clusters and converged in the series expansion. This is well-separated from a residual amount of spectral weight of $0.47 k_B \ln 2$ per spin which has not been released even at $T\approx0.25 J$, the limit of convergence of our results. This is a large value compared to a number of highly frustrated models and materials, such as spin ice or the kagome $S=1/2$ Heisenberg antiferromagnet. We also find a non-monotonic dependence on $T$ of the magnetisation at low magnetic fields, reflecting the dominantly non-magnetic character of the low-energy spectral weight. A detailed comparison of our results to measurements for the $S=1$ material NaCaNi$_2$F$_7$ yields rough agreement of the functional form of the specific heat maximum, which in turn differs from the sharper maximum of the heat capacity of the spin ice material Dy$_2$Ti$_2$O$_7$, all of which are yet qualitatively distinct from conventional, unfrustrated magnets.
△ Less
Submitted 10 March, 2020;
originally announced March 2020.
-
Stability of edge magnetism against disorder in zigzag MoS$_2$ nanoribbons
Authors:
Péter Vancsó,
Imre Hagymási,
Pauline Castenetto,
Philippe Lambin
Abstract:
Molybdenum disulfide nanoribbons with zigzag edges show ferromagnetic and metallic properties based on previous \emph{ab-initio} calculations. The investigation of the role of disorder on the magnetic properties is, however, still lacking due to the computational costs of these methods. In this work we fill this gap by studying the magnetic and electronic properties of several nanometer long MoS…
▽ More
Molybdenum disulfide nanoribbons with zigzag edges show ferromagnetic and metallic properties based on previous \emph{ab-initio} calculations. The investigation of the role of disorder on the magnetic properties is, however, still lacking due to the computational costs of these methods. In this work we fill this gap by studying the magnetic and electronic properties of several nanometer long MoS$_2$ zigzag nanoribbons using tight-binding and Hubbard Hamiltonians. Our results reveal that proper tight-binding parameters for the edge atoms are crucial to obtain quantitatively the metallic states and the magnetic properties of MoS$_2$ nanoribbons. With the help of the fine-tuned parameters, we perform large-scale calculations and predict the spin domain-wall energy along the edges, which is found to be significantly lower compared to that of the zigzag graphene nanoribbons. The tight-binding approach allows us to address the effect of edge disorder on the magnetic properties. Our results open the way for investigating electron-electron effects in realistic-size nanoribbon devices in MoS$_2$ and also provide valuable information for spintronic applications.
△ Less
Submitted 3 June, 2019;
originally announced June 2019.
-
Dynamical topological quantum phase transitions in nonintegrable models
Authors:
I. Hagymási,
C. Hubig,
Ö. Legeza,
U. Schollwöck
Abstract:
We consider sudden quenches across quantum phase transitions in the $S=1$ XXZ model starting from the Haldane phase. We demonstrate that dynamical phase transitions may occur during these quenches that are identified by nonanalyticities in the rate function for the return probability. In addition, we show that the temporal behavior of the string order parameter is intimately related to the subsequ…
▽ More
We consider sudden quenches across quantum phase transitions in the $S=1$ XXZ model starting from the Haldane phase. We demonstrate that dynamical phase transitions may occur during these quenches that are identified by nonanalyticities in the rate function for the return probability. In addition, we show that the temporal behavior of the string order parameter is intimately related to the subsequent dynamical phase transitions. We furthermore find that the dynamical quantum phase transitions can be accompanied by enhanced two-site entanglement.
△ Less
Submitted 1 July, 2019; v1 submitted 1 April, 2019;
originally announced April 2019.
-
Interaction quench and thermalization in a one-dimensional topological Kondo insulator
Authors:
I. Hagymási,
C. Hubig,
U. Schollwöck
Abstract:
We study the nonequilibrium dynamics of a one-dimensional topological Kondo insulator, modelled by a $p$-wave Anderson lattice model, following a quantum quench of the on-site interaction strength. Our goal is to examine how the quench influences the topological properties of the system, therefore our main focus is the time evolution of the string order parameter, entanglement spectrum and the top…
▽ More
We study the nonequilibrium dynamics of a one-dimensional topological Kondo insulator, modelled by a $p$-wave Anderson lattice model, following a quantum quench of the on-site interaction strength. Our goal is to examine how the quench influences the topological properties of the system, therefore our main focus is the time evolution of the string order parameter, entanglement spectrum and the topologically-protected edge states. We point out that postquench local observables can be well captured by a thermal ensemble up to a certain interaction strength. Our results demonstrate that the topological properties after the interaction quench are preserved. Though the absolute value of the string order parameter decays in time, the analysis of the entanglement spectrum, Loschmidt echo and the edge states indicates the robustness of the topological properties in the time-evolved state. These predictions could be directly tested in state-of-the-art cold-atom experiments.
△ Less
Submitted 26 February, 2019; v1 submitted 23 October, 2018;
originally announced October 2018.
-
Entanglement and magnetism in high-spin graphene nanodisks
Authors:
I. Hagymasi,
O. Legeza
Abstract:
We investigate the ground-state properties of triangular graphene nanoflakes with zigzag edge configurations. The description of zero-dimensional nanostructures requires accurate many-body techniques since the widely used density-functional theory with local density approximation or Hartree-Fock methods cannot handle the strong quantum fluctuations. Applying the unbiased density-matrix renormaliza…
▽ More
We investigate the ground-state properties of triangular graphene nanoflakes with zigzag edge configurations. The description of zero-dimensional nanostructures requires accurate many-body techniques since the widely used density-functional theory with local density approximation or Hartree-Fock methods cannot handle the strong quantum fluctuations. Applying the unbiased density-matrix renormalization group algorithm we calculate the magnetization and entanglement patterns with high accuracy for different interaction strengths and compare them to the mean-field results. With the help of quantum information analysis and subsystem density matrices we reveal that the edges are strongly entangled with each other. We also address the effect of electron and hole doping and demonstrate that the magnetic properties of triangular nanoflakes can be controlled by electric field, which reveals features of flat-band ferromagnetism. This may open up new avenues in graphene based spintronics.
△ Less
Submitted 14 November, 2017;
originally announced November 2017.
-
A magnetic phase-transition graphene transistor with tunable spin polarization
Authors:
Peter Vancso,
Imre Hagymasi,
Levente Tapaszto
Abstract:
Graphene nanoribbons (GNRs) have been proposed as potential building blocks for field effect transistor (FET) devices due to their quantum confinement bandgap. Here, we propose a novel GNR device concept, enabling the control of both charge and spin signals, integrated within the simplest three-terminal device configuration. In a conventional FET device, a gate electrode is employed to tune the Fe…
▽ More
Graphene nanoribbons (GNRs) have been proposed as potential building blocks for field effect transistor (FET) devices due to their quantum confinement bandgap. Here, we propose a novel GNR device concept, enabling the control of both charge and spin signals, integrated within the simplest three-terminal device configuration. In a conventional FET device, a gate electrode is employed to tune the Fermi level of the system in and out of a static bandgap. By contrast, in the switching mechanism proposed here, the applied gate voltage can dynamically open and close an interaction gap, with only a minor shift of the Fermi level. Furthermore, the strong interplay of the band structure and edge spin configuration in zigzag ribbons enables such transistors to carry spin polarized current without employing an external magnetic field or ferromagnetic contacts. Using an experimentally validated theoretical model, we show that such transistors can switch at low voltages and high speed, and the spin polarization of the current can be tuned from 0% to 50% by using the same back gate electrode. Furthermore, such devices are expected to be robust against edge irregularities and can operate at room temperature. Controlling both charge and spin signal within the simplest FET device configuration could open up new routes in data processing with graphene based devices.
△ Less
Submitted 8 February, 2017;
originally announced February 2017.
-
Interaction effects in a chaotic graphene quantum billiard
Authors:
Imre Hagymasi,
Peter Vancso,
Andras Palinkas,
Zoltan Osvath
Abstract:
We investigate the local electronic structure of a Sinai-like, quadrilateral graphene quantum billiard with zigzag and armchair edges using scanning tunneling microscopy at room temperature. It is revealed that besides the $(\sqrt{3}\times\sqrt{3})R30$° superstructure, which is caused by the intervalley scattering, its overtones also appear in the STM measurements, which are attributed to the Umkl…
▽ More
We investigate the local electronic structure of a Sinai-like, quadrilateral graphene quantum billiard with zigzag and armchair edges using scanning tunneling microscopy at room temperature. It is revealed that besides the $(\sqrt{3}\times\sqrt{3})R30$° superstructure, which is caused by the intervalley scattering, its overtones also appear in the STM measurements, which are attributed to the Umklapp processes. We point out that these results can be well understood by taking into account the Coulomb interaction in the quantum billiard, accounting for both the measured density of state values and the experimentally observed topography patterns. The analysis of the level-spacing distribution substantiates the experimental findings as well. We also reveal the magnetic properties of our system which should be relevant in future graphene based electronic and spintronic applications.
△ Less
Submitted 27 January, 2017;
originally announced January 2017.
-
Entanglement, excitations and correlation effects in narrow zigzag graphene nanoribbons
Authors:
I. Hagymasi,
O. Legeza
Abstract:
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, \emph{et al.} Nature \textbf{531}, 489 (2016)], and their descriptions require more sophisticated methods since in this regime conventional methods, like mean-field or density-…
▽ More
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, \emph{et al.} Nature \textbf{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.
△ Less
Submitted 4 October, 2016; v1 submitted 10 May, 2016;
originally announced May 2016.
-
Characterization of a correlated topological Kondo insulator in one dimension
Authors:
I. Hagymasi,
O. Legeza
Abstract:
We investigate the ground-state of a p-wave Kondo-Heisenberg model introduced by Alexandrov and Coleman with an Ising-type anisotropy in the Kondo interaction and correlated conduction electrons. Our aim is to understand how they affect the stability of the Haldane state obtained in the SU(2) symmetric case without the Hubbard interaction. By applying the density-matrix renormalization group algor…
▽ More
We investigate the ground-state of a p-wave Kondo-Heisenberg model introduced by Alexandrov and Coleman with an Ising-type anisotropy in the Kondo interaction and correlated conduction electrons. Our aim is to understand how they affect the stability of the Haldane state obtained in the SU(2) symmetric case without the Hubbard interaction. By applying the density-matrix renormalization group algorithm and calculating the entanglement entropy we show that in the anisotropic case a phase transition occurs and a Néel state emerges above a critical value of the Coulomb interaction. These findings are also corroborated by the examination of the entanglement spectrum and the spin profile of the system which clarify the structure of each phase.
△ Less
Submitted 6 April, 2016; v1 submitted 18 January, 2016;
originally announced January 2016.
-
Momentum distribution functions in a one-dimensional extended periodic Anderson model
Authors:
I. Hagymasi,
J. Solyom,
O. Legeza
Abstract:
We study the momentum distribution of the electrons in an extended periodic Anderson model, where the interaction, $U_{cf}$, between itinerant and localized electrons is taken into account. In the symmetric half-filled model, due to the increase of the interorbital interaction, the $f$ electrons become more and more delocalized, while the itinerancy of conduction electrons decreases. Above a certa…
▽ More
We study the momentum distribution of the electrons in an extended periodic Anderson model, where the interaction, $U_{cf}$, between itinerant and localized electrons is taken into account. In the symmetric half-filled model, due to the increase of the interorbital interaction, the $f$ electrons become more and more delocalized, while the itinerancy of conduction electrons decreases. Above a certain value of $U_{cf}$ the $f$ electrons become again localized together with the conduction electrons. In the less than half-filled case, we observe that $U_{cf}$ causes strong correlations between the $f$ electrons in the mixed valence regime.
△ Less
Submitted 8 May, 2015;
originally announced May 2015.
-
Competition between Hund's coupling and Kondo effect in a one-dimensional extended periodic Anderson model
Authors:
I. Hagymasi,
J. Solyom,
O. Legeza
Abstract:
We study the ground-state properties of an extended periodic Anderson model to understand the role of Hund's coupling between localized and itinerant electrons using the density-matrix renormalization group algorithm. By calculating the von Neumann entropies we show that two phase transitions occur and two new phases appear as the hybridization is increased in the symmetric half-filled case due to…
▽ More
We study the ground-state properties of an extended periodic Anderson model to understand the role of Hund's coupling between localized and itinerant electrons using the density-matrix renormalization group algorithm. By calculating the von Neumann entropies we show that two phase transitions occur and two new phases appear as the hybridization is increased in the symmetric half-filled case due to the competition between Kondo-effect and Hund's coupling. In the intermediate phase, which is bounded by two critical points, we found a dimerized ground state, while in the other spatially homogeneous phases the ground state is Haldane-like and Kondo-singlet-like, respectively. We also determine the entanglement spectrum and the entanglement diagram of the system by calculating the mutual information thereby clarifying the structure of each phase.
△ Less
Submitted 19 June, 2015; v1 submitted 19 March, 2015;
originally announced March 2015.
-
Second-order Peierls transition in the spin-orbital Kumar-Heisenberg model
Authors:
Wojciech Brzezicki,
Imre Hagymási,
Jacek Dziarmaga,
Örs Legeza
Abstract:
We add a Heisenberg interaction term $\proptoλ$ in the one-dimensional SU(2)$\otimes$XY spin-orbital model introduced by B. Kumar. At $λ=0$ the spin and orbital degrees of freedom can be separated by a unitary transformation leading to an exact solution of the model. We show that a finite $λ>0$ leads to spontaneous dimerization of the system which in the thermodynamic limit becomes a smooth phase…
▽ More
We add a Heisenberg interaction term $\proptoλ$ in the one-dimensional SU(2)$\otimes$XY spin-orbital model introduced by B. Kumar. At $λ=0$ the spin and orbital degrees of freedom can be separated by a unitary transformation leading to an exact solution of the model. We show that a finite $λ>0$ leads to spontaneous dimerization of the system which in the thermodynamic limit becomes a smooth phase transition at $λ\to 0$, whereas it remains discontinuous within the first order perturbation approach. We present the behavior of the entanglement entropy, energy gap and dimerization order parameter in the limit of $λ\to 0$ confirming the critical behavior. Finally, we show the evidence of another phase transition in the Heisenberg limit, $λ\to\infty$, and give a qualitative analytical explanation of the observed dimerized states both in the limit of small and large $λ$.
△ Less
Submitted 13 March, 2015;
originally announced March 2015.
-
Room temperature magnetic order on zigzag edges of narrow graphene nanoribbons
Authors:
Gabor Zsolt Magda,
Xiaozhan Jin,
Imre Hagymasi,
Peter Vancso,
Zoltan Osvath,
Peter Nemes-Incze,
Chanyong Hwang,
Laszlo P. Biro,
Levente Tapaszto
Abstract:
Magnetic order emerging in otherwise non-magnetic materials as carbon is a paradigmatic example of a novel type of s-p electron magnetism predicted to be of exceptional high-temperature stability. It has been demonstrated that atomic scale structural defects of graphene can host unpaired spins. However, it is still unclear under which conditions long-range magnetic order can emerge from such defec…
▽ More
Magnetic order emerging in otherwise non-magnetic materials as carbon is a paradigmatic example of a novel type of s-p electron magnetism predicted to be of exceptional high-temperature stability. It has been demonstrated that atomic scale structural defects of graphene can host unpaired spins. However, it is still unclear under which conditions long-range magnetic order can emerge from such defect-bound magnetic moments. Here we propose that in contrast to random defect distributions, atomic scale engineering of graphene edges with specific crystallographic orientation, comprising edge atoms only from one sub-lattice of the bipartite graphene lattice, can give rise to a robust magnetic order. We employ a nanofabrication technique based on Scanning Tunneling Microscopy to define graphene nanoribbons with nanometer precision and well-defined crystallographic edge orientations. While armchair ribbons display quantum confinement gap, zigzag ribbons narrower than 7 nm reveal a bandgap of about 0.2 - 0.3 eV, which can be identified as a signature of interaction induced spin ordering along their edges. Moreover, a semiconductor to metal transition is revealed upon increasing the ribbon width, indicating the switching of the magnetic coupling between opposite ribbon edges from antiferromagnetic to ferromagnetic configuration. We found that the magnetic order on graphene edges of controlled zigzag orientation can be stable even at room temperature, raising hope for graphene-based spintronic devices operating under ambient conditions.
△ Less
Submitted 5 November, 2014;
originally announced November 2014.
-
Interorbital interaction in the one-dimensional periodic Anderson model: A density-matrix renormalization-group study
Authors:
I. Hagymasi,
J. Solyom,
O. Legeza
Abstract:
We investigate the effect of the Coulomb interaction, $U_{cf}$, between the conduction and f electrons in the periodic Anderson model using the density-matrix renormalization-group algorithm. We calculate the excitation spectrum of the half-filled symmetric model with an emphasis on the spin and charge excitations. In the one-dimensional version of the model it is found that the spin gap is smalle…
▽ More
We investigate the effect of the Coulomb interaction, $U_{cf}$, between the conduction and f electrons in the periodic Anderson model using the density-matrix renormalization-group algorithm. We calculate the excitation spectrum of the half-filled symmetric model with an emphasis on the spin and charge excitations. In the one-dimensional version of the model it is found that the spin gap is smaller than the charge gap below a certain value of $U_{cf}$ and the reversed inequality is valid for stronger $U_{cf}$. This behavior is also verified by the behavior of the spin and density correlation functions. We also perform a quantum information analysis of the model and determine the entanglement map of the f and conduction electrons. It is revealed that for a certain $U_{cf}$ the ground state is dominated by the configuration in which the conduction and f electrons are strongly entangled, and the ground state is almost a product state. For larger $U_{cf}$ the sites are occupied alternatingly dominantly by two f electrons or by two conduction electrons.
△ Less
Submitted 24 September, 2014; v1 submitted 24 June, 2014;
originally announced June 2014.
-
Quantum criticality and first-order transitions in the extended periodic Anderson model
Authors:
I. Hagymasi,
K. Itai,
J. Solyom
Abstract:
We investigate the behavior of the periodic Anderson model in the presence of $d$-$f$ Coulomb interaction ($U_{df}$) using mean-field theory, variational calculation, and exact diagonalization of finite chains. The variational approach based on the Gutzwiller trial wave function gives a critical value of $U_{df}$ and two quantum critical points (QCPs), where the valence susceptibility diverges. We…
▽ More
We investigate the behavior of the periodic Anderson model in the presence of $d$-$f$ Coulomb interaction ($U_{df}$) using mean-field theory, variational calculation, and exact diagonalization of finite chains. The variational approach based on the Gutzwiller trial wave function gives a critical value of $U_{df}$ and two quantum critical points (QCPs), where the valence susceptibility diverges. We derive the critical exponent for the valence susceptibility and investigate how the position of the QCP depends on the other parameters of the Hamiltonian. For larger values of $U_{df}$, the Kondo regime is bounded by two first-order transitions. These first-order transitions merge into a triple point at a certain value of $U_{df}$. For even larger $U_{df}$ valence skipping occurs. Although the other methods do not give a critical point, they support this scenario.
△ Less
Submitted 29 March, 2013; v1 submitted 19 December, 2012;
originally announced December 2012.
-
Hubbard physics in the symmetric half-filled periodic Anderson-Hubbard model
Authors:
I. Hagymasi,
K. Itai,
J. Solyom
Abstract:
Two very different methods -- exact diagonalization on finite chains and a variational method -- are used to study the possibility of a metal-insulator transition in the symmetric half-filled periodic Anderson-Hubbard model. With this aim we calculate the density of doubly occupied $d$ sites as a function of various parameters. In the absence of on-site Coulomb interaction ($U_f$) between $f$ elec…
▽ More
Two very different methods -- exact diagonalization on finite chains and a variational method -- are used to study the possibility of a metal-insulator transition in the symmetric half-filled periodic Anderson-Hubbard model. With this aim we calculate the density of doubly occupied $d$ sites as a function of various parameters. In the absence of on-site Coulomb interaction ($U_f$) between $f$ electrons, the two methods yield similar results. The double occupancy of $d$ levels remains always finite just as in the one-dimensional Hubbard model. Exact diagonalization on finite chains gives the same result for finite $U_f$, while the Gutzwiller method leads to a Brinkman-Rice transition at a critical value ($U_d^c$), which depends on $U_f$ and $V$.
△ Less
Submitted 2 July, 2012;
originally announced July 2012.
-
Periodic Anderson model with correlated conduction electrons: variational and exact diagonalization study
Authors:
I. Hagymasi,
K. Itai,
J. Solyom
Abstract:
We investigate an extended version of the periodic Anderson model (the so-called periodic Anderson-Hubbard model) with the aim to understand the role of interaction between conduction electrons in the formation of the heavy-fermion and mixed-valence states. Two methods are used: (i) variational calculation with the Gutzwiller wave function optimizing numerically the ground-state energy and (ii) ex…
▽ More
We investigate an extended version of the periodic Anderson model (the so-called periodic Anderson-Hubbard model) with the aim to understand the role of interaction between conduction electrons in the formation of the heavy-fermion and mixed-valence states. Two methods are used: (i) variational calculation with the Gutzwiller wave function optimizing numerically the ground-state energy and (ii) exact diagonalization of the Hamiltonian for short chains. The f-level occupancy and the renormalization factor of the quasiparticles are calculated as a function of the energy of the f-orbital for a wide range of the interaction parameters. The results obtained by the two methods are in reasonably good agreement for the periodic Anderson model. The agreement is maintained even when the interaction between band electrons, U_d, is taken into account, except for the half-filled case. This discrepancy can be explained by the difference between the physics of the one- and higher dimensional models. We find that this interaction shifts and widens the energy range of the bare f-level, where heavy-fermion behavior can be observed. For large enough U_d this range may lie even above the bare conduction band. The Gutzwiller method indicates a robust transition from Kondo insulator to Mott insulator in the half-filled model, while U_d enhances the quasi-particle mass when the filling is close to half filling.
△ Less
Submitted 31 May, 2012; v1 submitted 24 June, 2011;
originally announced June 2011.
-
Periodic Anderson model with d-f interaction
Authors:
I. Hagymasi,
K. Itai,
J. Solyom
Abstract:
We investigate an extended version of the periodic Anderson model where an interaction is switched on between the doubly occupied d- and f-sites. We perform variational calculations using the Gutzwiller trial wave function. We calculate the f-level occupancy as a function of the f-level energy with different interaction strengths. It is shown that the region of valence transition is sharpened due…
▽ More
We investigate an extended version of the periodic Anderson model where an interaction is switched on between the doubly occupied d- and f-sites. We perform variational calculations using the Gutzwiller trial wave function. We calculate the f-level occupancy as a function of the f-level energy with different interaction strengths. It is shown that the region of valence transition is sharpened due to the new interaction.
△ Less
Submitted 30 May, 2012; v1 submitted 24 June, 2011;
originally announced June 2011.
-
Critical Josephson current in ballistic superconductor-graphene systems
Authors:
I. Hagymási,
A. Kormányos,
J. Cserti
Abstract:
We calculate the phase, the temperature and the junction length dependence of the supercurrent for ballistic graphene Josephson-junctions. For low temperatures we find non-sinusoidal dependence of the supercurrent on the superconductor phase difference for both short and long junctions. The skewness, which characterizes the deviation of the current-phase relation from a simple sinusoidal one, show…
▽ More
We calculate the phase, the temperature and the junction length dependence of the supercurrent for ballistic graphene Josephson-junctions. For low temperatures we find non-sinusoidal dependence of the supercurrent on the superconductor phase difference for both short and long junctions. The skewness, which characterizes the deviation of the current-phase relation from a simple sinusoidal one, shows a linear dependence on the critical current for small currents. We discuss the similarities and differences with respect to the classical theory of Josephson junctions, where the weak link is formed by a diffusive or ballistic metal. The relation to other recent theoretical results on graphene Josephson junctions is pointed out and the possible experimental relevance of our work is considered as well.
△ Less
Submitted 16 September, 2010; v1 submitted 16 June, 2010;
originally announced June 2010.
-
Graphene Andreev Billiards
Authors:
J. Cserti,
I. Hagymasi,
A. Kormanyos
Abstract:
We studied the energy levels of graphene based Andreev billiards consisting of a superconductor region on top of a monolayer graphene sheet. For the case of Andreev retro-reflection we show that the graphene based Andreev billiard can be mapped to the normal metal-superconducting billiards with the same geometry. We also derived a semiclassical quantization rule in graphene based Andreev billiar…
▽ More
We studied the energy levels of graphene based Andreev billiards consisting of a superconductor region on top of a monolayer graphene sheet. For the case of Andreev retro-reflection we show that the graphene based Andreev billiard can be mapped to the normal metal-superconducting billiards with the same geometry. We also derived a semiclassical quantization rule in graphene based Andreev billiards. The exact and the semiclassically obtained spectrum agree very well both for the case of Andreev retro-reflection and specular Andreev reflection.
△ Less
Submitted 23 April, 2009; v1 submitted 17 February, 2009;
originally announced February 2009.