Search

The structure, stability, and bonding characteristics of dimers and trimers involving BX 3 and AlX 3 (X = H, F, Cl) in the gas phase, many of them explored for the first time, were investigated using different DFT (B3LYP, B3LYP/D3BJ, and M06-2X) and ab initio (MP2 and G4) methods together with different energy decomposition formalisms, namely, many-body interaction-energy and localized molecular orbital energy decomposition analysis. The electron density of the clusters investigated was analyzed with QTAIM, electron localization function, NCIPLOT, and adaptive natural density partitioning approaches. Our results for triel hydride dimers and Al 2 X 6 (X = F, Cl) clusters are in good agreement with previous studies in the literature, but in contrast with the general accepted idea that B 2 F 6 and B 2 Cl 6 do not exist, we have found that they are predicted to be weakly bound systems if dispersion interactions are conveniently accounted for in the theoretical schemes used. Dispersion interactions are also dominant in both homo- and heterotrimers involving boron halide monomers. Surprisingly, B 3 F 9 and B 3 Cl 9 C 3 v cyclic trimers, in spite of exhibiting rather strong B-X (X = F, Cl) interactions, were found to be unstable with respect to the isolated monomers due to the high energetic cost of the rehybridization of the B atom, which is larger than the two- and three-body stabilization contributions when the cyclic is formed. Another important feature is the enhanced stability of both homo- and heterotrimers in which Al is the central atom because Al is systematically pentacoordinated, whereas this is not the case when the central atom is B, which is only tri- or tetra-coordinated.

Considering that a molecular-level understanding of the azeotropic ethanol-water system can contribute to the search of new methodologies and/or modifications of industrial separation methods, this study tries to provide some clues to understand why azeotropes should be expected for ethanol, but not for methanol. Our exploration of the potential energy surface of (ethanol) 6 -water heteroheptamers, carried out at the B3LYP-D3/6-311++G(d,p) level, shows these heteroclusters to exhibit a cyclic structure where the cooperativity effects between the OH···O HBs is a fundamental ingredient. An analysis of this cooperativity clearly indicates that ethanol-water systems will exhibit a similarly high stability as the heterocluster size approaches the azeotrope. However, a similar behavior should not be expected for the methanol-containing analogues. A comparison between (ethanol) 7 , (ethanol) 6 -water, (methanol) 7 , and (methanol) 6 -water shows the ethanol-containing systems to be significantly more stable than the methanol-containing analogues. This result is probably due to the fact that the OH···O HBs are weaker than those found between ethanol molecules. However, our atoms in molecule (AIM) and noncovalent interaction (NCI) analyses unambiguously show that important contributors to the enhanced stability of the ethanol-containing clusters are the secondary van der Waals interactions between ethyl groups, which are not observed between methyl groups. Hence, while the formation of stable azeotropes is expected for the case of ethanol, for the methanol-containing analogues, the relative stability of the clusters is significantly smaller, and its formation is accompanied by an increase of the free energy.

Pnicogen bonds, which are weak noncovalent interactions (NCIs), can be significantly modified by the presence of beryllium bonds, one of the strongest NCIs known. We demonstrate the importance of this influence by studying ternary complexes in which both NCIs are present, that is, the ternary complexes formed by a nitrogen base (NH 3 , NHCH 2 , and NCH), a phosphine (fluorophosphane, PH 2 F) and a beryllium derivative (BeH 2 , BeF 2 , BeCl 2 , BeCO 3 , and BeSO 4 ). Energies, structures, and nature of the chemical bonding in these complexes are studied by means of ab initio computational methods. The pnicogen bond between the nitrogen base and the phosphine and the beryllium bond between the fluorine atom of fluorophosphane and the beryllium derivative show large cooperativity effects both on energies and geometries, with dissociation energies up to 296 kJ mol -1 and cooperativity up to 104 kJ mol -1 in the most strongly bound complex, CH 2 HN:PH 2 F:BeSO 4 . In the complexes between the strongest nitrogen bases and the strongest beryllium donors, phosphorus-shared and phosphorus-transfer bonds are found.

