Next Article in Journal
Metabolic Potential of Halophilic Filamentous Fungi—Current Perspective
Next Article in Special Issue
DFT Study of Molecular Structure, Electronic and Vibrational Spectra of Tetrapyrazinoporphyrazine, Its Perchlorinated Derivative and Their Al, Ga and In Complexes
Previous Article in Journal
The Tumor Microenvironment Drives Intrahepatic Cholangiocarcinoma Progression
Previous Article in Special Issue
Vibrational Spectra and Molecular Vibrational Behaviors of Dibenzyl Disulfide, Dibenzyl Sulphide and Bibenzyl
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 8
10.3390/ijms23084190
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Protonation of Borylated Carboxonium Derivative [2,6-B10H8O2CCH3]−: Theoretical and Experimental Investigation

by
Ilya N. Klyukin
1,*,
Anastasia V. Kolbunova
1,2,
Alexander S. Novikov
3,
Aleksey V. Nelyubin
1,
Nikita A. Selivanov
1,
Alexander Yu. Bykov
1,
Alexandra A. Klyukina
4,
Andrey P. Zhdanov
1,
Konstantin Yu. Zhizhin
1,5 and
Nikolay T. Kuznetsov
1
1
Kurnakov Institute of General and Inorganic Chemistry, Russian Academy of Sciences, Leninskii pr. 31, 117907 Moscow, Russia
2
Faculty of Chemistry, Higher School of Economics, Myasnitskaya yl., 101000 Moscow, Russia
3
Institute of Chemistry, Saint Petersburg State University, Universitetskaya Nab. 7-9, 199034 Saint Petersburg, Russia
4
Research Centre of Biotechnology, Winogradsky Institute of Microbiology, Russian Academy of Sciences, 60 let Oktjabrja pr-t, 7, bld. 2, 117312 Moscow, Russia
5
Lomonosov Institute of Fine Chemical Technologies, MIREA Russian Technological University, 78 Vernadskiy Avenue, 119571 Moscow, Russia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(8), 4190; https://doi.org/10.3390/ijms23084190
Submission received: 5 March 2022 / Revised: 7 April 2022 / Accepted: 8 April 2022 / Published: 10 April 2022
(This article belongs to the Special Issue Molecular Structure, Electronic and Vibrational Spectra Theoretical Calculations in Materials Sciences)
Browse Figures
Versions Notes

Abstract

:
The process of protonation of [2,6-B10H8O2CCH3]− was investigated both theoretically and experimentally. The most suitable conditions for protonation of the derivative [2,6-B10H8O2CCH3]− were found. The process of protonation was carried out in the presence of an excess of trifluoromethanesulfonic acid CF3SO3H at room temperature in dichloromethane solution. The structure of the resulting complex [2,6-B10H8O2CCH3*Hfac]0 was established using NMR data and the results of DFT calculations. An additional proton atom Hfac was found to be localized on one of the facets that was opposite the boron atom in a substituted position, and which bonded mainly with one apical boron atom. The main descriptors of the B-Hfac bond were established theoretically using QTAIM and NBO approaches. In addition, the mechanism of [2,6-B10H8O2CCH3]− protonation was investigated.

1. Introduction

The investigation of covalent and noncovalent interactions is one of the main tasks of modern inorganic chemistry [1,2]. Such studies provide an opportunity to better understand the structures of chemical substances and their properties. By examining chemical bonds, it is possible to establish factors that influence their breaking and formation [3,4]. This fact can be used to produce new substances with given properties. A combined theoretical and experimental approach is the best way to investigate such phenomena [5,6].
There are several driving forces behind bonding formation: orbital, electrostatic, and the van der Waals interactions [7,8,9]. Information about the energy characteristics of such interactions is one of the most important descriptors. Theoretical methods allow exploration of the nature of chemical interactions and estimation of their energies in a simple and intuitive way. The application of such popular and well-established methods as QTAIM (Quantum Theory of Atoms in Molecules), ELF (Electron Localization Function), and NBO (Natural Bond Orbitals) analysis provides the opportunity to find essential information about the chemical structure and bonding of target compounds [10,11,12,13,14].
Closo-borates and related compounds have always attracted attention due to their unusual structure. Their electron structures cannot be described in terms of the classic Lewis approach [11,15]. For a correct description of the structure of closo-borate compounds, the terms of 3D-aromacity are used. Closo-borates possess various interesting properties, such as thermodynamic and oxidative stability and biological activity [16,17,18]. These properties determine the potential application of closo-borate systems for the preparation of solid-state batteries [19], magnetic and hydrogen storage materials [20,21], and drugs for BNCT (boron neutron capture therapy) [22,23,24]. Such systems can form various covalent and noncovalent interactions. It is impossible not to mention the tremendous contributions of the pioneers of boron cluster chemistry, such as Lipscomb [25,26,27,28], Hoffman [29,30,31], Wade [32,33,34], Mingos [35,36,37], and Muetterties [38,39,40,41]. Their works have significantly expanded ideas about the structure and properties of such systems. Currently, there is a great deal of interest in the chemistry of closo-borate anions [42,43]. The phenomena of the exo-polyhedral B-X (X = H, C, N, O, F) bonds have been extensively studied [44,45,46], for example, research into the main characteristics of these bonds for hydroxy-derivatives [BnHn-1OH]2− (n = 6, 10, 12) was carried out [47]. For the exo-polyhedral B-O bonds, it has been established that, for the closo-borate anion, covalent interactions increase with an increase in the boron cluster size. Noncovalent contacts of closo-borate anions and related compounds have also been investigated [48]. This type of interaction is weaker than covalent bonds but also has a significant effect on the structure and properties of the systems under study. For example, the interactions of nitrile derivatives with various nucleophiles have previously been studied [34,49,50,51,52]. It was shown that the structure of the products obtained was the result of noncovalent interactions between protons of nucleophile and hydrogen atoms of the cluster anion.
Closo-borate anions possess an excess of electron density and can act as an electron donor for an electrophilic system. They can form a complex with different Lewis acids: protons, AlCl3, MLn, and others [42,53]. Previously, the structure of the [B10H11]− anion was studied both theoretically and experimentally [54,55]. It was shown that Hfacmoves freely in the equatorial belt, with a very low energy barrier (approximately 1 kcal/mol) [20].
The present study focused on [2,6-B10H8O2CCH3]−. This anion is a member of the class of borylated heterocycles. Analogous derivatives based on borylated boroxazole were also investigated for closo-decaborate and dodecaborate anions [56,57,58,59]. Several approaches for the preparation of this class of closo-borate derivatives have been explored previously [60,61]. In addition, the disclosure reaction of borylated heterocycle has been established [62]. [2,6-B10H8O2CCH3]− is a convenient model system for the investigation of the protonation process. This system is more symmetrical than mono-substituted derivatives and the number of possible isomers is reduced, compared to them. The protonation process of a given anion can be readily investigated. It is important to assess the differences and similarities between the protonation of the original anion [B10H10]2− and the derivatives based on it. In addition, the protonation of a given derivative enables the preparation of a derivative with a degree of substitution equal to three. Thus, in the present research, a combined theoretical and experimental study of the protonation of [2,6-B10H8O2CCH3]− was carried out.

2. Results and Discussion

2.1. Experimental Protonation

