Alkyl Chains Tune Molecular Orientations to Enable Dual Passivation in Inverted Perovskite Solar Cells.

Nonradiative recombination losses occurring at the interface pose a significant obstacle to achieve high-efficiency perovskite solar cells (PSCs), particularly in inverted PSCs. Passivating surface defects using molecules with different functional groups represents one of the key strategies for enhancing PSCs efficiency. However, a lack of insight into the passivation orientation of molecules on the surface is a challenge for rational molecular design. In this study, aminothiol hydrochlorides with different alkyl chains but identical electron-donating (-SH) and electron-withdrawing (-NH3+) groups were employed to investigate the interplay between molecular structure, orientation, and interaction on perovskite surface. The 2-Aminoethane-1-thiol hydrochloride with shorter alkyl chains exhibited a preference of parallel orientations, which facilitating stronger interactions with the surface defects through strong coordination and hydrogen bonding. The resultant perovskite films following defect passivation demonstrate reduced ion migration, inhibition of nonradiative recombination, and more n-type characteristics for efficient electron transfer. Consequently, an impressive power conversion efficiency of 25% was achieved, maintaining 95% of its initial efficiency after 500 hours of continuous maximum power point tracking.

Understanding cellulose pyrolysis via ab initio deep learning potential field.

Comprehensive and dynamic studies of cellulose pyrolysis reaction mechanisms are crucial in designing experiments and processes with enhanced safety, efficiency, and sustainability. The details of the pyrolysis mechanism are not readily available from experiments but can be better described via molecular dynamics (MD) simulations. However, the large size of cellulose molecules challenges accurate ab initio MD simulations, while existing reactive force field parameters lack precision. In this work, precise ab initio deep learning potentials field (DPLF) are developed and applied in MD simulations to facilitate the study of cellulose pyrolysis mechanisms. The formation mechanism and production rate of both valuable and greenhouse products from cellulose at temperatures larger than 1073 K are comprehensively described. This study underscores the critical role of advanced simulation techniques, particularly DLPF, in achieving efficient and accurate understanding of cellulose pyrolysis mechanisms, thus promoting wider industrial applications.

A-D-A Type Nonfullerene Acceptors Synthesized by Core Segmentation and Isomerization for Realizing Organic Solar Cells with Low Nonradiative Energy Loss.

Reducing non-radiative recombination energy loss (ΔEnonrad ) in organic solar cells (OSCs) has been considered an effective method to improve device efficiency. In this study, the backbone of PTBTT-4F/4Cl is divided into D1-D2-D3 segments and reconstructed. The isomerized TPBTT-4F/4Cl obtains stronger intramolecular charge transfer (ICT), thus leading to elevated highest occupied molecular orbital (HOMO) energy level and reduced bandgap (Eg ). According to ELoss  = Eg- qVOC , the reduced Eg and enhanced open circuit voltage (VOC ) result in lower ELoss , indicating that ELoss has been effectively suppressed in the TPBTT-4F/4Cl based devices. Furthermore, compared to PTBTT derivatives, the isomeric TPBTT derivatives exhibit more planar molecular structure and closer intermolecular stacking, thus affording higher crystallinity of the neat films. Therefore, the reduced energy disorder and corresponding lower Urbach energy (Eu ) of the TPBTT-4F/4Cl blend films lead to low ELoss and high charge-carrier mobility of the devices. As a result, benefitting from synergetic control of molecular stacking and energetic offsets, a maximum power conversion efficiency (PCE) of 15.72% is realized from TPBTT-4F based devices, along with a reduced ΔEnonrad of 0.276 eV. This work demonstrates a rational method of suppressing VOC loss and improving the device performance through molecular design engineering by core segmentation and isomerization.

Epitaxial Growth of α-FAPbI3 at a Well-Matched Heterointerface for Efficient Perovskite Solar Cells and Solar Modules.

Although the FAPbI3 perovskite system exhibits an impressive optoelectronic characteristic and thermal stability because of its energetically unstable black phase at room temperature, it is considerably challenging to attain a controllable and oriented nucleation of α-FAPbI3 . To overcome this challenge, a 2D perovskite with a released inorganic octahedral distortion designed by weakening the hydrogen interactions between the organic interlayer and [PbI6 ]4- octahedron is presented in this study. A highly matched heterointerface can be formed between the (002) facet of the 2D structure and the (100) crystal plane of the cubic α-FAPbI3 , thereby lowering the crystallization energy and inducing a heterogeneous nucleation of α-FAPbI3 . This "epitaxial growth" mechanism results form the highly preferred crystallographic orientation of the (100) facets, improved crystal quality and film uniformity, substantially increased charge transporting characteristics, and suppressed nonradiative recombination losses. An impressive power conversion efficiency (PCE) of 25.4% (certified 25.2%) is achieved using target PSCs, which demonstrates outstanding ambient and operational stability. The feasibility of this strategy is proved for the scalable deposition of homogeneous and high-quality perovskite thin films by demonstrating the remarkably increased PCE of the large-area perovskite solar module, from 18.2% to 20.1%.

