IUPAC definitions, features and characteristic properties that can be used to identify non-covalent interactions, such as hydrogen bonding, chalcogen bonding and halogen bonding in chemical systems have been promulgated in 2011 
, 2013 
and 2019 
, respectively. A task group has been charged with categorizing tetrel bonding, pnictogen bonding, and other non-covalent interactions, involving the elements of Groups 14–16 of the periodic table 
. The proposed characteristics of these bonding interactions have emerged from observations of a variety of engineered chemical systems in the crystalline, liquid and gas phases, and examining signatures arising from IR [5,6]
, Raman [7,8]
, UV/vis [9,10]
, and NMR [10,11]
methods, and from ab initio and density functional theory calculations [12,13,14]
. The characteristic features vary from system to system because they are dependent on the nature of the electron density distribution associated with the electron donor and electron acceptor fragments driving the non-covalent interactions.
In this featured resviearch, researchersw, we show that nitrogen in molecular entities can act as an electrophile when covalently bonded with appropriate electron-withdrawing atomic domain(s). This electrophilic feature of N can form nitrogen-centered pnictogen bonding interactions (or simply nitrogen bonding interactions) when in close proximity with a negative site on the same or a neighboring molecule, and hence is responsible in part for the stability of the resulting supermolecular or intramolecular entity.
The nitrogen 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.
A nitrogen bond in a molecular entity may be regarded as a σ-hole centered pnictogen bond, especially when the covalently bound nitrogen features a positive σ-hole along the R–N bond extension, where R is the remainder of the molecular entity and is in attractive engagement with a negative site in the same or a neighboring molecule. A σ-hole is defined as a charge density deficient region on an atom A that appears along the outer extension of the R–A covalent bond [15,16]
. By contrast, a nitrogen bond may be regarded as a π-hole centered pnictogen bond when the nitrogen in molecular entities features a π-hole 
on its electrostatic surface and has the ability to engage attractively with a negative site in a neighboring molecule, or a site that has an electron density different to that of the π-hole, thus providing stability to the geometry of the resulting structure. A σ-hole interaction in a chemical system is generally observed to be directional, whereas a π-hole is non-directional. There are many studies on σ-hole and π-hole 
interactions in chemical systems, including similarities, differences, controversies and misconceptions [18,19,20,21,22,23,24,25]
Nevertheless, nitrogen is a crucial constituent in the development of high-density materials, and, of course, is crucial for ammonia synthesis [26,27,28,29]
. Apart from its well-known use as a feedstock for the production of fertilizers 
may be on the brink of a viable ammonia economy 
. “Green ammonia” is a carbon-free hydrogen-containing compound with immense interest due to its high density and high hydrogen storage capability [29,32,33,34]
considered nitrogen because a fundamental understanding of its modes of interaction in forming complex chemical systems would surely elevate theour
current knowledge of nitrogen-centered pnictogen bonding (i.e., nitrogen as a pnictogen bond donor). The role of pnictogen bonding containing heavier members of the pnictogen family (viz. P, As, Sb and Bi) in catalysis has also been demonstrated on several occasions [35,36,37,38]
, and in other areas, such as anion transport and recognition chemistry 
, solution and gas-phase chemistry [40,41,42]
, computational chemistry 
, supramolecular chemistry 
, coordination chemistry [44,45]
, medicinal chemistry 
, crystallography [47,48]
and crystal engineering [49,50]
, among other research fields.
Nitrogen is the lightest member of the pnictogen family (Group 15), and the third most electronegative element after fluorine and oxygen. When it is present in a molecular entity, it is well known to be a site for interaction with electrophiles 
. HCN, N2
, NO, HCCCN, CH3
CN, and NO−3
are a few examples where nitrogen is entirely negative, and hence, serves as a site to interact with electrophiles to form molecular complexes in the gas phase or adducts in the solid state. Nitrogen in molecules, such as NH3
and in amines, is often found, both experimentally and theoretically [52,53]
, to be negative. This is a consequence of its strong electronegative character. It is also less polarizable than the other members of the pnictogen family. Its lone-pair electrons, if visualized from a Lewis structure viewpoint, dominate along the extension of the R–N covalent bonds in compounds in which it is covalently bonded and is an electron density donor towards electrophiles, forming hydrogen bonds, halogen bonds or chalcogen bonds, and acts as a Lewis base to metal ions 
. This is the probable reason why fluorine, oxygen and nitrogen—being highly electronegative but with low polarizability—often have a negative σ-hole 
. However, under certain circumstances, nitrogen can act as an acceptor of electron density, for example, in hypervalent non-covalent interactions, as demonstrated by Chandra and coworkers 
. They observed an intermolecular interaction between nitromethane, CH3
at low temperature within an inert gas matrix, which was characterized by IR spectroscopy and supported by first-principles calculations. The nitrogen in CH3
acts as an electron density acceptor from the nitrogen in NH3
, forming a CH3
dimer and is stabilized by an N···N pnictogen bond. The directional prevalence of this interaction over the C–H···N and N–H···O hydrogen bonding interactions was shown by ab initio calculations, demonstrating that σ-hole/π-hole driven interactions lead to the formation of the dimer, despite nitrogen’s low polarizability.
2. The Solid-State Structure of Dinitrogen, N2
begin by looking at dinitrogen itself. Solid N2
displays exceptional polymorphism under extreme conditions. Some of its known phases include α, β, γ, δ, ε, ι and θ [115,116,117]
. Because the phases of the crystal depend on an external agency (temperature and pressure), the packing is different in each phase. Consequently, the intermolecular interactions between the N2
molecules in these crystals are different.
The α-phase (third allotrope) of solid N2
, with Z
= 4, space group Pa-3
, Figure 1
a, determined using electron diffraction measurements 
, is mainly stabilized by π···N (lone-pair) interactions between the N2
molecules. The centers of the N2
molecules are on an f.c.c. lattice with each molecule pointing in a different direction. The intermolecular distance between N in one molecule and the midpoint of the π bond in a neighboring molecule is 3.565 Å. RW
attribute this interaction to π-centered pnictogen/nitrogen bonding.
Figure 1. The nature of bonding interactions between the N2 molecules in the crystals of (a) α-N2, (b) β-N2, (c) γ-N2, and (d) ε-N2. (e) The 0.001 a.u. isodensity mapped MESP plot of an isolated N2 molecule, obtained at the MP2(full)/aug-cc-pVTZ level of theory. The local most minima and maxima of potential are marked by tiny circles in blue and red, respectively. Selected bond lengths and bond angles are in Å and degree, respectively. Crystallographic axes are not shown in (a,c) for clarity. The crystal symmetry is shown for each case, together with ICSD references.
The β-phase of hexagonal N2
= 2), Figure 1
, comprises two types of intermolecular bonding modes, including N···N and π···π. The first of these contacts occur between the N2
molecules along the crystallographic c
-axis, as shown in Figure 1
b, with r
(N···N) = 3.015 Å and ∠N≡N···N = 132.3°. It is perhaps a Type-Ib topology of pnictogen bonding. The second type of contact is the result of a slip parallel arrangement between the N2
molecules in the crystal along the crystallographic a
axes, with an intermolecular distance of 3.861 Å between the centroids of a pair of triple bonds on two neighboring N2
(π···π) = 3.861 Å).
The γ-phase of tetragonal N2
= 2), Figure 1
, consists of N···N and π···N(lone-pair) interactions. The former ones are longer than the latter. Similar to the β-phase, the N···N contacts follow a Type-Ia topology of bonding, with r
(N···N) = 3.292 Å and ∠N≡N···N = 121.8°, but are weaker. The π···N (lone-pair) pnictogen bonded interactions have r
(N···π) = 3.415 Å and are stronger than the π···π interactions in β-N2
. Regardless of the nature of the intermolecular interactions, they are all marginally shorter or somewhat longer than twice the vdW radius of N, 3.32 Å (rvdW
(N) = 1.66 Å 
The ε-phase of rhombohedral N2
(Z = 24), Figure 1
, has two types of N≡N···N bond distances, r
(N···N) = 2.727 Å and 2.801 Å, corresponding to ∠N≡N···N of 133.1° and 120.9°, respectively, and are less than twice the vdW radius of N, 3.32 Å. Additionally, N(lone-pair)···π interactions may be present between the N2
molecules in the crystal, with r
(N···π) = 2.958 Å.
The crystal structures of the high-pressure δ and δ* phases of nitrogen were also investigated using single-crystal X-ray diffraction (not shown) 
. The structure of the δ phase is isostructurally similar to that of γ-O2
. Thus, it comprises spherically disordered molecules, with a preference for avoiding pointing along the cubic ⟨100⟩ directions, and disk-like molecules with a uniform distribution of orientations. The structure of the δ* phase is tetragonal, and the space group was identified unambiguously as P
at 14.5 GPa.
ResWearchers did not observe any potentially directional interactions in either of the pressure-induced phases of N2. In all cases, the N···N interactions follow a Type-Ia/Type-Ib topology of bonding and occur between interacting sites of dissimilar electron density. The possibility of the N···N and π···N (lone-pair) interaction can be inferred from the MESP model of an isolated molecule, Figure 1e. The delocalized bonding region in N≡N is equipped with a belt of positive potential characterized by near equivalent local most maxima (tiny red spheres). In contrast to what might be expected on the surface of a covalent bound halogen in a molecular entity, such as, for example, in HBr and HCl, rwesearcher observed that each N along the N≡N bond extension is accompanied by a local most minima of potential (tiny blue spheres in Figure 1e) and is very negative. This signifies that there is a buildup on charge density on the surface of N along the outer N≡N and the buildup is significant compared to the lateral sites on the same atoms.
3. The Nitrogen Trihalides, NX3 and Their Crystal Structures
ResWearchers now discuss a set of other examples in which covalently bound nitrogen acts as an electrophile, as observed in crystalline systems reported over many years. In these systems, covalently bound nitrogen, as in NX3 (X = F, Cl, Br, I), has a positive region on the R–N bond extension that is capable of attracting a negative site on an identical, or in a different molecule. This can be appreciated by looking at the MESPs shown in Figure 2, obtained using ωB97XD/Jorge-ATZP. The top view of Figure 2 shows nitrogen having a positive electrostatic potential along the X–N bond extensions, regardless of the nature of the halogen derivative attached to it. It is largest (16.6 kcal·mol−1) along the F–N bond extensions in NF3. As one passes from NF3 through NCl3 to NBr3, there is a decrease in VS,max from 16.6 kcal·mol−1 to 5.4 kcal·mol−1 to 2.7 kcal·mol−1, and the stability of the σ-hole on N in NX3 occurs in the order X = F > Cl > Br. By contrast, the lone-pair dominated region on N is negative, and the negativity increases in the order NF3 (−4.1 kcal·mol−1) > NCl3 (−11.8 kcal·mol−1) > NBr3 (−14.2 kcal·mol−1). These results demonstrate the amphoteric character of N in NX3.
Comparison of the ωB97XD/Jorge-ATZP calculated 0.001 a.u. isodensity envelope mapped potential on the electrostatic surfaces of NX3
(X = F, Cl, Br, I) molecules. Selected VS,max
values in kcal·mol−1
are shown, which are the local most maximum and minimum of potential (red and blue circles), respectively. Two views of each MESP graph are shown. In the top view, bonded N faces the reader, whereas, in the bottom view, it is the three halogens that face the reader. The Quantum Theory of Atoms in Molecules (QTAIM 
)-based molecular graph is shown for each entity, with the circles in green representing bond critical points and bonds in atom color.
The MESP results for NI3 are also included in Figure 2, but clearly do not follow the same trend noted above. This may be an artifact of the basis set rwesearcher used for these calculations, and we clarify this further below. Whilst the numerical values of potential may be out of line, the reactive nature of the N site along and around the I–N bond extensions in NI3 behaves qualitatively similar to that found for the other three NX3 entities.
As shown in the bottom view of Figure 2, VS,max in NX3 increases monotonically from X = F through I along the N–X bond extension, indicating that the magnitude of the σ-holes associated with these potentials increases in the same order and the σ-holes are surrounded by belt-like negative potentials. In other words, the least polarizable F has the smallest σ-hole on N–F bond extensions and the most polarizable I has the largest σ-hole on the N–I bond extensions. The difference in VS,max between NBr3 and NI3 is not large and, as we commented above, the MESP values for NI3 with the chosen basis set used may not be very reliable.
In all cases, the central region of the triangular face formed by the three X atoms in each NX3 is found to be positive. The positive and negative nature of VS,max and VS,min on the surface of each constituent atomic domain in NX3 suggest that each of them not only has the ability to form complexes with another identical (or different) molecule but also has the ability to act both as an acid and a base.
In order to determine whether or not the inconsistency in the nature of the negative and positive regions on the surface of specific atoms in the series of four molecules NX3
is a basis set artifact, resw
examined the same properties using another basis set, def2-TZVPD, available in the basis set exchange library 
have also adopted a higher level of theory, MP2, for the same calculation as it is one of the simplest and most useful levels of theory beyond Hartree–Fock that accounts for effects arising from electron–electron correlation. The results are summarized in Table 1
Table 1. ωB97XD/def2-TZVPD and MP2(full)/def2-TZVPD computed 0.001 a.u. isodensity envelope mapped the local most minima and maxima of potential on the electrostatic surfaces of NX3 (X = F, Cl, Br, I) molecules. VS,max and VS,min values in kcal·mol−1.
|Extrema of potential
|N–X (one on each X)
|X–N (one on each extension)
|Center of the triangular X3 face
|On lone-pair region on N
|On F (within the X3 face)
|On X (opposite of the X3 face)
As can be seen from the data in Table 1, and regardless of the correlated method used in conjunction with the def2-TZVPD basis set, the N has a positive σ-hole along the X–N bond extensions. These σ-holes become less positive as the size of the halogen in NX3 increases from F down to I, and there are three equivalent holes on N in each NX3. By contrast, the halogen derivative has a single σ-hole on its surface along the N–X bond extension, so there are also three positive σ-holes on the surfaces of the three X atoms that are also equivalent (only one is listed in Table 1). They systematically increase with the increasing size of the halogen in NX3, in agreement with what might be expected from their polarizabilities.
3.1. Nitrogen Trifluoride, NF3
The crystal structure of NF3
was reported only recently, and corresponds to the low-temperature α-phase 
. Since powder neutron diffraction measurements were performed, it is expected that the intermolecular geometry is more accurate than might have been obtained from X-ray diffraction measurements.
The α-phase of NF3 crystallizes in the orthorhombic space group Pnma, with lattice parameters a = 6.71457(13) Å, b = 7.30913(14) Å, c = 4.55189(8) Å, cell volume (V) = 223.396(7) Å3, and Z = 4 at T = 6 K. The β-phase of NF3 corresponds to a high-temperature phase and was observed to be a plastic crystal (space group P4/mnm) with lattice parameters a = 15.334(6) Å, c = 7.820(3) Å, V = 1838.6(12) Å3, and Z = 30 at T = 60 K. It was suggested that the crystal structure of this latter phase is closely related to that of the Frank–Kasper sigma phase, but the one deposited in the ICSD (ref. code 1891641) does not contain the geometry of β-NF3 (the fluorine atoms are missing).
The unit-cell of the α-phase comprises four units of NF3, as shown in Figure 3a. It shows the intermolecular bonding modes between the building blocks responsible for the unit-cell of the NF3 crystal; these are consistent with the attraction between the positive and negative sites of electrostatic potential localized on different atomic domains (Figure 3b). In particular, Type-IIa N–F···F and F–N···F non-covalent interactions are observed and are inferred from the intermolecular distances and angles of interaction. The former is less directional but shorter than the latter (Figure 3a).
Figure 3. (a) The unit-cell of α-NF3 (CSD ref. code 1891640), showing selected N···F and F···F σ-hole-centered intermolecular interactions and the intermolecular angles for the approach of the electrophile. (b) Illustration of selected local maximum and minimum of electrostatic potential (red and blue spheres, respectively) mapped on the 0.001 a.u. isoelectron density surface of NF3 obtained at the ωB97XD/aug-cc-pVTZ level of theory. (c) The network of N···F and F···F σ-hole centered pnictogen and halogen bonding interactions in 2 × 2 supercell geometry of α-NF3. Bond lengths and bond angles in Å and degrees, respectively. Dotted lines represent intermolecular interactions.
The real nature of the intermolecular interactions between the molecular units in the crystal may not be apparent in inspecting the unit-cell alone; a periodic extension of the unit-cell is necessary. The 2 × 2 supercell structure of α-NF3 is shown in Figure 3c. Although there are many more interactions between the molecules of NF3 in the crystal, rwesearcher have highlighted in Figure 3c only the prominent bonding modes formed by the N site. As revealed by the MESP model (Figure 2a), the three positive sites along the F–N bond extensions do indeed donate pnictogen-centered σ-hole bonds to the lateral portions of the F sites in the nearest NF3 units (Figure 3a). Of the three, two are equivalent and the other is longer (3.169 Å vs. 3.368 Å). They are all directional since ∠F–N···F for each of the two equivalent interactions is 172.3° and that for the longer bond is 174.5° (Figure 3a).
In crystalline α-NF3, the nitrogen in the NF3 molecule acts as a hexa-furcated center in donating σ-hole bonds. It donates three highly directional σ-hole bonds (vide supra) and two relatively less directional equivalent σ-hole bonds (r(N···F) = 3.301 Å and ∠F–N···F = 143.4°) that are caused by the attraction between the positive site on the central region of the surface formed by the three fluorine atoms in one molecule and the negative region dominated by the lone-pair of a neighboring molecule, thereby forming five pnictogen bonds (Figure 4a,b). Because nitrogen has a negative site (readily appreciated in its Lewis structure) evidenced by the VS,min = −3.6 kcal·mol−1 shown in Figure 2a, it shows a tendency to accept a σ-hole bond from the covalently bound fluorine of another interacting NF3 molecule, responsible for the formation of an F···N halogen bond (Figure 4b). The back-to-back arrangement between the NF3 molecules, Figure 4a, causing the additional intermolecular interactions noted above, is a result of a π(N)···π(N) interaction. The behavior of N in forming so many non-covalent interactions appears to be unique to this system.
Figure 4. (a) Nature of N···F pnictogen bonds in the α-phase of the NF3 crystal. (b–d) Nature of the local topology of intermolecular bonding interactions around N and F in NF3. Bond lengths and bond angles in Å and degrees, respectively. The symbol “σ” on N/F in (b–d) refers to the covalently bonded atom donating the σ-hole.
Our investigation shows that fluorine has the ability to act both as a hexa- and tetra-furcated center, as shown in Figure 4
c,d, respectively. It is evident from Figure 4
c that fluorine is an acceptor of three σ-hole bonds, thus forming one F–N···F and two N–F···F halogen bonds; it is a donor of a single σ-hole bond, forming an N–F···N halogen bond; and it shows capacity to form two or more N–F···F–N halogen-bonded Type-Ia and/or Type-Ib interactions—all within a distance of 2.90–3.20 Å. The Type-Ia N–F···F–N halogen-bonded contacts, such as those marked in Figure 4
d, are commonly observed in fully fluorinated compounds, for example, the crystal of C6
, and in similar compounds [126,127,128]
. Clearly, the complex topology of bonding between the NF3
molecules in the crystal structure requires a variety of theoretical studies to detail the nature and the strength of the interactions involved.
3.2. Nitrogen Trichloride, NCl
The crystal structure of NCl3
was reported in 1975 
. It crystalizes in the orthorhombic space group Pnma
. There are 12 molecules in the unit-cell, as shown in Figure 5
a. The mean N–Cl bond distance is 1.75(1) Å and the mean bond angle is 107(2)°. The insights gained from the values of VS,max
on the surface of an isolated NCl3
molecule above (Figure 2
b) can now be used to understand its chemical reactivity.
Figure 5. (a) The unit-cell of the NCl3 crystal (ICSD ref. code: 4034). (b) Illustration of the nature of Cl···Cl intermolecular interactions between the molecular building blocks. (c) Illustration of the nature of the N···Cl interactions. (d) Nature of the local topology of bonding around N in NCl3, showing N is a penta-furcated center. (e) Nature of the reactive sites on the electrostatic surface of the NCl3 molecule responsible for the N···Cl, Cl···N and Cl···Cl interactions, computed with ωB97XD/aug-cc-pVTZ. (f) The nature of various N···Cl and Cl···N interactions in the crystal. (g) The nature of Cl···Cl and Cl···N halogen-bonded interactions in the crystal. Values of bond lengths and bond angles in Å and degrees, respectively.
The N–Cl···N and N–Cl···Cl intermolecular interactions holding the NCl3 molecules together in the crystal can be rationalized by referring to Figure 5b,c, respectively. The extent to which covalently bound N in the NCl3 molecule is engaged with the negative sites on neighboring NCl3 molecules, and producing the observed solid-state structure, is illustrated in Figure 5d. The bonding features are as expected when positive and negative sites on molecular surfaces are in close proximity, inferred from the MESP of NCl3 (Figure 5e).
The pattern of intermolecular interactions between bonded N and Cl in the NCl3 crystal is illustrated in Figure 5f. It is evident that each N donates three σ-holes. Two of these, the attraction between a covalently bonded N in an NCl3 molecule and chlorines on neighboring molecules, generate a zig-zag chain-like pattern. They are inequivalent, in contrast to the analogous interactions observed in the NF3 crystal described above. The r(N···Cl) (∠Cl–N···Cl) for these two contacts are 3.960 Å (160.8°) and 3.894 Å (171.6°), respectively, suggesting a Type-IIa nitrogen bonding topology. The remaining σ-hole bond formed by N, Figure 5d,f, is longer and quasi-linear (r(N···Cl) = 4.358 Å and ∠Cl–N···Cl = 173.2°), and is a result of the packing of molecules in the crystal. As shown in Figure 5d, each N site serves as a penta-furcated center in accepting and donating σ-hole bonds.
The three pnictogen bonds are augmented by σ-centered N–Cl···N and N–Cl···Cl halogen bonds in stabilizing the crystal. These Type-IIa N–Cl···N bonds are highly directional (r(Cl···N) = 3.190 Å and ∠N–Cl···N = 177.1°), and are formed between the σ-hole on the N–Cl bond extension in one molecule and the lone-pair dominated region on N in another interacting molecule (Figure 5g).
The three N···Cl contacts highlighted in yellow in Figure 5
f are repeated throughout the crystal. Two of them are equivalent (r
(N···Cl) = 3.298 Å) and the other slightly longer (r
(N···Cl) = 3.360 Å). They are significantly bent, with ∠Cl–N···Cl of 140.3° and 140.1°, respectively. We assign them to be Type-IIb based on the classification of bonding topology provided in Scheme 1.
They are enforced by the attraction between the positive site in the central region of the triangular face formed by three Cl atoms in the NCl3
molecule and the lone-pair dominated negative site on the N of a neighboring molecule. One might characterize this as an N(π)···N(π) pnictogen bond, consistent with the surface extrema revealed by the MESP model (Figure 5
Figure 5g shows the pattern of occurrence of Type-IIa N–Cl···Cl halogen bonds. Those having a Cl···Cl bond distance (and ∠N–Cl···Cl) of 3.580 Å (168.4°) are more directional than those with the corresponding values of 3.548 Å (157.6°). Indeed, they are shorter and less directional than Cl···N halogen bonds with bond distances (angles) of 4.305 Å (173.8°), Figure 5g. There are numerous Cl···Cl contacts present in the crystal (not shown).