High-level G4 ab initio calculations allowed us to show that C 4 H 4 (BeX) 4 (X = H, Cl) derivatives behave as rather efficient electron capturers due to their ability to trap the extra electron through the formation of a four-membered beryllium ring. This finding is in agreement with previous work showing the ability of highly electron-deficient atoms, such as beryllium, to lead to multicenter one-electron bonds. In our particular case, the formation of the four-center bond is characterized, in very good harmony, by different topological methods such as quantum theory of atoms in molecules (QTAIM), the electron localization function (ELF), and the noncovalent interactions (NCI) approach and is accompanied by large electron affinity values, around 300 kJ·mol -1 , in the gas phase. Preliminary results may anticipate that the ability of groups of beryllium atoms to trap electrons decays on going to bigger systems.

Weak and strong noncovalent interactions such as tetrel bonds and alkaline-earth bonds, respectively, cooperate and get reinforced when acting together in ternary complexes of general formula RN··· SiH 3 F···MY, where MY is a Be or Mg derivative and RN is a N-containing Lewis base with different hybridization patterns. Cooperativity has been studied in the optimized MP2/aug'-cc-pVTZ ternary complexes by looking at changes on geometries, binding energies, 29 Si NMR chemical shifts, and topological features according to the atoms in molecules theoretical framework. Our study shows that cooperativity in terms of energy is in general significant: more than 40 kJ/mol, and up to 83.6 kJ/mol in the most favorable case. The weakest the isolated interaction, the strongest the reinforcement in the ternary complex; in this sense, the tetrel bond is shortened enormously, between 0.3 and 0.6 Å. This dramatic reinforcement of the tetrel bond is also nicely reflected in the positive variations of the 29 Si chemical shifts in all the ternary complexes. At the same time the ternary complexes are characterized by the presence of totally planar silyl group, due to the pentacoordination of the Si atom. Both the hybridization of the N base and the geometry imposed by the alkaline-earth ligands have a strong influence on the binding energies, as they modify the donor ability of N and the Lewis acid character of the alkaline-earth metal.

The association of Ca 2+ to serine and the subsequent gas-phase unimolecular reactivity of the [Ca(Ser)] 2+ (Ser = Serine) adduct was investigated throughout the use of tandem mass spectrometry techniques and B3LYP/6-311+G(3df,2p)//B3LYP/6-311+G(d,p) density functional theory calculations. In a first step, the structure and relative stability of all possible conformers of serine were obtained and analyzed, as well as the most stable [serine-Ca] 2+ adducts. For the analysis of the different potential energy surfaces associated with the gas-phase unimolecular reactivity of these adducts, only those that differ by less than 100 kJ·mol -1 from the global minimum were taken into account. In agreement with previous studies, the serine-Ca 2+ global minimum corresponds to a charge-solvated structure in which Ca is tricoordinated to neutral serine. The major peaks observed in the nanoelectrospray-MS/MS spectrum of [Ca(Ser)] 2+ adduct correspond to both Coulomb explosions, yielding either CaOH + + [C 3 ,H 6 ,N,O 2 ] + or [C 2 ,H 4 ,O,N] + + [Ca(C,H 3 ,O 2 )] + , and to the loss of neutrals, namely, CH 2 O and H 2 O. Our theoretical survey of the energy profile allow us to conclude that, although all the aforementioned fragmentation processes can have their origin at the global minimum, similar fragmentations involving low-lying conformers, both zwitterionic and nonzwitterionic, compete and should be considered to account for the observed reactivity. We have also found that in some specific cases post-transition state dynamics similar to the ones described before in the literature for formamide-Ca 2+ reactions, may also play a role.

Structures of beryllium subphthalocyanines and beryllium subporphyrazines complexes with different substituents are explored for the first time. Their photochemical properties are studied using time-dependent density functional theory calculations and compared to boron-related compounds for which their photochemical activity is already known. These beryllium compounds were found to be thermodynamically stable in a vacuum and present features similar to those of boron-containing analogues, although the nature of bonding between the cation and the macrocycle presents subtle differences. Most important contributions to the main peak in the Q-band region arise from HOMO to LUMO transitions in the case of subphthalocyanines and alkyl subporphyrazine complexes, whereas a mixture of that contribution and a HOMO-2 to LUMO contribution are present in the case of thioalkyl subporphyrazines. The absorption in the visible region could make these candidates suitable for photochemical devices if combined with appropriate donor groups.

