Low-loss contacts on textured substrates for inverted perovskite solar cells.

Inverted perovskite solar cells (PSCs) promise enhanced operating stability compared to their normal-structure counterparts1-3. To improve efficiency further, it is crucial to combine effective light management with low interfacial losses4,5. Here we develop a conformal self-assembled monolayer (SAM) as the hole-selective contact on light-managing textured substrates. Molecular dynamics simulations indicate that cluster formation during phosphonic acid adsorption leads to incomplete SAM coverage. We devise a co-adsorbent strategy that disassembles high-order clusters, thus homogenizing the distribution of phosphonic acid molecules, and thereby minimizing interfacial recombination and improving electronic structures. We report a laboratory-measured power conversion efficiency (PCE) of 25.3% and a certified quasi-steady-state PCE of 24.8% for inverted PSCs, with a photocurrent approaching 95% of the Shockley-Queisser maximum. An encapsulated device having a PCE of 24.6% at room temperature retains 95% of its peak performance when stressed at 65 °C and 50% relative humidity following more than 1,000 h of maximum power point tracking under 1 sun illumination. This represents one of the most stable PSCs subjected to accelerated ageing: achieved with a PCE surpassing 24%. The engineering of phosphonic acid adsorption on textured substrates offers a promising avenue for efficient and stable PSCs. It is also anticipated to benefit other optoelectronic devices that require light management.

Buried-Interface Engineering Enables Efficient and 1960-Hour ISOS-L-2I Stable Inverted Perovskite Solar Cells.

High-performance perovskite solar cells (PSCs) typically require interfacial passivation, yet this is challenging for the buried interface, owing to the dissolution of passivation agents during the deposition of perovskites. Here, this limitation is overcome with in situ buried-interface passivation-achieved via directly adding a cyanoacrylic-acid-based molecular additive, namely BT-T, into the perovskite precursor solution. Classical and ab initio molecular dynamics simulations reveal that BT-T spontaneously may self-assemble at the buried interface during the formation of the perovskite layer on a nickel oxide hole-transporting layer. The preferential buried-interface passivation results in facilitated hole transfer and suppressed charge recombination. In addition, residual BT-T molecules in the perovskite layer enhance its stability and homogeneity. A power-conversion efficiency (PCE) of 23.48% for 1.0 cm2 inverted-structure PSCs is reported. The encapsulated PSC retains 95.4% of its initial PCE following 1960 h maximum-power-point tracking under continuous light illumination at 65 °C (i.e., ISOS-L-2I protocol). The demonstration of operating-stable PSCs under accelerated ageing conditions represents a step closer to the commercialization of this emerging technology.

From Chalcogen Bonding to S-π Interactions in Hybrid Perovskite Photovoltaics.

The stability of hybrid organic-inorganic halide perovskite semiconductors remains a significant obstacle to their application in photovoltaics. To this end, the use of low-dimensional (LD) perovskites, which incorporate hydrophobic organic moieties, provides an effective strategy to improve their stability, yet often at the expense of their performance. To address this limitation, supramolecular engineering of noncovalent interactions between organic and inorganic components has shown potential by relying on hydrogen bonding and conventional van der Waals interactions. Here, the capacity to access novel LD perovskite structures that uniquely assemble through unorthodox S-mediated interactions is explored by incorporating benzothiadiazole-based moieties. The formation of S-mediated LD structures is demonstrated, including one-dimensional (1D) and layered two-dimensional (2D) perovskite phases assembled via chalcogen bonding and S-π interactions, through a combination of techniques, such as single crystal and thin film X-ray diffraction, as well as solid-state NMR spectroscopy, complemented by molecular dynamics simulations, density functional theory calculations, and optoelectronic characterization, revealing superior conductivities of S-mediated LD perovskites. The resulting materials are applied in n-i-p and p-i-n perovskite solar cells, demonstrating enhancements in performance and operational stability that reveal a versatile supramolecular strategy in photovoltaics.

Temporal disconnectivity of the energy landscape in glassy systems.

An alternative graphical representation of the potential energy landscape (PEL) has been developed and applied to a binary Lennard-Jones glassy system, providing insight into the unique topology of the system's potential energy hypersurface. With the help of this representation one is able to monitor the different explored basins of the PEL, as well as how--and mainly when--subsets of basins communicate with each other via transitions in such a way that details of the prior temporal history have been erased, i.e., local equilibration between the basins in each subset has been achieved. In this way, apart from detailed information about the structure of the PEL, the system's temporal evolution on the PEL is described. In order to gather all necessary information about the identities of two or more basins that are connected with each other, we consider two different approaches. The first one is based on consideration of the time needed for two basins to mutually equilibrate their populations according to the transition rate between them, in the absence of any effect induced by the rest of the landscape. The second approach is based on an analytical solution of the master equation that explicitly takes into account the entire explored landscape. It is shown that both approaches lead to the same result concerning the topology of the PEL and dynamical evolution on it. Moreover, a "temporal disconnectivity graph" is introduced to represent a lumped system stemming from the initial one. The lumped system is obtained via a specially designed algorithm [N. Lempesis, D. G. Tsalikis, G. C. Boulougouris, and D. N. Theodorou, J. Chem. Phys. 135, 204507 (2011)]. The temporal disconnectivity graph provides useful information about both the lumped and the initial systems, including the definition of "metabasins" as collections of basins that communicate with each other via transitions that are fast relative to the observation time. Finally, the two examined approaches are compared to an "on the fly" molecular dynamics-based algorithm [D. G. Tsalikis, N. Lempesis, G. C. Boulougouris, and D. N. Theodorou, J. Chem. Theory Comput. 6, 1307 (2010)].

