Abstract

The structural, relative stability, electronic, IR vibrational, and thermodynamic properties of asymmetric clusters (CH3FBN3)n (n = 1–6) are systematically investigated using density functional theory (DFT) method. Results show that clusters (CH3FBN3)n (n = 2–6) form a cyclic structure with a B atom and a Nα atom binding together. Five main characteristic regions are observed and assigned for the calculated IR spectra. The size-dependent second-order energy difference shows that clusters (CH3FBN3)3 and (CH3FBN3)5 have relatively higher stability and enhanced chemical inertness compared with the neighboring clusters. These two clusters may serve as the cluster-assembled materials. The variations of thermodynamic properties with temperature T or cluster size n are analyzed, respectively. Based on enthalpies in the range of 200–800 K, the formations of the most stable clusters (CH3FBN3)n (n = 2–6) from monomer are thermodynamically favorable. These data are helpful to design and synthesize other asymmetric boron azides.

1. Introduction

The extensive studies on boron azides are developed as these molecules have been shown to be good precursors for boron nitride (BN) film deposition [1]. The chemistry of boron azides commenced in 1954 with the synthesis of boron triazide B(N3)3 [2]; since then, it developed slowly. In 1963, Paetzold produced trimeric (Cl2BN3)3 from the reaction of LiN3 and BCl3 in CH2Cl2 solution [3]. The molecular structure of dichloroboron azide containing six-membered boron nitrogen heterocycle with diazo groups bounded to nitrogen atoms was determined by Müller in 1971 [4]. Boron dihalide azides (BX2N3)3 (X = F, Cl, or Br) from the reaction of BCl3 with trimethylsilyl azide in CH2Cl2 solution were also recorded by Wiberg et al. in 1972 [5]. Dehnicke reported the infrared spectra of the monoazide products, I2MN3 (M = B, Al, Ga), and observed the formation of oligomers of these species in 1978 [6]. Several alkyl and arylboron azides have been synthesized using similar methods. Oligomerization for Me2BN3 was established by 11B NMR spectroscopy in 1966 [7]. The chemistry of alkylboron azides was further investigated, and the development of this field has been reviewed by Paetzold, Fraenk, and coworkers [816].

Although the chemistry of boron azides in the condensed phase has been extensively discussed in the literature, research studies of these species in the gas phase are rather limited. With the goal to use boron azides as the single-source precursors (SSP), knowledge of their gas phase stability and thermodynamic properties becomes essential. The structure of monomers Cl2BN3 and (CH3)2BN3 has been explored using ab initio calculation, and the thermodynamic stabilities of X2BN3 with respect to its dimerization and trimerization have been gained and discussed [1720] from the theoretical view. Previous theoretical studies of the boron azides mostly focused on the symmetric boron azides.

The effect for the break of the symmetry of substituted boron azides is hardly considered. Here, in accordance with previous theoretical studies on the asymmetric clusters of inorganic boron azides [21, 22], and experimental studies on the asymmetric clusters of organic gallane azides and aluminum azides (RR’MN3)n (M = Al, Ga; R = CH3; R’ = H, Cl, Br; n = 3−4) [23, 24], the structure, stability, IR spectra, and thermodynamic characteristics of the asymmetric clusters of organic boron azides (CH3FBN3)n (n = 1 – 6) were theoretically discussed in detail. These discussions will provide fundamental data and references for experimentalist to design and synthesize the novel boron azides in the future.

2. Computational Methods

As is well known, it is important to choose an appropriate basis set to give the accurate description of clusters’ structures and energies. Usually, a substantial size of basis set is required. However, the size of clusters studied in this work excluded the use of a very large basis set, and hence, all calculations were performed using the DFT-B3LYP method with the 6-31G basis set [25, 26] via the Gaussian 09 program package with the default convergence thresholds [27]. To ensure the adequacy of this basis set, we also optimized the clusters with the 6-311 + G basis set. As shown later on, results obtained from these two basis sets were similar except for slightly numerical differences. The energies of all clusters were also evaluated using the 6-311 + G basis set. The initial configurations are searched by two ways: (1) by considering the numbers of known clusters (HClBN3)n (n = 1–6) in previous works [21] and (2) by placing -CH3 groups and F atoms at various substitutional sites (H or Cl) on the basis of optimized (HClBN3) geometries. No symmetry constraints were applied. Lots of rationally initial one-, two-, and three-dimensional configurations were built to seek the most stable structures. In this way, a large number of optimized isomers for the asymmetric clusters (CH3FBN3)n (n = 1–6) were obtained. Frequencies as well as their respective IR intensities were then calculated, and the lack of imaginary frequencies confirmed that the true minimum was obtained in each case. According to the previous study, the calculated frequencies were scaled uniformly by 0.96 to approximately correct the systematic overestimation [28].

