Dynamic Local Structure in Caesium Lead Iodide: Spatial Correlation and Transient Domains.

Metal halide perovskites are multifunctional semiconductors with tunable structures and properties. They are highly dynamic crystals with complex octahedral tilting patterns and strongly anharmonic atomic behavior. In the higher temperature, higher symmetry phases of these materials, several complex structural features are observed. The local structure can differ greatly from the average structure and there is evidence that dynamic 2D structures of correlated octahedral motion form. An understanding of the underlying complex atomistic dynamics is, however, still lacking. In this work, the local structure of the inorganic perovskite CsPbI3 is investigated using a new machine learning force field based on the atomic cluster expansion framework. Through analysis of the temporal and spatial correlation observed during large-scale simulations, it is revealed that the low frequency motion of octahedral tilts implies a double-well effective potential landscape, even well into the cubic phase. Moreover, dynamic local regions of lower symmetry are present within both higher symmetry phases. These regions are planar and the length and timescales of the motion are reported. Finally, the spatial arrangement of these features and their interactions are investigated and visualized, providing a comprehensive picture of local structure in the higher symmetry phases.

Direct Correlation of Nanoscale Morphology and Device Performance to Study Photocurrent Generation in Donor-Enriched Phases of Polymer Solar Cells.

The nanoscale morphology of polymer blends is a key parameter to reach high efficiency in bulk heterojunction solar cells. Thereby, research typically focusing on optimal blend morphologies while studying nonoptimized blends may give insight into blend designs that can prove more robust against morphology defects. Here, we focus on the direct correlation of morphology and device performance of thieno[3,4-b]-thiophene-alt-benzodithiophene (PTB7):[6,6]phenyl C71 butyric acid methyl ester (PC71BM) bulk heterojunction (BHJ) blends processed without additives in different donor/acceptor weight ratios. We show that while blends of a 1:1.5 ratio are composed of large donor-enriched and fullerene domains beyond the exciton diffusion length, reducing the ratio below 1:0.5 leads to blends composed purely of polymer-enriched domains. Importantly, the photocurrent density in such blends can reach values between 45 and 60% of those reached for fully optimized blends using additives. We provide here direct visual evidence that fullerenes in the donor-enriched domains are not distributed homogeneously but fluctuate locally. To this end, we performed compositional nanoscale morphology analysis of the blend using spectroscopic imaging of low-energy-loss electrons using a transmission electron microscope. Charge transport measurement in combination with molecular dynamics simulations shows that the fullerene substructures inside the polymer phase generate efficient electron transport in the polymer-enriched phase. Furthermore, we show that the formation of densely packed regions of fullerene inside the polymer phase is driven by the PTB7:PC71BM enthalpy of mixing. The occurrence of such a nanoscale network of fullerene clusters leads to a reduction of electron trap states and thus efficient extraction of photocurrent inside the polymer domain. Suitable tuning of the polymer-acceptor interaction can thus introduce acceptor subnetworks in polymer-enriched phases, improving the tolerance for high-efficiency BHJ toward morphological defects such as donor-enriched domains exceeding the exciton diffusion length.

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.

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.

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.

Tuning the thermal conductivity of methylammonium lead halide by the molecular substructure.

By using state-of-the-art atomistic methods we provide an accurate estimate of thermal conductivity of methylammonium lead halide as a function of sample size and temperature, in agreement with experimental works. We show that the thermal conductivity of methylammonium lead halide is intrinsically low, due to the low sound velocity of the PbI lattice. Furthermore, by selectively analyzing the effect of different molecular degrees of freedom, we clarify the role of the molecular substructure by showing that the internal modes above 150 cm(-1) (in addition to rotations) are effective in reducing the thermal conductivity of hybrid perovskites. This analysis suggests strategies to tailor the thermal conductivity by modifying the internal structure of organic cations.

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.

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).