Showing 1–50 of 114 results for author: Franchini, C

Automated Workflow for Accurate High-Throughput GW Calculations

Authors: Lorenzo Varrassi, Florian Ellinger, Espen Flage-Larsen, Michael Wolloch, Georg Kresse, Nicola Marzari, Cesare Franchini

Abstract: The GW approximation represents the state-of-the-art ab-initio method for computing excited-state properties. Its execution requires control over a larger number of (often interdependent) parameters, and therefore its application in high-throughput studies is hindered by the intricate and time-consuming convergence process across a multi-dimensional parameter space. To address these challenges, he… ▽ More The GW approximation represents the state-of-the-art ab-initio method for computing excited-state properties. Its execution requires control over a larger number of (often interdependent) parameters, and therefore its application in high-throughput studies is hindered by the intricate and time-consuming convergence process across a multi-dimensional parameter space. To address these challenges, here we develop a fully-automated open-source workflow for G$_0$W$_0$ calculations within the AiiDA-VASP plugin architecture. The workflow is based on an efficient estimation of the errors on the quasi-particle (QP) energies due to basis-set truncation and the pseudo-potential norm violation, which allows a reduction of the dimensionality of the parameter space and avoids the need for multi-dimensional convergence searches. Protocol validation is conducted through a systematic comparison against established experimental and state-of-the-art GW data. To demonstrate the effectiveness of the approach, we construct a database of QP energies for a diverse dataset of over 320 bulk structures. The openly accessible workflow and resulting dataset can serve as a valuable resource and reference for conducting accurate data-driven research. △ Less

Submitted 21 October, 2024; originally announced October 2024.

Comments: 11 pages, 6 figures, 2 tables

Machine Learning Small Polaron Dynamics

Authors: Viktor C. Birschitzky, Luca Leoni, Michele Reticcioli, Cesare Franchini

Abstract: Polarons are crucial for charge transport in semiconductors, significantly impacting material properties and device performance. The dynamics of small polarons can be investigated using first-principles molecular dynamics. However, the limited timescale of these simulations presents a challenge for adequately sampling infrequent polaron hopping events. Here, we introduce a message-passing neural n… ▽ More Polarons are crucial for charge transport in semiconductors, significantly impacting material properties and device performance. The dynamics of small polarons can be investigated using first-principles molecular dynamics. However, the limited timescale of these simulations presents a challenge for adequately sampling infrequent polaron hopping events. Here, we introduce a message-passing neural network that learns the polaronic potential energy surface by encoding the polaronic state, allowing for simulations of polaron hopping dynamics at the nanosecond scale. By leveraging the statistical significance of the long timescale, our framework can accurately estimate polaron (anisotropic) mobilities and activation barriers in prototypical polaronic oxides across different scenarios (hole polarons in rocksalt MgO and electron polarons in pristine and F-doped rutile TiO$_2$) in excellent agreement with experimental observations. △ Less

Submitted 24 September, 2024; originally announced September 2024.

Electron-beam-induced adatom-vacancy-complexes in mono- and bilayer phosphorene

Authors: Carsten Speckmann, Andrea Angeletti, Lukáš Kývala, David Lamprecht, Felix Herterich, Clemens Mangler, Lado Filipovic, Christoph Dellago, Cesare Franchini, Jani Kotakoski

Abstract: Phosphorene, a puckered two-dimensional allotrope of phosphorus, has sparked considerable interest in recent years due to its potential especially for optoelectronic applications with its layer-number-dependant direct band gap and strongly bound excitons. However, detailed experimental characterization of its intrinsic defects as well as its defect creation characteristics under electron irradiati… ▽ More Phosphorene, a puckered two-dimensional allotrope of phosphorus, has sparked considerable interest in recent years due to its potential especially for optoelectronic applications with its layer-number-dependant direct band gap and strongly bound excitons. However, detailed experimental characterization of its intrinsic defects as well as its defect creation characteristics under electron irradiation are scarce. Here, we report on the creation and stability of a variety of defect configurations under 60 kV electron irradiation in mono- and bilayer phosphorene including the first experimental reports of stable adatom-vacancy-complexes. Displacement cross section measurements in bilayer phosphorene yield a value of 7.7 +- 1.4 barn with an estimated lifetime of adatom-vacancy-complexes of 19.9 +- 0.7 s, while some are stable for up to 68 s under continuous electron irradiation. Surprisingly, ab initio-based simulations indicate that the complexes should readily recombine, even in structures strained by up to 3 %. The presented results will help to improve the understanding of the wide variety of defects in phosphorene, their creation, and their stability, which may enable new pathways for defect engineered phosphorene devices. △ Less

Submitted 17 September, 2024; originally announced September 2024.

Comments: 11 pages, 7 figures

Structure and dynamics of the magnetite(001)/water interface from molecular dynamics simulations based on a neural network potential

Authors: Salvatore Romano, Pablo Montero de Hijes, Matthias Meier, Georg Kresse, Cesare Franchini, Christoph Dellago

Abstract: The magnetite/water interface is commonly found in nature and plays a crucial role in various technological applications. However, our understanding of its structural and dynamical properties at the molecular scale remains still limited. In this study, we develop an efficient Behler-Parrinello neural network potential (NNP) for the magnetite/water system, paying particular attention to the accurat… ▽ More The magnetite/water interface is commonly found in nature and plays a crucial role in various technological applications. However, our understanding of its structural and dynamical properties at the molecular scale remains still limited. In this study, we develop an efficient Behler-Parrinello neural network potential (NNP) for the magnetite/water system, paying particular attention to the accurate generation of reference data with density functional theory. Using this NNP, we performed extensive molecular dynamics simulations of the magnetite (001) surface across a wide range of water coverages, from the single molecule to bulk water. Our simulations revealed several new ground states of low coverage water on the Subsurface Cation Vacancy (SCV) model and yielded a density profile of water at the surface that exhibits marked layering. By calculating mean square displacements, we obtained quantitative information on the diffusion of water molecules on the SCV for different coverages, revealing significant anisotropy. Additionally, our simulations provided qualitative insights into the dissociation mechanisms of water molecules at the surface. △ Less

Submitted 6 September, 2024; v1 submitted 21 August, 2024; originally announced August 2024.

Evolution of the Coulomb interactions in correlated transition-metal perovskite oxides from the constrained random phase approximation

Authors: Liang Si, Peitao Liu, Cesare Franchini