3. Results and Discussion

3.1. Structures and Charge Distribution

A number of isomers are calculated at each size. Here, structures of clusters (CH3FBN3)n (n = 1 − 6) having the lowest energy were focused because all properties of clusters (CH3FBN3)n (n = 1–6) were calculated based on the lowest energy. Two stable structures of the monomer CH3FBN3 were obtained (connectivity: CH3FB−Nα−Nβ−Nγ) with a slight difference in geometrical parameters. The clusters (CH3FBN3)n (n = 2 – 6) are produced by the head-to-tail oligomerization of the CH3FBN3 monomers, which is a starting point for the oligomerizations. Two dimers, seventeen trimers, sixty-four tetramers, two hundred and fifty-six pentamers, and five hundred and thirty-six hexamers are obtained in this manner. Judged by the total energies, the most stable isomers at each size are labeled as 1, 2, 3, 4, 5, and 6 and shown in Figure 1. Obviously, B and Nα atoms easily bond together, and B–B and Nα–Nα bonds are not formed in the clusters (CH3FBN3)n (n = 2 − 6). In other words, the clusters (CH3FBN3)n (n = 2 − 6) contain cyclic (BNα)2n structures with alternating boron and α-nitrogen atoms.

The corresponding geometrical parameters are collected in Table 1, and the data in parentheses are calculated results from the 6-311 + G basis set. The results obtained from different basis sets are generally consistent. In detail, bond lengths (except for B−F bonds) from the 6-311 + G basis set are slightly shorter than those from the 6-31G basis set. This whole agreement shows that it is appropriate to choose the 6-31G basis set to compute the titled systems. Therefore, in this work, we only report the geometrical parameters obtained with the 6-31G basis set. For n = 1, the azide group in the monomer CH3FBN3 is slightly bent with a Nα–Nβ–Nγ angle of 173.1°. The calculated B–N bond length of 1.447 Å is between B–N single (1.58 Å) and double bonds (1.37 Å). The Nβ−Nγ bond length at 1.135 Å is shorter than the Nα−Nβ bond length at 1.240 Å. The Nβ−Nγ bond length is between N–N double (1.25 Å) and N–N triple bonds (1.10 Å). For n = 2 – 6, the computed B–N bond length of 1.601–1.638 Å possesses typical character of B–N single bond, which is similar to the B–N length of other covalent boron azides previously determined, such as (F2BN3)3 (1.616 Å) [19] and (C6F5B(N3)2)3 (1.60 Å) [12]. The B−C and B−F bond lengths are in the range of 1.585–1.596 Å and 1.363–1.381 Å, respectively. The computed structural parameters of the azide units in the clusters (CH3FBN3)n (n = 2 – 6) are 1.240–1.255 Å for Nα−Nβ bonds, 1.128–1.133 Å for Nβ−Nγ bonds, and 177.1–179.4° for Nα−Nβ−Nγ bond angle. It is obvious that the azide group is nearly linear. These results are perfectly consistent with those found in (Cl2BN3)3 (1.25 Å, 1.09 Å, 178°) [4], (F2BN3)3 (1.249 Å, 1.129 Å, 179.8°) [19], and (C6F5B(N3)2)3 (1.27 Å, 1.11 Å, 179°) [12]. Therefore, the calculation method in this work is reliable, and it gives credence to our computed structural parameters of clusters (CH3FBN3)n (n = 2 – 6).

