Cooperativity between hydrogen bonds and beryllium bonds in (H2O)(n)BeX2 (n = 1-3, X = H, F) complexes. A new perspective

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.

The dubious origin of beryllium toxicity

Four mechanisms have been proposed in the literature to explain beryllium toxicity; they can be divided in two groups of two mechanisms: (i) replacement type: models 1 and 2; (ii) addition type: models 3 and 4. At this moment is not possible to select the best model not even to establish if one of these models will be the ultimate mechanism of beryllium toxicity. However, it is important to know the still open discussion about something so important associated with one of the simplest elements of the periodic table.

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 BeF₂ 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 BeF₂ molecule and the (BeF₂)₂ 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 (BeF₂)₂ dimer, that is, a Lewis acid much stronger than the monomer. Although the azoles investigated have neighboring N-basic sites, their interaction with the (BeF₂)₂ 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-(BeF₂)₂ complexes investigated become stronger nitrogenous acids than oxyacids such as perchloric acid.

Malonaldehyde-like Systems: BeF2 Clusters—A Subtle Balance between Hydrogen Bonds, Beryllium Bonds, and Resonance

The stability of malonaldehyde is governed by intramolecular hydrogen bonds (IMHBs) as well as in malonaldehyde-like systems where oxygen is replaced by N or S at any of the basic sites. As beryllium bonds have been shown to strongly cooperate with hydrogen bonds, this work explores at the high level ab initio G4 level of theory the effect of including this non-covalent interaction in the system through its association with BeF2. Although malonaldehyde follows the expected trends, where the formation of a pseudocyclic form is favored also when IMHB and Be bonds are present, the subtle balance between both non-covalent interactions leads to some surprising results when other heteroatoms are involved, to the point that interaction energies can be much larger than expected or even cyclization is not favored. A complete analysis using different computational tools gives an answer to those cases escaping the predictable trends.

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 Be₄ 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 Be₂ dimer. However, whereas for Be₂ the first dissociation process involves a low activation barrier which is above the reactants, for Be₄ 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 Be₄ cluster. As a consequence, in the global minimum the basic atoms (N, O and F) behave as hyper-coordinated centers. Also importantly, the Be₄ cluster and its complexes present RHF-UHF instabilities (not reported before for Be₄), 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 Mg₄ 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 Mg₄ cluster through electron density holes in the neighborhood of each metallic center. After the adduct formation, the metallic bonding of the Mg₄ 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 Mg₄ 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 Mg₄ clusters suggest that nanostructures of this metal might have an interesting catalytic behavior.

Towards developing a criterion to characterize non-covalent bonds: a quantum mechanical study

Chemical bonds are central to chemistry, biology, and allied fields, but still, the criterion to characterize an interaction as a non-covalent bond has not been studied rigorously.

Some interesting features of the rich chemistry around electron-deficient systems

In this short review, different phenomena that are triggered by the interaction of different compounds or clusters of compounds with electron-deficient systems, in particular beryllium and boron compounds, have been discussed in some detail. Particular attention was devoted to the huge acidity enhancements that can be induced through the interaction of conventional bases with B or Be containing compounds, which change these conventional bases in extremely strong proton donors. We have paid also attention to the cooperativity between Be bonds with other weak interactions, which results in a substantial increase of their strength, that can lead in some specific cases to the spontaneous formation of ion-pairs in the gas phase. Finally, the behavior of different Be derivatives as electron and anion sponges is discussed as well as the conditions needed to have clusters exhibiting rather strong Be–Be bonds, even though the Be–Be interaction in Be2 dimer is extremely weak. Finally, some attention was paid to systems with extremely short Be–Be distances but without a bond.

Beryllium bonding: insights from the σ- and π-hole analysis

Beryllium bonding is actually a subclass of secondary bonding. Similar to the case of halogen bonding, the σ- and π-holes on the Be atom of the monomers give in zeroth approximation the direction of electrophilic attack favorable to the formation of beryllium bonds. The nature of beryllium bonding is purely electrostatic so that the symmetry-adapted perturbation theory energy decomposition perfectly explains the relevance of the polarization and dispersion contribution on the formation of the beryllium bond.

Modulating the intrinsic reactivity of molecules through non-covalent interactions

Non-covalent interactions as tools for modifying molecular properties.

Enhancement of Thermodynamic Gas-Phase Acidity and Basicity of Water by Means of Secondary Interactions

A series of A⋅water, B⋅water complexes (A=acid, B=base) are studied at the G4 level of theory to show that water acidity or basicity can be modulated by non-covalent interactions. Protic and non-protic acids interacting with water form hydrogen bonds or other kinds of non-covalent interactions, respectively, that may dramatically change the acidity of water up to almost 360 kJ ⋅ mol^-1 in terms of enthalpy. Similarly, hydrogen bonds responsible for the interaction between typical small nitrogen-containing Lewis bases and water can enhance the proton affinity of water by almost 300 kJ ⋅ mol^-1 . Our results reveal that these large enhancements are linearly related with the binding energy of the charged complexes, and are determined by the Lewis acid-base properties of the molecule involved in the interaction, allowing a quite precise modulation of the corresponding acid-base properties of water.

