Source author record

Nicola Marzari

Nicola Marzari appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.mtrl-sci cond-mat.mes-hall physics.comp-ph physics.chem-ph cond-mat.str-el cond-mat.soft cond-mat.stat-mech physics.app-ph quant-ph Software Engineering

Catalog footprint

What is connected

61works

10topics

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.

preprint2026arXiv

Novel fast Li-ion conductors for solid-state electrolytes from first-principles

We present a high-throughput computational screening for fast lithium-ion conductors to identify promising materials for application in all solid-state electrolytes. Starting from more than 30,000 Li-containing experimental structures sourced from Crystallography Open Database, Inorganic Crystal Structure Database and Materials Platform for Data Science, we perform highly automated calculations to identify electronic insulators. On these ~1000 structures, we use molecular dynamics simulations to estimate Li-ion diffusivities using the pinball model, which describes the potential energy landscape of diffusing lithium with accuracy similar to density functional theory while being 200-500 times faster. Then we study the ~60 most promising and previously unknown fast conductors with full first-principles molecular dynamics simulations at several temperatures to estimate their activation barriers. The results are discussed in detail for the 9 fastest conductors, including $Li_7NbO_6$ which shows a remarkable ionic conductivity of ~5 mS/cm at room temperature. We further present the entire screening protocol, including the workflows where the accuracy of the pinball model is improved self-consistently, necessary to automatically running the required calculations and analysing their results.

preprint2025arXiv

Extreme anharmonicity and thermal contraction of 1D wires

Ultrathin nanowires could play a central role in next-generation downscaled electronics. Here, we explore some of the most promising candidates identified from previous high-throughput screening: CuC$_2$, TaSe$_3$, and AuSe$_2$, to gain insight into the thermodynamic and anharmonic behaviors of nanowires that could be exfoliated from weakly-bonded three-dimensional materials. We analyze thermal stability, linear thermal expansion, and anharmonic heat capacity using the stochastic self-consistent harmonic approximation. Notably, our work unveils exotic features common among all the 1D wires: a colossal record negative thermal expansion and very large deviations from the Dulong-Petit law due to strong anharmonicity.

preprint2022arXiv

Bloch's theorem in orbital-density-dependent functionals: Band structures from Koopmans spectral functionals

Koopmans-compliant functionals provide an orbital-density-dependent framework for an accurate evaluation of spectral properties; they are obtained by imposing a generalized piecewise-linearity condition on the total energy of the system with respect to the occupation of any orbital. In crystalline materials, due to the orbital-density-dependent nature of the functionals, minimization of the total energy to a ground state provides a set of minimizing variational orbitals that are localized and thus break the periodicity of the underlying lattice. Despite this, we show that Bloch symmetry can be preserved and it is possible to describe the electronic states with a band-structure picture, thanks to the Wannier-like character of the variational orbitals. We also present a method to unfold and interpolate the electronic bands from supercell ($Γ$-point) calculations, which enables us to calculate full band structures with Koopmans-compliant functionals. The results obtained for a set of benchmark semiconductors and insulators show very good agreement with state-of-the-art many-body perturbation theory and experiments, underscoring the reliability of these spectral functionals in predicting band structures.

preprint2022arXiv

HP -- A code for the calculation of Hubbard parameters using density-functional perturbation theory

We introduce HP, an implementation of density-functional perturbation theory, designed to compute Hubbard parameters (on-site $U$ and inter-site $V$) in the framework of DFT+$U$ and DFT+$U$+$V$. The code does not require the use of computationally expensive supercells of the traditional linear-response approach; instead, unit cells are used with monochromatic perturbations that significantly reduce the computational cost of determining Hubbard parameters. HP is an open-source software distributed under the terms of the GPL as a component of Quantum ESPRESSO. As with other components, HP is optimized to run on a variety of different platforms, from laptops to massively parallel architectures, using native mathematical libraries (LAPACK and FFTW) and a hierarchy of custom parallelization layers built on top of MPI. The effectiveness of the code is showcased by computing Hubbard parameters self-consistently for the phospho-olivine Li$_x$Mn$_{1/2}$Fe$_{1/2}$PO$_4$ ($x=0, 1/2, 1$) and by highlighting the accuracy of predictions of the geometry and Li intercalation voltages.

preprint2022arXiv

Infrared-active phonons in one-dimensional materials and their spectroscopic signatures

Dimensionality provides a clear fingerprint on the dispersion of infrared-active, polar-optical phonons. For these phonons, the local dipoles parametrized by the Born effective charges drive the LO-TO splitting of bulk materials; this splitting actually breaks down in two-dimensional materials. Here, we extend the existing theory to the one-dimensional (1D) case. Combining an analytical model with the implementation of density-functional perturbation theory in 1D boundary conditions, we show that the dielectric splitting in the dispersion relations collapses logarithmically at the zone center. The dielectric properties and the radius of the 1D materials are linked by the present work to these red shifts, opening novel IR and Raman characterization avenues.

preprint2022arXiv

Koopmans spectral functionals in periodic-boundary conditions

Koopmans spectral functionals aim to describe simultaneously ground state properties and charged excitations of atoms, molecules, nanostructures and periodic crystals. This is achieved by augmenting standard density functionals with simple but physically motivated orbital-density-dependent corrections. These corrections act on a set of localized orbitals that, in periodic systems, resemble maximally localized Wannier functions. At variance with the original, direct supercell implementation [Phys. Rev. X 8, 021051 (2018)], we discuss here i) the complex but efficient formalism required for a periodic-boundary code using explicit Brillouin zone sampling, and ii) the calculation of the screened Koopmans corrections with density-functional perturbation theory. In addition to delivering improved scaling with system size, the present development makes the calculation of band structures with Koopmans functionals straightforward. The implementation in the open-source Quantum ESPRESSO distribution and the application to prototypical insulating and semiconducting systems are presented and discussed.

preprint2022arXiv

Many-body Green's function approach to lattice thermal transport

Recent progress in understanding thermal transport in complex crystals has highlighted the prominent role of heat conduction mediated by interband tunneling processes, which emerge between overlapping phonon bands (i.e. with energy differences smaller than their broadenings). These processes have recently been described in different ways, relying on the Wigner or Green-Kubo formalism, leading to apparently different results which question the definition of the heat-current operator. Here, we implement a full quantum approach based on the Kubo formula, elucidating analogies and differences with the recently introduced Wigner or Green-Kubo formulations, and extending the description of thermal transport to the overdamped regime of atomic vibrations, where the phonon quasiparticle picture breaks down. We rely on first-principles calculations on complex crystals with ultralow conductivity to compare numerically the thermal conductivity obtained within the aforementioned approaches, showing that at least in the quasiparticle regime the differences are negligible for practical applications.

preprint2022arXiv

Microscopic picture of paraelectric perovskites from structural prototypes

We show with first-principles molecular dynamics the persistence of intrinsic $\langle111\rangle$ Ti off-centerings for BaTiO$_3$ in its cubic paraelectric phase. Intriguingly, these are inconsistent with the Pm$\bar 3$m space group often used to atomistically model this phase using density functional theory or similar methods. Therefore we deploy a systematic symmetry analysis to construct representative structural models in the form of supercells that satisfy a desired point symmetry but are built from the combination of lower-symmetry primitive cells. We define as structural prototypes the smallest of these that are both energetically and dynamically stable. Remarkably, two 40-atom prototypes can be identified for paraelectric BaTiO$_3$; these are also common to many other $AB$O$_3$ perovskites. These prototypes can offer structural models of paraelectric phases that can be used for the computational engineering of functional materials. Last, we show that the emergence of B-cation off-centerings and the primitive-cell phonon instabilities are controlled by the equilibrium volume, in turn dictated by the filler A cation.

preprint2022arXiv

Phonon-assisted luminescence in defect centers from many-body perturbation theory

Phonon-assisted luminescence is a key property of defect centers in semiconductors, and can be measured to perform the readout of the information stored in a quantum bit, or to detect temperature variations. The investigation of phonon-assisted luminescence usually employs phenomenological models, such as that of Huang and Rhys, with restrictive assumptions that can fail to be predictive. In this work, we predict luminescence and study exciton-phonon couplings within a rigorous many-body perturbation theory framework, an analysis that has never been performed for defect centers. In particular, we study the optical emission of the negatively-charged boron vacancy in 2D hexagonal boron nitride, which currently stands out among defect centers in 2D materials thanks to its promise for applications in quantum information and quantum sensing. We show that phonons are responsible for the observed luminescence, which otherwise would be dark due to symmetry. We also show that the symmetry breaking induced by the static Jahn-Teller effect is not able to describe the presence of the experimentally observed peak at 1.5 eV.

preprint2022arXiv

Pulay forces in density-functional theory with extended Hubbard functionals: From nonorthogonalized to orthogonalized manifolds