First, the protonation process of [2,6-B10H8O2CCH3]− was conducted experimentally. As described previously [63], systems of the general form [B10H9L]− have a lesser possibility of forming a protonated complex than the [B10H10]2− anion. A possible reason for this is that the introduction of a positively charged group reduces the ability of the cluster cage to donate its electron density. Anion [B10H10]2− can be easily transferred into [B10H11]− in the presence of trifluoroacetic acid CF3COOH. For the preparation of a protonated complex of the general form [B10H9OR1R2*Hfac]0, however, trifluoromethanesulfonic acid CF3SO3H was used. In the present work, trifluoroacetic acid CF3COOH was also used as a proton donor, but the protonation process did not take place. As in the case of oxonium derivatives [B10H9OR1R2]−, the protonation of the derivative [2,6-B10H8O2CCH3]− was carried out in the presence of CF3SO3H at room temperature (Scheme 1). Dichloromethane CH2Cl2 was used as a solvent. The application of another solvent was impossible due to the interaction of [2,6-B10H8O2CCH3*Hfac]0 with the molecules of the solvents. [2,6-B10H8O2CCH3*Hfac]0 was formed immediately after the addition of the acid. The distinctive feature of this process is that the protonation of the cluster anion requires an excess of CF3SO3H acid. For an investigation of these phenomena, the experiment with sequential addition of CF3SO3H acid to a dichloromethane solution of the initial [2,6-B10H8O2CCH3]− derivative was conducted. The process was monitored by 11B-NMR spectroscopy (Figure S1). Detailed information about the NMR spectra of the initial and target substances and the correlation of all signals are given below in the NMR spectra analysis section. When one equivalent of acid was added to the system, the overall structure of the spectrum characteristic of the original derivative was preserved, but one of the signals was broadened. This signal corresponds to boron atoms in apical positions. The given pattern indicates the possibility of the formation of intermolecular contacts between the carboxonium derivative and the trifluoromethanesulfonic acid molecule. When one more acid equivalent is added to the system, a more significant change in the spectrum was observed, which may indicate the formation of molecular complexes between the [2,6-B10H8O2CCH3]− derivative and the acid molecule. The formation of the target substance was also observed, as indicated by the presence of a broadened signal in the region of 20 ppm. However, the yield of the target product did not exceed approximately 10 percent. A more accurate estimation of the target product yield could not be made due to the overlapping signals of the target product and the cluster-acid complex. At ratios [2,6-B10H8O2CCH3]−:CF3SO3H equal to 1:3 and 1:4, an increase in the yield of the target product was observed. In the case of the 1:4 ratio, the yield of the target product reached approximately 40 percent. Finally, in the presence of 5 equivalents of CF3SO3H, complete protonation of [2,6-B10H8O2CCH3]− with formation of [2,6-B10H8O2CCH3*Hfac]0 occurred. Further addition of the acid excess did not change the form of the spectrum, indicating the completeness of the protonation process. It is noteworthy that [2,6-B10H8O2CCH3*Hfac]0 is stable only as a solvent at room temperature and without air.

2.2. Fukui Function

In the next step, the potential position in the boron cluster for the coordination of an additional proton atom was investigated. The protonation process of [2,6-B10H8O2CCH3]− can be considered to be an electrophilic attack on the boron cluster. The Fukui functions approach is the best tool for investigation of the most likely position in the molecule for an electrophilic or nucleophilic attack [64,65,66]. Previously, it has been shown, with the help of Fukui functions, that electrophilic attack on closo-borate anions is performed predominantly in the apical position [67,68]. In the present investigation, Fukui functions were used with different approaches employed for their calculation (Table S2). As in the case of [BnHn]2− systems [68], the Hirshfeld approach is the best for calculating the electrophilic attack positions. In the present case, electrophilic attack on closo-borate anions was performed predominantly in the apical position (the optimized structures of [2,6-B10H8O2CCH3]− and [2,6-B10H8O2CCH3*Hfac]0 are shown in Figure 1). Thus, the data from the Fukui functions analysis indicated that the Hfac was localized near apical positions (Table S3).

2.3. Protonation Mechanism

The mechanism of [2,6-B10H8O2CCH3]− protonation was investigated. Initially, complex [2,6-B10H8O2CCH3*CF3SO3H]− (Comp) formed endoergonically (by 23 kJ·mol−1, in terms of Gibbs free energies). The main driving force behind the formation of this complex is the dihydrogen bond between the proton atom from CF3SO3H and the equatorial hydrogen atom from the boron polyhedron. The contact length of Heq–H was equal to 1.60 Å. The presence of the dihydrogen bond was proved with the help of QTAIM analysis of the Comp structure (Figure 2). In addition, the C–F⋯H–B contact was detected. According to the main topological descriptors of electron density C–F⋯H–B, contact can be characterized as being very weak.
In the next step, the proton migrated to a boron polyhedral through the formation of the transition state (TS). The preferred location for the proton attack was the face of the polyhedron opposite the substituted position. The proton from CF3SO3H acid in TS connected mainly with the boron atom in the apical position, which was demonstrated with the help of QTAIM analysis (Figure 3). In the case of TS, the contact length between the apical boron atom and the proton atom from CF3SO3H was equal to 1.48 Å. The overall energy barrier of proton migration was 69 kJ·mol−1.
Finally, [B10H8O2CCH3*Hfac]0 formed. The overall process of boron cluster protonation is endoergonical (by 18.5 kJ·mol−1, in terms of Gibbs free energies) (Figure 4). Thus, due to ΔG > 0 overall, a reaction did not occur in the presence of 1 equivalent of CF3SO3H. This finding correlates well with experimental data and the main reason it occurred is the weak stabilization of CF3SO3− in the dichloromethane solution. Additional proton donors are required to stabilize this anion. The use of an excess of trifluoromethanesulfonic acid can provide this. The CF3SO3− can be stabilized by hydrogen bonds between oxygen atoms of this anion and hydrogen atoms of CF3SO3H. For example, complex between CF3SO3− anion and CF3SO3H was observed (Figure S8). This complex possesses strong hydrogen bonds and the formation of a given system is an exergonic process. It is worth noting that the addition of the acid excess also significantly changed the properties of the medium, which in turn can change the solvation energy of the model system under consideration. This aspect may also contribute to the fact that the protonation process becomes exergonic when an excess of the acid is supplied.

2.4. NMR Spectra Analysis