Fostering the Basic Instinct of Boron in Boron-Beryllium Interactions

A set of complexes L₂HB···BeX₂ (L = CNH, CO, CS, N₂, NH₃, NCCH₃, PH₃, PF₃, PMe₃, OH₂; X = H, F) containing a boron-beryllium bond is described at the M06-2X/6-311+G(3df,2pd)//M062-2X/6-31+G(d) level of theory. In this quite unusual bond, boron acts as a Lewis base and beryllium as a Lewis acid, reaching binding energies up to -283.3 kJ/mol ((H₂O)₂HB···BeF₂). The stabilization of these complexes is possible thanks to the σ-donor role of the L ligands in the L₂HB···BeX₂ structures and the powerful acceptor nature of beryllium. According to the topology of the density, these B-Be interactions present positive laplacian values and negative energy densities, covering different degrees of electron sharing. ELF calculations allowed measuring the population in the interboundary B-Be region, which varies between 0.20 and 2.05 electrons upon switching from the weakest ((CS)₂HB···BeH₂) to the strongest complex ((H₂O)₂HB···BeF₂). These B-Be interactions can be considered as beryllium bonds in most cases.

Cooperative and anticooperative effects in resonance assisted hydrogen bonds in merged structures of malondialdehyde

We analyzed non-additive effects in resonance assisted hydrogen bonds (RAHBs) in different β-enolones, which are archetypal compounds of these types of interactions. For this purpose, we used (i) potential energy curves to compute the formation energy, ΔE, of the RAHBs of interest in different circumstances along with (ii) tools offered by quantum chemical topology, namely, the Quantum Theory of Atoms In Molecules (QTAIM) and the Interacting Quantum Atoms (IQA) electronic energy partition. We established the effect that a given H-bond exerts over ΔE associated with another RAHB, determining in this way the cooperativity or the anticooperativity of these interactions. The mesomeric structures and the QTAIM delocalisation indices are consistent with the determined cooperative or anticooperative character of two given RAHBs. The HB cooperativity and anticooperativity studied herein are directly reflected in the IQA interaction energy E, but they are modulated by the surrounding hydrocarbon chain. The IQA decomposition of ΔE_coop, a measure of the cooperativity between a pair of interacting RAHBs, indicates that the analyzed H-bond cooperative/anticooperative effects are associated with greater/smaller (i) strengthening of the pseudo-bicyclic structure of the compounds of interest and (ii) electron localisations with its corresponding changes in the intra and intermolecular exchange-correlation contributions to ΔE. Overall, we expect that this investigation will provide valuable insights into the interplay among hydrogen bonded atoms and the π system in RAHBs contributing in this way to the understanding of the general features of H-bonds.

Large Proton-Affinity Enhancements Triggered by Noncovalent Interactions

High affinity: The proton affinity (PA) of the OH group of YH_x OH compounds is always increased by noncovalent interaction (NCI) with a Lewis base (LB; see figure). This PA enhancement depends on the type of NCI. Weak NCIs can give rise to PA enhancements equal to or even larger than strong NCIs. The binding energies of protonated species play a major role in the case of sigma-hole interactions.

Viable aromatic BenHn stars enclosing a planar hypercoordinate boron or late transition metal

Similar to B_n rings, star-like Be_nH_n rings can serve as the n-electron σ-donors for designing species with planar hypercoordinate atom.

The effect of anions on noncovalent interactions in model clusters of chalcogen-containing (CH3)2X (X = O, S, Se) molecules

A computational study of F⁻⋯(CH₃)₂O⋯CH₃F with F⁻ bound to the protons of the two methyl groups, found significant enhancement of the O⋯C interaction relative to the neutral (CH₃)₂O⋯CH₃F dyad.

Modeling cooperative effects in halogen-bonded infinite linear chains

Non-additivity in noncovalent interactions is an important aspect of complex systems that can lead to stronger (cooperative) interactions when three or more molecular units influence each other. The halogen bond (XB) is a highly-directional noncovalent interaction that has been found to be cooperative. Here the strength and nature of cooperativity arising in X-bonded infinite linear chains of cyanogen halides and 4-halopyridines are investigated by means of density functional theory calculations. It is found that cyanogen halide chains are highly cooperative (up to 77%), whereas pyridines are only slightly cooperative (below 21%). It is demonstrated that XB and its non-additivity can be modeled as the sum of a local term, which depends on first nearest-neighbors only, and long-range effective dipole-dipole attractions. It is shown that the local term in cyanogen halides primarily accounts for repulsive short-range screened Coulomb interactions, whereas in 4-halopyridines such a term also includes attractive contributions, which are particularly sizeable in some elongated XB conformations. This outcome reveals differences in the nature of the XBs formed in these molecular systems. Nevertheless, it is shown that both systems behave as effective point dipoles regarding cooperative effects, at any point of the XB dissociation path. As such, these results are useful contributions for the understanding and modeling of non-additive effects of noncovalent interactions.