We present a derivation of the exact expression for Pulay forces in density-functional theory calculations augmented with extended Hubbard functionals, and arising from the use of orthogonalized atomic orbitals as projectors for the Hubbard manifold. The derivative of the inverse square root of the orbital overlap matrix is obtained as a closed-form solution of the associated Lyapunov (Sylvester) equation. The expression for the resulting contribution to the forces is presented in the framework of ultrasoft pseudopotentials and the projector-augmented-wave method, and using a plane wave basis set. We have benchmarked the present implementation with respect to finite differences of total energies for the case of NiO, finding excellent agreement. Owing to the accuracy of Hubbard-corrected density-functional theory calculations - provided the Hubbard parameters are computed for the manifold under consideration - the present work paves the way for systematic studies of solid-state and molecular transition-metal and rare-earth compounds.

preprint2022arXiv

Wigner formulation of thermal transport in solids

Two different heat-transport mechanisms are discussed in solids: in crystals, heat carriers propagate and scatter particle-like as described by Peierls' formulation of the Boltzmann transport equation for phonon wavepackets. In glasses, instead, carriers behave wave-like, diffusing via a Zener-like tunneling between quasi-degenerate vibrational eigenstates, as described by the Allen-Feldman equation. Recently, it has been shown that these two conduction mechanisms emerge from a Wigner transport equation, which unifies and extends the Peierls-Boltzmann and Allen-Feldman formulations, allowing to describe also complex crystals where particle-like and wave-like conduction mechanisms coexist. Here, we discuss the theoretical foundations of such transport equation as is derived from the Wigner phase-space formulation of quantum mechanics, elucidating how the interplay between disorder, anharmonicity, and the quantum Bose-Einstein statistics of atomic vibrations determines thermal conductivity. This Wigner formulation argues for a preferential phase convention for the dynamical matrix in the reciprocal Bloch representation and related off-diagonal velocity operator's elements; such convention is the only one yielding a conductivity which is invariant with respect to the non-unique choice of the crystal's unit cell and is size-consistent. We rationalize the conditions determining the crossover from particle-like to wave-like heat conduction, showing that phonons below the Ioffe-Regel limit (i.e. with a mean free path shorter than the interatomic spacing) contribute to heat transport due to their wave-like capability to interfere and tunnel. Finally, we show that the present approach overcomes the failures of the Peierls-Boltzmann formulation for materials with ultralow or glass-like thermal conductivity.

preprint2021arXiv

Self-consistent Hubbard parameters from density-functional perturbation theory in the ultrasoft and projector-augmented wave formulations

The self-consistent evaluation of Hubbard parameters using linear-response theory is crucial for quantitatively predictive calculations based on Hubbard-corrected density-functional theory. Here, we extend a recently-introduced approach based on density-functional perturbation theory (DFPT) for the calculation of the on-site Hubbard $U$ to also compute the inter-site Hubbard $V$. DFPT allows to reduce significantly computational costs, improve numerical accuracy, and fully automate the calculation of the Hubbard parameters by recasting the linear response of a localized perturbation into an array of monochromatic perturbations that can be calculated in the primitive cell. In addition, here we generalize the entire formalism from norm-conserving to ultrasoft and projector-augmented wave formulations, and to metallic ground states. After benchmarking DFPT against the conventional real-space Hubbard linear response in a supercell, we demonstrate the effectiveness of the present extended Hubbard formulation in determining the equilibrium crystal structure of Li$_x$MnPO$_4$ (x=0,1) and the subtle energetics of Li intercalation.

preprint2020arXiv

2-D materials for ultra-scaled field-effect transistors: hundred candidates under the ab initio microscope

Thanks to their unique properties single-layer 2-D materials appear as excellent candidates to extend Moore's scaling law beyond the currently manufactured silicon FinFETs. However, the known 2-D semiconducting components, essentially transition metal dichalcogenides, are still far from delivering the expected performance. Based on a recent theoretical study that predicts the existence of more than 1,800 exfoliable 2-D materials, we investigate here the 100 most promising contenders for logic applications. Their "current vs. voltage" characteristics are simulated from first-principles, combining density-functional theory and advanced quantum transport calculations. Both n- and p-type configurations are considered, with gate lengths ranging from 15 down to 5 nm. From this unprecedented collection of electronic materials, we identify 13 compounds with electron and hole currents potentially much higher than in future Si FinFETs. The resulting database widely expands the design space of 2-D transistors and provides original guidelines to the materials and device engineering community.

preprint2020arXiv

Automated high-throughput Wannierisation

Maximally-localised Wannier functions (MLWFs) are routinely used to compute from first-principles advanced materials properties that require very dense Brillouin zone integration and to build accurate tight-binding models for scale-bridging simulations. At the same time, high-throughput (HT) computational materials design is an emergent field that promises to accelerate the reliable and cost-effective design and optimisation of new materials with target properties. The use of MLWFs in HT workflows has been hampered by the fact that generating MLWFs automatically and robustly without any user intervention and for arbitrary materials is, in general, very challenging. We address this problem directly by proposing a procedure for automatically generating MLWFs for HT frameworks. Our approach is based on the selected columns of the density matrix method (SCDM) and we present the details of its implementation in an AiiDA workflow. We apply our approach to a dataset of 200 bulk crystalline materials that span a wide structural and chemical space. We assess the quality of our MLWFs in terms of the accuracy of the band-structure interpolation that they provide as compared to the band-structure obtained via full first-principles calculations. Finally, we provide a downloadable virtual machine that can be used to reproduce the results of this paper, including all first-principles and atomistic simulations as well as the computational workflows.

preprint2020arXiv

Electronic structure of pristine and Ni-substituted LaFeO$_3$ from near edge x-ray absorption fine structure experiments and first-principles simulations

We present a joint theoretical and experimental study of the oxygen $K$-edge spectra for LaFeO$_3$ and homovalent Ni-substituted LaFeO$_3$ (LaFe$_{0.75}$Ni$_{0.25}$O$_3$), using first-principles simulations based on density-functional theory with extended Hubbard functionals and x-ray absorption near edge structure (XANES) measurements. Ground-state and excited-state XANES calculations employ Hubbard on-site $U$ and inter-site $V$ parameters determined from first principles and the Lanczos recursive method to obtain absorption cross sections, which allows for a reliable description of XANES spectra in transition-metal compounds in a very broad energy range, with an accuracy comparable to that of hybrid functionals but at a substantially lower cost. We show that standard gradient-corrected exchange-correlation functionals fail in capturing accurately the electronic properties of both materials. In particular, for LaFe$_{0.75}$Ni$_{0.25}$O$_3$ they do not reproduce its semiconducting behaviour and provide a poor description of the pre-edge features at the O $K$ edge. The inclusion of Hubbard interactions leads to a drastic improvement, accounting for the semiconducting ground state of LaFe$_{0.75}$Ni$_{0.25}$O$_3$ and for a good agreement between calculated and measured XANES spectra. We show that the partial substitution of Fe for Ni affects the conduction-band bottom by generating a strongly hybridized O($2p$)-Ni($3d$) minority-spin empty electronic state. The present work, based on a consistent correction of self-interaction errors, outlines the crucial role of extended Hubbard functionals to describe the electronic structure of complex transition-metal oxides such as LaFeO$_3$ and LaFe$_{0.75}$Ni$_{0.25}$O$_3$ and paves the way to future studies on similar systems.

preprint2020arXiv

Equipartition of Energy Defines the Size-Thickness Relationship in Liquid-Exfoliated Nanosheets

Liquid phase exfoliation is a commonly used method to produce 2D nanosheets from a range of layered crystals. However, such nanosheets display broad size and thickness distributions and correlations between area and thickness, issues that limit nanosheet application potential. To understand the factors controlling the exfoliation process, we have liquid-exfoliated 11 different layered materials, size-selecting each into fractions before using AFM to measure the nanosheet length, width, and thickness distributions for each fraction. The resultant data show a clear power-law scaling of nanosheet area with thickness for each material. We have developed a simple nonequilibrium thermodynamics-based model predicting that the power-law prefactor is proportional to both the ratios of in-plane-tearing/out-of-plane-peeling energies and in-plane/out-of-plane moduli. By comparing the experimental data with the modulus ratio calculated from first-principles, we find close agreement between experiment and theory. This supports our hypothesis that energy equipartition holds between nanosheet tearing and peeling during sonication-assisted exfoliation.

preprint2020arXiv

Hubbard-corrected density functional perturbation theory with ultrasoft pseudopotentials

