Gabriela
Castillo-Toraya
a,
Mesías
Orozco-Ic
a,
Eugenia
Dzib
a,
Ximena
Zarate
b,
Filiberto
Ortíz-Chi
c,
Zhong-hua
Cui
d,
Jorge
Barroso
*a and
Gabriel
Merino
*a
aDepartamento de Física Aplicada, Centro de Investigación y de Estudios Avanzados, Unidad Mérida, Mérida, Yucatán, México. E-mail: gmerino@cinvestav.mx; jorge.barroso@cinvestav.mx
bInstituto de Ciencias Químicas Aplicadas, Facultad de Ingeniería, Universidad Autónoma de Chile, Santiago, Chile
cCONACYT-Universidad Juárez Autónoma de Tabasco, Centro de Investigación de Ciencia y Tecnología Aplicada de Tabasco, Cunduacán 86690, Tabasco, México
dInstitute of Atomic and Molecular Physics, Key Laboratory of Physics and Technology for Advanced Batteries (Ministry of Education), Jilin University, Changchun 130012, China
First published on 8th April 2021
Among the list of planar tetracoordinate atoms, fluorine is missing. So far, there are no theoretical or experimental reports suggesting their existence. Herein, we introduce the first six combinations (FIn4+, FTl4+, FGaIn3+, FIn2Tl2+, FIn3Tl+, and FInTl3+) whose global minima contain a planar tetracoordinate fluorine. The bonding analyses indicate that the interactions between the fluorine and the peripheral atoms are significantly electrostatic, which is also reflected in the electronic delocalization. As opposed to other planar tetracoordinate systems with carbon, nitrogen, or oxygen atoms, the fluorine in the ptFs does not act as a σ-acceptor, restraining any back-donation. On the other hand, σ-electrons show a diatropic response, which would characterize these clusters as σ-aromatic.
So, the theory has played a fundamental role in predicting and understanding these peculiar species.20–23 Although the list is long, it does not include a theoretically or experimentally reported planar tetracoordinate fluorine (ptF) atom, which does not make sense since the D4h CAl42−, NAl4−, and OAl4 structures are the global minima on their corresponding potential energy surfaces (PESs). Why is the isoelectronic FAl4+ not planar? Is it possible to introduce fluorine into the flat world? Herein, we systematically explored the PESs of all viable pentaatomic clusters containing a fluorine atom combined with main group elements, totaling 18ve and charges from +2 to −2. This amounts to 440 combinations and around 100000 optimizations. Gratifyingly, we found that the global minima of FIn4+, FTl4+, FGaIn3+, FIn2Tl2+, FIn3Tl+, and FInTl3+ contain a ptF atom.
The Wiberg Bond Index (WBI)32 and Natural Population Analysis (NPA)33 were estimated following the NBO scheme.34 The AdNDP method, developed by Zubarev and Boldyrev,35 was employed for better insight into the bonding of these systems. The AdNDP leads to a partitioning of the charge density into elements with the lowest possible number of atomic centers per electron pair: n-center-two-electron (nc-2e) bonds, including the Lewis bonding elements as core electrons, lone pairs, and 2c–2e bonds.
Also, the analysis of the induced magnetic field36–38 (Bind) and induced current density (Jind)39–41 was analyzed using the programs Aromagnetic42 and GIMIC,39–41 respectively. The molecular magnetic response to a homogeneous external magnetic field is computed on a rectangular grid with the molecule in its center. These calculations were performed at the TPSS/TZP-DKH43 level.
![]() | ||
Fig. 1 TPSS-D3(BJ)/def2-TZVP structures of the FX4+ (X = Al, Ga, In, Tl) series. Bond lengths are in Å. |
With this in mind, we explored the possibility of combining two elements of group 13 (except boron because it creates a small cavity) to produce more clusters with a ptF. Notably, there is no ptF as the global minimum if the structure contains aluminum. The only ptFs, but as local minima, are FIn3Al+ and FTl3Al+. It might be assumed that replacing a Ga atom in FGa4+ with a heavier element should create a cavity large enough to stabilize a ptF, but it is not so obvious (this only occurs for FGaIn3+). Interestingly, while the putative global minimum of FGaIn3+ (rF–Ga = 2.180 Å) is a ptF, the ptF isomer for FGaTl3+ (rF–Ga = 2.110 Å) is 0.6 kcal mol−1 higher in energy than the structure with the tricoordinate fluorine atom. In fact, the F–Ga bond lengths are shortened in clusters containing Tl, reducing the cavity space. So, using larger atoms does not necessarily lead to larger cavities. That is why the three combinations including In and Tl (FIn2Tl2+, FIn3Tl+, and FInTl3+) are the only ones with perfect cavities for stabilizing a true ptF.
Why are ptFs so elusive compared to their isoelectric analogues CAl42−, NAl4−, and OAl4? Table 1 reveals particular trends among planar tetracoordinate atoms surrounded by Al. First, our computations show that going from CAl42− to FAl4+, there is a charge decrease on the central atom (from −2.62 |e| to −0.88 |e|) with a charge increment on the peripheral atoms (from +0.16 |e| to +0.47 |e|). So, the charge transfer from the Al4 skeleton to the central atom is more effective in the fluorine case. In other words, the ionic character between the hypercoordinate center and the metal fragment increases from C to F.
q | Configuration | WBItotal | Δ H–L (eV) | ν min | |
---|---|---|---|---|---|
FAl4+ | −0.88 | 2s1.97 2px1.97 2py1.97 2pz1.97 | 0.24 | 1.15 | −47 |
FGa4+ | −0.86 | 2s1.97 2px1.96 2py1.96 2pz1.97 | 0.28 | 1.13 | −71 |
FIn4+ | −0.88 | 2s1.98 2px1.96 2py1.96 2pz1.97 | 0.24 | 0.95 | 33 |
FTl4+ | −0.87 | 2s1.98 2px1.96 2py1.96 2pz1.97 | 0.25 | 0.86 | 28 |
FGa3In+ | −0.85 | 2s1.97 2px1.96 2py1.95 2pz1.97 | 0.30 | 0.84 | 26 |
FGa3Tl+ | −0.84 | 2s1.97 2px1.95 2py1.95 2pz1.96 | 0.32 | 0.79 | 24 |
FGa2In2+ | −0.84 | 2s1.97 2px1.96 2py1.95 2pz1.96 | 0.31 | 0.80 | 30 |
FGa2Tl2+ | −0.83 | 2s1.97 2px1.96 2py1.95 2pz1.96 | 0.33 | 0.79 | 25 |
FGaIn3+ | −0.87 | 2s1.97 2px1.96 2py1.97 2pz1.97 | 0.26 | 0.97 | 42 |
FGaTl3+ | −0.86 | 2s1.97 2px1.96 2py1.96 2pz1.97 | 0.27 | 0.89 | 33 |
FIn2Tl2+ | −0.88 | 2s1.98 2px1.96 2py1.96 2pz1.97 | 0.25 | 0.89 | 10 |
FIn3Tl+ | −0.88 | 2s1.98 2px1.96 2py1.96 2pz1.97 | 0.25 | 0.92 | 25 |
FInTl3+ | −0.87 | 2s1.98 2px1.96 2py1.96 2pz1.97 | 0.25 | 0.87 | 31 |
CAl42− | −2.62 | 2s1.60 2px1.76 2py1.76 2pz1.46 | 2.29 | 1.55 | 86 |
NAl4− | −2.27 | 2s1.79 2px1.88 2py1.88 2pz1.71 | 1.33 | 1.79 | 95 |
OAl4 | −1.68 | 2s1.91 2px1.94 2py1.94 2pz1.88 | 0.63 | 1.64 | 94 |
This is also supported by the significant reduction in the total values of the WBI32 of the central atom from 2.29 in CAl42− to only 0.24 in FAl4+, which, of course, is a consequence of the electronegativity difference between the peripheral atoms and the corresponding planar tetracoordinate atom. The electron transfer mechanism is less obvious, while carbon acts as σ-acceptor, and this is compensated by the back-donation from its 2pz electrons to the π-bonding (2s1.60 2px1.76 2py1.76 2pz1.46), the back-donation is null in the fluorine compounds. This is the main reason why these species are so rare. While several ptCs have been predicted and even detected experimentally, only six from 440 compositions contain a planar fluorine atom within the established parameters. For all ptFs, global minimum or not, there is a negligible variation of the fluorine charges and orbital occupancies. The only distinction is a subtle increase in the total WBI values of the gallium clusters due to their relatively short Ga–F bonds. However, the charge transfer is more effective for the global minima, i.e., higher than 0.86 |e| on F. The exception is FAl4+ (−0.88 |e|), but the cavity formed by the aluminum atoms is relatively small. In other words, both cavity size and charge transfer are essential to stabilize the ptF.
AdNDP is a powerful tool for understanding the distribution of electron pairs in the system.35 However, in some cases, several solutions are possible. Our strategy for identifying the most viable ones was to determine those schemes with the highest occupation number (ON). Let us compare two bonding schemes, or solutions, for CAl42− and FAl4+. In the first solution (Fig. 2a and b), there are four 2c–2e E-Al σ-bonds (E = C or F), four 3c–2e Al-E-Al σ-bonds, and one π-bond. In both clusters, the σ-bonds remain almost identical, but in the carbon molecule, the π-bond involves the four ligands, whereas, in FAl4+, it is located at the central atom, i.e., it can be reclassified as a lone-pair. The occupation numbers in all cases are higher than 1.95.
![]() | ||
Fig. 2 (a) First AdNDP solution for CAl42− and (b) FAl4+. (c) Second AdNDP solution for CAl42− and (d) FAl4+. ON stands for occupation number. |
However, a second solution suggests the absence of back-donation from the central to the peripheral atoms (Fig. 2c and d). This bonding layout can be described by a delocalized 4c–2e σ-bond involving the external ring, one lone-pair on each peripheral atom, and four lone-pairs on the central fluorine, i.e., no overlap between fluoride and the X42+ skeleton. This is entirely consistent with the charge analysis and WBI values. The residual electrons for the carbon system are 1.47, and in the case of fluorine, only 0.65. These values, along with the corresponding analysis for N and O (Tables S1–S2†), reveal that, along the row, the “ionic” scheme becomes quite competitive. Besides, it should be emphasized that both systems identify the lone pair on the fluorine atom, and most importantly, the number of residual electrons in the ionic scheme varies from 0.65 to 0.35 as it goes from Al to Tl, indicating a proclivity for the scheme bonding with no orbital overlap, i.e., an electrostatic nature between the fluorine and the peripheral atoms.
The back-donation absence also has consequences on the magnetic response. Let us compare the induced magnetic field36–38 of FIn4+ and CAl42−. In both cases, the magnetic response of the σ- and π-electrons is diatropic, but the first one much stronger compared to that caused by the π-cloud (almost negligible), as Fig. 3 shows. Note that the σ-electrons response of FIn4+ is stronger than in the carbon structure. So, both clusters could be classified as aromatic, predominantly σ-aromatic.44 More details about the magnetic response are provided by the induced current density (Jind).39–41 For CAl42−, there is a strong diatropic ring current around the Al4 skeleton (in deep blue) and a paratropic (counterclockwise) current around the carbon vicinity (in red, see Fig. 3). In contrast, in FIn4+, a local diatropic current is entirely localized at the fluorine atom (light blue). In order to quantify the ring-current strength, the Jind is integrated into a specific area. For the typical reference system, benzene, the integration plane starts in the center of the ring, intersects the C–C bond ending about 4 Å away. This results in a ring-current strength of 12 nA/T. A ring-current strength of 20.5 (FIn4+) and 13.7 (CAl42−) nA/T are obtained using the same integration domain. However, these values include the currents located at the central atoms. So, it is mandatory to avoid these local ring-current effects to quantify only the peripheral ring current. The selected integration planes start at 1.2 (FIn4+) and 0.7 Å (CAl42−) from the central atom, which are inflection Jind points from the local pathway to a ring current pathway. This provides a ring-current strength of 17.1 nA/T for FIn4+, which is smaller than the corresponding one for CAl42− (22.6 nA/T).
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1sc01325d |
This journal is © The Royal Society of Chemistry 2021 |