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.

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.

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 BX₃ and AlX₃ (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₂X₆ (X = F, Cl) clusters are in good agreement with previous studies in the literature, but in contrast with the general accepted idea that B₂F₆ and B₂Cl₆ 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₃F₉ and B₃Cl₉ C_3v 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(η⁵-C₅H₄)(η⁵-C₅H₅)) 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-(CH₂)₂OH (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 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.

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.

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, N₂, NCH, CNH, OH₂, SH₂, NH₃, and PH₃) 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 (BeCl₂, CO₃Be and SO₄Be) 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 SO₄Be complexes, in particular for N₂ 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 (N₂ 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 (NH₃, NHCH₂, and NCH), a phosphine (fluorophosphane, PH₂F) and a beryllium derivative (BeH₂, BeF₂, BeCl₂, BeCO₃, and BeSO₄). 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₂HN:PH₂F:BeSO₄. In the complexes between the strongest nitrogen bases and the strongest beryllium donors, phosphorus-shared and phosphorus-transfer bonds are found.

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.

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

High-level G4 ab initio calculations allowed us to show that C₄H₄(BeX)₄ (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.

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.

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 BF₃ affinity (BF₃ A) respectively, were examined. The BF₃ adducts allow the use of the high-accuracy G4 approach, whereas for SbCl₅ 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 BF₃ 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 SbCl₅ enthalpies for phosphoramides are in error. In particular the DN for HMPA should be revised.

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··· SiH₃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, ²⁹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 ²⁹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.

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

The association of Ca²⁺ to serine and the subsequent gas-phase unimolecular reactivity of the [Ca(Ser)]²⁺ (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]²⁺ 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²⁺ 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)]²⁺ adduct correspond to both Coulomb explosions, yielding either CaOH⁺ + [C₃,H₆,N,O₂]⁺ or [C₂,H₄,O,N]⁺ + [Ca(C,H₃,O₂)]⁺, and to the loss of neutrals, namely, CH₂O and H₂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²⁺ reactions, may also play a role.

Modulating the intrinsic reactivity of molecules through non-covalent interactions

Non-covalent interactions as tools for modifying molecular properties.

Alkaline-earth (Be, Mg and Ca) bonds at the origin of huge acidity enhancements

The interaction between alkaline-earth derivatives with the general formula X₂M (X = H, F and Cl; M = Be, Mg and Ca) and a set of Lewis bases, including first and second-row hydrides, namely YH_n (Y = O, N, F, S, P and Cl) hydrides, as well as other typical cyclic organic bases, such as aniline, 1H-1,2,3-triazole, 1H-tetrazole and phenylphosphine, was investigated using the G4 ab initio composite method. Contrary to what was expected, it was found that the interactions involving Mg and Ca derivatives were not necessarily weaker than those between beryllium bonds. The origin is two-fold: larger deformation of the interacting systems when Be-derivatives are involved and appearance of secondary non-covalent interactions in the formation of some of the Mg- and Ca-containing complexes. Hence, the dissociation of the latter complexes may require higher enthalpies than that of the Be complexes. These deformations are triggered by a significant redistribution of electron density of the two interacting moieties, which also result in dramatic changes in the reactivity of the interacting compounds and in particular in the intrinsic basicity of the Lewis bases investigated, to the point that conventional bases, such as ammonia or aniline, upon complexation with MCl₂ (M = Be, Mg and Ca), become stronger Brønsted acids than phosphoric acid, whereas other bases, such as 1H-tetrazole, become stronger acids than perchloric acid.

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.

Complexes between neutral oxyacid beryllium salts and dihydrogen: a possible way for hydrogen storage?

Very stable 1 : 1 and 1 : 2 complexes between oxyacid beryllium salts and H₂ are found.

Formation of unexpected silicon- and disiloxane-bridged multiferrocenyl derivatives bearing Si-O-CH[double bond, length as m-dash]CH2 and Si-(CH2)2C(CH3)3 substituents via cleavage of tetrahydrofuran and trapping of its ring fragments

