Phonon-assisted magnetic Mott-insulating state in the charge density wave phase of single-layer 1TNbSe2

We study the structural, electronic and vibrational properties of single-layer 1TNbSe$_2$ from first principles. Within the generalized gradient approximation, the 1T polytype is highly unstable with respect to the 2H. The DFT+U method improves the stability of the 1T phase, explaining its detection in experiments. A charge density wave occurs with a $\sqrt{13}\times\sqrt{13}~R30^{\circ}$ periodicity, in agreement with STM data. At $U=0$, the David-star reconstruction displays a flat band below the Fermi level with a marked d$_{z^2-r^2}$ orbital character of the central Nb. The Hubbard interaction induces a magnetic Mott insulating state. Magnetism distorts the lattice around the central Nb atom in the star, reduces the hybridization between the central Nb d$_{z^2-r^2}$ orbital and the neighbouring Se p-states and lifts in energy the flat band becoming non-bonding. This cooperative lattice and magnetic effect amplifies the Mott gap. Single-layer 1TNbSe$_2$ is then a phonon-assisted spin-$1/2$ Magnetic Mott insulator.Comment: 6 pages, 9 picture

Chemically-exfoliated single-layer MoS2_2 : stability, lattice dynamics and catalytic adsorption from first principles

Chemically and mechanically exfoliated MoS$_2$ single-layer samples have substantially different properties. While mechanically exfoliated single-layers are mono-phase (1H polytype with Mo in trigonal prismatic coordination), the chemically exfoliated samples show coexistence of three different phases, 1H, 1T (Mo in octahedral coordination) and 1T$^{'}$ (a distorted $2\times 1$ 1T-superstructure). By using first-principles calculations, we investigate the energetics and the dynamical stability of the three phases. We show that the 1H phase is the most stable one, while the metallic 1T phase, strongly unstable, undergoes a phase transition towards a metastable and insulating 1T$^{'}$ structure composed of separated zig-zag chains. We calculate electronic structure, phonon dispersion, Raman frequencies and intensities for the 1T$^{'}$ structure. We provide a microscopical description of the J$_1$, J$_2$ and J$_3$ Raman features first detected more then $20$ years ago, but unexplained up to now. Finally, we show that H adsorbates, that are naturally present at the end of the chemical exfoliation process, stabilize the 1T$^{\prime}$ over the 1H one.Comment: 7 Pages, 8 Pictures, To appear on Phys. Rev.

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

Charge density wave and spin 1/21/2 insulating state in single layer 1T-NbS2_2

In bulk samples and few layer flakes, the transition metal dichalcogenides NbS$_2$ and NbSe$_2$ assume the H polytype structure with trigonal prismatic coordination of the Nb atom. Recently, however, single and few layers of 1T-NbSe$_2$ with octahedral coordination around the transition metal ion were synthesized. Motivated by these experiments and by using first-principles calculations, we investigate the structural, electronic and dynamical properties of single layer 1T-NbS$_2$. We find that single-layer 1T-NbS$_2$ undergoes a $\sqrt{13}\times\sqrt{13}$ star-of-David charge density wave. Within the generalized gradient approximation, the weak interaction between the stars leads to an ultraflat band at the Fermi level isolated from all other bands. The spin-polarized generalized gradient approximation stabilizes a total spin $1/2$ magnetic state with opening of a $0.15$ eV band gap and a $0.21\mu_B$ magnetic moment localized on the central Nb in the star. Within GGA+U, the magnetic moment on the central Nb is enhanced to $0.41\mu_{B}$ and a larger gap occurs. Most important, this approximation gives a small energy difference between the 1T and 1H polytypes (only $+0.5$ mRy/Nb), suggesting that the 1T-polytype can be synthesized in a similar way as done for single layer 1T-NbSe$_2$. Finally we compute first and second nearest neighbors magnetic inter-star exchange interactions finding $J_1$=9.5~K and $J_2$=0.4~K ferromagnetic coupling constants

Projector augmented wave calculation of x-ray absorption spectra at the L2,3 edges

We develop a technique based on density functional theory and the projector augmented wave method in order to obtain the x-ray absorption cross section at a general edge, both in the electric dipole and quadrupole approximations. The method is a generalization of Taillefumier et al., PRB 66, 195107 (2002). We apply the method to the calculation of the Cu L2,3 edges in fcc copper and cuprite (Cu2O), and to the S L2,3 edges in molybdenite (2H-MoS2). The role of core-hole effects, modeled in a supercell approach, as well as the decomposition of the spectrum into different angular momentum channels are studied in detail. In copper we find that the best agreement with experimental data is obtained when core-hole effects are neglected. On the contrary, core-hole effects need to be included both in Cu2O and 2H-MoS2. Finally we show that a non-negligible component of S L2,3 edges in 2H-MoS2 involves transition to states with s character at all energy scales. The inclusion of this angular momentum channel is mandatory to correctly describe the angular dependence of the measured spectra. We believe that transitions to s character states are quantitatively significant at the L2,3 edges of third row elements from Al to Ar.Comment: 12 pages, 10 picture

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

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

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