Tautomeric mixture coordination enables efficient lead-free perovskite LEDs.

Lead halide perovskite light-emitting diodes (PeLEDs) have demonstrated remarkable optoelectronic performance1-3. However, there are potential toxicity issues with lead4,5 and removing lead from the best-performing PeLEDs-without compromising their high external quantum efficiencies-remains a challenge. Here we report a tautomeric-mixture-coordination-induced electron localization strategy to stabilize the lead-free tin perovskite TEA2SnI4 (TEAI is 2-thiopheneethylammonium iodide) by incorporating cyanuric acid. We demonstrate that a crucial function of the coordination is to amplify the electronic effects, even for those Sn atoms that aren't strongly bonded with cyanuric acid owing to the formation of hydrogen-bonded tautomeric dimer and trimer superstructures on the perovskite surface. This electron localization weakens adverse effects from Anderson localization and improves ordering in the crystal structure of TEA2SnI4. These factors result in a two-orders-of-magnitude reduction in the non-radiative recombination capture coefficient and an approximately twofold enhancement in the exciton binding energy. Our lead-free PeLED has an external quantum efficiency of up to 20.29%, representing a performance comparable to that of state-of-the-art lead-containing PeLEDs6-12. We anticipate that these findings will provide insights into the stabilization of Sn(II) perovskites and further the development of lead-free perovskite applications.

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.

Pseudo-halide anion engineering for α-FAPbI3 perovskite solar cells.

Metal halide perovskites of the general formula ABX3-where A is a monovalent cation such as caesium, methylammonium or formamidinium; B is divalent lead, tin or germanium; and X is a halide anion-have shown great potential as light harvesters for thin-film photovoltaics1-5. Among a large number of compositions investigated, the cubic α-phase of formamidinium lead triiodide (FAPbI3) has emerged as the most promising semiconductor for highly efficient and stable perovskite solar cells6-9, and maximizing the performance of this material in such devices is of vital importance for the perovskite research community. Here we introduce an anion engineering concept that uses the pseudo-halide anion formate (HCOO-) to suppress anion-vacancy defects that are present at grain boundaries and at the surface of the perovskite films and to augment the crystallinity of the films. The resulting solar cell devices attain a power conversion efficiency of 25.6 per cent (certified 25.2 per cent), have long-term operational stability (450 hours) and show intense electroluminescence with external quantum efficiencies of more than 10 per cent. Our findings provide a direct route to eliminate the most abundant and deleterious lattice defects present in metal halide perovskites, providing a facile access to solution-processable films with improved optoelectronic performance.

MiMiC: A high-performance framework for multiscale molecular dynamics simulations.

MiMiC is a framework for performing multiscale simulations in which loosely coupled external programs describe individual subsystems at different resolutions and levels of theory. To make it highly efficient and flexible, we adopt an interoperable approach based on a multiple-program multiple-data (MPMD) paradigm, serving as an intermediary responsible for fast data exchange and interactions between the subsystems. The main goal of MiMiC is to avoid interfering with the underlying parallelization of the external programs, including the operability on hybrid architectures (e.g., CPU/GPU), and keep their setup and execution as close as possible to the original. At the moment, MiMiC offers an efficient implementation of electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) that has demonstrated unprecedented parallel scaling in simulations of large biomolecules using CPMD and GROMACS as QM and MM engines, respectively. However, as it is designed for high flexibility with general multiscale models in mind, it can be straightforwardly extended beyond QM/MM. In this article, we illustrate the software design and the features of the framework, which make it a compelling choice for multiscale simulations in the upcoming era of exascale high-performance computing.

Oxidative Defect Detection Within Free and Packed DNA Systems: A Quantum Mechanical/Molecular Mechanics (QM/ MM) Approach.