We present in full detail a newly developed formalism enabling density functional perturbation theory (DFPT) calculations from a DFT+$U$ ground state. The implementation includes ultrasoft pseudopotentials and is valid for both insulating and metallic systems. It aims at fully exploiting the versatility of DFPT combined with the low-cost DFT+$U$ functional. This allows to avoid computationally intensive frozen-phonon calculations when DFT+$U$ is used to eliminate the residual electronic self-interaction from approximate functionals and to capture the localization of valence electrons e.g. on $d$ or $f$ states. In this way, the effects of electronic localization (possibly due to correlations) are consistently taken into account in the calculation of specific phonon modes, Born effective charges, dielectric tensors and in quantities requiring well converged sums over many phonon frequencies, as phonon density of states and free energies. The new computational tool is applied to two representative systems, namely CoO, a prototypical transition metal monoxide and LiCoO$_2$, a material employed for the cathode of Li-ion batteries. The results show the effectiveness of our formalism to capture in a quantitatively reliable way the vibrational properties of systems with localized valence electrons.

preprint2020arXiv

Self-consistent DFT+$U$+$V$ study of oxygen vacancies in SrTiO$_3$

Contradictory theoretical results for oxygen vacancies in SrTiO$_3$ (STO) were often related to the peculiar properties of STO, which is a $d^0$ transition metal oxide with mixed ionic-covalent bonding. Here, we apply, for the first time, density functional theory (DFT) within the extended Hubbard DFT+$U$+$V$ approach, including on-site as well as inter-site electronic interactions, to study oxygen-deficient STO with Hubbard $U$ and $V$ parameters computed self-consistently via density-functional perturbation theory. Our results demonstrate that the extended Hubbard functional is a promising approach to study defects in materials with electronic properties similar to STO. Indeed, DFT+$U$+$V$ provides a better description of stoichiometric STO compared to standard DFT or DFT+$U$, the band gap and crystal field splitting being in good agreement with experiments. In turn, also the description of the electronic properties of oxygen vacancies in STO is improved, with formation energies in excellent agreement with experiments as well as results obtained with the most frequently used hybrid functionals, however at a fraction of the computational cost. While our results do not fully resolve the contradictory findings reported in literature, our systematic approach leads to a deeper understanding of their origin, which stems from different cell sizes, STO phases, the exchange-correlation functional, and the treatment of structural relaxations and spin-polarization.

preprint2020arXiv

Two-dimensional materials from high-throughput computational exfoliation of experimentally known compounds

We search for novel two-dimensional materials that can be easily exfoliated from their parent compounds. Starting from 108423 unique, experimentally known three-dimensional compounds we identify a subset of 5619 that appear layered according to robust geometric and bonding criteria. High-throughput calculations using van-der-Waals density-functional theory, validated against experimental structural data and calculated random-phase-approximation binding energies, allow to identify 1825 compounds that are either easily or potentially exfoliable, including all that are commonly exfoliated experimentally. In particular, the subset of 1036 easily exfoliable cases---layered materials held together mostly by dispersion interactions and with binding energies up to $30-35$ meV$\cdot\textÅ^{-2}$---provides a wealth of novel structural prototypes and simple ternary compounds, and a large portfolio to search materials for optimal properties. For the 258 compounds with up to 6 atoms per primitive cell we comprehensively explore vibrational, electronic, magnetic, and topological properties, identifying in particular 56 ferromagnetic and antiferromagnetic systems, including half-metals and half-semiconductors.

preprint2019arXiv

Abundance of $\mathbb{Z}_2$ topological order in exfoliable two-dimensional insulators

Quantum spin Hall insulators are a class of two-dimensional materials with a finite electronic band gap in the bulk and gapless helical edge states. In the presence of time-reversal symmetry, $\mathbb{Z}_2$ topological order distinguishes the topological phase from the ordinary insulating one. Some of the phenomena that can be hosted in these materials, from one-dimensional low-dissipation electronic transport to spin filtering, could be very promising for many technological applications in the fields of electronics, spintronics and topological quantum computing. Nevertheless, the rarity of two-dimensional materials that can exhibit non-trivial $\mathbb{Z}_2$ topological order at room temperature hinders development. Here, we screen a comprehensive database we recently created of 1825 monolayers that can be exfoliated from experimentally known compounds, to search for novel quantum spin Hall insulators. Using density-functional and many-body perturbation theory simulations, we identify 13 monolayers that are candidates for quantum spin Hall insulators, including high-performing materials such as AsCuLi$_2$ and jacutingaite (Pt$_2$HgSe$_3$). We also identify monolayer Pd$_2$HgSe$_3$ as a novel Kane-Mele quantum spin Hall insulator, and compare it with jacutingaite. Such a handful of promising materials are mechanically stable and exhibit $\mathbb{Z}_2$ topological order, either unpertubed or driven by a small amount of strain. Such screening highlights a relative abundance of $\mathbb{Z}_2$ topological order of around 1%, and provides an optimal set of candidates for experimental efforts.

preprint2019arXiv

Emergent dual topology in the three-dimensional Kane-Mele Pt$_2$HgSe$_3$

Recently, the very first large-gap Kane-Mele quantum spin Hall insulator was predicted to be monolayer jacutingaite (Pt$_2$HgSe$_3$), a naturally-occurring exfoliable mineral discovered in Brazil in 2008. The stacking of quantum spin Hall monolayers into a van-der-Waals layered crystal typically leads to a (0;001) weak topological phase, which does not protect the existence of surface states on the (001) surface. Unexpectedly, recent angle-resolved photoemission spectroscopy experiments revealed the presence of surface states dispersing over large areas of the 001-surface Brillouin zone of jacutingaite single crystals. The 001-surface states have been shown to be topologically protected by a mirror Chern number $C_M=-2$, associated with a nodal line gapped by spin-orbit interactions. Here, we extend the two-dimensional Kane-Mele model to bulk jacutingaite and unveil the microscopic origin of the gapped nodal line and the emerging crystalline topological order. By using maximally-localized Wannier functions, we identify a large non-trivial second nearest-layer hopping term that breaks the standard paradigm of weak topological insulators. Complemented by this term, the predictions of the Kane-Mele model are in remarkable agreement with recent experiments and first-principles simulations, providing an appealing conceptual framework also relevant for other layered materials made of stacked honeycomb lattices.

preprint2019arXiv

Generalization of Fourier's law into viscous heat equations

Heat conduction in dielectric crystals originates from the propagation of atomic vibrations, whose microscopic dynamics is well described by the linearized phonon Boltzmann transport equation. Recently, it was shown that thermal conductivity can be resolved exactly and in a closed form as a sum over relaxons, $\mathit{i.e.}$ collective phonon excitations that are the eigenvectors of Boltzmann equation's scattering matrix [Cepellotti and Marzari, PRX $\mathbf{6}$ (2016)]. Relaxons have a well-defined parity, and only odd relaxons contribute to the thermal conductivity. Here, we show that the complementary set of even relaxons determines another quantity --- the thermal viscosity --- that enters into the description of heat transport, and is especially relevant in the hydrodynamic regime, where dissipation of crystal momentum by Umklapp scattering phases out. We also show how the thermal conductivity and viscosity parametrize two novel viscous heat equations --- two coupled equations for the temperature and drift-velocity fields --- which represent the thermal counterpart of the Navier-Stokes equations of hydrodynamics in the linear, laminar regime. These viscous heat equations are derived from a coarse-graining of the linearized Boltzmann transport equation for phonons, and encompass both limits of Fourier's law and of second sound, taking place, respectively, in the regimes of strong or weak momentum dissipation. Last, we introduce the Fourier deviation number as a descriptor that captures the deviations from Fourier's law due to hydrodynamic effects. We showcase these findings in a test case of a complex-shaped device made of graphite, obtaining a remarkable agreement with the very recent experimental demonstration of hydrodynamic transport in this material. The present findings also suggest that hydrodynamic behavior can appear at room temperature in micrometer-sized diamond crystals.

preprint2019arXiv

Intrinsic edge excitons in two-dimensional MoS$_2$

Using accurate first-principles calculations based on many-body perturbation theory we predict that two-dimensional MoS$_2$ hosts edge excitons with universal character, intrinsic to the existence of edges and lying well below the onset of bulk features. These excitons are largely insensitive to edge terminations or orientation, persisting even in the presence of metallic screening at zigzag edges, with large binding energies of $\sim$0.4 eV. Additional excitons can also emerge in ultranarrow ribbons, or as a function of the chemical nature of the termination. The chemical, structural, and electronic similarities with Se- or W-based transition-metal dichalcogenides suggest that these optical features could be common in this class of materials.

preprint2019arXiv

Koopmans Meets Bethe-Salpeter: Excitonic Optical Spectra without GW

The Bethe-Salpeter Equation (BSE) can be applied to compute from first-principles optical spectra that include the effects of screened electron-hole interactions. As input, BSE calculations require single-particle states, quasiparticle energy levels and the screened Coulomb interaction, which are typically obtained with many-body perturbation theory, whose cost limits the scope of possible applications. This work tries to address this practical limitation, instead deriving spectral energies from Koopmans-compliant functionals and introducing a new methodology for handling the screened Coulomb interaction. The explicit calculation of the $W$ matrix is bypassed via a direct minimization scheme applied on top of a maximally localised Wannier function basis. We validate and benchmark this approach by computing the low-lying excited states of the molecules in Thiel's set, and the optical absorption spectrum of a $\text{C}_{60}$ fullerene. The results show the same trends as quantum chemical methods and are in excellent agreement with previous simulations carried out at the TD-DFT or $G_{0}W_{0}$-\text{BSE} level. Conveniently, the new framework reduces the parameter space controlling the accuracy of the calculation, thereby simplifying the simulation of charge-neutral excitations, offering the potential to expand the applicability of first-principles spectroscopies to larger systems of applied interest.

