Paolo Umari

QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials

QUANTUM ESPRESSO is an integrated suite of computer codes for electronic-structure calculations and materials modeling, based on density-functional theory, plane waves, and pseudopotentials (norm-conserving, ultrasoft, and projector-augmented wave). The acronym ESPRESSO stands for opEn Source Package for Research in Electronic Structure, Simulation, and Optimization. It is freely available to researchers around the world under the terms of the GNU General Public License. QUANTUM ESPRESSO builds upon newly-restructured electronic-structure codes that have been developed and tested by some of the original authors of novel electronic-structure algorithms and applied in the last twenty years by some of the leading materials modeling groups worldwide. Innovation and efficiency are still its main focus, with special attention paid to massively parallel architectures, and a great effort being devoted to user friendliness. QUANTUM ESPRESSO is evolving towards a distribution of independent and interoperable codes in the spirit of an open-source project, where researchers active in the field of electronic-structure calculations are encouraged to participate in the project by contributing their own codes or by implementing their own ideas into existing codes.

Relativistic GW calculations on CH3NH3PbI3 and CH3NH3SnI3 perovskites for solar cell applications.

Hybrid AMX3 perovskites (A=Cs, CH3NH3; M=Sn, Pb; X=halide) have revolutionized the scenario of emerging photovoltaic technologies. Introduced in 2009 by Kojima et al., a rapid evolution very recently led to 15% efficient solar cells. CH3NH3PbI3 has so far dominated the field, while the similar CH3NH3SnI3 has not been explored for photovoltaic applications, despite the reduced band-gap. Replacement of Pb by the more environment-friendly Sn would facilitate the large uptake of perovskite-based photovoltaics. Despite the extremely fast progress, the materials electronic properties which are key to the photovoltaic performance are relatively little understood. Here we develop an effective GW method incorporating spin-orbit coupling which allows us to accurately model the electronic, optical and transport properties of CH3NH3SnI3 and CH3NH3PbI3, opening the way to new materials design. The different CH3NH3SnI3 and CH3NH3PbI3 properties are discussed in light of their exploitation for solar cells, and found to be entirely due to relativistic effects.Hybrid AMX3 perovskites (A = Cs, CH3NH3; M = Sn, Pb; X = halide) have revolutionized the scenario of emerging photovoltaic technologies, with very recent results demonstrating 15% efficient solar cells. The CH3NH3PbI3/MAPb(I1−xClx)3 perovskites have dominated the field, while the similar CH3NH3SnI3 has not been exploited for photovoltaic applications. Replacement of Pb by Sn would facilitate the large uptake of perovskite-based photovoltaics. Despite the extremely fast progress, the materials electronic properties which are key to the photovoltaic performance are relatively little understood. Density Functional Theory electronic structure methods have so far delivered an unbalanced description of Pb- and Sn-based perovskites. Here we develop an effective GW method incorporating spin-orbit coupling which allows us to accurately model the electronic, optical and transport properties of CH3NH3SnI3 and CH3NH3PbI3, opening the way to new materials design. The different CH3NH3SnI3 and CH3NH3PbI3 electronic properties are discussed in light of their exploitation for solar cells, and found to be dominantly due to relativistic effects. These effects stabilize the CH3NH3PbI3 material towards oxidation, by inducing a deeper valence band edge. Relativistic effects, however, also increase the material band-gap compared to CH3NH3SnI3, due to the valence band energy downshift (~0.7 eV) being only partly compensated by the conduction band downshift (~0.2 eV).

Cation-Induced Band-Gap Tuning in Organohalide Perovskites: Interplay of Spin-Orbit Coupling and Octahedra Tilting

Organohalide lead perovskites have revolutionized the scenario of emerging photovoltaic technologies. The prototype MAPbI3 perovskite (MA = CH3NH3(+)) has dominated the field, despite only harvesting photons above 750 nm (∼1.6 eV). Intensive research efforts are being devoted to find new perovskites with red-shifted absorption onset, along with good charge transport properties. Recently, a new perovskite based on the formamidinium cation ((NH2)2CH(+) = FA) has shown potentially superior properties in terms of band gap and charge transport compared to MAPbI3. The results have been interpreted in terms of the cation size, with the larger FA cation expectedly delivering reduced band-gaps in Pb-based perovskites. To provide a full understanding of the interplay among size, structure, and organic/inorganic interactions in determining the properties of APbI3 perovskites, in view of designing new materials and fully exploiting them for solar cells applications, we report a fully first-principles investigation on APbI3 perovskites with A = Cs(+), MA, and FA. Our results evidence that the tetragonal-to-quasi cubic structural evolution observed when moving from MA to FA is due to the interplay of size effects and enhanced hydrogen bonding between the FA cations and the inorganic matrix altering the covalent/ionic character of Pb-I bonds. Most notably, the observed cation-induced structural variability promotes markedly different electronic and optical properties in the MAPbI3 and FAPbI3 perovskites, mediated by the different spin-orbit coupling, leading to improved charge transport and red-shifted absorption in FAPbI3 and in general in pseudocubic structures. Our theoretical model constitutes the basis for the rationale design of new and more efficient organohalide perovskites for solar cells applications.

Advanced capabilities for materials modelling with Quantum ESPRESSO

