Methylammonium Tetrel Halide Perovskite Ion Pairs and Their Dimers: The Interplay between the Hydrogen-, Pnictogen- and Tetrel-Bonding Interactions.

The structural stability of the extensively studied organic-inorganic hybrid methylammonium tetrel halide perovskite semiconductors, MATtX3 (MA = CH3NH3+; Tt = Ge, Sn, Pb; X = Cl, Br, I), arises as a result of non-covalent interactions between an organic cation (CH3NH3+) and an inorganic anion (TtX3-). However, the basic understanding of the underlying chemical bonding interactions in these systems that link the ionic moieties together in complex configurations is still limited. In this study, ion pair models constituting the organic and inorganic ions were regarded as the repeating units of periodic crystal systems and density functional theory simulations were performed to elucidate the nature of the non-covalent interactions between them. It is demonstrated that not only the charge-assisted N-H···X and C-H···X hydrogen bonds but also the C-N···X pnictogen bonds interact to stabilize the ion pairs and to define their geometries in the gas phase. Similar interactions are also responsible for the formation of crystalline MATtX3 in the low-temperature phase, some of which have been delineated in previous studies. In contrast, the Tt···X tetrel bonding interactions, which are hidden as coordinate bonds in the crystals, play a vital role in holding the inorganic anionic moieties (TtX3-) together. We have demonstrated that each Tt in each [CH3NH3+•TtX3-] ion pair has the capacity to donate three tetrel (σ-hole) bonds to the halides of three nearest neighbor TtX3- units, thus causing the emergence of an infinite array of 3D TtX64- octahedra in the crystalline phase. The TtX44- octahedra are corner-shared to form cage-like inorganic frameworks that host the organic cation, leading to the formation of functional tetrel halide perovskite materials that have outstanding optoelectronic properties in the solid state. We harnessed the results using the quantum theory of atoms in molecules, natural bond orbital, molecular electrostatic surface potential and independent gradient models to validate these conclusions.

The Tetrel Bond and Tetrel Halide Perovskite Semiconductors.

The ion pairs [Cs+•TtX3-] (Tt = Pb, Sn, Ge; X = I, Br, Cl) are the building blocks of all-inorganic cesium tetrel halide perovskites in 3D, CsTtX3, that are widely regarded as blockbuster materials for optoelectronic applications such as in solar cells. The 3D structures consist of an anionic inorganic tetrel halide framework stabilized by the cesium cations (Cs+). We use computational methods to show that the geometrical connectivity between the inorganic monoanions, [TtX3-]∞, that leads to the formation of the TtX64- octahedra and the 3D inorganic perovskite architecture is the result of the joint effect of polarization and coulombic forces driven by alkali and tetrel bonds. Depending on the nature and temperature phase of these perovskite systems, the Tt···X tetrel bonds are either indistinguishable or somehow distinguishable from Tt-X coordinate bonds. The calculation of the potential on the electrostatic surface of the Tt atom in molecular [Cs+•TtX3-] provides physical insight into why the negative anions [TtX3-] attract each other when in close proximity, leading to the formation of the CsTtX3 tetrel halide perovskites in the solid state. The inter-molecular (and inter-ionic) geometries, binding energies, and charge density-based topological properties of sixteen [Cs+•TtX3-] ion pairs, as well as some selected oligomers [Cs+•PbI3-]n (n = 2, 3, 4), are discussed.

An Oxysulfide Photocatalyst Evolving Hydrogen with an Apparent Quantum Efficiency of 30 % under Visible Light.