By combining chemical dynamics simulations and RRKM statistical theory we have characterized collision induced dissociation (CID) mechanisms of [M(formamide)](2+) ions (M = Ca, Sr) at different timescales, from few femtoseconds to microseconds. Chemical dynamics simulations account for the short-time and dynamically driven reactivity, such as impulsive collision mechanism for formamide neutral loss. From the simulations, we also got the amounts of energy transferred during the collision and, especially important, the vibrational and rotational energy distributions of the ions that did not react during the simulation time length of 2.5 ps. These internal energy distributions were in turn used in combination with RRKM theory to estimate the rate constants of the possible reactive pathways. Hence, we performed a statistical analysis of the CID dynamics accounting for the long-time and statistical reactivity (i.e., through an IVR mechanism). This multiscale approach allowed us to account for all the products observed in the CID experimental spectra of [Ca(formamide)](2+) and [Sr(formamide)](2+) doubly charged cations, as well as the differences between them.

A significant acidity enhancement and changes on aromaticity were previously observed in squaric acid and its derivatives when beryllium bonds are present in those systems. In order to know if these changes on the chemical properties could be considered a general behavior of carboxylic acids upon complexation with beryllium compounds, complexes between a set of representative carboxylic acids RCOOH (formic acid, acetic acid, propanoic acid, benzoic acid, and oxalic acid) and beryllium compounds BeX2 (X = H, F, Cl) were studied by means of density functional theory calculations. Complexes that contain a dihydrogen bond or a OH···X interaction are the most stable in comparison with other possible BeX2 complexation patterns in which no other weak interactions are involved apart from the beryllium bond. Formic, acetic, propanoic, benzoic, and oxalic acid complexes with BeX2 are much stronger acids than their related free forms. The analysis of the topology of the electron density helps to clarify the reasons behind this acidity enhancement. Importantly, when the halogen atom is replaced by hydrogen in the beryllium compound, the dihydrogen bond complex spontaneously generates a new neutral complex [RCOO:BeH] in which a hydrogen molecule is lost. This seems to be a trend for carboxylic acids on complexing BeX2 compounds.

The mutual interaction between beryllium bonds and halogen bonds within H2Be···FCl···Base complexes, where Base includes a wide set of N- and O-containing Lewis bases, has been studied at the M06-2X/6-31+G(d,p) level of theory. The reliability of this theoretical model was assessed by comparison with ab initio CCSD/aug-cc-pVTZ reference calculations. Cooperative effects were investigated within the framework of the atoms in molecules theory (AIM) by analyzing the topology of the electron density and the changes in the atomic energy components. The decomposition of the total stabilization energy into atomic components is found to be a very reliable tool to describe halogen bond interactions. Both the topological analysis of the electron density and the changes in the atomic energy components of the binding energy show the existence of strong cooperative effects between beryllium and halogen bonds, which are in some cases very intense. In general, there is a correlation between the intrinsic basicity of the Lewis base participating in the halogen bond and the resulting cooperativity in the sense that the stronger the base, the larger the cooperative effects.

The complexes formed by the interaction between a series of phosphines R-PH(2) (R = H, CH(3), c-C(3)H(5), C(6)H(5)) and AlH(3) have been investigated through the use of high-level G4 ab initio calculations. These very stable complexes behave as much stronger acids than the isolated phosphines. This dramatic acidity enhancement, which can be as high as 174 kJ mol(-1), results from a much greater stabilization of the anionic deprotonated species with respect to the neutral one, upon AlH(3) association. This effect depends quantitatively on the nature of the substituent R and is smaller for R = C(6)H(5) because of the conjugation of the P lone pair with the aromatic system. More unexpectedly, however, the phosphine-alane complexes, RPH(2):AlH(3), are more acidic than the corresponding phosphine-borane RPH(2):BH(3) analogues. This unexpected result is due to the enhanced stability of the anionic deprotonated species for complexes involving AlH(3), because the delocalization of the newly created P lone pair with the P-Al bonding density is more favorable when the Lewis acid is aluminum trihydride than when it is borane.