Structure and stability of clusters of β-alanine in the gas phase: importance of the nature of intermolecular interactions

We present a theoretical study of neutral clusters of β-alanine molecules in the gas phase, (β-ala)_nn ≤ 5.

Modulating the strength of tetrel bonding through beryllium bonding

Quantum chemical calculations were performed to investigate the stability of the ternary complexes BeH2···XMH3···NH3 (X = F, Cl, and Br; M = C, Si, and Ge) and the corresponding binary complexes at the atomic level. Our results reveal that the stability of the XMH3···BeH2 complexes is mainly due to both a strong beryllium bond and a weak tetrel-hydride interaction, while the XMH3···NH3 complexes are stabilized by a tetrel bond. The beryllium bond with a halogen atom as the electron donor has many features in common with a beryllium bond with an O or N atom as the electron donor, although they do exhibit some different characteristics. The stability of the XMH3···NH3 complex is dominated by the electrostatic interaction, while the orbital interaction also makes an important contribution. Interestingly, as the identities of the X and M atoms are varied, the strength of the tetrel bond fluctuates in an irregular manner, which can explained by changes in electrostatic potentials and orbital interactions. In the ternary systems, both the beryllium bond and the tetrel bond are enhanced, which is mainly ascribed to increased electrostatic potentials on the corresponding atoms and charge transfer. In particular, when compared to the strengths of the tetrel and beryllium bonds in the binary systems, in the ternary systems the tetrel bond is enhanced to a greater degree than the beryllium bond. Graphical Abstract A tetrel bond can be strengthened greatly by a beryllium bond.

The nature of resonance-assisted hydrogen bonds: a quantum chemical topology perspective

Resonance Assisted Hydrogen Bonds (RAHBs) are particularly strong H-Bonds (HBs) which are relevant in several fields of chemistry. The traditional explanation for the occurrence of these HBs is built on mesomeric structures evocative of electron delocalisation in the system. Nonetheless, there are several theoretical studies which have found no evidence of such electron delocalisation. We considered the origin of RAHBs by employing Quantum Chemical Topology tools, more specifically, the Quantum Theory of Atoms in Molecules (QTAIM) and the Interacting Quantum Atoms energy partition. Our results indicate that the π-conjugated bonds allow for a larger adjustment of electron density throughout the H-bonded system as compared with non-conjugated carbonyl molecules. This rearrangement of charge distribution is a response to the electric field due to the H atom involved in the hydrogen bonding of the considered compounds. As opposed to the usual description of RAHB interactions, these HBs lead to a larger electron localisation in the system, and concomitantly to larger QTAIM charges which in turn lead to stronger electrostatic, polarization and charge transfer components of the interaction. Overall, the results presented here offer a new perspective on the cause of strengthening of these important interactions.

Hydrogen bonds in methane–water clusters

Water–methane clusters are stable at low temperatures as those found in Mars. Water cages enveloping methane are stable, although they present small probability to occur.

Fullerene and corannulene derivatives acting as insulators of Cl−and BeH2

Corannulene and its benzo-derivatives C_xH₁₀(x= 20–60) as prototypes of non-planar π-aromatic systems have been studied as insulators of BeH₂and Cl⁻.

On the existence and characteristics of π-beryllium bonds

The existence of π-beryllium bonds explains the stability of the complexes between ethylene and acetylene and BeX2 (X = H, F, and Cl) derivatives. These linkers involve a significant charge transfer from the π(CC) bonding orbitals into the empty p orbitals of Be and to a much smaller degree into the σ(BeH)* antibonding orbitals. The significant deformation of the BeX2 moiety and the slight deformation of the unsaturated hydrocarbon result in distortion energies as high as the dissociation energy of the complex. The π-beryllium bonds are about four times stronger than conventional π-hydrogen bonds and even stronger than the strongest π-hydrogen bond reported to date in the literature. The topology of their electron density is characterized as being very flat in the bonding region between the π-system and Be, which leads to topologically unstable structures close to catastrophe points. Among the functionals considered in our study M06 is the one that leads to values in better agreement with CCSD(T)/aug-cc-pVTZ calculations used as a reference. B3LYP underestimates some interactions, whereas M06-2X overestimates all of them. MP2 also yields good agreement with the CCSD(T) method.

Using beryllium bonds to change halogen bonds from traditional to chlorine-shared to ion-pair bonds

Dramatic synergistic cooperative effects between Be⋯F beryllium bonds and Cl⋯N halogen bonds in XYBe:FCl:N-base ternary complexes lead to changes in the halogen-bond type from traditional to chlorine-shared to ion-pair bonds.