Claudia Caddeo

Thermally Activated Point Defect Diffusion in Methylammonium Lead Trihalide: Anisotropic and Ultrahigh Mobility of Iodine

We study the diffusion of point defects in crystalline methylammonium lead halide (MAPI) at finite temperatures by using all-atoms molecular dynamics. We find that, for what concerns intrinsic defects, iodine diffusion is by far the dominant mechanism of ionic transport in MAPI, with diffusivities as high as 7.4 × 10(-7) and 4.3 × 10(-6) cm(2) s(-1) at 300 K and single activation energies of 0.24 and 0.10 eV, for interstitials and vacancies, respectively. The comparison with common covalent and oxide crystals reveals the ultrahigh mobility of defects in MAPI. Though at room temperature the vacancies are about 1 order of magnitude more diffusive, the anisotropic interstitial dynamics increases more rapidly with temperature, and it can be dominant at high temperatures. Present results are fully consistent with the involvement of iodide ions in hysteresis and have implications for improvement of the material quality by better control of defect diffusion.

Temperature Evolution of Methylammonium Trihalide Vibrations at the Atomic Scale

The temperature evolution of vibrations of CH3NH3PbI3 (MAPI) is studied by combining first principles and classical molecular dynamics and compared to available experimental data. The work has a fundamental character showing that it is possible to reproduce the key features of the vibrational spectrum by the simple physical quantities included in the classical model, namely the ionic-dispersive hybrid interactions and the mass difference between organic and inorganic components. The dynamics reveals a sizable temperature evolution of the MAPI spectrum along with the orthorhombic-to-tetragonal-to-cubic transformation and a strong dependence on molecular confinement and order. The thermally induced weakening of the H-I interactions and the anharmonic mixing of modes give two vibrational peaks at 200-250 cm(-1) that are not present at zero temperature and are expected to have detectable infrared activity. The infrared inactive vibrational peak at ∼140 cm(-1) due to molecular spinning disappears abruptly at the orthorhombic-to-tetragonal transition and forms a broad molecular band red-shifting progressively with temperature. This trend is correlated to the reduced confinement of the rotating cations due to thermal expansion of the lattice.

Modeling hybrid perovskites by molecular dynamics

The topical review describes the recent progress in the modeling of hybrid perovskites by molecular dynamics simulations. Hybrid perovskites and in particular methylammonium lead halide (MAPI) have a tremendous technological relevance representing the fastest-advancing solar material to date. They also represent the paradigm of an organic-inorganic crystalline material with some conceptual peculiarities: an inorganic semiconductor for what concerns the electronic and absorption properties with a hybrid and solution processable organic-inorganic body. After briefly explaining the basic concepts of ab initio and classical molecular dynamics, the model potential recently developed for hybrid perovskites is described together with its physical motivation as a simple ionic model able to reproduce the main dynamical properties of the material. Advantages and limits of the two strategies (either ab initio or classical) are discussed in comparison with the time and length scales (from pico to microsecond scale) necessary to comprehensively study the relevant properties of hybrid perovskites from molecular reorientations to electrocaloric effects. The state-of-the-art of the molecular dynamics modeling of hybrid perovskites is reviewed by focusing on a selection of showcase applications of methylammonium lead halide: molecular cations disorder; temperature evolution of vibrations; thermally activated defects diffusion; thermal transport. We finally discuss the perspectives in the modeling of hybrid perovskites by molecular dynamics.

Optoelectronic properties of (ZnO)60 isomers

We studied the optoelectronic properties of six possible structures of the (ZnO)(60) cluster using density functional theory (DFT). Vertical ionization energies and electron affinities are calculated through total energy differences, while the optical absorption spectra are obtained by using hybrid time-dependent DFT. The (ZnO)(60) cluster has been proven to be particularly stable and it is of potential interest for future applications in nanoelectronics, but its ground-state configuration has been unknown to date. Since the relative stability inferred from total energy calculations suffers from a strong dependence on the computational scheme adopted, we combined it with optical spectroscopy to identify the most abundant geometrical structure of this cluster. The calculated optical spectra are different for each isomer and they could be thus used in comparison with experimental data to explain the ground state of (ZnO)(60).

Atomistic simulations of P(NDI2OD-T2) morphologies: from single chain to condensed phases.