Nucleation Regulation and Mesoscopic Dielectric Screening in α-FAPbI3.

While significant advancements in power conversion efficiencies (PCEs) of α-FAPbI3perovskite solar cells (PSCs) have been made, attaining controllable perovskite crystallization is still a considerable hurdle. This challenge stems from the initial formation of δ-FAPbI3, a more energetically stable phase than the desired black α-phase, during film deposition. This disrupts the heterogeneous nucleation of α-FAPbI3, causing the formation of mixed phases and defects. To this end, polarity engineering using molecular additives, specifically ((methyl-sulfonyl)phenyl)ethylamines (MSPEs) are introduced. The findings reveal that the interaction of PbI2-MSPEs-FAI intermediates is enhanced with the increased polarity of MSPEs, which in turn expedites the nucleation of α-FAPbI3. This leads to the development of high-quality α-FAPbI3 films, characterized by vertical crystal orientation and reduced residual stresses. Additionally, the increased dipole moment of MSPE at perovskite grain boundaries attenuates Coulomb attractions among charged defects and screens carrier capture process, thereby diminishing non-radiative recombination. Utilizing these mechanisms, PSCs treated with highly polar 2-(4-MSPE) achieve an impressive PCE of 25.2% in small-area devices and 20.5% in large-area perovskite solar modules (PSMs) with an active area of 70 cm2. These results demonstrate the effectiveness of this strategy in achieving controllable crystallization of α-FAPbI3, paving the way for scalable-production of high-efficiency PSMs.

Molecular-Level Anion and Li+ Co-Regulation by Amphoteric Polymer Separator for High-Rate Stable Lithium Metal Anode.

Regulating ion transport is a prevailing strategy to suppress lithium dendrite growth, in which the distribution of ion regulatory sites plays an important role. Here a hyperbranched polyamidoamine (HBPA) grafted polyethylene (PE) composite separator (HBPA-g-PE) is reported. The densely and uniformly distributed positive -NH2 and negative -CHNO- groups efficiently restrict the anion migration and promote Li+ transport at the surface of the lithium metal anode. The obtained Li foil symmetric cell delivers a stable cycle performance with a low-voltage hysteresis of 130 mV for over 1500 h (3000 cycles) at an ultrahigh current density of 20 mA cm-2 and a practical areal capacity of 5 mAh cm-2. Moreover, HBPA-g-PE separator enables a practical lithium-sulfur battery to achieve over 200-cycle stable performance with initial and retained capacity of 700 and 455 mAh g-1, at a high sulfur loading of 4 mg cm-2 and a low electrolyte content/sulfur loading ratio of 8 µL mg-1.

Asymmetric Non-Fullerene Acceptor Derivatives Incorporated Ternary Organic Solar Cells.

Incorporating ITIC derivatives as guest acceptors into binary host systems is an effective strategy for constructing high-performance ternary organic solar cells (TOSCs). In this work, we introduced A-D-A type ITIC derivatives PTBTT-4F (asymmetric) and PTBTP-4F (symmetric) into the PM6:BTP-BO-4F (Y6-BO) binary blend and investigated the impacts of two guest acceptors on the performance of TOSCs. Differentiated device performance was observed, although PTBTT-4F and PTBTP-4F presented similar chemical structures and comparable absorptions. The PTBTT-4F ternary devices exhibited an improved power conversion efficiency (PCE) of 17.67% with increased open circuit (VOC) and current density (JSC), whereas the PTBTP-4F-based ternary devices yielded a relatively lower PCE of 16.34%. PTBTT-4F showed much better compatibility with the host acceptor BTP-BO-4F, so that they formed a well-mixed alloy phase state; more precise phase separation and increased crystallinity were thus induced in the ternary blends, leading to reduced molecular recombination and improved charge mobilities, which contributed to improved fill factors of the ternary devices. In addition, the optimized PTBTT-4F devices exhibited good performance tolerance of the photoactive layer thickness, as they even delivered a PCE of 15.25% when the active layer was as thick as up to ∼300 nm.

Harnessing Heteropolar Lithium Polysulfides by Amphoteric Polymer Binder for Facile Manufacturing of Practical Li-S Batteries.