Abstract: Determining the strength of electronic correlations of correlated electrons plays important roles in accurately describing the electronic structures and physical properties of transition-metal (TM) perovskite oxides. Here, we study the evolution of electronic interaction parameters as a function of $d$-electron occupancy in an extended class of TM perovskite oxides $AB$O$_3$ ($A$=Sr, Ca, and $B$=3… ▽ More Determining the strength of electronic correlations of correlated electrons plays important roles in accurately describing the electronic structures and physical properties of transition-metal (TM) perovskite oxides. Here, we study the evolution of electronic interaction parameters as a function of $d$-electron occupancy in an extended class of TM perovskite oxides $AB$O$_3$ ($A$=Sr, Ca, and $B$=3$d$-5$d$ TM elements) using the constrained random-phase-approximation method adopting two distinct models: $t_{2g}$-$t_{2g}$ and $d$-$dp$. For Sr$B$O$_3$ with $B$=Fe, Ru, and Ir, the $t_{2g}$-$t_{2g}$ model faces critical challenges, as the low-energy Hamiltonian spanning $t_{2g}$ manifolds is ill-defined. The $t_{2g}$-$t_{2g}$ model suggests that, for early $AB$O$_3$ series ($B$=$d^1$-$d^3$), the bare Coulomb interaction parameters $V$ remain nearly constant due to the competition between extended $t_{2g}$ Wannier orbitals and bandwidth reduction. As the $d$-electron filling increases, both partially screened Coulomb interaction parameters $U$ and fully screened Coulomb interaction parameters $W$ decrease, which are attributed to enhanced $e_g$-$t_{2g}$ and $e_g$-$p$ screenings. In contrast to the $t_{2g}$-$t_{2g}$ model, the $d$-$dp$ model effectively handles both early and late $AB$O$_3$ perovskites and reveals different trends. Specifically, $V$ varies inversely with the spreads of $d$-orbitals. $W$ reaches its minimum at the $d^3$ occupancy due to an interplay between increasing $d$-orbital localization and increasing screening effects. An unusual trend is observed for $U$, with local maxima at both $d^1$ and $d^4$ occupations. This can be understood from two aspects: (1) the increasing full screening effects from $d^1$ to $d^3$ and (2) the strongest $d$-$d$ and the weakest $d$-$p$ screening effects near $d^4$ for Sr$B$O$_3$. △ Less

Submitted 19 August, 2024; originally announced August 2024.

Comments: 21 pages, 10 figures, 3 tables

Hydrogen Diffusion in Magnesium Using Machine Learning Potentials

Authors: Andrea Angeletti, Luca Leoni, Dario Massa, Luca Pasquini, Stefanos Papanikolaou, Cesare Franchini

Abstract: Understanding and accurately predicting hydrogen diffusion in materials is challenging due to the complex interactions between hydrogen defects and the crystal lattice. These interactions span large length and time scales, making them difficult to address with standard ab initio techniques. This work addresses this challenge by employing accelerated machine learning (ML) molecular dynamics simulat… ▽ More Understanding and accurately predicting hydrogen diffusion in materials is challenging due to the complex interactions between hydrogen defects and the crystal lattice. These interactions span large length and time scales, making them difficult to address with standard ab initio techniques. This work addresses this challenge by employing accelerated machine learning (ML) molecular dynamics simulations through active learning. We conduct a comparative study of different ML-based interatomic potential schemes, including VASP, MACE, and CHGNet, utilizing various training strategies such as on-the-fly learning, pre-trained universal models, and fine-tuning. We obtain an optimal hydrogen diffusion coefficient value of \(2.1 \cdot 10^{-8}\) m\(^2\)/s at 673 K in MgH\(_{0.06}\), which aligns exceptionally well with experimental results, underlining the efficacy and accuracy of ML-assisted methodologies in the context of diffusive dynamics. Particularly, our procedure significantly reduces the computational effort associated with traditional transition state calculations or ad-hoc designed interatomic potentials. The results highlight the limitations of pre-trained universal solutions for defective materials and how they can be improved by fine-tuning. Specifically, fine-tuning the models on a database produced during on-the-fly training of VASP ML force field allows the retrieving of DFT-level accuracy at a fraction of the computational cost. △ Less

Submitted 30 July, 2024; originally announced July 2024.

Electron Mobilities in SrTiO$_3$ and KTaO$_3$: Role of Phonon Anharmonicity, Mass Renormalization and Disorder

Authors: Luigi Ranalli, Carla Verdi, Marios Zacharias, Jacky Evens, Feliciano Giustino, Cesare Franchini

Abstract: Accurately predicting carrier mobility in strongly anharmonic solids necessitates a precise characterization of lattice dyndamics as a function of temperature. We achieve consistency with experimental electron mobility data for bulk KTaO$_3$ and SrTiO$_3$ above 150 K by refining the Boltzmann transport equations. This refinement includes incorporating temperature-dependent anharmonic phonon eigenf… ▽ More Accurately predicting carrier mobility in strongly anharmonic solids necessitates a precise characterization of lattice dyndamics as a function of temperature. We achieve consistency with experimental electron mobility data for bulk KTaO$_3$ and SrTiO$_3$ above 150 K by refining the Boltzmann transport equations. This refinement includes incorporating temperature-dependent anharmonic phonon eigenfrequencies and eigenmodes into the electron-phonon interaction tensor, while maintaining the derivatives of the Kohn-Sham potential as computed in density functional perturbation theory. Using efficient machine-learned force fields and the stochastic self-consistent harmonic approximation, we accurately compute the dynamical matrices. At room temperature, the calculated mobility for SrTiO$_3$ exceeds experimental values by an order of magnitude, whereas the overestimation for KTaO$_3$ is way less pronounced. This discrepancy is explained through the more significant electron mass renormalization near the conduction band bottom due to anharmonic electron-phonon coupling and the presence of local disorder in SrTiO$_3$. △ Less

Submitted 26 July, 2024; originally announced July 2024.

Digging its own Site: Linear Coordination Stabilizes a Pt1/Fe2O3 Single-Atom Catalyst

Authors: Ali Rafsanjani-Abbasi, Florian Buchner, Faith J. Lewis, Lena Puntscher, Florian Kraushofer, Panukorn Sombut, Moritz Eder, Jiri Pavelec, Erik Rheinfrank, Giada Franceschi, Viktor C. Birschitzky, Michele Riva, Cesare Franchini, Michael Schmid, Ulrike Diebold, Matthias Meier, Georg K. H. Madsen, Gareth S. Parkinson

Abstract: Determining the local coordination of the active site is a pre-requisite for the reliable modeling of single-atom catalysts (SACs). Obtaining such information is difficult on powder-based systems, so much emphasis is placed on density functional theory-based computations based on idealized low-index surfaces of the support. In this work, we investigate how Pt atoms bind to the (1-102) facet of Fe2… ▽ More Determining the local coordination of the active site is a pre-requisite for the reliable modeling of single-atom catalysts (SACs). Obtaining such information is difficult on powder-based systems, so much emphasis is placed on density functional theory-based computations based on idealized low-index surfaces of the support. In this work, we investigate how Pt atoms bind to the (1-102) facet of Fe2O3, a common support material in SAC. Using a combination of scanning tunneling microscopy (STM), x-ray photoelectron spectroscopy (XPS), and an extensive computational evolutionary search, we find that Pt atoms significantly reconfigure the support lattice to facilitate a pseudo-linear coordination to surface oxygen atoms. Despite breaking three surface Fe-O bonds, this geometry is favored by 0.84 eV over the best configuration involving an unperturbed support. We suggest that the linear O-Pt-O configuration is common in reactive Pt-based SAC systems because it balances thermal stability with the ability to adsorb reactants from the gas phase, and that extensive structural searches are likely necessary to determine realistic active site geometry in single-atom catalysis. △ Less

Submitted 26 June, 2024; originally announced June 2024.

A Multi-Technique Study of C2H4 Adsorption on a Model Single-Atom Rh1 Catalyst

Authors: Chunlei Wang, Panukorn Sombut, Lena Puntscher, Manuel Ulreich, Jiri Pavelec, David Rath, Jan Balajka, Matthias Meier, Michael Schmid, Ulrike Diebold, Cesare Franchini, Gareth S. Parkinson

Abstract: Single-atom catalysts are potentially ideal model systems to investigate structure-function relationships in catalysis, if the active sites can be uniquely determined. In this work, we study the interaction of C2H4 with a model Rh/Fe3O4(001) catalyst that features 2-, 5-, and 6-fold coordinated Rh adatoms, as well as Rh clusters. Using multiple surface-sensitive techniques in combination with calc… ▽ More Single-atom catalysts are potentially ideal model systems to investigate structure-function relationships in catalysis, if the active sites can be uniquely determined. In this work, we study the interaction of C2H4 with a model Rh/Fe3O4(001) catalyst that features 2-, 5-, and 6-fold coordinated Rh adatoms, as well as Rh clusters. Using multiple surface-sensitive techniques in combination with calculations of density functional theory (DFT), we follow the thermal evolution of the system and disentangle the behavior of the different species. C2H4 adsorption is strongest at the 2-fold coordinated Rh1 with a DFT-determined adsorption energy of -2.26 eV. However, desorption occurs at lower temperatures than expected because the Rh migrates into substitutional sites within the support, where the molecule is more weakly bound. Adsorption at the 5-fold coordinated Rh sites is predicated to -1.49 eV, but the superposition of this signal with that from small Rh clusters and additional heterogeneity leads to a broad C2H4 desorption shoulder in TPD above room temperature. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Analytic method for quadratic polarons in nonparabolic bands

Authors: Serghei N. Klimin, Jacques Tempere, Matthew Houtput, Stefano Ragni, Thomas Hahn, Cesare Franchini, Andrey S. Mishchenko

Abstract: Including the effect of lattice anharmonicity on electron-phonon interactions has recently garnered attention due to its role as a necessary and significant component in explaining various phenomena, including superconductivity, optical response, and temperature dependence of mobility. This study focuses on analytically treating the effects of anharmonic electron-phonon coupling on the polaron sel… ▽ More Including the effect of lattice anharmonicity on electron-phonon interactions has recently garnered attention due to its role as a necessary and significant component in explaining various phenomena, including superconductivity, optical response, and temperature dependence of mobility. This study focuses on analytically treating the effects of anharmonic electron-phonon coupling on the polaron self-energy, combined with numerical Diagrammatic Monte Carlo data. Specifically, we incorporate a quadratic interaction into the method of squeezed phonon states, which has proven effective for analytically calculating the polaron parameters. Additionally, we extend this method to nonparabolic finite-width conduction bands while maintaining the periodic translation symmetry of the system. Our results are compared with those obtained from Diagrammatic Monte Carlo, partially reported in a recent study [S. Ragni et al., Phys. Rev. B 107, L121109(2023)], covering a wide range of coupling strengths for the nonlinear interaction. Remarkably, our analytic method predicts the same features as the Diagrammatic Monte Carlo simulation. △ Less

Submitted 10 August, 2024; v1 submitted 26 March, 2024; originally announced March 2024.

Comments: 32 pages, 8 figures

Journal ref: Phys. Rev. B 110, 075107 (2024)

Ultrastrong coupling between polar distortion and optical properties in ferroelectric MoBr$_2$O$_2$

Authors: Zhaojun Li, Lorenzo Varrassi, Yali Yang, Cesare Franchini, Laurent Bellaiche, Jiangang He

Abstract: Tuning the properties of materials using external stimuli is crucial for developing versatile smart materials. A strong coupling among order parameters within a single-phase material constitutes a potent foundation for achieving precise property control. However, cross-coupling is pretty weak in most single materials. Leveraging first principles calculations, we demonstrate the layered mixed anion… ▽ More Tuning the properties of materials using external stimuli is crucial for developing versatile smart materials. A strong coupling among order parameters within a single-phase material constitutes a potent foundation for achieving precise property control. However, cross-coupling is pretty weak in most single materials. Leveraging first principles calculations, we demonstrate the layered mixed anion compound MoBr$_2$O$_2$ exhibits electric-field switchable spontaneous polarization and ultrastrong coupling between polar distortion and electronic structures as well as optical properties. It offers feasible avenues of achieving tunable Rashba spin-splitting, electrochromism, thermochromism, photochromism, and nonlinear optics by applying an external electric field to a single domain sample, heating, as well as intense light illumination. Additionally, it exhibits an exceptionally large photostrictive effect. These findings not only showcase the feasibility of achieving multiple order parameter coupling within a single material, but also pave the way for comprehensive applications based on property control, such as energy harvesting, information processing, and ultrafast control. △ Less

Submitted 24 February, 2024; originally announced February 2024.

Interplay of superexchange and vibronic effects in the hidden order of Ba$_2$MgReO$_6$ unravelled from first principles

Authors: Dario Fiore Mosca, Cesare Franchini, Leonid V. Pourovskii

Abstract: The origin of the "hidden" quadrupolar and unconventional magnetic low-temperature orders observed in the spin-orbit double perovskite Ba$_2$MgReO$_6$ defies explanation through standard experimental and theoretical techniques. Here we address this problem by deriving and solving an ab initio low-temperature effective Hamiltonian including inter-site electronic exchange and vibronic (electron-latt… ▽ More The origin of the "hidden" quadrupolar and unconventional magnetic low-temperature orders observed in the spin-orbit double perovskite Ba$_2$MgReO$_6$ defies explanation through standard experimental and theoretical techniques. Here we address this problem by deriving and solving an ab initio low-temperature effective Hamiltonian including inter-site electronic exchange and vibronic (electron-lattice) couplings between $J_{eff}=3/2$ Jahn-Teller-active Rhenium states. Our findings disclose the nature of these elusive states, attributing it to intertwined exchange and electron-lattice couplings, thus diverging from the conventional dichotomy of purely electronic or lattice driving mechanisms. Our results indicate the resilience of the quadrupolar hidden order under pressure, yet its rapid suppression under uniaxial strain suggests that external or lattice-induced distortions play a pivotal role in determining the relative stability of competing phases in Ba$_2$MgReO$_6$ and similar $d^1$ double perovskites. △ Less

Submitted 27 February, 2024; v1 submitted 23 February, 2024; originally announced February 2024.

Comments: 6 pages plus 15 pages Supplementary. References in Supplementary fixed

Global sampling of Feynman's diagrams through Normalizing Flow

Authors: Luca Leoni, Cesare Franchini

Abstract: Normalizing Flows (NF) are powerful generative models with increasing applications in augmenting Monte Carlo algorithms due to their high flexibility and expressiveness. In this work we explore the integration of NF in Diagrammatic Monte Carlo (DMC), presenting an architecture designed to sample the intricate multidimensional space of Feynman's diagrams through dimensionality reduction. By decoupl… ▽ More Normalizing Flows (NF) are powerful generative models with increasing applications in augmenting Monte Carlo algorithms due to their high flexibility and expressiveness. In this work we explore the integration of NF in Diagrammatic Monte Carlo (DMC), presenting an architecture designed to sample the intricate multidimensional space of Feynman's diagrams through dimensionality reduction. By decoupling the sampling of diagram order and interaction times, the flow focuses on one interaction at a time. This enables constructing a general diagram by employing the same unsupervised model iteratively, dressing a zero-order diagram with interactions determined by the previously sampled order. The resulting NF-augmented DMC is tested on the widely used single-site Holstein polaron model in the entire electron-phonon coupling regime. The obtained data show that the model accurately reproduces the diagram distribution reducing sample correlation and statistical error for observables such as the polaron binding energy and the interacting Green's function. △ Less

Submitted 25 June, 2024; v1 submitted 1 February, 2024; originally announced February 2024.

Comments: 8 pages, 7 figures

Machine Learning Based Prediction of Polaron-Vacancy Patterns on the TiO$_2$(110) Surface

Authors: Viktor C. Birschitzky, Igor Sokolovic, Michael Prezzi, Krisztian Palotas, Martin Setvin, Ulrike Diebold, Michele Reticcioli, Cesare Franchini

Abstract: The multifaceted physics of oxides is shaped by their composition and the presence of defects, which are often accompanied by the formation of polarons. The simultaneous presence of polarons and defects, and their complex interactions, pose challenges for first-principles simulations and experimental techniques. In this study, we leverage machine learning and a first-principles database to analyze… ▽ More The multifaceted physics of oxides is shaped by their composition and the presence of defects, which are often accompanied by the formation of polarons. The simultaneous presence of polarons and defects, and their complex interactions, pose challenges for first-principles simulations and experimental techniques. In this study, we leverage machine learning and a first-principles database to analyze the distribution of surface oxygen vacancies (V$_{\rm O}$) and induced small polarons on rutile TiO$_2$(110), effectively disentangling the interactions between polarons and defects. By combining neural-network supervised learning and simulated annealing, we elucidate the inhomogeneous V$_{\rm O}$ distribution observed in scanning probe microscopy (SPM). Our innovative approach allows us to understand and predict defective surface patterns at previously inaccessible length scales, identifying the specific role of individual types of defects. Specifically, surface-polaron-stabilizing V$_{\rm O}$-configurations are identified, which could have consequences for surface reactivity. △ Less

Submitted 22 January, 2024; originally announced January 2024.

The origin of magnetism in a supposedly nonmagnetic osmium oxide

Authors: S. Agrestini, F. Borgatti, P. Florio, J. Frassineti, D. Fiore Mosca, Q. Faure, B. Detlefs, C. J. Sahle, S. Francoual, J. Choi, M. Garcia-Fernandez, K. -J. Zhou, V. F. Mitrovic, P. M. Woodward, G. Ghiringhelli, C. Franchini, F. Boscherini, S. Sanna, and M. Moretti Sala

Abstract: A supposedly nonmagnetic 5d$^1$ double perosvkite oxide is investigated by a combination of spectroscopic and theoretical methods, namely resonant inelastic X-ray scattering, X-ray absorption spectroscopy, magnetic circular dichroism, and multiplet ligand field calculations. We found that the large spin-orbit coupling admixes the 5d $t_{2g}$ and $e_g$ orbitals, covalency raises the 5d population w… ▽ More A supposedly nonmagnetic 5d$^1$ double perosvkite oxide is investigated by a combination of spectroscopic and theoretical methods, namely resonant inelastic X-ray scattering, X-ray absorption spectroscopy, magnetic circular dichroism, and multiplet ligand field calculations. We found that the large spin-orbit coupling admixes the 5d $t_{2g}$ and $e_g$ orbitals, covalency raises the 5d population well above the nominal value, and the local symmetry is lower than $O_h$. The obtained electronic interactions account for the finite magnetic moment of Os in this compound and, in general, of 5d$^1$ ions. Our results provide direct evidence of elusive Jahn-Teller distortions, hinting at a strong electron-lattice coupling. △ Less

Submitted 22 January, 2024; originally announced January 2024.

Journal ref: Phys. Rev. Lett. 133, 066501 (2024)

Pb$_9$Cu(PO$_4$)$_6$O is a charge-transfer semiconductor

Authors: Lorenzo Celiberti, Lorenzo Varrassi, Cesare Franchini

Abstract: By means of density functional theory and constrained random phase approximation we analyze the bandstructure of Pb$_9$Cu(PO$_4$)$_6$O (named LK-99). Our data show that the lead-phosphate apatite LK-99 in the proposed Cu-doped structure is a semiconductor with predominant charge-transfer nature, a result incompatible with a superconducting behaviour. In order to understand the interesting electron… ▽ More By means of density functional theory and constrained random phase approximation we analyze the bandstructure of Pb$_9$Cu(PO$_4$)$_6$O (named LK-99). Our data show that the lead-phosphate apatite LK-99 in the proposed Cu-doped structure is a semiconductor with predominant charge-transfer nature, a result incompatible with a superconducting behaviour. In order to understand the interesting electronic and magnetic properties of this compound, it will be necessary to study the actual response to doping, the possibility of alternative structural and stoichiometric (dis)orders and to clarify the magnetic interactions as well as their impact on the electronic structure. △ Less

Submitted 22 August, 2023; originally announced October 2023.

Oxygen-Terminated (1x1) Reconstruction of Reduced Magnetite Fe$_3$O$_4$(111)

Authors: Florian Kraushofer, Matthias Meier, Zdeněk Jakub, Johanna Hütner, Jan Balajka, Jan Hulva, Michael Schmid, Cesare Franchini, Ulrike Diebold, Gareth S. Parkinson

Abstract: The (111) facet of magnetite (Fe$_3$O$_4$) has been studied extensively by experimental and theoretical methods, but controversy remains regarding the structure of its low-energy surface terminations. Using density functional theory (DFT) computations, we demonstrate three reconstructions that are more favorable than the accepted Fe$_{\rm oct2}$ termination in reducing conditions. All three struct… ▽ More The (111) facet of magnetite (Fe$_3$O$_4$) has been studied extensively by experimental and theoretical methods, but controversy remains regarding the structure of its low-energy surface terminations. Using density functional theory (DFT) computations, we demonstrate three reconstructions that are more favorable than the accepted Fe$_{\rm oct2}$ termination in reducing conditions. All three structures change the coordination of iron in the kagome Fe$_{\rm oct1}$ layer to tetrahedral. With atomically-resolved microscopy techniques, we show that the termination that coexists with the Fe$_{\rm tet1}$ termination consists of tetrahedral iron capped by three-fold coordinated oxygen atoms. This structure explains the inert nature of the reduced patches. △ Less

Submitted 23 August, 2023; originally announced August 2023.

Journal ref: Journal of Physical Chemistry C Letters, 14 (2023) 3258 - 3265

A Multi-Technique Study of C2H4 Adsorption on Fe3O4(001)

Authors: Lena Puntscher, Panukorn Sombut, Chunlei Wang, Manuel Ulreich, Jiri Pavelec, Ali Rafsanjani-Abbasi, Matthias Meier, Adam Lagin, Martin Setvin, Ulrike Diebold, Cesare Franchini, Michael Schmid, Gareth S. Parkinson

Abstract: The adsorption/desorption of ethene (C2H4), also commonly known as ethylene, on Fe3O4(001) was studied under ultrahigh vacuum conditions using temperature programmed desorption (TPD), scanning tunneling microscopy, x-ray photoelectron spectroscopy, and density functional theory (DFT) based computations. To interpret the TPD data, we have employed a new analysis method based on equilibrium thermody… ▽ More The adsorption/desorption of ethene (C2H4), also commonly known as ethylene, on Fe3O4(001) was studied under ultrahigh vacuum conditions using temperature programmed desorption (TPD), scanning tunneling microscopy, x-ray photoelectron spectroscopy, and density functional theory (DFT) based computations. To interpret the TPD data, we have employed a new analysis method based on equilibrium thermodynamics. C2H4 adsorbs intact at all coverages and interacts most strongly with surface defects such as antiphase domain boundaries and Fe adatoms. On the regular surface, C2H4 binds atop surface Fe sites up to a coverage of 2 molecules per (rt2xrt2)R45° unit cell, with every second Fe occupied. A desorption energy of 0.36 eV is determined by analysis of the TPD spectra at this coverage, which is approximately 0.1-0.2 eV lower than the value calculated by DFT + U with van der Waals corrections. Additional molecules are accommodated in between the Fe rows. These are stabilized by attractive interactions with the molecules adsorbed at Fe sites. The total capacity of the surface for C2H4 adsorption is found to be close to 4 molecules per (rt2xrt2)R45° unit cell. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Evidence of Molecular Hydrogen in the N-doped LuH3 System: a Possible Path to Superconductivity?

Authors: Cesare Tresca, Pietro Maria Forcella, Andrea Angeletti, Luigi Ranalli, Cesare Franchini, Michele Reticcioli, Gianni Profeta

Abstract: The discovery of ambient superconductivity would mark an epochal breakthrough long-awaited for over a century, potentially ushering in unprecedented scientific and technological advancements. The recent findings on high-temperature superconducting phases in various hydrides under high pressure have ignited optimism, suggesting that the realization of near-ambient superconductivity might be on the… ▽ More The discovery of ambient superconductivity would mark an epochal breakthrough long-awaited for over a century, potentially ushering in unprecedented scientific and technological advancements. The recent findings on high-temperature superconducting phases in various hydrides under high pressure have ignited optimism, suggesting that the realization of near-ambient superconductivity might be on the horizon. However, the preparation of hydride samples tends to promote the emergence of various metastable phases, marked by a low level of experimental reproducibility. Identifying these phases through theoretical and computational methods entails formidable challenges, often resulting in controversial outcomes. In this paper, we consider N-doped LuH3 as a prototypical complex hydride: By means of machine-learning-accelerated force-field molecular dynamics, we have identified the formation of H2 molecules stabilized at ambient pressure by nitrogen impurities. Importantly, we demonstrate that this molecular phase plays a pivotal role in the emergence of a dynamically stable, low-temperature, experimental-ambient-pressure superconductivity. The potential to stabilize hydrogen in molecular form through chemical doping opens up a novel avenue for investigating disordered phases in hydrides and their transport properties under near-ambient conditions. △ Less

Submitted 22 May, 2024; v1 submitted 7 August, 2023; originally announced August 2023.

Journal ref: Nat. Comm. 2024

Spin-orbital Jahn-Teller bipolarons

Authors: Lorenzo Celiberti, Dario Fiore Mosca, Giuseppe Allodi, Leonid V. Pourovskii, Anna Tassetti, Paola Caterina Forino, Rong Cong, Erick Garcia, Phuong M. Tran, Roberto De Renzi, Patrick M. Woodward, Vesna F. Mitrović, Samuele Sanna, Cesare Franchini

Abstract: Polarons and spin-orbit (SO) coupling are distinct quantum effects that play a critical role in charge transport and spin-orbitronics. Polarons originate from strong electron-phonon interaction and are ubiquitous in polarizable materials featuring electron localization, in particular $\mathrm{3d}$ transition metal oxides (TMOs). On the other hand, the relativistic coupling between the spin and orb… ▽ More Polarons and spin-orbit (SO) coupling are distinct quantum effects that play a critical role in charge transport and spin-orbitronics. Polarons originate from strong electron-phonon interaction and are ubiquitous in polarizable materials featuring electron localization, in particular $\mathrm{3d}$ transition metal oxides (TMOs). On the other hand, the relativistic coupling between the spin and orbital angular momentum is notable in lattices with heavy atoms and develops in $\mathrm{5d}$ TMOs, where electrons are spatially delocalized. Here we combine ab initio calculations and magnetic measurements to show that these two seemingly mutually exclusive interactions are entangled in the electron-doped SO-coupled Mott insulator $\mathrm{Ba_2Na_{1-x}Ca_xOsO_6}$ ($0< x < 1$), unveiling the formation of spin-orbital bipolarons. Polaron charge trapping, favoured by the Jahn-Teller lattice activity, converts the Os $\mathrm{5d^1}$ spin-orbital $\mathrm{J_{eff}=3/2}$ levels, characteristic of the parent compound $\mathrm{Ba_2NaOsO_6}$ (BNOO), into a bipolaron $\mathrm{5d^2}$ $\mathrm{J_{eff}=2}$ manifold, leading to the coexistence of different J-effective states in a single-phase material. The gradual increase of bipolarons with increasing doping creates robust in-gap states that prevents the transition to a metal phase even at ultrahigh doping, thus preserving the Mott gap across the entire doping range from $\mathrm{d^1}$ BNOO to $\mathrm{d^2}$ $\mathrm{Ba_2CaOsO_6}$ (BCOO). △ Less

Submitted 27 June, 2023; originally announced June 2023.

Real-space investigation of polarons in hematite Fe2O3

Authors: Jesus Redondo, Michele Reticcioli, Vit Gabriel, Dominik Wrana, Florian Ellinger, Michele Riva, Giada Franceschi, Erik Rheinfrank, Igor Sokolovic, Zdenek Jakub, Florian Kraushofer, Aji Alexander, Laerte L. Patera, Jascha Repp, Michael Schmid, Ulrike Diebold, Gareth S. Parkinson, Cesare Franchini, Pavel Kocan, Martin Setvin

Abstract: In polarizable materials, electronic charge carriers interact with the surrounding ions, leading to quasiparticle behaviour. The resulting polarons play a central role in many materials properties including electrical transport, optical properties, surface reactivity and magnetoresistance, and polaron properties are typically investigated indirectly through such macroscopic characteristics. Here,… ▽ More In polarizable materials, electronic charge carriers interact with the surrounding ions, leading to quasiparticle behaviour. The resulting polarons play a central role in many materials properties including electrical transport, optical properties, surface reactivity and magnetoresistance, and polaron properties are typically investigated indirectly through such macroscopic characteristics. Here, noncontact atomic force microscopy (nc-AFM) is used to directly image polarons in Fe2O3 at the single quasiparticle limit. A combination of Kelvin probe force microscopy (KPFM) and kinetic Monte Carlo (KMC) simulations shows that Ti doping dramatically enhances the mobility of electron polarons, and density functional theory (DFT) calculations indicate that a metallic transition state is responsible for the enhancement. In contrast, hole polarons are significantly less mobile and their hopping is hampered further by the introduction of trapping centres. △ Less

Submitted 31 March, 2023; originally announced March 2023.

The Mott transition in the 5d$^1$ compound Ba$_2$NaOsO$_6:$ a DFT+DMFT study with PAW spinor projectors

Authors: Dario Fiore Mosca, Hermann Schnait, Lorenzo Celiberti, Markus Aichhorn, Cesare Franchini

Abstract: Spin-orbit coupling has been reported to be responsible for the insulating nature of the 5d$^1$ osmate double perovskite Ba$_2$NaOsO$_6$ (BNOO). However, whether spin-orbit coupling indeed drives the metal-to-insulator transition (MIT) in this compound is an open question. In this work we investigate the impact of relativistic effects on the electronic properties of BNOO via density functional the… ▽ More Spin-orbit coupling has been reported to be responsible for the insulating nature of the 5d$^1$ osmate double perovskite Ba$_2$NaOsO$_6$ (BNOO). However, whether spin-orbit coupling indeed drives the metal-to-insulator transition (MIT) in this compound is an open question. In this work we investigate the impact of relativistic effects on the electronic properties of BNOO via density functional theory plus dynamical mean-field theory calculations in the paramagnetic regime, where the insulating phase is experimentally observed. The correlated subspace is modeled with spinor projectors of the projector augumented wave method (PAW) employed in the Vienna Ab Initio Simulation Package (VASP), suitably interfaced with the TRIQS package. The inclusion of PAW spinor projectors in TRIQS enables the treatment of spin-orbit coupling effects fully ab-initio within the dynamical mean-field theory framework. In the present work, we show that spin-orbit coupling, although assisting the MIT in BNOO, is not the main driving force for its gapped spectra, placing this material in the Mott insulator regime. Relativistic effects primarily impact the correlated states' character, excitations, and magnetic ground-state properties. △ Less

Submitted 15 January, 2024; v1 submitted 29 March, 2023; originally announced March 2023.

Quasiparticle and excitonic properties of monolayer SrTiO$_3$

Authors: Lorenzo Varrassi, Peitao Liu, Cesare Franchini

Abstract: Strontium titanate SrTiO$_3$ is one of the most studied and paradigmatic transition metal oxides. Recently, a breakthrough has been achieved with the fabrication of freestanding SrTiO$_3$ ultrathin films down to the monolayer limit. However, the many-body effects on the quasiparticle and optical properties of monolayer SrTiO$_3$ remain unexplored. Using state-of-the-art many-body perturbation theo… ▽ More Strontium titanate SrTiO$_3$ is one of the most studied and paradigmatic transition metal oxides. Recently, a breakthrough has been achieved with the fabrication of freestanding SrTiO$_3$ ultrathin films down to the monolayer limit. However, the many-body effects on the quasiparticle and optical properties of monolayer SrTiO$_3$ remain unexplored. Using state-of-the-art many-body perturbation theory in the GW approximation combined with the Bethe-Salpeter equation, we study the quasiparticle band structure, optical and excitonic properties of monolayer SrTiO$_3$. We show that quasiparticle corrections significantly alter the band structure topology; however, the widely used diagonal $G_0W_0$ approach yields unphysical band dispersions. The correct band dispersions are restored only by taking into account the off-diagonal elements of the self-energy. The optical properties are studied both in the optical limit and for finite momenta by computing the electron energy loss spectra. We find that the imaginary part of dielectric function at the long wavelength limit is dominated by three strongly bound excitonic peaks and the direct optical gap is associated to a bright exciton state with a large binding energy of 0.93 eV. We discuss the character of the excitonic peaks via the contributing interband transitions, and reveal that the lowest bound excitonic state becomes optical inactive for finite momenta along $Γ$-M, while the other two excitonic peaks disperse to higher energies and eventually merge for momenta close to M. △ Less

Submitted 26 March, 2023; originally announced March 2023.

Comments: 10 pages, 4 figures

Polaron with Quadratic Electron-phonon Interaction

Authors: Stefano Ragni, Thomas Hahn, Zhongjin Zhang, Nikolay Prokof'ev, Anatoly Kuklov, Serghei Klimin, Matthew Houtput, Boris Svistunov, Jacques Tempere, Naoto Nagaosa, Cesare Franchini, Andrey S. Mishchenko

Abstract: We present the first numerically exact study of a polaron with quadratic coupling to the oscillator displacement, using two alternative methodological developments. Our results cover both anti-adiabatic and adiabatic regimes and the entire range of electron-phonon coupling $g_2$, from the system's stability threshold at attractive $g_2=-1$ to arbitrary strong repulsion at $g_2 \gg 1$. Key properti… ▽ More We present the first numerically exact study of a polaron with quadratic coupling to the oscillator displacement, using two alternative methodological developments. Our results cover both anti-adiabatic and adiabatic regimes and the entire range of electron-phonon coupling $g_2$, from the system's stability threshold at attractive $g_2=-1$ to arbitrary strong repulsion at $g_2 \gg 1$. Key properties of quadratic polarons prove dramatically different from their linear counterparts. They (i) are insensitive even to large quadratic coupling except in the anti-adiabatic limit near the threshold of instability at attraction; (ii) depend only on the adiabatic ratio but are insensitive to the electron dispersion and dimension of space; (iii) feature weak lattice deformations even at the instability point. Our results are of direct relevance to properties of electrons at low densities in polar materials, including recent proposals for their superconducting states. △ Less

Submitted 2 December, 2022; originally announced December 2022.

Comments: 6 pages, 6 figures + supplemental material

Journal ref: Phys. Rev. B 107, L121109 (2023)

Quantum paraelectricity and structural phase transitions in strontium titanate beyond density-functional theory

Authors: Carla Verdi, Luigi Ranalli, Cesare Franchini, Georg Kresse

Abstract: We demonstrate an approach for calculating temperature-dependent quantum and anharmonic effects with beyond density-functional theory accuracy. By combining machine-learned potentials and the stochastic self-consistent harmonic approximation, we investigate the cubic to tetragonal transition in strontium titanate and show that the paraelectric phase is stabilized by anharmonic quantum fluctuations… ▽ More We demonstrate an approach for calculating temperature-dependent quantum and anharmonic effects with beyond density-functional theory accuracy. By combining machine-learned potentials and the stochastic self-consistent harmonic approximation, we investigate the cubic to tetragonal transition in strontium titanate and show that the paraelectric phase is stabilized by anharmonic quantum fluctuations. We find that a quantitative understanding of the quantum paraelectric behavior requires a higher-level treatment of electronic correlation effects via the random phase approximation. This approach enables detailed studies of emergent properties in strongly anharmonic materials beyond density-functional theory. △ Less

Submitted 1 March, 2023; v1 submitted 17 November, 2022; originally announced November 2022.

Effects of charge doping on Mott insulator with strong spin-orbit coupling, Ba$_2$NaOsO$_6$

Authors: E. Garcia, R. Cong, P. C. Forino, A. Tassetti, G. Allodi, A. P. Reyes, P. M. Tran, P. M. Woodward, C. Franchini, S. Sanna, V. F. Mitrović

Abstract: The effects of doping on the electronic evolution of the Mott insulating state have been extensively studied in efforts to understand mechanisms of emergent quantum phases of materials. The study of these effects becomes ever more intriguing in the presence of entanglement between spin and orbital degrees of freedom. Here, we present a comprehensive investigation of charge doping in the double per… ▽ More The effects of doping on the electronic evolution of the Mott insulating state have been extensively studied in efforts to understand mechanisms of emergent quantum phases of materials. The study of these effects becomes ever more intriguing in the presence of entanglement between spin and orbital degrees of freedom. Here, we present a comprehensive investigation of charge doping in the double perovskite Ba$_2$NaOsO$_6$, a a complex Mott insulator where such entanglement plays an important role. We establish that the insulating magnetic ground state evolves from canted antiferromagnet (cAF)to Néel order for dopant levels exceeding ~ 10 %. Furthermore, we determine that a broken local point symmetry (BLPS) phase, precursor to the magnetically ordered state, occupies an extended portion of the (H-T) phase diagram with increased doping. This finding reveals that the breaking of the local cubic symmetry is driven by a multipolar order, most-likely of the antiferro-quadrupolar type. △ Less

Submitted 23 November, 2022; v1 submitted 10 October, 2022; originally announced October 2022.

Temperature-dependent anharmonic phonons in quantum paraelectric KTaO$_3$ by first principles and machine-learned force fields

Authors: Luigi Ranalli, Carla Verdi, Lorenzo Monacelli, Matteo Calandra, Georg Kresse, Cesare Franchini

Abstract: Understanding collective phenomena in quantum materials from first principles is a promising route toward engineering materials properties on demand and designing new functionalities. This work examines the quantum paraelectric state, an elusive state of matter characterized by the smooth saturation of the ferroelectric instability at low temperature due to quantum fluctuations associated with anh… ▽ More Understanding collective phenomena in quantum materials from first principles is a promising route toward engineering materials properties on demand and designing new functionalities. This work examines the quantum paraelectric state, an elusive state of matter characterized by the smooth saturation of the ferroelectric instability at low temperature due to quantum fluctuations associated with anharmonic phonon effects. The temperature-dependent evolution of the soft ferroelectric phonon mode in the quantum paraelectric KTaO$_3$ in the range 0-300 K is modelled by combining density functional theory (DFT) calculations with the stochastic self-consistent harmonic approximation assisted by an on-the-fly machine-learned force field. The calculated data show that including anharmonic terms is essential to stabilize the spurious imaginary ferroelectric phonon predicted by DFT, in agreement with experiments. Augmenting the DFT workflow with machine-learned force fields allows for efficient stochastic sampling of the configurational space using large supercells in a broad and dense temperature range, inaccessible by conventional ab initio protocols. This work proposes a robust computational workflow capable of accounting for collective behaviors involving different degrees of freedom and occurring at large time/length scales, paving the way for precise modeling and control of quantum effects in materials. △ Less

Submitted 24 September, 2022; originally announced September 2022.

Comments: 10 pages, 6 figures. To be published in Advanced Quantum Technologies. Editor: Luigi Ranalli

CO oxidation by Pt2/Fe3O4: metastable dimer and support configurations facilitate lattice oxygen extraction

Authors: Matthias Meier, Jan Hulva, Zdenek Jakub, Florian Kraushofer, Mislav Bobić, Roland Bliem, Martin Setvin, Michael Schmid, Ulrike Diebold, Cesare Franchini, Gareth S. Parkinson

Abstract: Heterogeneous catalysts based on sub-nanometer metal clusters often exhibit strongly size-dependent properties, and the addition or removal of a single atom can make all the difference. Identifying the most active species and deciphering the reaction mechanism is extremely difficult, however, because it is often not clear how the catalyst evolves in operando. Here, we utilize a combination of atom… ▽ More Heterogeneous catalysts based on sub-nanometer metal clusters often exhibit strongly size-dependent properties, and the addition or removal of a single atom can make all the difference. Identifying the most active species and deciphering the reaction mechanism is extremely difficult, however, because it is often not clear how the catalyst evolves in operando. Here, we utilize a combination of atomically resolved scanning probe microscopies, spectroscopic techniques, and density functional theory (DFT)-based calculations to study CO oxidation by a model Pt/Fe3O4(001) single-atom catalyst. We demonstrate that (PtCO)2 dimers, formed dynamically through the agglomeration of mobile Pt-carbonyl species, catalyse a reaction involving the oxide support to form CO2. Pt2 dimers produce one CO2 molecule before falling apart into two adatoms, releasing the second CO. Interestingly, Olattice extraction only becomes facile when both the Pt-dimer and the Fe3O4 support can access metastable configurations, suggesting that substantial, concerted rearrangements of both cluster and support must be considered for reactions occurring at elevated temperature. △ Less

Submitted 9 September, 2022; originally announced September 2022.

Journal ref: Science Advances 8 (13), eabn4580, 2022

Rapid oxygen exchange between hematite and water vapor

Authors: Zdenek Jakub, Matthias Meier, Florian Kraushofer, Jan Balajka, Jiri Pavelec, Michael Schmid, Cesare Franchini, Ulrike Diebold, Gareth S. Parkinson

Abstract: Oxygen exchange at oxide/liquid and oxide/gas interfaces is important in technology and environmental studies, as it is closely linked to both catalytic activity and material degradation. The atomic-scale details are mostly unknown, however, and are often ascribed to poorly defined defects in the crystal lattice. Here we show that even thermodynamically stable, well-ordered surfaces can be surpris… ▽ More Oxygen exchange at oxide/liquid and oxide/gas interfaces is important in technology and environmental studies, as it is closely linked to both catalytic activity and material degradation. The atomic-scale details are mostly unknown, however, and are often ascribed to poorly defined defects in the crystal lattice. Here we show that even thermodynamically stable, well-ordered surfaces can be surprisingly reactive. Specifically, we show that all the 3-fold coordinated lattice oxygen atoms on a defect-free single-crystalline r-cut (1-102) surface of hematite (α-Fe2O3) are exchanged with oxygen from surrounding water vapor within minutes at temperatures below 70 °C, while the atomic-scale surface structure is unperturbed by the process. A similar behavior is observed after liquid water exposure, but the experimental data clearly show most of the exchange happens during desorption of the final monolayer, not during immersion. Density functional theory computations show that the exchange can happen during on-surface diffusion, where the cost of the lattice oxygen extraction is compensated by the stability of an HO-HOH-OH complex. Such insights into lattice oxygen stability are highly relevant for many research fields ranging from catalysis and hydrogen production to geochemistry and paleoclimatology. △ Less

Submitted 9 September, 2022; originally announced September 2022.

Journal ref: Nature communications 12 (1), 1-8, 2021

Automated Real-Space Lattice Extraction for Atomic Force Microscopy Images

Authors: Marco Corrias, Lorenzo Papa, Igor Sokolović, Viktor Birschitzky, Alexander Gorfer, Martin Setvin, Michael Schmid, Ulrike Diebold, Michele Reticcioli, Cesare Franchini

Abstract: Analyzing atomically resolved images is a time-consuming process requiring solid experience and substantial human intervention. In addition, the acquired images contain a large amount of information such as crystal structure, presence and distribution of defects, and formation of domains, which need to be resolved to understand a material's surface structure. Therefore, machine learning techniques… ▽ More Analyzing atomically resolved images is a time-consuming process requiring solid experience and substantial human intervention. In addition, the acquired images contain a large amount of information such as crystal structure, presence and distribution of defects, and formation of domains, which need to be resolved to understand a material's surface structure. Therefore, machine learning techniques have been applied in scanning probe and electron microscopies during the last years, aiming for automatized and efficient image analysis. This work introduces a free and open source tool (AiSurf: Automated Identification of Surface Images) developed to inspect atomically resolved images via Scale-Invariant Feature Transform (SIFT) and Clustering Algorithms (CA). AiSurf extracts primitive lattice vectors, unit cells, and structural distortions from the original image, with no pre-assumption on the lattice and minimal user intervention. The method is applied to various atomically resolved non-contact atomic force microscopy (AFM) images of selected surfaces with different levels of complexity: anatase TiO2(101), oxygen deficient rutile TiO2(110) with and without CO adsorbates, SrTiO3(001) with Sr vacancies and graphene with C vacancies. The code delivers excellent results and has proved to be robust against atom misclassification and noise, thereby facilitating the interpretation scanning probe microscopy images. △ Less

Submitted 28 August, 2022; originally announced August 2022.

Journal ref: 2023 Mach. Learn.: Sci. Technol. 4 015015

Small Polaron Formation on the Nb-doped SrTiO$_\textbf{3}$(001) Surface

Authors: Florian Ellinger, Muhammad Shafiq, Iftikhar Ahmad, Michele Reticcioli, Cesare Franchini

Abstract: The cubic perovsike strontium titanate SrTiO$_3$ (STO) is one of the most studied, polarizable transition metal oxides. When excess charge is introduced to this material, e.g., through doping or atomic defects, STO tends to host polarons: Quasi-particles formed by excess charge carriers coupling with the crystal phonon field. Their presence alters the materials properties, and is a key for many ap… ▽ More The cubic perovsike strontium titanate SrTiO$_3$ (STO) is one of the most studied, polarizable transition metal oxides. When excess charge is introduced to this material, e.g., through doping or atomic defects, STO tends to host polarons: Quasi-particles formed by excess charge carriers coupling with the crystal phonon field. Their presence alters the materials properties, and is a key for many applications. Considering that polarons form preferentially on or near surfaces, we study small polaron formation at the TiO$_2$ termination of the STO(001) surface via density functional theory calculations. We model several supercell slabs of Nb-doped and undoped STO(001) surfaces with increasing size, also considering the recently observed as-cleaved TiO$_2$ terminated surface hosting Sr-adatoms. Our findings suggest that small polarons become less stable at low concentrations of Nb-doping, in analogy with polarons localized in the bulk. Further, we inspect the stability of different polaron configurations with respect to Nb- and Sr-impurities, and discuss their spectroscopic properties. △ Less

Submitted 12 September, 2022; v1 submitted 22 August, 2022; originally announced August 2022.

Comments: 8 pages, 7 figures, 1 table

Competing electronic states emerging on polar surfaces

Authors: Michele Reticcioli, Zhichang Wang, Michael Schmid, Dominik Wrana, Lynn A. Boatner, Ulrike Diebold, Martin Setvin, Cesare Franchini

Abstract: Excess charge on polar surfaces of ionic compounds is commonly described by the two-dimensional electron gas (2DEG) model, a homogeneous distribution of charge, spatially-confined in a few atomic layers. Here, by combining scanning probe microscopy with density functional theory calculations, we show that excess charge on the polar TaO$_2$ termination of KTaO$_3$(001) forms more complex electronic… ▽ More Excess charge on polar surfaces of ionic compounds is commonly described by the two-dimensional electron gas (2DEG) model, a homogeneous distribution of charge, spatially-confined in a few atomic layers. Here, by combining scanning probe microscopy with density functional theory calculations, we show that excess charge on the polar TaO$_2$ termination of KTaO$_3$(001) forms more complex electronic states with different degrees of spatial and electronic localization: charge density waves (CDW) coexist with strongly-localized electron polarons and bipolarons. These surface electronic reconstructions, originating from the combined action of electron-lattice interaction and electronic correlation, are energetically more favorable than the 2DEG solution. They exhibit distinct spectroscopy signals and impact on the surface properties, as manifested by a local suppression of ferroelectric distortions. Controlling the degree of charge ordering and the transition from ferroelectric to paraelectric states could be of great benefit for the generation and transport of carriers in electronic applications. △ Less

Submitted 1 July, 2022; originally announced July 2022.

Suppressed fluctuations as the origin of the static magnetic order in strained Sr$_2$RuO$_4$

Authors: Bongjae Kim, Sergii Khmelevskyi, Cesare Franchini, I. I. Mazin

Abstract: Combining first principle density functional calculations and Moriya's self-consistent renormalization theory, we explain the recently reported counterintuitive appearance of an ordered magnetic state in uniaxially strained Sr$_2$RuO$_4$ beyond the Lifshits transition. We show that strain weakens the quantum spin fluctuations, which destroy the static order, more strongly than the tendency to magn… ▽ More Combining first principle density functional calculations and Moriya's self-consistent renormalization theory, we explain the recently reported counterintuitive appearance of an ordered magnetic state in uniaxially strained Sr$_2$RuO$_4$ beyond the Lifshits transition. We show that strain weakens the quantum spin fluctuations, which destroy the static order, more strongly than the tendency to magnetism. A different rate of decrease of the spin fluctuations vs. magnetic stabilization energy promotes the onset of a static magnetic order beyond a critical strain. △ Less

Submitted 29 November, 2022; v1 submitted 28 June, 2022; originally announced June 2022.

Role of Polarons in Single-Atom Catalysts: Case Study of Me1 [Au1, Pt1, and Rh1] on TiO2(110)

Authors: Panukorn Sombut, Lena Puntscher, Marlene Atzmueller, Zdenek Jakub, Michele Reticcioli, Matthias Meier, Gareth S. Parkinson, Cesare Franchini

Abstract: The local environment of metal-oxide supported single-atom catalysts plays a decisive role in the surface reactivity and related catalytic properties. The study of such systems is complicated by the presence of point defects on the surface, which are often associated with the localization of excess charge in the form of polarons. This can affect the stability, the electronic configuration, and the… ▽ More The local environment of metal-oxide supported single-atom catalysts plays a decisive role in the surface reactivity and related catalytic properties. The study of such systems is complicated by the presence of point defects on the surface, which are often associated with the localization of excess charge in the form of polarons. This can affect the stability, the electronic configuration, and the local geometry of the adsorbed adatoms. In this work, through the use of density functional theory and surface-sensitive experiments, we study the adsorption of Rh1, Pt1, and Au1 metals on the reduced TiO2(110) surface; a prototypical polaronic material. A systematic analysis of the adsorption configurations and oxidation states of the adsorbed metals reveals different types of couplings between adsorbates and polarons. As confirmed by scanning tunneling microscopy measurements, the favored Pt1 and Au1 adsorption at oxygen vacancy sites is associated with a strong electronic charge transfer from polaronic states to adatom orbitals, which results in a reduction of the adsorbed metal. In contrast, the Rh1 adatoms interact weakly with the excess charge, which leaves the polarons largely unaffected. Our results show that an accurate understanding of the properties of single-atom catalysts on oxide surfaces requires a careful account of the interplay between adatoms, vacancy sites, and polarons. △ Less

Submitted 14 April, 2022; originally announced April 2022.

Electronic correlations and universal long-range scaling in kagome metals

Authors: Domenico Di Sante, Bongjae Kim, Werner Hanke, Tim Wehling, Cesare Franchini, Ronny Thomale, Giorgio Sangiovanni

Abstract: We investigate the real-space profile of effective Coulomb interactions in correlated kagome materials. By particularizing to KV$_3$Sb$_5$, Co$_3$Sn$_2$S$_2$, FeSn, and Ni$_3$In, we analyze representative cases that exhibit a large span of correlation-mediated phenomena, and contrast them to prototypical prevoskite transition metal oxides. From our constrained random phase approximation studies we… ▽ More We investigate the real-space profile of effective Coulomb interactions in correlated kagome materials. By particularizing to KV$_3$Sb$_5$, Co$_3$Sn$_2$S$_2$, FeSn, and Ni$_3$In, we analyze representative cases that exhibit a large span of correlation-mediated phenomena, and contrast them to prototypical prevoskite transition metal oxides. From our constrained random phase approximation studies we find that the on-site interaction strength in kagome metals not only depends on the screening processes at high energy, but also on the low-energy hybriziation profile of the electronic density of states. Our results indicate that rescaled by the onsite interaction amplitude, all kagome metals exhibit a universal long-range Coulomb behaviour. △ Less

Submitted 21 March, 2023; v1 submitted 9 March, 2022; originally announced March 2022.

Comments: 7 pages, 3 figures

Journal ref: Phys. Rev. Research 5, L012008 (2023)

Deep Learning the Functional Renormalization Group

Authors: Domenico Di Sante, Matija Medvidović, Alessandro Toschi, Giorgio Sangiovanni, Cesare Franchini, Anirvan M. Sengupta, Andrew J. Millis

Abstract: We perform a data-driven dimensionality reduction of the scale-dependent 4-point vertex function characterizing the functional Renormalization Group (fRG) flow for the widely studied two-dimensional $t - t'$ Hubbard model on the square lattice. We demonstrate that a deep learning architecture based on a Neural Ordinary Differential Equation solver in a low-dimensional latent space efficiently lear… ▽ More We perform a data-driven dimensionality reduction of the scale-dependent 4-point vertex function characterizing the functional Renormalization Group (fRG) flow for the widely studied two-dimensional $t - t'$ Hubbard model on the square lattice. We demonstrate that a deep learning architecture based on a Neural Ordinary Differential Equation solver in a low-dimensional latent space efficiently learns the fRG dynamics that delineates the various magnetic and $d$-wave superconducting regimes of the Hubbard model. We further present a Dynamic Mode Decomposition analysis that confirms that a small number of modes are indeed sufficient to capture the fRG dynamics. Our work demonstrates the possibility of using artificial intelligence to extract compact representations of the 4-point vertex functions for correlated electrons, a goal of utmost importance for the success of cutting-edge quantum field theoretical methods for tackling the many-electron problem. △ Less

Submitted 28 March, 2023; v1 submitted 26 February, 2022; originally announced February 2022.

Comments: 6 pages, 5 figures

Journal ref: Phys. Rev. Lett. 129, 136402 (2022)

Modelling Magnetic Multipolar Phases in Density Functional Theory

Authors: Dario Fiore Mosca, Leonid V. Pourovskii, Cesare Franchini

Abstract: Multipolar magnetic phases in correlated insulators represent a great challenge for Density Functional Theory (DFT) due to the coexistence of intermingled interactions, typically spin-orbit coupling, crystal field and complex non-collinear and high-rank inter-site exchange, creating a complected configurational space with multiple minima. Though the +U correction to DFT allows, in principle, the m… ▽ More Multipolar magnetic phases in correlated insulators represent a great challenge for Density Functional Theory (DFT) due to the coexistence of intermingled interactions, typically spin-orbit coupling, crystal field and complex non-collinear and high-rank inter-site exchange, creating a complected configurational space with multiple minima. Though the +U correction to DFT allows, in principle, the modelling of such magnetic ground states, its results strongly depend on the initially symmetry breaking, constraining the nature of order parameter in the converged DFT+U solution. As a rule, DFT+U calculations starting from a set of initial on-site magnetic moments result in a conventional dipolar order. A more sophisticated approach is clearly needed in the case of magnetic multipolar ordering, which is revealed by a null integral of the magnetization density over spheres centered on magnetic atoms, but with non-zero local contributions. Here we show how such phases can be efficiently captured using an educated constrained initialisation of the onsite density matrix, which is derived from the multipolar-ordered ground state of an ab initio effective Hamiltonian. Various properties of such exotic ground states, like their one-electron spectra, become therefore accessible by all-electron DFT+U methods. We assess the reliability of this procedure on the Ferro-Octupolar ground state recently predicted in Ba$_2$MOsO$_6$ (M = Ca, Mg, Zn) [Phys. Rev. Lett. 127, 237201 (2021)] △ Less

Submitted 25 February, 2022; originally announced February 2022.

Machine Learning for Exploring Small Polaron Configurational Space

Authors: Viktor C. Birschitzky, Florian Ellinger, Ulrike Diebold, Michele Reticcioli, Cesare Franchini

Abstract: Polaron defects are ubiquitous in materials and play an important role in many processes involving carrier mobility, charge transfer and surface reactivity. Determining the spatial distribution of small polarons is essential to understand materials properties and functionalities. This requires an exploration of the configurational space, which is computationally demanding when using standard first… ▽ More Polaron defects are ubiquitous in materials and play an important role in many processes involving carrier mobility, charge transfer and surface reactivity. Determining the spatial distribution of small polarons is essential to understand materials properties and functionalities. This requires an exploration of the configurational space, which is computationally demanding when using standard first principles methods, and technically prohibitive for many-polaron systems. Here, we propose a machine-learning (ML) accelerated search that compares the energy stability of different polaron patterns and determines the ground state configuration. The kernel-regression based ML model is trained on databases generated by density functional theory (DFT) calculations on a minimal set of initial polaron patterns, obtained by using either molecular dynamics simulations or a random sampling approach. To establish an efficient mapping between training data and configuration stability we designed simple descriptors that model the interactions among polarons and charged point defects. The proposed DFT+ML protocol is used here to explore millions of polaron configurations for two different systems, oxygen defective rutile TiO$_2$(110) and Nb-doped SrTiO$_3$(001). Our data shows that the ML-aided search correctly individuates the ground-state polaron patterns, proposes polaronic configurations not visited in the training and can be used to efficiently determine the optimal distribution of polarons at any charge concentration. △ Less

Submitted 2 February, 2022; originally announced February 2022.

Unravelling CO adsorption on model single-atom catalysts

Authors: Jan Hulva, Matthias Meier, Roland Bliem, Zdenek Jakub, Florian Kraushofer, Michael Schmid, Ulrike Diebold, Cesare Franchini, Gareth S. Parkinson

Abstract: Understanding how the local environment of a single-atom catalyst affects stability and reactivity remains a significant challenge. We present an in-depth study of Cu1, Ag1, Au1, Ni1, Pd1, Pt1, Rh1, and Ir1 species on Fe3O4(001); a model support where all metals occupy the same 2-fold coordinated adsorption site upon deposition at room temperature. Surface science techniques revealed that CO adsor… ▽ More Understanding how the local environment of a single-atom catalyst affects stability and reactivity remains a significant challenge. We present an in-depth study of Cu1, Ag1, Au1, Ni1, Pd1, Pt1, Rh1, and Ir1 species on Fe3O4(001); a model support where all metals occupy the same 2-fold coordinated adsorption site upon deposition at room temperature. Surface science techniques revealed that CO adsorption strength at single metal sites differs from the respective metal surfaces and supported clusters. Charge transfer into the support modifies the d-states of the metal atom and the strength of the metal-CO bond. These effects could strengthen the bond (as for Ag1-CO) or weaken it (as for Ni1-CO), but CO-induced structural distortions reduce adsorption energies from those expected based on electronic structure alone. The extent of the relaxations depends on the local geometry and could be predicted by analogy to coordination chemistry. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Journal ref: Science 371, 375-379 (2021)

Ferro-octupolar order and low-energy excitations in d$^2$ double perovskites of Osmium

Authors: Leonid V. Pourovskii, Dario Fiore Mosca, Cesare Franchini

Abstract: Conflicting interpretations of experimental data preclude the understanding of the quantum magnetic state of spin-orbit coupled $d^2$ double perovskites. Whether the ground state is a Janh-Teller-distorted order of quadrupoles or the hitherto elusive octupolar order remains debated. We resolve this uncertainty through direct calculations of all-rank inter-site exchange interactions and inelastic n… ▽ More Conflicting interpretations of experimental data preclude the understanding of the quantum magnetic state of spin-orbit coupled $d^2$ double perovskites. Whether the ground state is a Janh-Teller-distorted order of quadrupoles or the hitherto elusive octupolar order remains debated. We resolve this uncertainty through direct calculations of all-rank inter-site exchange interactions and inelastic neutron scattering cross-section for the $d^2$ double perovskite series Ba$_2M$OsO$_6$ ($M$= Ca, Mg, Zn). Using advanced many-body first principles methods we show that the ground state is formed by ferro-ordered octupoles coupled by superexchange interactions within the ground-state $E_g$ doublet. Computed ordering temperature of the single second-order phase-transition is consistent with experimentally observed material-dependent trends. Minuscule distortions of the parent cubic structure are shown to qualitatively modify the structure of gaped magnetic excitations. △ Less

Submitted 1 December, 2021; v1 submitted 9 July, 2021; originally announced July 2021.

Comments: 6 pages + 12 pages of Supplementary. Final published version

Journal ref: Phys. Rev. Lett. 127, 237201 (2021)

Optical and excitonic properties of transition metal oxide perovskites by the Bethe-Salpeter equation

Authors: Lorenzo Varrassi, Peitao Liu, Zeynep Ergönenc Yavas, Menno Bokdam, Georg Kresse, Cesare Franchini

Abstract: We present a systematic investigation of the role and importance of excitonic effects on the optical properties of transitions metal oxide perovskites. A representative set of fourteen compounds has been selected, including 3$d$ (SrTiO$_3$, LaScO$_3$, LaTiO$_3$, LaVO$_3$, LaCrO$_3$, LaMnO$_3$, LaFeO$_3$ and SrMnO$_3$), 4$d$ (SrZrO$_3$, SrTcO$_3$ and Ca$_2$RuO$_4$) and 5$d$ (SrHfO$_3$, KTaO$_3$ and… ▽ More We present a systematic investigation of the role and importance of excitonic effects on the optical properties of transitions metal oxide perovskites. A representative set of fourteen compounds has been selected, including 3$d$ (SrTiO$_3$, LaScO$_3$, LaTiO$_3$, LaVO$_3$, LaCrO$_3$, LaMnO$_3$, LaFeO$_3$ and SrMnO$_3$), 4$d$ (SrZrO$_3$, SrTcO$_3$ and Ca$_2$RuO$_4$) and 5$d$ (SrHfO$_3$, KTaO$_3$ and NaOsO$_3$) perovskites, covering a band gap ranging from 0.1 eV to 6.1 eV and exhibiting different electronic, structural and magnetic properties. Optical conductivities and optical transitions including electron-hole interactions are calculated through the solution of the Bethe-Salpeter equation (BSE) with quasi-particle energies evaluated by single-shot $G_0W_0$ approximation. The exciton binding energies are computed by means of a model-BSE (mBSE), carefully benchmarked against the full BSE method, in order to obtain well-converged results in terms of k-point sampling. The predicted results are compared with available measured data, with an overall satisfactory agreement between theory and experiment. △ Less

Submitted 28 June, 2021; originally announced June 2021.

Pressure-induced excitations in the out-of-plane optical response of the nodal-line semimetal ZrSiS

Authors: J. Ebad-Allah, S. Rojewski, M. Vöst, G. Eickerling, W. Scherer, E. Uykur, Raman Sankar, L. Varrassi, C. Franchini, K. Ahn, J. Kuneš, C. A. Kuntscher

Abstract: The anisotropic optical response of the layered, nodal-line semimetal ZrSiS at ambient and high pressure is investigated by frequency-dependent reflectivity measurements for the polarization along and perpendicular to the layers. The highly anisotropic optical conductivity is in very good agreement with results from density functional theory calculations and confirms the anisotropic character of Z… ▽ More The anisotropic optical response of the layered, nodal-line semimetal ZrSiS at ambient and high pressure is investigated by frequency-dependent reflectivity measurements for the polarization along and perpendicular to the layers. The highly anisotropic optical conductivity is in very good agreement with results from density functional theory calculations and confirms the anisotropic character of ZrSiS. Whereas the in-plane optical conductivity shows only modest pressure-induced changes, we found strong effects on the out-of-plane optical conductivity spectrum of ZrSiS, with the appearance of two prominent excitations. These pronounced pressure-induced effects can neither be attributed to a structural phase transition according to our single-crystal x-ray diffraction measurements, nor can they be explained by electronic correlation and electron-hole pairing effects, as revealed by theoretical calculations. Our findings are discussed in the context of the recently proposed excitonic insulator phase in ZrSiS. △ Less

Submitted 22 June, 2021; originally announced June 2021.

Comments: 6 pages, 3 figures, accepted for publication in Phys. Rev. Lett

Diagrammatic quantum Monte Carlo study of an acoustic lattice polaron

Authors: Thomas Hahn, Naoto Nagaosa, Cesare Franchini, Andrey S. Mishchenko

Abstract: We present the first approximation free diagrammatic Monte Carlo study of a lattice polaron interacting with an acoustic phonon branch through the deformation potential. Weak and strong coupling regimes are separated by a self-trapping region where quantum resonance between various possible lattice deformations is seen in the ground state properties, spectral function, and optical conductivity. Th… ▽ More We present the first approximation free diagrammatic Monte Carlo study of a lattice polaron interacting with an acoustic phonon branch through the deformation potential. Weak and strong coupling regimes are separated by a self-trapping region where quantum resonance between various possible lattice deformations is seen in the ground state properties, spectral function, and optical conductivity. The unique feature of such polaron is the interplay between long- and short wavelength acoustic vibrations creating a composite phonon cloud and leading to persistent self-trapping due to the existence of multiple quasi-stable states. This results in a spectral response whose structure is much more complex than in any of the previously considered polaron models. △ Less

Submitted 27 April, 2021; originally announced April 2021.

Comments: 6 pages, 4 figures + supplemental material

Journal ref: Phys. Rev. B 104, L161111 (2021)

Electronic state unfolding for plane waves: energy bands, Fermi surfaces and spectral functions

Authors: David Dirnberger, Georg Kresse, Cesare Franchini, Michele Reticcioli

Abstract: Modern computing facilities grant access to first-principles density-functional theory study of complex physical and chemical phenomena in materials, that require large supercell to properly model the system. However, supercells are associated to small Brillouin zones in the reciprocal space, leading to folded electronic eigenstates that make the analysis and interpretation extremely challenging.… ▽ More Modern computing facilities grant access to first-principles density-functional theory study of complex physical and chemical phenomena in materials, that require large supercell to properly model the system. However, supercells are associated to small Brillouin zones in the reciprocal space, leading to folded electronic eigenstates that make the analysis and interpretation extremely challenging. Various techniques have been proposed and developed in order to reconstruct the electronic band structures of super cells, unfolded into the reciprocal space of an ideal primitive cell. Here, we propose an efficient unfolding scheme embedded directly in the Vienna Ab-initio Simulation Package (VASP), that requires modest computational resources and allows for an automatized mapping from the reciprocal space of the supercell to primitive cell Brillouin zone. This algorithm can computes band structures, Fermi surfaces and spectral functions, by using an integrated post-processing tool (bands4vasp). The method is here applied to a selected variety of complex physical situations: the effect of doping on the band dispersion in the BaFe$_{\rm 2(1-x)}$Ru$_{\rm 2x}$As$_2$ superconductor, the interaction between adsorbates and polaronic states on the TiO$_2$(110) surface, and the band splitting induced by non-collinear spin fluctuations in EuCd$_2$As$_2$. △ Less

Submitted 26 May, 2021; v1 submitted 17 March, 2021; originally announced March 2021.

Interplay between multipolar spin interactions, Jahn-Teller effect and electronic correlation in a $J_{eff}=\frac{3}{2}$ insulator

Authors: Dario Fiore Mosca, Leonid V. Pourovskii, Beom Hyun Kim, Peitao Liu, Samuele Sanna, Federico Boscherini, Sergii Khmelevskyi, Cesare Franchini

Abstract: In this work we study the complex entanglement between spin interactions, electron correlation and Janh-Teller structural instabilities in the 5d$^1$ $J_{eff}=\frac{3}{2}$ spin-orbit coupled double perovskite $\rm Ba_2NaOsO_6$ using first principles approaches. By combining non-collinear magnetic calculations with multipolar pseudospin Hamiltonian analysis and many-body techniques we elucidate the… ▽ More In this work we study the complex entanglement between spin interactions, electron correlation and Janh-Teller structural instabilities in the 5d$^1$ $J_{eff}=\frac{3}{2}$ spin-orbit coupled double perovskite $\rm Ba_2NaOsO_6$ using first principles approaches. By combining non-collinear magnetic calculations with multipolar pseudospin Hamiltonian analysis and many-body techniques we elucidate the origin of the observed quadrupolar canted antifferomagnetic. We show that the non-collinear magnetic order originates from Jahn-Teller distortions due to the cooperation of Heisenberg exchange, quadrupolar spin-spin terms and both dipolar and multipolar Dzyaloshinskii-Moriya interactions. We find a strong competition between ferromagnetic and antiferromagnetic canted and collinear quadrupolar magnetic phases: the transition from one magnetic order to another can be controlled by the strength of the electronic correlation ($U$) and by the degree of Jahn-Teller distortions. △ Less

Submitted 22 February, 2021; originally announced February 2021.

Comparative ab initio study of the structural, electronic, magnetic, and dynamical properties of LiOsO$_3$ and NaOsO$_3$

Authors: Peitao Liu, Jiangang He, Bongjae Kim, Sergii Khmelevskyi, Alessandro Toschi, Georg Kresse, Cesare Franchini

Abstract: Despite similar chemical compositions, LiOsO$_3$ and NaOsO$_3$ exhibit remarkably distinct structural, electronic, magnetic, and spectroscopic properties. At low temperature, LiOsO$_3$ is a polar bad metal with a rhombohedral $R3c$ structure without the presence of long-range magnetic order, whereas NaOsO$_3$ is a $G$-type antiferromagnetic insulator with an orthorhombic $Pnma$ structure. By means… ▽ More Despite similar chemical compositions, LiOsO$_3$ and NaOsO$_3$ exhibit remarkably distinct structural, electronic, magnetic, and spectroscopic properties. At low temperature, LiOsO$_3$ is a polar bad metal with a rhombohedral $R3c$ structure without the presence of long-range magnetic order, whereas NaOsO$_3$ is a $G$-type antiferromagnetic insulator with an orthorhombic $Pnma$ structure. By means of comparative first-principles DFT+$U$ calculations with the inclusion of the spin-orbit coupling, we ($i$) identify the origin of the different structural ($R3c$ vs. $Pnma$) properties using a symmetry-adapted soft mode analysis, ($ii$) provide evidence that all considered exchange-correlation functionals (LDA, PBE, PBEsol, SCAN, and HSE06) and the spin disordered polymorphous descriptions are unsatisfactory to accurately describe the electronic and magnetic properties of both systems simultaneously, and ($iii$) clarify that the distinct electronic (metallic vs. insulating) properties originates mainly from a cooperative steric and magnetic effect. Finally, we find that although at ambient pressure LiOsO$_3$ with a $Pnma$ symmetry and NaOsO$_3$ with a $R\bar{3}c$ symmetry are energetically unfavorable, they do not show soft phonons and therefore are dynamically stable. A pressure-induced structural phase transition from $R3c$ to $Pnma$ for LiOsO$_3$ is predicted, whereas for NaOsO$_3$ no symmetry change is discerned in the considered pressure range. △ Less

Submitted 14 February, 2020; originally announced February 2020.

Comments: 9 pages, 9 figures

Journal ref: Phys. Rev. Materials 4, 045001 (2020)

SrRuO$_3$-SrTiO$_3$ heterostructure as a possible platform for studying unconventional superconductivity in Sr$_{2}$RuO$_{4}$

Authors: Bongjae Kim, Sergii Khmelevskyi, Cesare Franchini, I. I. Mazin, Kyoo Kim

Abstract: There is intense controversy around the unconventional superconductivity in strontium ruthenate, where the various theoretical and experimental studies suggest diverse and mutually exclusive pairing symmetries. Currently, the investigation is solely focused on only one material, Sr2RuO4, and the field suffers from the lack of comparison targets. Here, employing a density functional theory based an… ▽ More There is intense controversy around the unconventional superconductivity in strontium ruthenate, where the various theoretical and experimental studies suggest diverse and mutually exclusive pairing symmetries. Currently, the investigation is solely focused on only one material, Sr2RuO4, and the field suffers from the lack of comparison targets. Here, employing a density functional theory based analysis, we show that the heterostructure composed of SrRuO3 and SrTiO3 is endowed with all the key characteristics of Sr2RuO4, and, in principle, can host superconductivity. Furthermore, we show that competing magnetic phases and associated frustration, naturally affecting the superconducting state, can be tuned by epitaxial strain engineering. This system thus offers an excellent platform for gaining more insight into superconductivity in ruthenates. △ Less

Submitted 3 June, 2020; v1 submitted 27 January, 2020; originally announced January 2020.

Journal ref: Phys. Rev. B 101, 220502 (2020)

Electron-Phonon Interactions Using the PAW Method and Wannier Functions

Authors: Manuel Engel, Martijn Marsman, Cesare Franchini, Georg Kresse

Abstract: We present an ab-initio density-functional-theory approach for calculating electron-phonon interactions within the projector augmented-wave method. The required electron-phonon matrix elements are defined as the second derivative of the one-electron energies in the PAW method. As the PAW method leads to a generalized eigenvalue problem, the resulting electron-phonon matrix elements lack some symme… ▽ More We present an ab-initio density-functional-theory approach for calculating electron-phonon interactions within the projector augmented-wave method. The required electron-phonon matrix elements are defined as the second derivative of the one-electron energies in the PAW method. As the PAW method leads to a generalized eigenvalue problem, the resulting electron-phonon matrix elements lack some symmetries that are usually present for simple eigenvalue problems and all-electron formulations. We discuss the relation between our definition of the electron-phonon matrix element and other formulations. To allow for efficient evaluation of physical properties, we introduce a Wannier-interpolation scheme, again adapted to generalized eigenvalue problems. To explore the method's numerical characteristics, the temperature-dependent band-gap renormalization of diamond is calculated and compared with previous publications. Furthermore, we apply the method to selected binary compounds and show that the obtained zero-point renormalizations agree well with other values found in literature and experiments. △ Less

Submitted 23 January, 2020; originally announced January 2020.

Journal ref: Phys. Rev. B 101, 184302 (2020)

Optical response of an interacting polaron gas in strongly polar crystals

Authors: S. N. Klimin, J. Tempere, J. T. Devreese, C. Franchini, G. Kresse

Abstract: Optical conductivity of an interacting polaron gas is calculated within an extended random phase approximation which takes into account mixing of collective excitations of the electron gas with LO phonons. This mixing is important for the optical response of strongly polar crystals where the static dielectric constant is rather high: strontium titanate is the case. The present calculation sheds li… ▽ More Optical conductivity of an interacting polaron gas is calculated within an extended random phase approximation which takes into account mixing of collective excitations of the electron gas with LO phonons. This mixing is important for the optical response of strongly polar crystals where the static dielectric constant is rather high: strontium titanate is the case. The present calculation sheds light on unexplained features of experimentally observed optical conductivity spectra in $n$-doped SrTiO$_{3}$. These features appear to be due to dynamic screening of the electron-electron interaction by polar optical phonons and hence do not require additional mechanisms for the explanation. △ Less

Submitted 6 December, 2019; originally announced December 2019.

Comments: 8 pages, 1 figure

Osmates on the verge of a Hund's-Mott transition: The different fates of NaOsO$_3$ and LiOsO$_3$

Authors: Daniel Springer, Bongjae Kim, Peitao Liu, Sergii Khmelevskyi, Severino Adler, Massimo Capone, Giorgio Sangiovanni, Cesare Franchini, Alessandro Toschi

Abstract: We clarify the origin of the strikingly different spectroscopic properties of the chemically similar compounds NaOsO$_3$ and LiOsO$_3$. Our first-principle many-body analysis demonstrates that the highly sensitive physics of these two materials is controlled by their proximity to an adjacent Hund's-Mott insulating phase. In mildly correlated $5d$ oxides, this physics is triggered by the {\sl coope… ▽ More We clarify the origin of the strikingly different spectroscopic properties of the chemically similar compounds NaOsO$_3$ and LiOsO$_3$. Our first-principle many-body analysis demonstrates that the highly sensitive physics of these two materials is controlled by their proximity to an adjacent Hund's-Mott insulating phase. In mildly correlated $5d$ oxides, this physics is triggered by the {\sl cooperative} action of intraorbital repulsion and Hund's exchange in half-filled $t_{2g}$ outer-shells. Small material-specific details hence result in an extremely sharp change of the electronic mobility, explaining the surprisingly different properties of the paramagnetic high-temperature phases of the two compounds. △ Less

Submitted 24 September, 2020; v1 submitted 11 October, 2019; originally announced October 2019.

Journal ref: Phys. Rev. Lett. 125, 166402 (2020)