11B NMR spectra were used for the confirmation suggestion based on the Fukui index results and to establish the structure of the resulting complex. It is useful to consider the spectra of the [2,6-B10H8O2CCH3]− anion. Due to the high symmetry of the given system, the two boron atoms in the substituted positions and the two apical atoms are pairwise equivalent. Atoms in equatorial positions are split into two groups of signals (Figure 5). Thus, as described previously, the 11B NMR spectra of the initial anion were characterized by four signals, as follows: a signal at 0.0 ppm with integral intensities I = 2, corresponding to the boron atoms from equivalent substituted positions. This signal did not split in the absence of broadband decoupling; a signal at −7.1 with integral intensities I = 2, corresponding to boron atoms from equivalent apical positions; and signals from equatorial positions appeared at −17.6 (I = 2) and −30.0 (I = 4). This interpretation of the spectrum is based on previously obtained spectral data for several closo-decaborate derivatives with B-O exo-polyhedral bonds. This approach does not, however, enable a complete correlation of signals and leaves some room for speculation.
The theoretical calculation of chemical shielding enabled the full correlation of the signals in the 11B NMR spectra with the boron positions in a cluster cage. Several computational levels for calculating NMR shifts were tested: wB97X-D3, B3LYP, and B97 DFT-functionals with 6–31++G(d,p), IGLO-III, and EPR-III basis sets (Table S5). To evaluate the difference between the theoretical and experimental data, the values of the root-mean-square deviation (RMSD) were used. The wB97X-D3 level of theory provided the worst results, with all basis sets. The values of RMSD were quite large (the range of RMSD values lay in the range 46–47) and applying this method did not enable the deciphering of the 11B spectra of given compounds. Improved results were obtained with the application of the hybrid functionals B3LYP and B97. In the case of the B97 functional, the best values of RMSD were obtained (Table S5). The RMSD value for the B97/IGLO-III level of theory was equal to 3.7. Thus, applying the results of the theoretical calculation, the assumption about the chemical shifts of atoms in the substituted and apical positions proved to be correct. In addition, the signal at −17.6 ppm corresponded to boron atoms from the equatorial belt in the B3 and B9 positions. The signal from −30.0 ppm corresponded to atoms in the B4, B5, B7, and B8 positions.
The 1H NMR spectra of [2,6-B10H8O2CCH3]− were considered. For a simple view of the spectrum and the correct correlation of signals, it is best to use 1H-11B NMR spectra (Figure S2). Some of the signals appeared, however, as broadened lines, making it difficult to identify them. Signals from the methyl substituent appeared at 2.18 ppm, and signals from apical positions at 3.16 ppm were well represented. Signals from equatorial positions appeared at 1.79, 0.39 and 0.12 ppm.
The addition of a proton reduced the symmetry of the resulting system and the spectra of [2,6-B10H8O2CCH3*Hfac]0 were more complicated than those of [2,6-B10H8O2CCH3]−. Boron atoms were split into separate signals. Thus, it was impossible to elucidate the structure of [2,6-B10H8O2CCH3*Hfac]0 by only applying the 1D 11B spectra. To decipher these spectra, the 11B-11B COSY spectra and theoretical calculations were used (Figure S5). This type of spectrum indicates that, contrary to the case of the [B10H11]− anion, the proton in this anion does not migrate along one of the equatorial belts but is localized on one of the edges. Based on the calculated data, the most stable isomer Iso1 was found. In this isomer, the proton was localized on the edge opposite the substituted position. The structure of this isomer is described in more detail below. Similar to the initial anion [2,6-B10H8O2CCH3]−, the B97 functional worked well for the prediction of NMR spectra. The RMSD value for the B97/IGLO-III level of theory was equal to 3.5 (Table S5).
As in the case of [B10H11]−, the anion [2,6-B10H8O2CCH3*Hfac]0 has two non-equivalent apical positions. The first of them is bonded with additional proton atoms Hfac and the second one is not. Signal from the boron atom of the apical position B10, which was not bonded with Hfac, shifted into a low field and appeared at 19.4 ppm. Similar observations were made in the case of the apical position B10 of the [B10H11]− anion. In addition, the existence of this shift was proved by theoretical calculation. According to theoretical calculation, the signal of the apical position B10 appeared at 13.4 ppm. The chemical shift from the apical position B1, which bonded with Hfac, shifted in the high field and appeared at −26.9 ppm. The boron atoms in the substituted positions were also not equivalent in the case of a protonated cluster. The signal at 0.2 ppm correlated with the B10 position of the boron cluster in the 11B-11B COSY NMR spectra, thus it corresponded to the B6-substituted position, which is opposite to the apical position with Hfac. Another substituted position was found at −1.6 ppm. At −15.4 ppm, the signal from the B4,5 positions appeared. Signals from the B3 and B7 positions were relatively close, appearing at −19.0 and −20.1 ppm, respectively. Signals from the B9 and B8 positions had the greatest negative chemical shift value at −26.1 and −27.8 ppm, respectively.
In the 1H NMR spectra of [2,6-B10H8O2CCH3*Hfac]0, the signal of the methyl group was observed at 2.45 ppm (Figure S3). Signals from the exo-polyhedral cluster were in the range 5.16–0.33 ppm. It is noteworthy that the signal at 5.16 ppm corresponded to a hydrogen atom connected with a boron atom in the B10 position. This assumption is based on the result of 1H NMR spectra modeling. According to the results of the calculations, the proton in the downfield region belonged to the given atom. The signal from Hfac, according to theoretical calculations, was in the downfield region at 0.00 ppm. In the experimental spectra, the most upfield signal was at 0.33 ppm. In the 13C NMR spectra of [2,6-B10H8O2CCH3*Hfac]0, the signal of the carbonyl group carbon atom was observed at 194.0 ppm and the signal of the methyl group carbon atom was observed at 20.10 ppm (Figure S4). In the IR-spectra, the B-H bands were shifted from 2490 cm−1 in the initial [2,6-B10H8O2CCH3]− to 2540 cm−1 in [2,6-B10H8O2CCH3*Hfac]0 (Figure S6).

2.5. Structure Elucidation

The structure of [2,6-B10H8O2CCH3*Hfac]0 was examined using DFT calculation. Several possible isomers were calculated (Figure 6 and Figure S7, Table S4). Calculations were performed in the gas phase and considering the solvation effect. In both cases, the isomers differed significantly in terms of energy. The Iso1 and Iso2 isomers were the most thermodynamically stable (Figure 6). In these isomers, an extra proton, in addition to the equatorial boron atoms, was localized on the apical boron atom. The Iso1 in which the proton was localized on the edge opposite the substituted position had the highest negative Gibbs energy among all isomers. The energetic barrier for the migration of the proton atom from Iso1 to Iso2 was equal to 30 kJ·mol−1in the gas phase and 25 kJ·mol−1 in dichloromethane solution. Isomers with the localized proton only on equatorial boron atoms were the least thermodynamically stable and had the smallest negative Gibbs energy values. The energy barrier between Iso1 and Iso3 was 50 kJ·mol−1 for both the gas and dichloromethane phase. This energy barrier was less than that of the anion [B10H11]−. Thus, Iso1 was the most stable isomer among the other possible structures of [2,6-B10H8O2CCH3*Hfac]0.
The geometric parameters of this isomer were discussed in detail (Table S1). The bond length between Hfac and the apical boron atom was 1.28 Å, whereas the analogous parameter for the bond between Hfac and the equatorial boron atom was 1.42 Å. The bond length between equatorial boron atoms and the apical boron atom that bonded with the proton increased compared to the same parameter in the original cluster. This parameter increased from 1.72 to 1.84 Å. The phenomena of B-Hfac contacts were investigated using Wiberg bond order indices and QTAIM analysis (Figure 7). Applying the QTAIM method, a molecular graph of the electron density distribution in [2,6-B10H8O2CCH3*Hfac]0 was obtained. The molecular graph indicated that Hfac has interactions with apical and equatorial boron atoms. The Bap–Hfac interaction was characterized by a greater value of the Wiberg bond order index compared to that of Beq–Hfac. The Bap–Hfac interaction was also characterized by greater values of ρ(r), total energy at the bond critical point, and the delocalization index. These data indicate that Hfac is predominantly bonded with the boron atom in the apical position. An analogous investigation was carried out for the case of the [B10H11]− anion. For this anion, Hfac only bonded with the boron atom in the apical position. This contact was characterized by a longer Bap-Hfac bond and smaller values of the main bond descriptors. Thus, in the case of [2,6-B10H8O2CCH3*Hfac]0 Bap–Hfac has greater covalent interaction than the [B10H11]− anion.
The values of the main bond descriptors indicated that the Bap–Hfac interactions were significantly weaker than those of the covalent interactions between boron atoms and exo-polyhedral hydrogen atoms. It is noteworthy that the values of the main bond descriptors of B-H and B-O bonds were greater in the case of [2,6-B10H8O2CCH3*Hfac]0 than in [2,6-B10H8O2CCH3]−. This indicates that the strength of these contacts increases. A possible reason for these phenomena may be the reduction in the total charge in the system compared to that of the initial anion [2,6-B10H8O2CCH3]−. The atomic charges on the boron atoms are reduced and the electron repulsion in the exo-polyhedral bonds is also reduced.