Ab initio equation-of-motion coupled cluster (EOM-CCSD) one-bond spin-spin coupling constants (1)J(B-N), (1)J(B-H), and (1)J(B-F) have been evaluated for complexes X:BH(n)F(3-n) with X = N(2), NCH, NCLi, H(2)CNH, NF(3), and NH(3), for n = 0-3. These complexes can be classified as either covalent or van der Waals complexes, on the basis of their binding energies and B-N distances. (1)J(B-N) for covalent complexes varies significantly from -19 to +9 Hz, whereas (1)J(B-N) is less than 2 Hz for van der Waals complexes. An absolute value of (1)J(B-N) of 3 Hz or greater indicates that the complex is covalently bonded, but a small value of this coupling constant does not necessarily mean that it is a van der Waals complex, in view of the variation among these complexes found for (1)J(B-N) as a function of the B-N distance. Deformation of the boron acid upon complex formation and electron donation by the nitrogen base has opposing effects on both (1)J(B-H) and (1)J(B-F). These effects are relatively small in van der Waals complexes. In covalent complexes, electron donation has the dominant effect on (1)J(B-H), and on (1)J(B-F) in complexes with BH(2)F and BHF(2), but acid deformation has the dominant effect on (1)J(B-F) in complexes with BF(3). Values of both (1)J(B-H) and (1)J(B-F) reflect the van der Waals or covalent nature of the B-N bond.

Cyanoethenethiol (NC-CH=CH-SH) was obtained in a 8:1 Z:E ratio by flash vacuum thermolysis of the t-butylsulfide derivative. Density functional theory (DFT) and G3 ab initio calculations predict the existence of Z-and E-isomers, each of which exhibits two rotamers as a function of the relative position of the SH group. All these rotameric forms are planar (C(s) symmetry) and correspond to synperiplanar and antiperiplanar conformations. Calculations indicate that the synperiplanar Z-isomer is the more stable. In pure form, the cyanoethenethiol rapidly decomposes at room temperature, even at low pressure and partially condenses on the wall of the cell. To record its spectrum, a long optical path of 136 m was necessary, and several successive fillings of the cell were required. On the basis of the calculated harmonic and anharmonic vibrational frequencies, a complete and unambiguous assignment of the experimental spectrum has been carried out.

The interaction of 3,4 dinitrophenol (DNP) with cyclic ketones, lactones, and lactams was investigated by UV-visible spectroscopy and density functional theory (DFT) methods. Equilibrium constants K(HB) for 1:1 hydrogen bonded complexes were determined in solution in CCl(4) and C(6)H(12). For the entire range of studied bases, the pK(HB) scale, varying between 2.92 for beta-propiolactone to 5.53 for 1-methyl-epsilon-caprolactam, indicates that the basicity increases with the ring size. Geometries, energies, and vibrational characteristics of complexes were obtained by means of DFT calculations. For lactones and lactams, the energy difference between the two most stable conformers, cis and trans, with respect to the ring oxygen (nitrogen) atom, is relatively small, suggesting that the complex observed in solution is probably an equilibrium mixture of both forms. The good correlation between Gibbs free energies in solution and in the gas phase, computed at the B3LYP/6-311++G(3df,2p) level of theory, confirms the reliability of our results. The electron density of the complexes has been analyzed by means of the atoms in molecules (AIM) theory and the natural bond orbital (NBO) method have been used to characterize the orbital interactions. Our theoretical survey shows that the 1:1 complexes are stabilized by a network of conventional and/or nonconventional intermolecular hydrogen bonds.

The infrared spectra (3500-500 cm(-1)) of gaseous vinyl selenol and vinyl tellurol have been recorded at 0.1 cm(-1) resolution. For the latter the spectra were obtained at room temperature, but for the former a temperature of -40 degrees C was required because of the chemical instability of vinyl selenol at room temperature. To compensate the very weak vapor pressure of vinyl tellurol at room temperature, a long optical path up to 136 m was necessary to record its spectrum. B3LYP density functional theory (DFT) calculations have been performed to assign the different absorption bands. Since an unambiguous assignment of the absorption bands requires a precise knowledge on the relative abundance of the syn and gauche rotamers of these compounds, their relative energies and their anharmonic vibrational frequencies were obtained using a very extended Def2-QZVP basis set. Two rotamers, the syn, which is planar, and a nonplanar gauche, were found to be local minima for both compounds. The gauche rotamer presents two degenerate conformers, which differ by the position of the SeH (TeH) hydrogen atom above or below the molecular plane. Our theoretical results are in good agreement with the main features of the experimental spectra. Fundamental bands and some combination bands of vinyl selenol and vinyl tellurol were assigned and compared with those of vinyl alcohol and vinyl thiol, whose spectra had been reported previously in the literature.