Photocatalytic water splitting is a simple means of converting solar energy into storable hydrogen energy. Narrow-band gap oxysulfide photocatalysts have attracted much attention in this regard owing to the significant visible-light absorption and relatively high stability of these compounds. However, existing materials suffer from low efficiencies due to difficulties in synthesizing these oxysulfides with suitable degrees of crystallinity and particle sizes, and in constructing effective reaction sites. The present work demonstrates the production of a Gd2 Ti2 O5 S2 (λ<650 nm) photocatalyst capable of efficiently driving photocatalytic reactions. Single-crystalline, plate-like Gd2 Ti2 O5 S2 particles with atomically ordered surfaces were synthesized by flux and chemical etching methods. Ultrafine Pt-IrO2 cocatalyst particles that promoted hydrogen (H2 ) and oxygen (O2 ) evolution reactions were subsequently loaded on the Gd2 Ti2 O5 S2 while ensuring an intimate contact by employing a microwave-heating technique. The optimized Gd2 Ti2 O5 S2 was found to evolve H2 from an aqueous methanol solution with a remarkable apparent quantum efficiency of 30 % at 420 nm. This material was also stable during O2 evolution in the presence of a sacrificial reagent. The results presented herein demonstrates a highly efficient narrow-band gap oxysulfide photocatalyst with potential applications in practical solar hydrogen production.

Insights into the carbonate effect on water oxidation over metal oxide photocatalysts/photoanodes.

Photocatalytic/photoelectrochemical water splitting using metal oxide semiconductors is a promising technology for direct and simple solar-energy conversion. The addition of carbonate salts to an aqueous reaction solution has been known to promote stoichiometric O2 evolution and H2O2 production via H2O oxidation. To elucidate the effect of carbonates, density functional theory calculations are performed to study the photoinduced H2O and H2CO3 oxidation mechanisms on TiO2 and BiVO4. The oxidation reactions proceeded via peroxide intermediates, such as H2O2 for H2O, H2C2O6 for H2CO3, and H2CO4 for the coexistence of H2O and H2CO3 molecules. Regardless of the reactant and metal oxide, the free energy changes in the four proton-coupled electron-transfer (PCET) steps of the oxidation mechanism indicate that the first PCET requires the highest energy input and is the rate-limiting step. All PCET steps of the H2O oxidation, except the second one, are more endergonic than those of the H2CO3 oxidation. The H2O reactant requires a larger energy barrier at the highest energy profile, as well as at the final state, than the H2CO3 reactant. The computational results verify that the adsorbed H2CO3 molecule is easily photo-oxidized compared with the adsorbed H2O molecule, facilitating the formation of the peroxide intermediate and improving O2 evolution and H2O2 production.

The Pnictogen Bond, Together with Other Non-Covalent Interactions, in the Rational Design of One-, Two- and Three-Dimensional Organic-Inorganic Hybrid Metal Halide Perovskite Semiconducting Materials, and Beyond.

The pnictogen bond, a somewhat overlooked supramolecular chemical synthon known since the middle of the last century, is one of the promising types of non-covalent interactions yet to be fully understood by recognizing and exploiting its properties for the rational design of novel functional materials. Its bonding modes, energy profiles, vibrational structures and charge density topologies, among others, have yet to be comprehensively delineated, both theoretically and experimentally. In this overview, attention is largely centered on the nature of nitrogen-centered pnictogen bonds found in organic-inorganic hybrid metal halide perovskites and closely related structures deposited in the Cambridge Structural Database (CSD) and the Inorganic Chemistry Structural Database (ICSD). Focusing on well-characterized structures, it is shown that it is not merely charge-assisted hydrogen bonds that stabilize the inorganic frameworks, as widely assumed and well-documented, but simultaneously nitrogen-centered pnictogen bonding, and, depending on the atomic constituents of the organic cation, other non-covalent interactions such as halogen bonding and/or tetrel bonding, are also contributors to the stabilizing of a variety of materials in the solid state. We have shown that competition between pnictogen bonding and other interactions plays an important role in determining the tilting of the MX6 (X = a halogen) octahedra of metal halide perovskites in one, two and three-dimensions. The pnictogen interactions are identified to be directional even in zero-dimensional crystals, a structural feature in many engineered ordered materials; hence an interplay between them and other non-covalent interactions drives the structure and the functional properties of perovskite materials and enabling their application in, for example, photovoltaics and optoelectronics. We have demonstrated that nitrogen in ammonium and its derivatives in many chemical systems acts as a pnictogen bond donor and contributes to conferring stability, and hence functionality, to crystalline perovskite systems. The significance of these non-covalent interactions should not be overlooked, especially when the focus is centered on the rationale design and discovery of such highly-valued materials.