Base excision repair enzymes (BERs) detect and repair oxidative DNA damage with efficacy despite the small size of the defects and their often only minor structural impact. A charge transfer (CT) model for rapid scanning of DNA stretches has been evoked to explain the high detection rate in the face of numerous, small lesions. The viability of CT DNA defect detection is explored via hybrid QM/MM computational studies that leverage the accuracy of quantum mechanics (QM) for a region of interest and the descriptive power of molecularmechanics (MM) for the remainder of the system. We find that the presence of an oxidative lesion lowers theredox free energy of oxidation by approximately 1.0 eV regardless of DNA compaction (free DNA versus packed DNA in nucleosome core particles) and damage location indicating the high feasibility of a CT-based process for defect detection in DNA.

Recent Advances in First-Principles Based Molecular Dynamics.

First-principles molecular dynamics (FPMD) and its quantum mechanical-molecular mechanical (QM/MM) extensions are powerful tools to follow the real-time dynamics of a broad variety of systems in their ground as well as electronically excited states. The continued advances in computational power have enabled simulations of QM regions of larger sizes for more extended time scales. In addition, development of the parallel algorithms has boosted the performance of QM/MM methods even on existing computer architectures. In the case of density functional-based FPMD, systems of several hundreds to thousands of atoms can now be customarily simulated for tens to hundreds of picoseconds. In spite of this progress, the time scale limitations remain severe, especially when high-rung exchange-correlation functionals or high-level wave function based quantum mechanical methods are used. To ameliorate this, a large number of enhanced sampling methods have been introduced but most of the approaches that have been developed to increase the efficiency of FPMD based simulations sacrifice the real-time dynamics in favor of enhancing sampling. Here, we present some recent advances in boosting the efficiency of FPMD based simulations while keeping the full dynamic information. These include a highly efficient recent implementation of FPMD-based QM/MM simulations that not only enables fully flexible combinations of different electronic structure methods and force fields via a highly efficient communication library, it also fully exploits parallelism for both quantum and classical descriptions. The second type of acceleration methods we discuss is a large family of specially devised multiple-time-step algorithms that make use of suitable breakups of the total nuclear forces into fast components that can be calculated via lower level methods and slowly varying correction forces evaluated with a high-level method at long time intervals. The computational gain of this scheme mostly depends on the cost difference between the two methods and advantageous combinations can yield large speedups without compromising the accuracy of the high-level method. And finally, the third class of FPMD acceleration methods presented here are machine learning models to accelerated FPMD and their powerful combinations with multiple-time-step techniques. The combination of all the approaches enables substantial speedups of FPMD simulations of several orders of magnitude while fully preserving the real-time dynamics and accuracy.

MiMiCPy: An Efficient Toolkit for MiMiC-Based QM/MM Simulations.

MiMiC is a highly flexible, extremely scalable multiscale modeling framework. It couples the CPMD (quantum mechanics, QM) and GROMACS (molecular mechanics, MM) codes. The code requires preparing separate input files for the two programs with a selection of the QM region. This can be a tedious procedure prone to human error, especially when dealing with large QM regions. Here, we present MiMiCPy, a user-friendly tool that automatizes the preparation of MiMiC input files. It is written in Python 3 with an object-oriented approach. The main subcommand PrepQM can be used to generate MiMiC inputs directly from the command line or through a PyMOL/VMD plugin for visually selecting the QM region. Many other subcommands are also provided for debugging and fixing MiMiC input files. MiMiCPy is designed with a modular structure that allows seamless extensions to new program formats depending on the requirements of MiMiC.

A theoretical perspective of the ultrafast transient absorption dynamics of CsPbBr3.

Transient absorption spectra (TAS) of lead halide perovskites can provide important insights into the nature of the photoexcited state dynamics of this prototypical class of materials. Here, we perform ground and excited state molecular dynamics (MD) simulations within a restricted open shell Kohn-Sham (ROKS) approach in order to interpret the characteristic features of the TAS of CsPbBr3 . Our results reveal that properties such as the finite temperature band gap, the Stokes shift, and therefore, also the TAS are strongly size-dependent. Our TAS simulations show an early positive red-shifted feature on the fs scale that can be explained by geometric relaxation in the excited state. As excited-state processes can crucially affect the electronic properties of this class of photoactive materials, our observations are an important ingredient for further optimization of lead halide based optoelectronic devices.

Hydration dynamics and IR spectroscopy of 4-fluorophenol.