The reaction of aminoacrylonitrile, a species of astrochemical interest, with Ni(+)((2)D(5/2)) was investigated by means of mass spectrometry techniques and density functional theory calculations. The dominant fragmentations in the MIKE spectrum correspond to the loss of [C2,N,H3], HCN, and NH3, the loss of H2 being very minor. The structure and bonding of the different aminoacrylonitrile-Ni(+) complexes were investigated at the B3LYP/6-311G(d,p) level of theory. The same approach was employed in our survey of the corresponding potential energy surface. This survey indicates that the [C2,N,H3] neutral product can be formed either as ketenimine (CH2CNH) or acetonitrile. The formation of the latter is significantly more exothermic but involves slightly higher activation barriers; so very likely, both isomers are produced along the reaction process. The lost of HNC is not competitive with the loss of HCN, because when the former is formed the products lie higher in energy and the corresponding mechanisms involve energy barriers above the entrance channel. The loss of NH3 is associated with the formation of a complex between cyanoacetylene, HCCCN, which is very abundant in the interstellar media, and Ni(+).

The enthalpies of combustion, heat capacities, enthalpies of sublimation and enthalpies of formation of 2-tert-butylbenzimidazole (2tBuBIM) and 2-phenylimidazole (2PhIM) are reported and the results compared with those of benzene derivatives and a series of azoles (imidazoles, pyrazoles, benzimidazoles and indazoles). Theoretical estimates of the enthalpies of formation were obtained through the use of atom equivalent schemes. The necessary energies were obtained in single-point calculations at the B3LYP/6-311++G(d,p) on B3LYP/6-31G optimized geometries. The comparison of experimental and calculated values of all studied compounds bearing H (unsubstituted), methyl (Me) ethyl (Et), propyl (Pr), isopropyl (iPr), tert-butyl (tBu), benzyl (Bn) and phenyl (Ph) groups show remarkable homogeneity. The remarkable consistency of both the calculated and experimental results allows us to predict with reasonable certainty the missing experimental values. The crystal and molecular structure of the 2-benzylbenzimidazole (2BnBIM) has been determined by X-ray analysis. The observed molecular conformation permits the crystal being built up through N-H...N hydrogen bonds and van der Waals contacts between the molecules. An attempt has been made to relate the crystal structure to the enthalpies of sublimation.

The deprotonation of Cu2+ complexes with uracil, 2-thiouracil, 4-thiouracil, and 2,4-dithiouracil has been investigated by means of B3LYP/ 6-311+G(2df,2p)//6-31G(d) calculations. The most stable [(uracil-H)Cu]+ and [(thiouracil-H)Cu]+ complexes correspond to bidentate structures in which Cu interacts with the deprotonated ring-nitrogen atom and with the oxygen or the sulfur atom of the adjacent carbonyl or thiocarbonyl group. For 2- and 4-thiouracil derivatives, the structures in which the metal cation interacts with the thiocarbonyl group are clearly favored with respect to those in which Cu interacts with the carbonyl group. This is at variance with what was found to be the most stable structure of the corresponding Cu2+ complexes, where association to the carbonyl oxygen was always preferred over the association to the thiocarbonyl group. The [(uracil-H)Cu]+ and [(thiouracil-H)Cu]+ complexes can be viewed as the result of Cu+ attachment to the uracil-H and thiouracil-H radicals formed by the deprotonation of the corresponding uracil+* and thiouracil+* radical cations. As a matter of fact their relative stability is dictated by the intrinsic stability of the corresponding uracil-H and thiouracil-H radical and by the fact that, in general, the N3-deprotonated site is a better electron donor than the N1. In all complexes, the bonding of Cu both to nitrogen and sulfur and to nitrogen and oxygen has a significantly large covalent character.

The gas-phase basicities of a representative set of hydroxy- and methoxycarbonyl compounds (hydroxyacetone, 1, 3-hydroxybutanone, 2, 3-hydroxy-3-methylbutanone, 3, 1-hydroxy-2-butanone, 4, 4-hydroxy-2-butanone, 5, 5-hydroxy-2-pentanone, 6, methoxyacetone, 7, 3-methoxy-2-butanone, 8, 4-methoxy-2-butanone, 9, and 5-methoxy-2-pentanone, 10) were experimentally determined by the equilibrium method using Fourier transform ion cyclotron resonance and high-pressure mass spectrometry techniques. The latter method allows the measurement of proton transfer equilibrium constants at various temperatures and thus the estimate of both the proton affinities and the protonation entropies of the relevant species. Quantum chemical calculations at the G3 and the B3LYP/6-311+G(3df,2p)//6-31G(d) levels of theory were undertaken in order to find the most stable structures of the neutrals 1-10 and their protonated forms. Conformational and vibrational analyses have been done with the aim of obtaining a theoretical estimate of the protonation entropies.