The Stibium Bond or the Antimony-Centered Pnictogen Bond: The Covalently Bound Antimony Atom in Molecular Entities in Crystal Lattices as a Pnictogen Bond Donor.

A stibium bond, i.e., a non-covalent interaction formed by covalently or coordinately bound antimony, occurs in chemical systems when there is evidence of a net attractive interaction between the electrophilic region associated with an antimony atom and a nucleophile in another, or the same molecular entity. This is a pnictogen bond and are likely formed by the elements of the pnictogen family, Group 15, of the periodic table, and is an inter- or intra-molecular non-covalent interaction. This overview describes a set of illustrative crystal systems that were stabilized (at least partially) by means of stibium bonds, together with other non-covalent interactions (such as hydrogen bonds and halogen bonds), retrieved from either the Cambridge Structure Database (CSD) or the Inorganic Crystal Structure Database (ICSD). We demonstrate that these databases contain hundreds of crystal structures of various dimensions in which covalently or coordinately bound antimony atoms in molecular entities feature positive sites that productively interact with various Lewis bases containing O, N, F, Cl, Br, and I atoms in the same or different molecular entities, leading to the formation of stibium bonds, and hence, being partially responsible for the stability of the crystals. The geometric features, pro-molecular charge density isosurface topologies, and extrema of the molecular electrostatic potential model were collectively examined in some instances to illustrate the presence of Sb-centered pnictogen bonding in the representative crystal systems considered.

Chalcogen Bonding in the Molecular Dimers of WCh2 (Ch = S, Se, Te): On the Basic Understanding of the Local Interfacial and Interlayer Bonding Environment in 2D Layered Tungsten Dichalcogenides.

Layered two-dimensional transition metal dichalcogenides and their heterostructures are of current interest, owing to the diversity of their applications in many areas of materials nanoscience and technologies. With this in mind, we have examined the three molecular dimers of the tungsten dichalcogenide series, (WCh2)2 (Ch = S, Se, Te), using density functional theory to provide insight into which interactions, and their specific characteristics, are responsible for the interfacial/interlayer region in the room temperature 2H phase of WCh2 crystals. Our calculations at various levels of theory suggested that the Te···Te chalcogen bonding in (WTe2)2 is weak, whereas the Se···Se and S···S bonding interactions in (WSe2)2 and (WS2)2, respectively, are of the van der Waals type. The presence and character of Ch···Ch chalcogen bonding interactions in the dimers of (WCh2)2 are examined with a number of theoretical approaches and discussed, including charge-density-based approaches, such as the quantum theory of atoms in molecules, interaction region indicator, independent gradient model, and reduced density gradient non-covalent index approaches. The charge-density-based topological features are shown to be concordant with the results that originate from the extrema of potential on the electrostatic surfaces of WCh2 monomers. A natural bond orbital analysis has enabled us to suggest a number of weak hyperconjugative charge transfer interactions between the interacting monomers that are responsible for the geometry of the (WCh2)2 dimers at equilibrium. In addition to other features, we demonstrate that there is no so-called van der Waals gap between the monolayers in two-dimensional layered transition metal tungsten dichalcogenides, which are gapless, and that the (WCh2)2 dimers may be prototypes for a basic understanding of the physical chemistry of the chemical bonding environments associated with the local interfacial/interlayer regions in layered 2H-WCh2 nanoscale systems.

The Pnictogen Bond: The Covalently Bound Arsenic Atom in Molecular Entities in Crystals as a Pnictogen Bond Donor.