Halogenated groups are relevant in pharmaceutical applications and potentially useful spectroscopic probes for infrared spectroscopy. In this work, the structural dynamics and infrared spectroscopy of para-fluorophenol (F-PhOH) and phenol (PhOH) is investigated in the gas phase and in water using a combination of experiment and molecular dynamics (MD) simulations. The gas phase and solvent dynamics around F-PhOH and PhOH is characterized from atomistic simulations using empirical energy functions with point charges or multipoles for the electrostatics, Machine Learning (ML) based parametrizations and with full ab initio (QM) and mixed Quantum Mechanical/Molecular Mechanics (QM/MM) simulations with a particular focus on the CF- and OH-stretch region. The CF-stretch band is heavily mixed with other modes whereas the OH-stretch in solution displays a characteristic high-frequency peak around 3600 cm-1 most likely associated with the -OH group of PhOH and F-PhOH together with a characteristic progression below 3000 cm-1 due to coupling with water modes which is also reproduced by several of the simulations. Solvent and radial distribution functions indicate that the CF-site is largely hydrophobic except for simulations using point charges which renders them unsuited for correctly describing hydration and dynamics around fluorinated sites. The hydrophobic character of the CF-group is particularly relevant for applications in pharmaceutical chemistry with a focus on local hydration and interaction with the surrounding protein.

A multiple time step algorithm for trajectory surface hopping simulations.

A multiple time step (MTS) algorithm for trajectory surface hopping molecular dynamics has been developed, implemented, and tested. The MTS scheme is an extension of the ab initio implementation for Born-Oppenheimer molecular dynamics presented in the work of Liberatore et al. [J. Chem. Theory Comput. 14, 2834 (2018)]. In particular, the MTS algorithm has been modified to enable the simulation of non-adiabatic processes with the trajectory surface hopping (TSH) method and Tully's fewest switches algorithm. The specificities of the implementation lie in the combination of Landau-Zener and Tully's transition probabilities during the inner MTS time steps. The new MTS-TSH method is applied successfully to the photorelaxation of protonated formaldimine, showing that the important characteristics of the process are recovered by the MTS algorithm. A computational speed-up between 1.5 and 3 has been obtained compared to standard TSH simulations, which is close to the ideal values that could be obtained with the computational settings considered.

Nanoscale Phase Segregation in Supramolecular π-Templating for Hybrid Perovskite Photovoltaics from NMR Crystallography.

The use of layered perovskites is an important strategy to improve the stability of hybrid perovskite materials and their optoelectronic devices. However, tailoring their properties requires accurate structure determination at the atomic scale, which is a challenge for conventional diffraction-based techniques. We demonstrate the use of nuclear magnetic resonance (NMR) crystallography in determining the structure of layered hybrid perovskites for a mixed-spacer model composed of 2-phenylethylammonium (PEA+) and 2-(perfluorophenyl)ethylammonium (FEA+) moieties, revealing nanoscale phase segregation. Moreover, we illustrate the application of this structure in perovskite solar cells with power conversion efficiencies that exceed 21%, accompanied by enhanced operational stability.

All-atom simulations disentangle the functional dynamics underlying gene maturation in the intron lariat spliceosome.

The spliceosome (SPL) is a majestic macromolecular machinery composed of five small nuclear RNAs and hundreds of proteins. SPL removes noncoding introns from precursor messenger RNAs (pre-mRNAs) and ligates coding exons, giving rise to functional mRNAs. Building on the first SPL structure solved at near-atomic-level resolution, here we elucidate the functional dynamics of the intron lariat spliceosome (ILS) complex through multi-microsecond-long molecular-dynamics simulations of ∼1,000,000 atoms models. The ILS essential dynamics unveils (i) the leading role of the Spp42 protein, which heads the gene maturation by tuning the motions of distinct SPL components, and (ii) the critical participation of the Cwf19 protein in displacing the intron lariat/U2 branch helix. These findings provide unprecedented details on the SPL functional dynamics, thus contributing to move a step forward toward a thorough understanding of eukaryotic pre-mRNA splicing.