preprint2019arXiv

Operando XANES from first-principles and its application to iridium oxide

Efficient electro-catalytic water-splitting technologies require suitable catalysts for the oxygen evolution reaction (OER). The development of novel catalysts could benefit from the achievement of a complete understanding of the reaction mechanism on iridium oxide (IrO$_2$), an active catalyst material that is, however, too scarce for large-scale applications. Considerable insight has already been provided by \emph{operando} X-ray absorption near-edge structure (XANES) experiments, which paved the way towards an atomistic description of the catalyst's evolution in a working environment. We combine here first-principles simulations augmented with a continuum description of the solvent and electrolyte to investigate the electrochemical stability of various IrO$_2$ interfaces and to predict the XANES cross-section for selected terminations under realistic conditions of applied potential. The comparison of computed O K-edge XANES spectra to corresponding experiments supports the formation of electron-deficient surface oxygen species in the OER-relevant voltage regime. Furthermore, surface hydroxyl groups that are found to be stable up to $\sim$1 V are suggested to be progressively oxidized at larger potentials, giving rise to a shift in the Ir L$_3$-edge cross-section that qualitatively agrees with measurements.

preprint2019arXiv

Photorealistic modelling of metals from first principles

The colours of metals have attracted the attention of humanity since ancient times, and coloured metals, in particular gold compounds, have been employed for tools and objects symbolizing the aesthetics of power. In this work we develop a comprehensive framework to obtain the reflectivity and colour of metals, and show that the trends in optical properties and the colours can be predicted by straightforward first-principles techniques based on standard approximations. We apply this to predict reflectivity and colour of several elemental metals and of different types of metallic compounds (intermetallics, solid solutions and heterogeneous alloys), considering mainly binary alloys based on noble metals. We validate the numerical approach through an extensive comparison with experimental data and the photorealistic rendering of known coloured metals.

preprint2019arXiv

Wannier90 as a community code: new features and applications

Wannier90 is an open-source computer program for calculating maximally-localised Wannier functions (MLWFs) from a set of Bloch states. It is interfaced to many widely used electronic-structure codes thanks to its independence from the basis sets representing these Bloch states. In the past few years the development of Wannier90 has transitioned to a community-driven model; this has resulted in a number of new developments that have been recently released in Wannier90 v3.0. In this article we describe these new functionalities, that include the implementation of new features for wannierisation and disentanglement (symmetry-adapted Wannier functions, selectively-localised Wannier functions, selected columns of the density matrix) and the ability to calculate new properties (shift currents and Berry-curvature dipole, and a new interface to many-body perturbation theory); performance improvements, including parallelisation of the core code; enhancements in functionality (support for spinor-valued Wannier functions, more accurate methods to interpolate quantities in the Brillouin zone); improved usability (improved plotting routines, integration with high-throughput automation frameworks), as well as the implementation of modern software engineering practices (unit testing, continuous integration, and automatic source-code documentation). These new features, capabilities, and code development model aim to further sustain and expand the community uptake and range of applicability, that nowadays spans complex and accurate dielectric, electronic, magnetic, optical, topological and transport properties of materials.

preprint2016arXiv

Band-like Electron Transport with Record-High Mobility in the TCNQ family

In highest quality organic single-crystal field-effect transistors, electron transport occurs in the band-like regime, with the carrier mobility increasing upon lowering temperature. Neither the microscopic nature of this regime, nor why it occurs only in a small number of materials is currently understood. Here, comparative studies of closely related materials, exhibiting high-quality reproducible transport properties are needed. We performed a study of electron transport in single-crystals of different TCNQ (tetracyanoquinodimethane) molecules, combined with band structure calculations. We show that F2-TCNQ devices exhibit very high electron mobility and an unprecedented increase in mobility upon cooling, whereas in TCNQ and F4-TCNQ the mobility is substantially lower and decreases upon cooling. We analyze the crystal and electronic structures of these materials and find that F2-TCNQ crystals are indeed ideal to achieve outstanding transport properties. Our analysis also shows that to understand the difference between the three materials, studying their band structure is not sufficient, and that the electron-phonon coupling needs to be investigated as well. Besides the outstanding transport properties of F2-TCNQ, a key result of our work is the identification of the Fx-TCNQ family as a paradigm to investigate the most fundamental aspects of electronic transport in organic crystals.

preprint2016arXiv

Performance of arsenene and antimonene double-gate MOSFETs from first principles

In the race towards high-performance ultra-scaled devices, two-dimensional materials offer an alternative paradigm thanks to their atomic thickness suppressing short-channel effects. It is thus urgent to study the most promising candidates in realistic configurations, and here we present detailed multiscale simulations of field-effect transistors based on arsenene and antimonene monolayers as channels. The accuracy of first-principles approaches in describing electronic properties is combined with the efficiency of tight-binding Hamiltonians based on maximally-localised Wannier functions to compute the transport properties of the devices. These simulations provide for the first time estimates on the upper limits for the electron and hole mobilities in the Takagi's approximation, including spin-orbit and multi-valley effects, and demonstrate that ultra-scaled devices in the sub-10 nm scale show a performance that is compliant with industry requirements.

preprint2016arXiv

The electronic thermal conductivity of graphene

Graphene, as a semimetal with the largest known thermal conductivity, is an ideal system to study the interplay between electronic and lattice contributions to thermal transport. While the total electrical and thermal conductivity have been extensively investigated, a detailed first-principles study of its electronic thermal conductivity is still missing. Here, we first characterize the electron-phonon intrinsic contribution to the electronic thermal resistivity of graphene as a function of doping using electronic and phonon dispersions and electron-phonon couplings calculated from first principles at the level of density-functional theory and many-body perturbation theory (GW). Then, we include extrinsic electron-impurity scattering using low-temperature experimental estimates. Under these conditions, we find that the in-plane electronic thermal conductivity of doped graphene is ~300 W/mK at room temperature, independently of doping. This result is much larger than expected, and comparable to the total thermal conductivity of typical metals, contributing ~10 % to the total thermal conductivity of bulk graphene. Notably, in samples whose physical or domain sizes are of the order of few micrometers or smaller, the relative contribution coming from the electronic thermal conductivity is more important than in the bulk limit, since lattice thermal conductivity is much more sensitive to sample or grain size at these scales. Last, when electron-impurity scattering effects are included, we find that the electronic thermal conductivity is reduced by 30 to 70 %. We also find that the Wiedemann-Franz law is broadly satisfied at low and high temperatures, but with the largest deviations of 20-50 % around room temperature.

preprint2016arXiv

Thermal transport in crystals as a kinetic theory of relaxons

Thermal conductivity in dielectric crystals is the result of the relaxation of lattice vibrations described by the phonon Boltzmann transport equation. Remarkably, an exact microscopic definition of the heat carriers and their relaxation times is still missing: phonons, typically regarded as the relevant excitations for thermal transport, cannot be identified as the heat carriers when most scattering events conserve momentum and do not dissipate heat flux. This is the case for two-dimensional or layered materials at room temperature, or three-dimensional crystals at cryogenic temperatures. In this work we show that the eigenvectors of the scattering matrix in the Boltzmann equation define collective phonon excitations, termed here relaxons. These excitations have well defined relaxation times, directly related to heat flux dissipation, and provide an exact description of thermal transport as a kinetic theory of the relaxon gas. We show why Matthiessen's rule is violated, and construct a procedure for obtaining the mean free paths and relaxation times of the relaxons. These considerations are general, and would apply also to other semiclassical transport models, such as the electronic Boltzmann equation. For heat transport, they remain relevant even in conventional crystals like silicon, but are of the utmost importance in the case of two-dimensional materials, where they can revise by several orders of magnitude the relevant time- and length-scales for thermal transport in the hydrodynamic regime.

preprint2015arXiv

AiiDA: Automated Interactive Infrastructure and Database for Computational Science