We investigate theoretically the structure, crystallinity, and solubility of a high-mobility n-type semiconducting copolymer, P(NDI2OD-T2), and we propose a set of new force field parameters. The force field is reparametrized against density functional theory (DFT) calculations, with the aim to reproduce the correct torsional angles that govern the polymer chain flexibility and morphology. We simulate P(NDI2OD-T2) oligomers in different environments, namely, in vacuo, in the bulk phase, and in liquid toluene and chloronaphthalene solution. The choice of these solvents is motivated by the fact that they induce different kinds of molecular preaggregates during the casting procedures, resulting in variable device performances. Our results are in good agreement with the available experimental data; the polymer bulk structure, in which the chains are quite planar, is correcly reproduced, yet the isolated chains are flexible enough to fold in vacuo. We also calculate the solubility of P(NDI2OD-T2) in toluene and chloronaphthalene, predicting a much better solubility of the polymer in the latter, also in accordance to experimental observations. Different morphologies and dynamics of the oligomers in the two solvents have been observed. The proposed parameters make it possible to obtain the description of P(NDI2OD-T2) in different environments and can serve as a basis for extensive studies of this polymer semiconductor, such as, for example, the dynamics of aggregation in solvent.

Collective Molecular Mechanisms in the CH3NH3PbI3 Dissolution by Liquid Water

The origin of the dissolution of methylammonium lead trihalide (MAPI) crystals in liquid water is clarified by finite-temperature molecular dynamics by developing a MYP-based force field (MYP1) for water-MAPI systems. A thermally activated process is found with an energy barrier of 0.36 eV consisting of a layer-by-layer degradation with generation of inorganic PbI2 films and solvation of MA and I ions. We rationalize the effect of water on MAPI by identifying a transition from a reversible absorption and diffusion in the presence of vapor to the irreversible destruction of the crystal lattice in liquid due to a cooperative action of water molecules. A strong water-MAPI interaction is found with a binding energy of 0.41 eV/H2O and wetting energy of 0.23 N/m. The water vapor absorption is energetically favored (0.29 eV/H2O), and the infiltrated molecules can migrate within the crystal with a diffusion coefficient D = 1.7 × 10-8 cm2/s and activation energy of 0.28 eV.

Thermal boundary resistance from transient nanocalorimetry: A multiscale modeling approach

The Thermal Boundary Resistance at the interface between a nanosized Al film and an Al_{2}O_{3} substrate is investigated at an atomistic level. A room temperature value of 1.4 m^{2}K/GW is found. The thermal dynamics occurring in time-resolved thermo-reflectance experiments is then modelled via macro-physics equations upon insertion of the materials parameters obtained from atomistic simulations. Electrons and phonons non-equilibrium and spatio-temporal temperatures inhomo- geneities are found to persist up to the nanosecond time scale. These results question the validity of the commonly adopted lumped thermal capacitance model in interpreting transient nanocalorimetry experiments. The strategy adopted in the literature to extract the Thermal Boundary Resistance from transient reflectivity traces is revised at the light of the present findings. The results are of relevance beyond the specific system, the physical picture being general and readily extendable to other heterojunctions.

Photoacoustic Sensing of Trapped Fluids in Nanoporous Thin Films: Device Engineering and Sensing Scheme

Accessing fluid infiltration in nanogranular coatings is an outstanding challenge, of relevance for applications ranging from nanomedicine to catalysis. A sensing platform, allowing quantifying the amount of fluid infiltrated in a nanogranular ultrathin coating, with thickness in the 10-40 nm range, is here proposed and theoretically investigated by multiscale modeling. The scheme relies on impulsive photoacoustic excitation of hypersonic mechanical breathing modes in engineered gas-phase-synthesized nanogranular metallic ultrathin films and time-resolved acousto-optical read-out of the breathing modes frequency shift upon liquid infiltration. A superior sensitivity, exceeding 26 × 103 cm2/g, is predicted upon equivalent areal mass loading of a few ng/mm2. The capability of the present scheme to discriminate among different infiltration patterns is discussed. The platform is an ideal tool to investigate nanofluidics in granular materials and naturally serves as a distributed nanogetter coating, integrating fluid sensing capabilities. The proposed scheme is readily extendable to other nanoscale and mesoscale porous materials.

Photoluminescence, optical gain, and lasing threshold in CH3NH3PbI3 methylammonium lead-halide perovskites obtained by ab initio calculations

Using ab initio band energy calculations and van Roosbroeck–Shockley recombination theory, we model transient photoluminescence for lead-iodide CH3NH3PbI3 perovskites. We provide clear evidence that the most important features of the photoluminescence process, i.e. strong absorption, low recombination rates, and long lifetimes, can be all coherently derived from the band-to-band recombination process, and the sometimes invoked contradiction between strong absorption and low recombination rate in a direct band-gap material has no reason to exist; optical gain and charge carrier lasing threshold are also reproduced in satisfactory agreement with the observed values, at least within a significant charge-injection range. Our description provides a solid theoretical assessment to the widespread but debated idea that photoluminescence in these materials is predominantly produced by unbound charge recombination.