Molecular Basis of CLC Antiporter Inhibition by Fluoride.

CLC channels and transporters conduct or transport various kinds of anions, with the exception of fluoride, which acts as an effective inhibitor. Here, we performed sub-nanosecond DFT-based QM/MM simulations of the E. coli anion/proton exchanger ClC-ec1 and observed that fluoride binds incoming protons within the selectivity filter, with excess protons shared with the gating glutamate E148. Depending on E148 conformation, the competition for the proton can involve either a direct F-/E148 interaction or the modulation of water molecules bridging the two anions. The direct interaction locks E148 in a conformation that does not allow for proton transport, and thus inhibits protein function.

Crown Ether Modulation Enables over 23% Efficient Formamidinium-Based Perovskite Solar Cells.

The use of molecular modulators to reduce the defect density at the surface and grain boundaries of perovskite materials has been demonstrated to be an effective approach to enhance the photovoltaic performance and device stability of perovskite solar cells. Herein, we employ crown ethers to modulate perovskite films, affording passivation of undercoordinated surface defects. This interaction has been elucidated by solid-state nuclear magnetic resonance and density functional theory calculations. The crown ether hosts induce the formation of host-guest complexes on the surface of the perovskite films, which reduces the concentration of surface electronic defects and suppresses nonradiative recombination by 40%, while minimizing moisture permeation. As a result, we achieved substantially improved photovoltaic performance with power conversion efficiencies exceeding 23%, accompanied by enhanced stability under ambient and operational conditions. This work opens a new avenue to improve the performance and stability of perovskite-based optoelectronic devices through supramolecular chemistry.

Multidisciplinary Preclinical Investigations on Three Oxamniquine Analogues as New Drug Candidates for Schistosomiasis*.

Schistosomiasis is a disease of poverty affecting millions of people. Praziquantel (PZQ), with its strengths and weaknesses, is the only treatment available. We previously reported findings on three lead compounds derived from oxamniquine (OXA), an old antischistosomal drug: ferrocene-containing (Fc-CH2 -OXA), ruthenocene-containing (Rc-CH2 -OXA) and benzene-containing (Ph-CH2 -OXA) OXA derivatives. These derivatives showed excellent in vitro activity against both Schistosoma mansoni larvae and adult worms and S. haematobium adult worms, and were also active in vivo against adult S. mansoni. Encouraged by these promising results, we conducted additional in-depth preclinical studies and report in this investigation on metabolic stability studies, in vivo studies on S. haematobium and juvenile S. mansoni, computational simulations, and formulation development. Molecular dynamics simulations supported the in vitro results on the target protein. Though all three compounds were poorly stable within an acidic environment, they were only slightly cleared in the in vitro liver model. This is likely the reason why the promising in vitro activity did not translate into in vivo activity on S. haematobium. This limitation could not be overcome by the formulation of lipid nanocapsules as a way to improve the in vivo activity. Further studies should focus on increasing the compound's bioavailability, to reach an active concentration in the microenvironment of the parasite.

Regulation of adenylyl cyclase 5 in striatal neurons confers the ability to detect coincident neuromodulatory signals.

Long-term potentiation and depression of synaptic activity in response to stimuli is a key factor in reinforcement learning. Strengthening of the corticostriatal synapses depends on the second messenger cAMP, whose synthesis is catalysed by the enzyme adenylyl cyclase 5 (AC5), which is itself regulated by the stimulatory Gαolf and inhibitory Gαi proteins. AC isoforms have been suggested to act as coincidence detectors, promoting cellular responses only when convergent regulatory signals occur close in time. However, the mechanism for this is currently unclear, and seems to lie in their diverse regulation patterns. Despite attempts to isolate the ternary complex, it is not known if Gαolf and Gαi can bind to AC5 simultaneously, nor what activity the complex would have. Using protein structure-based molecular dynamics simulations, we show that this complex is stable and inactive. These simulations, along with Brownian dynamics simulations to estimate protein association rates constants, constrain a kinetic model that shows that the presence of this ternary inactive complex is crucial for AC5's ability to detect coincident signals, producing a synergistic increase in cAMP. These results reveal some of the prerequisites for corticostriatal synaptic plasticity, and explain recent experimental data on cAMP concentrations following receptor activation. Moreover, they provide insights into the regulatory mechanisms that control signal processing by different AC isoforms.

