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.

Remarkably Weak Anisotropy in Thermal Conductivity of Two-Dimensional Hybrid Perovskite Butylammonium Lead Iodide Crystals.

Two-dimensional (2D) hybrid organic-inorganic perovskites consisting of alternating organic and inorganic layers are a new class of layered structures. They have attracted increasing interest for photovoltaic, optoelectronic, and thermoelectric applications, where knowing their thermal transport properties is critical. We carry out both experimental and computational studies on thermal transport properties of 2D butylammonium lead iodide crystals and find their thermal conductivity is ultralow (below 0.3 W m-1 K-1) with very weak anisotropy (around 1.5) among layered crystals. Further analysis reveals that the unique structure with the preferential alignment of organic chains and complicated energy landscape leads to moderately smaller phonon lifetimes in the out-of-plane direction and comparable phonon group velocities in in-plane and out-of-plane directions. These new findings may guide the future design of novel hybrid materials with desired thermal conductivity for various applications.

Long-lived electrets and lack of ferroelectricity in methylammonium lead bromide CH3NH3PbBr3 ferroelastic single crystals.

Hybrid lead halides CH3NH3PbX3 (X = I, Br, and Cl) have emerged as a new class of semiconductors for low-cost optoelectronic devices with superior performance. Since their perovskite crystal structure may have lattice instabilities against polar distortions, they are also being considered as potential photo-ferroelectrics. However, so far, research on their ferroelectricity has yielded inconclusive results and the subject is far from being settled. Here, we investigate, using a combined experimental and theoretical approach, the possible presence of electric polarization in tetragonal and orthorhombic CH3NH3PbBr3 (T-MAPB and O-MAPB). We found that T-MAPB does not sustain spontaneous polarization but, under an external electric field, it is projected into a metastable, ionic space-charge electret state. The electret can be frozen on cooling, producing a large and long-lasting polarization in O-MAPB. Molecular dynamics simulations show that the ferroelastic domain boundaries are able to trap charges and segregate ionic point defects, thus playing a favorable role in the stabilization of the electret. At lower temperatures, the lack of ferroelectric behavior is explained using first principles calculations as the result of the tight competition among many metastable states with randomly oriented polarization; this large configurational entropy does not allow a single polar state to dominate at any significant temperature range.

Ultralow Thermal Conductivity of Two-Dimensional Metal Halide Perovskites.

We report on the thermal conductivities of two-dimensional metal halide perovskite films measured by time domain thermoreflectance. Depending on the molecular substructure of ammonium cations and owing to the weaker interactions in the layered structures, the thermal conductivities of our two-dimensional hybrid perovskites range from 0.10 to 0.19 W m-1 K-1, which is drastically lower than that of their three-dimensional counterparts. We use molecular dynamics simulations to show that the organic component induces a reduction of the stiffness and sound velocities along with giving rise to vibrational modes in the 5-15 THz range that are absent in the three-dimensional counterparts. By systematically studying eight different two-dimensional hybrid perovskites, we show that the thermal conductivities of our hybrid films do not depend on the thicknesses of the organic layers and instead are highly dependent on the relative orientation of the organic chains sandwiched between the inorganic constituents.

Crystal-Size-Induced Band Gap Tuning in Perovskite Films.

A comprehensive picture explaining the effect of the crystal size in metal halide perovskite films on their opto-electronic characteristics is currently lacking. We report that perovskite nanocrystallites exhibit a wider band gap due to concurrent quantum confinement and size dependent structural effects, with the latter being remarkably distinct and attributed to the perturbation from the surface of the nanocrystallites affecting the structure of their core. This phenomenon might assist in the photo-induced charge separation within the perovskite in devices employing mesoporous layers as they restrict the size of nanocrystallites present in them. We demonstrate that the crystal size effect is widely applicable as it is ubiquitous in different compositions and deposition methods employed in the fabrication of state-of-the-art perovskite solar cells. This effect is a convenient and effective way to tune the band gap of perovskites.

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.

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

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.

Confinement of vibrational modes within crystalline lattices using thin amorphous layers.

It is possible to confine vibrational modes to a crystal by encapsulating it within thin disordered layers with the same average properties as the crystal. This is not due to an impedance mismatch between materials but, rather, to higher order moments in the distribution of density and stiffness in the disordered phase-i.e. it is a result of material substructure. The concept is elucidated in an idealized one-dimensional setting and then demonstrated for a realistic nanocrystalline geometry. This offers the prospect of specifically engineering higher order property distributions as an alternate means of managing phonons.

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.

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.

Origin of unusual bandgap shift and dual emission in organic-inorganic lead halide perovskites.

