Electronic structure of heavily-doped graphene: the role of foreign atom states

Using density functional theory calculations we investigate the electronic structure of graphene doped by deposition of foreign atoms. We demonstrate that, as the charge transfer to the graphene layer increases, the band structure of the pristine graphene sheet is substantially affected. This is particularly relevant when Ca atoms are deposed on graphene at CaC$_{6}$ stoichiometry. Similarly to what happens in superconducting graphite intercalated compounds, a Ca bands occurs at the Fermi level. Its hybridization with the C states generates a strong non-linearity in one of the $\pi^{*}$ bands below the Fermi level, at energies comparable to the graphene E$_{2g}$ phonon frequency. This strong non-linearity, and not manybody effects as previously proposed, explains the large and anisotropic values of the apparent electron-phonon coupling measured in angular resolved photoemission.Comment: 4 pages, 2 figures, see also M. Calandra and F. Mauri,arXiv:0707.146

Possibility of superconductivity in graphite intercalated with alkaline earths investigated with density functional theory

Using density functional theory we investigate the occurrence of superconductivity in AC$_6$ with A=Mg,Ca,Sr,Ba. We predict that at zero pressure, Ba and Sr should be superconducting with critical temperatures (T$_c$) 0.2 K and 3.0 K, respectively. We study the pressure dependence of T$_c$ assuming the same symmetry for the crystal structures at zero and finite pressures. We find that the SrC$_6$ and BaC$_6$ critical temperatures should be substantially enhanced by pressure. On the contrary, for CaC$_6$ we find that in the 0 to 5 GPa region, T$_c$ weakly increases with pressure. The increase is much smaller than what shown in several recent experiments. Thus we suggest that in CaC$_6$ a continous phase transformation, such as an increase in staging, occurs at finite pressure. Finally we argue that, although MgC$_6$ is unstable, the synthesis of intercalated systems of the kind Mg$_x$Ca$_{1-x}$C$_y$ could lead to higher critical temperatures.Comment: 9 page

Charge density wave and superconducting dome in TiSe2 from electron-phonon interaction

At low temperature TiSe2 undergoes a charge density wave instability. Superconductivity is stabilized either by pressure or by Cu intercalation. We show that the pressure phase diagram of TiSe2 is well described by first-principles calculations. At pressures smaller than 4 GPa charge density wave ordering occurs, in agreement with experiments. At larger pressures the disappearing of the charge density wave is due to a stiffening of the short-range force-constants and not to the variation of nesting with pressure. Finally we show that the behavior of Tc as a function of pressure is entirely determined by the electron-phonon interaction without need of invoking excitonic mechanisms. Our work demonstrates that phase-diagrams with competing orders and a superconducting dome are also obtained in the framework of the electron-phonon interaction.Comment: 4 pages, 7 picture

Anharmonic phonon frequency shift in MgB2

We compute the anharmonic shift of the phonon frequencies in MgB2, using density functional theory. We explicitly take into account the scattering between different phonon modes at different q-points in the Brillouin zone. The shift of the E2g mode at the Gamma point is +5 % of the harmonic frequency. This result comes from the cancellation between the contributions of the four- and three-phonon scattering, respectively +10 % and -5 %. A similar shift is predicted at the A point, in agreement with inelastic X-ray scattering phonon-dispersion measurements. A smaller shift is observed at the M point.Comment: 4 pages, 1 figur

Local and non-local electron-phonon couplings in K3Picene and the effect of metallic screening

We analyze the properties of electron-phonon couplings in K3Picene by exploiting a molecular orbital representation derived in the maximally localized Wannier function formalism. This allows us to go beyond the analysis done in Phys. Rev. Lett. 107, 137006 (2011), and separate not only the intra- and intermolecular phonon contributions but also the local and non-local electronic states in the electron-phonon matrix elements. Despite the molecular nature of the crystal, we find that the purely molecular contributions (Holstein-like couplings where the local deformation potential is coupled to intramolecular phonons) account for only 20% of the total electron-phonon interaction lambda. In particular, the Holstein-like contributions to lambda in K3Picene are four times smaller than those computed for an isolated neutral molecule, as they are strongly screened by the metallic bands of the doped crystal. Our findings invalidate the use of molecular electron-phonon calculations to estimate the total electron-phonon coupling in metallic picene, and possibly in other doped metallic molecular crystals. The major contribution (80%) to lambda in K3Picene comes from non-local couplings due to phonon modulated hoppings. We show that the crystal geometry together with the molecular picene structure leads to a strong 1D spatial anisotropy of the non-local couplings. Finally, based on the parameters derived from our density functional theory calculations, we propose a lattice modelization of the electron-phonon couplings in K3Picene which gives 90% of ab-initio lambda.Comment: 13 pages, 8 figures, 3 table

