Source author record

Mikhail I. Katsnelson

Mikhail I. Katsnelson 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

81works
18topics
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

81 published item(s)

preprint2024arXiv

Strong electron-phonon coupling and phonon-induced superconductivity in tetragonal C$_3$N$_4$ with hole doping

C$_3$N$_4$ is a recently discovered phase of carbon nitrides with the tetragonal crystal structure [D.Laniel $\textit{et al.}$, Adv. Mater. 2023, 2308030] that is stable at ambient conditions. C$_3$N$_4$ is a semiconductor exhibiting flat-band anomalies in the valence band, suggesting the emergence of many-body instabilities upon hole doping. Here, using state-of-the-art first-principles calculations we show that hole-doped C$_3$N$_4$ reveals strong electron-phonon coupling, leading to the formation of a gapped superconducting state. The phase transition temperatures turn out to be strongly dependent on the hole concentration. We propose that holes could be injected into C$_3$N$_4$ via boron doping which induces, according to our results, a rigid shift of the Fermi energy without significant modification of the electronic structure. Based on the electron-phonon coupling and Coulomb pseudopotential calculated from first principles, we conclude that the boron concentration of 6 atoms per nm$^3$ would be required to reach the critical temperature of $\sim$36 K at ambient pressure.

preprint2023arXiv

Dielectric tunability of magnetic properties in orthorhombic ferromagnetic monolayer CrSBr

Monolayer CrSBr is a recently discovered semiconducting spin-3/2 ferromagnet with a Curie temperature around 146 K. Unlike many other known two-dimensional (2D) magnets, CrSBr has an orthorhombic lattice, giving rise, for instance, to spatial anisotropy of the magnetic excitations within the 2D plane. Theoretical description of CrSBr within a spin Hamiltonian approach turns out to be nontrivial due to the triaxial magnetic anisotropy as well as due to magnetic dipolar interactions, comparable to spin-orbit effects in CrSBr. Here, we employ a Green's function formalism combined with first-principles calculations to systematically study the magnetic properties of monolayer CrSBr in different regimes of surrounding dielectric screening. We find that the magnetic anisotropy and thermodynamical properties of CrSBr depend significantly on the Coulomb interaction and its external screening. In the free-standing limit, the system turns out to be close to an easy-plane magnet, whose long-range ordering is partially suppressed. On the contrary, in the regime of large external screening, monolayer CrSBr behaves like an easy-axis ferromagnet with more stable magnetic ordering. Despite being relatively large, the magnetic dipolar interactions have only little effect on the magnetic properties. Our findings suggests that 2D CrSBr is suitable platform for studying the effects of substrate screening on magnetic ordering in low dimensions.

preprint2022arXiv

Dynamical correlations in single-layer CrI$_3$

Chromium triiodide is an intrinsically magnetic van der Waals material down to the single-layer limit. Here, we provide a first-principles description of finite-temperature magnetic and spectral properties of monolayer (ML) CrI$_3$ based on fully charge self-consistent density functional theory (DFT) combined with dynamical mean-field theory, revealing a formation of local moments on Cr from strong local Coulomb interactions. We show that the presence of local dynamical correlations leads to a modification of the electronic structure of ferromagnetically ordered CrI$_3$. In contrast to conventional DFT+$U$ calculations, we find that the top of the valence band in ML CrI$_3$ demonstrates essentially different orbital character for minority and majority spin states, which is closer to the standard DFT results. This leads to a strong spin polarization of the optical conductivity upon hole doping, which could be verified experimentally.

preprint2022arXiv

Electronic and optical properties of crystalline nitrogen versus black phosphorus: A comparative first-principles study

Crystalline black nitrogen (BN) is an allotrope of nitrogen with the black phosphorus (BP) structure recently synthesized at high pressure by two independent research groups [Ji et al., Sci. Adv. 6, eaba9206 (2020); Laniel et al., Phys. Rev. Lett. 124, 216001 (2020)]. Here, we present a systematic study of the electronic and optical properties of BN focusing on its comparison with BP. To this end, we use the state-of-the-art quasiparticle self-consistent $GW$ approach with vertex corrections in both the electronic and optical channels. Despite many similarities, the properties of BN are found to be considerably different. Unlike BP, BN exhibits a larger optical gap (2.5 vs 0.26 eV), making BN transparent in the visible spectral region with a highly anisotropic optical response. This difference can be primarily attributed to a considerably reduced dielectric screening in BN, leading to enhancement of the effective Coulomb interaction. Despite relatively strong Coulomb interaction, exciton formation is largely suppressed in both materials. Our analysis of the elastic properties shows exceptionally high stiffness of BN, comparable to that of diamond.

preprint2022arXiv

Gating orbital memory with an atomic donor

Orbital memory is defined by two stable valencies that can be electrically switched and read-out. To explore the influence of an electric field on orbital memory, we studied the distance-dependent influence of an atomic Cu donor on the state favorability of an individual Co atom on black phosphorus. Using low temperature scanning tunneling microscopy/spectroscopy, we characterized the electronic properties of individual Cu donors, corroborating this behavior with ab initio calculations based on density functional theory. We studied the influence of an individual donor on the charging energy and stochastic behavior of an individual Co atom. We found a strong impact on the state favorability in the stochastic limit. These findings provide quantitative information about the influence of local electric fields on atomic orbital memory.

preprint2022arXiv

Orbital memory from individual Fe atoms on black phosphorus

Bistable valency in individual atoms presents a new approach toward single-atom memory, as well as a building block to create tunable and stochastic multi-well energy landscapes. Yet, this concept of orbital memory has thus far only been observed for cobalt atoms on the surface of black phosphorus, which are switched using tip-induced ionization. Here, we show that individual iron atoms on the surface of black phosphorus exhibit orbital memory using a combination of scanning tunneling microscopy and spectroscopy with ab initio calculations based on density functional theory. Unlike cobalt, the iron orbital memory can be switched in its non-ionized ground state. Based on calculations, we confirm that each iron valency has a distinct magnetic moment that is characterized by a distinguishable charge distribution due to the different orbital population. By studying the stochastic switching of the valency with varying tunneling conditions, we propose that the switching mechanism is based on a two-electron tunneling process.

preprint2022arXiv

Perturbative renormalization and thermodynamics of quantum crystalline membranes

We analyze the statistical mechanics of a free-standing quantum crystalline membrane within the framework of a systematic perturbative renormalization group (RG). A power-counting analysis shows that the leading singularities of correlation functions can be analyzed within an effective renormalizable model in which the kinetic energy of in-plane phonons and subleading geometrical nonlinearities in the expansion of the strain tensor are neglected. For membranes at zero temperature, governed by zero-point motion, the RG equations of the effective model provide a systematic derivation of logarithmic corrections to the bending rigidity and to the elastic Young modulus derived in earlier investigations. In the limit of a weakly applied external tension, the stress-strain relation at $T = 0$ is anomalous: the linear Hooke's law is replaced with a singular law exhibiting logarithmic corrections. For small, but finite temperatures, we use techniques of finite-size scaling to derive general relations between the zero-temperature RG flow and scaling laws of thermodynamic quantities such as the thermal expansion coefficient $α$, the entropy S, and the specific heat C. A combination of the scaling relations with an analysis of thermal fluctuations shows that, for small temperatures, the thermal expansion coefficient $α$ is negative and logarithmically dependent on $T$, as predicted in an earlier work. Although the requirement $\lim_{T \to 0} α= 0$, expected from the third law of thermodynamics is formally satisfied, $α$ is predicted to exhibit such a slow variation to remain practically constant down to unaccessibly small temperatures.

preprint2021arXiv

Distribution of ripples in graphene membrane

Intrinsic ripples with various configurations and sizes were reported to affect the physical and chemical properties of 2D materials. By performing molecular dynamics simulations and theoretical analysis, we use two geometric models of the ripple shape to explore numerically the distribution of ripples in graphene membrane. We focus on the ratio of ripple height to its diameter (t/D) which was recently shown to be the most relevant for chemical activity of graphene membranes. Our result demonstrates that the ripple density decreases as the coefficient t/D increases, in a qualitative agreement with the Boltzmann distribution derived analytically from the bending energy of the membrane. Our theoretical study provides also specific quantitative information on the ripple distribution in graphene and gives new insights applicable to other 2D materials.

preprint2021arXiv

Electronic and Optical properties of transition metal dichalcogenides under symmetric and asymmetric field-effect doping

Doping via electrostatic gating is a powerful and widely used technique to tune the electron densities in layered materials. The microscopic details of how these setups affect the layered material are, however, subtle and call for careful theoretical treatments. Using semiconducting monolayers of transition metal dichalcogenides (TMDs) as prototypical systems affected by electrostatic gating, we show that the electronic and optical properties change indeed dramatically when the gating geometry is properly taken into account. This effect is implemented by a self-consistent calculation of the Coulomb interaction between the charges in different sub-layers within the tight-binding approximation. Thereby we consider both, single- and double-sided gating. Our results show that, at low doping levels of $10^{13}$ cm$^{-2}$, the electronic bands of monolayer TMDs shift rigidly for both types of gating, and subsequently undergo a Lifshitz transition. When approaching the doping level of $10^{14}$ cm$^{-2}$, the band structure changes dramatically, especially in the case of single-sided gating where we find that monolayer \ce{MoS2} and \ce{WS2} become indirect gap semiconductors. The optical conductivities calculated within linear response theory also show clear signatures of these doping-induced band structure renormalizations. Our numerical results based on light-weighted tight-binding models indicate the importance of electronic screening in doped layered structures, and pave the way for further understanding gated super-lattice structures formed by mutlilayers with extended Moiré pattern.