Through the above discussions, it is obvious that Nβ−Nγ bond is shorter than Nα−Nβ bond in the clusters (CH3FBN3)n (n = 1 – 6). This can be interpreted as a higher bond order for the terminal N−N bond, showing a preformation of the N2 molecule. Moreover, the cluster size n has an important effect on the geometries. An obvious increase is shown for the length of B-Nα bond with cluster size n increasing from 1 to 2∼6 owing to the tension of the ring. However, B−Nα bond length shows little change with cluster size in the range of 2∼6, which fluctuates in the range of 1.601–1.638 Å. The Nα−Nβ bond lengths increase as the cluster size n increases from 1 to 5 and show little change from n = 5 to n = 6. Similarly, when cluster size n increases from 1 to 3, the B−C and B−F bond lengths increase. However, the change of B−C and B−F bond lengths is not obvious from n = 4 to n = 6. The Nβ−Nγ bond length shortens from 1.135 Å in cluster size with 1 to 1.129 Å in cluster size with 6 and tends to be constant around 1.130 Å with n = 3 – 6. Compared with monomer 1, the increasing bond lengths of Nα−Nβ, B−C and B−F bonds that are outside of rings show that N2 (Nβ−Nγ), -CH3, and F- groups can be easily removed to yield BN material.

To give a deeper understanding on these structures, we calculate the charge distribution of the clusters (CH3FBN3)n (n = 1 – 6) as shown in Table 2. It can be seen that the charge distribution of CH3FBN3 molecule exhibits zwitterionic characteristics with the charge centers of N and B atoms, respectively. As for (CH3FBN3)n (n = 2 – 6), charge distributions are similar to those of CH3FBN3. However, because B atoms in these molecules all reach the high coordination, the characteristic of zwitterionic ion for these clusters is a problem needing more studies.

3.2. Relative Stabilities

All energies are displayed in Table 3, and the data in parentheses are from the 6-311 + G basis set. Obviously, the energies obtained at the B3LYP/6-31G levels are similar to those obtained with the 6-311 + G basis set. The use of larger basis sets has no significant influence on the binding energies, which again shows that the 6-31G basis set is suitable for the clusters studied here. Thus, in this work, we only report the relative stability of the clusters (CH3FBN3)n with the 6-31G basis set. The uncorrected binding energies (Eb), corrected binding energies (Eb-c), average corrected binding energy (Eb-c-ave), and the second-order energy difference (Δ2E) are calculated using the following equations:where E (CH3FBN3)n and E (CH3FBN3) are the total energies of the clusters (CH3FBN3)n (n = 2 – 6) and CH3FBN3, respectively; ZPE (CH3FBN3)n and ZPE (CH3FBN3) represent the zero point energies of the most stable clusters (CH3FBN3)n (n = 2 – 6) and CH3FBN3, respectively; and the 0.96 is a scaling factor [28].

From Table 3, the ratios of ZPE corrections to their binding energies Eb are large for the clusters (CH3FBN3)n (n = 2 − 6). Thus, it is necessary to carry out the ZPE corrections for the binding energies. For intuitive presentation, the calculated energy is also plotted as a function of cluster size n as shown in Figure 2. The curves of Eb (Figure 2(a)) and Eb-c (Figure 2(b)) are seen to increase monotonically with the augment of cluster size n, which means that clusters can continue to gain energy during the clusters’ growth process. The Eb-c-ave values increase sharply with the clusters size n from 2 to 3; then, it approaches to be stable around 13–14 kJ·mol−1 when cluster size n ≥ 3. Therefore, the geometrical structures of clusters (CH3FBN3)n (n = 1 − 6) tend to be stable when cluster size n ≥ 3.

The relative stability of clusters can be also estimated through the Δ2E (n). According to the definition, clusters with positive Δ2E are more stable than those with negative Δ2E. A clear odd-even oscillation behavior is presented in Figure 2(d). The clusters (CH3FBN3)n (n = 1 − 6) at n = 3, 5 correspond to local maximal on the curve of Δ2E, indicating that they are relatively more stable than other clusters. This trend is in excellent agreement with that of the asymmetric clusters of the inorganic boron azides (HClBN3)n (n = 1 – 6) [21]. In particular, the trimer (CH3FBN3)3 possesses the highest stability among all clusters in terms of the largest value of Δ2E. The stable structure of six-member ring of trimer plays a vital role.

