Source author record

M. I. Katsnelson

M. 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

207works

29topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Revised Tolmachev-Morel-Anderson pseudopotential for layered conventional superconductors with nonlocal Coulomb interaction

We study the effects of static nonlocal Coulomb interactions in layered conventional superconductors and show that they generically suppress superconductivity and reduce the critical temperature. Although the nonlocal Coulomb interaction leads to a significant structure in the superconducting gap function, we find that most properties can be effectively described by means of an appropriately revised local Coulomb pseudopotential $\tildeμ_C^*$, which is larger than the commonly adopted retarded Tolmachev-Morel-Anderson pseudopotential $μ^*_C$. To understand this, we analyze the Bethe-Salpeter equation describing the screening of Coulomb interaction in the superconducting state and obtain an expression for $\tildeμ_C^*$, which is valid in the presence of nonlocal Coulomb interactions. This analysis also reveals how the structure of the nonlocal Coulomb interaction weakens the screening effects from high-energy pair fluctuations and therefore yields larger values of the pseudopotential. Our findings are especially important for layered conventional superconductors with small Fermi energies and can be easily taken into account ab initio studies.

preprint2022arXiv

Spin dynamics of itinerant electrons: local magnetic moment formation and Berry phase

The state-of-the-art theoretical description of magnetic materials relies on solving effective Heisenberg spin problems or their generalizations to relativistic or multi-spin-interaction cases that explicitly assume the presence of local magnetic moments in the system. We start with a general interacting fermionic model that is often obtained in ab initio electronic structure calculations and show that the corresponding spin problem can be introduced even in the paramagnetic regime, which is characterized by a zero average value of the magnetization. Further, we derive a physical criterion for the formation of the local magnetic moment and confirm that the latter exists already at high temperatures well above the transition to the ordered magnetic state. The use of path-integral techniques allows us to disentangle spin and electronic degrees of freedom and to carefully separate rotational dynamics of the local magnetic moment from Higgs fluctuations of its absolute value. It also allows us to accurately derive the topological Berry phase and relate it to a physical bosonic variable that describes dynamics of the spin degrees of freedom. As the result, we demonstrate that the equation of motion in the case of a large magnetic moment takes a conventional Landau-Lifshitz form that explicitly accounts for the Gilbert damping due to itinerant nature of the original electronic model.

preprint2021arXiv

Certification of quantum states with hidden structure of their bitstrings

The rapid development of quantum computing technologies already made it possible to manipulate a collective state of several dozen of qubits. This success poses a strong demand on efficient and reliable methods for characterization and verification of large-scale many-body quantum states. Traditional methods, such as quantum tomography, which require storing and operating wave functions on classical computers, become problematic to use in the regime of large number of degrees of freedom. In this paper, we propose a numerically cheap procedure to describe and distinguish quantum states which is based on a limited number of simple projective measurements in at least two different bases and computing inter-scale dissimilarities of the resulting bit-string patterns via coarse-graining. The information one obtains through this procedure can be viewed as a "hash function" of quantum state -- a simple set of numbers which is specific for a concrete many-body wave function and can be used for certification. By studying a number of archetypal examples, we show that it is enough to characterize quantum states with different structure of entanglement, including the chaotic quantum states. The connection of the dissimilarity to standard measures of quantum correlations such as von Neumann entropy is discussed. We also demonstrate that our approach can be employed to detect phase transitions of different nature in many-body quantum magnetic systems.

preprint2021arXiv

Coexisting charge density wave and ferromagnetic instabilities in monolayer InSe

Recently fabricated InSe monolayers exhibit remarkable characteristics that indicate the potential of this material to host a number of many-body phenomena. Here, we consistently describe collective electronic effects in hole-doped InSe monolayers using advanced many-body techniques. To this end, we derive a realistic electronic-structure model from first principles that takes into account the most important characteristics of this material, including a flat band with prominent van Hove singularities in the electronic spectrum, strong electron-phonon coupling, and weakly-screened long-ranged Coulomb interactions. We calculate the temperature-dependent phase diagram as a function of band filling and observe that this system is in a regime with coexisting charge density wave and ferromagnetic instabilities that are driven by strong electronic Coulomb correlations. This regime can be achieved at realistic doping levels and high enough temperatures, and can be verified experimentally. We find that the electron-phonon interaction does not play a crucial role in these effects, effectively suppressing the local Coulomb interaction without changing the qualitative physical picture.

preprint2021arXiv

Local Plaquette Physics as Key Ingredient of High-Temperature Superconductivity in Cuprates

A major pathway towards understanding complex systems is given by exactly solvable reference systems that contain the essential physics of the system. For the $t-t'-U$ Hubbard model, the four-site plaquette is known to have a quantum critical point in the $U-μ$ space where states with electron occupations $N=2, 3, 4$ per plaquette are degenerate [Phys. Rev. B {\bf 94}, 125133 (2016)]. We show that such a critical point in the lattice causes an instability in the particle-particle singlet d-wave channel and manifests some of the essential elements of the cuprate superconductivity. For this purpose we design an efficient superperturbation theory -- based on the dual fermion approach -- with the critical plaquette as the reference system. Thus, the perturbation theory already contains the relevant d-wave fluctuations from the beginning via the two-particle correlations of the plaquette. We find that d-wave superconductivity remains a leading instability channel under reasonably broad range of parameters. The next-nearest-neighbour hopping $t'$ is shown to play a crucial role in a formation of strongly bound electronic bipolarons whose coherence at lower temperature results in superconductivity. The physics of the pseudogap within the developed picture is also discussed.

preprint2020arXiv

From local to non-local correlations: the Dual Boson perspective

Extended dynamical mean-field theory (EDMFT) is insufficient to describe non-local effects in strongly correlated systems, since corrections to the mean-field solution are generally large. We present an efficient scheme for the construction of diagrammatic extensions of EDMFT that avoids usual double counting problem by using an exact change of variables (the dual boson formalism) to distinguish the correlations included in the mean-field solution and those beyond. With a computational efficiency comparable to EDMFT+GW approach, our scheme significantly improves on the charge order transition phase boundary in the extended Hubbard model.

preprint2020arXiv

Limits on gas impermeability of graphene

Despite being only one-atom thick, defect-free graphene is considered to be completely impermeable to all gases and liquids. This conclusion is based on theory and supported by experiments that could not detect gas permeation through micrometre-size membranes within a detection limit of 10^5 to 10^6 atoms per second. Here, using small monocrystalline containers tightly sealed with graphene, we show that defect-free graphene is impermeable with an accuracy of eight to nine orders of magnitude higher than in the previous experiments. We could discern permeation of just a few helium atoms per hour, and this detection limit is also valid for all other tested gases (neon, nitrogen, oxygen, argon, krypton and xenon), except for hydrogen. Hydrogen shows noticeable permeation, even though its molecule is larger than helium and should experience a higher energy barrier. The puzzling observation is attributed to a two-stage process that involves dissociation of molecular hydrogen at catalytically active graphene ripples, followed by adsorbed atoms flipping to the other side of the graphene sheet with a relatively low activation energy of about 1.0 electronvolt, a value close to that previously reported for proton transport. Our work provides a key reference for the impermeability of two-dimensional materials and is important from a fundamental perspective and for their potential applications.

preprint2020arXiv

Magnetic Two-Dimensional Chromium Trihalides: A Theoretical Perspective

The discovery of ferromagnetic order in monolayer 2D crystals has opened a new venue in the field of two dimensional (2D) materials. 2D magnets are not only interesting on their own, but their integration in van der Waals heterostructures allows for the observation of new and exotic effects in the ultrathin limit. The family of Chromium trihalides, CrI$_3$, CrBr$_3$ and CrCl$_3$, is, so far, the most studied among magnetic 2D crystals. In this mini-review, we provide a perspective of the state of the art of the theoretical understanding of magnetic 2D trihalides, most of which will also be relevant for other 2D magnets, such as vanadium trihalides. We discuss both the well-established facts, such as the origin of the magnetic moment and magnetic anisotropy and address as well open issues such as the nature of the anisotropic spin couplings and the magnitude of the magnon gap. Recent theoretical predictions on Moir\' e magnets and magnetic skyrmions are also discussed. Finally, we give some prospects about the future interest of these materials and possible device applications.

preprint2020arXiv

Nonlinear optical study of commensurability effects in graphene-hBN heterostructures

Second-order nonlinear optical response allows to detect different properties of the system associated with the inversion symmetry breaking. Here, we use a second harmonic generation effect to investigate the alignment of a graphene/hexagonal Boron Nitride heterostructure. To achieve that, we activate a commensurate-incommensurate phase transition by a thermal annealing of the sample. We find that this structural change in the system can be directly observed through a strong modification of a nonlinear optical signal. This result reveals the potential of a second harmonic generation technique for probing structural properties of layered systems.

preprint2020arXiv

Origin of the vortex displacement field in twisted bilayer graphene

Model description of patterns of atomic displacements in twisted bilayer systems has been proposed. The model is based on the consideration of several dislocation ensembles, employing a language that is widely used for grain boundaries and film/substrate systems. We show that three ensembles of parallel screw dislocations are sufficient both to describe the rotation of the layers as a whole, and for the vortex-like displacements resulting from elastic relaxation. The results give a clear explanation of the observed features of the structural state such as vortices, accompanied by alternating stacking.

preprint2020arXiv

Probing the topology of the quantum analog of a classical skyrmion

In magnetism, skyrmions correspond to classical three-dimensional spin textures characterized by a topological invariant that keeps track of the winding of the magnetization in real space, a property that cannot be easily generalized to the quantum case since the orientation of a quantum spin is in general ill-defined. Moreover, as we show, the quantum skyrmion state cannot be directly observed in modern experiments that probe the local magnetization of the system. However, we show that this novel quantum state can still be identified and fully characterized by a special local three-spin correlation function defined on neighbouring lattice sites -- the scalar chirality -- which reduces to the classical topological invariant for large systems, and which is shown to be nearly constant in the quantum skyrmion phase.

preprint2020arXiv

Temperature-Induced Lifshitz Transition and Possible Excitonic Instability in ZrSiSe

The nodal-line semimetals have attracted immense interest due to the unique electronic structures such as the linear dispersion and the vanishing density of states as the Fermi energy approaching the nodes. Here, we report temperature-dependent transport and scanning tunneling microscope (spectroscopy) (STM[S]) measurements on nodal-line semimetal ZrSiSe.Our experimental results and theoretical analyses consistently demonstrate that the temperature induces Lifshitz transitions at 80 and 106 K in ZrSiSe, which results in the transport anomalies at the same temperatures. More strikingly, we observe a V-shaped dip structure around Fermi energy from the STS spectrum at low temperature,which can be attributed to co-effect of the spin-orbit coupling and excitonic instability. Our observations indicate the correlation interaction may play an important role in ZrSiSe, which owns the quasi-two-dimensional electronic structures.

preprint2020arXiv

Uncommon clustering in dilute Ti-Fe alloys

We present the results of ab initio modeling of structure of dilute Ti-Fe, a typical representative of quenched Ti-based transition-metal alloys. We have demonstrated that beyond the solubility limit this alloy cannot be described in common terms of substitutional and interstitial alloys. Instead, very stable local clusters are formed in both low-temperature hcp and high-temperature bcc phases of alloys, with almost identical local structures. This gives an example of geometrically frustrated state and explains unusual concentration behavior of Mössbauer spectra discovered long ago for this system.

preprint2019arXiv

Giant anisotropy of Gilbert damping in a Rashba honeycomb antiferromagnet

Giant Gilbert damping anisotropy is identified as a signature of strong Rashba spin-orbit coupling in a two-dimensional antiferromagnet on a honeycomb lattice. The phenomenon originates in spin-orbit induced splitting of conduction electron subbands that strongly suppresses certain spin-flip processes. As a result, the spin-orbit interaction is shown to support an undamped non-equilibrium dynamical mode that corresponds to an ultrafast in-plane Néel vector precession and a constant perpendicular-to-the-plane magnetization. The phenomenon is illustrated on the basis of a two dimensional $s$-$d$ like model. Spin-orbit torques and conductivity are also computed microscopically for this model. Unlike Gilbert damping these quantities are shown to reveal only a weak anisotropy that is limited to the semiconductor regime corresponding to the Fermi energy staying in a close vicinity of antiferromagnetic gap.

preprint2019arXiv

Magnetic polaron and antiferro-ferromagnetic transition in doped bilayer CrI$_3$

Gate-induced magnetic switching in bilayer CrI$_3$ has opened new ways for the design of novel low-power magnetic memories based on van der Waals heterostructures. The proposed switching mechanism seems to be fully dominated by electrostatic doping. Here we explain, by first-principle calculations, the ferromagnetic transition in doped bilayer CrI$_3$. For the case of a very small electron doping, our calculations predict the formation of magnetic polarons ("ferrons", "fluctuons") where the electron is self-locked in a ferromagnetic droplet in an antiferromagnetic insulating matrix. The self-trapping of holes is impossible, at least, within our approximation.

preprint2019arXiv

Orbitally-resolved ferromagnetism of monolayer CrI$_3$

Few-layer CrI$_3$ is the most known example among two-dimensional (2D) ferromagnets, which have attracted growing interest in recent years. Despite considerable efforts and progress in understanding the properties of 2D magnets both from theory and experiment, the mechanism behind the formation of in-plane magnetic ordering in chromium halides is still under debate. Here, we propose a microscopic orbitally-resolved description of ferromagnetism in monolayer CrI$_3$. Starting from first-principles calculations, we construct a low-energy model for the isotropic Heisenberg exchange interactions. We find that there are two competing contributions to the long-range magnetic ordering in CrI$_3$: (i) Antiferromagnetic Anderson's superexchange between half-filled $t_{2g}$ orbitals of Cr atoms; and (ii) Ferromagnetic exchange governed by the Kugel-Khomskii mechanism, involving the transitions between half-filled $t_{2g}$ and empty $e_g$ orbitals. Using numerical calculations, we estimate the exchange interactions in momentum-space, which allows us to restore the spin-wave spectrum, as well as estimate the Curie temperature. Contrary to the nearest-neighbor effective models, our calculations suggest the presence of sharp resonances in the spin-wave spectrum at 5--7 meV, depending on the vertical bias voltage. Our estimation of the Curie temperature in monolayer CrI$_3$ yields 55--65 K, which is in good agreement with experimental data.

preprint2016arXiv

Capturing non-local interaction effects in the Hubbard model: optimal mappings and limits of applicability

We investigate the Peierls-Feynman-Bogoliubov variational principle to map Hubbard models with nonlocal interactions to effective models with only local interactions. We study the renormalization of the local interaction induced by nearest-neighbor interaction and assess the quality of the effective Hubbard models in reproducing observables of the corresponding extended Hubbard models. We compare the renormalization of the local interactions as obtained from numerically exact determinant Quantum Monte Carlo to approximate but more generally applicable calculations using dual boson, dynamical mean field theory, and the random phase approximation. These more approximate approaches are crucial for any application with real materials in mind. Furthermore, we use the dual boson method to calculate observables of the extended Hubbard models directly and benchmark these against determinant Quantum Monte Carlo simulations of the effective Hubbard model.

preprint2016arXiv

Chirality-dependent transmission of spin waves through domain walls

Spin-wave technology (magnonics) has the potential to further reduce the size and energy consumption of information processing devices. In the submicrometer regime (exchange spin waves), topological defects such as domain walls may constitute active elements to manipulate spin waves and perform logic operations. We predict that spin waves that pass through a domain wall in an ultrathin perpendicular-anisotropy film experience a phase shift that depends on the orientation of the domain wall (chirality). The effect, which is absent in bulk materials, originates from the interfacial Dzyaloshinskii-Moriya interaction and can be interpreted as a geometric phase. We demonstrate analytically and by means of micromagnetic simulations that the phase shift is strong enough to switch between constructive and destructive interference. The two chirality states of the domain wall may serve as a memory bit or spin-wave switch in magnonic devices.

preprint2016arXiv

Correlation effects and orbital magnetism of Co clusters

Recent experiments on isolated Co clusters have shown huge orbital magnetic moments in comparison with their bulk and surface counterparts. These clusters hence provide the unique possibility to study the evolution of the orbital magnetic moment with respect to the cluster size and how competing interactions contribute to the quenching of orbital magnetism. We investigate here different theoretical methods to calculate the spin and orbital moments of Co clusters, and assess the performances of the methods in comparison with experiments. It is shown that density functional theory in conventional local density or generalized gradient approximations, or even with a hybrid functional, severely underestimates the orbital moment. As natural extensions/corrections we considered the orbital polarization correction, the LDA+U approximation as well as the LDA+DMFT method. Our theory shows that of the considered methods, only the LDA+DMFT method provides orbital moments in agreement with experiment, thus emphasizing the importance of dynamic correlations effects for determining fundamental magnetic properties of magnets in the nano-size regime.

preprint2016arXiv

Decoherence wave in magnetic systems and creation of Néel antiferromagnetic state by measurement

The interplay between the singlet ground state of the antiferromagnetic Heisenberg model and the experimentally measured Néel state of antiferromagnets is studied. To verify the hypothesis [M. I. Katsnelson et al., Phys. Rev. B 63, 212404 (2001)] that the latter can be considered to be a result of local measurements destroying the entanglement of the quantum ground state, we have performed systematic simulations of the effects of von Neumannmeasurements for the case of a one-dimensional antiferromagnetic spin-1/2 system for various types and degrees of magnetic anisotropies. It is found that in the ground state, a magnetization measurement can create decoherence waves [M. I. Katsnelson et al. Phys. Rev. A 62, 022118 (2000)] in the magnetic sublattices, and that a symmetry breaking anisotropy does not lead to alignment of the spins in a particular direction. However, for an easy-axis anisotropy of the same order magnitude as the exchange constant, a measurement on the singlet ground state can create Néel-ordering in finite systems of experimentally accessible size.

preprint2016arXiv

Friedel oscillations at the surfaces of rhombohedral $N$-layer graphene

The low-energy physics of rhombohedral $N$-layer graphene mainly arises on the external layers, where most of the π electrons are located. Their Bloch band structure defines a two-band semimetal; the dispersion relation scales as $\pm q^{N}$ with the momentum norm $q$ in the vicinity of two nonequivalent valleys. In this paper, we address the problem of elastic scattering through a localized impurity located either on the surface of the material or within the bulk, and focus on the quantum interferences it induces on the two external layers. It is apprehended in the framework of a $T$-matrix approach, both numerically and analytically, regardless of the impurity magnitude, which enables the description of realistic scatters. In rhombohedral multilayer graphene, the impurity induces Friedel oscillations that always decay as $1/r$. As a result, monolayer graphene is the only material of the rhombohedral class that exhibits $1/r^{2}$-decaying Friedel oscillations. The interference patterns are subsequently analyzed in momentum space. This analysis enables a clear distinction between monolayer graphene and multilayer graphene. It also shows that the interference pattern reveals the whole Bloch band structure, and highlights the number of layers stacked in the material, as well as the $π$-quantized Berry phases that characterize the existence of nodal points in the semimetallic spectrum. Experimentally, these features may be probed from scanning tunneling microscopy, when imaging the local density of states at the surfaces of suspended rhombohedral $N$-layer graphene.

preprint2016arXiv

Interaction-driven Lifshitz transition with dipolar fermions in optical lattices

Anisotropic dipole-dipole interactions between ultracold dipolar fermions break the symmetry of the Fermi surface and thereby deform it. Here we demonstrate that such a Fermi surface deformation induces a topological phase transition -- so-called Lifshitz transition -- in the regime accessible to present-day experiments. We describe the impact of the Lifshitz transition on observable quantities such as the Fermi surface topology, the density-density correlation function, and the excitation spectrum of the system. The Lifshitz transition in ultracold atoms can be controlled by tuning the dipole orientation and -- in contrast to the transition studied in crystalline solids -- is completely interaction-driven.

preprint2016arXiv

Intrinsic Charge Carrier Mobility in Single-Layer Black Phosphorus

We present a theory for single- and two-phonon charge carrier scattering in anisotropic two-dimensional semiconductors applied to single-layer black phosphorus (BP). We show that in contrast to graphene, where two-phonon processes due to the scattering by flexural phonons dominate at any practically relevant temperatures and are independent of the carrier concentration $n$, two-phonon scattering in BP is less important and can be considered negligible at $n\gtrsim10^{13}$ cm$^{-2}$. At smaller $n$, however, phonons enter in the essentially anharmonic regime. Compared to the hole mobility, which does not exhibit strong anisotropy between the principal directions of BP ($μ_{xx}/μ_{yy}\sim1.4$ at $n=10^{13}$ cm$^{-2}$ and $T=300$ K), the electron mobility is found to be significantly more anisotropic ($μ_{xx}/μ_{yy}\sim6.2$). Absolute values of $μ_{xx}$ do not exceed 250 (700) cm$^2$V$^{-1}$s$^{-1}$ for holes (electrons), which can be considered as an upper limit for the mobility in BP at room temperature.