Enabling efficient and durable charge storage under high sulfur loading and lean electrolyte remains a paramount challenge for Li-S battery technology to truly demonstrate its commercial viability. This work reports an amphoteric polymer binder, whose negatively and positively charged moieties allow for coregulation of both lithium cations and heteropolar lithium polysulfides through multiple intermolecular interactions. These interactions and the physical properties lead to simultaneously improved Li+ transport, polysulfide adsorption and catalysis, cathode robustness and anode stability. Therefore, this multifunctional binder endows Li-S batteries with compelling overall performances even under rigorous conditions. At low sulfur loading and copious electrolyte, the cell shows a low capacity-fading rate of 0.056% cycle-1 upon 700 cycles. At sulfur loading of 6.8 mg cm-2 and low E/S of 6 µL mg-1 , the cell still delivers stable areal capacities between 4.2 and 4.8 mAh cm-2 in 50 cycles without obvious decay at 0.2 C. The commercial feasibility of this work is further manifested by its zero added weight, low material cost, and ease of manufacturing and scale-up. The efficacy and simplicity of this work symbolize an example of lab-scale battery research aiming at improved technology and manufacturing readiness level.

Near-Ultraviolet Circular Dichroism and Two-Dimensional Spectroscopy of Polypeptides.

A fully quantitative theory of the relationship between protein conformation and optical spectroscopy would facilitate deeper insights into biophysical and simulation studies of protein dynamics and folding. In contrast to intense bands in the far-ultraviolet, near-UV bands are much weaker and have been challenging to compute theoretically. We report some advances in the accuracy of calculations in the near-UV, which were realised through the consideration of the vibrational structure of the electronic transitions of aromatic side chains.

Force Fields for Macromolecular Assemblies Containing Diketopyrrolopyrrole and Thiophene.

Utilizing a force-matching procedure, we parametrize new force fields systematically for large conjugated systems. We model both conjugated polymers and molecular crystals that contain diketopyrrolopyrrole, thiophene, and thieno[3,2-b]thiophene units. These systems have recently been found to have low band gaps, which exhibit high efficiency for photovoltaic devices. The equilibrium structures, forces, and energies of the building block chromophores, diketopyrrolopyrrole, thiophene, and thieno[3,2-b]thiophene computed using our parameters are comparable to those computed using the reference electronic structure method. We assess the suitability of this new force field for electronic property calculations by comparing the electronic excitation properties computed along classical and ab initio molecular dynamics trajectories. For both trajectories, we find similar distributions of TDDFT-calculated excitation energies and oscillator strengths for the building block chromophore diketopyrrolopyrrole-thieno[3,2-b]thiophene. The structural, dynamic, and electronic properties of the macromolecular assemblies built upon these chromophores are characterized. For both polymers and molecular crystals, pronounced peaks around 0° or 180° are observed for the torsions between chromophores under ambient conditions. The high planarity in these systems can promote local ordering and π-π stacking, thereby potentially facilitating charge transport across these materials. For the model conducting polymers, we found that the fluctuations in the density of states per chain per monomer is negligibly small and does not vary significantly with chains comprising 20-40 monomers. Analysis of the electron-hole distributions and the transition density matrices indicates that the delocalized length is approximately 4-6 monomers, which is in good agreement with other theoretical and experimental studies of different conducting polymers. For the molecular crystals, our investigation of the characteristic time scale of the fluctuation in the excitonic couplings shows that a low-frequency vibration below 100 cm-1 is observed for the nearest neighbors. These observations are in line with previous studies on other molecular crystals, in which low-frequency vibrations are believed to be responsible for the large modulation of the excitonic coupling. Thus, our approach and the new force fields provide a direct route for studying the structure-property relations and the molecular level origins of the high efficiency of these classes of materials.

Cluster integrals and virial coefficients for realistic molecular models.

We present a concise, general, and efficient procedure for calculating the cluster integrals that relate thermodynamic virial coefficients to molecular interactions. The approach encompasses nonpairwise intermolecular potentials generated from quantum chemistry or other sources; a simple extension permits efficient evaluation of temperature and other derivatives of the virial coefficients. We demonstrate with a polarizable model of water. We argue that cluster-integral methods are a potent yet underutilized instrument for the development and application of first-principles molecular models and methods.

A comparative study of mechanisms of the adsorption of CO2 confined within graphene-MoS2 nanosheets: a DFT trend study.