Quantum ESPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudo-potential and projector-augmented-wave approaches. Quantum ESPRESSO owes its popularity to the wide variety of properties and processes it allows to simulate, to its performance on an increasingly broad array of hardware architectures, and to a community of researchers that rely on its capabilities as a core open-source development platform to implement theirs ideas. In this paper we describe recent extensions and improvements, covering new methodologies and property calculators, improved parallelization, code modularization, and extended interoperability both within the distribution and with external software.Quantum EXPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the-art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudopotential and projector-augmented-wave approaches. Quantum EXPRESSO owes its popularity to the wide variety of properties and processes it allows to simulate, to its performance on an increasingly broad array of hardware architectures, and to a community of researchers that rely on its capabilities as a core open-source development platform to implement their ideas. In this paper we describe recent extensions and improvements, covering new methodologies and property calculators, improved parallelization, code modularization, and extended interoperability both within the distribution and with external software.

GW quasiparticle spectra from occupied states only

We introduce a method that allows for the calculation of quasi-particle spectra in the GW approximation, yet avoiding any explicit reference to empty one-electron states. This is achieved by expressing the irreducible polarizability operator and the self-energy operator through a set of linear response equations, which are solved using a Lanczos-chain algorithm. We first validate our approach by calculating the vertical ionization energies of the benzene molecule and then show its potential by addressing the spectrum of a large molecule such as free-base tetraphenylporphyrin.

Electronic and optical properties of mixed Sn–Pb organohalide perovskites: a first principles investigation

Organohalide lead perovskites have attracted considerable interest among emerging photovoltaic technologies, delivering highly efficient solid-state solar cells. Despite the huge potential of this class of materials, the use of Pb-containing systems will likely hamper the wide spread take off of perovskite solar cells. The development of lead-free hybrid perovskites represents thus an important step in the exploitation of this promising technology. Very recently, the use of new mixed Pb–Sn MASnxPb(1−x)I3 perovskites has been reported, with a considerable increase in the extension of the solar spectrum absorption with respect to lead-halide perovskites, shifting the absorption onset down to the near-IR. In light of the anticipated potential of Sn-based organohalide perovskites in replacing lead-based materials, here we apply a recently developed computational approach to the description of mixed Sn–Pb compounds showing a range of compositions comparable to experimentally characterized compounds. For the investigated series of MASnxPb(1−x)I3 perovskites we find a continuous and monotonic variation of the energy levels, shifting at lower potentials; and band-gaps, which shift towards the near-IR, as the Sn content in the perovskite is increased. Notably, while we find slightly unbalanced electron/hole transport in the pure phases, Pb (Sn) materials being better electron (hole) transporters, for intermediate compositions an almost perfectly balanced charge carrier transport can be achieved, in line with recent experimental observations.

Electronic and optical properties of MAPbX3 perovskites (X = I, Br, Cl): a unified DFT and GW theoretical analysis

Materials engineering is a key for the enhancement of photovoltaics technology. This is particularly true for the novel class of perovskite solar cells. Accurate theoretical modelling can help establish general trends of behavior when addressing structural changes. Here, we consider the effects due to halide substitution in organohalide CH3NH3PbX3 perovskites exploring the halide series with X = Cl, Br, I. For this task, we use accurate DFT and GW methods including spin-orbit coupling. We find the expected band gap increase when moving from X = I to Cl, in line with the experimental data. Most notably, the calculated absorption coefficients for I, Br and Cl are nicely reproducing the behavior reported experimentally. A common feature of all the simulated band structures is a significant Rashba effect. This is similar for MAPbI3 and MAPbBr3 while MAPbCl3 shows in general a reduced Rashba interaction coefficient. Finally, a monotonic increase of the exciton reduced masses is calculated when moving from I to Br to Cl, in line with the stronger excitonic character of the lighter perovskite halides.

Communication: electronic band gaps of semiconducting zig-zag carbon nanotubes from many-body perturbation theory calculations.

Electronic band gaps for optically allowed transitions are calculated for a series of semiconducting single-walled zig-zag carbon nanotubes of increasing diameter within the many-body perturbation theory GW method. The dependence of the evaluated gaps with respect to tube diameters is then compared with those found from previous experimental data for optical gaps combined with theoretical estimations of exciton binding energies. We find that our GW gaps confirm the behavior inferred from experiment. The relationship between the electronic gap and the diameter extrapolated from the GW values is also in excellent agreement with a direct measurement recently performed through scanning tunneling spectroscopy.

Modeling of the Raman spectrum of vitreous silica: concentration of small ring structures

To calculate the Raman spectrum of vitreous silica, we use a model structure consisting of a disordered network of cornersharing tetrahedra, for which the vibrational frequencies and eigenmodes were obtained previously by first principles. The Raman cross-section is evaluated using a bond polarizability model with parameters derived from first-principles results for a-quartz. This model reproduces the major features of the spectrum. By comparing measured and calculated intensities of the defect lines D-1 and D-2, weestimate that similar to0.2% and similar to3% of oxygen atoms in vitreous silica belong to three- and four-membered rings, respectively

Chlorine Incorporation in the CH3NH3PbI3 Perovskite: Small Concentration, Big Effect

The role of chlorine doping in CH3NH3PbI3 represents an important open issue in the use of hybrid perovskites for photovoltaic applications. In particular, even if a positive role of chlorine doping on perovskite film formation and on material morphology has been demonstrated, an inherent positive effect on the electronic and photovoltaic properties cannot be excluded. Here we carried out periodic density functional theory and Car-Parrinello molecular dynamics simulations, going down to ∼1% doping, to investigate the effect of chlorine on CH3NH3PbI3. We found that such a small doping has important effects on the dynamics of the crystalline structure, both with respect to the inorganic framework and with respect to the cation libration motion. Together, we observe a dynamic spatial localization of the valence and conduction states in separated spatial material regions, which takes place in the 10-1 ps time scale and which could be the key to ease of exciton dissociation and, likely, to small charge recombination in hybrid perovskites. Moreover, such localization is enhanced by chlorine doping, demonstrating an inherent positive role of chlorine doping on the electronic properties of this class of materials.