The formation of a family of silicon- and siloxane-bridged multiferrocenyl derivatives carrying different functional groups attached to silicon, including Fc₂(CH₃)₃C(CH₂)₂SiCH[double bond, length as m-dash]CH₂ (5), Fc₂(CH₂[double bond, length as m-dash]CH-O)SiCH[double bond, length as m-dash]CH₂ (6), Fc₂(OH)SiCH[double bond, length as m-dash]CH₂ (7), Fc₂(CH₂[double bond, length as m-dash]CH-O)Si-O-Si(O-CH[double bond, length as m-dash]CH₂)Fc₂ (8) and Fc₂(CH₂[double bond, length as m-dash]CH-O)Si-O-SiFc₃ (9) is described. Silyl vinyl ether molecules 6, 8 and 9 and the heteroleptic vinylsilane 5 resulted from the competing metathesis reaction of lithioferrocene (FcLi), CH₂[double bond, length as m-dash]CH-OLi or (CH₃)₃C(CH₂)₂Li with the corresponding multifunctional chlorosilane, Cl₃SiCH[double bond, length as m-dash]CH₂ or Cl₃Si-O-SiCl₃. The last two organolithium species have been likely formed in situ by fragmentation of the tetrahydrofuran solvent. Diferrocenylvinyloxyvinylsilane 6 is noteworthy since it represents a rare example of a redox-active silyl mononomer in which two different C[double bond, length as m-dash]C polymerisable groups are directly connected to silicon. The molecular structures of the silicon-containing multiferrocenyl species 5, 6, 8 and 9 have been investigated by single-crystal X-ray diffraction studies, demonstrating the capture and storage processes of two ring fragments resulting from the cleavage of cyclic THF in redox-active and stable crystalline organometallic compounds. From electrochemical studies we found that by changing the anion of the supporting electrolyte from [PF₆]^- to [B(C₆F₅)₄]^-, the redox behaviour of tetrametallic disiloxane 8 can be switched from a poorly resolved multistep redox process to four consecutive well-separated one-electron oxidations, corresponding to the sequential oxidation of the four ferrocenyl moieties.

Beryllium-based fluorenes as efficient anion sponges

The F⁻, Cl⁻, CN⁻, NO₂⁻, NO₃⁻, and SO₄²⁻ anion affinities of 4,5-bis(BeX)-fluorene (X = H, F, Cl, CN, NC, and OCH₃) derivatives have been calculated at the B3LYP/6-311+G(3df,2p)//B3LYP/6-31+G(d,p) level of theory.

Beryllium subphthalocyanines self-assembling properties

Beryllium subphthalocyanines have been recently shown to be suitable candidates for photochemical devices if combined with appropriate donor systems. The ability of beryllium subphthalocyanines to self-assemble is explored for the first time by means of density functional theory calculations. Free dimers of beryllium subphtalocyanine and their corresponding complexes with water and pyridine are computed at the wB97X-D/6-311+G(d,p) level of theory. In contrast with the behavior reported for beryllium phthalocyanines, for beryllium subphtalocyanines, beryllium–aza–nitrogen intermolecular interactions are observed, suggesting that these species are likely to self-assemble. Aggregates of related structures such as beryllium subporphyrazines with axial groups confirm the importance of hydrogen bonds in the stacking.

On the existence of intramolecular one-electron Be–Be bonds

The electron attachment to 1,8-diBeX-naphthalene derivatives leads to rather stable radical anions through the formation of one-electron Be–Be bonds.

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

Ga(+) Basicity and Affinity Scales Based on High-Level Ab Initio Calculations

The structure, relative stability and bonding of complexes formed by the interaction between Ga(+) and a large set of compounds, including hydrocarbons, aromatic systems, and oxygen-, nitrogen-, fluorine and sulfur-containing Lewis bases have been investigated through the use of the high-level composite ab initio Gaussian-4 theory. This allowed us to establish rather accurate Ga(+) cation affinity (GaCA) and Ga(+) cation basicity (GaCB) scales. The bonding analysis of the complexes under scrutiny shows that, even though one of the main ingredients of the Ga(+) -base interaction is electrostatic, it exhibits a non-negligible covalent character triggered by the presence of the low-lying empty 4p orbital of Ga(+) , which favors a charge donation from occupied orbitals of the base to the metal ion. This partial covalent character, also observed in AlCA scales, is behind the dissimilarities observed when GaCA are compared with Li(+) cation affinities, where these covalent contributions are practically nonexistent. Quite unexpectedly, there are some dissimilarities between several Ga(+) -complexes and the corresponding Al(+) -analogues, mainly affecting the relative stability of π-complexes involving aromatic compounds.