3. Materials and Methods

3.1. Computational Details

The full geometry optimization of all model structures was carried out at the ωB97X-D3/6–31++G(d,p) level of theory [69,70] with the help of the ORCA 4.2.1 program package (Mülheim an der Ruhr, Germany) [71] (the atom-pairwise dispersion correction with the zero-damping scheme was utilized [70]). The convergence tolerances for the geometry optimization procedure were as follows: energy change = 5.0 × 10−6 Eh, maximal gradient = 3.0 × 10−4 Eh/Bohr, RMS gradient = 1.0 × 10−4 Eh/Bohr, maximal displacement = 4.0 × 10−3 Bohr, and RMS displacement = 2.0 × 10−3 Bohr. Spin-restricted approximation for the model structures with closed electron shells was used. Symmetry operations were not applied during the geometry optimization procedure for any of the model structures. The Hessian matrices were calculated numerically for all optimized model structures in order to prove the location of correct minima on the potential energy surfaces (no imaginary frequencies for all reactants, intermediates, and final products; only one imaginary frequency for transition states). The connectivity of each reaction step was also confirmed using the intrinsic reaction coordinate (IRC) calculation from the transition states [72,73,74]. Solvent effects were considered using the Solvation Model based on Density (SMD) [75]. The natural bond orbital (NBO) method was emplyed, using the NBO7 program package (Madison, WI, USA) [76]. Topological analysis of the electron density distribution, using the Quantum Theory of Atoms in Molecules (QTAIM) formalism developed by Bader [12], was employed with the Multiwfn program, version 3.7 (Beijing, China) [77]. The Cartesian atomic coordinates for all optimized equilibrium model structures are presented in the Supplementary Materials. Visualization of the optimized structures was carried out with the help of ChemCraft program version 1.7 (Ivanovo, Russia) [78]. In the case of the molecular graph showing the results of the topological analysis of the electron density distribution visualization, the Multiwfn program (version 3.7) was employed [77].

3.2. IR Spectra

The IR spectra of the prepared compounds were recorded on an Infralyum FT 02 Fourier transform spectrometer (Lumex Instruments Research and Production Company, Fraserview Place, Vancouver, BC, Canada) in the region of 300–4000 cm–1 and with a resolution of 1 cm–1. Samples were prepared as dichloromethane CH2Cl2 solution.

3.3. NMR Spectra

The NMR (1H, 11B, 13C) spectra of the solutions of the studied compounds in CD3CN were recorded on a Bruker Avance II 300 spectrometer (Ettlingen, Germany) operating at 300.3, 96.32, and 75.49 MHz, respectively, using an internal deuterium lock. Tetramethylsilane and boron trifluoride etherate were used as external references.

3.4. Protonation of ((C6H5)4P)[2,6-B10H8O2CCH3]

CF3SO3OH (0.018 mL, 1.5 mmol) was added to a solution of ((C6H5)4P)[2,6-B10H8O2CCH3] (20 mg, 0.3 mmol) in deuterated dichloromethane CD2Cl2 (0.5 mL), in a dry argon atmosphere. The resulting mixture was kept at room temperature for 10 min. IR-spectra (dichloromethane solution): ʋ(BH) 2590 cm−1, 2555 cm−1. 11B-{1H} NMR (CD2Cl2, ppm): 19.2 (s, B10, I = 1), 1.4 (d, B6, I = 1), −0.2 (d, B2, I = 1), −15.2 (d, B4, I = 1), −15.7 (d, B3, I = 1), −19.1 (d, B5, I = 1), −20.1 (d, B9, I = 1), −24.4 (d, B7, I = 1), −26.1 (d, B1, I = 1), −27.8 (d, B8, I = 1). 1H NMR (CD2Cl2, ppm): 11.75 (CF3SO3H), 7.97 (m, Ph4P+, I = 4), 7.80 (m, Ph4P+, I = 8), 7.66 (m, Ph4P+, I = 8), 5.17 (s, B(10)H, I = 1), 2.68 (s, B(3)H, I = 1), 2.50 (s, CH3, I = 3), 2.03 (s, B(9)H, I = 1), 1.60 (s, B(4)H, I = 1), 1.32 (s, B(5)H, I = 1), 0.93 (s, B(7)H, I = 1), 0.62 (s, B(8)H, I = 1), 0.36 (s, Hfac, I = 1). 13C NMR (CD2Cl2, ppm): 194.7 (O2CCH3), 140.6 (Ph4P+), 139.8 (Ph4P+), 135.8 (Ph4P+), 122.7 (Ph4P+), 118.1 (CF3SO3H), 20.1 (O2CCH3).

4. Conclusions

The process of protonation of the carboxonium derivative [2,6-B10H8O2CCH3]− was investigated. In the case of the carboxonium derivatives of closo-decaborate anions, trifluoromethanesulfonic acid CF3SO3H was used as a proton donor. By considering the reaction mechanisms, the reason for the excess of CF3SO3H was established. The excess acid was required to stabilize the anion CF3SO3−, which in turn caused the total protonation process to become exergonic. In contrast to the anion [B10H11]−, in the case of [2,6-B10H8O2CCH3*Hfac]0, the proton did not migrate along the equatorial belt and was localized on a facet opposite the B atoms bonded with the exo-polyhedral substituent. Based on theoretical modeling data, it was shown that the proton was predominantly bound to the apical boron atom. In addition, the 11B-1H NMR spectra of the final compound were quite complicated and difficult to decipher. Using theoretical modeling data, all signals were correlated with the positions of the boron atoms in the cluster. Theoretical calculations at the B97/IGLO-III level of theory corresponded perfectly with the experimental 11B-1H NMR spectra data. [2,6-B10H8O2CCH3*Hfac]0 can be considered to be a promising synthon for the preparation of trisubstituted derivatives of closo-borate anion.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23084190/s1.

Author Contributions

Manuscript conception, I.N.K. and A.S.N.; writing and original draft preparation, I.N.K., A.S.N. and A.P.Z.; experimental protonation of [2,6-B10H8O2CCH3]−, A.V.K. and I.N.K.; NMR analysis, N.A.S. and A.Y.B.; IR analysis, A.V.N.; Fukui function calculation and structure elucidation, I.N.K., A.V.K. and A.S.N.; mechanism investigation, visualization, A.A.K. and I.N.K.; editing, data analysis, and interpretation, I.N.K., A.P.Z., A.S.N. and K.Y.Z.; project administration I.N.K., K.Y.Z. and N.T.K.; supervision N.T.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Russian Science Foundation (project no. 20-73-00326).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