preprint2016arXiv

Logical inference approach to relativistic quantum mechanics: derivation of the Klein-Gordon equation

The logical inference approach to quantum theory, proposed earlier [Ann. Phys. 347 (2014) 45-73], is considered in a relativistic setting. It is shown that the Klein-Gordon equation for a massive, charged, and spinless particle derives from the combination of the requirements that the space-time data collected by probing the particle is obtained from the most robust experiment and that on average, the classical relativistic equation of motion of a particle holds.

preprint2016arXiv

Macroscopic self-reorientation of interacting two-dimensional crystals

Microelectromechanical systems, which can be moved or rotated with nanometre precision, already find applications in such fields like radio-frequency electronics, micro-attenuators, sensors and many others. Especially interesting are those which allow fine control over the motion on atomic scale due to self-alignment mechanisms and forces acting on the atomic level. Such machines can produce well-controlled movements as a reaction to small changes of the external parameters. Here we demonstrate that, for the system of graphene on hexagonal boron nitride, the interplay between the van der Waals and elastic energies results in graphene mechanically self-rotating towards the hexagonal boron nitride crystallographic directions. Such rotation is macroscopic (for graphene flakes of tens of micrometres the tangential movement can be on hundreds of nanometres) and can be used for reproducible manufacturing of aligned van der Waals heterostructures.

preprint2016arXiv

Many-body effects on Cr(001) surfaces: An LDA+DMFT study

The electronic structure of the Cr(001) surface with its sharp resonance at the Fermi level is a subject of controversial debate of many experimental and theoretical works. To date, it is unclear whether the origin of this resonance is an orbital Kondo or an electron-phonon coupling effect. We have combined ab initio density functional calculations with dynamical mean-field simulations to calculate the orbitally resolved spectral function of the Cr(001) surface. The calculated orbital character and shape of the spectrum is in agreement with data from (inverse) photoemission experiments. We find that dynamic electron correlations crucially influence the surface electronic structure and lead to a low energy resonance in the $d_{z^2}$ and $d_{xz/yz}$ orbitals. Our results help to reconvene controversial experimental results from (I)PES and STM measurements.

preprint2016arXiv

Many-body effects on graphene conductivity: Quantum Monte Carlo calculations

Optical conductivity of graphene is studied using Quantum Monte Carlo calculations. We start from Euclidean current-current correlator and extract $σ(ω)$ from Green-Kubo relations using Backus-Gilbert method. Calculations were performed both for long-range interactions and taking into account only contact term. In both cases we vary interaction strength and study its influence on optical conductivity. We compare our results with previous theoretical calculations choosing $ω\approx κ$ thus working in the region of the plateau in $σ(ω)$ which corresponds to optical conductivity of Dirac quasiparticles. No dependence of optical conductivity on interaction strength is observed unless we approach antiferromagnetic phase transition in case of artificially enhanced contact term. Our results strongly support previous theoretical studies claimed very weak regularization of graphene conductivity.

preprint2016arXiv

Midgap states and band gap modification in defective graphene/h-BN heterostructures

The role of defects in van der Waals heterostructures made of graphene and hexagonal boron nitride (h-BN) is studied by a combination of ab initio and model calculations. Despite the weak van der Waals interaction between layers, defects residing in h-BN, such as carbon impurities and antisite defects, reveal a hybridization with graphene p$_{\rm z}$ states, leading to midgap state formation. The induced midgap states modify the transport properties of graphene and can be reproduced by means of a simple effective tight-binding model. In contrast to carbon defects, it is found that oxygen defects do not strongly hybridize with graphene's low-energy states. Instead, oxygen drastically modifies the band gap of graphene, which emerges in a commensurate stacking on h-BN lattices.

preprint2016arXiv

Monte-Carlo study of Dirac semimetals phase diagram

In this paper the phase diagram of Dirac semimetals is studied within lattice Monte-Carlo simulation. In particular, we concentrate on the dynamical chiral symmetry breaking which results in semimetal/insulator transition. Using numerical simulation we determined the values of the critical coupling constant of the semimetal/insulator transition for different values of the anisotropy of the Fermi velocity. This measurement allowed us to draw tentative phase diagram for Dirac semimetals. It turns out that within the Dirac model with Coulomb interaction both Na$_3$Bi and Cd$_3$As$_2$ known experimentally to be Dirac semimetals would lie deeply in the insulating region of the phase diagram. It probably shows a decisive role of screening of the interelectron interaction in real materials, similar to the situation in graphene.

preprint2016arXiv

Plaquette Valence Bond Theory of High-Temperature Superconductivity

We present a strong-coupling approach to the theory of high-temperature superconductivity based on the observation of a quantum critical point in the plaquette within the t,t' Hubbard model. The crossing of ground state energies in the N=2,3,4 sectors occurs for parameters close to the optimal doping. The theory predicts the maximum of the d(x2-y2)-wave order parameter at the border between localized and itinerant electron behavior and gives a natural explanation for the pseudo-gap formation via soft-fermion mode related to local singlet states of the plaquette in the environment. Our approach follows the general line of resonating valence bond theory stressing a crucial role of singlets in the physics of high-Tc superconductors, but focuses on the formation of local singlets, similar to phenomena observed in frustrated one-dimensional quantum spin models.

preprint2016arXiv

Polaronic effects in monolayer black phosphorus on polar substrates

We investigate the effect of charge carrier interaction with surface optical phonons on the band properties of monolayer black phosphorus induced by polar substrates. We develop an analytical method based on the Lee-Low-Pines theory to calculate the spectrum of Fröhlich type continuum Hamiltonian in the long-wavelength limit. We examine the modification of a band gap and renormalization of effective masses due to the substrate-related polaronic effect. Our results show that an energy gap in supported monolayer black phosphorus is enlarged depending on a particular substrate and the interlayer distance, $z$. Among the substrate considered, the largest gap broadening at $ z=2.5$ Å is observed for the Al$_{2} $O$_{3} $ substrate, which is found to be $ \sim 50$ meV. Carrier-phonon coupling also renormalizes the effective masses which is more pronounced along the zigzag direction. Anisotropy of the effective masses becomes stronger by the influence of the polaronic effect corresponding to direction-dependent carrier-phonon coupling. We conclude that substrate phonons have a non-negligible effect on the static band properties of monolayer black phosphorus, which may be further exploited in its experimental and theoretical studies.

preprint2016arXiv

Quantum elasticity of graphene: Thermal expansion coefficient and specific heat

We explore thermodynamics of a quantum membrane, with a particular application to suspended graphene membrane and with a particular focus on the thermal expansion coefficient. We show that an interplay between quantum and classical anharmonicity-controlled fluctuations leads to unusual elastic properties of the membrane. The effect of quantum fluctuations is governed by the dimensionless coupling constant, $g_0 \ll 1$, which vanishes in the classical limit ($\hbar \to 0$) and is equal to $\simeq 0.05$ for graphene. We demonstrate that the thermal expansion coefficient $α_T$ of the membrane is negative and remains nearly constant down to extremely low temperatures, $T_0\propto \exp (-2/g_0)$. We also find that $α_T$ diverges in the classical limit: $α_T \propto - \ln(1/g_0)$ for $g_0 \to 0$. For graphene parameters, we estimate the value of the thermal expansion coefficient as $α_T \simeq - 0.23\:{\rm eV}^{-1}$, which applies below the temperature $T_{\rm uv} \sim g_0 \varkappa_0 \sim 500$\:K (where $\varkappa_0 \sim 1$\:eV is the bending rigidity) down to $T_0 \sim 10^{-14}$\:K. For $T<T_0$, the thermal expansion coefficient slowly (logarithmically) approaches zero with decreasing temperature. This behavior is surprising since typically the thermal expansion coefficient goes to zero as a power-law function. We discuss possible experimental consequences of this anomaly. We also evaluate classical and quantum contributions to the specific heat of the membrane and investigate the behavior of the Grüneisen parameter.

preprint2016arXiv

Quantum theory as plausible reasoning applied to data obtained by robust experiments

We review recent work that employs the framework of logical inference to establish a bridge between data gathered through experiments and their objective description in terms of human-made concepts. It is shown that logical inference applied to experiments for which the observed events are independent and for which the frequency distribution of these events is robust with respect to small changes of the conditions under which experiments are carried out yields, without introducing any concept of quantum theory, the quantum theoretical description in terms of the Schrödinger or the Pauli equation, the Stern-Gerlach or Einstein-Podolsky-Rosen-Bohm experiments. The extraordinary descriptive power of quantum theory then follows from the fact that it is plausible reasoning, that is common sense, applied to reproducible and robust experimental data.

preprint2016arXiv

Role of direct exchange and Dzyaloshinskii-Moriya interactions in magnetic properties of graphene derivatives: C$_2$F and C$_2$H

According to the Lieb's theorem the ferromagnetic interaction in graphene-based materials with bipartite lattice is a result of disbalance between the number of sites available for $p_z$ electrons in different sublattices. Here, we report on another mechanism of the ferromagnetism in functionalized graphene that is the direct exchange interaction between spin orbitals. By the example of the single-side semihydrogenated (C$_2$H) and semifluorinated (C$_2$F) graphene we show that such a coupling can partially or even fully compensate antiferromagnetic character of indirect exchange interactions reported earlier [Phys. Rev. B {\bf 88}, 081405(R) (2013)]. As a result, C$_2$H is found to be a two-dimensional material with the isotropic ferromagnetic interaction and negligibly small magnetic anisotropy, which prevents the formation of the long-range magnetic order at finite temperature in accordance with the Mermin-Wagner theorem. This gives a rare example of a system where direct exchange interactions play a crucial role in determining a magnetic structure. In turn, C$_2$F is found to be at the threshold of the antiferromagnetic-ferromagnetic instability, which in combination with the Dzyaloshinskii-Moriya interaction can lead to a skyrmion state.

preprint2016arXiv

Self-consistent Dual Boson approach to single-particle and collective excitations in correlated systems

We propose an efficient dual boson scheme, which extends the DMFT paradigm to collective excitations in correlated systems. The theory is fully self-consistent both on the one- and on the two-particle level, thus describing the formation of collective modes as well as the renormalization of electronic and bosonic spectra on equal footing. The method employs an effective impurity model comprising both fermionic and bosonic hybridization functions. Only single- and two-electron Green's functions of the reference problem enter the theory, due to the optimal choice of the self-consistency condition for the effective bosonic bath. We show that the theory is naturally described by a dual Luttinger-Ward functional and obeys the relevant conservation laws.

preprint2016arXiv

Spin-orbit coupling and magnetic interactions in Si(111):{C,Si,Sn,Pb}

We study the magnetic properties of the adatom systems on a semiconductor surface Si(111):\{C,Si,Sn,Pb\} - ($\sqrt{3} \times \sqrt{3}$). On the basis of all-electron density functional theory calculations we construct effective low-energy models taking into account spin-orbit coupling and electronic correlations. In the ground state the surface nanostructures are found to be insulators with the non-collinear 120$^{\circ}$ Néel (for C, Si, Sn monolayer coverages) and 120$^{\circ}$ row-wise (for Pb adatom) antiferromagnetic orderings. The corresponding spin Hamiltonians with anisotropic exchange interactions are derived by means of the superexchange theory and the calculated Dzyaloshinskii-Moriya interactions are revealed to be very strong and compatible with the isotropic exchange couplings in the systems with Sn and Pb adatoms. To simulate the excited magnetic states we solve the constructed spin models by means of the Monte Carlo method. At low temperatures and zero magnetic field we observe complex spin spiral patterns in Sn/Si(111) and Pb/Si(111). On this basis the formation of antiferromagnetic skyrmion lattice states in adatom $sp$ electron systems in strong magnetic fields is discussed.

preprint2016arXiv

Standard model of the rare-earths, analyzed from the Hubbard I approximation

In this work we examine critically the electronic structure of the rare-earth elements by use of the so-called Hubbard I approximation. From the theoretical side all measured features of both occupied and unoccupied states are reproduced, without significant deviations between observations and theory. We also examine cohesive properties like the equilibrium volume and bulk modulus, where we find, in general, a good agreement between theory and measurements. In addition we have reproduced the spin and orbital moments of these elements, as they are reflected from measurements of the saturation moment. We have also employed the Hubbard I approximation to extract the interatomic exchange parameters of an effective spin Hamiltonian for the heavy rare earths. We show that the Hubbard I approximation gives results which are consistent with calculations where $4f$ electrons are treated as core states for Gd. The latter approach was also used to address the series of the heavy/late rare-earths. Via Monte Carlo simulations we obtained ordering temperatures which reproduce measurements within about $20\%$. We have further illustrated the accuracy of these exchange parameters by comparing measured and calculated magnetic configurations for the heavy rare earths and the magnon dispersion for Gd. The Hubbard I approximation is compared to other theories of the electronic structure, and we argue that it is superior. We discuss the relevance of our results in general, and how this makes it possible to treat the electronic structure of materials containing rare-earth elements, such as permanent magnets, magnetostrictive compounds, photovoltaics, optical fibers, topological insulators, and molecular magnets.

preprint2016arXiv

Treatment of 4f states of the rare-earths: the case study of TbN

The lattice constant, bulk modulus and shear constant of TbN are calculated by means of density functional theory (DFT) in the local density approximation (LDA) and generalized gradient approximation (GGA), with 4f-states treated as valence electrons or core electrons. In addition, local Coulomb repulsions U are treated both statically as in the LDA+U approach and dynamically as in the dynamical mean-field theory (DMFT) in the Hubbard-I approximation. It is shown that all methods, except DFT-LDA with 4f electrons treated as either valence states, produce lattice constants and bulk moduli in good agreement with experiment. In the LDA+U approach multiple minima are found, and we focus on the competition between a state with cubic symmetry and a state obtained from atomic Hund's rules. We find the state with cubic symmetry to be 0.59 eV lower in energy than the Hund's rules state, while the opposite was obtained in previous literature. The shear constant is shown to be rather sensitive to the theoretical method used, and the Hund's rules state obtained in LDA+U is found to be unstable towards tetragonal shear. As to the magnetism, we find that the calculation based on the Hubbard-I approximation reproduces observations with the best accuracy. Finally, the spectral properties of TbN are discussed, together with the general applicability of the different methods in describing rare-earth elements and compounds.

preprint2015arXiv

Chemical modifications and stability of phosphorene with impurities: A first principles study

We perform a systematic first-principles study of phosphorene in the presence of typical monovalent (hydrogen, fluorine) and divalent (oxygen) impurities. The results of our modeling suggest a decomposition of phosphorene into weakly bonded one-dimensional (1D) chains upon single- and double-side hydrogenation and fluorination. In spite of a sizable quasiparticle band gap (2.29 eV), fully hydrogenated phosphorene found to be dynamically unstable. In contrast, full fluorination of phosphorene gives rise to a stable structure, being an indirect gap semiconductor with the band gap of 2.27 eV. We also show that fluorination of phosphorene from the gas phase is significantly more likely than hydrogenation due to the relatively low energy barrier for the dissociative adsorption of F2 (0.19 eV) compared to H2 (2.54 eV). At low concentrations, monovalent impurities tend to form regular atomic rows phosphorene, though such patterns do not seem to be easily achievable due to high migration barriers (1.09 and 2.81 eV for H2 and F2, respectively). Oxidation of phosphorene is shown to be a qualitatively different process. Particularly, we observe instability of phosphorene upon oxidation, leading to the formation of disordered amorphous-like structures at high concentrations of impurities.

preprint2015arXiv

Chiral tunneling through generic one-dimensional potential barriers in bilayer graphene

We study tunneling of charge carriers in single- and bilayer graphene. We propose an explanation for non-zero "magic angles" with 100% transmission for the case of symmetric potential barrier, as well as for their almost-survival for slightly asymmetric barrier in the bilayer graphene known previously from numerical simulations. Most importantly, we demonstrate that these magic angles are not protected in the case of bilayer and give an explicit example of a barrier with very small electron transmission probability for any angles. This means that one can lock charge carriers by a p-n-p (or n-p-n) junction without opening energy gap. This creates new opportunities for the construction of graphene transistors.

preprint2015arXiv

Dirac electrons and domain walls: a realization in junctions of ferromagnets and topological insulators

We study a system of Dirac electrons with finite density of charge carriers coupled to an external electromagnetic field in two spatial dimensions, with a domain wall (DW) mass term. The interface between a thin-film ferromagnet and a three-dimensional topological insulator provides a condensed-matter realization of this model, when an out-of-plane domain wall magnetization is coupled to the TI surface states. We show how, for films with very weak intrinsic in-plane anisotropies, the torque generated by the edge electronic current flowing along the DW competes with an effective in-plane anisotropy energy, induced by quantum fluctuations of the chiral electrons bound to the wall, in a mission to drive the internal angle of the DW from a Bloch configuration towards a Néel configuration. Both the edge current and the induced anisotropy contribute to stabilize the internal angle, so that for weak intrinsic in-plane anisotropies DW motion is still possible without suffering from an extremely early Walker breakdown.

preprint2015arXiv

Effect of structural relaxation on the electronic structure of graphene on hexagonal boron nitride

We performed calculations of electronic, optical and transport properties of graphene on hBN with realistic moiré patterns. The latter are produced by structural relaxation using a fully atomistic model. This relaxation turns out to be crucially important for electronic properties. We describe experimentally observed features such as additional Dirac points and the "Hofstadter butterfly" structure of energy levels in a magnetic field. We find that the electronic structure is sensitive to many-body renormalization of the local energy gap.

preprint2015arXiv

Effective Hamiltonians for fastly driven many-body lattice systems

We consider 1D lattices described by Hubbard or Bose-Hubbard models, in the presence of periodic high-frequency perturbations, such as uniform ac force or modulation of hopping coefficients. Effective Hamiltonians for interacting particles are derived using an averaging method resembling classical canonical perturbation theory. As is known, a high-frequency force may renormalize hopping coefficients, causing interesting phenomena such as coherent destruction of tunnelling and creation of artificial gauge fields. We find explicitly additional corrections to the effective Hamiltonians due to interactions, corresponding to non-trivial processes such as single-particle density-dependent tunnelling, correlated pair hoppings, nearest neighbour interactions, etc. Some of these processes arise also in multiband lattice models, and are capable to give rise to a rich variety of quantum phases. The apparent contradiction with other methods, e.g. Floquet-Magnus expansion, is explained. The results may be useful for designing effective Hamiltonian models in experiments with ultracold atoms, as well as in the field of ultrafast nonequilibrium magnetism. An example of manipulating exchange interaction in a Mott-Hubbard insulator is considered, where our corrections play an essential role.

preprint2015arXiv

Exchange parameters of strongly correlated materials: extraction from spin-polarised density functional theory plus dynamical mean field theory

In this paper we present an accurate numerical scheme for extracting inter-atomic exchange parameters ($J_{ij}$) of strongly correlated systems, based on first-principles full-potential electronic structure theory. The electronic structure is modelled with the help of a full-potential linear muffin-tin orbital method. The effects of strong electron correlations are considered within the charge self-consistent density functional theory plus dynamical mean-field theory (DFT+DMFT). The exchange parameters are then extracted using the magnetic force theorem, hence all the calculations are performed within a single computational framework. The method allows to investigate how the $J_{ij}$-parameters are affected by dynamical electron correlations. In addition to describing the formalism and details of the implementation, we also present magnetic properties of a few commonly discussed systems, characterised by different degrees of electron localisation. In bcc Fe we found a minor renormalisation of the $J_{ij}$ interactions once the dynamical correlations are introduced. However, generally, if the magnetic coupling has several competing contributions from different orbitals, the redistribution of the spectral weight and changes in the exchange splitting of these states can lead to a dramatic modification of the total interaction parameter. In NiO we found that both static and dynamical mean-field results provide an adequate description of the exchange interactions, which is somewhat surprising given the fact that these two methods result in quite different electronic structures. By employing Hubbard-I approximation for the treatment of the $4f$ states in hcp Gd we reproduce the experimentally observed multiplet structure. The calculated exchange parameters result to be rather close to the ones obtained by treating the $4f$ electrons as non-interacting core states.

preprint2015arXiv

Fingerprints of disorder in graphene

We present a systematic study of the electronic, transport and optical properties of disordered graphene including the next-nearest-neighbor hopping. We show that this hopping has a non-negligible effect on resonant scattering but is of minor importance for long-range disorder such as charged impurities, random potentials or hoppings induced by strain fluctuations. Different types of disorders can be recognized by their fingerprints appearing in the profiles of dc conductivity, carrier mobility, optical spectroscopy and Landau level spectrum. The minimum conductivity $4e^{2}/h$ found in the experiments is dominated by long-range disorder and the value of $4e^{2}/πh$ is due to resonant scatterers only.

preprint2015arXiv

