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.


Introduction
Atomic nitrogen, N, the first member of the pnictogen family, is one of the major components of Earth's atmosphere. Its molecular analogue, N 2 , is a crucial constituent in feedstocks leading to the production of fertilizers and mining explosives, and high-density materials [1][2][3][4][5]. It is likely to play an ever more important role in the economy of the future feedstocks leading to the production of fertilizers and mining explosives, and high-density materials [1][2][3][4][5]. It is likely to play an ever more important role in the economy of the future as the era of the ammonia economy and "green ammonia" approaches [4,[6][7][8][9]. Nitrogen in numerous molecules has been recognized widely as a Lewis base, and thus can serve as an electron density donor D for the formation of various types of intermolecular interactions, including hydrogen bonds (HBs), tetrel bonds (TrBs), pnictogen bonds (PnBs), chalcogen bonds (ChBs) and halogen bonds (XBs), among others (Scheme 1a).
In the crystalline phase, N2 molecules bond to each other through π···N(lone-pair) interactions, forming a variety of polymorphs referred to as the α-, β-, γ-, δ-, and ε-phases. The occurrence of these non-covalent interactions is unsurprising, given that the bonding region in N2 carries a positive potential, and the electron density accumulates along and around the extension of the NN triple bond [10]. Scheme 1. (a) Schematic representation of five different types of Type-II (140° < θ < 180°) non-covalent interactions formed by covalently bonded hydrogen, tetrel, pnictogen, chalcogen, and halogen atoms in molecular entities, where R is the remainder part of the molecular entity and θ is the angular approach of the electrophile centered on the H/Tr/Ch/Pn/X atoms. (b) A ball-and-stick model of dodecane-1,12-diamine (H2N-(CH2)12-NH2), in which the ammine moiety acts both as a hydrogen bonded acceptor and a hydrogen bond donor (CSD ref: UJONUD [11]). The H···N hydrogen bonds in b) are depicted as dotted lines in green and red, with the latter are hanging contacts. Color codes: N-blue; C-gray; H-white. N forms single, double and triple bonds, and, given its significant electronegativity (χPauling = 3.04; χAllen = 3.066), often carries a negative charge and its ability to act as a Lewis base is well known. Ammines can act both as hydrogen bond donors and acceptors. Scheme 1b, for instance, shows the way the nitrogen in dodecane-1,12-diamine bonds noncovalently with its neighbors through hydrogen bonds (r(H···N) = 2.233 Å , N-H···N = 170.9°), resulting in the formation of a crystal lattice [11]. N in NH3 and ammines is a ligand in thousands of transition metal complexes, and the coordination of NH3 and its derivatives by metal ions has been extensively studied for many years, and often used to explore trends in stability constants of the transition metal ions (for example [12][13][14][15][16][17][18][19][20][21]).
There are many ammonium, diammonium and their derivatives found in solid state structures, including those that may be suitable in the development of photovoltaic and Scheme 1. (a) Schematic representation of five different types of Type-II (140 • < θ < 180 • ) non-covalent interactions formed by covalently bonded hydrogen, tetrel, pnictogen, chalcogen, and halogen atoms in molecular entities, where R is the remainder part of the molecular entity and θ is the angular approach of the electrophile centered on the H/Tr/Ch/Pn/X atoms. (b) A ball-and-stick model of dodecane-1,12-diamine (H 2 N-(CH 2 ) 12 -NH 2 ), in which the ammine moiety acts both as a hydrogen bonded acceptor and a hydrogen bond donor (CSD ref: UJONUD [11]). The H···N hydrogen bonds in b) are depicted as dotted lines in green and red, with the latter are hanging contacts. Color codes: N-blue; C-gray; H-white. N forms single, double and triple bonds, and, given its significant electronegativity (χ Pauling = 3.04; χ Allen = 3.066), often carries a negative charge and its ability to act as a Lewis base is well known. Ammines can act both as hydrogen bond donors and acceptors. Scheme 1b, for instance, shows the way the nitrogen in dodecane-1,12-diamine bonds non-covalently with its neighbors through hydrogen bonds (r(H···N) = 2.233 Å, ∠N-H···N = 170.9 • ), resulting in the formation of a crystal lattice [11]. N in NH 3 and ammines is a ligand in thousands of transition metal complexes, and the coordination of NH 3 and its derivatives by metal ions has been extensively studied for many years, and often used to explore trends in stability constants of the transition metal ions (for example [12][13][14][15][16][17][18][19][20][21]).
There are many ammonium, diammonium and their derivatives found in solid state structures, including those that may be suitable in the development of photovoltaic and optoelectronic materials (Scheme 2). The presence of cations such as these, often referred to as spacers, or additives, stabilizes the inorganic frameworks, and hence plays a structuredetermining role during the synthesis of organic-inorganic hybrid materials and is responsible for their functionality. Lower-dimensional organic-inorganic hybrid metal trihalide perovskite semiconductors are prominent examples. In these, bifunctional organoammonium cations X(CH 2 ) 2 NH 3 + (X = OH, Cl, Br, I, CN) control the in-and out-of-plane distortions of X(CH 2 ) 2 NH 3 ) 2 PbI 4 perovskites [22]. The family of 〈100〉-oriented layered perovskites with general formula (RNH 3 ) 2 A n-1 B n X 3n+1 are obtained by taking n layers along the 〈100〉 direction of the parent structure [22]. A search for the fragment "-CH 2 -CH 2 -CH 2 -NH 3 " in the Cambridge Structural Database (CSD [23,24]) led to 3656 hits; that for the fragment "-X-C-NH 3 " or "X-CH 2 -NH 3 " (X = any atom) led to 16,834 and 9509 hits, respectively, illustrating the wide occurrence of structures featuring an RNH 3 + entity. Although a few studies have discussed the importance of pnictogen bonding in the metal trihalide perovskite semiconductors (for example [25,26]), in addition to other non-covalent interactions such as hydrogen bonding and tetrel bonding [27][28][29], our exploration of the geometry of a variety of crystal structures suggests that the occurrence of nitrogen-centered pnictogen bonding in these systems is very common and that it plays an important role in their structural integrity.
This overview is therefore focused on exploring representative crystal structures deposited in the CSD and the Inorganic Chemistry Structural Database (ICSD) [30,31] in which covalently bound positively-charged nitrogen contributes in part to their structural stability through inter-or intra-molecular non-covalent interactions. In many cases, we attribute the attraction between N and the negative site as a charge-assisted nitrogencentered pnictogen bond (or simply, nitrogen bond [10]). This could be either a σ-hole centered or a π-hole centered nitrogen bonding interaction. A σ-hole on an atom A is recognized along the extension of the R-A covalent bond, and is deficient in electron density, where R is the remainder part of the molecular entity [32,33]. A π-hole can be recognized on an atom (for example, N in NO 3 − ), or on a molecular fragment (for example, the central portion of the C≡C bond in acetylene), or at the center of a delocalized system (such as an arene moiety) that has the ability to accept electron density from a lone-pair on a Lewis base (such as O in H 2 O and N in NH 3 ). Examples of this have been discussed elsewhere [10,[34][35][36]. Our survey included only single crystals in the CSD and ICSD that were free of errors and distortions, and that had an R-factor ≤ 0.1. A statistical analysis was performed using the geometric data obtained from this survey, which is presented before the conclusions section.
The nitrogen bond, or a covalently bound nitrogen-centered pnictogen bond, in chemical systems occurs when there is evidence of a net attractive interaction between the electrophilic region associated with a covalently or coordinately bound nitrogen atom in a molecular entity and a nucleophile in another, or the same, molecular entity. It is the first member of the family of pnictogen bonds formed by the first atom of the pnictogen family, Group 15, of the periodic table, and is an inter-or intra-molecular non-covalent interaction [10]. The possible occurrence of pnictogen bonds in many crystal lattices, formed by covalently bound nitrogen [10,33], phosphorous [33,35], arsenic [33,34], antimony [36] and bismuth [37], has already been discussed recently. the putative interaction [42,43]; (3) the positive nature of the electrostatic potential on the surface of N in molecular entities and its engagement with a negative site [10]; and (4) the existence of promolecular charge density-based isosurface volumes between the bonded atomic basins [44,45]. While our discussion is largely focused on the utilization of 1 and 2 in identifying pnictogen bonding in the illustrative crystal systems, the latter two properties were computed for some chosen systems to confirm the occurrence of such interactions between molecular entities that play any appreciable role in the overall stability of the crystal lattice, and hence in the functionality of these materials. Application of the four features above has been informative in rationalizing inter-and intra-molecular interactions of various kinds [10,[34][35][36]38,42] in a variety of chemical systems, and hence further demonstration is unnecessary. We note that the appearance of N-centered pnictogen bonding in the illustrative crystals is usually accompanied by hydrogen bonds. In the majority of the crystals explored, we found that these are stronger than the nitrogen bonds; it is therefore reasonable to conclude that they are usually the cause, and that the emergence of nitrogen bonding is probably an effect.
This overview is structured as follows. Following the Computational Details section, we provide examples of illustrative crystal systems featuring the wide occurrence of N-centered pnictogen bonds. This section comprises subsections that highlight 0D, 1D, 2D and 3D crystals, the dimension of the inorganic framework in the organic-organic frameworks. Prior to the Conclusion section, we include a section that shows the frequency of the various N-center pnictogen binding interactions that occur in crystals. Electron density donors of different types are considered in this section.

Computational Details
Ammonium and its derivatives (see Scheme 2) are entirely positive. This means that both the H and N atoms in NH 4 + , as well as the -CH 2 -and -CH 3 fragments in the derivatives of mono-, di-, tri-and multi-ammonium molecular entities are all positive. Depending on the nature of the (R-CH 2 -) n fragment in [R-(CH 2 ) n -NH 3 ] n+ (R = H, X, NH 3 and CH 3 , etc.), the potential on the electrostatic surface of each constituent species differs from other species. Therefore, it is not necessary to compute the electrostatic potential on the surface of the organic cation contributing to the structural integrity of the organicinorganic crystal lattices that are shown in the following section. Even so, we have chosen a few organic cations to show that the electrostatic potential on their surfaces differs from one atomic domain to the another that constitute the molecular cation.
Ab initio calculations at the MP2(FC) level theory [46,47] were performed on a selected number of organic cations in combination with Dunning's Aug-CC-pVTZ basis set system [48,49] to obtain their equilibrium geometries; FC refers to the frozen core approximation. Frequency calculations were also performed for each of them, and positive eigenvalues were found. Using fully optimized geometries, we used the resulting wavefunctions to calculate the electrostatic surface potential of the molecular entities. The Gaussian 16 suit of programs was used [50]. The 0.001 a.u. (electrons Bohr −3 ) isoelectron density envelope that arbitrarily defines the van der Waals surface of a molecular entity was used on which to calculate the potential.
The details of the Molecular Electrostatic Surface Potential (MESP) have been exhaustively described [32,33,38,[51][52][53][54][55] and will therefore not be repeated here. In essence, the MESP model features two extrema, V S,min and V S,max , called the local most minimum and the local most maximum of potential, respectively, that can be found on the electrostatic surface of a molecular entity. The sign and magnitude enable us to gain insight into the electrophilic and nucleophilic character and strength of specific regions on the surface of a molecular entity, respectively [54,55]. In particular, when V S,min or V S,max is positive (i.e., V S,min > 0 or V S,max > 0), they each represent regions that are electrophilic [10,[34][35][36]42]. Similarly, when both are negative (i.e., V S,min < 0 or V S,max < 0), they represent regions that are nucleophilic. When a positive V S,max appears along the outer extension of an R-A covalent or coordinate bond, it is called a σ-hole (as on C-X bond extensions in C 6 X 6 (X = F, Cl, Br, I)) [56][57][58][59]. A σ-hole on atom A can be positive, negative or neutral, and is an electron density deficient region, along the extension of an R-A covalent, or coordinate bond, that is more positive relative to the lateral portions of A [32,55]. However, when a V S,min > 0 is found on an R-A covalent, or coordinate, bond extension, it may represent an electron deficient region that could be positive [34]; it would not be regarded as a π-hole since its location is not perpendicular to the bond axis. A V S,min > 0 occurs, for instance, along the extension of the bond in a P≡P molecule, as well as in an As 2 , Sb 2 or Bi 2 molecule [34][35][36]. When a positive region with V S,min < 0 or V S,max < 0 is close to a region described by V S,min > 0 or V S,max > 0, the formation of a non-covalent interaction is likely to occur. There have been instances found in which interacting regions with negative V S,max [60,61] (or positive V S,max [62]) along bond extensions can engage in forming non-covalent interactions.
A π-hole is an electron density deficient region, which is generally associated with a region on the surface of a molecular entity that lies orthogonal to the bond axis (such as on the centroid region of a C 6 H 6 or C 6 X 6 molecule).
In this overview, we are interested in the characterization of pnictogen bonds between interacting entities in some illustrative crystals. To this end, we selected monomeric entities in several crystals. We analyzed the promolecular charge density-based isosurface calculated by the Independent Gradient Model, IGM-δg promolecular framework [44,45] to strengthen our identification of chemical bonded interactions as revealed using intermolecular geometry and the conceptual framework of the MESP model as described above. IGM-δg has already been shown in several studies [10,35,36,42] to be useful in providing an insightful, local description of the source of the atomic or fragmentary domain that is involved in an interaction. Within the framework of IGM, the atomic electron densities of pro-molecules are summed, and the associated atomic gradients do not interfere. This is achieved by using absolute values for the sum of the atomic gradients and eliminating any electron gradient contragradience feature. The resulting total gradient |∇ρ IGM | is the upper limit of the true gradient, and the difference δg quantifies the collapse of the net electron density gradient due to the interaction. In other words, the δg descriptor identifies the existence of the opposite sign in the components of the total electron density gradient |∇ρ(r)| due to the interaction. IGM automatically separates intramolecular and interfragment interactions in a molecular entity, which can be plotted in two or three dimensions (isosurface volumes) to determine the presence of intramolecular or intramolecular interactions. From the shape of the isosurface volume, we can infer the localized or delocalized nature of the interactions. The colors of these volumes, blue and green, represent strong and weak attraction, respectively, and red represents a repulsive interaction.
To clarify the nature of the interaction energy between the interacting units responsible for the crystal geometry, and which are the building blocks in the crystals, we have used two theoretical methods, MP2(FC) and ωB97X-D [63], in conjunction with the def2-TZVPPD pseudopotential basis set [64], to optimize the geometries of the ion-pairs [NH 3 (1) and (2), respectively. E T (ion-pair), E T (ion1), E T (ion2) and E(BSSE) are the electronic total energies of the ion-pair, ion1, ion2 and energy due to BSSE, respectively.
The quantum theory of atoms in molecules (QTAIM) has been widely used to provide insight into the chemical interactions present in molecules, molecular complexes and crystals [65][66][67]. It relies on the zero-flux boundary condition for space partitioning of atomic basins in molecular entities. The bond path topology is the graphical representation of atomic domains bonded to each other. Properties such as the charge density (ρ b ), the sign and magnitude of the Laplacian of the charge density (∇ 2 ρ b ), the total energy density (H b ) the bond critical points have been shown to be valid indicators in the process of recognition of inter-and intra-molecular interactions. ∇ 2 ρ b < 0 or ∇ 2 ρ b > 0 gives an indication of whether an interaction between atomic basins has an open-shell (region with charge concentration) or a closed-shell (region with charge depletion) character [65], although ∇ 2 ρ b > 0 has been identified for polar covalent bonds [68]. We have applied QTAIM to a few cases such as those mentioned above to clarify the assumed (misleading) perception that pnictogen bonding in the building blocks responsible for crystals is negligible and hence they do not make a significant contribution to the overall stability. Since the interaction energy of ion-pairs is electrostatically dominant, we used the empirical relationship, E b = 1 ⁄2V b , to calculate the energy of the pnictogen bond in the simplified ion-pairs used for illustrative purposes [69], where V b is the potential energy density at the bond critical point (bcp); this relationship has been utilized in many studies to calculate the interaction energy of hydrogen bonds in chemical systems (for example [70,71]).
Natural Bond Orbital (NBO)'s Second Order Perturbation Theory Analysis of the Fock Matrix in NBO Basis [72][73][74][75], which is valid within the Restricted Hartree-Fock (RHF) approximation, provides insight into the type of charge transfer delocalization energy between the donor NBO (i) and the acceptor NBO (j). We have used a few systems to estimate the second order stabilization energy E (2) associated with delocalization i → j given by Equation (3).
In Equation (3), q i is the donor orbital occupancy, ε i and ε j are diagonal elements (orbital energies) and F(i,j) is the off-diagonal NBO Fock matrix element. For this purpose, the Restricted Hartree-Fock (RHF) method in conjunction with the split valence 6-311G* basis set was used. The complexation energy, QTAIM-based interaction energy and second order charge transfer energy were evaluated on the MP2(FC), ωB97X-D and crystal geometries of the chosen ion-pairs to demonstrate the persistency of results and conclusions achieved.
Analysis and drawing of geometries of various molecular entities and crystals were performed using the Mercury 4.0 [76] and VMD [77] suite of programs. AIMAll [78] and MultiWfn [79] codes were used for calculation and analysis of MESP and QTAIM graphs, and VMD [77] was used for drawing of IGM-δg based isosurfaces.

Illustrative Crystal Systems
In what follows, we examine the bonding in representative simple ammonium and dihydrazinium (diazanediium) crystal structures, and in organic-inorganic hybrid metal halide perovskites [80]. Figure 1a, is probably one of the simplest chemical adducts known in the crystalline phase. It was deposited in the CSD by Bolte in 2021 [81]; it is a cubic crystal (space group Pa3), with lattice constant a = 7.8436 Å and Z = 4.
halide perovskites [80]. Figure 1a, is probably one of the simplest chemical adducts known in the crystalline phase. It was deposited in the CSD by Bolte in 2021 [81]; it is a cubic crystal (space group Pa3 ̅ ), with lattice constant a = 7.8436 Å and Z = 4.   [81]) and the MP2(fc)/def2-TZVPPD optimized geometry, respectively. Values (in a.u.) marked in these molecular graphs represent the charge density (ρ b ) and the Laplacian of the charge density (∇ 2 ρ b ) at N···Cl bond critical point. Bond paths are in atom color and bond critical points between atomic basins are shown as tiny spheres in green. Bond distances and bond angles are in Å and degrees, respectively. Intermolecular contacts in (a,b,e) are colored in cyan and hanging contacts colored in red. N-H···Cl hydrogen bonds are not shown between the interacting ions in the crystal in (a). Figure 1c is the geometry of the pnictogen bonds between the ion pairs [NH 3 NH 3 ] 2+ and 2Cl − in the crystal, and the space-filling model showing the overlap between the N and Cl atomic basins is depicted in Figure 1d. For clarity, the H···Cl hydrogen bonds between the triangular face of the NH 3 fragments in [NH 3 NH 3 ] and the two chloride ions are not shown. They are abundant, and the H···Cl hydrogen bonds are clearly the secondary interactions while N···Cl pnictogen bonds are primary interactions, especially when the ion pairs [NH 3 NH 3 ] 2+ and 2Cl − are in a geometrical arrangement shown in Figure 1c.

Shown in
The nature of the pnictogen bonding in crystalline [NH 3  are strongly positive, where the former are more positive than the latter. This is clearly evident in the MESP model ( Figure 2); the V S,max associated with the two N atoms are equivalent due to symmetry (V S,max = 294.6 kcal mol −1 each).
There are two types of intermolecular interactions that occur in the crystal lattice of [NH3NH3] [Cl]2: H···Cl quasi-linear hydrogen bonds (N-H···Cl = 170.5°) and linear N···Cl (N-N···Cl = 180°), Figure 1b. The H···Cl hydrogen bonds are likely to be stronger than the N···Cl pnictogen bonds based on their respective bond distances (r(H···Cl) = 2.220 Å ; r(N···Cl) = 3.061 Å ). In both cases the intermolecular distances are directional and less than the sum of the vdW radii of the respective atomic basins (rvdW(H) = 1.20 Å ; rvdW(N) = 1.66 Å ; rvdW(Cl) = 1.82 Å ). All van der Waals (vdW) radii quoted in this work are from [41]). Shown in Figure 1c is the geometry of the pnictogen bonds between the ion pairs [NH3NH3] 2+ and 2Cl − in the crystal, and the space-filling model showing the overlap between the N and Cl atomic basins is depicted in Figure 1d. For clarity, the H···Cl hydrogen bonds between the triangular face of the NH3 fragments in [NH3NH3] and the two chloride ions are not shown. They are abundant, and the H···Cl hydrogen bonds are clearly the secondary interactions while N···Cl pnictogen bonds are primary interactions, especially when the ion pairs [NH3NH3] 2+ and 2Cl − are in a geometrical arrangement shown in Figure  1c.
The  2+ are strongly positive, where the former are more positive than the latter. This is clearly evident in the MESP model ( Figure 2); the VS,max associated with the two N atoms are equivalent due to symmetry (VS,max = 294.6 kcal mol −1 each).  Although the MESP model provides a useful guide to identify the two types of intermolecular bonding interactions between the interacting molecular entities of the ion-pair [NH 3 NH 3 ][Cl] 2 , the results of the IGM-δg analysis, presented in Figure 1e, strengthens the conclusion. As expected, there are two greenish isosurface volumes appearing between interacting units [NH 3 NH 3 ] 2+ and [Cl] − that signify the possibility of pnictogen bonding; the flat nature of the isosurface volume between N and Cl indicates the delocalized nature of this non-covalent interaction and possible formation of weak H···Cl hydrogen bonds. However, the QTAIM-based molecular graphs presented in Figure 1f,g, which were obtained using the geometry of the ion-pair as in the crystal and that optimized at the MP2(FC)/def2-TZVPPD level, respectively, provide evidence of the formation of the charge assisted N···Cl pnictogen bonds in [NH 3 [Cl] − in the ion-pair in the crystal geometry is weaker than in the gas phase geometry where the intermolecular distance is markedly shorter. This explains why the bond path between N and C in Figure 1f is presented as a dotted line, and that in Figure 1g as a solid line, and the charge densities at the bcps of the corresponding bonds differ markedly from each other (ρ b 0.0144, 0.0347 and 0.0327 a.u. on the crystal, MP2(FC)/def2-TZVPPD, and ωB97XD/def2-TZVPPD geometries of the ion-pair, respectively). Similarly, ∇ 2 ρ b at the N···Cl bcps of the corresponding geometries were 0.0388, 0.1422 and 0.1359 a.u., respectively. Both properties indicate the closed-shell nature of the interaction, evidenced by ∇ 2 ρ b > 0 at the N···Cl bcp, and small value of ρ b at the same bcp [65,67].
We have calculated the interaction energy between the ion pairs [NH 3 [Cl] − as another fragment to assess the interaction of a single charge assisted N···Cl pnictogen bond. The uncorrected and basis-set-corrected interaction energies, E int and E int (BSSE), respectively, were calculated to be -145.64 and -144.10 kcal mol −1 , respectively. The effect of BSSE on energy is marginal; the large interaction energy is largely due to cation-anion electrostatic interaction. This can be effectively rationalized based on the interaction energy of each N···Cl pnictogen bond in the ion-pair evaluated using the empirical relationship, E b = 1 ⁄2V b , which was -11.45, -10.10 and -2.60 kcal mol −1 when the geometry of the ion-pair optimized with MP2(FC)/def2-TZVPPD, ωB97XD/def2-TZVPPD and that in the crystal was used, respectively. These energies correlate well with the N···Cl intermolecular bond distance of 2.572, 2.595 and 3.065 Å in the corresponding ion-pair geometry, respectively.
Second order NBO analysis at RHF/6-311G* enabled us to demonstrate that the N···Cl pnictogen bond in NH 3 NH 3 Cl 2 is indeed stronger than the N-H···Cl hydrogen bonds. Because the latter ones were very weak, they were unidentifiable via the bond path and bcp topologies of QTAIM (Figure 1f When the crystal geometry of the ion-pair was used, the two types of non-covalent interactions remain identifiable. The only difference is that the magnitude of E (2) is greatly reduced. This is understandable since the N···Cl intermolecular distance in the gas phase geometry is shorter than that in the crystal geometry. The E (2) associated with [LP(4)Cl → BD*(1) N-N] and [LP(4)Cl → RY*(1)N], where RY refers to a Rydberg type orbital, which accompany each of the two pnictogen bonds, were 4.33 and 3.03 kcal mol −1 , respectively. The charge transfer interactions associated with the N-H···Cl hydrogen bonds were about 0.1 kcal mol −1 . These results show that the pnictogen bonds in the [NH 3 NH 3 ][Cl] 2 ion-pair are stronger than the hydrogen bonds. Since this study was primarily centered on pnictogen bonds, we did not explore the conformational details of hydrogen bonds, such as those shown in Figure 1b, that also play an important role in the effective stabilization of the crystal lattice.

Methylammonium Lead Halide Perovskites
Crystals containing the methylammonium action (MA, CH 3 -NH 3 + ) have been extensively studied, especially in the field of metal halide perovskites for applications in optoelectronics [26,83,84]. Monovalent MA is entirely electrophilic. A simple CSD search for this organic cation yielded 931 hits, showing its prominent role in the design of materials. Along with its derivatives, it appears prominently in the lead trihalide perovskites, viz. CH 3 NH 3 MX 3 (M = Pb, Sn, Ge; X = Cl, Br, I). These materials have a band gap in the range 1.2-3.5 eV [26,84] in the far infrared to visible region, making some of them suitable for application in solar cells. However, the actual role of the organic cation in CH 3 NH 3 MX 3 perovskites is not yet fully understood [85], despite numerous experimental and computational studies on this topic that have been reported.
It has been suggested that the stability of the crystal structure of CH 3 NH 3 MX 3 single perovskites is driven by H···X hydrogen bonds formed between the ammonium guest and the halide of the MX 6 octahedra [25,26]. However, we have recently argued [25,26] that this is incorrect for orthorhombic CH 3 NH 3 MX 3 (X = Br, I), since the overall geometric stability of these crystals arises from an environment consisting of an interplay between hydrogen, carbon and pnictogen bonded interactions. These are collectively responsible for the octahedral tilt in the low-temperature Pnma orthorhombic (o) phase of the CH 3 NH 3 PbI 3 system. We show below that this feature also holds for other members of the same family.
The structures of o-CH 3 ND 3 PbCl 3 , o-CH 3 ND 3 PbBr 3 , and o-CD 6 ND 6 PbCl 3 are shown in Figure 3a-c. They are stable at low temperature and adopt an orthorhombic crystal symmetry with tilting of the PbX 6 octahedra with the organic cation trapped inside the inorganic cage. There are a variety of N-H···X and C-H···X hydrogen bonds in each of these crystals; these are probably the primary interactions.   The primary intermolecular interactions in these crystals are charge-assisted, and strong. There are also secondary interactions that appear as pnictogen bonds and tetrel bonds (N···X and C···X (X = Cl, Br, I), respectively). These are somewhat longer, and therefore weaker, than the primary interactions. Based upon the bond distances, the stability of these interactions is in the order (N)D···X > (C)H/D···X > N···X > C···X (X = Br, I). For instance, in o-CD 3 ND 3 PbI 3 [88] [87] and o-CD 3 ND 3 PbI 3 [88], this is not the case in o-CH 3 ND 3 PbCl 3 [86] because the N···Cl and C···Cl pnictogen-and tetrel-bonds in the crystal appear alternately; hence the two orientations of the MA cation in the crystal have different types of N···Cl and C···Cl bond distances (cf. Figure 3a). Although the C···Cl tetrel bond in both the orientations are somewhat more directional (∠N-C···Cl = 165.4 • or 162.9 • ) than the N···Cl pnictogen bonds (∠C-N···Cl = 157.1 • or 152.9 • ), the N···Cl and C···Cl bond distances associated with one orientation are 3.460 and 3.571 Å, and those with the other orientation are 3.635 Å and 3.421 Å (Figure 3a). From this and, given that directional interactions are stronger than weakly directional interactions, it is difficult to determine whether the N···Cl pnictogen bonds are stronger than C···Cl tetrel bonds in o-CH 3 ND 3 PbCl 3 .
The presence of hydrogen, carbon and pnictogen bonding interactions in crystalline o-CH 3 NH 3 PbX 3 has been confirmed by QTAIM and RDG analyses reported elsewhere [25,26]. The results are shown in Figure 4 for o-CH 3 NH 3 PbBr 3 and o-CH 3 NH 3 PbI 3 . Although the principal bonding features are evident from an analysis of the intermolecular geometry alone, the bond path and bcp, as well as the isosurface topologies of QTAIM and RDG suggested that each H atom in the ammonium/methyl head in MA could be involved with one/two Br atoms to form one/two (N)H···Br/(C)H···Br hydrogen bonds in o-CH 3 NH 3 PbBr 3 described by the bcps (tiny spheres in green) between the bonded atomic basins. Bond paths between N and C in MA and the Br sites are also apparent, indicating the presence of N···Br and C···Br pnictogen-and tetrel-bonds in the crystal.   [25,26]. The PBE relaxed geometry for o-CH 3 NH 3 PbBr 3 and the neutron diffraction geometry (T ≈ 100 K [89]) for o-CH 3 NH 3 PbI 3 were used. Atoms: iodine-purple, carbondark-grey, nitrogen-deep-blue, hydrogen-white-grey; bond paths are shown as solid and dotted lines in atom color; (3, -1) bond critical points are shown as tiny blue and red spheres between bonded atomic basins. Dotted lines (bond paths) indicate the presence of non-covalent interactions between atoms that are liked with each other. Values of ρ b are in a.u.
Because the (N)H···I and (C)H···I hydrogen bonds in o-CH 3 NH 3 PbI 3 are longer than those in o-CH 3 NH 3 PbBr 3 , the QTAIM results suggest the formation of a single (N)H···I hydrogen bond by each H atom of the -NH 3 group and at least two (C)H···I bonds formed by each of the H atoms of -CH 3 group of MA (Figure 4b). The inward curvature of the bond paths between N or C in MA and the nearest I sites indicates that the N···I and C···I pnictogen-and tetrel-bonds in o-CH 3 NH 3 PbI 3 are dispersive and weak. The nature of the bond path topologies associated with these interactions is very similar to that found for o-CH 3 NH 3 PbBr 3 (see Figure 4a,b). The charge density values depicted in Figure 4 suggest that they are very small at the bcps and are typical of non-covalent interactions (ρ b < 0.05 au). It is also apparent that ρ b at the N···I and C···I bcps are smaller than those at the (N)H···I and (C)H···I bcps, confirming that first two are weaker than the last two. These analyses highlight the co-existence of several types of non-covalent interactions in these halide perovskite structures, including I···I ( Figure 4a) and Br···Br (Figure 4a) halogen-halogen interactions.

Methylammonium and Deutero-Methylammonium Halides
Two different views of the crystal structure of perdeuteromethylammonium iodide, [CD 3 ND 3 ][I] (CSD ref. KUBNOK [90]), reported using neutron diffraction measurements on a powder sample, are shown in Figure 5a,b. It crystallizes in the Pbma space group with Z = 4. The P 4 /nmm tetragonal crystal of the same compound is known, with Z = 2, but the D atom positions were not accurately determined (CSD ref. KUBNOK01 [90] (Figure 1), it appears that the D···I and H···I deuterium and hydrogen bonds, N···I pnictogen bonds and C···I tetrel bonds play a structure-determining role in [CD 3  From Figure 5c, it is apparent that the N···I pnictogen bond and the C···I tetrel bond in the perdeuteromethylammonium iodide crystal (CSD ref: KUBNOK, space group Pbma) are near-linear interactions, with ∠C-N···I = 179.3 • and ∠N-C···I = 178.7 • . The former are stronger than the latter (r(N···I) = 3.513 Å vs. r(C···I) = 3.838 Å), a bond distance which is smaller than the sum of the vdW radii N and I, 3.70 Å; the value of r(C···I) is marginally longer than the sum of the vdW radii of C and I, 3.81 Å, which suggests they are dispersive in nature. By contrast, the (N)D···I and (C)D···I deuterium bonds are inequivalent because the -ND 3 group is more acidic than the -CD 3 group in d 6  The attractive interactions between CH 3 NH 3 + and I − accord with the results of the MESP model, Figure 6. N is more electropositive than C in CH 3 NH 3 + (V S,max = 158.3 kcal mol −1 for N and 119.7 kcal mol −1 for C), and N interacts more strongly than C with I − . Both N···I and C···I are directional interactions. Their presence is evident in the IGM-δg plot (Figure 5d), in which the isosurfaces are thicker and wider between N and I than that between C and I.
[NH3NH3] [Cl]2 crystal (Figure 1), it appears that the D···I and H···I deuterium and hydrogen bonds, N···I pnictogen bonds and C···I tetrel bonds play a structure-determining role in [CD3ND3][I] and [CH3NH3][I]. The dotted lines in cyan between atomic basins in Figure  5b,c indicate the presence of the latter interactions in [CD3ND3][I], and are substantially longer than those found in [NH3NH3] [Cl]2; this is expected since the sum of the vdW radii of C and I, as well as that of N and I, is greater than that of N and Cl (rvdW(N) = 1.66 Å; rvdW(C) = 1.77 Å ; rvdW(Cl) = 1.82 Å ; rvdW(I) = 2.04 Å ).    Table 1. Clearly, the interaction energy is very large regardless of the geometry of the ion-pair used, and that the ionpair with the longer intermolecular distances is predicted to have the smaller interaction energy, as expected. Specifically, the interaction energy of the pnictogen bonded ion-pair H3C-H3N···I was found be larger than that of the carbon bonded ion-pair H3N-H3C···I present in [CH3NH3][I], and is typical of a Mulliken inner complex [94]. The strength of the interaction energy found for these systems is comparable to that in analogous systems reported in previous studies [85].
The strength of the carbon bond in the conformer shown in Figure 5f based on Eb = ½ Vb gave values of −2.79 and −0.75 kcal mol −1 when evaluated on the ωB97X-D and crystal geometries of the ion-pair, respectively. Similarly, the strength of the pnictogen bond in the conformer shown in Figure 5e was −4.64 and −1.79 kcal mol −1 when evaluated on the ωB97X-D and crystal geometries of the ion-pair, respectively. Since there was no bond path between the H atoms of -NH3 and -CH3 fragments of the CH3NH3 cation and the interacting iodide anion, the strength of hydrogen bonds between them could not be evaluated. These results suggest that pnictogen bond is the strongest of all interactions investigated.    [I] was appreciably shorter than that in the crystal. However, they are indeed less than the typical value of ρ b (ρ b < 0.05 au) often reported for non-covalent interactions [92,93]. The Laplacian of the charge density (∇ 2 ρ b ) values at the N···I (or C···I) bcp are also unequal and are positive (∇ 2 ρ b > 0), and hence attest to the closed-shell nature of these interactions [65][66][67]. Although the IGM-δg plot indicated N···I and C···I pnictogen and tetrel bonding interactions (as well as (N)D···I and (C)D···I deuterium bonds) between the interacting monomers in the ion-pair [CD 3 Table 1. Clearly, the interaction energy is very large regardless of the geometry of the ion-pair used, and that the ion-pair with the longer intermolecular distances is predicted to have the smaller interaction energy, as expected. Specifically, the interaction energy of the pnictogen bonded ion-pair H 3 C-H 3 N···I was found be larger than that of the carbon bonded ion-pair H 3

N-H 3 C···I present in [CH 3 NH 3 ][I]
, and is typical of a Mulliken inner complex [94]. The strength of the interaction energy found for these systems is comparable to that in analogous systems reported in previous studies [85].   The strength of the carbon bond in the conformer shown in Figure 5f based on E b = 1 ⁄2V b gave values of −2.79 and −0.75 kcal mol −1 when evaluated on the ωB97X-D and crystal geometries of the ion-pair, respectively. Similarly, the strength of the pnictogen bond in the conformer shown in Figure 5e was −4.64 and −1.79 kcal mol −1 when evaluated on the ωB97X-D and crystal geometries of the ion-pair, respectively. Since there was no bond path between the H atoms of -NH 3 and -CH 3 fragments of the CH 3 NH 3 cation and the interacting iodide anion, the strength of hydrogen bonds between them could not be evaluated. These results suggest that pnictogen bond is the strongest of all interactions investigated.
Zeroth-order symmetry-adapted perturbation theory (SPAT0 [95]) calculations were performed for the [CH 3 NH 3 ][I] adduct, on the structure extracted from the crystal and on the ωB97X-D and MP2(FC) optimized geometry, in conjunction with the def2-TZVPPD-RI basis set. This approach dissects the interaction energy into four major components arising from electrostatics (E els ), exchange-repulsion (E exch ), induction (E ind ) and dispersion (E disp ). The results are summarized in Table 2. As can be seen, the total interaction energies E int (SAPT0) obtained using SAPT0 are in reasonable agreement with those computed using the counterpoise method of the SCF (Self-Consistent Field) approaches (ωB97X-D and MP2(FC)). Inspection of the dissected components suggests that the energy due to electrostatics is significant as expected given that [CH 3 NH 3 ][I] is an ion-pair and hence coulombic attraction plays a very crucial role in stabilizing the intermolecular interaction. Since the intermolecular distances in the adducts in the crystal geometry are longer than those in the ωB97X-D and MP2(FC) optimized geometries, the interaction between the monomers becomes stronger in the latter case than in the former. Concomitant with this is the increase in each of the dissected components (see Table 2).
The collective contribution due to induction and dispersion cannot be overlooked since they contribute appreciably to the net interaction energy E(SAPT0). Whereas E exch is repulsive, and its magnitude is comparable with the induction energy (viz. 4.4 kcal mol −1 for H 3

N-H 3 C···I) when the crystal geometries of the H 3 C-H 3 N···I and H 3 N-H 3 C···I ion-
pairs were used, the net interaction energy E(SAPT0) is then due to energies arising from electrostatics and dispersion. However, when the ωB97X-D and MP2(FC) geometries of the H-containing ion-pairs were used, each of the four components responsible for E int (SAPT0) are significantly affected. This is not very surprising since the intermolecular distances in H 3 C-H 3 N···I and H 3 N-H 3 C···I were significantly reduced in the gas phase structure. This has caused an increase in E els , E ind and E disp , with concomitant increase in E exch . Even so, the E int (SAPT0) value calculated on the ωB97X-D and MP2(FC) geometries of the two ion-pairs are in reasonable agreement with the corresponding E int (ωB97X-D) and E int ([MP2(FC)]), respectively (see Table 2). Table 2. Comparison of total and dissected energy SAPT0 components with the SCF based interaction energies of H 3 C-H 3 N···I and H 3 N-H 3 C···I ion-pairs, obtained on the fixed crystal geometry, and the MP2(FC)/TZVPPD and ωB97X-D/TZVPPD optimized geometries. Values are in kcal mol −1 . The strengths of the pnictogen bond in H 3 C-H 3 N···I and the tetrel bond in H 3 N-H 3 C···I are significant when compared to the N-H···I and C-H···I hydrogen bonds in the corresponding ion-pairs. This is also evidenced by our second order analysis of the charge transfer delocalization energy [72][73][74][75]96]. This was carried out at the RHF/6-311G* level using the MP2(FC) and ωB97X-D optimized geometries. Our results are summarized in Table 3. The data show that for the formation of the N···I pnictogen bond and N-H···I hydrogen bonds, the lone-pair(4) on the iodide anion, LP(4)I, interacts attractively via a charge transfer process with the σ* type C-N and N-H anti-bonding orbitals BD* of the organic cation NH 3 CH 3 + , respectively. The magnitude of second order energy E (2) associated with the two types of charge transfer delocalizations is 4.51 and 1.14 kcal mol −1 , respectively. In addition, we have also identified that the same lone-pair on I is involved in a charge transfer delocalization with a Rydberg type orbital RY*(1) on N with an E (2) of 2.86 kcal mol −1 . There are, in fact, three weak N-H···I hydrogen bonds, in which LP(2), LP(3) and LP(4) are simultaneously involved in hyperconjugative charge transfer interactions with the three BD*(N-H) orbitals. Table 3. RHF/6-311G* level second-order perturbative estimates of 'donor-acceptor' (bond-antibond) interaction energies, E (2) (in kcal mol −1 ), in the NBO basis, obtained for H 3 C-H 3 N···I and H 3 N-H 3 C···I using their MP2(FC)/TZVPPD and ωB97X-D/TZVPPD level optimized geometries.

Interaction
Type Donor NBO Acceptor NBO The nature of the hyperconjugative charge transfer delocalization is very similar between the interacting units CH 3 NH 3 + and I − , leading to the formation of the ion-pair H 3 N-H 3 C···I. However, in this case, the three C-H···I hydrogen bonds are markedly weaker (compared to the N-H···I hydrogen bonds), with E (2) energies close to 0.12 kcal mol −1 . By contrast, the charge transfer delocalization is very pronounced between LP(4)I and BD*(1)C- The results given above were obtained on the MP2(FC) geometries; they are also evident when the ωB97X-D geometries of the ion-pair was used, and E (2) values are slightly smaller, in accordance with the longer intermolecular distance between CH 3 NH 3 + and I − . These results suggest that the triangular face in -NH 3 or -CH 3 is simultaneously engaged with the iodide anion, leading to the stability of the H 3 C-H 3 N···I and H 3 N-H 3 C···I ionpairs, respectively. While we did not explore the conformational space, it is no doubt that the crystal also contains a significant number of C-H···I and N-H···I hydrogen bonds of another variety.

Ammonium Cyanate
The transformation of ammonium cyanate into urea was first studied over 170 years ago by Wöhler and Liebig [97]. Knowledge of the crystal structure of ammonium cyanate is a prerequisite for a basic understanding of the nature of its solid-state reaction. The neutron powder diffraction data have enabled MacLean et al. [97] to demonstrate the formation of N−H···N hydrogen bonds between an NH 4 + ion to four cyanate N atoms at alternate corners of a distorted cube, overturning their initial own view of the presence of N−H···O hydrogen bonds to cyanate O atoms at the other four corners.
We investigated the intermolecular geometry of the NH 4 CNO crystal with MESP and IGM-δg calculations. The results we obtained support both proposals of the authors [97]. Three different views of the 2 × 2 × 2 supercell crystal structure of tetragonal ammonium cyanate (space group P 4 /nmm) are shown in Figure 7a Figure 8c,d are the IGM-δg based promolecular charge density plots between NH 4 + and CNO − ions, obtained using a geometry extracted from the crystal. From the size of the IGM-δg volumes obtained using two different isovalues (0.015 and 0.005 a.u.), it is clear that the charge-assisted N···O pnictogen bond is moderately strong given its presence even at the higher isovalue of 0.015 a.u. The isosurface corresponding to the H···O close contacts that appears at low isovalues indicates that the H···O hydrogen bonds are weaker than the N···O pnictogen bonds (see Figure 8c,d). On the other hand, the isosurfaces between N and H are bluish as a result of more electron density concentration in the bonding region (cf. Figure 8e). The presence of (C)π···π(C) interactions in the crystal is also confirmed by the IGM-δg analysis, with a contact distance of 3.644 Å. However, these interactions are probably the weakest of all the interactions identified in the crystal, as evidenced by the flat greenish isosurface between a pair of CNO − ions (see Figure 8e).

Shown in
To further clarify the presence of pnictogen bonding in NH 4 CNO, we have used the geometry of the ion-pair found in the crystal and performed QTAIM calculations on it at the MP2(FC)/def2-TZVPPD level. Similar to other ion-pairs discussed above and below, QTAIM did not indicate the formation of hydrogen bonding in the conformation shown in Figure 7g. However, it indeed revealed the presence of an N···O pnictogen bonding between the N and O atomic basins of the two interacting ion pairs. There is a bond path, and bond critical point between the bonded atomic basins, shown as dotted line in atom color between N and O. The ρ b and ∇ 2 ρ b at the N···O bcp is as small as 0.0074 and 0.0315 au. The empirical expression, E b = 1 ⁄2V b , gave values of −1.51 kcal mol −1 in the corresponding geometries, respectively, suggesting the possibility of a weak interaction between N and O in NH 4 CNO. Of course, this interaction energy is substantially smaller than the complexation energy found at the MP2(fc)/def2-TZVPPD level: the uncorrected and BSSE corrected energies of −88.16 and −87.80 kcal mol −1 , respectively, arising primarily from the electrostatic interaction between two ionic moieties that cause the two ionic domains to come into close proximity. NBO's second order analysis with RHF/6-311G* gave an indication of the occurrence of a LP(3)O → BD(1)* N-H hyperconjugative interaction, with an E (2) of 0.56 kcal mol −1 ; however, there was no such appreciable hyperconjugative interaction revealed that could represent hydrogen bonding between NH 4 + and CNO − for the conformer shown in Figure 7g.

The Ammonium Halides
Referring to their 2007 work [98], Nelyubina et al. have shown in a subsequent 2014 study [99] that N···Cl close contacts could be present in crystalline NH 4 Cl. Experimental electron density data derived from X-ray diffraction measurements enabled them to provide an accurate description of the interionic bonding in NH 4 Cl, and gave results that were consistent with quantum chemistry based-calculations [98].
The crystal structure of NH 4 Cl is shown in Figure 9a-c, and features H···Cl hydrogen bonding and N···Cl pnictogen bonding. The first of these interactions is the stronger, as revealed from the bond distances (r(H···Cl) = 2.309 Å and r(N···Cl) = 3.318 Å), and both are linear (∠N-H···Cl = 180.0 • and ∠H-N···Cl = 180.0 • ). From the local geometries shown in Figure 9d,e, it is clear that N is a donor of four σ-hole-centered pnictogen bonds, and simultaneously four hydrogen bond donors; this means that each NH 4 + unit is surrounded by eight nearest-neighbor Cl − anions in the crystal lattice. The former are engaged attractively with four nearest-neighbor chloride ions in forming four σ-hole centered charge-assisted H···Cl hydrogen bonds. Similarly, the four σ-holes on N in NH 4 + are linked with four nearest-neighbor Cl ions, thereby forming four σ-hole-centered pnictogen bonds. This shows that the Cl − ion is the versatile acceptor of eight non-covalent bonds (four H-bonds and four N-bonds), Figure 9f. Our RHF/6-311G* level second order NBO analysis suggested that delocalization of charge transfer promotes attraction between the interacting units. In particular, this analysis suggests two possible interactions, H···Cl and N···Cl, that stabilize the ion-pair configuration shown in Figure 9g. There are three equivalent close contacts of the former type, each with an E (2)   Nelyubina et al. [99] have estimated the energy of the pnictogen bonds to vary between 1.2 and 1.5 kcal mol −1 , so each of them accounts for 20 to 25% of the interaction energy per NH 4 + moiety in NH 4 Cl. As mentioned above, our view is consistent with that of these authors: the hydrogen bonds in the crystal are the primary interactions. They are the cause of the occurrence of the pnictogen bonds with the ammonium cation. Clearly, there is a straightforward interplay between these two types of bonding interactions that not only provides stability to the crystal structure, but also assists in the determination of phase transition behavior, ferroelastoelectric and other physical properties.
QTAIM calculations for the NH 4 Cl ion-pair were performed for both the crystal geometry and the MP2(FC)/def2-TZVPPD and ωB97X-D/def2-TZVPPD optimized geometries. As for the CH 3 NH 3 I case, we also found that the intermolecular distance between N and Cl in NH 4  that are associated with E (2) of 4.88, 0.30 and 2.66 kcal mol −1 , respectively. These insights collected using the MP2/def2-TZVPPD gas phase optimized geometry are consistent with that gained when the ion-pair geometry of the crystal was used. In this case, the strength of the hyperconjugative interaction is reduced because the interaction between the ions is longer ranged. Therefore, the strength of the hydrogen bond interactions was predicted to be of the van der Waals type, with an E (2)  Similarly, as found for NH4Cl (see above), the RH/6-311G* level NBO's second order analysis performed on the ωB97X-D/def2-TZVPPD level fully-relaxed geometry of the ion-pair NH4Br suggested the possible occurrence of both N···Br pnictogen bonding and H···Br hydrogen bonding. It provided evidence of three equivalent hydrogen bonds formed between the three H atom on the triangular face of NH4 + and the Br anion (see A number of structures of the H-and D-derivatives of ammonium iodide NH 4 I (P43m, Pm3m, Fm3m and P 4 /nmm), have also been deposited in the ICSD. Ammonium iodide undergoes a I-II-IV-V sequence of phase transitions at room temperature under pressure. The I-II phase transition occurs at P I-II = 0.5 kbar [104], while there is no reliable information on the II-IV transition pressure; it was speculated that P II-IV = 27 kbar on the basis of a subtle change observed in Raman spectra. Phase transition into the phase V was observed at P = 5 4 kbar, which was also observed in the Raman spectra [103]. The tetragonal crystal HP-IV (high pressure phase, 8600 MPa) contains accurate information on the D atom positions (ICSD ref. 89094) [103]. Balagurov et al. found that ND 4 I(V) has the same structure as the low-temperature phase ND 4 I(III)-a tetragonal structure with antiparallel ordering of ammonium ions, space group P4/nmm [100]. As in the P 4 /nmm structure of ND 4 Br, we found that the angles of interaction in the crystal structure of polymorphic phase IV of NH 4 I, Figure 11a  The IGM-δg plots shown in Figure 11g,h confirm the presence of all these interactions. The presence of the H···I hydrogen bonds in NH 4 I is indicated by the thick circular volumes in greenish-blue, whereas the N···I pnictogen bonds manifest as triangle-like isosurface volumes in deep-green; both appear at isovalues of 0.010 and 0.015 a.u. The flat isosurface volumes in green corresponding to the I···I contacts were obtained with low isovalues because they are weak (Figure 11h). That N and I atomic basins in NH 4 I are bonded can be inferred from the MESP model of NH 4 + (cf. Figure 8a). The gas-phase relaxation of the geometry of the NH 4 + I − ion-pair gave N···I inter-ionic distances of 3.103 and 3.055 Å with ωB97X-D/def2-TZVPPD and MP2(FC)/def2-TZVPPD, respectively. Indeed, these are very close to what was found in the geometry of the same ion-pair in the crystal [105] (r(H 4 N···I) = 3.374 Å). However, the gas-phase geometry was linear, whereas that in the crystal is quasi-linear. This is not very surprising since the H···I hydrogen bonds compete with pnictogen bonds in the crystal, and each ammonium cation is involved in forming eight non-bonded interactions with eight nearest-neighbor iodide ions (see Figure 11f Figure 11i shows the QTAIM-based molecular graph that provides evidence of the N···I pnictogen bonded interaction in ND 4 + I − , obtained using the crystal geometry of the ion-pair at MP2(FC)/def2-TZVPPD, and is reminiscent of the gasphase NH 4 + I − geometries as well (not shown). E b = 1 2 V b gave values of −4.71, −5.52 and −2.57 kcal mol −1 as the strength of the N···I pnictogen bond in the ωB97X-D/def2-TZVPPD, MP2(FC)/def2-TZVPPD and crystal geometries of the ion pair, respectively, the latter calculated on the ND 4 + I − geometry with MP2(FC)/def2-TZVPPD. NBO's second order analysis performed using the ωB97X-D/def2-TZVPPD geometry of the NH 4 + I − ion pair suggested the occurrence of both H···I hydrogen bonding and N···I pnictogen bonding interactions between the interacting ions in NH 4 + I − . The former is several times weaker than the latter. This conclusion is based on the charge transfer delocal-ization energy; E (2) corresponding to the delocalizations, LP(4)I → BD*(1) N-H], and [LP(4)I → RY*(1)N], associated with the N···I pnictogen bonding were 3.73 and 2.44 kcal mol −1 , respectively. There were three delocalizations of the type, LP(4)I → BD*(1) H-N], associated with three equivalent H···I hydrogen bonds, each with an E (2) of 1.1 kcal mol −1 . Other lone-pairs on I are also involved in hyperconjugative interactions, explaining both the interaction types noted above, but E (2) for them were < 0.75 kcal mol −1 . When the crystal geometry of the ion-pair was used, E (2) for two types of charge transfer delocalizations associated N···I pnictogen bonding were 2.20 and 1.52 kcal mol −1 , respectively, and that associated with each of the three hydrogen bonds was 0.13 kcal mol −1 .

The Crystal Structure of ([I 44 Pb 18 ][CH 3 NH 3 ]
Fateev and coworkers [106] reported that solvents such as DMSO, DMF and γ-butyrolactone (GBL) strongly influence the process of metal halide perovskite crystallization because of the formation of intermediate adducts with different structures and morphologies. Based on these observations, they were able to synthesize new adducts of leadbased halide perovskites and GBL with either an unusual cluster structure, [(MA) 8 a and b directions. The intermolecular distance between a pair of I atoms separating the inorganic frameworks is 4.081 Å and the interactions are linear (∠Pb-I···I = 180.0 • ). These I···I contacts are probably tertiary interactions, and are close to twice the sum of the vdW radius of the atoms, 4.08 Å. They arise as the consequence of primary (H···I hydrogen bonds) and secondary (N···I pnictogen bonds) interactions in the crystal lattice (Figure 13a,b).

Pnictogen Bond in 2D Functional Crystals of Metal Halide Perovskites: Derivatives of Ammonium as Pnictogen Bond Donors
Billing and Lemmerer [108] reported a series of layered crystal structures of the inorganic-organic perovskite-type hybrids, stable at room temperature, with general formula [(CnH2n+1NH3)2PbI4] (n = 12, 14, 16 and 18). Each of them displays two reversible phase transitions above room temperature, with phases labelled III, II and I. The transition from phase III to phase II was first-order and was related to a change in the conformation of the alkylammonium chains and a shift of the inorganic layers relative to each other. Depending on the orientation of the organic cation, a metal-coordinated iodide can act as a receptor of at least three pnictogen bonds, with r(N···I) = 3.894 Å and ∠C-N···I = 171.1 • . They are quasi-linear Type-IIa interactions since primary N-H···I hydrogen bonds formed by the ammonium fragment of the organic cation are simultaneously engaged attractively with the nearest-neighbor I − sites of the inorganic cage structure. They appear in different flavors, with intermolecular bond distances varying between 2.80 and 3.50 Å.

Pnictogen Bond in 2D Functional Crystals of Metal Halide Perovskites: Derivatives of Ammonium as Pnictogen Bond Donors
Billing and Lemmerer [108] reported a series of layered crystal structures of the inorganic-organic perovskite-type hybrids, stable at room temperature, with general formula [(C n H 2n+1 NH 3 ) 2 PbI 4 ] (n = 12, 14, 16 and 18). Each of them displays two reversible phase transitions above room temperature, with phases labelled III, II and I. The transition from phase III to phase II was first-order and was related to a change in the conformation of the alkylammonium chains and a shift of the inorganic layers relative to each other.
The crystal structure of one of these perovskites, with n = 16, is shown in Figure 14a.   Mixed polyhedral (inorganic) and ball-and-stick (organic) models of the crystals of (a) [(CH 3 (CH 2 ) 15  The pattern of intermolecular interactions in the system referred to above is similar to that in the crystal structure with odd n, as in [(C n H 2n+1 NH 3 ) 2 PbI 4 ] (n = 9), Figure 14d. In this case the N···I nitrogen bond formed by the terminal N atom of the ammonium moiety of CH 3 (CH 2 ) 8 NH 3 + is long-ranged and directional (Type-IIa), with r(N···I) and ∠C-N···I = 3.830 Å and 176.9 • , respectively. The crystal structure of a 2D metal bromide perovskite system, catena-(tetrakis(4carboxybutylammonium) tetrakis(m-chloro)-tetrachloro-di-lead), [C 5 H 12 NO 2 + ] 4 (Pb 2 Cl 8 ] n , is shown in Figure 14b. The two mono-cationic linkers connecting the two neighboring inorganic layers are chain-like organic dications. Each dication, [(HOOC(CH 2 ) 4 NH 3 )] 2+ , is linked to another dication through two non-equivalent C-H···O hydrogen bonds (r(H···O) 1.909 and 1.950 Å) (see the central region marked yellow in Figure 14b). The -NH 3 + fragment of the linker is bonded to the Cl atoms of a face of a PbCl 6 octahedron and involves N-H···Cl hydrogen bonds and C-N···Cl pnictogen bonds. The two types of pnictogen bonds have r(N···Cl) of 3.461 and 3.298 Å and are both quasi-linear Type-IIa interactions, with ∠C-N···Cl = 173.7 • and 161.8 • , respectively. The -NH 3 + head is also involved in the formation of three or four non-equivalent r(H···Cl) hydrogen bonds. The shorter Hbond distances of three of them are at 2.582, 2.474 and 2.322 Å, with angles of approach of the electrophile of 138.6 • , 152.9 • and 175.7 • , respectively. The first follows a Type-IIb pattern and the latter two follow a Type-IIa pattern of bonding. Similarly, the H sites of the other ammonium fragment involved in the other pnictogen bond (r(N···Cl) = 3.298 Å) is also involved in forming hydrogen bonds with distances of 2.474, 2.327 and 2.322 Å, corresponding to the angles of approach of the electrophile of 152.9 • , 170.7 • and 175.7 • , respectively. Clearly, these attractive interactions interplay with each other to determine octahedral tilting of the PbCl 6 octahedra of the inorganic frameworks of the metal halide perovskite. The octahedral tilting is a factor influencing the band gap and effective masses of charge carriers, a necessary perquisite for the functioning of perovskite compounds in photovoltaics and optoelectronics [26,112], and is driven by the competition between the N and H sites of the ammonium head in forming pnictogen and hydrogen bonding interactions.
Diammonium cations with the generic formula + H 3 N-R-NH 3 + , where R group may be a linear-chain alkyl group, or an even-membered derivative of ethane, or a fused aromatic ring such as naphthalene, are well known. As spacers, they have been used as template moieties to stabilize the inorganic framework of a large number of inorganic-organic hybrid nanocomposites (cf. Figures 14c and 15). Many two-dimensional layered perovskites with a RbAlF 4 type structure and with linear alkyl chains between inorganic layers have been reported.  6 ], which have edge-sharing twin-anionic chains of corner-sharing ribbons, have also been reported.
The structure of (NH 3 (CH 2 ) 12 NH 3 ) n PbI 4 ·n(PbI 4 ) (Figure 14c) has been reported by two different groups [110,113]. . The pnictogen bond in the first two of these systems is far from linear. The r(N···Br) in (H 3 N(CH 2 ) 10 NH 3 )PbBr 4 is 3.597 Å, which is less than the sum of the vdW radii of Br and N, 3.52 Å. This is a notable feature of all the 2D iodide-based perovskites (cf. Figures 14 and 15).
In cases where the templating agent is an are aromatic moiety, as in (H 3 NC 10 H 6 NH 3 )PbI 4 , shown in Figure 15d, the C-N···I pnictogen bonds are somewhat longer and linear. Both the ammonium fragments form equivalent pnictogen bonds, with ∠C-N···I = 178.6 • and r(N···I) = 3.930 Å. The three shortest N-H···I hydrogen bonds identified for all cases are also illustrated in Figure 15d. They are formed by the H atoms of the ammonium fragment; two of them followed a Type-IIa pattern and the other a Type-IIb pattern.
A similar system is [(NH 3 (CH 2 ) 3 NH 2 (CH 2 ) 3 NH 3 ) 2 ][Cd 2 Cl 10 ], in which the organic cation carries a formal charge of +3 (Figure 16a). The 2D framework manifests through H···Cl hydrogen bonds and N···Cl pnictogen bonds; the former are formed by the terminal H atoms of the two ammonium groups and chloride ions. The latter are directional, with ∠C-N···Cl = 176.   [110]. Bond lengths and bond angles associated with selected intermolecular interactions are in Å and degrees, respectively. Space groups are shown for each case. Intermolecular interactions are marked as dotted lines in green.
The structure of (NH3(CH2)12NH3)n PbI4·n(PbI4) (Figure 14c) has been reported by two different groups [110,113]. Both    The crystal structures shown in Figure 16b-e are also two-dimensional polymeric structures, with different organic linkers. Whereas the H and N atoms of the ammonium fragment in these linkers are involved both in the formation of H···Cl hydrogen bonds and N···Cl pnictogen bonds, the chain-length of the organic spacers between the inorganic frameworks in Figure 16b,c,e is reinforced by attractive intermolecular interactions between a pair of two organic cations. The reinforcement is driven predominantly by chargeassisted C···O tetrel bonds and H···O hydrogen bonds between the two 4-(methoxyphenyl)methanaminium cations in Figure 16b, by C···F tetrel bonds between the 2,2,2-trifluoroethan-1-aminium cations in Figure 16c and by C···C/C···H interactions between the POHFEA [114] [118]). Intermolecular interactions are marked as dotted lines in cyan, the hanging contacts are marked as dotted lines in red. Selected bond distances and bond angles are in Å and degrees, respectively.
The crystal structures shown in Figure 16b-e are also two-dimensional polymeric structures, with different organic linkers. Whereas the H and N atoms of the ammonium fragment in these linkers are involved both in the formation of H···Cl hydrogen bonds and N···Cl pnictogen bonds, the chain-length of the organic spacers between the inorganic frameworks in Figure 16b,c,e is reinforced by attractive intermolecular interactions between a pair of two organic cations. The reinforcement is driven predominantly by charge-assisted C···O tetrel bonds and H···O hydrogen bonds between the two 4-(methoxyphenyl)methanaminium cations in Figure 16b, by C···F tetrel bonds between the 2,2,2-trifluoroethan-1-aminium cations in Figure 16c and by C···C/C···H interactions between the N-ethylethylenediammonium cations in Figure 16e. In the case of the crystal of (CH 3 CH 2 NH 2 CH 2 CH 2 NH 3 ) 2n ·n(CdCl 4 ), Figure 16e, the H site of the acetylenic fragment of the cation forms a strong H···Cl hydrogen bond with the inorganic framework and is strong enough to compete with the H···Cl hydrogen bonds and N···Cl pnictogen bonds formed by its other end. As noted above, perovskite-type layer compounds with general formula (R-NH 3 ) 2 [MX 4 ] or (NH 3 -(CH 2 ) n -NH 3 )[MX 4 ], (M = Cd, Mn, Cu, Pb . . . ; X = Cl, Br, I; and R = a saturated or an unsaturated chain) have received considerable attention for some time as their crystal structures exhibit a quasi-two-dimensional arrangement with alternating inorganic (MX 6 )-organic (alkyl or alkylene chains) layers and undergo numerous structural phase transitions under temperature and/or pressure effects [119][120][121][122][123][124][125][126][127][128][129][130][131].
Another set of 2D metal halide perovskite systems in which H···N hydrogen bonds and directional pnictogen bonds formed between the organic and inorganic octahedra frameworks reinforce the layered structure, but in which they are not solely responsible for the crystalline structure, is shown in Figure 17. The pnictogen bonds appear as N···I contacts in Figure 17a,b, and as N···Br contacts in Figure 17c,d. The N···I and N···Br contacts, respectively, in Figure 17a,c,d are longer than the N···I contacts in Figure 17b, which are exceptionally short, suggesting a significant degree of covalency [132][133][134][135]. Other than the two types of interactions formed by the ammonium head mentioned above, we found that part of the stability of the system is driven by the attractive interaction of iodine in I−(CH2)2−NH3 + and the inorganic framework. As shown in Figures 17a  and 18a, these are nothing other than charge assisted C-I···I halogen bonds and C-H···I hydrogen bonds. The first are stronger than the second, inferred from the intermolecular Other than the two types of interactions formed by the ammonium head mentioned above, we found that part of the stability of the system is driven by the attractive interaction of iodine in I−(CH 2 ) 2 −NH 3 + and the inorganic framework. As shown in Figures 17a and 18a, these are nothing other than charge assisted C-I···I halogen bonds and C-H···I hydrogen bonds. The first are stronger than the second, inferred from the intermolecular distances and directionality features. In a study of some related structures, viz. (Cl-(CH 2 ) 2 -NH 3 ) 2 PbI 4 ( Figure 18b) and (Br-(CH 2 ) 2 -NH 3 ) 2 PbI 4 (Figure 18c)) [140], it was suggested that the structural features of these systems are the result of halogen and hydrogen bonding at the organic−inorganic interface; the hydrogen bonds between equatorial iodine atoms of the perovskite layer and ammonium parts are absent.
Although the characterization of intermolecular bonding interactions in these systems was based purely on geometric factors, our close inspection not only revealed various other interactions, as mentioned above, but the conclusions of the authors may be incorrect since the equatorial iodine atoms of the perovskite layer and ammonium fragments are not free from hydrogen bonding (cf. Figure 18a-c). It was further argued that the red salts, such as (Cl−(CH 2 ) 2 −NH 3 ) 2 PbI 4 , display a reduced band gap (experimental: 2.2 eV; calculation: 1.9 eV) compared to (I−(CH 2 ) 2 −NH 3 ) 2 PbI 4 in which tilting of the PbI 6 octahedra and a larger bandgap (experimental: 2.45 eV; calculation: 2.3 eV) occurs [140]. This may not be correct since it is not a straightforward matter to determine whether the reduced band gap for the red-salt perovskites of the series has emerged from hydrogen bonding or halogen bonding interactions, or a combination of the two, at the organic−inorganic interface.
We did not find such directional N···I pnictogen bonds between the organic and inorganic counterparts in the undistorted lead iodide perovskite systems with X −(CH 2 ) 2 −NH 3 + (X = Br, Cl) as spacer cations. It could be due to the organic spacer adopting a cis configuration as a result of its small size and reactivity of the lighter halogens. The geometry of these cations is accompanied by a charge-assisted intramolecular N-H···X hydrogen bonding interaction (r(H···X ) 2.898 Å in Br−(CH 2 ) 2 −NH 3 + and 2.778 Å in Cl−(CH 2 ) 2 −NH 3 + ), thus providing additional stability to the structure (cf. Figure 18b,c, respectively). It should be appreciated that the (Br−(CH 2 ) 2 −NH 3 ) 2 PbI 4 and (Cl−(CH 2 ) 2 −NH 3 ) 2 PbI 4 crystals are not free from C-H···I and N-H···I hydrogen bonds, as well as charge-assisted X ···I (X = Cl, Br) halogen bonds. The marked differences in intermolecular bonding features, as well as the inflexible nature of the conformational degree of freedom of the three bifunctional organic cations (adopts to a trans configuration in Figure 18a and to a cis configuration in Figure 18b,c), between the three crystal systems may explain why the degree of distortion of the perovskite layers, in terms of the degree of tilting of the octahedral related to the mean equatorial PbI 2 plane, occurs in (I−(CH 2 ) 2 −NH 3 ) 2 PbI 4 , but not in (Cl−(CH 2 ) 2 −NH 3 ) 2 PbI 4 and (Br−(CH 2 ) 2 −NH 3 ) 2 PbI 4 . From these structural features, it is clear that the bandgap reduction in the red-salts is a joint effect of the conformational freedom of the organic cation, the extent of various intermolecular interaction and the tilting of the octahedra.
The N···I pnictogen bonds in the crystal [2(C 9  In the case of (PPA) 2 (MA 0.5 FA 0.5 ) 2 Pb 3 I 10 (Figure 20b), the pnictogen bond is also short but less directional (∠C-N···I = 167.7 • and (r(N···I) = 2.954 Å). The change in the geometric property of the pnictogen bond is not unexpected since the conjugated counter-cation spacer adjusts its position to maximize its attractive interaction with the iodides of the PbI 6 octahedra of the inorganic frameworks. Both these crystals also show -HC-H···I type hydrogen bonds within a bond length range of 2.8-3.3 Å, which are significantly non-linear (∠C-H···I = 110-145 • ); they are weaker than the N-H···I hydrogen bonds. The N···I pnictogen bonds in the crystal [2(C9H12N),0.5(C2H6),0.5(Pb4I14)] (referred t by the authors as (PPA)2(MA0.5FA0.5)Pb2I7 (PPA = 3-phenyl-2-propenammonium; MA methylammonium; FA = formamidinium [137]), Figures 17b and 19a, are short and dire tional, r(N···I) = 2.814 Å and C-N···I = 175.3°, and are accompanied by non-linear N-H·· hydrogen bonds (r(H···I ) = 2.6-2.7 Å ; N-H···I values in the range 89-120.0°).  eties has led to the assembly of more efficient solar cells. The structure consists of 2D inorganic layers of corner sharing [PbI6] 4− octahedra which are separated by two layers of organic phenylpropylene ammonium cations (PPA + ), and adjacent inorganic layers are shifted by ½ unit-cell length. The recorded trends in the absorption spectra suggested a progressively reduced bandgap as n increases because of less quantum and dielectric confinement, with optical bandgaps of approximately 1.94 and 1.75 eV for (PPA)2(FA0.5MA0.5)n−1PbnI3n+1 when n = 2 and 3, respectively.  Polyhedral and ball-and-stick models of the crystal structure of (PPA) 2 (FA 0.5 MA 0.5 ) n−1 Pb n I 3n+1 materials in 2D/3D: (a) n = 2 and (b) n = 3. Neither of these structures, as deposited in the CSD, show the correct orientation of the organic cations within the 3D perovskite cages, as provided in the supplementary material to the study [137], and hence we replaced them here by blue ellipses. The pnictogen bonded contacts are shown as dotted lines in cyan and the hanging contacts are in red. Selected intermolecular N···I, C-H···I and N-H···I bond distances and bond angles in Å and degrees, respectively.
Both crystals referred to above adopt a Ruddlesden-Popper (RP) type structure characterized by the general formula (A ) 2 (A) n−1 Pb n X 3n+1 (A = spacer cation, A = cage cation, and X = halide anion). In the reported crystal systems, 3-phenyl-2-propenammonium (PPA) with conjugated backbone plays the role of the spacer cation between the inorganic counterparts, and FA and MA as the inorganic cage cations (shown as blue ellipses in Figure 20), at low-temperature. This arrangement between the organic and inorganic moieties has led to the assembly of more efficient solar cells. The structure consists of 2D inorganic layers of corner sharing [PbI 6 ] 4− octahedra which are separated by two layers of organic phenylpropylene ammonium cations (PPA + ), and adjacent inorganic layers are shifted by 1 2 unit-cell length. The recorded trends in the absorption spectra suggested a progressively reduced bandgap as n increases because of less quantum and dielectric confinement, with optical bandgaps of approximately 1.94 and 1.75 eV for (PPA) 2 (FA 0.5 MA 0.5 ) n−1 Pb n I 3n+1 when n = 2 and 3, respectively.
The yellow plate-like crystals (BA) 4 AgBiBr 8 and (BA) 2 CsAgBiBr 7 (BA = n-butylammonium = CH 3 (CH 2 ) 3 NH 3 + ) are also members of the family A 2 A n−1 B n X 3n+1 (n = 1 and 2). The structure of these organic-inorganic halide (001) perovskites were reported at 100 and 298 K [141]. However, the atomic positions of the organic cation, determined at the higher temperature, overlapped; hence this was not considered in this study. Although the overall geometry of each crystal lattices determined at the two different temperatures look similar, the local bonding environment is very different (Figure 21a,b). In the low temperature phase, the ammonium head of the organic cation is involved both in the formation of pnictogen bonding and hydrogen bonding, whereas the methyl head is involved in tetrel bonding. The pnictogen bond follows a Type-IIa pattern.
hanging contacts are in red. Selected intermolecular N···I, C-H···I and N-H···I bond distances and bond angles in Å and degrees, respectively.
The yellow plate-like crystals (BA)4AgBiBr8 and (BA)2CsAgBiBr7 (BA = n-butylammonium = CH3(CH2)3NH3 + ) are also members of the family A2An−1BnX3n+1 (n = 1 and 2). The structure of these organic-inorganic halide (001) perovskites were reported at 100 and 298 K [141]. However, the atomic positions of the organic cation, determined at the higher temperature, overlapped; hence this was not considered in this study. Although the overall geometry of each crystal lattices determined at the two different temperatures look similar, the local bonding environment is very different (Figure 21a,b). In the low temperature phase, the ammonium head of the organic cation is involved both in the formation of pnictogen bonding and hydrogen bonding, whereas the methyl head is involved in tetrel bonding. The pnictogen bond follows a Type-IIa pattern. Shown on the right of (a,b) are the illustration of local modes of various interlayer intermolecular interactions (viz. hydrogen bonds, tetral bonds, and pnictogen bonds) formed between the inorganic and organic frameworks. Selected bond distances and bond angles are in Å and degrees, respectively. Some of the organic cations in (b) have been deleted for clarity. The CSD references are shown in each case in uppercase letters. The purple and white-gray octahedra correspond to BiBr6 and AgBr6 octahedra, respectively, whereas the dark-purple balls inside the cage in (b) represent the cesium cations.   Figure 21b (right). The shortest C-H···Br hydrogen bonds formed by the two methyl groups of the organic cation chain are longer and not equivalent (see Figure 21b, right).
In the case of (BA) 4 AgBiBr 8 , all three types of non-covalent links persist, yet the N···Br pnictogen bonds are equivalent and are shorter than that found in (BA) 2 CsAgBiBr 7 . By contrast, the C-C···Br tetrel bonds formed by the terminal methyl carbon of the organic cation repeat in an alternate fashion in (BA) 4 (Figure 21a, right); however, they are equiv-alent and linear in (BA) 2 CsAgBiBr 7 , with r(C···Br) = 3.992 Å (and ∠C-C···Br = 179.1 • ) (Figure 21b, right). Clearly, the competition between N···Br pnictogen bonding, C···Br carbon bonding, and N-H···Br and -(H)C-H···Br hydrogen bonding interactions in the crystal causes tilting of the corner-shared PbB 6 /AgBr 6 octahedra in both the crystals. Because of the observed differences in the intermolecular bonding properties, the two perovskite systems in their powder and thin film forms display an onset of optical absorption at energies close to 2.7 eV ((BA) 2 CsAgBiBr 7 ) and 2.77 eV ((BA) 4 AgBiBr 8 ), highlighting the semiconducting nature of these systems in the solid state [141].
A very similar multiaxial layered halide double perovskite system, determined at 200, 298 and 400 K, was recently reported [139] (Figure 17d). It exhibits a fascinating ferroelectric feature with a Curie temperature of~273 K and unique semiconducting nature. In the low temperature phase of the system the N···Br pnictogen bond is more directional (r(N···Br) = 3.514 Å; ∠C-N···Br = 169.9 • ) and is accompanied by N-H···Br and -(H)C-H···Br hydrogen bonding interactions. As above, the [(CH 3 ) 2 CHCH 2 NH 3 ] 2 PbCl 4 lead halide 2D perovskite, Figure 19a, features an exceptionally large electrocaloric effect near room temperature [142]. It also exhibits a sharp first-order phase transition at 302 K, a superior spontaneous polarization > 4.8 µC cm −2 and relatively small coercive field (<15 kV cm −1 ). The low temperature monoclinic phase (285 K) of the crystal is a result of a marriage between the [(CH 3 ) 2 CHCH 2 NH 3 ] 2 and PbCl 4 ions, accompanied by hydrogen and pnictogen bonding interactions. The latter interactions, N···Cl, are shown in Figure 19a. They follow a quasi-linear Type-IIa bonding topology, and the N-H···Cl hydrogen bonds are in the range 2.35-2.90 Å.
McClure et al. [143] have reported the synthesis, structure, and bandgaps of four layered halide double perovskites, BA 2 Cu 0.5 In 0.5 Cl 4 , BA 2 Ag 0.5 In 0.5 Cl 4 , BA 2 Ag 0.5 Sb 0.5 Cl 4 and BA 2 Ag 0.5 Sb 0.5 Br 4 (BA = CH 3 (CH 2 ) 3 NH 3 + ). They have an RP structure with n = 1. They revisited the crystal structure of BA 2 Ag 0.5 Bi 0.5 Br 4 (Figure 19b) reported by Corner et al. [141] and also reported the crystal structure of BA 2 PbCl 4 . All these systems exhibited bandgaps in the range 2.65 eV (BA 2 Ag 0.5 Bi 0.5 Br 4 ) to 4.27 eV (for BA 2 Ag 0.5 In 0.5 Cl 4 ), with the bromide compounds featuring smaller bandgaps than the chloride compounds; this is expected since halide perovskites with heavier halogens generally exhibit smaller bandgaps [84]. The bandgaps of the bromide-containing RP phases, BA 2 Ag 0.5 M 0.5 Br 4 (M = Bi 3+ and Sb 3+ ), fall in the visible range, whereas those with chloride phases are wide bandgap materials (bandgap values between 3.2 and 4.3 eV); the former are shown to be suitable for visible light absorption. Similarly, the bandgaps of layered structures with composition BA 2 M 0.5 M 0.5 X 4 were observed to be between 0.5-0.8 eV larger than the analogous 3D Cs 2 MM X 6 cubic double perovskites, a consequence of dimensional reduction (3D → 2D), accompanied by distortions of the octahedral environment around the M/M ions, and octahedral tilting. In the 2D systems, the octahedral tilting is driven by various non-covalent interactions between the inorganic and organic layers, including pnictogen bonding, as noted in several crystal systems discussed above. For example, BA 2 Ag 0.5 Bi 0.5 Br 4 , Figure 19b, has equivalent N···Br pnictogen bonds, with r(N···Br) = 3.522 Å and ∠C-N···Br = 169.4 • . They are directional, and together with hydrogen bonds, are clearly responsible for the rotation and tilting of the SbBr 6 and AgBr 6 octahedra in the structure.
Mei et al. [144] were probably the first to show that the protonated 5-amino valeric acid cation, HOOC(CH 2 ) 4 NH 3 + , (5-AVA), can be used as an organic additive for the development of long-term stability of 2D perovskite solar cells. The cation has the advantage that it is not as toxic as organic diammonium cations, which are potentially deleterious to human health [109]. Subsequently, many such systems were reported, including, for instance, the quaternary representative 5-AVA 2 MA n-1 Pb n Br 3n-1 series of 2D bromide perovskite compounds [145] where n = 1-3. The value of n determines whether the compound comprises one, two or three layers. Similar 2D perovskite compounds containing chloride and iodide have also appeared. In the case of the iodide derivatives, n can reach values up to 6 or 7 [146]. However, the authors of this study were unable to synthesize pure samples of 5-AVA 2 MAPb 2 Br 7 and 5-AVA 2 MA 2 Pb 3 Br 10 [145]. The crystal structures of these 2D compounds demonstrate that they share some aspects of the structure of the cubic 3D perovskites with layers of corner-sharing PbBr 6 octahedrons in (100) orientation. The structural stability of these crystals, such as other perovskite systems discussed above, arises from a significant number of hydrogen bonding interactions between the inorganic and organic components of the structure.
The PbX 6 octahedra are rotated in these systems, and of course, the extent of this depends on the nature of halides in the inorganic frameworks and the intermolecular interactions involved. There is a monotonic decrease of the band gap from 3.02 eV through 2.80 eV to 2.70 eV as n increases from 1 to 2 and 3 in 5-AVA 2 MA n-1 Pb n Br 3n-1 . The 5-AVA 2 PbBr 4 perovskite system also displayed strong blue fluorescence (413 nm) with a remarkably small Stokes shift (47 nm). However, the stability of the systems with n = 2 and n = 3 was reduced compared to the simpler systems 5-AVA 2 PbBr 4 and MAPbBr 3 -a result that was not in line with what was reported for the corresponding systems featuring iodide in the inorganic layer. This may be a consequence of the intermolecular chemical bonding between the organic and inorganic moieties and merits a detailed study.
The crystal structures of [HOOC(CH 2 ) 4 NH 3 ] 2 PbBr 4 and [HOOC(CH 2 ) 4 NH 3 ] 2 PbCl 4 , showing the local nature of hydrogen-and pnictogen bonding interactions formed by the halides and the ammonium head of the cation, as well as between the halides and the two methyl groups the organic cation, are shown in Figure 22 [147]. The N-H···Br hydrogen bonds formed by ammonium moiety are significantly stronger than the -(H)C-H···Br hydrogen bonds formed by the methyl moieties ( Figure 22a). There is a network of pnictogen bonds in the crystal that appear in two varieties, with r(N···Br) of 3.602 and 3.498 Å; the first is less directional than the second (∠C-N···Br = 156.4 • vs. 161.3 • , Figure 22b). The large deviation from linearity is because there is competition between all three H atoms of the ammonium moiety for formation of hydrogen bonds of appreciable strengths (based on the H···Br bond distances). The structure of [HOOC(CH 2 ) 4 NH 3 ] 2 PbCl 4 features virtually similar intermolecular bonding at the interface region between the halides and the ammonium head (Figure 22c,d). The N-H···Cl and C-H···Cl hydrogen bonding and N···Cl pnictogen bonding interactions reported by Hillebrecht [148] are slightly different to those reported by Yang and coworkers [109]. Nevertheless, the overall nature of the local topology of non-covalent interactions is very such similar.
It was argued [109] that the intermolecular hydrogen bonding in AVA 2 PbCl 4 causes a larger cation penetration depth that enables larger structural deformation than that observed in the reference system PDAPbCl 4 (PDA = pentamethylenediammonium). The deformation of the PbCl 6 octahedra leads to ultrabroadband emission in AVA 2 PbCl 4 , which has an enhanced photoluminescence quantum yield (2.83%) compared to PDAPbCl 4 (0.4%). This led to the suggestion that the role played by intermolecular interactions results in the modulation of the photophysical properties of 2D perovskites, and this will be beneficial for the design of green perovskites for optoelectronic applications.  [109] and VUHVIG [148]) showing the nature of N-H···Cl and C-H···Cl hydrogen bonding and N···Cl pnictogen bonding interactions formed at the organic and inorganic interface. Selected bond distances and bond angles are in Å and degrees, respectively. It was argued [109] that the intermolecular hydrogen bonding in AVA2PbCl4 causes a larger cation penetration depth that enables larger structural deformation than that observed in the reference system PDAPbCl4 (PDA = pentamethylenediammonium). The deformation of the PbCl6 octahedra leads to ultrabroadband emission in AVA2PbCl4, which has an enhanced photoluminescence quantum yield (2.83%) compared to PDAPbCl4 (0.4%). This led to the suggestion that the role played by intermolecular interactions results in the modulation of the photophysical properties of 2D perovskites, and this will be beneficial for the design of green perovskites for optoelectronic applications.

Pnictogen Bonding in 1D (One-Dimensional) Perovskite Systems
Many low-dimensional halide-based organic-inorganic hybrid perovskite compounds are known [149][150][151][152]. They sometimes adopt a step-like (SL-type) structure which is effectively a 1D quantum wire with chains of corner-sharing octahedra "insulated" by blocks of face-sharing octahedra [150]. For instance, among others, Hoffman et al. [150] have reported a new series of 1D SL-type structures with the general formula (PA)2m+4(MA)m−2Pb2m+1I7m+4 (m = 2, 3, 4), where the PbI6 octahedra connect in a corner-and face-sharing motif that exhibit resistivity trends that are dominated by ionic transport and  [109] and VUHVIG [148]) showing the nature of N-H···Cl and C-H···Cl hydrogen bonding and N···Cl pnictogen bonding interactions formed at the organic and inorganic interface. Selected bond distances and bond angles are in Å and degrees, respectively.

Pnictogen Bonding in 1D (One-Dimensional) Perovskite Systems
Many low-dimensional halide-based organic-inorganic hybrid perovskite compounds are known [149][150][151][152]. They sometimes adopt a step-like (SL-type) structure which is effectively a 1D quantum wire with chains of corner-sharing octahedra "insulated" by blocks of face-sharing octahedra [150]. For instance, among others, Hoffman et al. [150] have reported a new series of 1D SL-type structures with the general formula (PA) 2m+4 (MA) m−2 Pb 2m+1 I 7m+4 (m = 2, 3, 4), where the PbI 6 octahedra connect in a corner-and face-sharing motif that exhibit resistivity trends that are dominated by ionic transport and no photoresponse. Yuan and coworkers [153] have reported the synthesis, crystal structure and photophysical properties of a 1D organic lead bromide perovskite, C 4 N 2 H 14 PbBr 4 . It consists of edge-sharing octahedral lead bromide chains [PbBr 4 2− ] n that are surrounded by the C 4 N 2 H 14 2+ organic cations to form the bulk assembly of core-shell quantum wires. The 1D structure has enabled the authors to demonstrate the presence of strong quantum confinement with the formation of self-trapped excited states that give efficient bluish white-light emissions with photoluminescence quantum efficiencies of approximately 20% for the bulk single crystals and 12% for the microscale crystals. Shown in Figure 23 are a few illustrative 1D crystal structures that are stabilized by the joint involvement of hydrogen-and pnictogen-bonding interactions between the organic and inorganic frameworks. temperature, Figure 23a [154]. The photoluminescence measurements showed a strong emission line at 3.64eV. The unaided-eye detectable blue luminescence emission came from the excitonic transition in the Sb2Br10 anions. The stability of the crystal structure arises from the intermolecular interaction between the bioctahedral [Sb2Br10] 4− dimers composed of two equivalent distorted octahedrons sharing one edge and the organic N,Ndiethylethylendiammonium cation. Although the N-H⋯Br hydrogen bonds were shown to be the sole synthon for holding the crystal network, we found that C-N⋯Br pnictogen bonds in the crystal should not be overlooked. The bond distance associated with the latter are less than the sum of the vdW radii of Br and N and are directional (r(N···Br) = 3.400 Å ; C-N···Br = 176.0°). The monoclinic hybrid compound [C 6 H 18 N 2 ] 2 Sb 2 Br 10 is monoclinically stable at room temperature, Figure 23a [154]. The photoluminescence measurements showed a strong emission line at 3.64 eV. The unaided-eye detectable blue luminescence emission came from the excitonic transition in the Sb 2 Br 10 anions. The stability of the crystal structure arises from the intermolecular interaction between the bioctahedral [Sb 2 Br 10 ] 4− dimers composed of two equivalent distorted octahedrons sharing one edge and the organic N,Ndiethylethylendiammonium cation. Although the N-H· · · Br hydrogen bonds were shown to be the sole synthon for holding the crystal network, we found that C-N· · · Br pnictogen bonds in the crystal should not be overlooked. The bond distance associated with the latter are less than the sum of the vdW radii of Br and N and are directional (r(N···Br) = 3.400 Å; ∠C-N···Br = 176.0 • ).
The crystal structure of CH 3 NH 2 -CH 3 NH 3 PbI 3 , Figure 23b, was also considered to involve H···N hydrogen bonds between the organic moieties CH 3 NH 2 and CH 3 NH 3 + ) [155]. The CH 3 NH 2 ···CH 3 NH 3 + dimers in the CH 3 NH 2 -CH 3 NH 3 PbI 3 intermediates were accompanied by 1D-PbI 3 − chains (δ-phase). The weakly hydrogen-bonded CH 3 NH 2 molecules could be easily released from the CH 3 NH 2 -CH 3 NH 3 PbI 3 intermediates, probably causing rapid, spontaneous phase transition from 1D-PbI 3 − (δ-phase) to 3D-PbI 3 − (α-phase). From the crystal structure shown in Figure 23b we found that, in addition to H···N hydrogen bonds formed by the methylamine and methylammonium moieties, there are also weak C-N· · · I pnictogen bonds in the crystal that are directional. The weak aspect of these bonds is evident from the N···I intermolecular distances (r(N···I) = 3.655 Å; ∠C-N···I = 176.5 • ) that are less than the vdW radii sum of N and I atoms, 3.70 Å, and are longer than N···Br pnictogen bonds found in the [C 6 H 18 N 2 ] 2 Sb 2 Br 10 crystal (Figure 23a).
Senior and coworkers [157] have reported the crystal structures of six halobismuth(III) salts of a variety of substituted aminopyridinium cations. These crystals feature discrete mononuclear [BiCl 6 ] 3− and binuclear [Bi 2 X 10 ] 4− anions (X = Cl or Br), and polymeric cis-double-halo-bridged [Bi n X 4n ] n− anionic chains (X = Br or I). One of these, bis(4-amino-3-ammoniopyridinium) diµ-chlorido-bis[tetrachloridobismuth(III)] dihydrate, (C 5 H 9 N 3 ) 2 [Bi 2 Cl 10 ]·2H 2 O, incorporates discrete [Bi 2 Cl 10 ] 4− anions. The Bi(III) centers in this and other crystals adopted a slightly distorted octahedral geometry and there was a correlation between the Bi-X bond lengths and the number of classic N-H···X hydrogen bonds that the X ligand accepts, with a greater number of interactions corresponding with slightly longer Bi-X distances. The supramolecular networks formed by classic N-H···X hydrogen bonds include ladders, bilayers and three-dimensional frameworks. As shown in Figure 23c, the geometric stability of (C 5 H 9 N 3 ) 2 [Bi 2 Cl 10 ]·2H 2 O not only emerges from the hydrogen bonds, but also C-N···Cl pnictogen bonds that are somewhat stronger than the C-N···Br and C-N···I found in the crystals shown Figure 23a,b, respectively.
The structure of catena-[methylammonium bis(µ-iodo)-iodo-lead monohydrate (CH 6 N + )(I 3 Pb − ) n ,(H 2 O) [158] is shown in Figure 23d. The crystal is a needle-shaped solvation intermediate similar to CH 3 NH3PbI 3 ·DMF that has been recognized as the main cause for the incomplete coverage of the resultant thin film. By avoiding these intermediates, the films crystallized at the gas-solid interface offer several beneficial features for device performance including high surface coverage, small surface roughness, as well as controllable grain size. The stability of the intermediate is organized by N-H···I and N-H···O hydrogen bonds, and directional C-N···I pnictogen bonds.
Another illustrative 1D perovskite system in 1D, [(C 8 H 12 N) 2 (PbI 4 )] n , that crystallizes as an organic-inorganic hybrid is shown in Figure 24 [159]. The structure consists of chains of [PbI 2 ] − units extending along the crystallographic b axis. In the asymmetric unit Pb is on a twofold axis and the two I atoms are in general positions. Each Pb atom is octahedrally coordinated to six I atoms, arranged as chains of edge-sharing octahedra. The N-H···I hydrogen-bonding interactions between the three H atoms of the ammonium group of the organic cation and the metal coordinate I atoms (r(H···I) 2.566, 2.835 and 3.147 Å); these are stronger than the C-H···I formed by a methyl group of the same cation and the coordinate I atoms (viz. r(H···I) 3.312 and 3.393 Å). There are also C-H···π and N-H···π hydrogen bonds between the ammonium and methyl groups of one organic cation and the arene ring system of another cation, as well as the formation of arene C-H···I hydrogen bonds and weakly bound Type-I I···I close contacts. Apart from these, the N site of the ammonium group is involved in forming pnictogen bonds with the nearest I sites of the inorganic framework, Figure 24  inorganic framework, Figure 24 (right). These pnictogen bonds are quasi-linear, (r(N···I) = 3.634 Å ; C-N···I = 177.7°. Figure 24. Ball-and-stick and polyhedral model of the 2 × 2 × 2 cell of the crystal of catena-(bis dimethylanilinium) bis(μ2-iodo)-diiodo-lead(II)) [(C8H12N + )2n, (I4Pb 2− )] [159]. Intermolecular con between bonded atomic basins are illustrated as dotted lines in cyan and hanging contacts in Selected bond distance and bond angle associated with the nitrogen-centered pnictogen bond Å and degrees, respectively. The CSD reference is shown in uppercase letters.

Pnictogen Bonding in Zero-Dimensional (0D) Crystal Systems
Several zero-dimensional (0D) crystal systems have been reported, wherein pn gen bonding seemingly plays an important role in determining the stability of the crys   [159]. Intermolecular contacts between bonded atomic basins are illustrated as dotted lines in cyan and hanging contacts in red. Selected bond distance and bond angle associated with the nitrogen-centered pnictogen bond is in Å and degrees, respectively. The CSD reference is shown in uppercase letters.
Intermolecular contacts between bonded atomic basins are illustrated as dotted lines in cyan and hanging contacts in red. Selected bond distances and bond angles associated with the nitrogen-centered pnictogen bonds are in Å and degrees, respectively.
(C 6 H 7 NCl) 3 (BiCl 6 )·H 2 O, Figure 25c [ 165], is a non-centrosymmetric organic-inorganic metallate complex and synthesized at room temperature. The cohesion between the organic and inorganic entities was suggested to be driven by N-H···Cl, N-H···O, C-H···O, C-H···Cl and O-H···Cl hydrogen bonding interactions between the 4-dichloroanilinium cations, the [BiCl 6 ] 3− anions and water molecules, thus forming a pseudo-3D network. Our investigation of the crystal, however, revealed that pnictogen bonds also contribute to the cohesion of the system; this was overlooked in the original study. The pnictogen bonds appear at different intermolecular distances, Figure 25c; one of them has r(N···Cl) = 3.231 Å and ∠C-N···Cl = 175.6 • , and the other has r(N···Cl) = 3.214 Å and ∠C-N···Cl = 179.3 • .
The orthorhombic organic-inorganic compound (C 2 H 10 N 2 2+ ) 2 ,(MnBr 2− )(Br − ) 2 , Figure 25d [166], is a 0D structure. The authors examined intermolecular interactions using single crystal X-ray diffraction and Hirshfeld surface analyses. It was found that this system exhibits an optical bandgap of 2.61 eV, a strong photoluminescence excitation line at 362 nm and a strong fluorescence at 537 nm and that the hybrid material features good thermal stability (decomposition temperature: 288-660 • C). Their structural analysis suggested that the inorganic anions are stabilized by [NH 3 CH 2 CH 2 NH 3 ] 2+ dications via an extensive number of H···Br hydrogen bonds. However, we also found the presence of N···Br pnictogen bonds in the crystal. These are formed between N of ammonium moiety and Br of four-coordinate MnBr 4 2− ; they are quasi-linear (∠C-N···Br = 174.8 • ) and r(N···Br) = 3.386 Å). Shown in Figure 26a,b are the crystal structures of bis(1,4-phenylenediammonium) hexachloro-lead(II) [2(C 12 H 14 N 2 2+ ) 2 (PbCl 6 4− ) reported by two different groups [167,168]. Whilst the nature of intermolecular interactions between the interacting organic and inorganic moieties are virtually similar, they differ in the intermolecular bond distances and bond angles. The PbCl 6 4− octahedra are not linked to each other; this is therefore not a 1D crystal structure. The three H atoms of the ammonium head of the diammonium cation are involved in four N-H···Cl hydrogen bonds with bond distances of 2.278, 2.255, 2.356 and 2.883 Å (Figure 26a) and 2.282, 2.315, 2.547 and 2.885 Å (Figure 26b). The first three short bonds are Type-II interactions and the fourth is bent. The H atoms of the arene moiety placed close to the PbCl 6 4− octahedra are also hydrogen bonded, with C-H···Cl in the range 2.8-3.0 Å. Apart from these, the N site in the ammonium heads of the diammonium cation are pnictogen bonded, with r(N···Cl) = 3.550 Å in Figure 26a [169]; (d) 1,3-propanediammonium tetraiodo-cadmium(II) dihydrate (C3H12N2 2+ )(CdI4 2− )·2(H2O) [170]. Intermolecular contacts between bonded atomic basins are illustrated as dotted lines in cyan and hanging contacts in red. Selected bond distances and bond angles associated with the nitrogen-centered pnictogen bonds are in Å and degrees, respectively.

Statistical Analysis of N-Centered Pnictogen Bond Distances and Angles in Crystals
Our search of the CSD for an R-(H/D) 3 N···X-R' geometric motif with N···X in the range 2.6-4.0 Å and ∠R-N···X in the range 150-180 • (X = any halogen derivative; R = any element; R' = any metal) produced 1452 hits with 2337 close contacts. When the angular range was changed to 140- [171]), but the N site in ammonia in this crystal system probably does not possess a positive region on its electrostatic surface with which to engage with the halogen (Cl) in a neighboring building block, causing an apparent N···Cl interaction. We did not categorize N···X (X = halogen) close contacts in the other crystal systems as false contacts since N in them was positive. For example, the crystal [C 4 H 16  The histogram plots, together with the normal distribution of r(N···X) and ∠R-N···X found in the remaining 81 hits, comprising of 89 close contacts, are shown in Figure 27a,b. Figure 27c,d are from 1690 crystals with 3122 close contacts that emerged when N···X was in the range 2.6-4.0 Å and ∠R-N···X was in the range 140-180 • . As can be inferred from the former, the formation of N-centered pnictogen bonding in crystals does not occur frequently in the bond distance and bond angle ranges of 2.6-3.7 Å and 175-180 • , respectively. This suggests that highly directional nitrogen bonds are rare. This is understandable given that H atoms of the ammonium head are frequently involved in hydrogen bonds with nearestneighbor halogens. From Figure 27c,d, it is clear that nitrogen bonds occur frequently in the distance range 3.2-3.8 Å, and angular range between 145-175 • ; these were originated especially when the electron density donors are the heavier halogens, Cl, Br and I. We only analyzed 82 hits from one of these searches to determine whether there were any false hits; only one false close contact was found. This was an N···Cl contact that appeared genuine in [Ru(NO)(NH3)4(SO4)](HSO4)·H2O (r(N···Cl) and [R-N···Cl] = 3.307 Å and [175.8°]) (CSD ref. AXOQAH [171]), but the N site in ammonia in this crystal system probably does not possess a positive region on its electrostatic surface with which to engage with the halogen (Cl) in a neighboring building block, causing an apparent N···Cl interaction. We did not categorize N···X (X = halogen) close contacts in the other crystal systems as false contacts since N in them was positive. For example, the crystal The histogram plots, together with the normal distribution of r(N···X) and R-N···X found in the remaining 81 hits, comprising of 89 close contacts, are shown in Figure 27a,b. Figure 27c,d are from 1690 crystals with 3122 close contacts that emerged when N···X was in the range 2.6-4.0 Å and R-N···X was in the range 140-180°. As can be inferred from the former, the formation of N-centered pnictogen bonding in crystals does not occur frequently in the bond distance and bond angle ranges of 2.6-3.7 Å and 175-180°, respectively. This suggests that highly directional nitrogen bonds are rare. This is understandable given that H atoms of the ammonium head are frequently involved in hydrogen bonds with nearest-neighbor halogens. From Figure 27c,d, it is clear that nitrogen bonds occur frequently in the distance range 3.2-3.8 Å , and angular range between 145-175°; these were originated especially when the electron density donors are the heavier halogens, Cl, Br and I.  We have also made an individual analysis of the N···X geometry for X in crystals obtained from a CSD search. We used the same geometry, R-(H/D) 3 N···X-R', in our search, but the angle ∠R-N···X (X = I) was varied in the ranges 140 [174], was found to be a false hit since N in [C 5 H 10 I 2 N 4 OPt], does not have a positive region on its electrostatic surface. Figure 28a,b shows the normal distribution of r(N···I) and ∠R-N···I found in the remaining 81 hits comprising 89 geometric instances. From the former, it is apparent that N···I close contacts are extremely rare in the 2.6-3.2 Å range, and tend to occur in the range from 3.3 to 4.0 Å. The maximum of the normal distribution peaks around 3.77 Å. Similarly, the angle ∠R-N···I is populated largely around 160 • , because the occurrence of pnictogen bond is largely influenced by hydrogen bonds formed by the ammonium head. There are a small number of N···I close contacts in the 175-180 • range. We have also made an individual analysis of the N···X geometry for X in crystals obtained from a CSD search. We used the same geometry, R-(H/D)3N···X-R', in our search, but the angle R-N···X (X = I) was varied in the ranges 140-180° (150-180°) [ [174], was found to be a false hit since N in [C5H10I2N4OPt], does not have a positive region on its electrostatic surface. Figure 28a,b shows the normal distribution of r(N···I) and R-N···I found in the remaining 81 hits comprising 89 geometric instances. From the former, it is apparent that N···I close contacts are extremely rare in the 2.6-3.2 Å range, and tend to occur in the range from 3.3 to 4.0 Å . The maximum of the normal distribution peaks around 3.77 Å . Similarly, the angle R-N···I is populated largely around 160°, because the occurrence of pnictogen bond is largely influenced by hydrogen bonds formed by the ammonium head. There are a small number of N···I close contacts in the 175-180° range. Our search of the CSD for the geometry R-(H/D) 3 N···X-R' with ∠R-N···X (X = Br, Cl) in the range 140-180 • and N···X in the range 2.6-4.0 Å gave 653 hits (691 close contacts) for N···Br, and 776 hits (1355 close contacts) for N···Cl. When we examined some of these crystals, we found that several of the N···X were bent (∠R-N···X between 140 and 150 • ). We then decided to search for the N···X close contacts that are within the angular range 150-180 • . The resulting geometric instances are plotted in Figure 28c,d for N···Br close contacts, and in Figure 28e,f for N···Cl close contacts. The former consists of 554 geometric instances from 309 hits, and the latter of 1000 geometric instances from 658 hits. As can be seen, the nature of nature of the normal distribution of the angle ∠R-N···X is very similar when X = Br and X = Cl. A significant population occurs between 155 • and 170 • both for X = Cl and X = Br. However, the normal distribution curve associated with X = Cl shows that N···Cl close contacts are largely in the range 3.2-3.6 Å with a peak around 3.4 Å. This was not the case when X = Br, where the N···Br close contacts are largely in the 3.4-3.8 Å with a peak around 3.5 Å. These results suggest that N···X intermolecular distance can be longer than the sum of the vdW radii of N and X in many crystals.

Conclusions
This overview focused on several crystal structures, formed between various organic and inorganic moieties, to reveal the frequent existence of nitrogen-centered charge-assisted pnictogen bonding interactions; these have been largely overlooked in the metal halide perovskite literature. Since our motivation was to draw attention to the occurrence of this type of interaction in crystals, we examined its chemical manifestations. Our review of crystals structures deposited in the CSD suggests that N in ammonium derivatives is a potential pnictogen bond donor, and is, in part at least, responsible for the stability and functionality of many organic-inorganic hybrid perovskite systems, and hence could be exploited in the design of novel materials for photovoltaics, photocatalysis and other optoelectronics.
Our analysis suggests that the ammonium parts of mono-, bi-and multi-functional organic cations interact attractively with halogen atoms of the inorganic counterparts of the metal halide perovskites in two ways. Firstly, the N atom forms directional pnictogen bonding interactions with the bridging halogen sites of the MX 6 n− octahedra of the inorganic counterpart. Secondly, the H atoms form hydrogen bonding interactions either with the bridging halogens and/or terminal halogen atoms, depending on the architectures of the organic cations and inorganic anions. Based on the spatial orientation of the organic cation, an additive, a single H-atom in the ammonium fragment could sustain an attractive engagement with at least two halogen atoms of the inorganic octahedra, thus forming a bior tri-furcated interaction bonding topology.
In cases where one of the terminal sites of the organic cation is -H, -X, or -CH 3 , the formation of charge-assisted H···X hydrogen bonding, or X ···X halogen bonding, or H···X/C···X hydrogen/tetrel bonding interactions between the organic cation and inorganic metal halide anion framework is a likely outcome. These, together with the attractive engagements driven by the ammonium part of the organic cation, act as forces that hold the layered inorganic framework together in low dimension and hence provide stability to the overall geometry of the organic-inorganic metal halide compounds.
The tilting of the MX 6 octahedra of the inorganic framework in halide-based perovskites is often found to be driven by the competition between N···X pnictogen bonding, H···X hydrogen bonding and together with, when they occur, other non-covalent interactions. This attribute is shown through several illustrative crystal systems, including, for instance, orthorhombic MAPbX 3 (X = Cl, Br, I) systems in 3D (Figure 4), and those shown in Figures 14-16, Figures 18-20 in 2D. What specific role does the nitrogen-centered pnictogen bonds play in determining the optoelectronic properties of the functional perovskite materials is not revealed in this study, and hence should be the topic of any future studies. We believe that the nitrogen-centered pnictogen bonding interactions revealed in this study will be beneficial to those are interested in the development of novel organic-inorganic hybrid halide-based perovskite materials and beyond.
While the complexation energy of several illustrative simplified ion-pairs examined in this study was very large owning to the dominance of coulombic interactions involved, the empirical relationship known between potential energy density at the bcp and interaction energy allowed us to suggest that pnictogen bonds explored in these models were locally of weak-to-medium strength, and hence their contribution to the crystal stability should be not overlooked. However, it should be kept in mind that the intermolecular interaction energies estimated from intermolecular bond critical point properties of electron densities could be inherently unreliable, since they typically underestimate, but sometimes overestimate, more reliable values. Conclusions drawn regarding the energetic importance of specific intermolecular interactions based on this empirical relationship can be misleading, especially in the context of crystal packing and crystal engineering [175]. These interactions appeared persistently regardless of whether the geometry of the ion-pair in the crystal, or optimized with MP2 and DFT, was used for QTAIM analysis; their feasibility in the ion-pair was also confirmed by NBO's second order analysis. Although it is not straightforward to quantify the energetic strength of each of the interaction types in each crystal system illustrated in this study, even with detailed computational studies, the simple model systems presented at the beginning of this overview may be considered for investigation using state of the art energy decomposition procedures, such as Interacting Quantum Atom Model, to quantify the strength of each of the interactions involved. In addition, charge density approaches (viz. ELF, RDG and QTAIM) may also be applicable to most of the crystal systems to explore the nature of various interlayer interactions that occur in the low dimensional crystal systems. These approaches allow one to visualize the nature of electron density localization, delocalization and depletion, among others. Future studies will address whether such computational and theoretical attempts would enable us to better understand the pnictogen-centered non-covalent bonding environments responsible not only for the stability of the crystals, but also their functionality in areas as diverse as optoelectronics.