NMR Studies were performed within the framework of the State Assignment of the Kurnakov Institute (IGIC RAS) in the field of fundamental scientific research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pino-Rios, R.; Inostroza, D.; Tiznado, W. Neither too Classic nor too Exotic: One-Electron Na⋅B Bond in NaBH 3− Cluster. Angew. Chem. Int. Ed. 2021, 60, 12747–12753. [Google Scholar] [CrossRef] [PubMed]
  2. Su, B.; Cao, Z.-C.; Shi, Z.-J. Exploration of Earth-Abundant Transition Metals (Fe, Co, and Ni) as Catalysts in Unreactive Chemical Bond Activations. Acc. Chem. Res. 2015, 48, 886–896. [Google Scholar] [CrossRef] [PubMed]
  3. Woodward, R.B.; Hoffmann, R. The Conservation of Orbital Symmetry. Angew. Chem. Int. Ed. 1969, 8, 781–853. [Google Scholar] [CrossRef]
  4. Abu-Youssef, M.A.M.; Langer, V.; Barakat, A.; Haukka, M.; Soliman, S.M. Molecular, Supramolecular Structures Combined with Hirshfeld and DFT Studies of Centrosymmetric M(II)-azido {M = Ni(II), Fe(II) or Zn(II)} Complexes of 4-Benzoylpyridine. Symmetry 2021, 13, 2026. [Google Scholar] [CrossRef]
  5. Gouda, M.; Ferjani, H.; El-Lateef, H.M.A.; Khalaf, M.M.; Shaaban, S.; Yousef, T.A. A Competition between Hydrogen, Stacking, and Halogen Bonding in N-(4-((3-Methyl-1,4-dioxo-1,4-dihydronaphthalen-2-yl)selanyl)phenyl)acetamide: Structure, Hirshfeld Surface Analysis, 3D Energy Framework Approach, and DFT Calculation. Int. J. Mol. Sci. 2022, 23, 2716. [Google Scholar] [CrossRef]
  6. Kolomeychuk, F.M.; Safonova, E.A.; Polovkova, M.A.; Sinelshchikova, A.A.; Martynov, A.G.; Shokurov, A.V.; Kirakosyan, G.A.; Efimov, N.N.; Tsivadze, A.Y.; Gorbunova, Y.G. Switchable Aromaticity of Phthalocyanine via Reversible Nucleophilic Aromatic Addition to an Electron-Deficient Phosphorus(V) Complex. J. Am. Chem. Soc. 2021, 143, 14053–14058. [Google Scholar] [CrossRef]
  7. Weinhold, F.; Landis, C.R. Discovering Chemistry with Natural Bond Orbitals; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2012; ISBN 9781118229101. [Google Scholar]
  8. Lepetit, C.; Fau, P.; Fajerwerg, K.; Kahn, M.L.; Silvi, B. Topological analysis of the metal-metal bond: A tutorial review. Coord. Chem. Rev. 2017, 345, 150–181. [Google Scholar] [CrossRef] [Green Version]
  9. Cortes-Guzman, F.; Bader, R. Complementarity of QTAIM and MO theory in the study of bonding in donor-acceptor complexes. Coord. Chem. Rev. 2005, 249, 633–662. [Google Scholar] [CrossRef]
  10. Kumar, P.S.V.; Raghavendra, V.; Subramanian, V. Bader’s Theory of Atoms in Molecules (AIM) and its Applications to Chemical Bonding. J. Chem. Sci. 2016, 128, 1527–1536. [Google Scholar] [CrossRef]
  11. Bader, R.F.W.; Legare, D.A. Properties of atoms in molecules: Structures and reactivities of boranes and carboranes. Can. J. Chem. 1992, 70, 657–676. [Google Scholar] [CrossRef]
  12. Bader, R.F.W. Atoms in Molecules: A Quantum Theory; Oxford University Press: Oxford, UK, 1990. [Google Scholar]
  13. Savin, P.-D.A.; Nesper, R.; Wengert, S.; Fässler, T.F. ELF: The Electron Localization Function. Angew. Chem. Int. Ed. 1997, 36, 1808–1832. [Google Scholar] [CrossRef]
  14. Morgante, P.; Peverati, R. The devil in the details: A tutorial review on some undervalued aspects of density functional theory calculations. Int. J. Quant. Chem. 2020, 120, e26332. [Google Scholar] [CrossRef]
  15. Shen, Y.-F.; Xu, C.; Cheng, L.-J. Deciphering chemical bonding in BnHn2−(n = 2–17): Flexible multicenter bonding. RSC Adv. 2017, 7, 36755–36764. [Google Scholar] [CrossRef] [Green Version]
  16. Teixidor, C.V. The uniqueness of boron as a novel challenging element for drugs in pharmacology, medicine and for smart biomaterials. Future Med. Chem. 2013, 5, 617–619. [Google Scholar] [CrossRef] [PubMed]
  17. Evamarie Hey-Hawkins, C.V.T. Boron-Based Compounds: Potential and Emerging Applications in Medicine; John Wiley & Sons Ltd.: Hoboken, NJ, USA, 2018; ISBN 978-1-119-27555-8. [Google Scholar]
  18. Poater, J.; Solà, M.; Viñas, C.; Teixidor, F. Hückel’s Rule of Aromaticity Categorizes Aromatic closo Boron Hydride Clusters. Chem. A Eur. J. 2016, 22, 7437–7443. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Duchêne, L.; Remhof, A.; Hagemann, H.; Battaglia, C. Status and prospects of hydroborate electrolytes for all-solid-state batteries. Energy Stor. Mater. 2020, 25, 782–794. [Google Scholar] [CrossRef]
  20. Kochnev, V.K.; Avdeeva, V.V.; Goeva, L.V.; Malinina, E.A.; Kuznetsov, N.T. The undecahydrodecaborate anion B10H11—As the starting reagent in exopolyhedral substitution and complexation: Theoretical and experimental prerequisites. Russ. J. Inorg. Chem. 2012, 57, 331–336. [Google Scholar] [CrossRef]
  21. Huang, Z.; Wang, S.; Dewhurst, R.D.; Ignat’Ev, N.V.; Finze, M.; Braunschweig, H. Boron: Its Role in Energy—Related Processes and Applications. Angew. Chem. Int. Ed. 2019, 59, 8800–8816. [Google Scholar] [CrossRef]
  22. Ali, F.; Hosmane, N.S.; Zhu, Y. Boron Chemistry for Medical Applications. Molecules 2020, 25, 828. [Google Scholar] [CrossRef] [Green Version]
  23. Hu, K.; Yang, Z.; Zhang, L.; Xie, L.; Wang, L.; Xu, H.; Josephson, L.; Liang, S.H.; Zhang, M.-R. Boron agents for neutron capture therapy. Coord. Chem. Rev. 2020, 405, 213139. [Google Scholar] [CrossRef]
  24. Vorobyeva, M.; Dymova, M.; Novopashina, D.; Kuligina, E.; Timoshenko, V.; Kolesnikov, I.; Taskaev, S.; Richter, V.; Venyaminova, A. Tumor Cell-Specific 2′-Fluoro RNA Aptamer Conjugated with Closo-Dodecaborate as A Potential Agent for Boron Neutron Capture Therapy. Int. J. Mol. Sci. 2021, 22, 7326. [Google Scholar] [CrossRef] [PubMed]
  25. Lipscomb, W.N. Framework Rearrangement in Boranes and Carboranes. Science 1966, 153, 373–378. [Google Scholar] [CrossRef] [PubMed]
  26. Dobrott, R.D.; Lipscomb, W.N. Structure of Cu2B10H10. J. Chem. Phys. 1962, 37, 1779–1784. [Google Scholar] [CrossRef]
  27. Schwalbe, C.H.; Lipscomb, W.N. Structures of B20H182- and B20H18NO3—And conformations of the triethylammonium ion. J. Am. Chem. Soc. 1969, 91, 194–196. [Google Scholar] [CrossRef]
  28. Lipscomb, W.N. The Boranes and Their Relatives. Science 1977, 196, 1047–1055. [Google Scholar] [CrossRef] [Green Version]
  29. Hoffmann, R.; Lipscomb, W.N. Theory of Polyhedral Molecules. I. Physical Factorizations of the Secular Equation. J. Chem. Phys. 1962, 36, 2179–2189. [Google Scholar] [CrossRef]
  30. Hoffmann, R.; Lipscomb, W.N. Theory of Polyhedral Molecules. III. Population Analyses and Reactivities for the Carboranes. J. Chem. Phys. 1962, 36, 3489–3493. [Google Scholar] [CrossRef]
  31. Hoffmann, R.; Lipscomb, W.N. Boron Hydrides: LCAO—MO and Resonance Studies. J. Chem. Phys. 1962, 37, 2872–2883. [Google Scholar] [CrossRef]
  32. Wade, K. Structural and Bonding Patterns in Cluster Chemistry. Adv. Inorg. Chem. Radiochem. 1976, 18, 1–66. [Google Scholar] [CrossRef]
  33. Fox, M.A.; Nervi, C.; Crivello, A.; Batsanov, A.S.; Howard, J.A.K.; Wade, K.; Low, P.J. Structural, spectroscopic, electrochemical and computational studies of C,C′-diaryl-ortho-carboranes, 1-(4-XC6H4)-2-Ph-1,2-C2B10H10 (X = H, F, OMe, NMe2, NH2, OH and O−). J. Solid State Electrochem. 2009, 13, 1483–1495. [Google Scholar] [CrossRef] [Green Version]
  34. Nelyubin, A.V.; Selivanov, N.A.; Bykov, A.Y.; Klyukin, I.N.; Novikov, A.S.; Zhdanov, A.P.; Karpechenko, N.Y.; Grigoriev, M.S.; Zhizhin, K.Y.; Kuznetsov, N.T. Primary Amine Nucleophilic Addition to Nitrilium Closo-Dodecaborate [B12H11NCCH3]−: A Simple and Effective Route to the New BNCT Drug Design. Int. J. Mol. Sci. 2021, 22, 13391. [Google Scholar] [CrossRef] [PubMed]
  35. Mingos, D.M.P. Polyhedral skeletal electron pair approach. Acc. Chem. Res. 1984, 17, 311–319. [Google Scholar] [CrossRef]
  36. Copley, R.C.B.; Mingos, D.M.P. The novel structure of the [Au11(PMePh2)10]3+ cation: Crystal structures of [Au11(PMePh2)10][C2B9H12]3·4thf and [Au11(PMePh2)10][C2B9H12]3 (thf = tetrahydrofuran). J. Chem. Soc. Dalt. Trans. 1996, 17, 479–489. [Google Scholar] [CrossRef]
  37. Yan, Y.-K.; Mingos, D.M.P.; Kurmoo, M.; Li, W.-S.; Scowen, I.J.; McPartlin, M.; Coomber, A.T.; Friend, R.H. Synthesis, structure and physical properties of tetrathiafulvalenium salts of the ferracarborane complex commo-[3,3′-Fe{1-(C4H3S)-1,2-C2B9H10}2]−. J. Chem. Soc. Dalt. Trans. 1995, 17, 2851–2860. [Google Scholar] [CrossRef]
  38. Hertler, W.R.; Knoth, W.H.; Muetterties, E.L. Chemistry of Boranes. XXIV. Carbonylation of Derivatives of B10H102- and B12Hl22- with Oxalyl Chloride. Inorg. Chem. 1965, 4, 288–293. [Google Scholar] [CrossRef]
  39. Miller, H.C.; Miller, N.E.; Muetterties, E.L. Chemistry of Boranes. XX. Syntheses of Polyhedral Boranes. Inorg. Chem. 1964, 3, 1456–1463. [Google Scholar] [CrossRef]
  40. Chamberland, B.L.; Muetterties, E.L. Chemistry of Boranes. XVIII. Oxidation of B10H10-2 and Its Derivatives. Inorg. Chem. 1964, 231, 1450–1456. [Google Scholar] [CrossRef]
  41. Miller, H.C.; Hertler, W.R.; Muetterties, E.L.; Knoth, W.H.; Miller, N.E. Chemistry of Boranes. XXV. Synthesis andc Chemistry of Base Derivatives of B10H102- and B12Hl22-. Inorg. Chem. 1965, 4, 1216–1221. [Google Scholar] [CrossRef]
  42. Golub, I.; Filippov, O.; Belkova, N.; Epstein, L.; Shubina, E. The Reaction of Hydrogen Halides with Tetrahydroborate Anion and Hexahydro-closo-hexaborate Dianion. Molecules 2021, 26, 3754. [Google Scholar] [CrossRef]
  43. Černý, R.; Brighi, M.; Murgia, F. The Crystal Chemistry of Inorganic Hydroborates. Chemistry 2020, 2, 805–826. [Google Scholar] [CrossRef]
  44. Zhao, X.; Yang, Z.; Chen, H.; Wang, Z.; Zhou, X.; Zhang, H. Progress in three-dimensional aromatic-like closo-dodecaborate. Coord. Chem. Rev. 2021, 444, 214042. [Google Scholar] [CrossRef]
  45. Kamin, A.A.; Juhasz, M.A. Exhaustive Cyanation of the Dodecaborate Dianion: Synthesis, Characterization, and X-ray Crystal Structure of [B12(CN)12]2−. Inorg. Chem. 2019, 59, 189–192. [Google Scholar] [CrossRef]
  46. Lacroix, M.R.; Liu, Y.; Strauss, S.H. Hydrated Metal Ion Salts of the Weakly Coordinating Fluoroanions PF6−, TiF62−, B12F122−, Ga(C2F5)4−, B(3,5-C6H3(CF3)2)4−, and Al(OC(CF3)3)4−. In Search of the Weakest HOH⋯F Hydrogen Bonds. Inorg. Chem. 2019, 58, 14900–14911. [Google Scholar] [CrossRef] [PubMed]
  47. Klyukin, I.N.; Novikov, A.S.; Zhdanov, A.P.; Zhizhin, K.Y.; Kuznetsov, N.T. QTAIM Analysis of Mono-Hydroxy Derivatives of closo-Borate Anions [BnHn–1OH]2 (n = 6, 10, 12). Russ. J. Inorg. Chem. 2019, 64, 1825–1828. [Google Scholar] [CrossRef]
  48. Kononova, E.; Klemenkova, Z. The electronic structure of nido-B10H14 and [6-Ph-nido-6-CB9H11] in terms of Bader’s theory (AIM). J. Mol. Struct. 2013, 1036, 311–317. [Google Scholar] [CrossRef]
  49. Mindich, A.L.; Bokach, N.A.; Dolgushin, F.M.; Haukka, M.; Lisitsyn, L.A.; Zhdanov, A.P.; Zhizhin, K.Y.; Miltsov, S.A.; Kuznetsov, N.T.; Kukushkin, V.Y. 1,3-Dipolar Cycloaddition of Nitrones to a Nitrile Functionality in closo-Decaborate Clusters: A Novel Reactivity Mode for the Borylated C≡N Group. Organometallics 2012, 31, 1716–1724. [Google Scholar] [CrossRef]
  50. Burianova, V.K.; Mikherdov, A.S.; Bolotin, D.S.; Novikov, A.S.; Mokolokolo, P.P.; Roodt, A.; Boyarskiy, V.P.; Suslonov, V.V.; Zhdanov, A.P.; Zhizhin, K.Y.; et al. Electrophilicity of aliphatic nitrilium closo-decaborate clusters: Hyperconjugation provides an unexpected inverse reactivity order. J. Organomet. Chem. 2018, 870, 97–103. [Google Scholar] [CrossRef]
  51. Bolotin, D.S.; Burianova, V.K.; Novikov, A.S.; Demakova, M.Y.; Pretorius, C.; Mokolokolo, P.P.; Roodt, A.; Bokach, N.A.; Suslonov, V.V.; Zhdanov, A.P.; et al. Nucleophilicity of Oximes Based upon Addition to a Nitrilium closo-Decaborate Cluster. Organometallics 2016, 35, 3612–3623. [Google Scholar] [CrossRef]
  52. Stogniy, M.Y.; Erokhina, S.A.; Sivaev, I.B.; Bregadze, V.I. Nitrilium derivatives of polyhedral boron compounds (boranes, carboranes, metallocarboranes): Synthesis and reactivity. Phosph. Sulf. Silicon Relat. Elem. 2019, 194, 983–988. [Google Scholar] [CrossRef]
  53. Golub, I.E.; Filippov, O.A.; Belkova, N.V.; Epstein, L.M.; Shubina, E.S. The Mechanism of Halogenation of Decahydro-closo-Decaborate Dianion by Hydrogen Chloride. Russ. J. Inorg. Chem. 2021, 66, 1639–1648. [Google Scholar] [CrossRef]
  54. Shore, S.G.; Hamilton, E.J.M.; Bridges, A.N.; Bausch, J.; Krause-Bauer, J.A.; Dou, D.; Liu, J.; Liu, S.; Bin, D.; Hall, H.; et al. The Solid State Structure of [B10H11]− and Its Dynamic NMR Spectra in Solution. Inorg. Chem. 2003, 42, 1175–1186. [Google Scholar] [CrossRef] [PubMed]
  55. Kochnev, V.K.; Avdeeva, V.V.; Malinina, E.A.; Kuznetsov, N.T. Theoretical study of H2 elimination from [BnHn+1] monoanions (n = 6–9, 11). Russ. J. Inorg. Chem. 2014, 59, 1268–1275. [Google Scholar] [CrossRef]
  56. Sun, Y.; Zhang, J.; Zhang, Y.; Liu, J.; Van Der Veen, S.; Duttwyler, S. The closo-Dodecaborate Dianion Fused with Oxazoles Provides 3D Diboraheterocycles with Selective Antimicrobial Activity. Chem. A Eur. J. 2018, 24, 10364–10371. [Google Scholar] [CrossRef] [PubMed]
  57. Cao, K.; Zhang, C.-Y.; Xu, T.-T.; Wu, J.; Wen, X.-Y.; Jiang, W.-J.; Chen, M.; Yang, J. Synthesis of Polyhedral Borane Cluster Fused Heterocycles via Transition Metal Catalyzed B-H Activation. Molecules 2020, 25, 391. [Google Scholar] [CrossRef] [Green Version]
  58. Zhu, T.-C.; Xing, Y.-Y.; Sun, Y.; Duttwyler, S.; Hong, X. Directed B–H functionalization of the closo-dodecaborate cluster via concerted iodination—deprotonation: Reaction mechanism and origins of regioselectivity. Org. Chem. Front. 2020, 7, 3648–3655. [Google Scholar] [CrossRef]
  59. Voinova, V.V.; Selivanov, N.A.; Plyushchenko, I.V.; Vokuev, M.F.; Bykov, A.Y.; Klyukin, I.N.; Novikov, A.S.; Zhdanov, A.P.; Grigoriev, M.S.; Rodin, I.A.; et al. Anion—Novel Members of Diboroheterocycle Class. Molecules 2021, 26, 248. [Google Scholar] [CrossRef]
  60. Klyukin, I.; Zhdanov, A.; Matveev, E.; Razgonyaeva, G.; Grigoriev, M.; Zhizhin, K.; Kuznetsov, N. Synthesis and reactivity of closo -decaborate anion derivatives with multiple carbon-oxygen bonds. Inorg. Chem. Commun. 2014, 50, 28–30. [Google Scholar] [CrossRef]
  61. Laskova, J.; Ananiev, I.; Kosenko, I.; Serdyukov, A.; Stogniy, M.Y.; Sivaev, I.B.; Grin, M.A.; Semioshkin, A.; Bregadze, V.I. Nucleophilic addition reactions to nitrilium derivatives [B12H11NCCH3]− and [B12H11NCCH2CH3]. Synthesis and structures of closo-dodecaborate-based iminols, amides and amidines. Dalt. Trans. 2022, 51, 3051–3059. [Google Scholar] [CrossRef]
  62. Avdeeva, V.V.; Polyakova, I.N.; Goeva, L.V.; Malinina, E.A.; Kuznetsov, N.T. [2,6(9)-B10H8(O)2CCH3]− and [2,7(8)- B10H8 (OC(O)CH3)2]2− derivatives in synthesis of position isomers of the [B10H8(OC(O)CH3)(OH)]2− anion with the 2,6(9)- and 2,7(8)-arrangement of functional groups. Russ. J. Inorg. Chem. 2014, 59, 1247–1258. [Google Scholar] [CrossRef]
  63. Klyukin, I.; Kubasov, A.; Limarev, I.; Zhdanov, A.; Matveev, E.; Polyakova, I.; Zhizhin, K.; Kuznetsov, N. The new approach to formation of exo boron-oxygen bonds from the decahydro-closo-decaborate(2-) anion. Polyhedron 2015, 101, 215–222. [Google Scholar] [CrossRef]
  64. Oller, J.; Pérez, P.; Ayers, P.W.; Vöhringer-Martinez, E. Global and local reactivity descriptors based on quadratic and linear energy models forα,β-unsaturated organic compounds. Int. J. Quant. Chem. 2018, 118, 1–13. [Google Scholar] [CrossRef]
  65. Gómez, T.; Fuentealba, P.; Robles-Navarro, A.; Cárdenas, C. Links among the Fukui potential, the alchemical hardness and the local hardness of an atom in a molecule. J. Comput. Chem. 2021, 42, 1681–1688. [Google Scholar] [CrossRef] [PubMed]
  66. Khanum, G.; Ali, A.; Shabbir, S.; Fatima, A.; Alsaiari, N.; Fatima, Y.; Ahmad, M.; Siddiqui, N.; Javed, S.; Gupta, M. Vibrational Spectroscopy, Quantum Computational and Molecular Docking Studies on 2-[(1H-Benzimidazol-1-yl)-methyl]benzoic Acid. Crystals 2022, 12, 337. [Google Scholar] [CrossRef]
  67. Klyukin, I.N.; Vlasova, Y.S.; Novikov, A.S.; Zhdanov, A.P.; Hagemann, H.R.; Zhizhin, K.Y.; Kuznetsov, N.T. B-F bonding and reactivity analysis of mono- and perfluoro-substituted derivatives of closo-borate anions (6, 10, 12): A computational study. Polyhedron 2022, 211, 115559. [Google Scholar] [CrossRef]
  68. Klyukin, I.; Vlasova, Y.; Novikov, A.; Zhdanov, A.; Zhizhin, K.; Kuznetsov, N. Theoretical Study of closo-Borate Anions [BnHn]2− (n = 5–12): Bonding, Atomic Charges, and Reactivity Analysis. Symmetry 2021, 13, 464. [Google Scholar] [CrossRef]
  69. Hehre, W.J.; Ditchfield, R.; Pople, J.A. Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules. J. Chem. Phys. 1972, 56, 2257–2261. [Google Scholar] [CrossRef]
  70. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132, 154104–154119. [Google Scholar] [CrossRef] [Green Version]
  71. Neese, F. The ORCA program system. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2012, 2, 73–78. [Google Scholar] [CrossRef]
  72. Gonzalez, C.; Schlegel, H.B. Improved algorithms for reaction path following: Higher-order implicit algorithms. J. Chem. Phys. 1991, 95, 5853–5860. [Google Scholar] [CrossRef] [Green Version]
  73. Gonzalez, C.; Schlegel, H.B. An improved algorithm for reaction path following. J. Chem. Phys. 1989, 90, 2154–2161. [Google Scholar] [CrossRef]
  74. Gonzalez, C.; Schlegel, H.B. Reaction path following in mass-weighted internal coordinates. J. Phys. Chem. 1990, 94, 5523–5527. [Google Scholar] [CrossRef]
  75. Marenich, A.V.; Cramer, C.J.; Truhlar, D.G. Universal Solvation Model Based on Solute Electron Density and on a Continuum Model of the Solvent Defined by the Bulk Dielectric Constant and Atomic Surface Tensions. J. Phys. Chem. B 2009, 113, 6378–6396. [Google Scholar] [CrossRef] [PubMed]
  76. Glendening, E.D.; Landis, C.R.; Weinhold, F. NBO 7.0: New vistas in localized and delocalized chemical bonding theory. J. Comput. Chem. 2019, 40, 2234–2241. [Google Scholar] [CrossRef] [PubMed]
  77. Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2011, 33, 580–592. [Google Scholar] [CrossRef] [PubMed]
  78. Adrienko, G.A. Chemcraft, Version 1.8 (Build 530). 2017. Available online: http://www.chemcraftprog.com (accessed on 9 April 2022).
