BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Nicola Marzari

Nicola Marzari contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.mtrl-sciphysics.comp-phcond-mat.mes-hallcond-mat.str-elphysics.chem-phcond-mat.stat-mechphysics.app-phquant-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
28works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access