Why choosing the right partner is important: stabilization of ternary CsyGUAxFA(1-y-x)PbI3 perovskites.

Lead halide perovskites with mixtures of monovalent cations have attracted wide attention due to the possibility of preferentially stabilizing the perovskite phase with respect to photovoltaically less suitable competing phases. Here, we present a theoretical analysis and interpretation of the phase stability of binary (CH6N3)x[HC(NH2)2](1-x)PbI3 = GUAxFA(1-x)PbI3 and ternary CsyGUAxFA(1-y-x)PbI3 mixtures. We first estimate if such mixtures are stable and if they lead to a stabilization of the perovskite phase based on static Density Functional Theory (DFT) calculations. In order to investigate the finite temperature stability of the phases, we also employ first-principles molecular dynamics (MD) simulations. It turns out that in contrast to the FA+-rich case of FA/Cs mixtures, although mixing of FA/GUA is possible, it is not sufficient to stabilize the perovskite phase at room temperature. In contrast, stable ternary mixtures that contain 17% of Cs+ can be formed that lead to a preferential stabilization of the perovskite phase. In such a way, the enthalpic destabilization due to the introduction of a too large/too small cation that lies outside the Goldschmidt tolerance range can be (partially) compensated through the introduction of a third cation with complementary size. This allows to suggest a new design principle for the preparation of stable perovskite structures at room temperature with cations that lie outside the Goldschmidt range through mixtures with size-complementary cations in such a way that the effective average cation radius of the mixture lies within the stability range.

Guanine-Stabilized Formamidinium Lead Iodide Perovskites.

Formamidinium (FA) lead iodide perovskite materials feature promising photovoltaic performances and superior thermal stabilities. However, conversion of the perovskite α-FAPbI3 phase to the thermodynamically stable yet photovoltaically inactive δ-FAPbI3 phase compromises the photovoltaic performance. A strategy is presented to address this challenge by using low-dimensional hybrid perovskite materials comprising guaninium (G) organic spacer layers that act as stabilizers of the three-dimensional α-FAPbI3 phase. The underlying mode of interaction at the atomic level is unraveled by means of solid-state nuclear magnetic resonance spectroscopy, X-ray crystallography, transmission electron microscopy, molecular dynamics simulations, and DFT calculations. Low-dimensional-phase-containing hybrid FAPbI3 perovskite solar cells are obtained with improved performance and enhanced long-term stability.

Association of Both Inhibitory and Stimulatory Gα Subunits Implies Adenylyl Cyclase 5 Deactivation.

Adenylyl cyclase (AC) generates cyclic AMP required for a variety of cellular functions, and its regulation plays a major role in cellular signal transduction in eukaryotes and prokaryotes. All membrane-bound AC isoforms in eukaryotes can be activated by stimulatory G-proteins, but only AC1, AC5, and AC6 can be both stimulated and inhibited by active Gα subunits, Gαs and Gαi, respectively. In principle, these Gαi-sensitive AC isoforms could form both binary and ternary complexes with Gα subunits due to the noncompetitive association of inhibitory and stimulatory Gα. However, the formation and possible catalytic activity of a putative ternary complex have not yet been experimentally confirmed due to its proposed short-lived nature. Here, the catalytic activity of such a ternary complex consisting of apo AC5, stimulatory Gαolf, and inhibitory Gαi1 is investigated via classical molecular dynamics simulations. Trajectories of inhibited and stimulated binary complexes, AC5:Gαi1 and AC5:Gαolf, respectively, as well as Gα-free AC5 were also obtained to compare the sampled AC5 conformation in the ternary complex to those sampled under different Gα conditions. This comparison suggests that association of both Gα subunits results in an AC5 conformation similar to that sampled by the AC5:Gαi1 complex, indicating that the ternary complex mainly samples an inactive conformation.