preprint2021arXiv

Excitons in Bulk and Layered Chromium Tri-Halides: From Frenkel to the Wannier-Mott Limit

Excitons with large binding energies $\sim$2-3 eV in CrX$_{3}$ are historically characterized as being localized (Frenkel) excitons that emerge from the atomic $d{-}d$ transitions between the Cr-3$d$-$t_{2g}$ and $e_{g}$ orbitals. The argument has gathered strength in recent years as the excitons in recently made monolayers are found at almost the same energies as the bulk. The Laporte rule, which restricts such parity forbidden atomic transitions, can relax if, at least, one element is present: spin-orbit coupling, odd-parity phonons or Jahn-Teller distortion. While what can be classified as a purely Frenkel exciton is a matter of definition, we show using an advanced first principles parameter-free approach that these excitons in CrX$_{3}$, in both its bulk and monolayer variants, have band-origin and do not require the relaxation of Laporte rule as a fundamental principle. We show that, the character of these excitons is mostly determined by the Cr-$d$ orbital manifold, nevertheless, they appear only as a consequence of X-p states hybridizing with the Cr-$d$. The hybridization enhances as the halogen atom becomes heavier, bringing the X-$p$ states closer to the Cr-$d$ states in the sequence Cl{\textrightarrow}Br{\textrightarrow}I, with an attendant increase in exciton intensity and decrease in binding energy. By applying a range of different kinds of perturbations, we show that, moderate changes to the two-particle Hamiltonian that essentially modifies the Cr-$d$-X-$p$ hybridization, can alter both the intensities and positions of the exciton peaks. A detailed analysis of several deep lying excitons, with and without strain, reveals that the exciton is most Frenkel like in CrCl$_{3}$ and acquires mixed Frenkel-Wannier character in CrI$_{3}$.

preprint2021arXiv

First principles vs second principles: Role of charge self-consistency in strongly correlated systems

First principles approaches have been successful in solving many-body Hamiltonians for real materials to an extent when correlations are weak or moderate. As the electronic correlations become stronger often embedding methods based on first principles approaches are used to better treat the correlations by solving a suitably chosen many-body Hamiltonian with a higher level theory. Such combined methods are often referred to as second principles approaches. At such level of the theory the self energy, i.e. the functional that embodies the stronger electronic correlations, is either a function of energy or momentum or both. The success of such theories is commonly measured by the quality of the self energy functional. However, self-consistency in the self-energy should, in principle, also change the real space charge distribution in a correlated material and be able to modify the electronic eigenfunctions, which is often undermined in second principles approaches. Here we study the impact of charge self-consistency within two example cases: TiSe$_{2}$, a three-dimensional charge-density-wave candidate material, and CrBr$_{3}$, a two-dimensional ferromagnet, and show how real space charge re-distribution due to correlation effects taken into account within a first principles Green's function based many-body perturbative approach is key in driving qualitative changes to the final electronic structure of these materials.

preprint2021arXiv

Interlayer hybridization in graphene quasicrystal and other bilayer graphene systems

