Discovering trends in the Lewis acidity of beryllium and magnesium hydrides and fluorides with increasing clusters size.

We have reported in the last years the strong effect that Be- and Mg-containing Lewis acids have on the intrinsic properties of typical bases, which become acids upon complexation. In an effort to investigate these changes when the Be and Mg derivatives form clusters of increasing size, we have examined the behavior of the (MX2)n (M = Be, Mg; X = H, F; n = 1, 2, 3) clusters when they interact with ammonia, methanimine, hydrogen cyanide and pyridine, and with their corresponding deprotonated forms. The complexes obtained at the M06-2X/aug-cc-pVTZ level were analyzed using the MBIE energy decomposition formalism, in parallel with QTAIM, ELF, NCIPLOT and AdNDP analyses of their electron density. For n = 1 the interaction enthalpy for the different families of monomers, Be (Mg) hydrides and Be (Mg) fluorides, follows the same trend as the intrinsic basicity of the base that interacts with them. This interaction is greatly reinforced after the deprotonation of the base, resulting in a significant enhancement of the intrinsic acidity of the corresponding MX2-Base complex. For (MX2)2 clusters a further reinforcement of the interaction with the base is observed, this reinforcement being again larger for the deprotonated complexes. However, the concomitant increase of their intrinsic acidity is one order of magnitude larger for hydrides than for fluorides. Unexpectedly, the cyclic conformers (MX2)3, which are more unstable than the linear ones, become the global minima after association with the base and the same is true for the deprotonated complex. Accordingly, a further increase of the intrinsic acidity of the (MX2)3-Base complexes with respect to the (MX2)2-Base ones is observed. This effect is maximum for (MgF2)3 clusters, to the point that the (MgF2)3-Base complexes become more acidic than nitric acid, the extreme case being the cluster (MgF2)3-NCH, whose acidity is higher than that of perchloric acid.

Ionization, intrinsic basicity, and intrinsic acidity of unsaturated diols of astrochemical interest: 1,1- and 1,2-ethenediol: A theoretical survey.

The structure, stability, and bonding characteristics of 1,1- and 1,2-ethenediol, their radical cations, and their protonated and deprotonated species were investigated using high-level ab initio G4 calculations. The electron density of all the neutral and charged systems investigated was analyzed using the QTAIM, ELF, and NBO approaches. The vertical ionization potential (IP) of the five stable tautomers of 1,2-ethenediol and the two stable tautomers of 1,1-ethenediol go from 11.81 to 12.27 eV, whereas the adiabatic ones go from 11.00 to 11.72 eV. The adiabatic ionization leads to a significant charge delocalization along the O-C-C-O skeleton. The most stable protonated form of (Z)-1,2-ethenediol can be reached by the protonation of both the anti-anti and the syn-anti conformers, whereas the most stable deprotonated form arises only from the syn-anti one. Both charged species are extra-stabilized by the formation of an O-H···O intramolecular hydrogen bond (IHB) which is not found in the neutral system. (Z)-1,2-ethenediol is predicted to be less stable, less basic, and more acidic than its cis-glycolaldehyde isomer. The most stable protonated species of (E)-1,2-ethenediol comes from its syn-syn conformer, although the anti-anti conformer is the most basic one. Contrarily, the three conformers yield a common deprotonated species, so their acidity follows exactly their relative stability. Again, the (E)-1,2-ethenediol is predicted to be less stable, less basic, and more acidic than its trans-glycolaldehyde isomer. Neither the neutral nor the protonated or the deprotonated forms of 1,1-ethenediol show the formation of any O-H···O IHB. The most stable protonated species is formed by the protonation of any of the two tautomers, but the most stable deprotonated form arises exclusively from the syn-anti neutral conformer. The conformers of 1,1-ethenediol are much less stable and significantly less basic than their isomer, acetic acid, and only slightly more acidic.

A Holistic View of the Interactions between Electron-Deficient Systems: Clustering of Beryllium and Magnesium Hydrides and Halides.