The structure and relative stability of methanol complexes with various cyclic ketones, lactones, lactams, and N-methyl lactams from three- to seven-membered rings have been investigated using the density functional theory method. The geometries, harmonic frequencies, and energies were calculated at the B3LYP/6-311+G(d,p) level. Three stable structures, cis-a, cis-b, and trans, with respect to the ring oxygen (nitrogen) atom, were found to be local minima of the potential energy surface. For lactones and N-methyl lactams, the most stable structure is trans; it is stabilized, as in cyclic ketones, through the conventional hydrogen bond (HB) interaction between the basic carbonyl oxygen and the acidic methanolic hydrogen and an unconventional HB interaction between the methanolic oxygen and the CH hydrogen, in the alpha position of the carbonyl group. For unsubstituted lactams, the cis-a structure, stabilized through a HB interaction between the NH group and the methanol oxygen in addition to the conventional HB interaction, is the most stable. The topological properties of the electron density ratify the existence of conventional (N,O-H. . .O) and unconventional (C-H. . .O) hydrogen bonding. A good correlation was found between the HB distances and the electron density at the HB critical point. The unsubstituted lactams yield more stable complexes with methanol than N-methyl lactams, lactones, and cyclic ketones. In the most stable complexes, both components behave simultaneously as a HB donor and as a HB acceptor.

The binding energies (BEs) of alkali metal monocations and alkaline-earth metal dications to a series of small oxygen and nitrogen bases have been evaluated by means of CCSD(T) calculations on B3-LYP optimized structures. These calculations were carried out both using all-electron basis sets, and additionally using an effective-core potential (ECP) to describe the inner electrons of the metal. A theoretical model aiming at analyzing the effects on the binding energy trends of electrostatic, polarization, and covalent contributions, as well as geometry distortion, was employed. From this analysis, we conclude that although the neutral-ion interaction energy for alkali and alkaline-earth metal cations is dominated by electrostatic contributions, in many cases the correct basicity trends are only attained once polarization effects are also included in the model. This is indeed the case when Ca2+ and K+ are bound to ammonia and formaldehyde. Geometry distortions triggered by polarization are also necessary, in some cases, to obtain the correct basicity trends. In addition, in particular for alkaline-earth dications, the energy associated with covalent interactions sometimes dictates the basicity trend. Our observations imply that simple models based on ion-dipole interactions, that are frequently used in the literature to explain affinity trends in ion-molecule reactions, are generally not likely to be reliable.

Bond dissociation energies (BDEs) of all possible A-X single bonds involving the first- and second-row atoms, from Li to Cl, where the free valences are saturated by hydrogens, have been estimated through the use of the G3-theory and at the B3LYP/6-311+G(3df,2pd)//B3LYP/6-31G(2df,p) DFT level of theory. BDEs exhibit a periodical behavior. The A-X (A = Li, Be, B, Na, Mg, Al, and Si) BDEs show a steady increase along the first and the second row of the periodic table as a function of the atomic number Z(X). For A-X bonds involving electronegative atoms (A = C, N, O, F, P, S, and Cl) the bond energies achieve a maximum around Z(X) = 5. The same behavior is observed when BDEs are plotted against the electronegativity chi(X) of the atom X. Thus, for A-X bonds (A = Li, Be, B, Na, Mg, Al, Si), the BDEs for a fixed A increases, grosso modo, as the electronegativity differences between X and A increase, with some exceptions, which reflect the differences in the relaxation energies of the radicals produced upon the bond cleavage. A similar trend, albeit less pronounced, is found for single A-X bonds, where A = C, N, O, F, P, S, and Cl. However, there is an additional feature embodied in the enhancement of the strength of the A-boron bonds due to the ability of boron to act as a strong electron acceptor. The trends in bond lengths and charge densities at the bond critical points are in line with the aforementioned behavior.