3.3. IR Spectrum

The IR spectrum is not only the basic property of compounds and effective method to analyze or identify substances, but also has a direct relationship with thermodynamic properties. However, there are no experimental data for the title compounds. Thus, it is necessary to use the theoretical method to calculate and predict IR spectra for both theoretical and practical reasons. Figure 3 provides the calculated IR spectra of the clusters (CH3FBN3)n (n = 1 − 6) at the B3LYP/6-31G level. Due to intrinsic complexity, it is difficult to assign all vibrational modes. Thus, we only analyze and discuss some typical vibrational modes.

There are five main characteristic regions for clusters (CH3FBN3)n (n = 1 – 6) that correspond to the –N3 asymmetric and symmetric stretching, –CH3 stretching and in-plane bending vibrations, and the fingerprint region. As for monomer CH3FBN3, the –N3 asymmetric and symmetric stretching vibrations with intense peak lie at 2215 and 1360 cm−1, respectively. The modes in the range of 2923–3022 cm−1 are classified to the characteristic –CH3 stretching vibrations. The peak around 1440 cm−1 corresponds to the –CH3 in-plane bending vibrations. The peaks at less than 1210 cm−1 belong to the fingerprint region, which are mainly composed of wagging and scissoring of –CH3 and the stretching of B−F and B−C bonds. As can be seen from Figure 3, similar vibrational modes are observed in clusters (CH3FBN3)n (n = 2 – 6). The regions of 2917–3019, 2202–2243, 1331–1458, 1217–1308, and 0–1204 cm−1 are individually assigned to -CH3 stretching vibrations, -N3 asymmetric stretching vibrations, –CH3 in-plane bending vibrations, -N3 symmetric stretching vibrations, and fingerprint region. In previous studies, the N–N asymmetrical stretching vibration for Cl2BN3, Br2BN3, Me2BN3, (C6F5)2BN3 polymer, and NaN3 crystal were individually located at 2210–2219 [15], 2200–2215 [15], 2128 [7, 14], 2135–2209 cm−1 [11, 12], and 2128–2151 cm−1 [29] from experiments. Our calculated results are in accordance with those in crystal, which shows that the structures of N3 groups in clusters and in crystal are identical. Moreover, the obtained result is similar to that found in the clusters (HClBN3)n (n = 1 – 6) (2217–2250 cm−1) of DFT calculations [21].

Furthermore, compared with monomer, the N3 asymmetric stretching vibration presents the blue-shifted phenomenon for clusters with n = 2 − 6. However, the N3 symmetric stretching vibration and the –CH3 stretching vibration show red-shifted trend. Meanwhile, in three characteristic regions of N3 asymmetric, N3 symmetric, and –CH3 stretching vibration, the number of vibration mode equals to that of azido groups and C–H bonds, respectively. For example, the monomer 1 has three bands at 2923, 2972, and 3022 cm−1 with C–H stretching vibration, and the trimer 3 has three bands at 2213, 2225, and 2240 cm−1 with N3 asymmetric stretching vibration, and the hexamer 6 has six bands at 1237, 1237, 1245, 1245, 1247, and 1255 cm−1 with N3 symmetric stretching vibration.

3.4. Thermodynamic Properties

Thermodynamic functions are important parameters for clusters to predict reactive property in a chemical reaction. According to vibrational analysis and statistical thermodynamic method, the standard molar heat capacity (), standard molar thermal entropy (), and standard molar thermal enthalpy () of the asymmetric clusters (CH3FBN3)n (n = 1 – 6) in the range of 200–800 K are evaluated and tabulated in Table 4. It is obvious that the calculated thermodynamic functions increase with the temperature T raising. For more clear intuitive, taking monomer 1 as an example, the temperature-dependent relationships for , , and are expressed in formulas (2)–(4) and all are plotted in Figure 4(a). These correlations approximate to linear equations due to the small coefficients of T2; namely, the , , and increase linearly with the increase of temperature T. This can be understood from the fact that these three thermodynamic functions are dominated by the weak translational and rotational motions of the clusters at lower temperature, whereas the vibrational motion is intensified and makes more contributions to , , and at higher temperature. The same linear relationships are found for clusters (CH3FBN3)n with n = 2 – 6:where R2 represents the square of correlation coefficient and SD represents standard deviation.