In the search for common bonding patterns in pure and mixed clusters of beryllium and magnesium derivatives, the most stable dimers and trimers involving BeX2 and MgX2 (X = H, F, Cl) have been studied in the gas phase using B3LYP and M06-2X DFT methods and the G4 ab initio composite procedure. To obtain some insight into their structure, stability, and bonding characteristics, we have used two different energy decomposition formalisms, namely MBIE and LMO-EDA, in parallel with the analysis of the electron density with the help of QTAIM, ELF, NCIPLOT, and AdNDP approaches. Some interesting differences are already observed in the dimers, where the stability sequence observed for the hydrides differs entirely from that of the fluorides and chlorides. Trimers also show some peculiarities associated with the presence of compact trigonal cyclic structures that compete in stability with the more conventional hexagonal and linear forms. As observed for dimers, the stability of the trimers changes significantly from hydrides to fluorides or chlorides. Although some of these clusters were previously explored in the literature, the novelty of this work is to provide a holistic approach to the entire series of compounds by using chemical bonding tools, allowing us to understand the stability trends in detail and providing insights for a significant number of new, unexplored structures.

Dispersion, Rehybridization, and Pentacoordination: Keys to Understand Clustering of Boron and Aluminum Hydrides and Halides.

The structure, stability, and bonding characteristics of dimers and trimers involving BX3 and AlX3 (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 Al2X6 (X = F, Cl) clusters are in good agreement with previous studies in the literature, but in contrast with the general accepted idea that B2F6 and B2Cl6 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, B3F9 and B3Cl9 C3v 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.

Thiol-yne chemistry of diferrocenylacetylene: from synthesis and electrochemistry to theoretical studies.

The thiol-yne coupling chemistry of diferrocenylacetylene (FcCCFc) 1, bearing two electron rich and redox-active ferrocenyl units (Fc = Fe(η5-C5H4)(η5-C5H5)) and an internal triple bond, has been investigated for the first time. In order to determine whether steric limitations might affect hydrothiolation, a model reaction using a functionalized monothiol was tested, namely 2-mercaptoethanol I. The thiol-diferrocenylacetylene reactions were initiated either thermally (in toluene with AIBN) or by UV light irradiation (in THF and in the presence of DMPA as the photoinitiator). The outcomes of these thiol-yne reactions showed a strong dependence on the initiation method used, with the thermally initiated one being the most efficient. These thiol-diferrocenylacetylene reactions mainly afforded the (Z)-stereoisomer of the newly obtained vinyl thioether sulfide FcCHC(Fc)S-(CH2)2OH (2), unlike the more common (E)-vinyl sulfides found in other additions to alkynes. The hydrothiolation of the internal -CC- bond in 1 was successfully extended to dithiol 2,2'-(ethylenedioxy)diethanethiol II, leading to the formation of the (ZZ)-isomer, with four ferrocenyl units, as the major product. According to the electrochemical studies, the new asymmetrical ferrocenyl-vinyl sulfides show iron-iron electronic and electrostatic interactions. Theoretical results for the (Z)-stereoisomer (2) suggest that adiabatic oxidation would lead to the loss of almost one electron on the ferrocenyl subunit closer to the thioether chain. Furthermore, the thiol-yne chemistry of the internal -CC- bond in diferrocenylacetylene has been compared to the external triple bond in ethynylferrocene, the theoretical results of which helped us to rationalize the very different reactivities observed in both metallocenes.

Disrupting bonding in azoles through beryllium bonds: Unexpected coordination patterns and acidity enhancement.

Although triazoles and tetrazole are amphoteric and may behave as weak acids, the latter property can be hugely enhanced by beryllium bonds. To explain this phenomenon, the structure and bonding characteristics of the complexes between triazoles and tetrazoles with one and two molecules of BeF2 have been investigated through the use of high-level G4 ab initio calculations. The formation of the complexes between the N basic sites of the azoles and the Be center of the BeF2 molecule and the (BeF2)2 dimer leads to a significant bonding perturbation of both interacting subunits. The main consequence of these electron density rearrangements is the above-mentioned increase in the intrinsic acidity of the azole subunit, evolving from a typical nitrogen base to a very strong nitrogenous acid. This effect is particularly dramatic when the interaction involves the (BeF2)2 dimer, that is, a Lewis acid much stronger than the monomer. Although the azoles investigated have neighboring N-basic sites, their interaction with the (BeF2)2 dimer yields a monodentate complex. However, the deprotonated species becomes extra-stabilized because a second N-Be bond is formed, leading to a new five-membered ring, with the result that the azole-(BeF2)2 complexes investigated become stronger nitrogenous acids than oxyacids such as perchloric acid.

Enthalpies of Adduct Formation between Boron Trifluoride and Selected Organic Bases in Solution: Toward an Accurate Theoretical Entry to Lewis Basicity.

The Lewis basicity of selected organic bases, modeled by the enthalpies of adduct formation between gaseous BF3 and bases in dichloromethane (DCM) solution, is critically examined. Although experimental enthalpies for a large number of molecules have been reported in the literature, it may be desirable to estimate missing or uncertain data for important Lewis bases. We decided to use high-level ab initio procedures, combined with a polarized continuum solvation model, in which the solvated species were the clusters formed by specific hydrogen bonding of DCM with the Lewis base and the Lewis base/BF3 adduct. This mode of interaction with DCM corresponds to a specific solvation model (SSM). The results essentially showed that the enthalpy of BF3 adduct formation in DCM solution was clearly influenced by specific interactions, with DCM acting as hydrogen-bonding donor (HBD) molecule in two ways: base/DCM and adduct/DCM, confirming that specific solvation is an important contribution to experimentally determined Lewis basicity scales. This analysis allowed us to conclude that there are reasons to suspect some gas-phase values to be in error by more than the stated experimental uncertainty. Some experimental values in DCM solution that were uncertain for identified reasons could be complemented by the computed values.

Perturbating Intramolecular Hydrogen Bonds through Substituent Effects or Non-Covalent Interactions.

An analysis of the effects induced by F, Cl, and Br-substituents at the α-position of both, the hydroxyl or the amino group for a series of amino-alcohols, HOCH2(CH2)nCH2NH2 (n = 0-5) on the strength and characteristics of their OH···N or NH···O intramolecular hydrogen bonds (IMHBs) was carried out through the use of high-level G4 ab initio calculations. For the parent unsubstituted amino-alcohols, it is found that the strength of the OH···N IMHB goes through a maximum for n = 2, as revealed by the use of appropriate isodesmic reactions, natural bond orbital (NBO) analysis and atoms in molecules (AIM), and non-covalent interaction (NCI) procedures. The corresponding infrared (IR) spectra also reflect the same trends. When the α-position to the hydroxyl group is substituted by halogen atoms, the OH···N IMHB significantly reinforces following the trend H < F < Cl < Br. Conversely, when the substitution takes place at the α-position with respect to the amino group, the result is a weakening of the OH···N IMHB. A totally different scenario is found when the amino-alcohols HOCH2(CH2)nCH2NH2 (n = 0-3) interact with BeF2. Although the presence of the beryllium derivative dramatically increases the strength of the IMHBs, the possibility for the beryllium atom to interact simultaneously with the O and the N atoms of the amino-alcohol leads to the global minimum of the potential energy surface, with the result that the IMHBs are replaced by two beryllium bonds.

Spontaneous bond dissociation cascades induced by Ben clusters (n = 2,4).

High-level single and multireference ab initio calculations show that the Be4 cluster behaves as a very efficient Lewis acid when interacting with conventional Lewis bases such as ammonia, water or hydrogen fluoride, to the point that the corresponding acid-base interaction triggers a sequential dissociation of all the bonds of the Lewis base. Notably, this behavior is already found for the simplest beryllium cluster, the Be2 dimer. However, whereas for Be2 the first dissociation process involves a low activation barrier which is above the reactants, for Be4 all the bond dissociation processes involve barriers below the entrance channel leading to a cascade of successive exothermic processes, which end up spontaneously in a global minimum in which the bonding patterns of both the base and the Lewis acid are completely destroyed. Indeed, the global minimum, in all cases, is stabilized by three-center Be-H-Be bonds and covalent interactions between the Be atoms and the basic center of the base, which replace the initial metallic bond stabilizing the Be4 cluster. As a consequence, in the global minimum the basic atoms (N, O and F) behave as hyper-coordinated centers. Also importantly, the Be4 cluster and its complexes present RHF-UHF instabilities (not reported before for Be4), which require the use of multireference methods to correctly describe them.

Significant bonding rearrangements triggered by Mg4 clusters.

The structure, stability, and bonding of the complexes formed by the interaction of Mg4 clusters and first row Lewis bases, namely, ammonia, water, and hydrogen fluoride, have been investigated through the use of high-level G4 single-reference and CASPT2 multireference formalisms. The adducts formed reflect the high electrophilicity of the Mg4 cluster through electron density holes in the neighborhood of each metallic center. After the adduct formation, the metallic bonding of the Mg4 moiety is not significantly altered so that the hydrogen shifts from the Lewis base toward the Mg atoms lead to new local minima with enhanced stability. For the particular case of ammonia and water, the global minima obtained when all the hydrogens of the Lewis base are shifted to the Mg4 moiety have in common a very stable scaffold with a N or an O center covalently tetracoordinated to the four Mg atoms, so the initial bonding arrangements of both reactants have completely disappeared. The reactivity features exhibited by these Mg4 clusters suggest that nanostructures of this metal might have an interesting catalytic behavior.

The Importance of Strain (Preorganization) in Beryllium Bonds.

In order to explore the angular strain role on the ability of Be to form strong beryllium bonds, a theoretical study of the complexes of four beryllium derivatives of orthocloso-carboranes with eight molecules (CO, N2, NCH, CNH, OH2, SH2, NH3, and PH3) acting as Lewis bases has been carried out at the G4 computational level. The results for these complexes, which contain besides Be other electron-deficient elements, such as B, have been compared with the analogous ones formed by three beryllium salts (BeCl2, CO3Be and SO4Be) with the same set of Lewis bases. The results show the presence of large and positive values of the electrostatic potential associated to the beryllium atoms in the isolated four beryllium derivatives of ortho-carboranes, evidencing an intrinsically strong acidic nature. In addition, the LUMO orbital in these systems is also associated to the beryllium atom. These features led to short intermolecular distances and large dissociation energies in the complexes of the beryllium derivatives of ortho-carboranes with the Lewis bases. Notably, as a consequence of the special framework provided by the ortho-carboranes, some of these dissociation energies are larger than the corresponding beryllium bonds in the already strongly bound SO4Be complexes, in particular for N2 and CO bases. The localized molecular orbital energy decomposition analysis (LMOEDA) shows that among the attractive terms associated with the dissociation energy, the electrostatic term is the most important one, except for the complexes with the two previously mentioned weakest bases (N2 and CO), where the polarization term dominates. Hence, these results contribute to further confirm the importance of bending on the beryllium environment leading to strong interactions through the formation of beryllium bonds.

Mutual Influence of Pnicogen Bonds and Beryllium Bonds: Energies and Structures in the Spotlight.

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 (NH3, NHCH2, and NCH), a phosphine (fluorophosphane, PH2F) and a beryllium derivative (BeH2, BeF2, BeCl2, BeCO3, and BeSO4). 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, CH2HN:PH2F:BeSO4. In the complexes between the strongest nitrogen bases and the strongest beryllium donors, phosphorus-shared and phosphorus-transfer bonds are found.

Are Anions of Cyclobutane Beryllium Derivatives Stabilized through Four-Center One-Electron Bonds?

High-level G4 ab initio calculations allowed us to show that C4H4(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.

Relativistic Effects on NMR Parameters of Halogen-Bonded Complexes.

Relativistic effects are found to be important for the estimation of NMR parameters in halogen-bonded complexes, mainly when they involve the heavier elements, iodine and astatine. A detailed study of 60 binary complexes formed between dihalogen molecules (XY with X, Y = F, Cl, Br, I and At) and four Lewis bases (NH3, H2O, PH3 and SH2) was carried out at the MP2/aug-cc-pVTZ/aug-cc-pVTZ-PP computational level to show the extent of these effects. The NMR parameters (shielding and nuclear quadrupolar coupling constants) were computed using the relativistic Hamiltonian ZORA and compared to the values obtained with a non-relativistic Hamiltonian. The results show a mixture of the importance of the relativistic corrections as both the size of the halogen atom and the proximity of this atom to the basic site of the Lewis base increase.

On the Lewis Basicity of Phosphoramides: A Critical Examination of Their Donor Number through Comparison of Enthalpies of Adduct Formation with SbCl5 and BF3.

The Lewis basicity of a series of phosphoryl compounds was examined using DFT and ab initio methods, including solvation effects. The enthalpies of adduct formation with two archetypal Lewis acids, antimony pentachloride and boron trifluoride, used to define the donor number DN and the BF3 affinity (BF3 A) respectively, were examined. The BF3 adducts allow the use of the high-accuracy G4 approach, whereas for SbCl5 adducts, three different DFT formalisms, including empirical dispersion corrections, were used because the G4 formalism is not available for third-row elements. For a comparison with experimental data, solvation effects were taken into account by using the polarizable continuum model. The experimental BF3 affinities were well reproduced by G4 calculations when including PCM solvation. Conversely, comparisons of our calculated values and experimental results reported in the literature show that SbCl5 enthalpies for phosphoramides are in error. In particular the DN for HMPA should be revised.

Combined Experimental and Theoretical Survey of the Gas-Phase Reactions of Serine-Ca2+ Adducts.

The association of Ca2+ 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-Ca2+ 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+ + [C3,H6,N,O2]+ or [C2,H4,O,N]+ + [Ca(C,H3,O2)]+, and to the loss of neutrals, namely, CH2O and H2O. 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-Ca2+ reactions, may also play a role.

Weak Interactions Get Strong: Synergy between Tetrel and Alkaline-Earth Bonds.

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··· SiH3F···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, 29Si 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 29Si 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.

Ternary Complexes Stabilized by Chalcogen and Alkaline-Earth Bonds: Crucial Role of Cooperativity and Secondary Noncovalent Interactions.

High-level G4 calculations show that the strength of chalcogen interactions is enhanced dramatically if chalcogen compounds simultaneously form alkaline-earth bonds. This phenomenon is studied by exploring binary YX2 ⋅⋅⋅N-Base complexes and two types of ternary MCl2 ⋅⋅⋅YX2 ⋅⋅⋅N-Base, YX2 ⋅⋅⋅N-Base⋅⋅⋅MCl2 complexes, in which YX2 is a chalcogen compound (Y=S, Se; X=F, Cl), the N-Bases are sp, sp2 , and sp3 bases (NCH, HN=CH2 , NH3 ), and MCl2 are alkaline-earth BeCl2 or MgCl2 derivatives. Starting from the chalcogen-bonded complexes YX2 ⋅⋅⋅NH3 and YX2 ⋅⋅⋅HN=CH2 , the binding site of a new incoming alkaline-earth bond is found, surprisingly, to depend on the nature of the halogen atom attached to the chalcogen. For the YF2 binary complexes the association site is the F atom of the YF2 subunit, whereas for YCl2 it is the N atom of the nitrogen base. Regarding YX2 ⋅⋅⋅NCH complexes, N is the most favorable site for an alkaline-earth interaction in ternary complexes, regardless of which YX2 derivative is used. The explanation relies on the interplay of all the noncovalent interactions involved: the strong cooperativity between chalcogen and alkaline-earth bonds, and the appearance of secondary noncovalent interactions in the form of hydrogen bonds.

Modulating the intrinsic reactivity of molecules through non-covalent interactions.

Non-covalent interactions unavoidably involve a certain disturbance of the electronic density of the interacting systems. Such perturbations are particularly strong when dealing with electron deficient systems such as boron, beryllium, magnesium (pre-p elements) or calcium (a pre-d element) derivatives. Indeed, these compounds have been shown to modify the intrinsic reactivity of the systems interacting with them. In the first part of this paper, we present an overview on (i) how electron deficient systems, acting as Lewis acids, modulate the intrinsic acidity of Lewis bases, explaining for instance why a typical base, such as aniline, can be converted by association with borane into an acid as strong as phosphoric acid; (ii) how other weak non-covalent interactions, such as halogen bonds, permit one to modulate the intrinsic basicity of typical oxyacids changing them into strong Brønsted bases; (iii) how cooperativity between different non-covalent interactions may lead to the spontaneous formation of ion-pairs in the gas phase; (iv) how non-covalent interactions generate σ-holes in systems where this feature is not present; and (v) how these interactions can induce exergonic and spontaneous formation of neutral radicals. In the second part of the paper, we show, by using G4 high-level ab initio calculations, that the acidity enhancement phenomenon is a general mechanism whenever a given base interacts with non-protic and protic acids. In the non-protic acid case, the underlying mechanism behind the enhancement is similar to the one reported for electron-deficient compounds, whereas the protic acid case appears in complexes stabilized through conventional hydrogen bonds. We also show that the former could be classified as an a priori mechanism, whereas the latter would be an a posteriori mechanism. This same a posteriori mechanism is behind the significant basicity enhancement of water and ammonia when interacting with conventional N-bases. Finally, we present a detailed analysis of the role that deformation can play in the intensity and nature of these enhancements.