The incommensurate 30$^{\circ}$ twisted bilayer graphene (BG) possesses both relativistic Dirac fermions and quasiperiodicity with 12-fold rotational symmetry arising from the interlayer interaction [\href{https://science.sciencemag.org/content/361/6404/782}{Ahn et al., Science \textbf{361}, 782 (2018)} and \href{https://www.pnas.org/content/115/27/6928}{Yao et al., Proc. Natl. Acad. Sci. \textbf{115}, 6928 (2018)}]. Understanding how the interlayer states interact with each other is of vital importance for identifying and subsequently engineering the quasicrystalline order in the layered structures. Herein, via symmetry and group representation theory we unravel the interlayer hybridization selection rules governing the interlayer coupling in both untwisted and twisted BG systems. Compared with the only allowed equivalent hybridization in $D_{6h}$ untwisted BG, $D_6$ twisted BG permits equivalent and mixed hybridizations, and $D_{6d}$ graphene quasicrystal allows both equivalent and non-equivalent hybridizations. The energy-dependent hybridization strengths in graphene quasicrystal and $D_6$ twisted BG show two remarkable characteristics: (i) near the Fermi level the weak hybridization owing to the relatively large energy difference between Dirac bands from top and bottom layers, and (ii) in high-energy regions the electron-hole asymmetry of hybridization strength with stronger interlayer coupling for holes, which arises from the non-nearest-neighbor interlayer hoppings and the wave-function phase difference between paring states. These hybridization-generated band structures and their hybridization strength characteristics are verified by the calculated optical conductivity spectra. Our theoretical study paves a way for revealing the interlayer hybridization in van der Waals layered systems.

preprint2021arXiv

Plasmonic Quantum Dots in Twisted Bilayer Graphene

We derive a material-realistic real-space many-body Hamiltonian for twisted bilayer graphene from first principles, including both single-particle hopping terms for $p_z$ electrons and long-range Coulomb interactions. By disentangling low- and high-energy subspaces of the electronic dispersion, we are able to utilize state-of-the-art constrained Random Phase Approximation calculations to reliably describe the non-local background screening from the high-energy $s$, $p_x$, and $p_y$ electron states for arbitrary twist angles. The twist-dependent low-energy screening from $p_z$ states is subsequently added to obtain a full screening model. We use this approach to study real-space plasmonic patterns in electron-doped twisted bilayer graphene supercells and find, next to classical dipole-like modes, also twist-angle-dependent plasmonic quantum-dot-like excitations with $s$ and $p$ symmetries. Based on their inter-layer charge modulations and their footprints in the electron energy loss spectrum, we can classify these modes into "bright" and "dark" states, which show different dependencies on the twist angle.

preprint2021arXiv

Thermally-induced magnetic order from glassiness in elemental neodymium

Temperature in thermodynamics is synonymous with disorder, and responsible for ultimately destroying ordered phases. Here, we show an unusual magnetic transition where, with increasing the temperature of elemental neodymium, long-range multi-Q magnetic order emerges from a self-induced spin glass. Using temperature-dependent spin-polarized scanning tunneling microscopy, we characterize the local Q order in the spin-Q glass phase and quantify the emergence of long-range multi-Q order with increasing temperature. We develop two distinct analysis tools, which enable the quantification of the glass transition temperature, based on measured spatially-dependent magnetization. We compare these observations with atomic spin dynamics simulations, which reproduce the qualitative observation of a phase transition from a low-temperature spin glass phase to an intermediate ordered multi-Q phase. These simulations trace the origin of the unexpected high temperature order in weakened frustration driven by temperature-dependent sublattice correlations. These findings constitute an example of order from disorder and provide a rich platform to study magnetization dynamics in a self-induced spin glass.

preprint2020arXiv

Control of magnetic interactions between surface adatoms via orbital repopulation

We propose a reversible mechanism for switching Heisenberg-type exchange interactions between deposited transition metal adatoms from ferromagnetic to antiferromagnetic. Using first-principles calculations, we show that this mechanism can be realized for cobalt atoms on the surface of black phosphorus by making use of electrically-controlled orbital repopulation, as recently demonstrated by scanning probe techniques [Nat. Commun. 9, 3904 (2018)]. We find that field-induced repopulation not only affects the spin state, but also causes considerable modification of exchange interaction between adatoms, including its sign. Our model analysis demonstrates that variable adatom-substrate hybridization is a key factor responsible for this modification. We perform quantum simulations of inelastic tunneling characteristics and discuss possible ways to verify the proposed mechanism experimentally.

preprint2020arXiv

Detecting quantum critical points in the $t-t'$ Fermi-Hubbard model via complex network theory

A considerable success in phenomenological description of high-T$_{\rm c}$ superconductors has been achieved within the paradigm of Quantum Critical Point (QCP) - a parental state of a variety of exotic phases that is characterized by dense entanglement and absence of well-defined quasiparticles. However, the microscopic origin of the critical regime in real materials remains an open question. On the other hand, there is a popular view that a single-band $t-t'$ Hubbard model is the minimal model to catch the main relevant physics of superconducting compounds. Here, we suggest that emergence of the QCP is tightly connected with entanglement in real space and identify its location on the phase diagram of the hole-doped $t-t'$ Hubbard model. To detect the QCP we study a weighted graph of inter-site quantum mutual information within a four-by-four plaquette that is solved by exact diagonalization. We demonstrate that some quantitative characteristics of such a graph, viewed as a complex network, exhibit peculiar behavior around a certain submanifold in the parametric space of the model. This method allows us to overcome difficulties caused by finite size effects and to identify the transition point even on a small lattice, where long-range asymptotics of correlation functions cannot be accessed.

preprint2020arXiv

Electron correlation effects on exchange interactions and spin excitations in 2D van der Waals materials

Despite serious effort, the nature of the magnetic interactions and the role of electron-correlation effects in magnetic two-dimensional (2D) van der Waals materials remain elusive. Using CrI$_3$ as a model system, we show that the calculated electronic structure including nonlocal electron correlations yields spin excitations consistent with inelastic neutron scattering measurements. Remarkably, this approach identifies an unreported correlation-enhanced interlayer super-superexchange, which rotates the magnon Dirac lines off, and introduces a gap along, the high-symmetry $Γ$-$K$-$M$ path. This discovery provides a different perspective on the gap opening mechanism observed in CrI$_3$, which was previously associated with spin-orbit coupling through the Dzyaloshinskii-Moriya interaction or Kitaev interaction. Our observation elucidates the critical role of electron correlations on the spin ordering and spin dynamics in magnetic van der Waals materials and demonstrates the necessity of explicit treatment of electron correlations in the broad family of 2D magnetic materials.

preprint2020arXiv

Electronic structure of 30° twisted double bilayer graphene

In this paper, the electronic properties of 30° twisted double bilayer graphene, which loses the translational symmetry due to the incommensurate twist angle, are studied by means of the tight-binding approximation. We demonstrate the interlayer decoupling in the low-energy region from various electronic properties, such as the density of states, effective band structure, optical conductivity and Landau level spectrum. However, at Q points, the interlayer coupling results in the appearance of new Van Hove singularities in the density of states, new peaks in the optical conductivity and importantly the 12-fold-symmetry-like electronic states. The k-space tight-binding method is adopted to explain this phenomenon. The electronic states at Q points show the charge distribution patterns more complex than the 30° twisted bilayer graphene due to the symmetry decrease. These phenomena appear also in the 30° twisted interface between graphene monolayer and AB stacked bilayer.

preprint2020arXiv

Linearized spectral decimation in fractals

In this article we study the energy level spectrum of fractals which have block-hierarchical structures. We develop a method to study the spectral properties in terms of linearization of spectral decimation procedure and verify it numerically. Our approach provides qualitative explanations for various spectral properties of self-similar graphs within the theory of dynamical systems, including power-law level-spacing distribution, smooth density of states and effective chaotic regime.

preprint2020arXiv

Nonequilibrium dual-boson approach

We develop nonequilibrium auxiliary quantum master equation dual boson method (aux-DB), and argue that it presents a convenient way to describe steady states of correlated impurity models (such as single molecule optoelectronic devices) where electron and energy transport should be taken into account. The aux-DB is shown to provide high accuracy with relatively low numerical cost. Theoretical analysis is followed by illustrative simulations within generic two-level junction model, where the new scheme is benchmarked against numerically exact results.

preprint2020arXiv

Phonon-mediated superconductivity in strongly correlated electron systems: a Luttinger-Ward functional approach

We use a Luttinger-Ward functional approach to study the problem of phonon-mediated superconductivity in electron systems with strong electron-electron interactions (EEIs). Our derivation does not rely on an expansion in skeleton diagrams for the EEI and the resulting theory is therefore nonperturbative in the strength of the latter. We show that one of the building blocks of the theory is the irreducible six-leg vertex related to EEIs. Diagrammatically, this implies five contributions (one of the Fock and four of the Hartree type) to the electronic self-energy, which, to the best of our knowledge, have never been discussed in the literature. Our approach is applicable to (and in fact designed to tackle superconductivity in) strongly correlated electron systems described by generic lattice models, as long as the glue for electron pairing is provided by phonons.

preprint2020arXiv

Pressure and electric field dependence of quasicrystalline electronic states in 30$^{\circ}$ twisted bilayer graphene

30$^{\circ}$ twisted bilayer graphene demonstrates the quasicrystalline electronic states with 12-fold symmetry. These states are however far away from the Fermi level, which makes conventional Dirac fermion behavior dominating the low energy spectrum in this system. By using tight-binding approximation, we study the effect of external pressure and electric field on the quasicrystalline electronic states. Our results show that by applying the pressure perpendicular to graphene plane one can push the quasicrystalline electronic states towards the Fermi level. Then, the electron or hole doping of the order of $\sim$ $4\times10^{14}$ $cm^{-2}$ is sufficient for the coincidence of the Fermi level with these quasicrystalline states. Moreover, our study indicates that applying the electric field perpendicular to the graphene plane can destroy the 12-fold symmetry of these states and break the energy degeneracy of the 12-wave states, and it is easier to reach this in the conduction band than in the valence band. Importantly, the application of the pressure can recover the 12-fold symmetry of these states to some extent against the electric field. We propose a hybridization picture which can explain all these phenomena.

preprint2020arXiv

Punctuated equilibrium as the default mode of evolution of large populations on fitness landscapes dominated by saddle points in the weak-mutation limit

Punctuated equilibrium is a mode of evolution in which phenetic change occurs in rapid bursts that are separated by much longer intervals of stasis during which mutations accumulate but no major phenotypic change occurs. Punctuated equilibrium has been originally proposed within the framework of paleobiology, to explain the lack of transitional forms that is typical of the fossil record. Theoretically, punctuated equilibrium has been linked to self-organized criticality (SOC), a model in which the size of avalanches in an evolving system is power-law distributed, resulting in increasing rarity of major events. We show here that, under the weak-mutation limit, a large population would spend most of the time in stasis in the vicinity of saddle points in the fitness landscape. The periods of stasis are punctuated by fast transitions, in lnNe time (Ne, effective population size), when a new beneficial mutation is fixed in the evolving population, which moves to a different saddle, or on much rarer occasions, from a saddle to a local peak. Thus, punctuated equilibrium is the default mode of evolution under a simple model that does not involve SOC or other special conditions.

preprint2020arXiv

Scaling behavior of crystalline membranes: an $ε$-expansion approach

We study the scaling behavior of two-dimensional (2D) crystalline membranes in the flat phase by a renormalization group (RG) method and an $ε$-expansion. Generalization of the problem to non-integer dimensions, necessary to control the $ε$-expansion, is achieved by dimensional continuation of a well-known effective theory describing out-of-plane fluctuations coupled to phonon-mediated interactions via a scalar composite field, equivalent for small deformations to the local Gaussian curvature. The effective theory, which will be referred to as Gaussian curvature interaction (GCI) model, is equivalent to theories of elastic $D$-dimensional manifolds fluctuating in a $(D + d_{c})$-dimensional embedding space in the physical case $D = 2$ for arbitrary $d_{c}$. For $D\neq 2$, instead, the GCI model is not equivalent to a direct dimensional continuation of elastic membrane theory and it defines an alternative generalization to generic internal dimension $D$. We calculate explicitly RG functions at two-loop order and determine the exponent $η$ characterizing the long-wavelength scaling of correlation functions to order $ε^{2}$ in an $ε=(4-D)$-expansion. The self-consistent screening approximation (SCSA) for the GCI model is shown to be exact to O($ε^{2}$). For $d_{c} = 1$, the O($ε^{2}$) correction is suppressed by a small numerical prefactor. As a result, despite the large value of $ε= 2$, extrapolation of the first and second order results to $D = 2$ leads to very close numbers, $η= 0.8$ and $η\simeq 0.795$. The calculated exponent values are close to earlier reference results obtained by non-perturbative RG, the SCSA and numerical simulations. These indications suggest that a perturbative analysis of the GCI model could provide an useful framework for accurate quantitative predictions of the scaling exponent even at $D = 2$.

preprint2020arXiv

Self-induced spin glass state in elemental and crystalline neodymium

Spin glasses are a highly complex magnetic state of matter, intricately linked to spin frustration and structural disorder. They exhibit no long-range order and exude aging phenomena, distinguishing them from quantum spin liquids. We report a new type of spin glass state, the spin-Q glass, observable in bulk-like crystalline metallic neodymium thick films. Using spin-polarized scanning tunneling microscopy combined with ab initio calculations and atomistic spin-dynamics simulations, we visualize the variations in atomic-scale non-collinear order and its response to magnetic field and temperature. We quantify the aging phenomena relating the glassiness to crystalline symmetry and the energy landscape. This result not only resolves the long-standing debate of the magnetism of neodymium, but also suggests that glassiness may arise in other magnetic solids lacking extrinsic disorder.

preprint2019arXiv

Anisotropic destruction of the Fermi surface in inhomogeneous holographic lattices

We analyze fermionic response of strongly correlated holographic matter in presence of inhomogeneous periodically modulated potential mimicking the crystal lattice. The modulation is sourced by a scalar operator that explicitly breaks the translational symmetry in one direction. We compute the fermion spectral function and show that it either exhibits a well defined Fermi surface with umklapp gaps opening on the Brillouin zone boundary at small lattice wave vector, or, when the wave vector is large, the Fermi surface is anisotropically deformed and the quasiparticles get significantly broadened in the direction of translation symmetry breaking. Making use of the ability of our model to smoothly extrapolate to the homogeneous Q-lattice like setup, we show that this novel effect is not due to the periodic modulation of the potential and Umklapp physics, but rather due to the anisotropic features of the holographic horizon. That means it encodes novel physics of strongly correlated critical systems which may be relevant for phenomenology of exotic states of electron matter.

preprint2019arXiv

Atom-by-atom construction of attractors in a tunable finite size spin array

We demonstrate that a two-dimensional finite and periodic array of Ising spins coupled via RKKY-like exchange can exhibit tunable magnetic states ranging from three distinct magnetic regimes: (1) a conventional ferromagnetic regime, (2) a glass-like regime, and (3) a new multi-well regime. These magnetic regimes can be tuned by one gate-like parameter, namely the ratio between the lattice constant and the oscillating interaction wavelength. We characterize the various magnetic regimes, quantifying the distribution of low energy states, aging relaxation dynamics, and scaling behavior. The glassy and multi-well behavior results from the competing character of the oscillating long-range exchange interactions. The multi-well structure features multiple attractors, each with a sizable basin of attraction. This may open the possible application of such atomic arrays as associative memories.

preprint2019arXiv

Exactly solvable model of strongly correlated $d$-wave superconductivity

We present an infinite-dimensional lattice of two-by-two plaquettes, the quadruple Bethe lattice, with Hubbard interaction and solve it exactly by means of the cluster dynamical mean-field theory. It exhibits a $d$-wave superconducting phase that is related to a highly degenerate point in the phase diagram of the isolated plaquette at that the groundstates of the particle number sectors $N=2,3,4$ cross. The superconducting gap is formed by the renormalized lower Slater peak of the correlated, hole-doped Mott insulator. We engineer parts of the interaction and find that pair hoppings between $X/Y$-momenta are the main two-particle correlations of the superconducting phase. The suppression of the superconductivity in the overdoped regime is caused by the diminishing of pair hopping correlations and in the underdoped regime by charge blocking. The optimal doping is $\sim 0.15$ at which the underlying normal state shows a Lifshitz transition. The model allows for different intra- and inter-plaquette hoppings that we use to disentangle superconductivity from antiferromagnetism as the latter requires larger inter-plaquette hoppings.

preprint2019arXiv

Hall conductivity of Sierpinski carpet

We calculate the Hall conductivity of a Sierpinski carpet using Kubo-Bastin formula. The quantization of Hall conductivity disappears when we increase the depth of the fractal. The Hall conductivity is no more proportional to the Chern number. Nevertheless, these quantities behave in a similar way showing some reminiscence of a topological nature of the Hall conductivity. We also study numerically the bulk-edge correspondence and find that the edge states become less manifested when the depth of Sierpinski carpet is increased.

preprint2019arXiv

Large-area, periodic, and tunable intrinsic pseudo-magnetic fields in low-angle twisted bilayer graphene

A properly strained graphene monolayer or bilayer is expected to harbour periodic pseudo-magnetic fields with high symmetry, yet to date, a convincing demonstration of such pseudo-magnetic fields has been lacking, especially for bilayer graphene. Here, we report the first definitive experimental proof for the existence of large-area, periodic pseudo-magnetic fields, as manifested by vortex lattices in commensurability with the moiré patterns of low-angle twisted bilayer graphene. The pseudo-magnetic fields are strong enough to confine the massive Dirac electrons into circularly localized pseudo-Landau levels, as observed by scanning tunneling microscopy/spectroscopy, and also corroborated by tight-binding calculations. We further demonstrate that the geometry, amplitude, and periodicity of the pseudo-magnetic field can be fine-tuned by both the rotation angle and heterostrain applied to the system. Collectively, the present study substantially enriches twisted bilayer graphene as a powerful enabling platform for exploration of new and exotic physical phenomena, including quantum valley Hall effects and quantum anomalous Hall effects.

preprint2016arXiv

Anomalous Magnetothermopower in a Metallic Frustrated Antiferromagnet

We report the temperature $T$ and magnetic field $H$ dependence of the thermopower $S$ of an itinerant triangular antiferromagnet PdCrO$_2$ in high magnetic fields up to 32 T. In the paramagnetic phase, the zero-field thermopower is positive with a value typical of good metals with a high carrier density. In marked contrast to typical metals, however, $S$ decreases rapidly with increasing magnetic field, approaching zero at the maximum field scale for $T >$ 70 K. We argue here that this profound change in the thermoelectric response derives from the strong interaction of the 4$d$ correlated electrons of the Pd ions with the short-range spin correlations of the Cr$^{3+}$ spins that persist beyond the Néel ordering temperature due to the combined effects of geometrical frustration and low dimensionality.

preprint2016arXiv

Double occupancy in DMFT and the Dual Boson approach

We discuss the calculation of the double occupancy using Dynamical Mean-Field Theory (DMFT) in finite dimensions. The double occupancy can be determined from the susceptibility of the auxiliary impurity model or from the lattice susceptibility. The former method underestimates, whereas the latter overestimates the double occupancy. We illustrate this for the square-lattice Hubbard model. We propose an approach for which both methods lead to identical results by construction and which resolves this ambiguity. This self-consistent dual boson scheme results in a double occupancy that is numerically close to benchmarks available in the literature.

preprint2016arXiv

Edge pseudo-magnetoplasmons

We study the properties of edge plasmons in two-component electron liquids in the presence of pseudomagnetic fields, which have opposite signs for the two different electronic populations and therefore preserve the time-reversal symmetry. The physical realizations of such systems are many. We discuss the cases of strained graphene and of electrons in proximity to a Skyrmion lattice, solving the problem with the Wiener-Hopf technique. We show (i) that two charged counter-propagating acoustic edge modes exist at the boundary and (ii) that, in the limit of large pseudomagnetic fields, each of them involves oscillations of only one of the two electronic components. We suggest that the edge pseudo-magnetoplasmons of graphene can be used to selectively address the electrons of one specific valley, a feature relevant for the emerging field of valleytronics. Conversely, the spin-polarized plasmons at the boundary of Skyrmion lattices can be exploited for spintronics applications. Our solution highlights new features missing in previous (similar) results obtained with uncontrolled approximations, namely a logarithmic divergence of the plasmon velocity, and the absence of gapped edge modes inside the bulk-plasmon gap.

preprint2016arXiv

Magnetism and exchange interaction of small rare-earth clusters; Tb as a representative

Here we follow, both experimentally and theoretically, the development of magnetism in Tb clusters from the atomic limit, adding one atom at a time. The exchange interaction is, surprisingly, observed to drastically increase compared to that of bulk, and to exhibit irregular oscillations as a function of the interatomic distance. From electronic structure theory we find that the theoretical magnetic moments oscillate with cluster size in exact agreement with experimental data. Unlike the bulk, the oscillation is not caused by the RKKY mechanism. Instead, the inter-atomic exchange is shown to be driven by a competition between wave-function overlap of the 5d shell and the on-site exchange interaction, which leads to a competition between ferromagnetic double-exchange and antiferromagnetic super-exchange. This understanding opens up new ways to tune the magnetic properties of rare-earth based magnets with nano-sized building blocks.

preprint2016arXiv

Non-equilibrium itinerant-electron magnetism: a time-dependent mean-field theory

We study the dynamical magnetic susceptibility of a strongly correlated electronic system in the presence of a time-dependent hopping field, deriving a generalized Bethe-Salpeter equation which is valid also out of equilibrium. Focusing on the single-orbital Hubbard model within the time-dependent Hartree-Fock approximation, we solve the equation in the non-equilibrium adiabatic regime, obtaining a closed expression for the transverse magnetic susceptibility. From this, we provide a rigorous definition of non-equilibrium (time-dependent) magnon frequencies and exchange parameters, expressed in terms of non-equilibrium single-electron Green functions and self-energies. In the particular case of equilibrium, we recover previously known results.

preprint2016arXiv

Quantum Hall effect and semiconductor to semimetal transition in biased black phosphorus

We study the quantum Hall effect of 2D electron gas in black phosphorus in the presence of perpendicular electric and magnetic fields. In the absence of a bias voltage, the external magnetic field leads to a quantization of the energy spectrum into equidistant Landau levels, with different cyclotron frequencies for the electron and hole bands. The applied voltage reduces the band gap, and eventually a semiconductor to semimetal transition takes place. This nontrivial phase is characterized by the emergence of a pair of Dirac points in the spectrum. As a consequence, the Landau levels are not equidistant anymore, but follow the $\varepsilon_n\propto \sqrt{nB}$ characteristic of Dirac crystals as graphene. By using the Kubo-Bastin formula in the context of the kernel polynomial method, we compute the Hall conductivity of the system. We obtain a $σ_{xy}\propto 2n$ quantization of the Hall conductivity in the gapped phase (standard quantum Hall effect regime), and a $σ_{xy}\propto 4(n+1/2)$ quantization in the semimetalic phase, characteristic of Dirac systems with non-trivial topology.

preprint2016arXiv

Self-induced glassiness and pattern formation in spin systems subject to long-range interactions

We study the glass formation in two- and three-dimensional Ising and Heisenberg spin systems subject to competing interactions and uniaxial anisotropy with a mean-field approach. In three dimensions, for sufficiently strong anisotropy the systems always modulates in a striped phase. Below a critical strength of the anisotropy, a glassy phase exists in a finite range of temperature, and it becomes more stable as the system becomes more isotropic. In two dimension the criticality is always avoided and the glassy phase always exists.

preprint2016arXiv

Versatile Approach to the Spin Dynamics in Correlated Electron Systems

Time-dependent spin phenomena in condensed matter are most often either described in the weakly correlated limit of metallic Stoner/Slater-like magnetism via band theory or in the strongly correlated limit of Heisenberg-like interacting spins in an insulator. However many experimental studies, e.g. of (de)magnetization processes, focus on itinerant local-moment materials such as transition metals and various of their compounds. We here present a general theoretical framework that is capable of addressing correlated spin dynamics, also in the presence of a vanishing charge gap. A real-space implementation of the time-dependent rotational-invariant slave boson methodology allows to treat non-equilibrium spins numerically fast and efficiently beyond linear response as well as beyond the band-theoretical or Heisenberg limit.

preprint2015arXiv

Dirac points with giant spin-orbit splitting in the electronic structure of two-dimensional transition-metal carbides

Two-dimensional (2D) materials, especially their most prominent member, graphene, have greatly influenced many scientific areas. Moreover, they have become a base for investigating the relativistic properties of condensed matter within the emerging field of Dirac physics. This has ignited an intense search for new materials where charge carriers behave as massless or massive Dirac fermions. Here, we theoretically show the existence of Dirac electrons in a series of 2D transition-metal carbides, known as MXenes. They possess twelve conical crossings in the 1st Brillouin zone with giant spin-orbit splitting. Our findings indicate that the 2D band structure of MXenes is protected against external perturbations and preserved even in multilayer phases. These results, together with the broad possibilities to engineer the properties of these materials phases, make Dirac MXenes a potential candidate for studying and developing novel Dirac-physics-based technologies.

preprint2015arXiv

Electronic Structures and Optical Properties of Partially and Fully Fluorinated Graphene

In this letter we study the electronic structures and optical properties of partially and fully fluorinated graphene by a combination of abinitio G0W0 calculations and large-scale multi-orbital tight-binding simulations. We find that for partially fluorinated graphene, the appearance of paired fluorine atoms is more favorable than unpaired atoms. We also show that different types of structural disorder, such as carbon vacancies, fluorine vacancies, fluorine vacancy-clusters and fluorine armchair- and zigzag-clusters, will introduce different types of midgap states and extra excitations within the optical gap. Furthermore we argue that the local formation of $sp^3$ bonds upon fluorination can be distinguished from other disorder inducing mechanisms which do not destroy the $sp^2$ hybrid orbitals by measuring the polarization rotation of passing polarized light.

preprint2015arXiv

Helicons in Weyl semimetals

Helicons are transverse electromagnetic waves propagating in three-dimensional (3D) electron systems subject to a static magnetic field. We present a theory of helicons propagating through a 3D Weyl semimetal. Our approach relies on the evaluation of the optical conductivity tensor from semiclassical Boltzmann transport theory, with the inclusion of certain Berry curvature corrections that have been neglected in the earlier literature (such as the one due to the orbital magnetic moment). We demonstrate that the axion term characterizing the electromagnetic response of Weyl semimetals dramatically alters the helicon dispersion with respect to that in nontopological metals. We also discuss axion-related anomalies that appear in the plasmon dispersion relation.

preprint2015arXiv

Magnetic interactions in strongly correlated systems: spin and orbital contributions

We present a technique to map an electronic model with local interactions (a generalized multi-orbital Hubbard model) onto an effective model of interacting classical spins, by requiring that the thermodynamic potentials associated to spin rotations in the two systems are equivalent up to second order in the rotation angles. This allows to determine the parameters of relativistic and non-relativistic magnetic interactions in the effective spin model in terms of equilibrium Green's functions of the electronic model. The Hamiltonian of the electronic system includes, in addition to the non-relativistic part, relativistic single-particle terms such as the Zeeman coupling to an external magnetic fields, spin-orbit coupling, and arbitrary magnetic anisotropies; the orbital degrees of freedom of the electrons are explicitly taken into account. We determine the complete relativistic exchange tensors, accounting for anisotropic exchange, Dzyaloshinskii-Moriya interactions, as well as additional non-diagonal symmetric terms (which may include dipole-dipole interaction). Our procedure provides the complete exchange tensors in a unified framework, including previously disregarded features such as the vertices of two-particle Green's functions and non-local self-energies. We do not assume any smallness in spin-orbit coupling, so our treatment is in this sense exact. Finally, we show how to distinguish and address separately the spin, orbital and spin-orbital contributions to magnetism.

preprint2015arXiv

Pressure induced crossing of the core-levels in 5d metals

Pressure induced interaction between core electrons, the core level crossing (CLC) transition has been observed in hcp Os at P~400 GPa [L. Dubrovinksy, et al., Nature 525, 226-229 (2015)]. In this work, we carry out a systematic study of the influence of pressure on the electronic structure in all metals of the 5d series (Hf,Ta,W,Re,Os,Ir,Pt,Au) using first-principles electronic structure calculations. We have found that CLC is a general effect for this series of metals. While in Pt it occurs at ~1500 GPa, at a pressure substantially higher than in Os, in Ir it occurs already at 80 GPa. Moreover, we predict that in Re the CLC transition may appear at ambient pressure. We analyze the shifts of the CLC transition pressure across the series within the Thomas-Fermi model, and show that the effect has many common features to the atomic collapse in the rare-earth elements.

preprint2015arXiv

Quantum theory as a description of robust experiments: derivation of the Pauli equation

It is shown that the Pauli equation and the concept of spin naturally emerge from logical inference applied to experiments on a charged particle under the conditions that (i) space is homogeneous (ii) the observed events are logically independent, and (iii) the observed frequency distributions are robust with respect to small changes in the conditions under which the experiment is carried out. The derivation does not take recourse to concepts of quantum theory and is based on the same principles which have already been shown to lead to e.g. the Schrödinger equation and the probability distributions of pairs of particles in the singlet or triplet state. Application to Stern-Gerlach experiments with chargeless, magnetic particles, provides additional support for the thesis that quantum theory follows from logical inference applied to a well-defined class of experiments.

preprint2015arXiv

Screening and plasmons in pure and disordered single- and bilayer black phosphorus

We study collective plasmon excitations and screening of disordered single- and bilayer black phosphorus beyond the low energy continuum approximation. The dynamical polarizability of phosphorene is computed using a tight-binding model that properly accounts for the band structure in a wide energy range. Electron-electron interaction is considered within the Random Phase Approximation. Damping of the plasmon modes due to different kinds of disorder, such as resonant scatterers and long-range disorder potentials, is analyzed. We further show that an electric field applied perpendicular to bilayer phosphorene can be used to tune the dispersion of the plasmon modes. For sufficiently large electric field, the bilayer BP enters in a topological phase with a characteristic plasmon spectrum, which is gaped in the armchair direction.

preprint2015arXiv

Spin and orbital exchange interactions from Dynamical Mean Field Theory

We derive a set of equations expressing the parameters of the magnetic interactions characterizing a strongly correlated electronic system in terms of single-electron Green's functions and self-energies. This allows to establish a mapping between the initial electronic system and a spin model including up to quadratic interactions between the effective spins, with a general interaction (exchange) tensor that accounts for anisotropic exchange, Dzyaloshinskii-Moriya interaction and other symmetric terms such as dipole-dipole interaction. We present the formulas in a format that can be used for computations via Dynamical Mean Field Theory algorithms.

preprint2015arXiv

Stripe glasses in ferromagnetic thin films

Domain walls in magnetic multilayered systems can exhibit a very complex and fascinating behavior. For example, the magnetization of thin films of hard magnetic materials is in general perpendicular to the thin-film plane, thanks to the strong out-of-plane anisotropy, but its direction changes periodically, forming an alternating spin-up and spin-down stripe pattern. The latter is stabilized by the competition between the ferromagnetic coupling and dipole-dipole interactions, and disappears when a moderate in-plane magnetic field is applied. It has been suggested that such a behavior may be understood in terms of a self-induced stripe glassiness. In this paper we show that such a scenario is compatible with the experimental findings. The strong out-of-plane magnetic anisotropy of the film is found to be beneficial for the formation of both the stripe-ordered and glassy phases. At zero magnetic field the system can form a glass only in a narrow interval of fairly large temperatures. An in-plane magnetic field, however, shifts the glass transition towards lower temperatures, therefore enabling it at or below room temperature. In good qualitative agreement with the experimental findings, we show that a moderate in-plane magnetic field of the order of $30~{\rm mT}$ can lead to the formation of defects in the stripe pattern, which sets the onset of the glass transition.

preprint2015arXiv

Thermodynamic consistency of the charge response in dynamical mean-field based approaches

We consider the thermodynamic consistency of the charge response function in the (extended) Hubbard model. In DMFT, thermodynamic consistency is preserved. We prove that the static, homogeneous DMFT susceptibility is consistent as long as vertex corrections obtained from the two-particle impurity correlation function are included. In presence of a nonlocal interaction, the problem may be treated within extended DMFT (EDMFT), or its diagrammatic extension, the dual boson approach. We show that here, maintaining thermodynamic consistency requires knowledge of three- and four-particle impurity correlation functions, which are typically neglected. Nevertheless, the dual boson approximation to the response is remarkably close to consistency. This holds even when two-particle vertex corrections are neglected. EDMFT is consistent only in the strongly correlated regime and near half-filling, where the physics is predominantly local.

preprint2015arXiv

Transport and optical properties of an electron gas in a Sierpinski carpet

Recent progress in the design and fabrication of artificial two-dimensional (2D) materials paves the way for the experimental realization of electron systems moving on plane fractals. In this work, we present the results of computer simulations for the conductance and optical absorption spectrum of a 2D electron gas roaming on a Sierpinski carpet, i.e. a plane fractal with Hausdorff dimension intermediate between one and two. We find that the conductance is sensitive to the spatial location of the leads and that it displays fractal fluctuations whose dimension is compatible with the Hausdorff dimension of the sample. Very interestingly, electrons in this fractal display a broadband optical absorption spectrum, which possesses sharp "molecular" peaks at low photon energies.

preprint2015arXiv

Ultralong-range order in the Fermi-Hubbard model with long-range interactions

We use the dual boson approach to reveal the phase diagram of the Fermi-Hubbard model with long-range dipole-dipole interactions. By using a large-scale finite-temperature calculation on a $64 \times 64$ square lattice we demonstrate the existence of a novel phase, possessing an `ultralong-range' order. The fingerprint of this phase -- the density correlation function -- features a non-trivial behavior on a scale of tens of the lattice sites. We study the properties and the stability of the ultralong-range ordered phase, and show that it is accessible in modern experiments with ultracold polar molecules and magnetic atoms.

preprint2014arXiv

Beyond extended dynamical mean-field theory: Dual boson approach to the two-dimensional extended Hubbard model

The dual boson approach [Ann. Phys. 327, 1320 (2012)] provides a means to construct a diagrammatic expansion around the extended dynamical mean-field theory (EDMFT). In this paper, we present the numerical implementation of the approach and apply it to the extended Hubbard model with nearest-neighbor interaction $V$. We calculate the EDMFT phase diagram and study the effect of diagrams beyond EDMFT on the transition to the charge-ordered phase. Including diagrammatic corrections to the EDMFT polarization shifts the EDMFT phase boundary to lower values of $V$. The approach interpolates between the random phase approximation in the weak coupling limit and EDMFT for strong coupling. Neglecting vertex corrections leads to results reminiscent of the EDMFT+$GW$ approximation. We however find significant deviations from the dual boson results already for small values of the interaction, emphasizing the crucial importance of fermion-boson vertex corrections.

preprint2014arXiv

Collective Charge Excitations of Strongly Correlated Electrons, Vertex Corrections and Gauge Invariance

We consider the collective, long-wavelength charge excitations in correlated media in presence of short- and long-range forces. As an example for the case of a short-range interaction, we examine the two-dimensional Hubbard model within dynamical mean-field theory (DMFT). It is shown explicitly that the DMFT susceptibility including vertex corrections respects the Ward identity and yields a manifestly gauge-invariant response in finite dimensions. For computing the susceptibility, we use a different expression and establish its formal equivalence to the standard DMFT formula. It allows for a more stable analytical continuation. We find a zero-sound mode expected for short-range forces. The relation between the vertex corrections, gauge invariance, and the appearance of the collective modes is discussed. Long-range forces are treated within extended dynamical mean-field theory. In order to obtain a gauge-invariant response, it is necessary to additionally incorporate some non local vertex corrections into the polarization. In doing so, we obtain plasmons in the three-dimensional Hubbard model. The plasma frequency is determined by the (single-particle) density distribution as a consequence of gauge invariance. We compare this result with the plasma frequency extracted from the analytical continuation of the susceptibility. It is in good agreement with the prediction from the gauge-invariance condition.

preprint2014arXiv

Fermi Condensation Near van Hove Singularities Within the Hubbard Model on the Triangular Lattice

The proximity of the Fermi surface to van Hove singularities drastically enhances interaction effects and leads to essentially new physics. In this work we address the formation of flat bands ("Fermi condensation") within the Hubbard model on the triangular lattice and provide a detailed analysis from an analytical and numerical perspective. To describe the effect we consider both weak-coupling and strong-coupling approaches, namely the renormalization group and dual fermion methods. It is shown that the band flattening is driven by correlations and is well pronounced even at sufficiently high temperatures, of the order of 0.1-0.2 of the hopping parameter. The effect can therefore be probed in experiments with ultracold fermions in optical lattices.

preprint2014arXiv

From Hubbard bands to spin-polaron excitations in the doped Mott material Na$_x$CoO$_2$

We investigate the excitation spectrum of strongly correlated sodium cobaltate within a realistic many-body description beyond dynamical mean-field theory (DMFT). At lower doping around $x$=0.3, rather close to Mott-critical half-filling, the single-particle spectral function of Na$_x$CoO$_2$ displays an upper Hubbard band which is captured within DMFT. Momentum-dependent self-energy effects beyond DMFT become dominant at higher doping. Around a doping level of $x\sim 0.67$, the incoherent excitations give way to finite-energy spin-polaron excitations in close agreement with optics experiments. These excitations are a direct consequence of the formation of bound states between quasiparticles and paramagnons in the proximity to in-plane ferromagnetic ordering.

preprint2014arXiv

Phonon-pump XUV-photoemission-probe in graphene: evidence for non-adiabatic heating of Dirac carriers by lattice deformation

We modulate the atomic structure of bilayer graphene by driving its lattice at resonance with the in-plane E1u lattice vibration at 6.3um. Using time- and angle-resolved photoemission spectroscopy (tr-ARPES) with extreme ultra-violet (XUV) pulses, we measure the response of the Dirac electrons near the K-point. We observe that lattice modulation causes anomalous carrier dynamics, with the Dirac electrons reaching lower peak temperatures and relaxing at faster rate compared to when the excitation is applied away from the phonon resonance or in monolayer samples. Frozen phonon calculations predict dramatic band structure changes when the E1u vibration is driven, which we use to explain the anomalous dynamics observed in the experiment.

preprint2014arXiv

Phonons and electron-phonon coupling in graphene-h-BN heterostructures

First principle calculations of the phonons of graphene-BN heterostructures are presented and compared to those of the constituents. We show that AA and AB' stacking are not only energetically less favoured than AB but also dynamically unstable. We have identified low energy flat phonon branches of BN character with out of plane displacement and evaluated their coupling to electrons in graphene.

preprint2014arXiv

Quantum capacitance and Landau parameters of massless Dirac fermions in graphene

We present extensive numerical results for the thermodynamic density of states (i.e. quantum capacitance) of a two-dimensional massless Dirac fermion fluid in a doped graphene sheet. In particular, by employing the random phase approximation, we quantify the impact of screening exerted by a metal gate located nearby a graphene flake. Finally, we calculate the spin- and circularly-symmetric Landau parameter, which can be experimentally extracted from independent measurements on the same setup of the quantum capacitance and quasiparticle velocity.

preprint2014arXiv

Quasiparticle band structure and tight-binding model for single- and bilayer black phosphorus

By performing \emph{ab initio} calculations for one- to four-layer black phosphorus within the $GW$ approximation, we obtain a significant difference in the band gap ($\sim$1.5 eV), which is in line with recent experimental data. The results are analyzed in terms of the constructed four-band tight-binding model, which gives accurate descriptions of the mono- and bilayer band structure near the band gap, and reveal an important role of the interlayer hoppings, which are largely responsible for the obtained gap difference.

preprint2014arXiv

Schottky barriers at hexagonal boron nitride/metal interfaces: a first principles study

The formation of a Schottky barrier at the interface between a metal and hexagonal boron nitride (h-BN) is studied using density functional theory. For metals whose work functions range from 4.2 to 6.0 eV, we find Schottky barrier heights for holes between 1.2 and 2.3 eV. A central role in determining the Schottky barrier height is played by a potential step of between 0.4 and 1.8 eV that is formed at the metal|h-BN interface and effectively lowers the metal work function. If h-BN is physisorbed, as is the case on fcc Cu, Al, Au, Ag and Pt(111) substrates, the interface potential step is described well by a universal function that depends only on the distance separating h-BN from the metal surface. The interface potential step is largest when h-BN is chemisorbed, which is the case for hcp Co and Ti (0001) and for fcc Ni and Pd (111) substrates.

preprint2014arXiv

Stiffening graphene by controlled defect creation

Graphene extraordinary strength, stiffness and lightness have generated great expectations towards its application in flexible electronics and as mechanical reinforcement agent. However, the presence of lattice defects, unavoidable in sheets obtained by scalable routes, might degrade its mechanical properties. Here we report a systematic study on the elastic modulus and strength of graphene with controlled density of defects. Counter intuitively, the in-plane Young modulus increases with increasing defect density up to almost twice the initial value for vacancy content of ~0.2%, turning it into the stiffest material ever reported. For higher density of vacancies, elastic modulus decreases with defect inclusion. The initial increase in Young modulus is explained in terms of a dependence of the elastic coefficients with the momentum of flexural modes predicted for 2D membranes. In contrast, the fracture strength decreases with defect density according to standard fracture continuum models. These quantitative structure-property relationships, measured in atmospheric conditions, are of fundamental and technological relevance and provide guidance for applications in which graphene mechanics represents a disruptive improvement.

preprint2013arXiv

Graphene as a Prototype Crystalline Membrane

The understanding of the structural and thermal properties of membranes, low-dimensional flexible systems in a space of higher dimension, is pursued in many fields from string theory to chemistry and biology. The case of a two-dimensional (2D) membrane in three dimensions is the relevant one for dealing with real materials. Traditionally, membranes are primarily discussed in the context of biological membranes and soft matter in general. The complexity of these systems hindered a realistic description of their interatomic structures based on a truly microscopic approach. Therefore, theories of membranes were developed mostly within phenomenological models. From the point of view of statistical mechanics, membranes at finite temperature are systems governed by interacting long-range fluctuations. Graphene, the first truly two-dimensional system consisting of just one layer of carbon atoms, provides a model system for the development of a microscopic description of membranes. In this Account, we review key results in the microscopic theory of structural and thermal properties of graphene and compare them with the predictions of phenomenological theories. The two approaches are in good agreement for the various scaling properties of correlation functions of atomic displacements. However, some other properties, such as the temperature dependence of the bending rigidity, cannot be understood based on phenomenological approaches. We also consider graphene at very high temperature and compare the results with existing models for two-dimensional melting. The melting of graphene presents a different scenario, and we describe that process as the decomposition of the graphene layer into entangled carbon chains.

preprint2013arXiv

Interaction phenomena in graphene seen through quantum capacitance

Capacitance measurements provide a powerful means of probing the density of states. The technique has proved particularly successful in studying 2D electron systems, revealing a number of interesting many-body effects. Here, we use large-area high-quality graphene capacitors to study behavior of the density of states in this material in zero and high magnetic fields. Clear renormalization of the linear spectrum due to electron-electron interactions is observed in zero field. Quantizing fields lead to splitting of the spin- and valley-degenerate Landau levels into quartets separated by interaction-enhanced energy gaps. These many-body states exhibit negative compressibility but the compressibility returns to positive in ultrahigh B. The reentrant behavior is attributed to a competition between field-enhanced interactions and nascent fractional states.

preprint2013arXiv

Origin of Anomalous Water Permeation through Graphene Oxide Membrane

Water inside the low dimensional carbon structures has been considered seriously owing to fundamental interest in its flow and structures as well as its practical impact. Recently, the anomalous perfect penetration of water through graphene oxide membrane was demonstrated although the membrane was impenetrable for other liquids and even gases. The unusual auxetic behavior of graphene oxide in the presence of water was also reported. Here, based on first-principles calculations, we establish atomistic models for hybrid systems composed of water and graphene oxides revealing the anomalous water behavior inside the stacked graphene oxides. We show that formation of hexagonal ice bilayer in between the flakes as well as melting transition of ice at the edges of flakes are crucial to realize the perfect water permeation across the whole stacked structures. The distance between adjacent layers that can be controlled either by oxygen reduction process or pressure is shown to determine the water flow thus highlighting a unique water dynamics in randomly connected two-dimensional spaces.

preprint2012arXiv

Atomically thin boron nitride: a tunnelling barrier for graphene devices

We investigate the electronic properties of heterostructures based on ultrathin hexagonal boron nitride (h-BN) crystalline layers sandwiched between two layers of graphene as well as other conducting materials (graphite, gold). The tunnel conductance depends exponentially on the number of h-BN atomic layers, down to a monolayer thickness. Exponential behaviour of I-V characteristics for graphene/BN/graphene and graphite/BN/graphite devices is determined mainly by the changes in the density of states with bias voltage in the electrodes. Conductive atomic force microscopy scans across h-BN terraces of different thickness reveal a high level of uniformity in the tunnel current. Our results demonstrate that atomically thin h-BN acts as a defect-free dielectric with a high breakdown field; it offers great potential for applications in tunnel devices and in field-effect transistors with a high carrier density in the conducting channel.

preprint2012arXiv

Electron-hole puddles in the absence of charged impurities

It is widely believed that carrier-density inhomogeneities ("electron-hole puddles") in single-layer graphene on a substrate like quartz are due to charged impurities located close to the graphene sheet. Here we demonstrate by using a Kohn-Sham-Dirac density-functional scheme that corrugations in a real sample are sufficient to determine electron-hole puddles on length scales that are larger than the spatial resolution of state-of-the-art scanning tunneling microscopy.

preprint2012arXiv

Equilibration and Thermalization of Classical Systems

It is demonstrated that the canonical distribution for a subsystem of a closed system follows directly from the solution of the time-reversible Newtonian equation of motion in which the total energy is strictly conserved. It is shown that this conclusion holds for both integrable or nonintegrable systems even though the whole system may contain as little as a few thousand particles. In other words, we demonstrate that the canonical distribution holds for subsystems of experimentally relevant sizes and observation times.

preprint2012arXiv

Lattice Expansion in Seamless Bi layer Graphene Constrictions at High Bias

Our understanding of sp2 carbon nanostructures is still emerging and is important for the development of high performance all carbon devices. For example, in terms of the structural behavior of graphene or bi-layer graphene at high bias, little to nothing is known. To this end we investigated bi-layer graphene constrictions with closed edges (seamless) at high bias using in situ atomic resolution transmission electron microscopy. We directly observe a highly localized anomalously large lattice expansion inside the constriction. Both the current density and lattice expansion increase as the bi-layer graphene constriction narrows. As the constriction width decreases below 10 nm, shortly before failure, the current density rises to 4 \cdot 109 A cm-2 and the constriction exhibits a lattice expansion with a uniaxial component showing an expansion approaching 5 % and an isotropic component showing an expansion exceeding 1 %. The origin of the lattice expansion is hard to fully ascribe to thermal expansion. Impact ionization is a process in which charge carriers transfer from bonding states to antibonding states thus weakening bonds. The altered character of C-C bonds by impact ionization could explain the anomalously large lattice expansion we observe in seamless bi-layer graphene constrictions. Moreover, impact ionization might also contribute to the observed anisotropy in the lattice expansion, although strain is probably the predominant factor.

preprint2012arXiv

Non-equilibrium magnetic interactions in strongly correlated systems

We formulate a low-energy theory for the magnetic interactions between electrons in the multi-band Hubbard model under non-equilibrium conditions determined by an external time-dependent electric field which simulates laser-induced spin dynamics. We derive expressions for dynamical exchange parameters in terms of non-equilibrium electronic Green functions and self-energies, which can be computed, e.g., with the methods of time-dependent dynamical mean-field theory. Moreover, we find that a correct description of the system requires, in addition to exchange, a new kind of magnetic interaction, that we name "twist exchange", which formally resembles Dzyaloshinskii-Moriya coupling, but is not due to spin-orbit, and is actually due to an effective three-spin interaction. Our theory allows the evaluation of the related time-dependent parameters as well.

preprint2012arXiv

Peculiarities of dynamics of Dirac fermions associated with zero-mass lines

Zero-mass lines result in appearance of linear dispersion modes for Dirac fermions. These modes play an important role in various physical systems. However, a Dirac fermion may not precisely follow a single zero-mass line, due to either tunneling between different lines or centrifugal forces. Being shifted from a zero-mass line the Dirac fermion acquires mass which can substantially influence its expected "massless" behavior. In the paper we calculate the energy gap caused by the tunneling between two zero-mass lines and show that its opening leads to the delocalization of linear dispersion modes. The adiabatic bending of a zero-mass line gives rise to geometric phases. These are the Berry phase, locally associated with a curvature, and a new phase resulting from the mass square asymmetry in the vicinity of a zero-mass line.

preprint2012arXiv

Polarization of graphene in a strong magnetic field beyond the Dirac cone approximation

In this paper we study the excitation spectrum of graphene in a strong magnetic field, beyond the Dirac cone approximation. The dynamical polarizability is obtained using a full $π$-band tight-binding model where the effect of the magnetic field is accounted for by means of the Peierls substitution. The effect of electron-electron interaction is considered within the random phase approximation, from which we obtain the dressed polarization function and the dielectric function. The range of validity of the Landau level quantization within the continuum approximation is studied, as well as the non-trivial quantization of the spectrum around the Van Hove singularity. We further discuss the effect of disorder, which leads to a smearing of the absorption peaks, and temperature, which activates additional inter-Landau level transitions induced by the Fermi distribution function.

preprint2012arXiv

Spectral flow for Aharonov--Bohm rings generated by zero-mass lines

We found that the spectrum of the linear dispersion mode localized along the zero-mass line for variable mass Dirac fermions shows the nonzero spectral flow as a function of the external magnetic field B. This happens due to the Aharonov--Bohm effect, leading to the linear dependence of every eigenvalue on B. The nonzero spectral flow allows to use the magnetic field slowly varying in time to control the energy of the linear dispersion mode.

preprint2011arXiv

Excitation spectrum and high energy plasmons in single- and multi-layer graphene

In this paper we study the excitation spectrum of single- and multi-layer graphene beyond the Dirac cone approximation. The dynamical polarizability of graphene is computed using a full $π$-band tight-binding model, considering the possibility of inter-layer hopping in the calculation. The effect of electron-electron interaction is considered within the random phase approximation. We further discuss the effect of disorder in the spectrum, which leads to a smearing of the absorption peaks. Our results show a redshift of the $π$-plasmon dispersion of single-layer graphene with respect to graphite, in agreement with experimental results. The inclusion of inter-layer hopping in the kinetic Hamiltonian of multi-layer graphene is found to be very important to properly capture the low energy region of the excitation spectrum.

preprint2011arXiv

Landau Level Spectrum of ABA- and ABC-stacked Trilayer Graphene

We study the Landau level spectrum of ABA- and ABC-stacked trilayer graphene. We derive analytic low energy expressions for the spectrum, the validity of which is confirmed by comparison to a π-band tight-binding calculation of the density of states on the honeycomb lattice. We further study the effect of a perpendicular electric field on the spectrum, where a zero-energy plateau appears for ABC stacking order, due to the opening of a gap at the Dirac point, while the ABA-stacked trilayer graphene remains metallic. We discuss our results in the context of recent electronic transport experiments. Furthermore, we argue that the expressions obtained can be useful in the analysis of future measurements of cyclotron resonance of electrons and holes in trilayer graphene.

preprint2011arXiv

Optical conductivity of disordered graphene beyond the Dirac cone approximation

In this paper we systemically study the optical conductivity and density of states of disorded graphene beyond the Dirac cone approximation. The optical conductivity of graphene is computed by using the Kubo formula, within the framework of a full π-band tight-binding model. Different types of non-correlated and correlated disorders are considered, such as random or Gaussian potentials, random or Gaussian nearest-neighbor hopping parameters, randomly distributed vacancies or their clusters, and random adsorbed hydrogen atoms or their clusters. For a large enough concentration of resonant impurities, a new peak in the optical conductivity is found, associated to transitions between the midgap states and the Van Hove singularities of the main π-band. We further discuss the effect of doping on the spectrum, and find that small amounts of resonant impurities are enough to obtain a background contribution to the conductivity in the infra-red part of the spectrum, in agreement with recent experiments.

preprint2011arXiv

Suppression of anharmonicities in crystalline membranes by external strain

In practice, physical membranes are exposed to a certain amount of external strain (tension or compression), due to the environment where they are placed. As a result, the behavior of the phonon modes of the membrane is modified. We show that anharmonic effects in stiff two-dimensional membranes are highly suppressed under the application of tension. For this, we consider the anharmonic coupling between bending and stretching modes in the self-consistent screening approximation (SCSA), and compare the obtained height-height correlation function in the SCSA to the corresponding harmonic propagator. The elasticity theory results are compared to atomistic Monte Carlo simulations for a graphene membrane under tension. We find that, while rather high values of strain are needed to avoid anharmonicity in soft membranes, strain fields less than 1% are enough to suppress all the anharmonic effects in stiff membranes, as graphene.

preprint2010arXiv

\textit{Ab initio} study of interacting lattice vibrations and stabilization of the $β$-phase in Ni-Ti shape-memory alloy

Lattice dynamical methods used to predict phase-transformations in crystals typically evaluate the harmonic phonon spectra and therefore do not work in frequent and important situations where the crystal structure is unstable in the harmonic approximation, such as the $β$ structure when it appears as a high-temperature phase of the shape memory alloy (SMA) NiTi. Here it is shown by self consistent {\it ab initio} lattice dynamical calculations (SCAILD) that the critical temperature for the pre-martensitic $R$ to $β$ phase-transformation in NiTi can be effectively calculated with good accuracy, and that the $β$-phase is a result primarily of the stabilizing interaction between different lattice vibrations.

preprint2010arXiv

Approach to equilbrium in nano-scale systems at finite temperatur

We study the time evolution of the reduced density matrix of a system of spin-1/2 particles interacting with an environment of spin-1/2 particles. The initial state of the composite system is taken to be a product state of a pure state of the system and a pure state of the environment. The latter pure state is prepared such that it represents the environment at a given finite temperature in the canonical ensemble. The state of the composite system evolves according to the time-dependent Schr{ö}dinger equation, the interaction creating entanglement between the system and the environment. It is shown that independent of the strength of the interaction and the initial temperature of the environment, all the eigenvalues of the reduced density matrix converge to their stationary values, implying that also the entropy of the system relaxes to a stationary value. We demonstrate that the difference between the canonical density matrix and the reduced density matrix in the stationary state increases as the initial temperature of the environment decreases. As our numerical simulations are necessarily restricted to a modest number of spin-1/2 particles ($<36$), but do not rely on time-averaging of observables nor on the assumption that the coupling between system and environment is weak, they suggest that the stationary state of the system directly follows from the time evolution of a pure state of the composite system, even if the size of the latter cannot be regarded as being close to the thermodynamic limit.

preprint2010arXiv

Complex edge effects in zigzag graphene nanoribbons due to hydrogen loading

We have performed density functional calculations as well as employed a tight-binding theory, to study the effect of passivation of zigzag graphene nanoribbons (ZGNR) by Hydrogen. We show that each edge C atom bonded with 2 H atoms open up a gap and destroys magnetism for small widths of the nanoribbon. However, a re-entrant magnetism accompanied by a metallic electronic structure is observed from 8 rows and thicker nanoribbons. The electronic structure and magnetic state are quite complex for this type of termination, with sp$^3$ bonded edge atoms being non-magnetic, whereas the nearest neighboring atoms are metallic and magnetic. We have also evaluated the phase stability of several thicknesses of ZGNR, and demonstrate that sp$^3$ bonded edge atoms, with 2 H atoms at the edge, should be stable at temperatures and pressures which are reachable in a laboratory environment.

preprint2010arXiv

Electronic Transport in Disordered Bilayer and Trilayer Graphene

We present a detailed numerical study of the electronic transport properties of bilayer and trilayer graphene within a framework of single-electron tight-binding model. Various types of disorder are considered, such as resonant (hydrogen) impurities, vacancies, short- or long-range Gaussian random potentials, and Gaussian random nearest neighbor hopping. The algorithms are based on the numerical solution of the time-dependent Schr ödinger equation and applied to calculate the density of states and conductivities (via the Kubo formula) of large samples containing millions of atoms. In the cases under consideration, far enough from the neutrality point, depending on the strength of disorders and the stacking sequence, a linear or sublinear electron-density dependent conductivity is found. The minimum conductivity $σ_{\min}\approx 2e^{2}/h$ (per layer) at the charge neutrality point is the same for bilayer and trilayer graphene, independent of the type of the impurities, but the plateau of minimum conductivity around the neutrality point is only observed in the presence of resonant impurities or vacancies, originating from the formation of the impurity band.

preprint2010arXiv

Modeling electronic structure and transport properties of graphene with resonant scattering centers

We present a detailed numerical study of the electronic properties of single-layer graphene with resonant ("hydrogen") impurities and vacancies within a framework of noninteracting tight-binding model on a honeycomb lattice. The algorithms are based on the numerical solution of the time-dependent Schrödinger equation and applied to calculate the density of states, \textit{quasieigenstates}, AC and DC conductivities of large samples containing millions of atoms. Our results give a consistent picture of evolution of electronic structure and transport properties of functionalized graphene in a broad range of concentration of impurities (from graphene to graphane), and show that the formation of impurity band is the main factor determining electrical and optical properties at intermediate impurity concentrations, together with a gap opening when approaching the graphane limit.

preprint2009arXiv

Origin of the Canonical Ensemble: Thermalization with Decoherence

We solve the time-dependent Schrodinger equation for the combination of a spin system interacting with a spin bath environment. In particular, we focus on the time development of the reduced density matrix of the spin system. Under normal circumstances we show that the environment drives the reduced density matrix to a fully decoherent state, and furthermore the diagonal elements of the reduced density matrix approach those expected for the system in the canonical ensemble. We show one exception to the normal case is if the spin system cannot exchange energy with the spin bath. Our demonstration does not rely on time-averaging of observables nor does it assume that the coupling between system and bath is weak. Our findings show that the canonical ensemble is a state that may result from pure quantum dynamics, suggesting that quantum mechanics may be regarded as the foundation of quantum statistical mechanics.