Emission characteristics of metal halide perovskites play a key role in the current widespread investigations into their potential uses in optoelectronics and photonics. However, a fundamental understanding of the molecular origin of the unusual blueshift of the bandgap and dual emission in perovskites is still lacking. In this direction, we investigated the extraordinary photoluminescence behavior of three representatives of this important class of photonic materials, that is, CH3NH3PbI3, CH3NH3PbBr3, and CH(NH2)2PbBr3, which emerged from our thorough studies of the effects of temperature on their bandgap and emission decay dynamics using time-integrated and time-resolved photoluminescence spectroscopy. The low-temperature (<100 K) photoluminescence of CH3NH3PbI3 and CH3NH3PbBr3 reveals two distinct emission peaks, whereas that of CH(NH2)2PbBr3 shows a single emission peak. Furthermore, irrespective of perovskite composition, the bandgap exhibits an unusual blueshift by raising the temperature from 15 to 300 K. Density functional theory and classical molecular dynamics simulations allow for assigning the additional photoluminescence peak to the presence of molecularly disordered orthorhombic domains and also rationalize that the unusual blueshift of the bandgap with increasing temperature is due to the stabilization of the valence band maximum. Our findings provide new insights into the salient emission properties of perovskite materials, which define their performance in solar cells and light-emitting devices.

Entropy-Suppressed Ferroelectricity in Hybrid Lead-Iodide Perovskites.

The actual nature of the electric polarization in hybrid lead-iodide perovskites is unveiled on the basis of ab initio and model results. A finite, albeit small electric polarization of few µC/cm(2) is found to be pervasive in this system, due to the polar-uncompensated alignment of methylammonium dimers, at least for temperature lower than the activation energy of dimer rotations; however, the presence of a large number of structural local minima corresponding to differently oriented polarization directions counteracts the stabilization of an ordered ferroelectric phase at the macroscale. According to our estimate, only for temperatures lower than 40-50 K a clear ferroelectric behavior is displayed. At higher temperature the polarization is progressively suppressed and the ferroelectric ordering hindered by the large configurational entropy, giving rise to a superparaelectric-like behavior at the macroscale.

Competing Forces in the Self-Assembly of Coupled ZnO Nanopyramids.

Self-assembly (SA) of nanostructures has recently gained increasing interest. A clear understanding of the process is not straightforward since SA of nanoparticles is a complex multiscale phenomenon including different driving forces. Here, we study the SA between aluminum doped ZnO nanopyramids into couples by combining inorganic chemistry and advanced electron microscopy techniques with atomistic simulations. Our results show that the SA of the coupled nanopyramids is controlled first by morphology, as coupling only occurs in the case of pyramids with well-developed facets of the basal planes. The combination of electron microscopy and atomistic modeling reveals that the coupling is further driven by strong ligand-ligand interaction between the bases of the pyramids as dominant force, while screening effects due to Al doping or solvent as well as core-core interaction are only minor contributions. Our combined approach provides a deeper understanding of the complex interplay between the interactions at work in the coupled SA of ZnO nanopyramids.

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.

Self-assembling of zinc phthalocyanines on ZnO (1010) surface through multiple time scales.

We adopt a hierarchic combination of theoretical methods to study the assembling of zinc phthalocyanines (ZnPcs) on a ZnO (1010) surface through multiple time scales. Atomistic simulations, such as model potential molecular dynamics and metadynamics, are used to study the energetics and short time evolution (up to ∼100 ns) of small ZnPc aggregates. The stability and the lifetime of large clusters is then studied by means of an atomistically informed coarse-grained model using classical molecular dynamics. Finally, the macroscopic time scale clustering phenomenon is studied by Metropolis Monte Carlo algorithms as a function of temperature and surface coverage. We provide evidence that at room temperature the aggregation is likely to occur at sufficiently high coverage, and we characterize the nature, morphology, and lifetime of ZnPc's clusters. We identify the molecular stripes oriented along [010] crystallographic directions as the most energetically stable aggregates.

Nonuniform growth of embedded silicon nanocrystals in an amorphous matrix.

The phase transformation of a metastable system occurs when islands of a second stable phase form and grow. The growth velocity of the islands controls the kinetics of the phase transformation. In this work we consider the amorphous-to-crystalline transformation in silicon as the prototype of a solid-to-solid transformation. The results of atomistic simulations are fit using an analytic model for the growth of [100]-oriented nanosized crystalline fibers embedded into an amorphous matrix. We demonstrate that the radius of the island does not grow, in general, at constant velocity. On the contrary, we identify a decelerated motion that is due to anisotropic effects of the crystal grain. Such a nonuniform growth should be taken into account in the modeling of solid-to-solid crystallization.