Lumping analysis for the prediction of long-time dynamics: from monomolecular reaction systems to inherent structure dynamics of glassy materials.

In this work we develop, test, and implement a methodology that is able to perform, in an automated manner, "lumping" of a high-dimensional, discrete dynamical system onto a lower-dimensional space. Our aim is to develop an algorithm which, without any assumption about the nature of the system's slow dynamics, is able to reproduce accurately the long-time dynamics with minimal loss of information. Both the original and the lumped systems conform to master equations, related via the "lumping" analysis introduced by Wei and Kuo [Ind. Eng. Chem. Fundam. 8, 114 (1969)], and have the same limiting equilibrium probability distribution. The proposed method can be used in a variety of processes that can be modeled via a first order kinetic reaction scheme. Lumping affords great savings in the computational cost and reveals the characteristic times governing the slow dynamics of the system. Our goal is to approach the best lumping scheme with respect to three criteria, in order for the lumped system to be able to fully describe the long-time dynamics of the original system. The criteria used are: (a) the lumping error arising from the reduction process; (b) a measure of the magnitude of singular values associated with long-time evolution of the lumped system; and (c) the size of the lumped system. The search for the optimum lumping proceeds via Monte Carlo simulation based on the Wang-Landau scheme, which enables us to overcome entrapment in local minima in the above criteria and therefore increases the probability of encountering the global optimum. The developed algorithm is implemented to reproduce the long-time dynamics of a glassy binary Lennard-Jones mixture based on the idea of "inherent structures," where the rate constants for transitions between inherent structures have been evaluated via hazard plot analysis of a properly designed ensemble of molecular dynamics trajectories.

On the role of inherent structures in glass-forming materials: II. Reconstruction of the mean square displacement by rigorous lifting of the inherent structure dynamics.

In a previous paper, we investigated the role of inherent structures in the vitrification process of glass-forming materials, showing that the dynamical transitions between inherent structures (states) can be well predicted by a first-order kinetic scheme based on infrequent-event theory at low temperatures. In this work, we utilize and extend that methodology in order to completely reconstruct the system dynamics in the form of the mean square atomic displacement as a function of time at finite temperatures on the basis of the succession of transitions in a network of states, the vibrational contribution being evaluated on the basis of short molecular dynamics runs artificially trapped within each one of the states. In order to do so, we provide the mathematical formulation for lifting the coarse-grained Poisson process model of transitions between states back to the atomistic level and thereby reproducing the full dynamics of the atomistic system within the Poisson approximation. Our result shows excellent agreement for temperatures around and below the glass transition temperature of our model Lennard-Jones two-component mixtures. Clearly, our approach is able to reproduce the full dynamics of the atomistic system at low temperatures, where the Poisson approximation is valid.

On the role of inherent structures in glass-forming materials: I. The vitrification process.

In this work, we investigate the role of inherent structures in the vitrification process of glass-forming materials by using a two-component Lennard-Jones mixture. We start with a simplified model that describes the dynamics of the atomistic system as a Poisson process consisting of a series of transitions from one potential energy minimum (inherent structure) to another, the rate of individual transitions being described by a first-order kinetic law. We investigate the validity of this model by comparing the mean square displacement resulting from atomistic molecular dynamics (MD) trajectories with the corresponding mean square displacement based on inherent structures. Furthermore, in the case of vitrification via stepwise cooling, we identify the role of the potential energy landscape in determining the properties of the resulting glass. Interestingly, the cooling rate is not sufficient to define the resulting glass in a stepwise cooling process, because the time spent by the system at different temperatures (length of the steps) has a highly nonlinear impact on the properties of the resulting glass. In contrast to previous investigations of supercooled liquids, we focus on a range of temperatures close to and below the glass transition temperature, where the use of MD is incapable of producing equilibrated samples of the metastable supercooled state. Our aim is to develop a methodology that enables mapping the dynamics under these conditions to a coarse-grained first-order kinetic model based on the Poisson process approximation. This model can be used in order to extend our dynamical sampling ability to much broader time scales and therefore allow us to create computer glasses with cooling rates closer to those used experimentally. In a continuation to this work, we provide the mathematical formulation for lifting the coarse-grained Poisson process model and reproducing the full dynamics of the atomistic system.

Temperature accelerated dynamics in glass-forming materials.

In this work we propose a methodology for improving dynamical sampling in molecular simulations via temperature acceleration. The proposed approach combines the novel methods of Voter for temperature-accelerated dynamics with the multiple histogram reweighting method of Ferrenberg and Swendsen, its dynamical extension by Nieto-Draghi et al., and with hazard plot analysis, allowing optimal sampling with small computational cost over time scales inaccessible to classical molecular dynamics simulations by utilizing the "inherent structure" idea. The time evolution of the system is viewed as a succession of transitions between "basins" in its potential energy landscape, each basin containing a local minimum of the energy (inherent structure). Applying the proposed algorithm to a glass-forming material consisting of a mixture of spherical atoms interacting via Lennard-Jones potentials, we show that it is possible to perform an exhaustive search and evaluate rate constants for basin-to-basin transitions that cover several orders of magnitude on the time scale, far beyond the abilities of any competitive dynamical study, revealing an extreme ruggedness of the potential energy landscape in the vicinity of the glass transition temperature. By analyzing the network of inherent structures, we find that there are characteristic distances and rate constants related to the dynamical entrapment of the system in a neighborhood of basins (a metabasin), whereas evidence to support a random energy model is provided. The multidimensional configurational space displays a self-similar character depicted by a fractal dimension around 2.7 (+/-0.5) for the set of sampled inherent structures. Only transitions with small Euclidean measure can be considered as localized.