In chemical systems, the arsenic-centered pnictogen bond, or simply the arsenic bond, occurs when there is evidence of a net attractive interaction between the electrophilic region associated with a covalently or coordinately bound arsenic atom in a molecular entity and a nucleophile in another or the same molecular entity. It is the third member of the family of pnictogen bonds formed by the third atom of the pnictogen family, Group 15 of the periodic table, and is an inter- or intramolecular noncovalent interaction. In this overview, we present several illustrative crystal structures deposited into the Cambridge Structure Database (CSD) and the Inorganic Chemistry Structural Database (ICSD) during the last and current centuries to demonstrate that the arsenic atom in molecular entities has a significant ability to act as an electrophilic agent to make an attractive engagement with nucleophiles when in close vicinity, thereby forming σ-hole or π-hole interactions, and hence driving (in part, at least) the overall stability of the system's crystalline phase. This overview does not include results from theoretical simulations reported by others as none of them address the signatory details of As-centered pnictogen bonds. Rather, we aimed at highlighting the interaction modes of arsenic-centered σ- and π-holes in the rationale design of crystal lattices to demonstrate that such interactions are abundant in crystalline materials, but care has to be taken to identify them as is usually done with the much more widely known noncovalent interactions in chemical systems, halogen bonding and hydrogen bonding. We also demonstrate that As-centered pnictogen bonds are usually accompanied by other primary and secondary interactions, which reinforce their occurrence and strength in most of the crystal structures illustrated. A statistical analysis of structures deposited into the CSD was performed for each interaction type As···D (D = N, O, S, Se, Te, F, Cl, Br, I, arene's π system), thus providing insight into the typical nature of As···D interaction distances and ∠R-As···D bond angles of these interactions in crystals, where R is the remainder of the molecular entity.

The Phosphorus Bond, or the Phosphorus-Centered Pnictogen Bond: The Covalently Bound Phosphorus Atom in Molecular Entities and Crystals as a Pnictogen Bond Donor.

The phosphorus bond in chemical systems, which is an inter- or intramolecular noncovalent interaction, occurs when there is evidence of a net attractive interaction between an electrophilic region associated with a covalently or coordinately bonded phosphorus atom in a molecular entity and a nucleophile in another, or the same, molecular entity. It is the second member of the family of pnictogen bonds, formed by the second member of the pnictogen family of the periodic table. In this overview, we provide the reader with a snapshot of the nature, and possible occurrences, of phosphorus-centered pnictogen bonding in illustrative chemical crystal systems drawn from the ICSD (Inorganic Crystal Structure Database) and CSD (Cambridge Structural Database) databases, some of which date back to the latter part of the last century. The illustrative systems discussed are expected to assist as a guide to researchers in rationalizing phosphorus-centered pnictogen bonding in the rational design of molecular complexes, crystals, and materials and their subsequent characterization.

Understanding the electrochemical properties of bulk phase and surface structures of Na3TMPO4CO3 (TM = Fe, Mn, Co, Ni) from first principles calculations.

First principles calculations were performed to investigate the electrochemical performance (voltage, cycling stability, electrical conductivity, mechanical properties and safety) of the bulk phase and surface structures of Na3TMPO4CO3 (TM = Fe, Mn, Co, Ni). Na3FePO4CO3 and Na3MnPO4CO3 are estimated to be promising candidates for the cathode materials of sodium ion batteries because of the moderate voltages, good stability and high safety during the cycling process of two sodium ions per formula unit. For the purpose of improving the rate performances, Na3MnPO4CO3 was chosen as an example to explore its surface performance. The surface energies, equilibrium morphology, redox potentials and electronic conductivities of surfaces are explored in detail. The results suggest that (010), (001), (111) and (110) orientations are the dominating surfaces in the Wulff shape, while the surfaces (010) and (001) possess high second surface redox potentials, corresponding to the unsatisfactory specific capacity and ionic conductivity. Moreover, low surface band gaps are discovered in all orientations, which gives a good explanation for the enhanced electronic conductivity as a consequence of decreasing particle size. In addition, the (110), (101) and (12-1) surfaces display significantly lower surface band gaps and comparatively lower second redox potentials, thus enlarging the relative surface areas of surfaces (110), (101) and (12-1) could be an efficient methodology to further improve the specific capacity and electronic conductivity of the Na3MnPO4CO3 material.

Laser-induced fluorescence of the trans-CHBrCHO radical.