Computational science has seen in the last decades a spectacular rise in the scope, breadth, and depth of its efforts. Notwithstanding this prevalence and impact, it is often still performed using the renaissance model of individual artisans gathered in a workshop, under the guidance of an established practitioner. Great benefits could follow instead from adopting concepts and tools coming from computer science to manage, preserve, and share these computational efforts. We illustrate here our paradigm sustaining such vision, based around the four pillars of Automation, Data, Environment, and Sharing. We then discuss its implementation in the open-source AiiDA platform (http://www.aiida.net), that has been tuned first to the demands of computational materials science. AiiDA's design is based on directed acyclic graphs to track the provenance of data and calculations, and ensure preservation and searchability. Remote computational resources are managed transparently, and automation is coupled with data storage to ensure reproducibility. Last, complex sequences of calculations can be encoded into scientific workflows. We believe that AiiDA's design and its sharing capabilities will encourage the creation of social ecosystems to disseminate codes, data, and scientific workflows.

preprint2015arXiv

Emergence of one-dimensional wires of free carriers in transition-metal-dichalcogenide nanostructures

We highlight the emergence of metallic states in two-dimensional transition-metal-dichalcogenide nanostructures -nanoribbons, islands, and inversion domain boundaries- as a widespread and universal phenomenon driven by the polar discontinuities occurring at their edges or boundaries. We show that such metallic states form one-dimensional wires of electrons or holes, with a free charge density that increases with the system size, up to complete screening of the polarization charge, and can also be controlled by the specific edge or boundary configurations, e.g. through chemisorption of hydrogen or sulfur atoms at the edges. For triangular islands, local polar discontinuities occur even in the absence of a total dipole moment for the island and lead to an accumulation of free carriers close to the edges, providing a consistent explanation of previous experimental observations. To further stress the universal character of these mechanisms, we show that polar discontinuities give rise to metallic states also at inversion domain boundaries. These findings underscore the potential of engineering transition-metal-dichalcogenide nanostructures for manifold applications in nano- and opto-electronics, spintronics, catalysis, and solar-energy harvesting.

preprint2015arXiv

First-principles investigation of organic photovoltaic materials C$_{60}$, C$_{70}$, [C$_{60}$]PCBM, and bis-[C$_{60}$]PCBM using a many-body $G_0W_0$-Lanczos approach

We present a first-principles investigation of the excited-state properties of electron acceptors in organic photovoltaics including C$_{60}$, C$_{70}$, [6,6]-phenyl-C$_{61}$-butyric-acid-methyl-ester ([C$_{60}$]PCBM), and bis-[C$_{60}$]PCBM using many-body perturbation theory within the Hedin's $G_0W_0$ approximation and an efficient Lanczos approach. Calculated vertical ionization potentials (VIP) and vertical electron affinities (VEA) of C$_{60}$ and C$_{70}$ agree very well with experimental values measured in gas phase. The density of states of all three molecules is also compared to photoemission and inverse photoemission spectra measured on thin-films, exhibiting a close agreement - a rigid energy-gap renormalization owing to intermolecular interactions in the thin-films. In addition, it is shown that the low-lying unoccupied states of [C$_{60}$]PCBM are all derived from the highest-occupied molecular orbitals and the lowest-unoccupied molecular orbitals of fullerene C$_{60}$. The functional side group in [C$_{60}$]PCBM introduces a slight electron transfer to the fullerene cage, resulting in small decreases of both VIP and VEA. This small change of VEA provides a solid justification for the increase of open-circuit voltage when replacing fullerene C$_{60}$ with [C$_{60}$]PCBM as the electron acceptor in bulk heterojunction polymer solar cells.

preprint2015arXiv

Thermoelastic properties of $α$-iron from first-principles

We calculate the thermomechanical properties of $α$-iron, and in particular its isothermal and adiabatic elastic constants, using first-principles total-energy and lattice-dynamics calculations, minimizing the quasi-harmonic vibrational free energy under finite strain deformations. Particular care is made in the fitting procedure for the static and temperature-dependent contributions to the free energy, in discussing error propagation for the two contributions separately, and in the verification and validation of pseudopotential and all-electron calculations. We find that the zero-temperature mechanical properties are sensitive to the details of the calculation strategy employed, and common semi-local exchange-correlation functionals provide only fair to good agreement with experimental elastic constants, while their temperature dependence is in excellent agreement with experiments in a wide range of temperature almost up to the Curie transition.

preprint2014arXiv

Electron-Phonon Interactions and the Intrinsic Electrical Resistivity of Graphene

We present a first-principles study of the temperature- and density-dependent intrinsic electrical resistivity of graphene. We use density-functional theory and density-functional perturbation theory together with very accurate Wannier interpolations to compute all electronic and vibrational properties and electron-phonon coupling matrix elements; the phonon-limited resistivity is then calculated within a Boltzmann-transport approach. An effective tight-binding model, validated against first-principles results, is also used to study the role of electron-electron interactions at the level of many-body perturbation theory. The results found are in excellent agreement with recent experimental data on graphene samples at high carrier densities and elucidate the role of the different phonon modes in limiting electron mobility. Moreover, we find that the resistivity arising from scattering with transverse acoustic phonons is 2.5 times higher than that from longitudinal acoustic phonons. Last, high-energy, optical, and zone-boundary phonons contribute as much as acoustic phonons to the intrinsic electrical resistivity even at room temperature and become dominant at higher temperatures.

preprint2014arXiv

Electrostatics of solvated systems in periodic boundary conditions

Continuum solvation methods can provide an accurate and inexpensive embedding of quantum simulations in liquid or complex dielectric environments. Notwithstanding a long history and manifold applications to isolated systems in open boundary conditions, their extension to materials simulations --- typically entailing periodic-boundary conditions --- is very recent, and special care is needed to address correctly the electrostatic terms. We discuss here how periodic-boundary corrections developed for systems in vacuum should be modified to take into account solvent effects, using as a general framework the self-consistent continuum solvation model developed within plane-wave density-functional theory [O. Andreussi et al. J. Chem. Phys. 136, 064102 (2012)]. A comprehensive discussion of real-space and reciprocal-space corrective approaches is presented, together with an assessment of their ability to remove electrostatic interactions between periodic replicas. Numerical results for zero-dimensional and two-dimensional charged systems highlight the effectiveness of the different suggestions, and underline the importance of a proper treatement of electrostatic interactions in first-principles studies of charged systems in solution.

preprint2014arXiv

Engineering polar discontinuities in honeycomb lattices

Unprecedented and fascinating phenomena have been recently observed at oxide interfaces between centrosymmetric cubic materials, such as LaAlO$_3$ and SrTiO$_3$, where a polar discontinuity across the boundary gives rise to polarization charges and electric fields that drive a metal-insulator transition, with the appearance of free carriers at the interface. Two-dimensional analogues of these systems are possible, and honeycomb lattices could offer a fertile playground, thanks to their versatility and the extensive on-going experimental efforts in graphene and related materials. Here we suggest different realistic pathways to engineer polar discontinuities across interfaces between honeycomb lattices, and support these suggestions with extensive first-principles calculations. Two broad approaches are discussed, that are based on (i) nanoribbons, where a polar discontinuity against the vacuum emerges, and (ii) selective functionalizations, where covalent ligands are used to engineer polar discontinuities by selective or total functionalization of the parent system. All the cases considered have the potential to deliver innovative applications in ultra-thin and flexible solar-energy devices and in micro- and nano-electronics.

preprint2014arXiv

Koopmans-compliant functionals and their performance against reference molecular data

Koopmans-compliant functionals emerge naturally from extending the constraint of piecewise linearity of the total energy as a function of the number of electrons to each fractional orbital occupation. When applied to approximate density-functional theory, these corrections give rise to orbital-density-dependent functionals and potentials. We show that the simplest implementations of Koopmans' compliance provide accurate estimates for the quasiparticle excitations and leave the total energy functional almost or exactly intact, i.e., they describe correctly electron removals or additions, but do not necessarily alter the electronic charge density distribution within the system. Additional functionals can then be constructed that modify the potential energy surface, including e.g. Perdew-Zunger corrections. These functionals become exactly one-electron self-interaction free and, as all Koopmans-compliant functionals, are approximately many-electron self-interaction free. We discuss in detail these different formulations, and provide extensive benchmarks for the 55 molecules in the reference G2-1 set, using Koopmans-compliant functionals constructed from local-density or generalized-gradient approximations. In all cases we find excellent performance in the electronic properties, comparable or improved with respect to that of many-body perturbation theories, such as G$_0$W$_0$ and self-consistent GW, at a fraction of the cost and in a variational framework that also delivers energy derivatives. Structural properties and atomization energies preserve or slightly improve the accuracy of the underlying density-functional approximations (Note: Supplemental Material is included in the source).

preprint2014arXiv

Phonon-limited resistivity of graphene by first-principle calculations: electron-phonon interactions, strain-induced gauge field and Boltzmann equation

Electron-phonon coupling in graphene is extensively modeled and simulated from first principles. We find that using an accurate model for the polarizations of the acoustic phonon modes is crucial to obtain correct numerical results. The interactions between electrons and acoustic phonon modes, the gauge field and deformation potential, are calculated at the DFT level in the framework of linear response. The zero-momentum limit of acoustic phonons is interpreted as a strain pattern, allowing the calculation of the acoustic gauge field parameter in the GW approximation. The role of electronic screening on the electron-phonon matrix elements is investigated. We then solve the Boltzmann equation semi-analytically in graphene, including both acoustic and optical phonon scattering. We show that, in the Bloch-Grüneisen and equipartition regimes, the electronic transport is mainly ruled by the unscreened acoustic gauge field, while the contribution due to the deformation potential is negligible and strongly screened. By comparing with experimental data, we show that the contribution of acoustic phonons to resistivity is doping- and substrate-independent. The DFT+GW approach underestimates this contribution to resistivity by about 30 %. Above 270K, the calculated resistivity underestimates the experimental one more severely, the underestimation being larger at lower doping. We show that, beside remote phonon scattering, a possible explanation for this disagreement is the electron-electron interaction that strongly renormalizes the coupling to intrinsic optical-phonon modes. Finally, after discussing the validity of the Matthiessen rule in graphene, we derive simplified analytical solutions of the Boltzmann equation to extract the coupling to acoustic phonons, related to the strain-induced gauge field, directly from experimental data.

preprint2014arXiv

Photoemission Spectroscopy and Orbital Imaging from Koopmans-Compliant Functionals

The determination of spectral properties from first principles can provide powerful connections between microscopic theoretical predictions and experimental data, but requires complex electronic-structure formulations that fall outside the domain of applicability of common approaches, such as density-functional theory. We show here that Koopmans-compliant functionals, constructed to enforce piecewise linearity in energy functionals with respect to fractional occupations - i.e. with respect to charged excitations - provide molecular photoemission spectra and momentum maps of Dyson orbitals that are in excellent agreement with experimental ultraviolet photoemission spectroscopy and orbital tomography data. These results highlight the role of Koopmans-compliant functionals as accurate and inexpensive quasiparticle approximations to the spectral potential.

preprint2014arXiv

Spin-resolved optical conductivity of two-dimensional group-VIB transition-metal dichalcogenides

We present an ab-initio study of the spin-resolved optical conductivity of two-dimensional (2D) group-VIB transition-metal dichalcogenides (TMDs). We carry out fully-relativistic density-functional-theory calculations combined with maximally localized Wannier functions to obtain band manifolds at extremely high resolutions and focus on the photo-response of 2D TMDs to circularly-polarized light in a wide frequency range. We present extensive numerical results for monolayer TMDs involving molybdenum and tungsten combined with sulphur and selenium. Our numerical approach allows us to locate with a high degree of accuracy the positions of the points in the Brillouin zone that are responsible for van Hove singularities in the optical response. Surprisingly, some of the saddle points do not occur exactly along high-symmetry directions in the Brillouin zone, although they happen to be in their close proximity.

preprint2013arXiv

BoltzWann: A code for the evaluation of thermoelectric and electronic transport properties with a maximally-localized Wannier functions basis

We present a new code to evaluate thermoelectric and electronic transport properties of extended systems with a maximally-localized Wannier function basis set. The semiclassical Boltzmann transport equations for the homogeneous infinite system are solved in the constant relaxation-time approximation and band energies and band derivatives are obtained via Wannier interpolations. Thanks to the exponential localization of the Wannier functions obtained, very high accuracy in the Brillouin zone integrals can be achieved with very moderate computational costs. Moreover, the analytical expression for the band derivatives in the Wannier basis resolves any issues that may occur when evaluating derivatives near band crossings. The code is tested on binary and ternary skutterudites CoSb_3 and CoGe_{3/2}S_{3/2}.

preprint2013arXiv

Defect ordering and defect-domain wall interactions in PbTiO$_3$: A first-principles study

The properties of ferroelectric materials, such as lead zirconate titanate (PZT), are heavily influenced by the interaction of defects with domain walls. These defects are either intrinsic, or are induced by the addition of dopants. We study here PbTiO$_3$ (the end member of a key family of solid solutions) in the presence of acceptor (Fe) and donor (Nb) dopants, and the interactions of the different defects and defect associates with the domain walls. For the case iron acceptors, the calculations point to the formation of defect associates involving an iron substitutional defect and a charged oxygen vacancy (Fe$^{'}_{Ti}$-V$^{^{\textbf{..}}}_O$). This associate exhibits a strong tendency to align in the direction of the bulk polarization; in fact, ordering of defects is also observed in pure PbTiO$_3$ in the form of lead-oxygen divacancies. Conversely, calculations on donor-doped PbTiO$_3$ do not indicate the formation of polar defect complexes involving donor substitutions. Last, it is observed that both isolated defects in donor-doped materials and defect associates in acceptor-doped materials are more stable at 180$^o$ domain walls. However, polar defect complexes lead to asymmetric potentials at domain walls due to the interaction of the defect polarization with the bulk polarization. The relative pinning characteristics of different defects are then compared, to develop an understanding of defect-domain wall interactions in both doped and pure PbTiO$_3$. These results may also help understanding hardening and softening mechanisms in PZT.

preprint2013arXiv

First-principles molecular dynamics simulations of proton diffusion in cubic BaZrO3 perovskite under strain conditions

First-principles molecular dynamics simulations have been employed to analyze the proton diffusion in cubic BaZrO3 perovskite at 1300K, and a non-linear effect of an applied isometric strain of 2% on the lattice parameter has been observed. The structural and electronic properties of BaZrO3 are analyzed, based on Density Functional Theory calculations, and after an analysis of the electronic structure, we provide a possible explanation for an enhanced ionic conductivity, that can be caused by the formation of a preferential path for proton diffusion under compressive strain conditions.

preprint2013arXiv

Thermodynamic stability of alkali metal/zinc double-cation borohydrides at low temperatures

We study the thermodynamic stability at low temperatures of a series of alkali metal/zinc double-cation borohydrides, including LiZn(BH$_4$)$_3$, LiZn$_2$(BH$_4$)$_5$, NaZn(BH$_4$)$_3$, NaZn$_2$(BH$_4$)$_5$, KZn(BH$_4$)$_3$, and KZn$_2$(BH$_4$)$_5$. While LiZn$_2$(BH$_4$)$_5$, NaZn(BH$_4$)$_3$, NaZn$_2$(BH$_4$)$_5$ and KZn(BH$_4$)$_3$ were recently synthesized, LiZn(BH$_4$)$_3$ and KZn$_2$(BH$_4$)$_5$ are hypothetical compounds. Using the minima-hopping method, we discover two new lowest-energy structures for NaZn(BH$_4$)$_3$ and KZn$_2$(BH$_4$)$_5$ which belong to the $C2/c$ and $P2$ space groups, respectively. These structures are predicted to be both thermodynamically stable and dynamically stable, implying that their existence may be possible. On the other hand, the lowest-energy $P1$ structure of LiZn(BH$_4$)$_3$ is predicted to be unstable, suggesting a possible reason elucidating why this compound has not been experimentally identified. In exploring the low-energy structures of these compounds, we find that their energetic ordering is sensitive to the inclusion of the van der Waals interactions. We also find that a proper treatment of these interactions, e.g., as given by a non-local density functional such as vdW-DF2, is necessary to address the stability of the low-energy structures of these compounds.

preprint2012arXiv

First Principles Quantum Transport with Electron-vibration Interactions: A Maximally Localized Wannier Function Approach

We present an ab initio inelastic quantum transport approach based on maximally localized Wannier functions. Electronic-structure properties are calculated with density-functional theory in a planewave basis, and electron-vibration coupling strengths and vibrational properties are determined with density-functional perturbation theory. Vibration-induced inelastic transport properties are calculated with non-equilibrium Green's function techniques, which are based on localized orbitals. For this purpose we construct maximally localized Wannier functions. Our formalism is applied to investigate inelastic transport in a benzene molecular junction connected to mono-atomic carbon chains. In this benchmark system the electron-vibration self-energy is calculated either in the self-consistent Born approximation or by lowest-order perturbation theory. It is observed that upward and downward conductance steps occur, which can be understood using multi-eigenchannel scattering theory and symmetry conditions. In a second example where the mono-atomic carbon chain electrode is replaced by a (3; 3) carbon nanotube, we focus on the non-equilibrium vibration populations driven by the conducting electrons using a semi-classical rate equation.

preprint2012arXiv

Maximally localized Wannier functions: Theory and applications

The electronic ground state of a periodic system is usually described in terms of extended Bloch orbitals, but an alternative representation in terms of localized "Wannier functions" was introduced by Gregory Wannier in 1937. The connection between the Bloch and Wannier representations is realized by families of transformations in a continuous space of unitary matrices, carrying a large degree of arbitrariness. Since 1997, methods have been developed that allow one to iteratively transform the extended Bloch orbitals of a first-principles calculation into a unique set of maximally localized Wannier functions, accomplishing the solid-state equivalent of constructing localized molecular orbitals, or "Boys orbitals" as previously known from the chemistry literature. These developments are reviewed here, and a survey of the applications of these methods is presented. This latter includes a description of their use in analyzing the nature of chemical bonding, or as a local probe of phenomena related to electric polarization and orbital magnetization. Wannier interpolation schemes are also reviewed, by which quantities computed on a coarse reciprocal-space mesh can be used to interpolate onto much finer meshes at low cost, and applications in which Wannier functions are used as efficient basis functions are discussed. Finally the construction and use of Wannier functions outside the context of electronic-structure theory is presented, for cases that include phonon excitations, photonic crystals, and cold-atom optical lattices.

preprint2012arXiv

Transport properties of room temperature ionic liquids from classical molecular dynamics

Room Temperature Ionic Liquids (RTILs) have attracted much of the attention of the scientific community in the past decade due the their novel and highly customizable properties. Nonetheless their high viscosities pose serious limitations to the use of RTILs in practical applications. To elucidate some of the physical aspects behind transport properties of RTILs, extensive classical molecular dynamics (MD) calculations are reported. Bulk viscosities and ionic conductivities of butyl-methyl-imidazole based RTILs are presented over a wide range of temperatures. The dependence of the properties of the liquids on simulation parameters, e.g. system size effects and choice of the interaction potential, is analyzed.

preprint2011arXiv

Automated quantum conductance calculations using maximally-localised Wannier functions

A robust, user-friendly, and automated method to determine quantum conductance in disordered quasi-one-dimensional systems is presented. The scheme relies upon an initial density- functional theory calculation in a specific geometry after which the ground-state eigenfunctions are transformed to a maximally-localised Wannier function (MLWF) basis. In this basis, our novel algorithms manipulate and partition the Hamiltonian for the calculation of coherent electronic transport properties within the Landauer-Buttiker formalism. Furthermore, we describe how short- ranged Hamiltonians in the MLWF basis can be combined to build model Hamiltonians of large (>10,000 atom) disordered systems without loss of accuracy. These automated algorithms have been implemented in the Wannier90 code[Mostofi et al, Comput. Phys. Commun. 178, 685 (2008)], which is interfaced to a number of electronic structure codes such as Quantum-ESPRESSO, AbInit, Wien2k, SIESTA and FLEUR. We apply our methods to an Al atomic chain with a Na defect, an axially heterostructured Si/Ge nanowire and to a spin-polarised defect on a zigzag graphene nanoribbon.

preprint2011arXiv

Berry Phase and Pseudospin Winding Number in Bilayer Graphene

Ever since the novel quantum Hall effect in bilayer graphene was discovered, and explained by a Berry phase of 2pi [K. S. Novoselov et al., "Unconventional quantum Hall effect and Berry's phase of 2pi in bilayer graphene", Nature Phys. 2, 177 (2006)], it has been widely accepted that the low-energy electronic wavefunction in this system is described by a non-trivial Berry phase of 2pi, different from the zero phase of a conventional two-dimensional electron gas. Here, we show that (i) the relevant Berry phase for bilayer graphene is not different from that for a conventional two-dimensional electron gas (as expected, given that Berry phase is only meaningful modulo 2pi) and that (ii) what is actually observed in the quantum Hall measurements is not the absolute value of the Berry phase but the pseudospin winding number.

preprint2011arXiv

Dominant Phonon Wavevectors and Strain-induced Splitting of the 2D Graphene Raman Mode

The dominant phonon wavevectors $q^{*}$ probed by the 2D Raman mode of graphene are highly anisotropic and rotate with the orientation of the polarizer:analyzer direction relative to the lattice. The corresponding electronic transitions connect the electronic equibandgap contours where the product of the ingoing and outgoing optical matrix elements is strongest, showing a finite component along $\bm{K}-\bmΓ$ that sensitively determines $q^{*}$. We revoke the notion of 'inner' and 'outer' processes. Our findings explain the splitting of the 2D mode of graphene under uniaxial tensile strain. The splitting originates from a strain-induced distortion of the phonon dispersion; changes in the electronic band structure and resonance conditions are negligeable for the 2D Raman spectrum.

preprint2011arXiv

Electronic, vibrational and transport properties of pnictogen substituted ternary skutterudites

First principles calculations are used to investigate electronic band structure and vibrational spectra of pnictogen substituted ternary skutterudites. We compare the results with the prototypical binary composition CoSb$_3$ to identify the effects of substitutions on the Sb site, and evaluate the potential of ternary skutterudites for thermoelectric applications. Electronic transport coefficients are computed within the Boltzmann transport formalism assuming a constant relaxation time, using a new methodology based on maximally localized Wannier function interpolation. Our results point to a large sensitivity of the electronic transport coefficients to carrier concentration and to scattering mechanisms associated with the enhanced polarity. The ionic character of the bonds is used to explain the detrimental effect on the thermoelectric properties.

preprint2011arXiv

Photoelectron properties of DNA and RNA bases from many-body perturbation theory

The photoelectron properties of DNA and RNA bases are studied using many-body perturbation theory within the GW approximation, together with a recently developed Lanczos-chain approach. Calculated vertical ionization potentials, electron affinities, and total density of states are in good agreement with experimental values and photoemission spectra. The convergence benchmark demonstrates the importance of using an optimal polarizability basis in the GW calculations. A detailed analysis of the role of exchange and correlation in both many-body and density-functional theory calculations shows that while self-energy corrections are strongly orbital-dependent, they nevertheless remain almost constant for states that share the same bonding character. Finally, we report on the inverse lifetimes of DNA and RNA bases, that are found to depend linearly on quasi-particle energies for all deep valence states. In general, our G0W0-Lanczos approach provides an efficient yet accurate and fully converged description of quasiparticle properties of five DNA and RNA bases.

preprint2011arXiv

Revised self-consistent continuum solvation in electronic-structure calculations

The solvation model proposed by Fattebert and Gygi [Journal of Computational Chemistry 23, 662 (2002)] and Scherlis et al. [Journal of Chemical Physics 124, 074103 (2006)] is reformulated, overcoming some of the numerical limitations encountered and extending its range of applicability. We first recast the problem in terms of induced polarization charges that act as a direct mapping of the self-consistent continuum dielectric; this allows to define a functional form for the dielectric that is well behaved both in the high-density region of the nuclear charges and in the low-density region where the electronic wavefunctions decay into the solvent. Second, we outline an iterative procedure to solve the Poisson equation for the quantum fragment embedded in the solvent that does not require multi-grid algorithms, is trivially parallel, and can be applied to any Bravais crystallographic system. Last, we capture some of the non-electrostatic or cavitation terms via a combined use of the quantum volume and quantum surface [Physical Review Letters 94, 145501 (2005)] of the solute. The resulting self-consistent continuum solvation (SCCS) model provides a very effective and compact fit of computational and experimental data, whereby the static dielectric constant of the solvent and one parameter allow to fit the electrostatic energy provided by the PCM model with a mean absolute error of 0.3 kcal/mol on a set of 240 neutral solutes. Two parameters allow to fit experimental solvation energies on the same set with a mean absolute error of 1.3 kcal/mol. A detailed analysis of these results, broken down along different classes of chemical compounds, shows that several classes of organic compounds display very high accuracy, with solvation energies in error of 0.3-0.4 kcal/mol, whereby larger discrepancies are mostly limited to self-dissociating species and strong hydrogen-bond forming compounds.

preprint2011arXiv

Variational Minimization of Orbital-dependent Density Functionals

Functionals that strive to correct for such self-interaction errors, such as those obtained by imposing the Perdew-Zunger self-interaction correction or the generalized Koopmans' condition, become orbital dependent or orbital-density dependent, and provide a very promising avenue to go beyond density-functional theory, especially when studying electronic, optical and dielectric properties, charge-transfer excitations, and molecular dissociations. Unlike conventional density functionals, these functionals are not invariant under unitary transformations of occupied electronic states, which leave the total charge density intact, and this added complexity has greatly inhibited both their development and their practical applicability. Here, we first recast the minimization problem for non-unitary invariant energy functionals into the language of ensemble density-functional theory, decoupling the variational search into an inner loop of unitary transformations that minimize the energy at fixed orbital subspace, and an outer-loop evolution of the orbitals in the space orthogonal to the occupied manifold. Then, we show that the potential energy surface in the inner loop is far from convex parabolic in the early stages of the minimization and hence minimization schemes based on these assumptions are unstable, and present an approach to overcome such difficulty. The overall formulation allows for a stable, robust, and efficient variational minimization of non-unitary-invariant functionals, essential to study complex materials and molecules, and to investigate the bulk thermodynamic limit, where orbitals converge typically to localized Wannier functions. In particular, using maximally localized Wannier functions as an initial guess can greatly reduce the computational costs needed to reach the energy minimum while not affecting or improving the convergence efficiency.

preprint2010arXiv

Koopmans' condition for density-functional theory

In approximate Kohn-Sham density-functional theory, self-interaction manifests itself as the dependence of the energy of an orbital on its fractional occupation. This unphysical behavior translates into qualitative and quantitative errors that pervade many fundamental aspects of density-functional predictions. Here, we first examine self-interaction in terms of the discrepancy between total and partial electron removal energies, and then highlight the importance of imposing the generalized Koopmans' condition -- that identifies orbital energies as opposite total electron removal energies -- to resolve this discrepancy. In the process, we derive a correction to approximate functionals that, in the frozen-orbital approximation, eliminates the unphysical occupation dependence of orbital energies up to the third order in the single-particle densities. This non-Koopmans correction brings physical meaning to single-particle energies; when applied to common local or semilocal density functionals it provides results that are in excellent agreement with experimental data -- with an accuracy comparable to that of GW many-body perturbation theory -- while providing an explicit total energy functional that preserves or improves on the description of established structural properties.

preprint2009arXiv

A first-principles characterization of the structure and electronic structure of $α$-S and Rh-S chalcogenides

We have used first-principles calculations to study the structural, electronic, and thermodynamic properties of the three known forms of Rh-S chalcogenides: Rh$_2$S$_3$, Rh$_3$S$_4$, and Rh$_{17}$S$_{15}$. Only the first of these materials of interest for catalysis had been studied previously within this approach. We find that Rh$_{17}$S$_{15}$ crystallizes in a {\em Pm}\={3}{\em m} centrosymmetric structure, as believed experimentally but never confirmed. We show band structures and densities of states for these compounds. Finally, we also investigated the ground state structure of solid sulfur ($α$-S), one of the few elements that represents a challenge for full first-principles calculations due to its demanding 128-atom unit cell and dispersion interactionis between S$_8$ units.

preprint2009arXiv

Dynamics and Thermodynamics of a Novel Phase of NaAlH4

We characterize a novel orthorhombic phase (gamma) of NaAlH4, discovered using first-principles molecular dynamics, and discuss its relevance to the dehydrogenation mechanism. This phase is close in energy to the known low-temperature structure and becomes the stabler phase above 320 K, thanks to a larger vibrational entropy associated with AlH4 rotational modes. The structural similarity of gamma-NaAlH4 to alpha-Na3AlH6 suggests it acts as a key intermediate during hydrogen release. Findings are consistent with recent experiments recording an unknown phase during dehydrogenation.

preprint2009arXiv

Linear and Non-linear Susceptibilities from Diffusion Quantum Monte Carlo: Application to Periodic Hydrogen Chains

We calculate the linear and non-linear susceptibilities of periodic longitudinal chains of hydrogen dimers with different bond-length alternations using a diffusion quantum Monte Carlo approach. These quantities are derived from the changes in electronic polarization as a function of applied finite electric field - an approach we recently introduced and made possible by the use of a Berry-phase, many-body electric-enthalpy functional. Calculated susceptibilities and hyper-susceptibilities are found to be in excellent agreement with the best estimates available from quantum chemistry - usually extrapolations to the infinite-chain limit of calculations for chains of finite length. It is found that while exchange effects dominate the proper description of the susceptibilities, second hyper-susceptibilities are greatly affected by electronic correlations. We also assess how different approximations to the nodal surface of the many-body wavefunction affect the accuracy of the calculated susceptibilities.

preprint2009arXiv

Surface energies, work functions, and surface relaxations of low index metallic surfaces from first-principles

We study the relaxations, surface energies, and work functions of low index metallic surfaces using pseudopotential plane-wave density-functional calculations within the generalized gradient approximation. We study here the (100), (110), and (111) surfaces of Al, Pd, Pt, and Au and the (0001) surface of Ti, chosen for their use as contact or lead materials in nanoscale devices. We consider clean, mostly non-reconstructed surfaces in the slab-supercell approximation. Particular attention is paid to the convergence of these quantities with respect to slab thickness; furthermore, different methodologies for the calculation of work functions and surfaces energies are compared. We find that the use of bulk references for calculations of surface energies and work functions can be detrimental to convergence unless numerical grids are closely matched, especially when surface relaxations are being considered. Our results and comparison show that calculated values often do not quantitatively match experimental values. This may be understandable for the surface relaxations and surface energies, where experimental values can have large error, but even for the work functions, neither local nor semi-local functionals emerge as an accurate choice for every case.

Nicola Marzari

What is connected

Connect this record

See the researcher in context

Building this map preview

61 published item(s)

Novel fast Li-ion conductors for solid-state electrolytes from first-principles

Extreme anharmonicity and thermal contraction of 1D wires

Bloch's theorem in orbital-density-dependent functionals: Band structures from Koopmans spectral functionals

HP -- A code for the calculation of Hubbard parameters using density-functional perturbation theory

Infrared-active phonons in one-dimensional materials and their spectroscopic signatures

Koopmans spectral functionals in periodic-boundary conditions

Many-body Green's function approach to lattice thermal transport

Microscopic picture of paraelectric perovskites from structural prototypes

Phonon-assisted luminescence in defect centers from many-body perturbation theory

Pulay forces in density-functional theory with extended Hubbard functionals: From nonorthogonalized to orthogonalized manifolds

Wigner formulation of thermal transport in solids

Self-consistent Hubbard parameters from density-functional perturbation theory in the ultrasoft and projector-augmented wave formulations

2-D materials for ultra-scaled field-effect transistors: hundred candidates under the ab initio microscope

Automated high-throughput Wannierisation

Electronic structure of pristine and Ni-substituted LaFeO$_3$ from near edge x-ray absorption fine structure experiments and first-principles simulations

Equipartition of Energy Defines the Size-Thickness Relationship in Liquid-Exfoliated Nanosheets

Hubbard-corrected density functional perturbation theory with ultrasoft pseudopotentials

Self-consistent DFT+$U$+$V$ study of oxygen vacancies in SrTiO$_3$

Two-dimensional materials from high-throughput computational exfoliation of experimentally known compounds

Abundance of $\mathbb{Z}_2$ topological order in exfoliable two-dimensional insulators

Emergent dual topology in the three-dimensional Kane-Mele Pt$_2$HgSe$_3$

Generalization of Fourier's law into viscous heat equations

Intrinsic edge excitons in two-dimensional MoS$_2$

Koopmans Meets Bethe-Salpeter: Excitonic Optical Spectra without GW

Operando XANES from first-principles and its application to iridium oxide

Photorealistic modelling of metals from first principles

Wannier90 as a community code: new features and applications

Band-like Electron Transport with Record-High Mobility in the TCNQ family

Performance of arsenene and antimonene double-gate MOSFETs from first principles

The electronic thermal conductivity of graphene

Thermal transport in crystals as a kinetic theory of relaxons

AiiDA: Automated Interactive Infrastructure and Database for Computational Science

Emergence of one-dimensional wires of free carriers in transition-metal-dichalcogenide nanostructures

First-principles investigation of organic photovoltaic materials C$_{60}$, C$_{70}$, [C$_{60}$]PCBM, and bis-[C$_{60}$]PCBM using a many-body $G_0W_0$-Lanczos approach

Thermoelastic properties of $α$-iron from first-principles

Electron-Phonon Interactions and the Intrinsic Electrical Resistivity of Graphene

Electrostatics of solvated systems in periodic boundary conditions

Engineering polar discontinuities in honeycomb lattices

Koopmans-compliant functionals and their performance against reference molecular data

Phonon-limited resistivity of graphene by first-principle calculations: electron-phonon interactions, strain-induced gauge field and Boltzmann equation

Photoemission Spectroscopy and Orbital Imaging from Koopmans-Compliant Functionals

Spin-resolved optical conductivity of two-dimensional group-VIB transition-metal dichalcogenides

BoltzWann: A code for the evaluation of thermoelectric and electronic transport properties with a maximally-localized Wannier functions basis

Defect ordering and defect-domain wall interactions in PbTiO$_3$: A first-principles study

First-principles molecular dynamics simulations of proton diffusion in cubic BaZrO3 perovskite under strain conditions

Thermodynamic stability of alkali metal/zinc double-cation borohydrides at low temperatures

First Principles Quantum Transport with Electron-vibration Interactions: A Maximally Localized Wannier Function Approach

Maximally localized Wannier functions: Theory and applications

Transport properties of room temperature ionic liquids from classical molecular dynamics

Automated quantum conductance calculations using maximally-localised Wannier functions

Berry Phase and Pseudospin Winding Number in Bilayer Graphene

Dominant Phonon Wavevectors and Strain-induced Splitting of the 2D Graphene Raman Mode

Electronic, vibrational and transport properties of pnictogen substituted ternary skutterudites

Photoelectron properties of DNA and RNA bases from many-body perturbation theory

Revised self-consistent continuum solvation in electronic-structure calculations

Variational Minimization of Orbital-dependent Density Functionals

Koopmans' condition for density-functional theory

A first-principles characterization of the structure and electronic structure of $α$-S and Rh-S chalcogenides

Dynamics and Thermodynamics of a Novel Phase of NaAlH4

Linear and Non-linear Susceptibilities from Diffusion Quantum Monte Carlo: Application to Periodic Hydrogen Chains

Surface energies, work functions, and surface relaxations of low index metallic surfaces from first-principles