Geometric, electronic and magnetic structure of Fe$_{x}$O$_{y}^{+}$ clusters

Correlation between geometry, electronic structure and magnetism of solids is both intriguing and elusive. This is particularly strongly manifested in small clusters, where a vast number of unusual structures appear. Here, we employ density functional theory in combination with a genetic search algorithm, GGA$+U$ and a hybrid functional to determine the structure of gas phase Fe$_{x}$O$_{y}^{+/0}$ clusters. For Fe$_{x}$O$_{y}$ cation clusters we also calculate the corresponding vibration spectra and compare them with experiments. We successfully identify Fe$_{3}$O$_{4}^{+}$, Fe$_{4}$O$_{5}^{+}$, Fe$_{4}$O$_{6}^{+}$, Fe$_{5}$O$_{7}^{+}$ and propose structures for Fe$_{6}$O$_{8}^{+}$. Within the triangular geometric structure of Fe$_{3}$O$_{4}^{+}$ a non-collinear, ferrimagnetic and ferromagnetic state are comparable in energy. Fe$_{4}$O$_{5}^{+}$ and Fe$_{4}$O$_{6}^{+}$ are ferrimagnetic with a residual magnetic moment of 1~\muB{} due to ionization. Fe$_{5}$O$_{7}^{+}$ is ferrimagnetic due to the odd number of Fe atoms. We compare the electronic structure with bulk magnetite and find Fe$_{4}$O$_{5}^{+}$, Fe$_{4}$O$_{6}^{+}$, Fe$_{6}$O$_{8}^{+}$ to be mixed valence clusters. In contrast, in Fe$_{3}$O$_{4}^{+}$ and Fe$_{5}$O$_{7}^{+}$ all Fe are found to be trivalent.

preprint2015arXiv

Landau levels of single layer and bilayer phosphorene

In this work we introduce a low-energy Hamiltonian for single layer and bilayer black phosphorus that describes the electronic states at the vicinity of the gamma point. The model is based on a recently proposed tight-binding description for electron and hole bands close to the Fermi level. We calculate expressions for the Landau level spectrum as function of magnetic field and in the case of bilayer black phosphorus we investigate the effect of an external bias on the electronic band gap. The results showcase the highly anisotropic character of black phosphorus and in particular for bilayer BP, the presence of bias allows for a field-induced semiconductor-metal transition.

preprint2015arXiv

Layer-resolved magnetic exchange interactions of surfaces of late 3d elements: effects of electronic correlations

We present the results of an ab initio study of magnetic properties of Fe, Co and Ni surfaces. In particular, we discuss their electronic structure and magnetic exchange interactions (Jij), as obtained by means of a combination of density functional theory and dynamical mean-field theory. All studied systems have a pronounced tendency to ferromagnetism both for bulk and surface atoms. The presence of narrow-band surface states is shown to enhance the magnetic moment as well as the exchange couplings. The most interesting results were obtained for the Fe surface where the atoms have a tendency to couple antiferromagnetically with each other. This interaction is relatively small, when compared to interlayer ferromagnetic interaction, and strongly depends on the lattice parameter. Local correlation effects are shown to lead to strong changes of the overall shape of the spectral functions. However, they seem to not play a decisive role on the overall picture of the magnetic couplings studied here. We have also investigated the influence of correlations on the spin and orbital moments of the bulk-like and surface atoms. We found that dynamical correlations in general lead to enhanced values of the orbital moment.

preprint2015arXiv

Magnetism and interaction-induced gap opening in graphene with vacancies or hydrogen adatoms: Quantum Monte Carlo study

We study electronic properties of graphene with finite concentration of vacancies or other resonant scatterers by a straightforward lattice Quantum Monte Carlo calculations. Taking into account realistic long-range Coulomb interaction we calculate distribution of spin density associated to midgap states and demonstrate antiferromagnetic ordering. Energy gaps are open due to the interaction effects, both in the bare graphene spectrum and in the vacancy/impurity bands. In the case of 5 % concentration of resonant scatterers the latter gap is estimated as 0.7 eV and 1.1 eV for graphene on boron nitride and freely suspended graphene, respectively.

preprint2015arXiv

Many-body effects in graphene beyond the Dirac model with Coulomb interaction

This paper is devoted to development of perturbation theory for studying the properties of graphene sheet of finite size, at nonzero temperature and chemical potential. The perturbation theory is based on the tight-binding Hamiltonian and arbitrary interaction potential between electrons, which is considered as a perturbation. One-loop corrections to the electron propagator and to the interaction potential at nonzero temperature and chemical potential are calculated. One-loop formulas for the energy spectrum of electrons in graphene, for the renormalized Fermi velocity and also for the dielectric permittivity are derived.

preprint2015arXiv

Melting temperature of graphene

We present an approach to the melting of graphene based on nucleation theory for a first order phase transition from the 2D solid to the 3D liquid via an intermediate quasi-2D liquid. The applicability of nucleation theory, supported by the results of systematic atomistic Monte Carlo simulations, provides an intrinsic definition of the melting temperature of graphene, $ T_m $, and allows us to determine it. We find $T_m \simeq 4510$ K, about 250 K higher than that of graphite using the same interatomic interaction model. The found melting temperature is shown to be in good agreement with the asymptotic results of melting simulations for finite disks and ribbons of graphene. Our results strongly suggest that graphene is the most refractory of all known materials.

preprint2015arXiv

Microscopic origin of Heisenberg and non-Heisenberg exchange interactions in ferromagnetic bcc Fe

By means of first principles calculations we investigate the nature of exchange coupling in ferromagnetic bcc Fe on a microscopic level. Analyzing the basic electronic structure reveals a drastic difference between the $3d$ orbitals of $E_g$ and $T_{2g}$ symmetries. The latter ones define the shape of the Fermi surface, while the former ones form weakly-interacting impurity levels. We demonstrate that, as a result of this, in Fe the $T_{2g}$ orbitals participate in exchange interactions, which are only weakly dependent on the configuration of the spin moments and thus can be classified as Heisenberg-like. These couplings are shown to be driven by Fermi surface nesting. In contrast, for the $E_g$ states the Heisenberg picture breaks down, since the corresponding contribution to the exchange interactions is shown to strongly depend on the reference state they are extracted from. Our analysis of the nearest-neighbour coupling indicates that the interactions among $E_g$ states are mainly proportional to the corresponding hopping integral and thus can be attributed to be of double-exchange origin.

preprint2015arXiv

Misfit stabilized embedded nanoparticles in metallic alloys

Nanoscale inhomogeneities are typical for numerous metallic alloys and crucially important for their practical applications. At the same time, stabilization mechanisms of such a state are poorly understood. We present a general overview of the problem, together with a more detailed discussion of the prototype example, namely, Guinier-Preston zones in Al-based alloys. It is shown that coherent strain due to a misfit between inclusion and host crystal lattices plays a decisive role in the emergence of the inhomogeneous state. We suggest a model explaining formation of ultrathin plates (with the thickness of a few lattice constants) typical for Al-Cu alloys. Discreteness of the array of misfit dislocations and long-ranged elastic interactions between them are the key ingredients of the model. This opens a way to a general understanding of the nature of (meta)stable embedded nanoparticles in practically important systems.

preprint2015arXiv

Modeling Klein tunneling and caustics of electron waves in graphene

We employ the tight-binding propagation method to study Klein tunneling and quantum interference in large graphene systems. With this efficient numerical scheme, we model the propagation of a wave packet through a potential barrier and determine the tunneling probability for different incidence angles. We consider both sharp and smooth potential barriers in n-p-n and n-n' junctions and find good agreement with analytical and semiclassical predictions. When we go outside the Dirac regime, we observe that sharp n-p junctions no longer show Klein tunneling because of intervalley scattering. However, this effect can be suppressed by considering a smooth potential. Klein tunneling holds for potentials changing on the scale much larger than the interatomic distance. When the energies of both the electrons and holes are above the Van Hove singularity, we observe total reflection for both sharp and smooth potential barriers. Furthermore, we consider caustic formation by a two-dimensional Gaussian potential. For sufficiently broad potentials we find a good agreement between the simulated wave density and the classical electron trajectories.

preprint2015arXiv

Moir{é} patterns as a probe of interplanar interactions: graphene on h-BN

By atomistic modeling of moir{é} patterns of graphene on a substrate with a small lattice mismatch, we find qualitatively different strain distributions for small and large misorientation angles, corresponding to the commensurate-incommensurate transition recently observed in graphene on hexagonal BN. We find that the ratio of C-N and C-B interactions is the main parameter determining the different bond lengths in the center and edges of the moir{é} pattern. Agreement with experimental data is obtained only by assuming that the C-B interactions are at least twice weaker than the C-N interactions. The correspondence between the strain distribution in the nanoscale moir{é} pattern and the potential energy surface at the atomic scale found in our calculations, makes the moir{é} pattern a tool to study details of dispersive forces in van der Waals heterostructures.

preprint2015arXiv

Phonon-assisted resonant tunneling of electrons in graphene-boron nitride transistors

We observe a series of sharp resonant features in the differential conductance of graphene-hexagonal boron nitride-graphene tunnel transistors over a wide range of bias voltages between $\sim$10 and 200 mV. We attribute them to electron tunneling assisted by the emission of phonons of well-defined energy. The bias voltages at which they occur are insensitive to the applied gate voltage and hence independent of the carrier densities in the graphene electrodes, so plasmonic effects can be ruled out. The phonon energies corresponding to the resonances are compared with the lattice dispersion curves of graphene-boron nitride heterostructures and are close to peaks in the single phonon density of states.

preprint2015arXiv

Quantum theory as a description of robust experiments: Application to Stern-Gerlach and Einstein-Podolsky-Rosen-Bohm experiments

We propose and develop the thesis that the quantum theoretical description of experiments emerges from the desire to organize experimental data such that the description of the system under scrutiny and the one used to acquire the data are separated as much as possible. Application to the Stern-Gerlach and Einstein-Podolsky-Rosen-Bohm experiments are shown to support this thesis. General principles of logical inference which have been shown to lead to the Schrödinger and Pauli equation and the probabilistic descriptions of the Stern-Gerlach and Einstein-Podolsky-Rosen-Bohm experiments, are used to demonstrate that the condition for the separation procedure to yield the quantum theoretical description is intimately related to the assumptions that the observed events are independent and that the data generated by these experiments is robust with respect to small changes of the conditions under which the experiment is carried out.

preprint2015arXiv

Racah materials: role of atomic multiplets in intermediate valence systems

We address the long-standing mystery of the nonmagnetic insulating state of the intermediate valence compound SmB$_6$. Within a combination of the local density approximation (LDA) and an exact diagonalization (ED) of an effective discrete Anderson impurity model, the intermediate valence ground state with the $f$-shell occupation $\langle n_{4f} \rangle=5.6$ is found for the Sm atom in SmB$_6$. This ground state is a singlet, and the first excited triplet state $\sim 3$ meV higher in the energy. SmB$_6$ is a narrow band insulator already in LDA, with the direct band gap of $\sim 10$ meV. The electron correlations increase the band gap which now becomes indirect. Thus, the many-body effects are relevant to form the indirect band gap, crucial for the idea of ``topological Kondo insulator" in SmB$_6$. Also, an actinide analog PuB$_6$ is considered, and the intermediate valence singlet ground state is found for the Pu atom. We propose that [Sm,Pu]B$_6$ belong to a new class of the intermediate valence materials with the multi-orbital ``Kondo-like" singlet ground-state. Crucial role of complex spin-orbital $f^n$-$f^{n+1}$ multiplet structure differently hybridized with ligand states in such Racah materials is discussed.

preprint2015arXiv

Relaxation of moiré patterns for slightly misaligned identical lattices: graphene on graphite

We study the effect of atomic relaxation on the structure of moiré patterns in twisted graphene on graphite and double layer graphene by large scale atomistic simulations. The reconstructed structure can be described as a superlattice of `hot spots' with vortex-like behaviour of in-plane atomic displacements and increasing out-of-plane displacements with decreasing angle. These lattice distortions affect both scalar and vector potential and the resulting electronic properties. At low misorientation angles (<$\sim$1$^\circ$) the optimized structures deviate drastically from the sinusoidal modulation which is often assumed in calculations of the electronic properties. The proposed structure might be verified by scanning probe microscopy measurements.

preprint2015arXiv

Scaling behavior and strain dependence of in-plane elastic properties of graphene

We show by atomistic simulations that, in the thermodynamic limit, the in-plane elastic moduli of graphene at finite temperature vanish with system size $ L $ as a power law $ ~ L^{-η_u} $ with $ η_u \simeq 0.325 $, in agreement with the membrane theory. Our simulations clearly reveal the size and strain dependence of graphene's elastic moduli, allowing comparison to experimental data. Although the recently measured difference of a factor 2 between the asymptotic value of the Young modulus for tensilely strained systems and the value from {\it ab initio} calculations remains unsolved, our results do explain the experimentally observed increase of more than a factor 2 for a tensile strain of only a few permille. We also discuss the scaling of the Poisson ratio, for which our simulations disagree with the predictions of the self-consistent screening approximation.

preprint2015arXiv

Toward a realistic description of multilayer black phosphorus: from $GW$ approximation to large-scale tight-binding simulations

We provide a tight-binding model parametrization for black phosphorus (BP) with an arbitrary number of layers. The model is derived from partially self-consistent $GW_0$ approach, where the screened Coulomb interaction $W_0$ is calculated within the random phase approximation on the basis of density functional theory. We thoroughly validate the model by performing a series of benchmark calculations, and determine the limits of its applicability. The application of the model to the calculations of electronic and optical properties of multilayer BP demonstrates good quantitative agreement with \emph{ab initio} results in a wide energy range. We also show that the proposed model can be easily extended for the case of external fields, yielding the results consistent with those obtained from first principles. The model is expected to be suitable for a variety of realistic problems related to the electronic properties of multilayer BP including different kinds of disorder, external fields, and many-body effects.

preprint2015arXiv

Transport and Optical Properties of Single- and Bilayer Black Phosphorus with Defects

We study the electronic and optical properties of single- and bilayer black phosphorus with shortand long-range defects by using the tight-binding propagation method. Both types of defect states are localized and induce a strong scattering of conduction states reducing significantly the charge carrier mobility. In contrast to properties of pristine samples, the anisotropy of defect-induced optical excitations is suppressed due to the isotropic nature of the defects. We also investigate the Landau level spectrum and magneto-optical conductivity, and find that the discrete Landau levels are sublinearly dependent on the magnetic field and energy level index, even at low defect concentrations.

preprint2014arXiv

Commensurate-incommensurate transition for graphene on hexagonal boron nitride

When a crystal is subjected to a periodic potential, under certain circumstances (such as when the period of the potential is close to the crystal periodicity; the potential is strong enough, etc.) it might adjust itself to follow the periodicity of the potential, resulting in a, so called, commensurate state. Such commensurate-incommensurate transitions are ubiquitous phenomena in many areas of condensed matter physics: from magnetism and dislocations in crystals, to vortices in superconductors, and atomic layers adsorbed on a crystalline surface. Of particular interest might be the properties of topological defects between the two commensurate phases: solitons, domain walls, and dislocation walls. Here we report a commensurate-incommensurate transition for graphene on top of hexagonal boron nitride (hBN). Depending on the rotational angle between the two hexagonal lattices, graphene can either stretch to adjust to a slightly different hBN periodicity (the commensurate state found for small rotational angles) or exhibit little adjustment (the incommensurate state). In the commensurate state, areas with matching lattice constants are separated by domain walls that accumulate the resulting strain. Such soliton-like objects present significant fundamental interest, and their presence might explain recent observations when the electronic, optical, Raman and other properties of graphene-hBN heterostructures have been notably altered.

preprint2014arXiv

Effect of Point Defects on the Optical and Transport Properties of MoS2 and WS2

Imperfections in the crystal structure, such as point defects, can strongly modify the optical and transport properties of materials. Here, we study the effect of point defects on the optical and DC conductivities of single layers of semiconducting transition metal dichalcogenides with the form $M$S$_2$, where $M$=Mo or W. The electronic structure is considered within a six bands tight-binding model, which accounts for the relevant combination of $d$ orbitals of the metal $M$ and $p$ orbitals of the chalcogen $S$. We use the Kubo formula for the calculation of the conductivity in samples with different distributions of disorder. We find that $M$ and/or S defects create mid-gap states that localize charge carriers around the defects and which modify the optical and transport properties of the material, in agreement with recent experiments. Furthermore, our results indicate a much higher mobility for $p$-doped WS$_2$ in comparison to MoS$_2$.

preprint2014arXiv

Electron correlations in Mn$_x$Ga$_{1-x}$As as seen by resonant electron spectroscopy and dynamical mean field theory

After two decades from the discovery of ferromagnetism in Mn-doped GaAs, its origin is still debated, and many doubts are related to the electronic structure. Here we report an experimental and theoretical study of the valence electron spectrum of Mn-doped GaAs. The experimental data are obtained through the differences between off- and on-resonance photo-emission data. The theoretical spectrum is calculated by means of a combination of density-functional theory in the local density approximation and dynamical mean-field theory (LDA+DMFT), using exact diagonalisation as impurity solver. Theory is found to accurately reproduce measured data, and illustrates the importance of correlation effects. Our results demonstrate that the Mn states extend over a broad range of energy, including the top of the valence band, and that no impurity band splits off from the valence band edge, while the induced holes seem located primarily around the Mn impurity.

preprint2014arXiv

Electronic correlations determine the phase stability of iron up to the melting temperature

We present theoretical results on the high-temperature phase stability and phonon spectra of paramagnetic bcc iron which explicitly take into account many-body effects. Several peculiarities, including a pronounced softening of the [110] transverse (T$_1$) mode and a dynamical instability of the bcc lattice in harmonic approximation are identified. We relate these features to the $α$-to-$γ$ and $γ$-to-$δ$ phase transformations in iron. The high-temperature bcc phase is found to be highly anharmonic and appears to be stabilized by the lattice entropy.

preprint2014arXiv

First-principles modelling of magnetic excitations in Mn12

We have developed a fully microscopic theory of magnetic properties of the prototype molecular magnet Mn12. First, the intra-molecular magnetic properties have been studied by means of first-principles density functional-based methods, with local correlation effects being taken into account within the local density approximation plus U (LDA+U) approach. Using the magnetic force theorem, we have calculated the interatomic isotropic and anisotropic exchange interactions and full tensors of single-ion anisotropy for each Mn ion. Dzyaloshinskii-Moriya (DM) interaction parameters turned out to be unusually large, reflecting a low symmetry of magnetic pairs in molecules, in comparison with bulk crystals. Based on these results we predict a distortion of ferrimagnetic ordering due to DM interactions. Further, we use an exact diagonalization approach allowing to work with as large Hilbert space dimension as 10^8 without any particular symmetry (the case of the constructed magnetic model). Based on the computational results for the excitation spectrum, we propose a distinct interpretation of the experimental inelastic neutron scattering spectra.

preprint2014arXiv

Inverse magneto-refraction as a mechanism for laser modification of spin-spin exchange parameters and subsequent terahertz emission from iron oxides

Ultrafast non-thermal manipulation of magnetization by light relies on either indirect coupling of the electric field component of the light with spins via spin-orbit interaction or direct coupling between the magnetic field component and spins. Here we propose a novel scenario for coupling between the electric field of light and spins via optical modification of the exchange interaction, one of the strongest quantum effects, the strength of which can reach 1000 Tesla. We demonstrate that this isotropic opto-magnetic effect, which can be called the inverse magneto-refraction, is allowed in a material of any symmetry. Its existence is corroborated by the experimental observation of THz emission by magnetic-dipole active spin resonances optically excited in a broad class of iron oxides with a canted spin configuration. From its strength we estimate that a sub-picosecond laser pulse with a moderate fluence of ~ 1 mJ/cm^2 acts as a pulsed effective magnetic field of 0.01 Tesla, arising from the optically perturbed balance between the exchange parameters. Our findings are supported by a low-energy theory for the microscopic magnetic interactions between non-equilibrium electrons subjected to an optical field which suggests a possibility to modify the exchange interactions by light over 1 %.

preprint2014arXiv

Metal-Insulator Transition in Graphene on Boron Nitride

Electrons in graphene aligned with hexagonal boron nitride are modelled by Dirac fermions in a correlated random-mass landscape subject to a scalar- and vector-potential disorder. We find that the system is insulating in the commensurate phase since the average mass deviates from zero. At the transition the mean mass is vanishing and electronic conduction in a finite sample can be described by a critical percolation along zero-mass lines. In this case graphene at the Dirac point is in a critical state with the conductivity $\sqrt{3}e^2/h$. In the incommensurate phase the system behaves as a symplectic metal.

preprint2014arXiv

Motion of domain walls and the dynamics of kinks in the magnetic Peierls potential

We study the dynamics of magnetic domain walls in the Peierls potential due to the discreteness of the crystal lattice. The propagation of a narrow domain wall (comparable to the lattice parameter) under the effect of a magnetic field proceeds through the formation of kinks in its profile. We predict that, despite the discreteness of the system, such kinks can behave like sine-Gordon solitons in thin films of materials such as yttrium iron garnets, and we derive general conditions for other materials. In our simulations we also observe long-lived breathers. We provide analytical expressions for the effective mass and limiting velocity of the kink in excellent agreement with our numerical results.

preprint2014arXiv

Plasmons in strongly correlated systems: spectral weight transfer and renormalized dispersion

We study the charge-density dynamics within the two-dimensional extended Hubbard model in the presence of long-range Coulomb interaction across the metal-insulator transition point. To take into account strong correlations we start from self-consistent extended dynamical mean-field theory and include non-local dynamical vertex corrections through a ladder approximation to the polarization operator. This is necessary to fulfill charge conservation and to describe plasmons in the correlated state. The calculated plasmon spectra are qualitatively different from those in the random-phase approximation: they exhibit a spectral density transfer and a renormalized dispersion with enhanced deviation from the canonical $\sqrt{q}$-behavior. Both features are reminiscent of interaction induced changes found in single-electron spectra of strongly correlated systems.

preprint2014arXiv

Probing of valley polarization in graphene via optical second-harmonic generation

Valley polarization in graphene breaks inversion symmetry and therefore leads to second-harmonic generation. We present a complete theory of this effect within a single-particle approximation. It is shown that this may be a sensitive tool to measure the valley polarization created, e.g., by polarized light and, thus, can be used for a development of ultrafast valleytronics in graphene.

preprint2014arXiv

Proton transport through one atom thick crystals

Graphene is impermeable to all gases and liquids, and even such a small atom as hydrogen is not expected to penetrate through graphene's dense electronic cloud within billions of years. Here we show that monolayers of graphene and hexagonal boron nitride (hBN) are unexpectedly permeable to thermal protons, hydrogen ions under ambient conditions. As a reference, no proton transport could be detected for a monolayer of molybdenum disulfide, bilayer graphene or multilayer hBN. At room temperature, monolayer hBN exhibits the highest proton conductivity with a low activation energy of about 0.3 eV but graphene becomes a better conductor at elevated temperatures such that its resistivity to proton flow is estimated to fall below 10^-3 Ohm per cm2 above 250 C. The proton barriers can be further reduced by decorating monolayers with catalytic nanoparticles. These atomically thin proton conductors could be of interest for many hydrogen-based technologies.

preprint2014arXiv

Quantum theory as the most robust description of reproducible experiments

It is shown that the basic equations of quantum theory can be obtained from a straightforward application of logical inference to experiments for which there is uncertainty about individual events and for which the frequencies of the observed events are robust with respect to small changes in the conditions under which the experiments are carried out.

preprint2014arXiv

Reply to 'Comment on "Thermodynamics of quantum crystalline membranes"'

In this note, we reply to the comment made by E.I.Kats and V.V.Lebedev [arXiv:1407.4298] on our recent work "Thermodynamics of quantum crystalline membranes" [Phys. Rev. B 89, 224307 (2014)]. Kats and Lebedev question the validity of the calculation presented in our work, in particular on the use of a Debye momentum as a ultra-violet regulator for the theory. We address and counter argue the criticisms made by Kats and Lebedev to our work.

preprint2014arXiv

Role of magnetic degrees of freedom in a scenario of phase transformations in steel

The diversity of mesostructures formed in steel at cooling from a high-temperature austenite ("gamma") phase is determined by the interplay of shear reconstructions of crystal lattice and diffusion of carbon. Combining first-principles calculations with large-scale phase-field simulations we demonstrate a decisive role of magnetic degrees of freedom in the formation of energy relief along the Bain path of "gamma"-"alpha" transformation and, thus, in this interplay. We show that there is the main factor, namely, the magnetic state of iron and its evolution with temperature which controls the change in character of the transformation. Based on the computational results we propose a simple model which reproduces, in good agreement with experiment, the most important curves of the phase transformation in Fe-C, namely, the lines relevant to a start of ferrite, bainite, and martensite transformations. Phase-field simulations within the model describe qualitatively typical patterns at these transformations.

preprint2014arXiv

Spin relaxation related to the edge scattering in graphene

We discuss the role of spin-flip scattering of electrons from the magnetized edges in graphene nanoribbons. The spin-flip scattering is associated with strong fluctuations of the magnetic moments at the edge. Using the Boltzmann equation approach, which is valid for not too narrow nanoribbons, we calculate the spin relaxation time in the case of Berry-Mondragon and zigzag graphene edges. We also consider the case of ballistic nanoribbons characterized by very long momentum relaxation time in the bulk, when the main source of momentum and spin relaxation is the spin-dependent scattering at the edges. We found that in the case of zigzag edges, an anomalous spin diffusion is possible, which is related to very weak spin-flip scattering of electrons gliding along the nanoribbon edge.

preprint2014arXiv

Thermodynamics of quantum crystalline membranes

We investigate the thermodynamic properties and the lattice stability of two-dimensional crystalline membranes, such as graphene and related compounds, in the low temperature quantum regime $T\rightarrow0$. A key role is played by the anharmonic coupling between in-plane and out-of plane lattice modes that, in the quantum limit, has very different consequences than in the classical regime. The role of retardation, namely of the frequency dependence, in the effective anharmonic interactions turns out to be crucial in the quantum regime. We identify a crossover temperature, $T^{*}$, between classical and quantum regimes, which is $\sim 70 - 90$ K for graphene. Below $T^{*}$, the heat capacity and thermal expansion coefficient decrease as power laws with decreasing temperature, tending to zero for $T\rightarrow0$ as required by the third law of thermodynamics.

preprint2014arXiv

Topological matter: graphene and superfluid 3He

Physics of graphene and physics of superfluid phases of 3He have many common features. Both systems are topological materials where quasiparticles behave as relativistic massless (Weyl, Majorana or Dirac) fermions. We formulate the points where these features are overlapping. This will allow us to use graphene for study the properties of superfluid 3He, to use superfluid 3He for study the properties of graphene, and to use the combination to study the physics of topological quantum vacuum. We suggest also some particular experiments with superfluid 3He using graphene as an atomically thin membrane impenetrable for He atoms but allowing for spin, momentum and energy transfer.

preprint2014arXiv

Twisted magnetic patterns: Exploring the Dzyaloshinskii--Moriya vector

Magnetism - the spontaneous alignment of atomic moments in a material - is driven by quantum-mechanical `exchange' interactions which operate over atomic distances as a result of the fundamental symmetry of electrons. Currently, one of the most active fields of condensed matter physics involves the study of magnetic interactions that cause, or are caused by a twisting of nearby atoms. This can lead to the magnetoelectric effect that couples electric and magnetic properties, and is predicted to play a prominent role in future technology. Here, we discuss the complex relativistic interplay between magnetism and atomic crystal structure in a class of materials called `weak ferromagnets'. The sign of the underpinning Dzyaloshinskii--Moriya interaction has been determined for the first time, by using synchrotron radiation to study iron borate (FeBO3). We present a novel experimental technique based on interference between two x-ray scattering processes (one acts as a reference wave) which we combine with a second unusual approach of turning the atomic antiferromagnetic motif with a small magnetic field. We show that the experimental results provide a clear validation of state-of-the-art theoretical calculations. These experimental and theoretical approaches open up new possibilities for exploring, modelling and exploiting novel magnetic and magnetoelectric materials.