A new laser-induced fluorescence spectrum was observed in the region of 350 nm-360 nm. The spectrum was observed in the reaction between the CHBrCHBr and OH radicals and in the reaction of CHBrCHBr and CH2CHBr with atomic oxygen O(3P). The spectrum was assigned to the B̃--X̃ transition of the trans-CHBrCHO (trans-2-bromovinoxy) radical. The B̃--X̃ electronic transition energy (T0) was 28 542 cm-1, which was 242 cm-1 lower than that of the unsubstituted vinoxy radical (CH2CHO). From an analysis of the laser-induced single vibronic level fluorescence aided by ab initio calculations, some of the vibrational frequencies were assigned to the ground electronic state ν3 (C-O str.) = 1581 cm-1, ν6 (C-C str.) = 1130 cm-1, and ν8 (C-C-O bend.) = 409 cm-1. The fluorescence lifetimes of the excited B̃ state were 35 ns-75 ns, depending on the excited vibrational modes, implying that predissociation had accelerated as the energy level (v') increased.

Application of accelerated long-range corrected exchange functional [LC-DFT(2Gau)] to periodic boundary condition systems: CO adsorption on Cu(111) surface.

Several different types of density functional theory (DFT) exchange correlation functionals were applied to a periodic boundary condition (PBC) system [carbon monoxide (CO) adsorbed on Cu(111): CO/Cu(111)] and the differences in the results calculated using these functionals were compared. The exchange correlation functionals compared were those of Perdew-Burke-Ernzerhof (PBE) and those of long-range corrected density functional theory (LC-DFT), such as LC-ωPBE(2Gau) and LC-BLYP(2Gau). Solid state properties such as the partial density of states were calculated in order to elucidate the detailed adsorption mechanisms and back-bonding peculiar to the CO/Cu(111) system. In addition, our benchmark analysis of the correlations among the orbitals of CO and Cu metal using LC-DFT reasonably was in line with the experimentally observed adsorption site. The computation time was reasonable, and other numerical results were found to agree well with the experimental results and also with the theoretical results of other researchers. This suggests that the long-range Hartree-Fock exchange integral should be included to correctly predict the electronic nature of PBC systems.

A Theoretical study on the charge and discharge states of Na-ion battery cathode material, Na1+x FePO4 F.

Na2 FePO4 F is a promising cathode material for a Na-ion battery because of its high electronic capacity and good cycle performance. In this work, first principle calculations combined with cluster expansion and the Monte Carlo method have been applied to analyze the charge and discharge processes of Na2 FePO4 F by examining the voltage curve and the phase diagram. As a result of the density functional theory calculation and experimental verification with structural analysis, we found that the most stable structure of Na1.5 FePO4 F has the P21 /b11 space group, which has not been reported to date. The estimated voltage curve has two clear plateaus caused by the two-phase structure composed of P21 /b11 Na1.5 FePO4 F and Pbcn Na2 FePO4 F or Na1 FePO4 F and separated along the c-axis direction. The phase diagram shows the stability of the phase-separated structure. Considering that Na2 FePO4 F has diffusion paths in the a- and c-axis directions, Na2 FePO4 F has both innerphase and interphase diffusion paths. We suggest that the stable two-phase structure and the diffusion paths to both the innerphase and interphases are a key for the very clear plateau. We challenge to simulate a nonequilibrium state at high rate discharge with high temperature by introducing a coordinate-dependent chemical potential. The simulation shows agreement with the experimental discharge curve on the disappearance of the two plateaus. © 2018 Wiley Periodicals, Inc.

Molecular QTAIM Topology Is Sensitive to Relativistic Corrections.

The topology of the molecular electron density of benzene dithiol gold cluster complex Au4 -S-C6 H4 -S'-Au'4 changed when relativistic corrections were made and the structure was close to a minimum of the Born-Oppenheimer energy surface. Specifically, new bond paths between hydrogen atoms on the benzene ring and gold atoms appeared, indicating that there is a favorable interaction between these atoms at the relativistic level. This is consistent with the observation that gold becomes a better electron acceptor when relativistic corrections are applied. In addition to relativistic effects, here, we establish the sensitivity of molecular topology to basis sets and convergence thresholds for geometry optimization.

Laser-induced fluorescence of the CHFCHO radical and reaction of OH radicals with halogenated ethylenes.

A new laser-induced fluorescence spectrum of the 2-fluorovinoxy (CHFCHO) radical was first observed around 335 nm. The radical was produced in the reaction of an OH radical with 1,2-difluoroethylene (CHF=CHF). A single weak band was observed, which was assigned to the 00 0 band of the B̃-X̃ transition of the trans-CHFCHO radical. The B̃←X̃ electronic transition energy (T0) for trans-CHFCHO was 29 871 cm-1, which was just 3 cm-1 lower than that of its isomer, the 1-fluorovinoxy (CH2CFO) radical. The fluorescence lifetime at 29 871 cm-1 was shorter than 20 ns. This means that strong predissociation is probable at v' = 0 in the excited B̃ state of trans-CHFCHO. From an analysis of the dispersed fluorescence spectrum, some of the vibrational frequencies can be assigned for the ground electronic state: ν3 = 1557 cm-1 (C-O stretch), ν7 = 1162 cm-1 (C-C stretch), and ν8 = 541 cm-1 (CCO bend). These vibrational assignments were supported by ab initio calculations. The structure of the C-C-O skeleton and the spectroscopic character of trans-CHFCHO were close to those of CHClCHO and CH2CHO than those of CH2CFO. For the reaction of CH2=CHF with O(3P), the formation of both the regioisomeric radicals, i.e., 1- and 2-fluorovinoxy radicals, was confirmed. The regioselectivity of the oxygen atom added to the double bond of monofluoroethylene is discussed.

Do surfaces of positive electrostatic potential on different halogen derivatives in molecules attract? like attracting like!

Coulomb's law states that like charges repel, and unlike charges attract. However, it has recently been theoretically revealed that two similarly charged conducting spheres will almost always attract each other when both are in close proximity. Using multiscale first principles calculations, we illustrate practical examples of several intermolecular complexes that are formed by the consequences of attraction between positive atomic sites of similar or dissimilar electrostatic surface potential on interacting molecules. The results of the quantum theory of atoms in molecules and symmetry adapted perturbation theory support the attraction between the positive sites, characterizing the F•••X (X = F, Cl, Br) intermolecular interactions in a series of 20 binary complexes as closed-shell type, although the molecular electrostatic surface potential approach does not (a failure!). Dispersion that has an r-6 dependence, where r is the equilibrium distance of separation, is found to be the sole driving force pushing the two positive sites to attract. © 2017 Wiley Periodicals, Inc.

Halogen in materials design: Chloroammonium lead triiodide perovskite (ClNH3 PbI3 ) a dynamical bandgap semiconductor in 3D for photovoltaics.

Methylammonium lead trihalides and their derivatives are photovoltaic materials. CH3 NH3 PbI3 is the most efficient light harvester among all the known halide perovskites (PSCs). It is regarded as unsuitable for long-term stable solar cells, thus it is necessary to develop other types of PSC materials to achieve stable PSCs (Wang et al., Nat. Energy 2016, 2, 16195). Because of this, various research efforts are on-going to discover novel lead-based or lead-free single/double PSCs, which can be stable, synthesizable, transportable, abundant and efficient in solar energy conversion. Keeping these factors in mind, we report here the electronic structures, energetic stabilities and some materials properties (viz. band structures, density of states spectra and photo-carrier masses) of the PSC chloroammonium lead triiodide (ClNH3 PbI3 ). This emerges through compositional engineering that often focuses on B- and Y-site substitutions within the domain of the BMY3 PSC stoichiometry. ClNH3 PbI3 is found to be stable as orthorhombic and pseudocubic polymorphs, which are analogous with the low and high temperature polymorphs of CH3 NH3 PbI3 . The bandgap of ClNH3 PbI3 (values between 1.28 and 1.60 eV) is found to be comparable with that of CH3 NH3 PbI3 , (1.58 eV), both obtained with periodic DFT at the PBE level of theory. Spin orbit coupling is shown to have a pronounced effect on both the magnitude and character of the bandgap. The computed results show that ClNH3 PbI3 may act as a competitor for CH3 NH3 PbI3 for photovoltaics. © 2018 Wiley Periodicals, Inc.

Revealing Factors Influencing the Fluorine-Centered Non-Covalent Interactions in Some Fluorine-Substituted Molecular Complexes: Insights from First-Principles Studies.

We examine the equilibrium structure and properties of six fully or partially fluorinated hydrocarbons and several of their binary complexes using computational methods. In the monomers, the electrostatic surface of the fluorine is predicted to be either entirely negative or weakly positive. However, its lateral sites are always negative. This enables the fluorine to display an anisotropic distribution of charge density on its electrostatic surface. While this is the electrostatic surface scenario of the fluorine atom, its negative sites in some of these monomers are shown to have the potential to engage in attractive engagements with the negative site(s) on the same atom in another molecule of the same type, or a molecule of a different type, to form bimolecular complexes. This is revealed by analyzing the results of current state-of-the-art computational approaches such as DFT, together with those obtained from the quantum theory of atoms in molecules, molecular electrostatic surface potential and symmetry adapted perturbation theories. We demonstrate that the intermolecular interaction energy arising in part from the universal London dispersion, which has been underappreciated for decades, is an essential factor in explaining the attraction between the negative sites, although energy arising from polarization strengthens the extent of the intermolecular interactions in these complexes.

First-principles investigation of the Lewis acid-base adduct formation at the methylammonium lead iodide surface.

We have here performed a campaign of ab initio calculations focusing on the anchoring mechanism and adduct formation of some Lewis bases, both aliphatic and aromatic, on a PbI2-rich flat (001) methylammonium lead iodide (MAPI) surface. Our goal is to provide theoretical support to the recently reported experimental techniques of MAPI surface passivation via Lewis acid-base neutralization and similarly of MAI·PbI2·(Lewis base) adduct formation. We tested several X-donor bases (X = :N, :O, :S), paying attention to the thermodynamic stability of the final MAPI·base adducts and to their electronic properties. Factors that impact on the passivation mechanism are the directionality of the Lewis base lone pair and its enhanced/reduced overlap with MAPI Pb p orbitals, the dipole moment of the base and, similarly, the electronegativity of the X donor atom. Also non-covalent interactions, both at the surface side (intra, MAPI) and at the very interface (inter, MAPI·Lewis base), seem to contribute to the stability of the final adducts. Here we show that the thermodynamic stability does not necessarily correspond to the most effective base → acid dative bond formation. Starting from a low coverage (12.5% of the undercoordinated Pb atoms available at the surface are passivated) this paper paves the way towards the study of cooperative and steric effects among Lewis bases at higher coverages representing, to the best of our knowledge, one of the very first studies focusing on the molecular anchoring on the surfaces of this very important class of materials.

A DFT assessment of some physical properties of iodine-centered halogen bonding and other non-covalent interactions in some experimentally reported crystal geometries.

A set of six binary complexes that feature iodine-centered halogen bonding, extracted from structures deposited in the Cambridge Structure Database, has been examined computationally using density functional theory calculations with the M06-2X global hybrid, and dispersion corrected B3LYP-D3 and B97-D3, to determine their equilibrium geometries, binding energies and electronic properties. The results show that gas phase calculations are very informative in evaluating what occurs in the solid state, even though these calculations ignore the importance of lattice packing and counter ion effects. The calculated binding energies for the non-covalent interactions responsible for these complexes lie between -4.15 and -7.48 kcal mol-1 (M06-2X), which enables us to characterize them as weak-to-moderate in strength. The basis set superposition error energies are calculated to vary between 0.60 and 2.42 kcal mol-1 for all the complexes examined, even though an all-electron QZP basis set used in the analysis was of quadrupole-ζ (plus polarization) quality. Dispersion is found to have a profound effect on the binding energy of some of these complexes, and was estimated to be as large as 5.0 kcal mol-1. For one complex, the crystal geometry could not be precisely reproduced using a gas phase calculation. While both halogen- and hydrogen-bonding interactions were found competitive, they cooperate with each other to determine the stable configuration of the binary complex. The molecular electrostatic surface potential, quantum theory of atoms in molecules, and reduced density gradient non-covalent Interaction models were utilized to arrive at a fundamental understanding of the various inter- and intra-molecular molecular interactions involved, as well as some other previously-overlooked non-covalent interactions that emerge in the modelling.