Moreover, the cluster size-dependent relationships for , , and are expressed in formulas (5)–(7) and all are shown in Figure 4(b). All thermodynamic functions increase monotonically with the enlarged cluster size n; in other words, when one more CH3FBN3 is bound, the , , and increase with the average of 108 J·mol−1 K−1, 133 J·mol−1 K−1, and 18 kJ·mol−1, separately. This means that the contribution of monomer CH3FBN3 to the thermodynamic properties of cluster matches the cluster additivity:

In addition, the theoretical enthalpies (ΔH) and Gibbs free energies (ΔG) of various oligomerizations from monomer CH3FBN3 in the range of 200–800 K are compiled in Table 5. The values of ΔH in the processes are negative except for the process 1 ⟶ 2 beyond 700 K, indicating that these oligomerization processes are exothermic. The ΔG at a given temperature is evaluated by the equation ΔG = ΔH − TΔS. The values of ΔG are all positive, which reveals the oligomerizations cannot occur spontaneously in the range of 200–800 K. There are no experimentally thermodynamic data available for the asymmetric clusters (CH3FBN3)n (n = 1 − 6), so the calculated thermodynamic functions and the obtained relationships of them with the temperature and cluster size n may help the experiment to further study the physical, chemical, and energetic properties of the asymmetric clusters (CH3FBN3)n (n = 1 − 6) or other asymmetric boron azides.

3.5. Aromatic Properties

Here, the aromaticity of the studied systems will be discussed. Because our systems are nonplanar, it is not suitable to use the criterion of nucleus-independent chemical shift (NICS) [30]. In order to find out whether a large π bond exists in our studied system, we use the newly proposed localized orbital locator (LOL)-π method [31] that can be used for studying nonplanar system to study the π bonds of our systems. The LOL-π is analyzed and visualized with the Multiwfn [32]. From Figure 5, there is π bond of NNNB in 1 and similar π bonds in 2-6. However, the π conjugation of 2–6 is not as good as that of 1. Meanwhile, Figure 5 also shows no large π bonds in all systems.

4. Conclusions

We have systematically studied the structure, relative stability, IR, and thermodynamic properties of the asymmetric clusters (CH3FBN3)n (n = 1 – 6) at the DFT-B3LYP/6-31G. The B and Nα atoms tend to bond together in clusters (CH3FBN3)n (n = 2 − 6). The variation trend of geometrical parameters shows N2 (Nβ−Nγ), -CH3, and F- groups are easily eliminated and BN material is yielded. The calculated Δ2E of the asymmetric clusters (CH3FBN3)n (n = 1 − 6) exhibits a pronounced odd-even alternation phenomenon with cluster size n increasing, indicating that the cluster at n = 3 is more stable than other clusters. From monomer (n = 1) to clusters (n = 2 – 6), the N3 asymmetric stretching vibration presents blue-shifted trend while the N3 symmetric stretching vibration and –CH3 stretching vibration are red-shifted. The thermodynamic functions (, , and ) of clusters (CH3FBN3)n (n = 1 – 6) show linear increase with the augment of both temperature T and cluster size n of clusters. Judged by the enthalpies in the range of 298.2–600 K, the calculated thermodynamic properties demonstrate that the formation of clusters (CH3FBN3)n (n = 2 − 6) from the monomer is thermodynamically favorable; however, the formation cannot occur spontaneously.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

Qi-Ying Xia conceptualized the study and was responsible for project administration; Deng-Xue Ma, Yao-Yao Wei, and Yun-Zhi Li were involved in data curation; Guo-Kui Liu was involved in resource management; Deng-Xue Ma and Yao-Yao Wei wrote the original draft; and Qi-Ying Xia and Guo-Kui Liu wrote and edited the manuscript. Deng-Xue Ma and Yao-Yao Wei contributed equally to the work.

Acknowledgments

This work was supported by the Natural Science Foundation of Shandong Province (grant no. ZR2017LB011) and the Undergraduate Education Reform Project of Shandong Province (grant no. Z2018S006).