preprint2013arXiv

Accurate determination of electron-hole asymmetry and next-nearest neighbor hopping in graphene

The next-nearest neighbor hopping term t' determines a magnitude and, hence, importance of several phenomena in graphene, which include self-doping due to broken bonds and the Klein tunneling that in the presence of t' is no longer perfect. Theoretical estimates for t' vary widely whereas a few existing measurements by using polarization resolved magneto-spectroscopy have found surprisingly large t', close or even exceeding highest theoretical values. Here we report dedicated measurements of the density of states in graphene by using high-quality capacitance devices. The density of states exhibits a pronounced electron-hole asymmetry that increases linearly with energy. This behavior yields t' approx -0.30 eV +-15%, in agreement with the high end of theory estimates. We discuss the role of electron-electron interactions in determining t' and overview phenomena which can be influenced by such a large value of t'.

preprint2013arXiv

Density functional theory analysis of flexural modes, elastic constants, and corrugations in strained graphene

Ab initio density functional theory has been used to analyze flexural modes, elastic constants, and atomic corrugations on single and bi-layer graphene. Frequencies of flexural modes are sensitive to compressive stress; its variation under stress can be related to the anomalous thermal expansion via a simple model based in classical Elasticity Theory [Phys. Rev. B 86, 144103]. Under compression, flexural modes are responsible for a long wavelength rippling with a large amplitude and a marked anharmonic behavior. This is compared with corrugations created by thermal fluctuations and the adsorption of a light impurity (hydrogen). Typical values for the later are in the sub-Angstrom regime, while maximum corrugations associated to bending modes quickly increase up to a few Angstroms under a compressive stress, due to the intrinsic instability of flexural modes

preprint2013arXiv

Doping Mechanisms in Graphene-MoS2 Hybrids

We present a joint theoretical and experimental investigation of charge doping and electronic potential landscapes in hybrid structures composed of graphene and semiconducting single layer MoS2. From first-principles simulations we find electron doping of graphene due to the presence of rhenium impurities in MoS2. Furthermore, we show that MoS2 edges give rise to charge reordering and a potential shift in graphene, which can be controlled through external gate voltages. The interplay of edge and impurity effects allows the use of the graphene-MoS2 hybrid as a photodetector. Spatially resolved photocurrent signals can be used to resolve potential gradients and local doping levels in the sample.

preprint2013arXiv

Dual origin of defect magnetism in graphene and its reversible switching by molecular doping

A possibility to control magnetic properties by using electric fields is one of the most desirable characteristics for spintronics applications. Finding a suitable material remains an elusive goal, with only a few candidates found so far. Graphene is one of them and offers a hope due to its weak spin-orbit interaction, the ability to control electronic properties by the electric field effect and the possibility to introduce paramagnetic centres such as vacancies and adatoms. Here we show that adatoms magnetism in graphene is itinerant and can be controlled by doping, so that magnetic moments can be switched on and off. The much-discussed vacancy magnetism is found to have a dual origin, with two approximately equal contributions: one coming from the same itinerant magnetism and the other due to broken bonds. Our work suggests that magnetic response of graphene can be controlled by the field effect, similar to its transport and optical properties, and that spin diffusion length can be significantly enhanced above a certain carrier density.

preprint2013arXiv

Edge scattering of electrons in graphene

We discuss the contribution of edge scattering to the conductance of graphene nanoribbons and nanoflakes. Using different possible types of the boundary conditions for the electron wave function at the edge, we found dependences of the momentum relaxation time and conductance on the geometric sizes and on the carrier density. We also consider the case of ballistic nanoribbon and nanodisc, for which the edge scaterring is the main mechanism of momentum relaxation.

preprint2013arXiv

Effect of magnetism on kinetics of "gamma"-"alpha" transformation and pattern formation in iron

Kinetics of polymorphous "gamma"-"alfa" transformation in Fe is studied numerically within a model taking into account both lattice and magnetic degrees of freedom, based on first-principle calculations of the total energy for different magnetic states. It is shown that magnetoelastic phenomena, namely, a strong sensitivity of the potential relief along the Bain deformation path to the magnetic state, are crucial for the picture of the transformation. With the temperature increase, a scenario of the phase transformation evolves from a homogeneous lattice instability at T<Ms (Ms is the temperature of the beginning of the martensitic transformation) to the growth and nucleation of embryos of the new phase at T>Ms. In the latter case, a stage of formation of a tweed-like structure occurs, with a strong short-range order and slow interphase fluctuations.

preprint2013arXiv

Effects of structural and chemical disorders on the visible/UV spectra of carbonaceous interstellar grains

The recent spectacular progress in the experimental and theoretical understanding of graphene, the basic constituent of graphite, is applied here to compute, from first principles, the UV extinction of nano-particles made of stacks of graphene layers. The theory also covers cases where graphene is affected by structural, chemical or orientation disorder, each disorder type being quantitatively defined by a single parameter. The extinction bumps carried by such model materials are found to have positions and widths falling in the same range as the known astronomical 2175 Åfeatures: as the disorder parameter increases, the bump width increases from 0.85 to 2.5 $μ$m$^{-1}$, while its peak position shifts from 4.65 to 4.75 $μ$m$^{-1}$. Moderate degrees of disorder are enough to cover the range of widths of the vast majority of observed bumps (0.75 to 1.3 $μ$m$^{-1}$). Higher degrees account for outliers, also observed in the sky. The introduction of structural or chemical disorder amounts to changing the initial $sp^{2}$ bondings into $sp^{3}$ or $sp^{1}$, so the optical properties of the model material become similar to those of the more or less amorphous carbon-rich materials studied in the laboratory: a-C, a-C:H, HAC, ACH, coals etc. The present treatment thus bridges gaps between physically different model materials.

preprint2013arXiv

Euler - Heisenberg effective action and magnetoelectric effect in multilayer graphene

The low energy effective field model for the multilayer graphene (at ABC stacking) is considered. We calculate the effective action in the presence of constant external magnetic field $B$ (normal to the graphene sheet). We also calculate the first two corrections to this effective action caused by the in-plane electric field $E$ at $E/B \ll 1$ and discuss the magnetoelectric effect. In addition, we calculate the imaginary part of the effective action in the presence of constant electric field $E$ and the lowest order correction to it due to the magnetic field ($B/E \ll 1$).

preprint2013arXiv

Exchange interactions and frustrated magnetism in single-side hydrogenated and fluorinated graphene

Magnetism in single-side hydrogenated (C$_2$H) and fluorinated (C$_2$F) graphene is analyzed in terms of the Heisenberg model with parameters determined from first principles. We predict a frustrated ground state for both systems, which means the instability of collinear spin structures and sheds light on the absence of a conventional magnetic ordering in defective graphene demonstrated in recent experiments. Moreover, our findings suggest a highly correlated magnetic behavior at low temperatures offering the possibility of a spin-liquid state.

preprint2013arXiv

Ferromagnetic two-dimensional crystals: Single layers of K2CuF4

The successful isolation of graphene ten years ago has evoked a rapidly growing scientific interest in the nature of two-dimensional (2D) crystals. A number of different 2D crystals has been produced since then, with properties ranging from superconductivity to insulating behavior. Here, we predict the possibility for realizing ferromagnetic 2D crystals by exfoliating atomically thin films of K2CuF4. From a first-principles theoretical analysis, we find that single layers of K2CuF4 form exactly 2D Kosterlitz-Thouless systems. The 2D crystal can form a free-standing membrane, and exhibits an experimentally accessible transition temperature and robust magnetic moments of 1 Bohr magneton per formula unit. 2D K2CuF4 unites ferromagnetic and insulating properties and is a demonstration of principles for nanoelectronics such as novel 2D-based heterostructures.

preprint2013arXiv

Many body renormalization of the minimal conductivity in graphene

The conductance of ballistic graphene at the neutrality point is due to coherent electron tunneling between the leads, the so called pseudodiffusive regime. The conductance scales as function of the sample dimensions in the same way as in a diffusive metal, despite the difference in the physical mechanisms involved. The electron-electron interaction modifies this regime, and plays a role similar to that of the environment in macroscopic quantum phenomena. We show that interactions, and the presence of external gates, change substantially the transport properties, and can lead to a diverging resistivity at the neutrality point.

preprint2013arXiv

Monte-Carlo study of the semimetal-insulator phase transition in monolayer graphene with realistic inter-electron interaction potential

We report on the results of the first-principle numerical study of spontaneous breaking of chiral (sublattice) symmetry in suspended monolayer graphene due to electrostatic interaction, which takes into account the screening of Coulomb potential by electrons on $σ$-orbitals. In contrast to the results of previous numerical simulations with unscreened potential, we find that suspended graphene is in the conducting phase with unbroken chiral symmetry. This finding is in agreement with recent experimental results by the Manchester group \cite{Elias et al, 2011}. Further, by artificially increasing the interaction strength we demonstrate that suspended graphene is quite close to the phase transition associated with spontaneous chiral symmetry breaking, which suggests that fluctuations of chirality and nonperturbative effects might still be quite important.

preprint2013arXiv

Non-local correlation effects and metal-insulator transition in the s-d exchange model

The metal-insulator transition within the s-d exchange model is studied by the Dynamical Mean Field Theory and Dual Fermion approaches. The latter takes into account nonlocal correlation effects which are shown to be essential. In particular, the critical values of the s-d exchange coupling constant for these two methods turn out to be different by a factor of more than two (for the case of square lattice and spin 1/2 localized electrons). The calculations were performed using a Continuous-Time Quantum Monte Carlo method. The difference of the quantum spin-1/2 case and the classical-spin case is discussed. For the quantum case the sign of s-d exchange coupling constant is relevant which demonstrates an importance of the Kondo effect.

preprint2013arXiv

Optimal Hubbard models for materials with nonlocal Coulomb interactions: graphene, silicene and benzene

To understand how nonlocal Coulomb interactions affect the phase diagram of correlated electron materials, we report on a method to approximate a correlated lattice model with nonlocal interactions by an effective Hubbard model with on-site interactions U* only. The effective model is defined by the Peierls-Feynman-Bogoliubov variational principle. We find that the local part of the interaction U is reduced according to U*=U-V', where V' is a weighted average of nonlocal interactions. For graphene, silicene and benzene we show that the nonlocal Coulomb interaction can decrease the effective local interaction by more than a factor of 2 in a wide doping range.

preprint2013arXiv

Screening and Collective Modes in Disordered Graphene Antidot Lattices

The excitation spectrum and the collective modes of graphene antidot lattices (GALs) are studied in the context of a $π$-band tight-binding model. The dynamical polarizability and dielectric function are calculated within the random phase approximation. The effect of different kinds of disorder, such as geometric and chemical disorder, are included in our calculations. We highlight the main differences of GALs with respect to single-layer graphene (SLG). Our results show that, in addition to the well-understood bulk plasmon in doped samples, inter-band plasmons appear in GALs. We further show that the static screening properties of undoped and doped GALs are quantitatively different from SLG.

preprint2013arXiv

Semiclassical theory of potential scattering for massless Dirac fermions

In this paper we study scattering of two-dimensional massless Dirac fermions by a potential that depends on a single Cartesian variable. Depending on the energy of the incoming particle and its angle of incidence, there are three different regimes of scattering. To find the reflection and transmission coefficients in these regimes, we apply the Wentzel-Kramers-Brillouin (WKB), also called semiclassical, approximation. We use the method of comparison equations to extend our prediction to nearly normal incidence, where the conventional WKB method should be modified due to the degeneracy of turning points. We compare our results to numerical calculations and find good agreement.

preprint2013arXiv

Structure and magnetism of disordered carbon

Motivated by the observation of ferromagnetism in carbon foams, a massive search for (meta)stable disorder structures of elemental carbon is performed by a generate and test approach. We use the Density Functional based program SIESTA to optimize the structures and calculate the electronic spectra and spin densities. About 1% of the 24000 optimized structures presents magnetic moments, a necessary but not sufficient condition for intrinsic magnetic order. We analyze the results using elements of graph theory. Although the relation between structure and the occurrence of magnetic moments is not yet fully clarified, we give some minimal requirements for this possibility, such as the existence of three-fold coordinated atoms surrounded by four-fold coordinated atoms. We discuss in detail the most promising structures.

preprint2013arXiv

Unified character of correlation effects in unconventional Pu-based superconductors and δ-Pu

Electronic structure calculations combining the local-density approximation with an exact diagonalization of the Anderson impurity model show an intermediate 5f^5-5f^6-valence ground state and delocalization of the 5f^5 multiplet of the Pu atom 5f-shell in PuCoIn_5, PuCoGa_5, and δ-Pu. The 5f-local magnetic moment is compensated by a moment formed in the surrounding cloud of conduction electrons. For PuCoGa_5 and δ-Pu the compensation is complete and the Anderson impurity ground state is a singlet. For PuCoIn_5 the compensation is partial and the Pu ground state is magnetic. We suggest that the unconventional d-wave superconductivity is likely mediated by the 5f-states antiferromagnetic fluctuations in PuCoIn_5, and by valence fluctuations in PuCoGa_5.

preprint2013arXiv

Unruh effect in vacua with anisotropic scaling: Applications to multilayer graphene

We extend the calculation of the Unruh effect to the universality classes of quantum vacua obeying topologically protected invariance under anisotropic scaling ${\bf r} \rightarrow b {\bf r}$, $t \rightarrow b^z t$. Two situations are considered. The first one is related to the accelerated detector which detects the electron - hole pairs. The second one is related to the system in external electric field, when the electron - hole pairs are created due to the Schwinger process. As distinct from the Unruh effect in relativistic systems (where $z=1$) the calculated radiation is not thermal, but has properties of systems in the vicinity of quantum criticality. The vacuum obeying anisotropic scaling can be realized, in particular, in multilayer graphene with the rhombohedral stacking. Opportunities of the experimental realization of Unruh effect in this situation are discussed.

preprint2013arXiv

Zero modes in magnetic systems: general theory and an efficient computational scheme

The presence of topological defects in magnetic media often leads to normal modes with zero frequency (zero modes). Such modes are crucial for long-time behavior, describing, for example, the motion of a domain wall as a whole. Conventional numerical methods to calculate the spin-wave spectrum in magnetic media are either inefficient or they fail for systems with zero modes. We present a new efficient computational scheme that reduces the magnetic normal-mode problem to a generalized Hermitian eigenvalue problem also in the presence of zero modes. We apply our scheme to several examples, including two-dimensional domain walls and Skyrmions, and show how the effective masses that determine the dynamics can be calculated directly. These systems highlight the fundamental distinction between the two types of zero modes that can occur in spin systems, which we call special and inertial zero modes. Our method is suitable for both conservative and dissipative systems. For the latter case, we present a perturbative scheme to take into account damping, which can also be used to calculate dynamical susceptibilities.

preprint2012arXiv

Adsorption of cobalt on graphene: Electron correlation effects from a quantum chemical perspective

In this work, we investigate the adsorption of a single cobalt atom (Co) on graphene by means of the complete active space self-consistent field approach, additionally corrected by the second-order perturbation theory. The local structure of graphene is modeled by a planar hydrocarbon cluster (C$_{24}$H$_{12}$). Systematic treatment of the electron correlations and the possibility to study excited states allow us to reproduce the potential energy curves for different electronic configurations of Co. We find that upon approaching the surface, the ground-state configuration of Co undergoes several transitions, giving rise to two stable states. The first corresponds to the physisorption of the adatom in the high-spin $3d^74s^2$ ($S=3/2$) configuration, while the second results from the chemical bonding formed by strong orbital hybridization, leading to the low-spin $3d^9$ ($S=1/2$) state. Due to the instability of the $3d^9$ configuration, the adsorption energy of Co is small in both cases and does not exceed 0.35 eV. We analyze the obtained results in terms of a simple model Hamiltonian that involves Coulomb repulsion ($U$) and exchange coupling ($J$) parameters for the 3$d$ shell of Co, which we estimate from first-principles calculations. We show that while the exchange interaction remains constant upon adsorption ($\simeq1.1$ eV), the Coulomb repulsion significantly reduces for decreasing distances (from 5.3 to 2.6$\pm$0.2 eV). The screening of $U$ favors higher occupations of the 3$d$ shell and thus is largely responsible for the interconfigurational transitions of Co. Finally, we discuss the limitations of the approaches that are based on density functional theory with respect to transition metal atoms on graphene, and we conclude that a proper account of the electron correlations is crucial for the description of adsorption in such systems.

preprint2012arXiv

Bending modes, anharmonic effects and thermal expansion coefficient in single layer and multilayer graphene

We present a simple analytical approach to study anharmonic effects in single layer, bilayer, and multilayer graphene. The coupling between in plane and out of plane modes leads to negative Grüneisen coefficients and negative thermal expansion. The value of the thermal expansion coefficient depends on the coupling to the substrate. The bending rigidity in bilayer graphene shows a crossover between a long wavelength regime where its value is determined by the in plane elastic properties and a short wavelength regime where its value approaches twice that of a single layer.

preprint2012arXiv

Dangling bonds and magnetism of grain boundaries in graphene

Grain boundaries with dangling bonds (DBGB) in graphene are studied by atomistic Monte Carlo and molecular dynamics simulations in combination with density functional (SIESTA) calculations. The most stable configurations are selected and their structure is analyzed in terms of grain boundary dislocations. It is shown that the grain boundary dislocation with the core consisting of pentagon, octagon and heptagon (5-8-7 defect) is a typical structural element of DBGB with relatively low energies. Electron energy spectrum and magnetic properties of the obtained DBGB are studied by density functional calculations. It is shown that the 5-8-7 defect is magnetic and that its magnetic moment survives after hydrogenation. The effects of hydrogenation and of out of plane deformations on the magnetic properties of DBGB are studied.

preprint2012arXiv

Effect of impurities on growth and morphology of cementite nanowires

Effects of doping on morphology of iron carbide (cementite) nanowires have been explored by first principle electronic structure calculations. We examined the role of several realistic impurities (Si, Mn, V, P and S) in the formation energies of the cementite nanowires with different sizes and morphologies. It is shown that the presence of the impurities decreases the formation energy and can switch the preferable axis of the cementite nanowire growth. The conditions of the formation and decomposition of cementite nanowires in steels are also discussed.

preprint2012arXiv

Electronic Properties of Disordered Graphene Antidot Lattices

Regular nanoscale perforations in graphene (graphene antidot lattices, GAL) are known to lead to a gap in the energy spectrum, thereby paving a possible way towards many applications. This theoretical prediction relies on a perfect placement of identical perforations, a situation not likely to occur in the laboratory. Here, we present a systematic study of the effects of disorder in GALs. We consider both geometric and chemical disorder, and evaluate the density-of-states as well as the optical conductivity of disordered GALs. The theoretical method is based on an efficient algorithm for solving the time-dependent Schr{ö}dinger equation in a tight-binding representation of the graphene sheet [S. Yuan et al., Phys. Rev. B 82, 115448 (2010)], which allows us to consider GALs consisting of 6400 $\times$ 6400 carbon atoms. The central conclusion for all kinds of disorder is that the gaps found for pristine GALs do survive at a considerable amount of disorder, but disappear for very strong disorder. Geometric disorder is more detrimental to gap formation than chemical disorder. The optical conductivity shows a low-energy tail below the pristine GAL band gap due to disorder-introduced transitions.

preprint2012arXiv

Electronic, Magnetic and Transport Properties of Graphene Ribbons Terminated by Nanotubes

We study by density functional and large scale tight-binding transport calculations the electronic structure, magnetism and transport properties of the recently proposed graphene ribbons with edges rolled to form nanotubes. Edges with armchair nanotubes present magnetic moments localized either in the tube or the ribbon and metallic or half-metallic character, depending on the symmetry of the junction. These properties have potential for spin valve and spin filter devices with advantages over other proposed systems. Edges with zigzag nanotubes are either metallic or semiconducting without affecting the intrinsic mobility of the ribbon. By varying the type and size of the nanotubes and ribbons offers the possibility to tailor the magnetic and transport properties, making these systems very promising for applications.

preprint2012arXiv

Enhanced Screening in Chemically Functionalized Graphene

Resonant scatterers such as hydrogen adatoms can strongly enhance the low energy density of states in graphene. Here, we study the impact of these impurities on the electronic screening. We find a two-faced behavior: Kubo formula calculations reveal an increased dielectric function $\varepsilon$ upon creation of midgap states but no metallic divergence of the static $\varepsilon$ at small momentum transfer $q\to 0$. This bad metal behavior manifests also in the dynamic polarization function and can be directly measured by means of electron energy loss spectroscopy. A new length scale $l_c$ beyond which screening is suppressed emerges, which we identify with the Anderson localization length.

preprint2012arXiv

Graphene in periodic deformation fields: dielectric screening and plasmons

We consider the effect of periodic scalar and vector potentials generated by periodic deformations of the graphene crystal lattice, on the energy spectrum of electrons. The dependence of electron velocity near the Dirac point on the periodic perturbations of different types is discussed. We also investigated the effect of screening of the scalar potential by calculating the dielectric function as a function of the wave length of the periodic potential. This calculation shows that the periodic scalar field is strongly suppressed by the screening. Using the dependence of electron velocity on the periodic field we also studied the variation of the plasmon spectra in graphene. We found that the spectrum of plasmon excitations can be effectively controlled by the periodic strain field.

preprint2012arXiv

GraXe, graphene and xenon for neutrinoless double beta decay searches

We propose a new detector concept, GraXe (to be pronounced as grace), to search for neutrinoless double beta decay in Xe-136. GraXe combines a popular detection medium in rare-event searches, liquid xenon, with a new, background-free material, graphene. In our baseline design of GraXe, a sphere made of graphene-coated titanium mesh and filled with liquid xenon (LXe) enriched in the Xe-136 isotope is immersed in a large volume of natural LXe instrumented with photodetectors. Liquid xenon is an excellent scintillator, reasonably transparent to its own light. Graphene is transparent over a large frequency range, and impermeable to the xenon. Event position could be deduced from the light pattern detected in the photosensors. External backgrounds would be shielded by the buffer of natural LXe, leaving the ultra-radiopure internal volume virtually free of background. Industrial graphene can be manufactured at a competitive cost to produce the sphere. Enriching xenon in the isotope Xe-136 is easy and relatively cheap, and there is already near one ton of enriched xenon available in the world (currently being used by the EXO, KamLAND-Zen and NEXT experiments). All the cryogenic know-how is readily available from the numerous experiments using liquid xenon. An experiment using the GraXe concept appears realistic and affordable in a short time scale, and its physics potential is enormous.

preprint2012arXiv

Importance of correlation effects in hcp iron revealed by a pressure-induced electronic topological transition

We discover that hcp phases of Fe and Fe0.9Ni0.1 undergo an electronic topological transition at pressures of about 40 GPa. This topological change of the Fermi surface manifests itself through anomalous behavior of the Debye sound velocity, c/a lattice parameter ratio and Mössbauer center shift observed in our experiments. First-principles simulations within the dynamic mean field approach demonstrate that the transition is induced by many-electron effects. It is absent in one-electron calculations and represents a clear signature of correlation effects in hcp Fe.

preprint2012arXiv

Quantum electrodynamics with anisotropic scaling: Heisenberg-Euler action and Schwinger pair production in the bilayer graphene

We discuss quantum electrodynamics emerging in the vacua with anisotropic scaling. Systems with anisotropic scaling were suggested by Horava in relation to the quantum theory of gravity. In such vacua the space and time are not equivalent, and moreover they obey different scaling laws, called the anisotropic scaling. Such anisotropic scaling takes place for fermions in bilayer graphene, where if one neglects the trigonal warping effects the massless Dirac fermions have quadratic dispersion. This results in the anisotropic quantum electrodynamics, in which electric and magnetic fields obey different scaling laws. Here we discuss the Heisenberg-Euler action and Schwinger pair production in such anisotropic QED

preprint2012arXiv

Scattering by flexural phonons in suspended graphene under back gate induced strain

We have studied electron scattering by out-of-plane (flexural) phonon modes in doped suspended graphene and its effect on charge transport. In the free-standing case (absence of strain) the flexural branch shows a quadratic dispersion relation, which becomes linear at long wavelength when the sample is under tension due to the rotation symmetry breaking. In the non-strained case, scattering by flexural phonons is the main limitation to electron mobility. This picture changes drastically when strains above $\bar{u}=10^{-4} n(10^{12}\,\text{cm}^{-2})$ are considered. Here we study in particular the case of back gate induced strain, and apply our theoretical findings to recent experiments in suspended graphene.

preprint2012arXiv

Strong Coulomb drag and broken symmetry in double-layer graphene

Spatially separated electron systems remain strongly coupled by electron-electron interactions even when they cannot exchange particles, provided that the layer separation d is comparable to a characteristic distance l between charge carriers within layers. One of the consequences of this remote coupling is a phenomenon called Coulomb drag, in which an electric current passed through one of the layers causes frictional charge flow in the other layer. Previously, only the regime of weak (d>>l) to intermediate (d ~ l) coupling could be studied experimentally. Here we use graphene-BN heterostructures with d down to 1 nm to probe interlayer interactions and Coulomb drag in the limit d<<l where the two Dirac liquids effectively nest within the same plane, but still can be tuned and measured independently. The strongly interacting regime reveals many unpredicted features that are qualitatively different from those observed previously. In particular, although drag vanishes because of electron-hole symmetry when either layer is neutral, we find that drag is at its strongest when /both/ layers are neutral. Under this circumstance, drag is positive in zero magnetic field but changes its sign and rapidly grows in strength with field. The drag remains strong at room temperature. The broken electron-hole symmetry is attributed to mutual polarization of closely-spaced interacting layers. The two-fluid Dirac system offers a large parameter space for further investigation of strong interlayer interaction phenomena.

preprint2012arXiv

The Aharonov-Bohm effect for massless Dirac fermions and the spectral flow of Dirac type operators with classical boundary conditions

We compute, in topological terms, the spectral flow of an arbitrary family of self-adjoint Dirac type operators with classical (local) boundary conditions on a compact Riemannian manifold with boundary under the assumption that the initial and terminal operators of the family are conjugate by a bundle automorphism. This result is used to study conditions for the existence of nonzero spectral flow of a family of self-adjoint Dirac type operators with local boundary conditions in a two-dimensional domain with nontrivial topology. Possible physical realizations of nonzero spectral flow are discussed.

preprint2012arXiv

Theory of Coulomb drag for massless Dirac fermions

Coulomb drag between two unhybridized graphene sheets separated by a dielectric spacer has recently attracted considerable theoretical interest. We first review, for the sake of completeness, the main analytical results which have been obtained by other authors. We then illustrate pedagogically the minimal theory of Coulomb drag between two spatially-separated two-dimensional systems of massless Dirac fermions which are both away from the charge-neutrality point. This relies on second-order perturbation theory in the screened interlayer interaction and on Boltzmann transport theory. In this theoretical framework and in the low-temperature limit, we demonstrate that, to leading (i.e. quadratic) order in temperature, the drag transresistivity is completely insensitive to the precise intralayer momentum-relaxation mechanism (i.e. to the functional dependence of the scattering time on energy). We also provide analytical results for the low-temperature drag transresistivity for both cases of "thick" and "thin" spacers and for arbitrary values of the dielectric constants of the media surrounding the two Dirac-fermion layers. Finally, we present numerical results for the low-temperature drag transresistivity in the case in which one of the media surrounding the Dirac-fermion layers has a frequency-dependent dielectric constant. We conclude by suggesting an experiment that can potentially allow for the observation of departures from the canonical Fermi-liquid quadratic-in-temperature behavior of the transresistivity.

preprint2011arXiv

Adhesion and electronic structure of graphene on hexagonal boron nitride substrates

We investigate the adsorption of graphene sheets on h-BN substrates by means of first-principles calculations in the framework of adiabatic connection fluctuation-dissipation theory in the random phase approximation. We obtain adhesion energies for different crystallographic stacking configurations and show that the interlayer bonding is due to long-range van der Waals forces. The interplay of elastic and adhesion energies is shown to lead to stacking disorder and moiré structures. Band structure calculations reveal substrate induced mass terms in graphene which change their sign with the stacking configuration. The dispersion, absolute band gaps and the real space shape of the low energy electronic states in the moiré structures are discussed. We find that the absolute band gaps in the moiré structures are at least an order of magnitude smaller than the maximum local values of the mass term. Our results are in agreement with recent STM experiments.

preprint2011arXiv

Chiral tunneling in single and bilayer graphene

We review chiral (Klein) tunneling in single-layer and bilayer graphene and present its semiclassical theory, including the Berry phase and the Maslov index. Peculiarities of the chiral tunneling are naturally explained in terms of classical phase space. In a one-dimensional geometry we reduced the original Dirac equation, describing the dynamics of charge carriers in the single layer graphene, to an effective Schrödinger equation with a complex potential. This allowed us to study tunneling in details and obtain analytic formulas. Our predictions are compared with numerical results. We have also demonstrated that, for the case of asymmetric n-p-n junction in single layer graphene, there is total transmission for normal incidence only, side resonances are suppressed.

preprint2011arXiv

Coexistence of electron and hole transport in graphene

When sweeping the carrier concentration in monolayer graphene through the charge neutrality point, the experimentally measured Hall resistivity shows a smooth zero crossing. Using a two- component model of coexisting electrons and holes around the charge neutrality point, we unambiguously show that both types of carriers are simultaneously present. For high magnetic fields up to 30 T the electron and hole concentrations at the charge neutrality point increase with the degeneracy of the zero-energy Landau level which implies a quantum Hall metal state at ν=0 made up by both electrons and holes.

preprint2011arXiv

Coulomb drag in graphene single layers separated by a thin spacer

Motivated by very recent studies of Coulomb drag in grahene-BN-graphene system we develop a theory of Coulomb drag for the Fermi liquid regime, for the case when the ratio of spacer thickness $d$ to the Fermi wavelength of electrons is arbitrary. The concentration ($n$) and thickness dependence of the drag resistivity is changed from $n^{-3}d^{-4}$ for the thick spacer to $n^{-1}|\ln{(nd^2)}|$ for the thin one.

preprint2011arXiv

Disorder induced metallicity in amorphous graphene

We predict a transition to metallicity when a sufficient amount of disorder is induced in graphene. Calculations were performed by means of a first principles stochastic quench method. The resulting amorphous graphene can be seen as nanopatches of graphene that are connected by a network of disordered small and large carbon rings. The buckling is minimal and we believe that it is a result of averaging of counteracting random in-plane stress forces. The linear response conductance is obtained by a model theory as function of lattice distortions. Such metallic behaviour is a much desired property for functionalisation of graphene to realize a transparent conductor, e.g. suitable for touch-screen devices.

preprint2011arXiv

Dual boson approach to collective excitations in correlated fermionic systems

We develop a general theory of a boson decomposition for both local and non-local interactions in lattice fermion models which allows us to describe fermionic degrees of freedom and collective charge and spin excitations on equal footing. An efficient perturbation theory in the interaction of the fermionic and the bosonic degrees of freedom is constructed in so-called dual variables in the path-integral formalism. This theory takes into account all local correlations of fermions and collective bosonic modes and interpolates between itinerant and localized regimes of electrons in solids. The zero-order approximation of this theory corresponds to extended dynamical mean-field theory (EDMFT), a regular way to calculate nonlocal corrections to EDMFT is provided. It is shown that dual ladder summation gives a conserving approximation beyond EDMFT. The method is especially suitable for consideration of collective magnetic and charge excitations and allows to calculate their renormalization with respect to "bare" RPA-like characteristics. General expression for the plasmonic dispersion in correlated media is obtained. As an illustration it is shown that effective superexchange interactions in the half-filled Hubbard model can be derived within the dual-ladder approximation.

preprint2011arXiv

Enhancement of the Curie temperature in small particles of weak itinerant ferromagnets

Self consistent renormalization theory of itinerant ferromagnets is used to calculate the Curie temperature of clusters down to approximately 100 atoms in size. In these clusters the electrons responsible for the magnetic properties are assumed to be (weakly) itinerant. It is shown that the Curie temperature can be larger than in the bulk. The effect originates from the phenomenon of level repulsion in chaotic quantum systems, which suppresses spin fluctuations. Since the latter destroy the magnetic order the resulting Curie temperature increases, contrary to expectations of the naive Stoner picture. The calculations are done assuming that the energy levels of the cluster are described by the Gaussian Orthogonal Ensemble of random matrix theory.

preprint2011arXiv

Field-effect tunneling transistor based on vertical graphene heterostructures

We report a bipolar field effect tunneling transistor that exploits to advantage the low density of states in graphene and its one atomic layer thickness. Our proof-of-concept devices are graphene heterostructures with atomically thin boron nitride acting as a tunnel barrier. They exhibit room temperature switching ratios ~50, a value that can be enhanced further by optimizing the device structure. These devices have potential for high frequency operation and large scale integration.

preprint2011arXiv

Gaps tunable by electrostatic gates in strained graphene

We show that when the pseudomagnetic fields created by long wavelength deformations are appropriately coupled with a scalar electric potential, a significant energy gap can emerge due to the formation of a Haldane state. Ramifications of this physical effect are examined through the study of various strain geometries commonly seen in experiments, such as strain superlattices and wrinkled suspended graphene. Of particular technological importance, we consider setup where this gap can be tunable through electrostatic gates, allowing for the design of electronic devices not realizable with other materials.

preprint2011arXiv

Gauge fields at the surface of topological insulators

We study the emergence of gauge couplings in the surface states of topological insulators. We show that gauge fields arise when a three dimensional strong topological insulator is coated with an easy-plane ferromagnet with magnetization parallel to the surface. We analyze the modification induced by the gauge fields on the surface spectrum for some specific magnetic configurations.

preprint2011arXiv

Giant Nonlocality near the Dirac Point in Graphene

Transport measurements have been a powerful tool for uncovering new electronic phenomena in graphene. We report nonlocal measurements performed in the Hall bar geometry with voltage probes far away from the classical path of charge flow. We observe a large nonlocal response near the Dirac point in fields as low as 0.1T, which persists up to room temperature. The nonlocality is consistent with the long-range flavor currents induced by lifting of spin/valley degeneracy. The effect is expected to contribute strongly to all magnetotransport phenomena near the neutrality point.

preprint2011arXiv

Graphene adhesion on mica: Role of surface morphology

We investigate theoretically the adhesion and electronic properties of graphene on a muscovite mica surface using the density functional theory (DFT) with van der Waals (vdW) interactions taken into account (the vdW-DF approach). We found that irregularities in the local structure of cleaved mica surface provide different mechanisms for the mica-graphene binding. By assuming electroneutrality for both surfaces, the binding is mainly of vdW nature, barely exceeding thermal energy per carbon atom at room temperature. In contrast, if potassium atoms are non uniformly distributed on mica, the different regions of the surface give rise to $n$- or $p$-type doping of graphene. In turn, an additional interaction arises between the surfaces, significantly increasing the adhesion. For each case the electronic states of graphene remain unaltered by the adhesion. It is expected, however, that the Fermi level of graphene supported on realistic mica could be shifted relative to the Dirac point due to asymmetry in the charge doping. Obtained variations of the distance between graphene and mica for different regions of the surface are found to be consistent with recent atomic force microscopy experiments. A relative flatness of mica and the absence of interlayer covalent bonding in the mica-graphene system make this pair a promising candidate for practical use.

preprint2011arXiv

Interaction-Driven Spectrum Reconstruction in Bilayer Graphene

The nematic phase transition in electronic liquids, driven by Coulomb interactions, represents a new class of strongly correlated electronic ground states. We studied suspended samples of bilayer graphene, annealed so that it achieves very high quasiparticle mobilities. Bilayer graphene is a truly two-dimensional material with complex chiral electronic spectra and the high quality of our samples allowed us to observe strong spectrum reconstructions and electron topological transitions that can be attributed to a nematic phase transition and a decrease in rotational symmetry. These results are especially surprising because no interaction effects have been observed so far in bilayer graphene in the absence of an applied magnetic field.

preprint2011arXiv

Interfacial interactions between local defects in amorphous SiO$_2$ and supported graphene

We present a density functional study of graphene adhesion on a realistic SiO$_2$ surface taking into account van der Waals (vdW) interactions. The SiO$_2$ substrate is modeled at the local scale by using two main types of surface defects, typical for amorphous silica: the oxygen dangling bond and three-coordinated silicon. The results show that the nature of adhesion between graphene and its substrate is qualitatively dependent on the surface defect type. In particular, the interaction between graphene and silicon-terminated SiO$_2$ originates exclusively from the vdW interaction, whereas the oxygen-terminated surface provides additional ionic contribution to the binding arising from interfacial charge transfer ($p$-type doping of graphene). Strong doping contrast for the different surface terminations provides a mechanism for the charge inhomogeneity of graphene on amorphous SiO$_2$ observed in experiments. We found that independent of the considered surface morphologies, the typical electronic structure of graphene in the vicinity of the Dirac point remains unaltered in contact with the SiO$_2$ substrate, which points to the absence of the covalent interactions between graphene and amorphous silica. The case of hydrogen-passivated SiO$_2$ surfaces is also examined. In this situation, the binding with graphene is practically independent of the type of surface defects and arises, as expected, from the vdW interactions. Finally, the interface distances obtained are shown to be in good agreement with recent experimental studies.

preprint2011arXiv

Magnetoelastic coupling in iron

Exchange interactions in α- and γ-Fe are investigated within an ab-initio spin spiral approach. We have performed total energy calculations for different magnetic structures as a function of lattice distortions, related with various cell volumes and the Bain tetragonal deformations. The effective exchange parameters in γ-Fe are very sensitive to the lattice distortions, leading to the ferromagnetic ground state for the tetragonal deformation or increase of the volume cell. At the same time, the magnetic-structure-independent part of the total energy changes very slowly with the tetragonal deformations. The computational results demonstrate a strong mutual dependence of crystal and magnetic structures in Fe and explain the observable "anti-Invar" behavior of thermal expansion coefficient in γ-Fe.

preprint2011arXiv

Melting of graphene: from two to one dimension

The high temperature behaviour of graphene is studied by atomistic simulations based on an accurate interatomic potential for carbon. We find that clustering of Stone-Wales defects and formation of octagons are the first steps in the process of melting which proceeds via the formation of carbon chains. The molten state forms a three-dimensional network of entangled chains rather than a simple liquid. The melting temperature estimated from the two-dimensional Lindemann criterion and from extrapolation of our simulation for different heating rates is about 4900 K.

preprint2011arXiv

Quantum transport via evanescent waves in undoped graphene

Charge carriers in graphene are chiral quasiparticles ("massless Dirac fermions"). Graphene provides therefore an amazing opportunity to study subtle quantum relativistic effects in condensed matter experiment. Here I review a theory of one of these unusual features of graphene, a "pseudodiffusive" transport in the limit of zero charge carrier concentration, which is related to existence of zero-modes of the Dirac operator and to the Zitterbewegung of unltrarelativistic particles. A conformal mapping technique is a powerful mathematical tool to study these phenomena, as demonstrated here, using the Aharonov-Bohm effect in graphene rings with Corbino geometry as an example.

preprint2011arXiv

sp-Electron Magnetic Clusters with a Large Spin in Graphene

Motivated by recent experimental data (Sepioni, M. et al. Phys. Rev. Lett. 2010, 105, 207205), we have studied the possibility of forming magnetic clusters with spin S> 1/2 on graphene by adsorption of hydrogen atoms or hydroxyl groups. Migration of hydrogen atoms and hydroxyl groups on the surface of graphene during the delamination of HOPG led to the formation of seven-atom or seven-OH-group clusters with S=5/2 that were of a special interest. The coincidence of symmetry of the clusters with the graphene lattice strengthens the stability of the cluster. For (OH)7 clusters that were situated greater than 3 nm from one another, the reconstruction barrier to a nonmagnetic configuration was approximately 0.4 eV, whereas for H7 clusters, there was no barrier and the high-spin state was unstable. Stability of the high-spin clusters increased if they were formed on top of ripples. Exchange interactions between the clusters were studied and we have shown that the ferromagnetic state is improbable. The role of the chemical composition of the solvent used for the delamination of graphite is discussed.

preprint2011arXiv

Strength of effective Coulomb interactions in graphene and graphite

To obtain an effective many-body model of graphene and related materials from first principles we calculate the partially screened frequency dependent Coulomb interaction. In graphene, the effective on-site (Hubbard) interaction is U_00 = 9.3 eV in close vicinity to the critical value separating conducting graphene from an insulating phase emphasizing the importance of non-local Coulomb terms. The nearest-neighbor Coulomb interaction strength is computed to U_01 = 5.5 eV. In the long wavelength limit, we find the effective background dielectric constant of graphite to be ε= 2.5 in very good agreement with experiment.

preprint2011arXiv

Temperature dependent resistivity in bilayer graphene due to flexural phonons

We have studied electron scattering by out-of-plane (flexural) phonons in doped suspended bilayer graphene. We have found the bilayer membrane to follow the qualitative behavior of the monolayer cousin. In the bilayer, different electronic structure combine with different electron-phonon coupling to give the same parametric dependence in resistivity, and in particular the same temperature $T$ behavior. In parallel with the single layer, flexural phonons dominate the phonon contribution to resistivity in the absence of strain, where a density independent mobility is obtained. This contribution is strongly suppressed by tension, and in-plane phonons become the dominant contribution in strained samples. Among the quantitative differences an important one has been identified: room $T$ mobility in bilayer graphene is substantially higher than in monolayer. The origin of quantitative differences has been unveiled.

preprint2011arXiv

Transition metal ad-atoms on graphene: Influence of local Coulomb interactions on chemical bonding and magnetic moments

We address the interaction of graphene with 3d transition metal adatoms and the formation of localized magnetic moments by means of first-principles calculations. By comparing calculations within the generalized gradient approximation (GGA) to GGA+U we find that the electronic configuration and the adsorption geometries can be very sensitive to effects of local Coulomb interactions U in the transition metal d-orbitals. We find high-spin configurations being favorable for Cr and Mn adatoms independent of the functional. For Fe, Co and Ni different electronic configurations are realized depending on the value of the local Coulomb interaction strength U. Chemical control over the spin of the adatoms by hydrogenation is demonstrated: NiH and CoH adsorbed to graphene exhibit spin S=1/2 and S=1, respectively.

preprint2011arXiv

Two-dimensional Mott-Hubbard electrons in an artificial honeycomb lattice

Electrons in artificial lattices enable explorations of the impact of repulsive Coulomb interactions in a tunable system. We have trapped two-dimensional electrons belonging to a gallium arsenide quantum well in a nanofabricated lattice with honeycomb geometry. We probe the excitation spectrum in a magnetic field identifying novel collective modes that emerge from the Coulomb interaction in the artificial lattice as predicted by the Mott-Hubbard model. These observations allow us to determine the Hubbard gap and suggest the existence of a novel Coulomb-driven ground state. This approach offers new venues for the study of quantum phenomena in a controllable solid-state system.

preprint2011arXiv

Two-Site Kondo Effect in Atomic Chains

Linear CoCu_nCo clusters on Cu(111) are fabricated by means of atomic manipulation. They represent a two-site Kondo system with tunable interaction. Scanning tunneling spectroscopy reveals oscillations of the Kondo temperature T_K with the number n of Cu atoms for n>=3. Density functional calculations show that the Ruderman-Kittel-Kasuya-Yosida interaction mediated by the Cu chains causes the oscillations. Calculations find ferromagnetic and antiferromagnetic interaction for n=1 and 2, respectively. Both interactions lead to a decrease of T_K as experimentally observed.

preprint2010arXiv

A new route towards uniformly functionalized single-layer graphene

It is shown, by DFT calculations, that the uniform functionalization of upper layer of graphite by hydrogen or fluorine does not change essentially its bonding energy with the underlying layers, whereas the functionalization by phenyl groups decreases the bonding energy by a factor of approximately ten. This means that the functionalized monolayer in the latter case can be easily separated by mild sonication. According to our computational results, such layers can be cleaned up to pure graphene, as well as functionalized further up to 25% coverage, without essential difficulties. The energy gap within the interval from 0.5 to 3 eV can be obtained by such one-side funtionalization using different chemical species.

preprint2010arXiv

A theoretical analysis of the chemical bonding and electronic structure of graphene interacting with Group IA and Group VIIA elements

We propose a new class of materials, which can be viewed as graphene derivatives involving Group IA or Group VIIA elements, forming what we refer to as graphXene. We show that in several cases large band gaps can be found to open up, whereas in other cases a semimetallic behavior is found. Formation energies indicate that under ambient conditions, sp$^3$ and mixed sp$^2$/sp$^3$ systems will form. The results presented allow us to propose that by careful tuning of the relative concentration of the adsorbed atoms, it should be possible to tune the band gap of graphXene to take any value between 0 and 6.4 eV.

preprint2010arXiv

Adsorption of diatomic halogen molecules on graphene: A van der Waals density functional study

The adsorption of fluorine, chlorine, bromine, and iodine diatomic molecules on graphene has been investigated using density functional theory with taking into account nonlocal correlation effects by means of vdW-DF approach. It is shown that the van der Waals interaction plays a crucial role in the formation of chemical bonding between graphene and halogen molecules, and is therefore important for a proper description of adsorption in this system. In-plane orientation of the molecules has been found to be more stable than the orientation perpendicular to the graphene layer. In the cases of F$_2$, Br$_2$ and I$_2$ we also found an ionic contribution to the binding energy, slowly vanishing with distance. Analysis of the electronic structure shows that ionic interaction arises due to the charge transfer from graphene to the molecules. Furthermore, we found that the increase of impurity concentration leads to the conduction band formation in graphene due to interaction between halogen molecules. In addition, graphite intercalation by halogen molecules has been investigated. In the presence of halogen molecules the binding between graphite layers becomes significantly weaker, which is in accordance with the results of recent experiments on sonochemical exfoliation of intercalated graphite.

preprint2010arXiv

Atomistic simulations of structural and thermodynamic properties of bilayer graphene

We study the structural and thermodynamic properties of bilayer graphene, a prototype two-layer membrane, by means of Monte Carlo simulations based on the empirical bond order potential LCBOPII. We present the temperature dependence of lattice parameter, bending rigidity and high temperature heat capacity as well as the correlation function of out-of-plane atomic displacements. The thermal expansion coefficient changes sign from negative to positive above $\approx 400$ K, which is lower than previously found for single layer graphene and close to the experimental value of bulk graphite. The bending rigidity is twice as large than for single layer graphene, making the out-of-plane fluctuations smaller. The crossover from correlated to uncorrelated out-of-plane fluctuations of the two carbon planes occurs for wavevectors shorter than $\approx 3$ nm$^{-1}$

preprint2010arXiv

Bending modes, elastic constants and mechanical stability of graphitic systems

The thermodynamic and mechanical properties of graphitic systems are strongly dependent on the shear elastic constant C44. Using state-of-the-art density functional calculations, we provide the first complete determination of their elastic constants and exfoliation energies. We show that stacking misorientations lead to a severe lowering of C44 of at least one order of magnitude. The lower exfoliation energy and the lower C44 (more bending modes) suggest that flakes with random stacking should be easier to exfoliate than the ones with perfect or rhombohedral stacking. We also predict ultralow friction behaviour in turbostratic graphitic systems.

preprint2010arXiv

Correlated band theory of spin and orbital contributions to Dzyaloshinskii-Moriya interactions

A new approach for calculations of Dzyaloshinskii-Moriya interactions in molecules and crystals is proposed. It is based on the exact perturbation expansion of total energy of weak ferromagnets in the canting angle with the only assumption of local Hubbard-type interactions. This scheme leads to a simple and transparent analytical expression for Dzyaloshinskii-Moriya vector with a natural separation into spin and orbital contributions. The main problem was transferred to calculations of effective tight-binding parameters in the properly chosen basis including spin-orbit coupling. Test calculations for La$_2$CuO$_4$ give the value of canting angle in a good agreement with experimental data.

preprint2010arXiv

Density of states and zero Landau level probed through capacitance of graphene

We report capacitors in which a finite electronic compressibility of graphene dominates the electrostatics, resulting in pronounced changes in capacitance as a function of magnetic field and carrier concentration. The capacitance measurements have allowed us to accurately map the density of states D, and compare it against theoretical predictions. Landau oscillations in D are robust and zero Landau level (LL) can easily be seen at room temperature in moderate fields. The broadening of LLs is strongly affected by charge inhomogeneity that leads to zero LL being broader than other levels.

preprint2010arXiv

Dual-Fermion approach to Non-equilibrium strongly correlated problems

We present a generalization of the recently developed dual fermion approach introduced for correlated lattices to non-equilibrium problems. In its local limit, the approach has been used to devise an efficient impurity solver, the superperturbation solver for the Anderson impurity model (AIM). Here we show that the general dual perturbation theory can be formulated on the Keldysh contour. Starting from a reference Hamiltonian system, in which the time-dependent solution is found by exact diagonalization, we make a dual perturbation expansion in order to account for the relaxation effects from the fermionic bath. Simple test results for closed as well as open quantum systems in a fermionic bath are presented.

preprint2010arXiv

Effect of ligand substitution on the exchange interactions in {Mn12}-type single-molecule magnets

We investigate how the ligand substitution affects the intra-molecular spin exchange interactions, studying a prototypal family of single-molecule magnets comprising dodecanuclear cluster molecules [Mn12O12(COOR)16]. We identify a simple scheme based on accumulated Pauling electronegativity numbers (a.e.n.) of the carboxylate ligand groups (R). The redistribution of the electron density, controlled by a.e.n. of a ligand, changes the degree of hybridization between 3d electrons of manganese and 2p electrons of oxygen atoms, thus changing the exchange interactions. This scheme, despite its conceptual simplicity, provides a strong correlation with the exchange energies associated with carboxylate bridges, and is confirmed by the electronic structure calculations taking into account the Coulomb correlations in magnetic molecules.

preprint2010arXiv

Electron density distribution and screening in rippled graphene sheets

Single-layer graphene sheets are typically characterized by long-wavelength corrugations (ripples) which can be shown to be at the origin of rather strong potentials with both scalar and vector components. We present an extensive microscopic study, based on a self-consistent Kohn-Sham-Dirac density-functional method, of the carrier density distribution in the presence of these ripple-induced external fields. We find that spatial density fluctuations are essentially controlled by the scalar component, especially in nearly-neutral graphene sheets, and that in-plane atomic displacements are as important as out-of-plane ones. The latter fact is at the origin of a complicated spatial distribution of electron-hole puddles which has no evident correlation with the out-of-plane topographic corrugations. In the range of parameters we have explored, exchange and correlation contributions to the Kohn-Sham potential seem to play a minor role.

preprint2010arXiv

Electron energy spectrum of the spin-liquid state in a frustrated Hubbard model

Non-local correlation effects in the half-filled Hubbard model on an isotropic triangular lattice are studied within a spin polarized extension of the dual fermion approach. A competition between the antiferromagnetic non-collinear and the spin liquid states is strongly enhanced by an incorporation of a k-dependent self-energy beyond the local dynamical mean-field theory. The dual fermion correc- tions drastically decrease the energy of a spin liquid state while leaving the non-collinear magnetic states almost non-affected. This makes the spin liquid to become a preferable state in a certain interval of interaction strength of an order of the magnitude of a bandwidth. The spectral function of the spin-liquid Mott insulator is determined by a formation of local singlets which results in the energy gap of about twice larger than that of the 120 degrees antiferromagnetic Neel state.

preprint2010arXiv

First-principles Theory of Nonlocal Screening in Graphene

Using the quasiparticle self-consistent GW (QSGW) and local-density (LD) approximations, we calculate the q-dependent static dielectric function, and derive an effective 2D dielectric function corresponding to screening of point charges. In the q-to-0 limit, the 2D function is found to scale approximately as the square root of the macroscopic dielectric function. Its value is ~4, a factor approximately 1.5 larger than predictions of Dirac model. Both kinds of dielectric functions depend strongly on q, in contrast with the Dirac model. The QSGW approximation is shown to describe QP levels very well, with small systematic errors analogous to bulk sp semiconductors. Local-field effects are rather more important in graphene than in bulk semiconductors.

preprint2010arXiv

Flexuron, a self-trapped state of electron in crystalline membranes

Self-trapping of an electron due to its interaction with bending fluctuations in a flexible crystalline membrane is considered. Due to the dependence of the electron energy on the corrugations of the membrane, the electron can create around itself an anomalously flat (or anomalously corrugated, depending on the sign of the interaction constant) region and be confined there. Using the Feynman path integral approach, the autolocalization energy and the size of the self-trapped state (flexuron) are estimated. It is shown that typically the size of the flexuron is of the order of the wavelength of fluctuations at the border between harmonic and anharmonic regimes. The flexuron states are connected with the fluctuation tail of the electron density-of-states, the asymptotic behavior of this tail being determined by the exponent of the renormalized bending rigidity.

preprint2010arXiv

Fluorographene: Two Dimensional Counterpart of Teflon

We report a stoichiometric derivative of graphene with a fluorine atom attached to each carbon. Raman, optical, structural, micromechanical and transport studies show that the material is qualitatively different from the known graphene-based nonstoichiometric derivatives. Fluorographene is a high-quality insulator (resistivity >10^12 Ohm per square) with an optical gap of 3 eV. It inherits the mechanical strength of graphene, exhibiting Young's modulus of 100 N/m and sustaining strains of 15%. Fluorographene is inert and stable up to 400C even in air, similar to Teflon.

preprint2010arXiv

Gauge fields in graphene

The physics of graphene is acting as a bridge between quantum field theory and condensed matter physics due to the special quality of the graphene quasiparticles behaving as massless two dimensional Dirac fermions. Moreover, the particular structure of the 2D crystal lattice sets the arena to study and unify concepts from elasticity, topology and cosmology. In this paper we analyze these connections combining a pedagogical, intuitive approach with a more rigorous formalism when required.

preprint2010arXiv

Limits on electron quality in suspended graphene due to flexural phonons

The temperature dependence of the mobility in suspended graphene samples is investigated. In clean samples, flexural phonons become the leading scattering mechanism at temperature $T \gtrsim 10\,\,$K, and the resistivity increases quadratically with $T$. Flexural phonons limit the intrinsic mobility down to a few $\text{m}^2/\text{Vs}$ at room $T$. Their effect can be eliminated by applying strain or placing graphene on a substrate.

preprint2010arXiv

Non-spherical shapes of capsules within a fourth-order curvature model

We minimize a discrete version of the fourth-order curvature based Landau free energy by extending Brakke's Surface Evolver. This model predicts spherical as well as non-spherical shapes with dimples, bumps and ridges to be the energy minimizers. Our results suggest that the buckling and faceting transitions, usually associated with crystalline matter, can also be an intrinsic property of non-crystalline membranes.

preprint2010arXiv

Orbitally controlled Kondo effect of Co ad-atoms on graphene

Based on ab-initio calculations we identify possible scenarios for the Kondo effect due to Co ad-atoms on graphene. General symmetry arguments show that for magnetic atoms in high-symmetry positions, the Kondo effect in graphene is controlled not only by the spin but also by the orbital degree of freedom. For a Co atom absorbed on top of a carbon atom, the Kondo effect is quenched by spin-orbit coupling below an energy scale of $\sim 15$\,K. For Co with spin $S=1/2$ located in the center of a hexagon, an SU(4) Kondo model describes the entanglement of orbital moment and spin at higher energies, while below \sim 60$\,meV spin-orbit coupling leads to a more conventional SU(2) Kondo effect. The interplay of the orbital Co physics and the peculiar band-structure of graphene is directly accessible in Fourier transform tunneling spectroscopy or in the gate-voltage dependence of the Kondo temperature displaying a very strong, characteristic particle-hole asymmetry.

preprint2010arXiv

Quantitative determination of spin-dependent quasiparticle lifetimes and electronic correlations in hcp cobalt

We report on a quantitative investigation of the spin-dependent quasiparticle lifetimes and electron correlation effects in ferromagnetic hcp Co(0001) by means of spin and angle-resolved photoemission spectroscopy. The experimental spectra are compared in detail to state-of-the-art many-body calculations within the dynamical mean field theory and the three-body scattering approximation, including a full calculation of the one-step photoemission process. From this comparison we conclude that although strong local many-body Coulomb interactions are of major importance for the qualitative description of correlation effects in Co, more sophisticated many-body calculations are needed in order to improve the quantitative agreement between theory and experiment, in particular concerning the linewidths. The quality of the overall agreement obtained for Co indicates that the effect of non-local correlations becomes weaker with increasing atomic number.

preprint2010arXiv

Resonant scattering by realistic impurities in graphene

We develop a first-principles theory of resonant impurities in graphene and show that a broad range of typical realistic impurities leads to the characteristic sublinear dependence of the conductivity on the carrier concentration. By means of density functional calculations various organic groups as well as ad-atoms like H absorbed to graphene are shown to create midgap states within +-0.03eV around the neutrality point. A low energy tight-binding description is mapped out. Boltzmann transport theory as well as a numerically exact Kubo formula approach yield the conductivity of graphene contaminated with these realistic impurities in accordance with recent experiments.

preprint2010arXiv

Self-Consistent Screening Approximation for Flexible Membranes: Application to Graphene

Crystalline membranes at finite temperatures have an anomalous behavior of the bending rigidity that makes them more rigid in the long wavelength limit. This issue is particularly relevant for applications of graphene in nano- and micro-electromechanical systems. We calculate numerically the height-height correlation function $G(q)$ of crystalline two-dimensional membranes, determining the renormalized bending rigidity, in the range of wavevectors $q$ from $10^{-7}$ Å$^{-1}$ till 10 Å$^{-1}$ in the self-consistent screening approximation (SCSA). For parameters appropriate to graphene, the calculated correlation function agrees reasonably with the results of atomistic Monte Carlo simulations for this material within the range of $q$ from $10^{-2}$ Å$^{-1}$ till 1 Å$^{-1}$. In the limit $q\rightarrow 0$ our data for the exponent $η$ of the renormalized bending rigidity $κ_R(q)\propto q^{-η}$ is compatible with the previously known analytical results for the SCSA $η\simeq 0.82$. However, this limit appears to be reached only for $q<10^{-5}$ Å$^{-1}$ whereas at intermediate $q$ the behavior of $G(q)$ cannot be described by a single exponent.

preprint2010arXiv

Stable and fast semi-implicit integration of the stochastic Landau-Lifshitz equation

We propose new semi-implicit numerical methods for the integration of the stochastic Landau-Lifshitz equation with built-in angular momentum conservation. The performance of the proposed integrators is tested on the 1D Heisenberg chain. For this system, our schemes show better stability properties and allow us to use considerably larger time steps than standard explicit methods. At the same time, these semi-implicit schemes are also of comparable accuracy to and computationally much cheaper than the standard midpoint implicit method. The results are of key importance for atomistic spin dynamics simulations and the study of spin dynamics beyond the macro spin approximation.

preprint2010arXiv

The Effect of Cluster Formation on Graphene Mobility

We investigate the effect of gold (Au) atoms in the form of both point-like charged impurities and clusters on the transport properties of graphene. Cryogenic deposition (18 K) of Au decreases the mobility and shifts the Dirac point in a manner that is consistent with scattering from point-like charged impurities. Increasing the temperature to room temperature promotes the formation of clusters, which is verified with atomic force microscopy. We find that for a fixed amount of Au impurities, the formation of clusters enhances the mobility and causes the Dirac point to shift back towards zero.

preprint2010arXiv

Theory of optically forbidden d-d transitions in strongly correlated crystals

A general multiband formulation of linear and non-linear optical response functions for realistic models of correlated crystals is presented. Dipole forbidden d-d optical transitions originate from the vertex functions, which we consider assuming locality of irreducible four-leg vertex. The unified formulation for second- and third-order response functions in terms of the three-leg vertex is suitable for practical calculations in solids. We illustrate the general approach by consideration of intraatomic spin-flip contributions, with the energy of 2J, where J is a Hund exchange, in the simplest two-orbital model.

preprint2009arXiv

About the strength of correlation effects in the electronic structure of iron

The strength of electronic correlation effects in the spin-dependent electronic structure of ferromagnetic bcc Fe(110) has been investigated by means of spin and angle-resolved photoemission spectroscopy. The experimental results are compared to theoretical calculations within the three-body scattering approximation and within the dynamical mean-field theory, together with one-step model calculations of the photoemission process. This comparison indicates that the present state of the art many-body calculations, although improving the description of correlation effects in Fe, give too small mass renormalizations and scattering rates thus demanding more refined many-body theories including non-local fluctuations.

preprint2009arXiv

Accurate electronic band gap of pure and functionalized graphane from GW calculations

Using the GW approximation, we study the electronic structure of the recently synthesized hydrogenated graphene, named graphane. For both conformations, the minimum band gap is found to be direct at the $Γ$ point, and it has a value of 5.4 eV in the stable chair conformation, where H atoms attach C atoms alternatively on opposite sides of the two dimensional carbon network. In the meta-stable boat conformation the energy gap is 4.9 eV. Then, using a supercell approach, the electronic structure of graphane was modified by introducing either an hydroxyl group or an H vacancy. In this last case, an impurity state appears at about 2 eV above the valence band maximum.

preprint2009arXiv

Aharonov-Bohm effect in undoped graphene: Magnetotransport via evanescent waves

Using conformal mapping technique, compact and general analytic expressions for the effects of magnetic fluxes on conductance and Fano factor of undoped graphene nanoflakes in pseudodiffusive regime are derived.

preprint2009arXiv

Analytical approximation for single-impurity Anderson model

We have applied the recently developed dual fermion technique to the spectral properties of single-band Anderson impurity problem (SIAM). In our approach a series expansion is constructed in vertices of the corresponding atomic Hamiltonian problem. This expansion contains a small parameter in two limiting cases: in the weak coupling case ($U/t \to 0$), due to the smallness of the irreducible vertices, and near the atomic limit ($U/t \to \infty$), when bare propagators are small. Reasonable results are obtained also for the most interesting case of strong correlations ($U \approx t$). The atomic problem of the Anderson impurity model has a degenerate ground state, so the application of the perturbation theory is not straightforward. We construct a special approach dealing with symmetry-broken ground state of the renormalized atomic problem. Formulae for the first-order dual diagram correction are obtained analytically in the real-time domain. Most of the Kondo-physics is reproduced: logarithmic contributions to the self energy arise, Kondo-like peak at the Fermi level appears, and the Friedel sum rule is fulfilled. Our approach describes also renormalization of atomic resonances due to hybridization with a conduction band. A generalization of the proposed scheme to a multi-orbital case can be important for the realistic description of correlated solids.

preprint2009arXiv

Correlated Electrons Step-by-Step: Itinerant-to-Localized Transition of Fe Impurities in Free-Electron Metal Hosts

High-resolution photoemission spectroscopy and realistic ab-initio calculations have been employed to analyze the onset and progression of d-sp hybridization in Fe impurities deposited on alkali metal films. The interplay between delocalization, mediated by the free-electron environment, and Coulomb interaction among d-electrons gives rise to complex electronic configurations. The multiplet structure of a single Fe atom evolves and gradually dissolves into a quasiparticle peak near the Fermi level with increasing the host electron density. The effective multi-orbital impurity problem within the exact diagonalization scheme describes the whole range of hybridizations.

preprint2009arXiv

Destruction of graphene by metal adatoms

The formation energies for mono- and bivacancies in graphene in the presence of adatoms of various metals and small metallic clusters have been calculated. It is shown that transition metal impurities, such as iron, nickel and, especially, cobalt reduce dramatically the vacancy formation energies whereas gold impurities have almost no effect on characteristics of the vacancies. This results highlight that special measures are required in order to protect graphene from damage by transition metal leads.

preprint2009arXiv

Effect of high-k environment on charge carrier mobility in graphene

It is widely assumed that the dominant source of scattering in graphene is charged impurities in a substrate. We have tested this conjecture by studying graphene placed on various substrates and in high-k media. Unexpectedly, we have found no significant changes in carrier mobility either for different substrates or by using glycerol, ethanol and water as a top dielectric layer. This suggests that Coulomb impurities are not the scattering mechanism that limits the mean free path currently attainable for graphene on a substrate.

preprint2009arXiv

Efficient perturbation theory for quantum lattice models

We present a novel approach to long-range correlations beyond dynamical mean-field theory through a ladder approximation to dual fermions. The new technique is applied to the two-dimensional Hubbard model. We demonstrate that the transformed perturbation series for the nonlocal dual fermions has superior convergence properties over standard diagrammatic techniques. The critical Neel temperature of the mean-field solution is suppressed in the ladder approximation, in accordance with quantum Monte-Carlo (QMC) results. An illustration of how the approach captures and allows to distinguish short- and long-range correlations is given.

preprint2009arXiv

Energy gaps, topological insulator state and zero-field quantum Hall effect in graphene by strain engineering

Among many remarkable qualities of graphene, its electronic properties attract particular interest due to a massless chiral character of charge carriers, which leads to such unusual phenomena as metallic conductivity in the limit of no carriers and the half-integer quantum Hall effect (QHE) observable even at room temperature [1-3]. Because graphene is only one atom thick, it is also amenable to external influences including mechanical deformation. The latter offers a tempting prospect of controlling graphene's properties by strain and, recently, several reports have examined graphene under uniaxial deformation [4-8]. Although the strain can induce additional Raman features [7,8], no significant changes in graphene's band structure have been either observed or expected for realistic strains of approx. 10% [9-11]. Here we show that a designed strain aligned along three main crystallographic directions induces strong gauge fields [12-14] that effectively act as a uniform magnetic field exceeding 10 T. For a finite doping, the quantizing field results in an insulating bulk and a pair of countercirculating edge states, similar to the case of a topological insulator [15-20]. We suggest realistic ways of creating this quantum state and observing the pseudo-magnetic QHE. We also show that strained superlattices can be used to open significant energy gaps in graphene's electronic spectrum.

preprint2009arXiv

gamma-Mn at the border between weak and strong correlations

We investigate the role of magnetic fluctuations in the spectral properties of paramagnetic gamma-Mn. Two methods are employed. The Local Density Approximation plus Dynamical Mean-Field Theory together with the numerically exact quantum Monte-Carlo solver is used as a reference for the spectral properties. Then the same scheme is used with the computationally less demanding perturbative spin-polarized fluctuation-exchange solver in combination with the Disordered Local Moment approach, and photoemission spectra are calculated within the one-step model. It is shown that the formation of local magnetic moments in gamma-Mn is very sensitive to the value of Hund's exchange parameter. Comparison with the experimental photoemission spectra demonstrates that gamma-Mn is a strongly correlated system, with the Hubbard band formation, which cannot be described by the perturbative approach. However, minor change of parameters would transform it into a weakly correlated system.

preprint2009arXiv

Generating quantizing pseudomagnetic fields by bending graphene ribbons

We analyze the mechanical deformations that are required to create uniform pseudomagnetic fields in graphene. It is shown that, if a ribbon is bent in-plane into a circular arc, this can lead to fields exceeding 10T, which is sufficient for the observation of pseudo-Landau quantization. The arc geometry is simpler than those suggested previously and, in our opinion, has much better chances to be realized experimentally soon. The effects of a scalar potential induced by dilatation in this geometry is shown to be negligible.

preprint2009arXiv

Half-metallicity in NiMnSb: a Variational Cluster Approach with ab-initio parameters

Electron correlation effects in the half-metallic ferromagnet NiMnSb are investigated within a combined density functional and many-body approach. Starting from a realistic multi-orbital Hubbard-model including Mn and Ni-d orbitals, the many-body problem is addressed via the Variational Cluster Approach. The density of states obtained in the calculation shows a strong spectral weight transfer towards the Fermi level in the occupied conducting majority spin channel with respect to the uncorrelated case, as well as states with vanishing quasiparticle weight in the minority spin gap. Although the two features produce competing effects, the overall outcome is a strong reduction of the spin polarisation at the Fermi level with respect to the uncorrelated case. This result emphasizes the importance of correlation in this material.

preprint2009arXiv

Majority-spin non-quasiparticle states in half-metallic ferrimagnet Mn$_2$VAl

The density of non-quasiparticle states in the ferrimagnetic full-Heuslers Mn$_2$VAl alloy is calculated from first principles upon appropriate inclusion of correlations. In contrast to most half-metallic compounds, this material displays an energy gap in the majority-spin spectrum. For this situation, non-quasiparticle states are located below the Fermi level, and should be detectable by spin-polarized photoemission. This opens a new way to study many-body effects in spintronic-related materials.

preprint2009arXiv

Relevance of complete Coulomb interaction matrix for the Kondo problem: Co impurity on Cu(111)

The electronic structure of a prototype Kondo system, a cobalt impurity in a copper host is calculated with accurate taking into account of correlation effects on the Co atom. Using the recently developed continuous-time QMC technique, it is possible to describe the Kondo resonance with a complete four-index Coulomb interaction matrix. This opens a way for completely first-principle calculations of the Kondo temperature. We have demonstrated that a standard practice of using a truncated Hubbard Hamiltonian to consider the Kondo physics can be quantitatively inadequate.

preprint2009arXiv

Scaling Properties of Flexible Membranes from Atomistic Simulations: Application to Graphene

Structure and thermodynamics of crystalline membranes are characterized by the long wavelength behavior of the normal-normal correlation function G(q). We calculate G(q) by Monte Carlo and Molecular Dynamics simulations for a quasi-harmonic model potential and for a realistic potential for graphene. To access the long wavelength limit for finite-size systems (up to 40000 atoms) we introduce a Monte Carlo sampling based on collective atomic moves (wave moves). We find a power-law behaviour $G(q)\propto q^{-2+η}$ with the same exponent $η\approx 0.85$ for both potentials. This finding supports, from the microscopic side, the adequacy of the scaling theory of membranes in the continuum medium approach, even for an extremely rigid material like graphene.

preprint2009arXiv

Topological defects and shape of aromatic self-assembled vesicles

We show that the stacking of flat aromatic molecules on a curved surface results in topological defects. We consider, as an example, spherical vesicles, self-assembled from molecules with 5- and 6-thiophene cores. We predict that the symmetry of the molecules influences the number of topological defects and the resulting equilibrium shape.

preprint2008arXiv

Correlation effects in the electronic structure of Mn$_4$ molecular magnet

We present joint theoretical-experimental study of the correlation effects in the electronic structure of (pyH)$_3$[Mn$_4$O$_3$Cl$_7$(OAc)$_3$]$\cdot$2MeCN molecular magnet (Mn$_4$). Describing the many-body effects by cluster dynamical mean-field theory, we find that Mn$_4$ is predominantly Hubbard insulator with strong electron correlations. The calculated electron gap (1.8 eV) agrees well with the results of optical conductivity measurements, while other methods, which neglect many-body effects or treat them in a simplified manner, do not provide such an agreement. Strong electron correlations in Mn$_4$ may have important implications for possible future applications.

preprint2008arXiv

Effect of magnetic state on the $γ-α$ transition in iron: First-principle calculations of the Bain transformation path

Energetics of the fcc ($γ$) - bcc ($α$) lattice transformation by the Bain tetragonal deformation is calculated for both magnetically ordered and paramagnetic (disordered local moment) states of iron. The first-principle computational results manifest a relevance of the magnetic order in a scenario of the $γ$ - $α$ transition and reveal a special role of the Curie temperature of $α$-Fe, $T_C$, where a character of the transformation is changed. At a cooling down to the temperatures $T < T_C$ one can expect that the transformation is developed as a lattice instability whereas for $T > T_C$ it follows a standard mechanism of creation and growth of an embryo of the new phase. It explains a closeness of $T_C$ to the temperature of start of the martensitic transformation, $M_s$.

preprint2008arXiv

Electron energy level statistics in graphene quantum dots

Motivated by recent experimental observations of size quantization of electron energy levels in graphene quantum dots \cite{ponomarenko} we investigate the level statistics in the simplest tight-binding model for different dot shapes by computer simulation. The results are in a reasonable agreement with the experiment which confirms qualitatively interpretation of observed level statistics in terms of "Dirac billiards" without taking into account many-body effects. It is shown that edge effects are in general sufficient to produce the observed level distribution and that even strong bulk disorder does not change the results drastically.

preprint2007arXiv

Anharmonic magnetic deformation of self-assembled molecular nanocapsules

High magnetic fields were used to deform spherical nanocapsules, self-assembled from bola-amphiphilic sexithiophene molecules. At low fields the deformation -- measured through linear birefringence -- scales quadratically with the capsule radius and with the magnetic field strength. These data confirm a long standing theoretical prediction (W. Helfrich, Phys. Lett. {\bf 43A}, 409 (1973)), and permits the determination of the bending rigidity of the capsules as (2.6$\pm$0.8)$\times 10^{-21}$ J. At high fields, an enhanced rigidity is found which cannot be explained within the Helfrich model. We propose a complete form of the free energy functional that accounts for this behaviour, and allows discussion of the formation and stability of nanocapsules in solution.

preprint2007arXiv

Cluster Dual Fermion Approach to Nonlocal Correlations

We formulate a general cluster Dual Fermion Approach to nonlocal correlations in crystals. The scheme allows the treatment of long-range correlations beyond cluster DMFT and nonlocal effects in realistic calculations of multiorbital systems. We show that the the simplest approximation exactly corresponds to free cluster DMFT. We further consider the relation between the two-particle Green functions in real and dual variables. We apply this approach by calculating the Green function of the Hubbard model in one dimension starting from the two-site cluster DMFT solution. The result agrees well with the Green function obtained from a DMRG calculation.

preprint2007arXiv

Conductance quantization in graphene nanoribbons: Adiabatic approximation

A theory of electron states for graphene nanoribbons with a smoothly varying width is developed. It is demonstrated that the standard adiabatic approximation allowing to neglect the mixing of different standing waves is more restrictive for the massless Dirac fermions in graphene than for the conventional electron gas. For the case of zigzag boundary conditions, one can expect a well-pronounced conductance quantization only for highly excited states. This difference is related to the relativistic Zitterbewegung effect in graphene.

preprint2007arXiv

Decoherence by a spin thermal bath: Role of the spin-spin interactions

We study the decoherence of two coupled spins that interact with a spin-bath environment. It is shown that the connectivity and the coupling strength between the spins in the environment are of crucial importance for the decoherence of the central system. Changing the connectivity or coupling strenghts changes the decoherence of the central system from Gaussian to exponential decay law.

preprint2007arXiv

Generalized kinetic equations for charge carriers in graphene

A system of generalized kinetic equations for the distribution functions of two-dimensional Dirac fermions scattered by impurities is derived in the Born approximation with respect to short-range impurity potential. It is proven that the conductivity following from classical Boltzmann equation picture, where electrons or holes have scattering amplitude reduced due chirality, is justified except for an exponentially narrow range of chemical potential near the conical point. When in this range, creation of infinite number of electron-hole pairs related to quasi-relativistic nature of electrons in graphene results in a renormalization of minimal conductivity as compared to the Boltzmann term and logarithmic corrections in the conductivity similar to the Kondo effect.

preprint2007arXiv

Graphene: new bridge between condensed matter physics and quantum electrodynamics

Graphene is the first example of truly two-dimensional crystals - it's just one layer of carbon atoms. It turns out to be a gapless semiconductor with unique electronic properties resulting from the fact that charge carriers in graphene demonstrate charge-conjugation symmetry between electrons and holes and possess an internal degree of freedom similar to ``chirality'' for ultrarelativistic elementary particles. It provides unexpected bridge between condensed matter physics and quantum electrodynamics (QED). In particular, the relativistic Zitterbewegung leads to the minimum conductivity of order of conductance quantum $e^2/h$ in the limit of zero doping; the concept of Klein paradox (tunneling of relativistic particles) provides an essential insight into electron propagation through potential barriers; vacuum polarization around charge impurities is essential for understanding of high electron mobility in graphene; index theorem explains anomalous quantum Hall effect.

preprint2007arXiv

Half-Metallic Ferromagnetism and the spin polarization in CrO$_2$

We present electronic structure calculations in combination with local and non-local many-body correlation effects for the half-metallic ferromagnet CrO$_2$. Finite-temperature Dynamical Mean Field Theory results show the existence of non-quasiparticle states, which were recently observed as almost currentless minority spin states near the Fermi energy in resonant scattering experients. At zero temperatures, Variational Cluster Approach calculations support the half-metallic nature of CrO$_2$ as seen in superconducting point contact spectroscopy. The combination of these two techniques allowed us to qualitatively describe the spin-polarization in CrO$_2$.

preprint2007arXiv

Intrinsic nanoscale inhomogeneity in ordering systems due to elastic-mediated interactions

Phase diagram and pattern formation in two-dimensional Ising model with coupling between order parameter and lattice vibrations is investigated by Monte-Carlo simulations. It is shown that if the coupling is strong enough (or phonons are soft enough) a short-range order exists in disordered phase for a broader temperature interval. Different types of this short-range order (stripe-like, checkboard-like, etc.) depending on the temperature and model parameters are investigated. With further increase of the coupling, a reconstruction of the ground state happens and new ordered phases appear at low enough temperatures.

preprint2007arXiv

Intrinsic ripples in graphene

The stability of two-dimensional (2D) layers and membranes is subject of a long standing theoretical debate. According to the so called Mermin-Wagner theorem, long wavelength fluctuations destroy the long-range order for 2D crystals. Similarly, 2D membranes embedded in a 3D space have a tendency to be crumpled. These dangerous fluctuations can, however, be suppressed by anharmonic coupling between bending and stretching modes making that a two-dimensional membrane can exist but should present strong height fluctuations. The discovery of graphene, the first truly 2D crystal and the recent experimental observation of ripples in freely hanging graphene makes these issues especially important. Beside the academic interest, understanding the mechanisms of stability of graphene is crucial for understanding electronic transport in this material that is attracting so much interest for its unusual Dirac spectrum and electronic properties. Here we address the nature of these height fluctuations by means of straightforward atomistic Monte Carlo simulations based on a very accurate many-body interatomic potential for carbon. We find that ripples spontaneously appear due to thermal fluctuations with a size distribution peaked around 70 Åwhich is compatible with experimental findings (50-100 Å) but not with the current understanding of stability of flexible membranes. This unexpected result seems to be due to the multiplicity of chemical bonding in carbon.

preprint2007arXiv

Molecular Doping of Graphene

Graphene, a one-atom thick zero gap semiconductor [1, 2], has been attracting an increasing interest due to its remarkable physical properties ranging from an electron spectrum resembling relativistic dynamics [3-12] to ballistic transport under ambient conditions [1-4]. The latter makes graphene a promising material for future electronics and the recently demonstrated possibility of chemical doping without significant change in mobility has improved graphene's prospects further [13]. However, to find optimal dopants and, more generally, to progress towards graphene-based electronics requires understanding the physical mechanism behind the chemical doping, which has been lacking so far. Here, we present the first joint experimental and theoretical investigation of adsorbates on graphene. We elucidate a general relation between the doping strength and whether or not adsorbates have a magnetic moment: The paramagnetic single NO2 molecule is found to be a strong acceptor, whereas its diamagnetic dimer N2O4 causes only weak doping. This effect is related to the peculiar density of states of graphene, which provides an ideal situation for model studies of doping effects in semiconductors. Furthermore, we explain recent results on its "chemical sensor" properties, in particular, the possibility to detect a single NO2 molecule [13].

preprint2006arXiv

Electron Correlations and the Minority-Spin Band Gap in Half-Metallic Heusler Alloys

Electron-electron correlations affect the band gap of half-metallic ferromagnets by introducing non-quasiparticle states just above the Fermi level. In contrast to the spin-orbit coupling, a large asymmetric non-quasiparticle spectral weight is present in the minority-spin channel, leading to a peculiar finite-temperature spin depolarization effects. Using recently developed first-principle dynamical mean field theory, we investigate these effects for the half-metallic ferrimagnetic Heusler compound FeMnSb. We discuss depolarization effects in terms of strength of local Coulomb interaction $U$ and temperature in FeMnSb. We propose Ni$_{1-x}$Fe$_{x}$MnSb alloys as a perspective materials to be used in spin-valve structures and for experimental search of non-quasiparticle states in half-metallic materials.

preprint2006arXiv

Impurity States in Graphene

Defects in graphene are of crucial importance for its electronic and magnetic properties. Here impurity effects on the electronic structure of surrounding carbon atoms are considered and the distribution of the local densities of states (LDOS) is calculated. As the full range from near field to the asymptotic regime is covered, our results are directly accessible by scanning tunnelling microscopy (STM). We also include exchange scattering at magnetic impurities and eludicate how strongly spin polarized impurity states arise.

preprint2005arXiv

Multiple scattering formalism for correlated systems: A KKR+DMFT approach

We present a charge and self-energy self-consistent computational scheme for correlated systems based on the Korringa-Kohn-Rostoker (KKR) multiple scattering theory with the many-body effects described by the means of dynamical mean field theory (DMFT). The corresponding local multi-orbital and energy dependent self-energy is included into the set of radial differential equations for the single-site wave functions. The KKR Green's function is written in terms of the multiple scattering path operator, the later one being evaluated using the single-site solution for the $t$-matrix that in turn is determined by the wave functions. An appealing feature of this approach is that it allows to consider local quantum and disorder fluctuations on the same footing. Within the Coherent Potential Approximation (CPA) the correlated atoms are placed into a combined effective medium determined by the dynamical mean field theory (DMFT) self-consistency condition. Results of corresponding calculations for pure Fe, Ni and Fe$_{x}$Ni$_{1-x}$ alloys are presented.

preprint2005arXiv

Role of the d-f Coulomb interaction in intermediate valence and Kondo systems: a numerical renormalization group study

Using numerical renormalization group method the temperature dependences of the magnetic susceptibility $χ(T)$ and specific heat $C(T)$ are obtained in the single-impurity Anderson model with inclusion of the $d-f$ Coulomb interaction. It is shown that the exciton effects owing to this interaction can change considerably the dependence $C(T)$ in comparison with the standard Anderson model at not too low temperatures, whereas the dependence $χ(T)$ remains universal. The renormalization of the effective hybridization parameter and $f$-level position, which is connected with the $d-f$ interaction, is calculated, a satisfactory agreement with the Hartree-Fock approximation being derived.

preprint2005arXiv

Solvent Driven Formation of Bolaamphiphilic Vesicles

We show that a spontaneous bending of single layer bolaamphiphiles results from the frustration due to the competition between core-core and tail-solvent interactions. We find that spherical vesicles are stable under rather general assumptions on these interactions described within the Flory-Huggins theory. We consider also the deformation of the vesicles in an external magnetic field that has been recently experimentally observed.

preprint2005arXiv

Two-Dimensional Gas of Massless Dirac Fermions in Graphene

Electronic properties of materials are commonly described by quasiparticles that behave as non-relativistic electrons with a finite mass and obey the Schroedinger equation. Here we report a condensed matter system where electron transport is essentially governed by the Dirac equation and charge carriers mimic relativistic particles with zero mass and an effective "speed of light" c* ~10^6m/s. Our studies of graphene - a single atomic layer of carbon - have revealed a variety of unusual phenomena characteristic of two-dimensional (2D) Dirac fermions. In particular, we have observed that a) the integer quantum Hall effect in graphene is anomalous in that it occurs at half-integer filling factors; b) graphene's conductivity never falls below a minimum value corresponding to the conductance quantum e^2/h, even when carrier concentrations tend to zero; c) the cyclotron mass m of massless carriers with energy E in graphene is described by equation E =mc*^2; and d) Shubnikov-de Haas oscillations in graphene exhibit a phase shift of pi due to Berry's phase.

preprint2004arXiv

Nonquasiparticle states in half-metallic ferromagnets

Anomalous magnetic and electronic properties of the half-metallic ferromagnets (HMF) have been discussed. The general conception of the HMF electronic structure which take into account the most important correlation effects from electron-magnon interactions, in particular, the spin-polaron effects, is presented. Special attention is paid to the so called non-quasiparticle (NQP) or incoherent states which are present in the gap near the Fermi level and can give considerable contributions to thermodynamic and transport properties. Prospects of experimental observation of the NQP states in core-level spectroscopy is discussed. Special features of transport properties of the HMF which are connected with the absence of one-magnon spin-flip scattering processes are investigated. The temperature and magnetic field dependences of resistivity in various regimes are calculated. It is shown that the NQP states can give a dominate contribution to the temperature dependence of the impurity-induced resistivity and in the tunnel junction conductivity. First principle calculations of the NQP-states for the prototype half-metallic material NiMnSb within the local-density approximation plus dynamical mean field theory (LDA+DMFT) are presented.

preprint2003arXiv

Nonquasiparticle states in half-metallic ferromagnet NiMnSb

Nonquasiparticle states above the Fermi energy are studied by first-principle dynamical mean field calculations for a prototype half-metallic ferromagnet, NiMnSb. We present a quantitative evaluation of the spectral weight of this characteristic feature and discuss the possible experimental investigation (BIS, NMR, STM and Andreev reflection) to clarify the existence of these states.

preprint2002arXiv

Destruction of the Kondo effect by a local measurement

We show that the local spin measurement which decoheres the localized spin in a Kondo system, suppresses the Abrikosov-Suhl resonance and destroys the Kondo effect. This happens due to elimination of the entanglement between the localized spin and the conduction electrons, and differs essentially from smearing of the resonance by dissipation. Considering decoherence by a spin bath, we predict that the Kondo effect disappears when the Kondo temperature becomes smaller than the coupling with a bath. This effect can be detected in experiments on ``quantum corrals'' or quantum dots doped by impurities with internal degrees of freedom.

preprint2002arXiv

Robustness of the van Hove scenario for high-Tc superconductors

The pinning of the Fermi level to the van Hove singularity and the formation of flat bands in the two-dimensional t-t' Hubbard model is investigated by the renormalization group technique. The "van-Hove" scenario of non-Fermi-liquid behavior for high-Tc compounds can take place in a broad enough range of the hole concentrations. The results are in qualitative agreement with the recent ARPES data on La2CuO4.

preprint2001arXiv

Effects of the Van Hove singularities on magnetism and superconductivity in the t-t' Hubbard model: a parquet approach

The t-t' Hubbard model for the Fermi level near the Van Hove singularity is considered within the renormalization group and parquet approaches. The interplay of ferromagnetic, antiferromagnetic, and superconducting channels is investigated, and the phase diagram of the model is constructed. In comparison with previous approaches, the account of ferromagnetic fluctuations suppresses superconducting pairing and, vice versa, the influence of the Cooper channel decreases the Curie temperature, so that the Stoner criterion is inapplicable even qualitatively.

preprint1999arXiv

On the self-consistent spin-wave theory of layered Heisenberg magnets

The versions of the self-consistent spin-wave theories (SSWT) of two-dimensional (2D) Heisenberg ferro- and antiferromagnets with a weak interlayer coupling and/or magnetic anisotropy, that are based on the non-linear Dyson-Maleev, Schwinger, and combined boson-pseudofermion representations, are analyzed. Analytical results for the temperature dependences of (sublattice) magnetization and short-range order parameter, and the critical points are obtained. The influence of external magnetic field is considered. Fluctuation corrections to SSWT are calculated within a random-phase approximation which takes into account correctly leading and next-leading logarithmic singularities. These corrections are demonstrated to improve radically the agreement with experimental data on layered perovskites and other systems. Thus an account of these fluctuations provides a quantitative theory of layered magnets.

preprint1997arXiv

1/N Expansion for Critical Exponents of Magnetic Phase Transitions in CP^{N-1} Model at 2<d<4

Critical exponents in the CP^{N-1} model, which describes localized-moment ferro- and antiferromagnets (N=2 in the Heisenberg model), are calculated from two-particle Green's functions to first order in 1/N. For d=2+εthe results agree with earlier renormalization group calculations. For d=3 the leading 1/N-corrections turn out to be very large at N=2. For d=4-εthe 1/N-corrections are small at any N and insufficient to describe correctly the magnetic phase transition.

preprint1997arXiv

On the phonon-induced superconductivity of disordered alloys

A model of alloy is considered which includes both quenched disorder in the electron subsystem (``alloy'' subsystem) and electron-phonon interaction. For given approximate solution for the alloy part of the problem, which is assumed to be conserving in Baym's sense, we construct the generating functional and derive the Eliashberg-type equations which are valid to the lowest order in the adiabatic parameter.The renormalization of bare electron-phonon interaction vertices by disorder is taken into account consistently with the approximation for the alloy self-energy. For the case of exact configurational averaging the same set of equations is established within the usual T-matrix approach. We demonstrate that for any conserving approximation for the alloy part of the self-energy the Anderson's theorem holds in the case of isotropic singlet pairing provided disorder renormalizations of the electron-phonon interaction vertices are neglected. Taking account of the disorder renormalization of the electron-phonon interaction we analyze general equations qualitatively and present the expressions for $T_{c}$ for the case of weak and intermediate electron-phonon coupling. Disorder renormalizations of the logarithmic corrections to the effective coupling, which arise when the effective interaction kernel for the Cooper channel has the second energy scale, as well as the renormalization of the dilute paramagnetic impurity suppression are discussed.

preprint1997arXiv

On the Self-Consistent Spin-Wave Theory of Two-Dimensional Magnets With Impurities

The self-consistent spin-wave theory is applied to investigate the magnetization distribution around the impurity in isotropic and easy-axis two-dimensional ferro- and antiferromagents. The temperature dependences of host magnetization disturbance and impurity magnetization are calculated. The short-range order in the isotropic case is investigated. Importance of dynamic and kinematic interactions of spin waves is demonstrated.

preprint1997arXiv

Peculiarities of phonon spectra and lattice heat capacity in Ir and Rh

A simple pseudopotential model is proposed, which allows the phonon spectra and temperature dependence of the lattice heat capacity of Ir and Rh be described with a high enough accuracy. A careful comparison of the calculated and experimental values of the lattice heat capacity is carried out, with the procedure of the identification of the phonon contribution to the heat capacity and determination of the characteristics (momenta) of the phonon density of states from the experimental values of the total heat capacity of metal at a constant pressure being described in detail. The results of the theoretical calculations explain, in particular, such peculiar feature of Ir and Rh, unusual for cubic metals, as a sharp (more than by a factor of 1.5) decrease in the effective Debye temperature with increasing termperature. The temperature dependence of the mean square amplitude of atomic displacements in Ir and Rh has been calculated. Basing on the band calculations the manifestation of the Kohn singularities in the phonon spectra of Ir are discussed.

M. I. Katsnelson

What is connected

Connect this record

See the researcher in context

Building this map preview

207 published item(s)

Revised Tolmachev-Morel-Anderson pseudopotential for layered conventional superconductors with nonlocal Coulomb interaction

Spin dynamics of itinerant electrons: local magnetic moment formation and Berry phase

Certification of quantum states with hidden structure of their bitstrings

Coexisting charge density wave and ferromagnetic instabilities in monolayer InSe

Local Plaquette Physics as Key Ingredient of High-Temperature Superconductivity in Cuprates

From local to non-local correlations: the Dual Boson perspective

Limits on gas impermeability of graphene

Magnetic Two-Dimensional Chromium Trihalides: A Theoretical Perspective

Nonlinear optical study of commensurability effects in graphene-hBN heterostructures

Origin of the vortex displacement field in twisted bilayer graphene

Probing the topology of the quantum analog of a classical skyrmion

Temperature-Induced Lifshitz Transition and Possible Excitonic Instability in ZrSiSe

Uncommon clustering in dilute Ti-Fe alloys

Giant anisotropy of Gilbert damping in a Rashba honeycomb antiferromagnet

Magnetic polaron and antiferro-ferromagnetic transition in doped bilayer CrI$_3$

Orbitally-resolved ferromagnetism of monolayer CrI$_3$

Capturing non-local interaction effects in the Hubbard model: optimal mappings and limits of applicability

Chirality-dependent transmission of spin waves through domain walls

Correlation effects and orbital magnetism of Co clusters

Decoherence wave in magnetic systems and creation of Néel antiferromagnetic state by measurement

Friedel oscillations at the surfaces of rhombohedral $N$-layer graphene

Interaction-driven Lifshitz transition with dipolar fermions in optical lattices

Intrinsic Charge Carrier Mobility in Single-Layer Black Phosphorus

Logical inference approach to relativistic quantum mechanics: derivation of the Klein-Gordon equation

Macroscopic self-reorientation of interacting two-dimensional crystals

Many-body effects on Cr(001) surfaces: An LDA+DMFT study

Many-body effects on graphene conductivity: Quantum Monte Carlo calculations

Midgap states and band gap modification in defective graphene/h-BN heterostructures

Monte-Carlo study of Dirac semimetals phase diagram

Plaquette Valence Bond Theory of High-Temperature Superconductivity

Polaronic effects in monolayer black phosphorus on polar substrates

Quantum elasticity of graphene: Thermal expansion coefficient and specific heat

Quantum theory as plausible reasoning applied to data obtained by robust experiments

Role of direct exchange and Dzyaloshinskii-Moriya interactions in magnetic properties of graphene derivatives: C$_2$F and C$_2$H

Self-consistent Dual Boson approach to single-particle and collective excitations in correlated systems

Spin-orbit coupling and magnetic interactions in Si(111):{C,Si,Sn,Pb}

Standard model of the rare-earths, analyzed from the Hubbard I approximation

Treatment of 4f states of the rare-earths: the case study of TbN

Chemical modifications and stability of phosphorene with impurities: A first principles study

Chiral tunneling through generic one-dimensional potential barriers in bilayer graphene

Dirac electrons and domain walls: a realization in junctions of ferromagnets and topological insulators

Effect of structural relaxation on the electronic structure of graphene on hexagonal boron nitride

Effective Hamiltonians for fastly driven many-body lattice systems

Exchange parameters of strongly correlated materials: extraction from spin-polarised density functional theory plus dynamical mean field theory

Fingerprints of disorder in graphene

Geometric, electronic and magnetic structure of Fe$_{x}$O$_{y}^{+}$ clusters

Landau levels of single layer and bilayer phosphorene

Layer-resolved magnetic exchange interactions of surfaces of late 3d elements: effects of electronic correlations

Magnetism and interaction-induced gap opening in graphene with vacancies or hydrogen adatoms: Quantum Monte Carlo study

Many-body effects in graphene beyond the Dirac model with Coulomb interaction

Melting temperature of graphene

Microscopic origin of Heisenberg and non-Heisenberg exchange interactions in ferromagnetic bcc Fe

Misfit stabilized embedded nanoparticles in metallic alloys

Modeling Klein tunneling and caustics of electron waves in graphene

Moir{é} patterns as a probe of interplanar interactions: graphene on h-BN

Phonon-assisted resonant tunneling of electrons in graphene-boron nitride transistors

Quantum theory as a description of robust experiments: Application to Stern-Gerlach and Einstein-Podolsky-Rosen-Bohm experiments

Racah materials: role of atomic multiplets in intermediate valence systems

Relaxation of moiré patterns for slightly misaligned identical lattices: graphene on graphite

Scaling behavior and strain dependence of in-plane elastic properties of graphene

Toward a realistic description of multilayer black phosphorus: from $GW$ approximation to large-scale tight-binding simulations

Transport and Optical Properties of Single- and Bilayer Black Phosphorus with Defects

Commensurate-incommensurate transition for graphene on hexagonal boron nitride

Effect of Point Defects on the Optical and Transport Properties of MoS2 and WS2

Electron correlations in Mn$_x$Ga$_{1-x}$As as seen by resonant electron spectroscopy and dynamical mean field theory

Electronic correlations determine the phase stability of iron up to the melting temperature

First-principles modelling of magnetic excitations in Mn12

Inverse magneto-refraction as a mechanism for laser modification of spin-spin exchange parameters and subsequent terahertz emission from iron oxides

Metal-Insulator Transition in Graphene on Boron Nitride

Motion of domain walls and the dynamics of kinks in the magnetic Peierls potential

Plasmons in strongly correlated systems: spectral weight transfer and renormalized dispersion

Probing of valley polarization in graphene via optical second-harmonic generation

Proton transport through one atom thick crystals

Quantum theory as the most robust description of reproducible experiments