Electrochemical doping of few layer ZrNCl from first-principles: electronic and structural properties in field-effect configuration

We develop a first-principles theoretical approach to doping in field-effect devices. The method allows for calculation of the electronic structure as well as complete structural relaxation in field-effect configuration using density-functional theory. We apply our approach to ionic-liquid-based field-effect doping of monolayer, bilayer, and trilayer ZrNCl and analyze in detail the structural changes induced by the electric field. We show that, contrary to what is assumed in previous experimental works, only one ZrNCl layer is electrochemically doped and that this induces large structural changes within the layer. Surprisingly, despite these structural and electronic changes, the density of states at the Fermi energy is independent of the doping. Our findings imply a substantial revision of the phase diagram of electrochemically doped ZrNCl and elucidate crucial differences with superconductivity in Li intercalated bulk ZrNCl.Comment: 15 pages, 14 figure

Field-effect-driven half-metallic multilayer graphene

Rhombohedral stacked multilayer graphene displays the occurrence of a magnetic surface state at low temperatures. Recent angular resolved photoemission experiments demonstrate the robustness of the magnetic state in long sequences of ABC graphene. Here, by using first-principles calculations, we show that field-effect doping of these graphene multilayers induces a perfect half-metallic behaviour with 100% of spin current polarization already at dopings attainable in conventional field effect transistors with solid state dielectrics. Our work demonstrates the realisability of a new kind of spintronic devices where the transition between the low resistance and the high resistance state is driven only by electric fields

Anharmonic free energies and phonon dispersions from the stochastic self-consistent harmonic approximation: application to platinum and palladium hydrides

Harmonic calculations based on density-functional theory are generally the method of choice for the description of phonon spectra of metals and insulators. The inclusion of anharmonic effects is, however, delicate as it relies on perturbation theory requiring a considerable amount of computer time, fast increasing with the cell size. Furthermore, perturbation theory breaks down when the harmonic solution is dynamically unstable or the anharmonic correction of the phonon energies is larger than the harmonic frequencies themselves.We present a stochastic implementation of the self-consistent harmonic approximation valid to treat anharmonicity at any temperature in the non-perturbative regime. The method is based on the minimization of the free energy with respect to a trial density matrix described by an arbitrary harmonic Hamiltonian. The minimization is performed with respect to all the free parameters in the trial harmonic Hamiltonian, namely, equilibrium positions, phonon frequencies and polarization vectors. The gradient of the free energy is calculated following a stochastic procedure. The method can be used to calculate thermodynamic properties, dynamical properties and anharmonic corrections to the Eliashberg function of the electron-phonon coupling. The scaling with the system size is greatly improved with respect to perturbation theory. The validity of the method is demonstrated in the strongly anharmonic palladium and platinum hydrides. In both cases we predict a strong anharmonic correction to the harmonic phonon spectra, far beyond the perturbative limit. In palladium hydrides we calculate thermodynamic properties beyond the quasiharmonic approximation, while in PtH we demonstrate that the high superconducting critical temperatures at 100 GPa predicted in previous calculations based on the harmonic approximation are strongly suppressed when anharmonic effects are included.Comment: 17 pages, 5 figure

Density-functional calculation of static screening in 2D materials: the long-wavelength dielectric function of graphene

We calculate the long-wavelength static screening properties of both neutral and doped graphene in the framework of density-functional theory. We use a plane-wave approach with periodic images in the third dimension and truncate the Coulomb interactions to eliminate spurious interlayer screening. We carefully address the issue of extracting two dimensional dielectric properties from simulated three-dimensional potentials. We compare this method with analytical expressions derived for two dimensional massless Dirac fermions in the random phase approximation. We evaluate the contributions of the deviation from conical bands, exchange-correlation and local-fields. For momenta smaller than twice the Fermi wavevector, the static screening of graphene within the density-functional perturbative approach agrees with the results for conical bands within random phase approximation and neglecting local fields. For larger momenta, we find that the analytical model underestimates the static dielectric function by $\approx 10$, mainly due to the conical band approximation