Ijms 23 04190 sch001 550
Scheme 1. Protonation of carboxonium derivative [2,6-B10H8O2CCH3]− anion.
Scheme 1. Protonation of carboxonium derivative [2,6-B10H8O2CCH3]− anion.
Ijms 23 04190 sch001
Ijms 23 04190 g001 550
Figure 1. Optimized structures of [2,6-B10H8O2CCH3]− and [2,6-B10H8O2CCH3*Hfac]0.
Figure 1. Optimized structures of [2,6-B10H8O2CCH3]− and [2,6-B10H8O2CCH3*Hfac]0.
Ijms 23 04190 g001
Ijms 23 04190 g002 550
Figure 2. Molecular graph showing the results of the topological analysis of the electron density distribution in the model structure of the [2,6-B10H8O2CCH3*CF3SO3H]− (Comp).
Figure 2. Molecular graph showing the results of the topological analysis of the electron density distribution in the model structure of the [2,6-B10H8O2CCH3*CF3SO3H]− (Comp).
Ijms 23 04190 g002
Ijms 23 04190 g003 550
Figure 3. Molecular graph showing the results of the topological analysis of the electron density distribution in the model structure of TS. Atomic designations are given in Figure 2.
Figure 3. Molecular graph showing the results of the topological analysis of the electron density distribution in the model structure of TS. Atomic designations are given in Figure 2.
Ijms 23 04190 g003
Ijms 23 04190 g004 550
Figure 4. Energy profile of protonation of [2,6-B10H8O2CCH3]−. TfOH–CF3SO3H, TfO−–CF3SO3−. The deciphering of other abbreviations is given in the main text.
Figure 4. Energy profile of protonation of [2,6-B10H8O2CCH3]−. TfOH–CF3SO3H, TfO−–CF3SO3−. The deciphering of other abbreviations is given in the main text.
Ijms 23 04190 g004
Ijms 23 04190 g005 550
Figure 5. 11B-1H NMR spectra of [B10H8O2CCH3]− (bottom) and [B10H8O2CCH3*Hfac]0 (top). The numbers are the positions of the boron atoms in the cluster polyhedral.
Figure 5. 11B-1H NMR spectra of [B10H8O2CCH3]− (bottom) and [B10H8O2CCH3*Hfac]0 (top). The numbers are the positions of the boron atoms in the cluster polyhedral.
Ijms 23 04190 g005
Ijms 23 04190 g006 550
Figure 6. Optimized structures of the main isomers of [2,6-B10H8O2CCH3*Hfac]0. The atomic designations are given in Figure 1. ΔGgp—relative Gibbs energy of isomer in gas phase. ΔGdcm—relative Gibbs energy of isomer in dichlormethane solution.
Figure 6. Optimized structures of the main isomers of [2,6-B10H8O2CCH3*Hfac]0. The atomic designations are given in Figure 1. ΔGgp—relative Gibbs energy of isomer in gas phase. ΔGdcm—relative Gibbs energy of isomer in dichlormethane solution.
Ijms 23 04190 g006
Ijms 23 04190 g007 550
Figure 7. Molecular graph showing the results of the topological analysis of the electron density distribution in the model structures of [B10H11]− and [2,6-B10H8O2CCH3*Hfac]0. The atomic designations are given in Figure 2.
Figure 7. Molecular graph showing the results of the topological analysis of the electron density distribution in the model structures of [B10H11]− and [2,6-B10H8O2CCH3*Hfac]0. The atomic designations are given in Figure 2.
Ijms 23 04190 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Klyukin, I.N.; Kolbunova, A.V.; Novikov, A.S.; Nelyubin, A.V.; Selivanov, N.A.; Bykov, A.Y.; Klyukina, A.A.; Zhdanov, A.P.; Zhizhin, K.Y.; Kuznetsov, N.T. Protonation of Borylated Carboxonium Derivative [2,6-B10H8O2CCH3]−: Theoretical and Experimental Investigation. Int. J. Mol. Sci. 2022, 23, 4190. https://doi.org/10.3390/ijms23084190

AMA Style

Klyukin IN, Kolbunova AV, Novikov AS, Nelyubin AV, Selivanov NA, Bykov AY, Klyukina AA, Zhdanov AP, Zhizhin KY, Kuznetsov NT. Protonation of Borylated Carboxonium Derivative [2,6-B10H8O2CCH3]−: Theoretical and Experimental Investigation. International Journal of Molecular Sciences. 2022; 23(8):4190. https://doi.org/10.3390/ijms23084190

Chicago/Turabian Style

Klyukin, Ilya N., Anastasia V. Kolbunova, Alexander S. Novikov, Aleksey V. Nelyubin, Nikita A. Selivanov, Alexander Yu. Bykov, Alexandra A. Klyukina, Andrey P. Zhdanov, Konstantin Yu. Zhizhin, and Nikolay T. Kuznetsov. 2022. "Protonation of Borylated Carboxonium Derivative [2,6-B10H8O2CCH3]−: Theoretical and Experimental Investigation" International Journal of Molecular Sciences 23, no. 8: 4190. https://doi.org/10.3390/ijms23084190

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop