Modulating the Strength of Hydrogen Bonds through Beryllium Bonds

Strong Be-N Interaction Induced Complementary Chemical Tuning to Design a Dual-gated Single Molecule Junction

The interaction between pyridines and the π-hole of BeH2 leads to the formation of strong beryllium-bonded complexes. Theoretical investigations demonstrate that the Be-N bonding interaction can effectively regulate the electronic current through a molecular junction. The electronic conductance exhibits distinct switching behavior depending on the substituent groups at the para position of pyridine, highlighting the role of Be-N interaction as a potent chemical gate in the proposed device. The complexes exhibit short intermolecular distances ranging from 1.724 to 1.752 Å, emphasizing their strong binding. Detailed analysis of electronic rearrangements and geometric perturbations upon complex formation provides insights into the underlying reasons for the formation of such strong Be-N bonds, with bond strengths varying from -116.25 to -92.96 kJ/mol. Moreover, the influence of chemical substituents on the local electronic transmission of the beryllium-bonded complex offers valuable insights for the implementation of a secondary chemical gate in single-molecule devices. This study paves the way for the development of chemically gateable, functional single-molecule transistors, advancing the design and fabrication of multifunctional single-molecule devices in the nanoscale regime.

Dismantlement of ammonia upon interaction with Ben (n ≤ 10) clusters

The interaction of ammonia with Be_n (n < 1-10) clusters has been investigated by density functional theory and ab initio calculations. The main conclusion is that, regardless of the size of the Be cluster, neither the structure of ammonia nor that of the Be clusters are preserved due to a systematic dissociation of its NH bonds and a spontaneous H-shift toward the available Be atoms. This H migration not only leads to rather stable BeH bonds, but dramatically enhances the strength of the BeN bonds as well. Accordingly, the maximum stability is found for the interaction with the beryllium trimer, leading to a complex with three NBe and three BeH bonds. Another maximum in stability, although lower than that reached for n = 3, is found for the Be heptamer, since from n = 6, a new NBe bond is formed, so that complexes from n = 6 to n = 10 are characterized by the formation of a NBe₄ moiety, whose stability reaches a maximum at n = 7. The bonding characteristics of the different species formed are analyzed by means of AIM, NBO, ELF and AdNDP approaches.

Triel Bond Formed by Malondialdehyde and Its Influence on the Intramolecular H-Bond and Proton Transfer

Malondialdehyde (MDA) engages in a triel bond (TrB) with TrX₃ (Tr = B and Al; X = H, F, Cl, and Br) in three modes, in which the hydroxyl O, carbonyl O, and central carbon atoms of MDA act as the electron donors, respectively. A H···X secondary interaction coexists with the TrB in the former two types of complexes. The carbonyl O forms a stronger TrB than the hydroxyl O, and both of them are better electron donors than the central carbon atom. The TrB formed by the hydroxyl O enhances the intramolecular H-bond in MDA and thus promotes proton transfer in MDA-BX₃ (X = Cl and Br) and MDA-AlX₃ (X = halogen), while a weakening H-bond and the inhibition of proton transfer are caused by the TrB formed by the carbonyl O. The TrB formed by the central carbon atom imposes little influence on the H-bond. The BH₂ substitution on the central C-H bond can also realise the proton transfer in the triel-bonded complexes between the hydroxyl O and TrH₃ (Tr = B and Al).

Promotion of TH3 (T = Si and Ge) group transfer within a tetrel bond by a cation-π interaction

The possibility of the transfer of the TH₃ group across a tetrel bond is considered by ab initio calculations. The TB is constructed by pairing PhTH₃ (Ph = phenyl; T = Si and Ge) with bases NH₃, NHCH₂, and the C₃N₂H₄ carbene. The TH₃ moves toward the base but only by a small amount in these dimers. However, when a Be²⁺ or Mg²⁺ dication is placed above the phenyl ring, the tetrel bond strength is greatly magnified reaching up to nearly 100 kcal mol^-1. This dication also induces a much higher degree of transfer which can be best categorized as half-transfer for the two N-bases and a near complete transfer for the carbene.

Enhancement of tetrel bond involving tetrazole-TtR3 (Tt = C, Si; R = H, F). Promotion of SiR3 transfer by a triel bond

The combination of a CR3 (R = H, F) with a tetrazole can result in a moderate carbon bond, which can be further strengthened by a triel bond. On the other hand, SiR3 group is half transferred between the two N atoms in these conditions.

Enhancement of the Tetrel Bond by the Effects of Substituents, Cooperativity, and Electric Field: Transition from Noncovalent to Covalent Bond

The T⋅⋅⋅N tetrel bond (TB) formed between TX₃ OH (T=C, Si, Ge; X=H, F) and the Lewis base N≡CM (M=H, Li, Na) is studied by ab initio calculations at the MP2/aug-cc-pVTZ level. Complexes involving TH₃ OH contain a conventional TB with interaction energy less than 10 kcal/mol. This bond is substantially strengthened, approaching 35 kcal/mol and covalent character, when fluorosubstituted TF₃ OH is combined with NCLi or NCNa. Along with this enhanced binding comes a near equalization of the TB T⋅⋅⋅N and the internal T-O bond lengths, and the associated structure acquires a trigonal bipyramidal shape, despite a high internal deformation energy. This structural transformation becomes more complete, and the TB is further strengthened upon adding an electron acceptor BeCl₂ to the Lewis acid and a base to the NCM unit. This same TB strengthening can be accomplished also by imposition of an external electric field.

Classification of So-Called Non-Covalent Interactions Based on VSEPR Model

The variety of interactions have been analyzed in numerous studies. They are often compared with the hydrogen bond that is crucial in numerous chemical and biological processes. One can mention such interactions as the halogen bond, pnicogen bond, and others that may be classified as σ-hole bonds. However, not only σ-holes may act as Lewis acid centers. Numerous species are characterized by the occurrence of π-holes, which also may play a role of the electron acceptor. The situation is complicated since numerous interactions, such as the pnicogen bond or the chalcogen bond, for example, may be classified as a σ-hole bond or π-hole bond; it ultimately depends on the configuration at the Lewis acid centre. The disadvantage of classifications of interactions is also connected with their names, derived from the names of groups such as halogen and tetrel bonds or from single elements such as hydrogen and carbon bonds. The chaos is aggravated by the properties of elements. For example, a hydrogen atom can act as the Lewis acid or as the Lewis base site if it is positively or negatively charged, respectively. Hence names of the corresponding interactions occur in literature, namely hydrogen bonds and hydride bonds. There are other numerous disadvantages connected with classifications and names of interactions; these are discussed in this study. Several studies show that the majority of interactions are ruled by the same mechanisms related to the electron charge shifts, and that the occurrence of numerous interactions leads to specific changes in geometries of interacting species. These changes follow the rules of the valence-shell electron-pair repulsion model (VSEPR). That is why the simple classification of interactions based on VSEPR is proposed here. This classification is still open since numerous processes and interactions not discussed in this study may be included within it.

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.

Perturbing the O-H…O Hydrogen Bond in 1-oxo-3-hydroxy-2-propene

Ab initio MP2/aug'-cc-pVTZ calculations have been carried out to identify and characterize equilibrium structures and transition structures on the 1-oxo-3-hydroxy-2-propene: Lewis acid potential energy surfaces, with the acids LiH, LiF, BeH2, and BeF2. Two equilibrium structures, one with the acid interacting with the C=O group and the other with the interaction occurring at the O-H group, exist on all surfaces. These structures are separated by transition structures that present the barriers to the interconversion of the two equilibrium structures. The structures with the acid interacting at the C=O group have the greater binding energies. Since the barriers to convert the structures with interaction occurring at the O-H group are small, only the isomers with interaction occurring at the C=O group could be experimentally observed, even at low temperatures. Charge-transfer energies were computed for equilibrium structures, and EOM-CCSD spin-spin coupling constants 2hJ(O-O), 1hJ(H-O), and 1J(O-H) were computed for equilibrium and transition structures. These coupling constants exhibit a second-order dependence on the corresponding distances, with very high correlation coefficients.

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.

High-energetic and low-sensitive 1,3,5-triamino 2,4,6-trinitrobenzene (TATB) crystal: first principles investigation and Hirshfeld surface analysis

The N–H⋯O hydrogen bonding between neighboring TATB molecules within a planar layer cooperates with π⋯π interaction between two adjacent molecules from two adjacent layers in the TATB crystal.

Molecular Hydrogen as a Lewis Base in Hydrogen Bonds and Other Interactions

The second-order Møller–Plesset perturbation theory calculations with the aug-cc-pVTZ basis set were performed for complexes of molecular hydrogen. These complexes are connected by various types of interactions, the hydrogen bonds and halogen bonds are most often represented in the sample of species analysed; most interactions can be classified as σ-hole and π-hole bonds. Different theoretical approaches were applied to describe these interactions: Quantum Theory of ‘Atoms in Molecules’, Natural Bond Orbital method, or the decomposition of the energy of interaction. The energetic, geometrical, and topological parameters are analysed and spectroscopic properties are discussed. The stretching frequency of the H-H bond of molecular hydrogen involved in intermolecular interactions is considered as a parameter expressing the strength of interaction.

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.

Not Only Hydrogen Bonds: Other Noncovalent Interactions

In this review, we provide a consistent description of noncovalent interactions, covering most groups of the Periodic Table. Different types of bonds are discussed using their trivial names. Moreover, the new name “Spodium bonds” is proposed for group 12 since noncovalent interactions involving this group of elements as electron acceptors have not yet been named. Excluding hydrogen bonds, the following noncovalent interactions will be discussed: alkali, alkaline earth, regium, spodium, triel, tetrel, pnictogen, chalcogen, halogen, and aerogen, which almost covers the Periodic Table entirely. Other interactions, such as orthogonal interactions and π-π stacking, will also be considered. Research and applications of σ-hole and π-hole interactions involving the p-block element is growing exponentially. The important applications include supramolecular chemistry, crystal engineering, catalysis, enzymatic chemistry molecular machines, membrane ion transport, etc. Despite the fact that this review is not intended to be comprehensive, a number of representative works for each type of interaction is provided. The possibility of modeling the dissociation energies of the complexes using different models (HSAB, ECW, Alkorta-Legon) was analyzed. Finally, the extension of Cahn-Ingold-Prelog priority rules to noncovalent is proposed.

Competition and cooperativity of hydrogen-bonding and tetrel-bonding interactions involving triethylene diamine (DABCO), H2O and CO2in air

Triethylene diamine (DABCO) can interact with H₂O and CO₂in air to form dimeric and trimeric complexesviahydrogen bond, tetrel bond as well as van der Waals interactions.

Bonding between electron-deficient atoms: strong Lewis-acid character preserved in X–Y–X (X = B, Al; Y = Be, Mg) bridges

Beryllium bis(diazaborolyl) derivatives and their Mg and Al-containing analogues are stable compounds stabilized through covalent bonds between electron-deficient atoms, and behave as good Lewis acids.

Evaluation of the aggregation process in a mixture of propofol and benzocaine

We report on a mass-resolved IR spectrosopic study on propofol-benzocaine aggregates. This is a complex system due to the several conformational isomers that both monomers may adopt and to the combination of functional groups they present, which allow the molecules to interact in many possible ways. However, our results demonstrate that a single conformation is favored for each stoichiometry. In the heterodimer, propofol acts as a proton donor to the ester group of benzocaine, while the whole cluster is stabilized by dispersive forces. These dispersive forces account for an important part of the system's stabilization energy as the calculations suggest. Propofol does not show any affinity for the amino group of benzocaine, even when a second molecule of propofol is introduced. These results demonstrate the difficulty in anticipating the aggregation preferences of even small organic molecules.

Coordination of Be and Mg Centres by HCN Ligands - Be…N and Mg…N Interactions

ωB97XD/aug-cc-pVTZ calculations were performed for clusters of Z2+ cations (Z=Be and Mg) and HCN molecules (up to six molecules). The clusters of Be(CH3 )2 and Mg(CH3 )2 with HCN species were also calculated to analyse the influence of the Be/Mg-C formally covalent bonds on interactions of Be or Mg centre with ligands. The beryllium and magnesium centres possess different areas of a positive electrostatic potential that depend on a number of HCN ligands in the cluster considered. Numerous correlations between geometrical, energetic and topological parameters of the clusters considered are discussed since various theoretical approaches are applied; Quantum Theory of 'Atoms in Molecules', Natural Bond Orbital method and decomposition of the energy of interaction. The Be/Mg…N interactions classified as beryllium and magnesium bonds possess numerous characteristics which are known for the hydrogen bonds. Different types of coordination of Be and Mg centres analysed here exist also in crystal structures.

Effect of Magnesium Bond on the Competition Between Hydrogen and Halogen Bonds and the Induction of Proton and Halogen Transfer

HOX (X=Cl, Br, I, and At) can engage in either a H-bond (HB) or halogen bond (XB) with a base-like HCN, NH₃ , and imidazole. Although the former is energetically preferred for X=Cl and Br, it is the XB that is more stable for At, with I showing little preference. MgY₂ forms a Mg-bond with the O atom of HOX, which grows stronger in the order X=Cl<Br<I<At and Y=F<Cl<Br. When all three molecules are combined, both the Mg and the H/X bonds are cooperatively strengthened to a large degree. Rather than causing a reversal in the HB/XB competition, the Mg-bond acts primarily to amplify the natural preference within the dimer. The Mg-bond induces a certain degree of transfer from O to N of the bridging atom in the H/X bond. Comparison is also made with the effects of a Be-bond.

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.

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.

Systematic Coupled Cluster Study of Noncovalent Interactions Involving Halogens, Chalcogens, and Pnicogens

The noncovalent interactions of 32 complexes involving pnicogens, chalcogens, and halogens atoms were investigated at the CCSD(T)/aug-cc-pVTZ level of theory. Two different types of complexes could be distinguished on the basis of geometric parameters, electron difference densities, and the charge transfer mechanisms associated with each type. In the type I conformation, the monomers adopt a skewed orientation allowing charge to be transfer between both monomers, whereas in the type II conformation the monomers adopt a linear arrangement, maximizing charge transfer in only one direction. Type I complexes involving the interaction between pnicogens and chalcogens cannot be unambiguously defined as chalcogen or pnicogen bonds, they are an admixture of both. The charge transfer dependence on the conformation adopted by the complexes described in this work can serve as a novel conformationally driven design concept for materials.

Hydrogen bonds, and σ-hole and π-hole bonds – mechanisms protecting doublet and octet electron structures

For various interactions electron charge shifts try to protect the former doublet or octet electronic structure of the Lewis acid centre.

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.

Fully Hydrogenated Beryllium Nanoclusters

We present the ground state and energetically low structures of BenH2n nanoclusters as predicted using density functional theory (DFT) and employing the M06 meta-hybrid exchange-correlation functional. Results using the M06 functional are benchmarked against high accuracy coupled-cluster CCSD(T) and found to be in excellent agreement. For small values of n, the linear or polymeric form is the lowest energy geometry, while for sizes larger, n > 9 ring type and link type structures are the energetically lowest configurations. This trend has also been observed through ab initio molecular dynamics (AIMD) simulations at finite temperatures. In addition to the binding energies of the structures we report on polymerization energies, Be-H bond energies with respect to coordination details, hydrogen desorption energies of saturated and oversaturated species, as well as computed infrared spectra of all the ground state and energetically low lying structures presented. Furthermore, we find that the saturated polymeric forms of the nanoclusters cannot retain molecular hydrogen, in contrast to what is expected when zero point energy corrections are not taken into account.

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⁻.

Interplay between Beryllium Bonds and Anion-π Interactions in BeR2:C6X6:Y- Complexes (R = H, F and Cl, X = H and F, and Y = Cl and Br)

A theoretical study of the beryllium bonds in BeR₂:C₆X₆ (R = H, F, Cl and X = H and F) has been carried out by means of MP2/aug′-cc-pVDZ computational methods. In addition, the ternary complexes BeR₂:C₆X₆:Y⁻ (Y = Cl and Br) have been analyzed. Geometric, energetic and electronic aspects of the complexes have been taken into account. All the parameters analyzed provide a clear indication of favorable cooperativity in both interactions observed, beryllium bond and aromatic ring:anion interaction.

Dual functions of Lewis acid and base of Se in F2C=Se and their interplay in F 2CSe•••NH 3•••HX

High-level quantum chemical calculations of the ternary systems F2CSe∙∙∙NH3∙∙∙HX (X=BeH, BH2, OH, CN, OCH3, Cl, and F) and the corresponding binary systems have been carried out in view of geometries, vibrational frequencies, interaction energies, orbital interactions, and electron densities. The molecular electrostatic potentials of F2CSe demonstrate that the Se atom could play a dual role of Lewis acid and base to form a chalcogen bond with NH3 and a hydrogen bond or a covalent interaction with HX, respectively. The chalcogen bond can compete with the hydrogen bond for the complexes involving F2CSe, but the covalent interaction is far stronger than the chalcogen bond. In the ternary complexes, both types of interactions are strengthened by each other, characterized by a shorter binding distance, a larger electron density, and a stronger orbital interaction. The covalent interaction has a greater enhancing effect on the chalcogen bond than the hydrogen bond does, resulting in a prominent shortening of ~0.23 Å distance for the Se∙∙∙N distance in F2CSe∙∙∙NH3∙∙∙BH3. The enhancement of both interactions in the ternary complexes has been understood with the electrostatic potentials and orbital interactions. Graphical Abstract The dual functions of Lewis acid and base of Se in F2CSe are enhanced each other in the ternary complexes.