The space within the interlayer of 2-dimensional (2D) nanosheets provides new and intriguing confined environments for molecular interactions. However, atomic level understanding of the adsorption mechanism of CO2 confined within the interlayer of 2D nanosheets is still limited. Herein, we present a comparative study of the adsorption mechanisms of CO2 confined within graphene-molybdenum disulfide (MoS2) nanosheets using density functional theory (DFT). A comprehensive analysis of CO2 adsorption energies (E AE) at various interlayer spacings of different multilayer structures comprising graphene/graphene (GrapheneB) and MoS2/MoS2 (MoS2B) bilayers as well as graphene/MoS2 (GMoS2) and MoS2/graphene (MoS2G) hybrids is performed to obtain the most stable adsorption configurations. It was found that 7.5 Å and 8.5 Å interlayer spacings are the most stable conformations for CO2 adsorption on the bilayer and hybrid structures, respectively. Adsorption energies of the multilayer structures decreased in the following trend: MoS2B > GrapheneB > MoS2G > GMoS2. By incorporating van der Waals (vdW) interactions between the CO2 molecule and the surfaces, we find that CO2 binds more strongly on these multilayer structures. Furthermore, there is a slight discrepancy in the binding energies of CO2 adsorption on the heterostructures (GMoS2, MoS2G) due to the modality of the atom arrangement (C-Mo-S-O and Mo-S-O-C) in both structures, indicating that conformational anisotropy determines to a certain degree its CO2 adsorption energy. Meanwhile, Bader charge analysis shows that the interaction between CO2 and these surfaces causes charge transfer and redistributions. By contrast, the density of states (DOS) plots show that CO2 physisorption does not have a substantial effect on the electronic properties of graphene and MoS2. In summary, the results obtained in this study could serve as useful guidance in the preparation of graphene-MoS2 nanosheets for the improved adsorption efficiency of CO2.

Reverse energy partitioning-An efficient algorithm for computing the density of states, partition functions, and free energy of solids.

A robust and model free Monte Carlo simulation method is proposed to address the challenge in computing the classical density of states and partition function of solids. Starting from the minimum configurational energy, the algorithm partitions the entire energy range in the increasing energy direction ("upward") into subdivisions whose integrated density of states is known. When combined with the density of states computed from the "downward" energy partitioning approach [H. Do, J. D. Hirst, and R. J. Wheatley, J. Chem. Phys. 135, 174105 (2011)], the equilibrium thermodynamic properties can be evaluated at any temperature and in any phase. The method is illustrated in the context of the Lennard-Jones system and can readily be extended to other molecular systems and clusters for which the structures are known.

Calculation of high-order virial coefficients for the square-well potential.

Accurate virial coefficients B_{N}(λ,ɛ) (where ɛ is the well depth) for the three-dimensional square-well and square-step potentials are calculated for orders N=5-9 and well widths λ=1.1-2.0 using a very fast recursive method. The efficiency of the algorithm is enhanced significantly by exploiting permutation symmetry and by storing integrands for reuse during the calculation. For N=9 the storage requirements become sufficiently large that a parallel algorithm is developed. The methodology is general and is applicable to other discrete potentials. The computed coefficients are precise even near the critical temperature, and thus open up possibilities for analysis of criticality of the system, which is currently not accessible by any other means.

Developing accurate molecular mechanics force fields for conjugated molecular systems.

A rapid method to parameterize the intramolecular component of classical force fields for complex conjugated molecules is proposed. The method is based on a procedure of force matching with a reference electronic structure calculation. It is particularly suitable for those applications where molecular dynamics simulations are used to generate structures that are therefore analysed by electronic structure methods, because it is possible to build force fields that are consistent with electronic structure calculations that follow classical simulations. Such applications are commonly encountered in organic electronics, spectroscopy of complex systems and photobiology (e.g. photosynthetic systems). We illustrate the method by parameterizing the force fields of a molecule used in molecular semiconductors (2,2-dicyanovinyl-capped S,N-heteropentacene or DCV-SN5), a polymeric semiconductor (thieno[3,2-b]thiophene-diketopyrrolopyrrole TT-DPP) and a chromophore embedded in a protein environment (15,16-dihydrobiliverdin or DBV) where several hundreds of parameters need to be optimized in parallel.

Calculation of the vibrational frequencies of carbon clusters and fullerenes with empirical potentials.

Vibrational frequencies for carbon clusters, fullerenes and nanotubes evaluated using empirical carbon-carbon potentials are presented. For linear and cyclic clusters, frequencies evaluated with the reactive empirical bond order (REBO) potential provide the closest agreement with experiment. The mean absolute deviation (MAD) between experiment and the calculated harmonic frequencies is 79 cm(-1) for the bending modes and 76 cm(-1) for the stretching modes. The effects of anharmonicity are included via second order vibrational perturbation theory and tend to increase the frequency of the bending modes while the stretching modes have negative shifts in the region of 20-60 cm(-1), with larger shifts for the higher frequency modes. This results in MADs for the bending and stretching modes of 84 cm(-1) and 58 cm(-1), respectively. For the fullerene molecule C60, the high frequency modes are predicted to have harmonic frequencies that are significantly higher than experiment, and this is not corrected by accounting for anharmonicity. This overestimation of experimental observed frequencies is also evident in the calculated frequencies of the G band in nanotubes. This suggests that the REBO potential is not optimal for these larger systems and it is shown that adjustment of the parameters within the potential leads to closer agreement with experiment, particularly if higher and lower frequency modes are considered separately.

Proton transfer or hemibonding? The structure and stability of radical cation clusters.

The basin hopping search algorithm in conjunction with second-order Møller-Plesset perturbation theory is used to determine the lowest energy structures of the radical cation clusters (NH3)n˙(+), (H2O)n˙(+), (HF)n˙(+), (PH3)n˙(+), (H2S)n˙(+) and (HCl)n˙(+), where n = 2-4. The energies of the most stable structures are subsequently evaluated using coupled cluster theory in conjunction with the aug-cc-pVTZ basis set. These cationic clusters can adopt two distinct structural types, with some clusters showing an unusual type of bonding, often referred to as hemibonding, while other clusters undergo proton transfer to give an ion and radical. It is found that proton transfer based structures are preferred by the (NH3)n˙(+), (H2O)n˙(+) and (HF)n˙(+) clusters while hemibonded structures are favoured by (PH3)n˙(+), (H2S)n˙(+) and (HCl)n˙(+). These trends can be attributed to the relative strengths of the molecules and molecular cations as Brønsted bases and acids, respectively, and the strength of the interaction between the ion and radical in the ion-radical clusters.

Structure and bonding in ionized water clusters.

The structure and bonding in ionized water clusters, (H2O)(n)(+) (n = 3­9), has been studied using the basin hopping search algorithm in combination with quantum chemical calculations. Initially candidate low energy isomers were generated using basin hopping in conjunction with density functional theory. Subsequently, the structures and energies were refined using second order Møller­Plesset perturbation theory and coupled cluster theory, respectively. The lowest energy isomers are found to involve proton transfer to give H(3)O(+) and a OH radical, which are more stable than isomers containing the hemibonded hydrazine-like fragment (H(2)O­OH(2)), with the calculated infrared spectra consistent with experimental data. For (H(2)O)(9)(+) the observation of a new structural motif comprising proton transfer to form H(3)O(+) and OH, but with the OH radical involved in hemibonding to another water molecule is discussed.

Computational study of the structure and electronic circular dichroism spectroscopy of blue copper proteins.

The calculation of the electronic circular dichroism (CD) spectra of the oxidized form of the blue copper proteins plastocyanin and cucumber basic protein and the relationship between the observed spectral features and the structure of the active site of the protein is investigated. Excitation energies and transition strengths are computed using multireference configuration interaction, and it is shown that computed spectra based on coordinates from the crystal structure or a single structure optimized in quantum mechanics/molecular mechanics (QM/MM) or ligand field molecular mechanics (LFMM) are qualitatively incorrect. In particular, the rotational strength of the ligand to metal charge transfer band is predicted to be too small or have the incorrect sign. By considering calculations on active site models with modified structures, it is shown that the intensity of this band is sensitive to the nonplanarity of the histidine and cysteine ligands coordinated to copper. Calculation of the ultraviolet absorption and CD spectra based upon averaging over many structures drawn from a LFMM molecular dynamics simulation are in good agreement with experiment, and superior to analogous calculations based upon structures from a classical molecular dynamics simulation. This provides evidence that the LFMM force field provides an accurate description of the molecular dynamics of these proteins.

Density of States Partitioning Method for Calculating the Free Energy of Solids.

We propose a new simulation method, which combines a cage model and a density of states partitioning technique, to compute the free energy of an arbitrary solid. The excess free energy is separated into two contributions, noninteracting and interacting. The excess free energy of the noninteracting solid is computed by partitioning its geometrical configuration space with respect to the ideal gas. This quantity depends on the lattice type and the number of molecules. The excess free energy of the interacting solid, with respect to the noninteracting solid, is calculated using density of states partitioning and a cage model. The cage model is better than the cell model in that it has a smaller configuration space and better represents the equilibrium distribution of solid configurations. Since the partition function (and hence free energy) is obtained from the density of states, which is independent of the temperature, equilibrium thermodynamic properties at any condition can be obtained by varying the density. We illustrate our method in the context of the free energy of dry ice.