Planar tetracoordinate beryllium compounds with a partially covalent Be-Ng bond.

An analysis of the thermodynamic and kinetic stability and the nature of the chemical bond in hypercoordinated compounds with the formula BeH3Ng+ (Ng = He-Rn) through high-level calculations is presented in this work. Thermochemical calculations show that, for the heavier noble gases (Ar-Rn), these systems are thermodynamically stable at room temperature; however, this stability decreases due to a weakening of the Be-H2 interaction, while the Be-Ng bond strengthens going down the periodic table. These results are complemented by Born Oppenheimer molecular dynamics simulations, in which the increasing tendency to dissociate the Be-H2 bond is evidenced. The nature of the chemical bonding depends on the analysis performed. On the one hand, the interacting quantum atoms method indicates that the covalent contribution is around 25 to 30%. On the other hand, the electron density topology indicates a covalent nature for compounds with Kr-Rn, while Hirshfeld population analysis in conjunction with Mayer's bond order establishes polar covalent behavior. The geometrical parameters and natural energy decomposition analysis (NEDA) indicate a covalent nature, allowing us to consider that the Be-Ng bond has a partially covalent character.

Breaking the plane: B5H5 is a three-dimensional structure.

In this study, we delved into the structure of B5H5 and questioned some of its accepted assumptions. By exploring the potential energy surface, we found a new three-dimensional structure as the global minimum. This finding is in contrast with the previously hypothesized planar and cage-like models. Our exploration extends to the kinetic stability of various B5H5 isomers, offering insights into the dynamic behavior of these molecules.

Li6 E5 Li6 : Tetrel Sandwich Complexes with 10-π-Electrons.

When (4n +2) π-electrons are located in single planar ring, it conventionally qualifies as aromatic. According Hückel's rule, systems possessing ten π-electrons should be aromatic. Herein we report a series of D5h  Li6 E5 Li6 sandwich structures, representing the first global minima featuring ten π-electrons E5 10- ring (E=Si-Pb). However, these π-electrons localize as five π-lone-pairs rather than delocalized orbitals. The high symmetry structure achieved is a direct consequence of σ-aromaticity, particularly favored in elements from Si to Pb, resulting in a pronounced diatropic ring current flow that contributes to the enhanced stability of these systems.

BH4 Ng+ (Ar-Rn): Viable Compounds with a B-Ng Covalent Bond.

In this work, we explore, using high-level calculations, the ability of BH4 + to interact with noble gases. The He system is energetically unstable, while the Ne system could only be observed at cryogenic temperatures. In the case of the Ar, Kr and Xe systems, all are energetically stable, even at room temperature. The different chemical bond descriptors reveal a covalent character between B and the noble gas from Ar to Rn. However, this interaction gradually weakens the multicentric bond between the boron atom and the H2 fragment. Thus, although BH4 Rn+ exhibits a strong covalent bond, it tends to dissociate at room temperature into BH2 Rn+ +H2 .

Strong carbon - noble gas covalent bond and fluxionality in hypercoordinate compounds.

Thermodynamic, kinetic, and chemical bonding analysis at the coupled cluster level has been carried out for a series of hypercoordinated carbon compounds with formula CH4Ng2+ (Ng = He-Rn). Results show that these compounds could be stable at room temperature and Born-Oppenheimer molecular dynamics simulations (BOMD) in conjunction with activation energies indicate high kinetic stability. In addition, all chemical bonding descriptors agree with a strong C-Ng covalent bond and a bonding pattern similar to that of CH5+. Finally, BOMD simulations showed that these compounds are fluxional, with a continuous formation/breaking of H-H and C-H bonds. To the best of the authors' knowledge, these results represent the first series of fluxional compounds possessing a covalent bond between a main group element and a noble gas atom.

Revisiting the potential-energy surface of CnBe3n+2H2n+22+ (n = 2-4) clusters: are planar pentacoordinate carbon structures the global minima?

Using various exploration strategies, in this study, we investigated the potential energy surfaces (PES) of CBe5H5+ and CnBe3n+2H2n+22+ (n = 2-4) clusters. Previous studies proposed that the planar pentacoordinate carbons (ppCs) were the global minima of these clusters. However, our study identified new putative global minima and competitive isomers, refuting some previous assignments. We employed several methods, including evolutive-inspired stochastic approaches guided by "chemical criteria", and ab initio molecular dynamics simulations at elevated temperatures. Our results showed that the size of the scanned population significantly affected the evolutive method and that constrained or guided procedures showed an advantage in identifying better minima for larger systems. This study demonstrated that using multiple complementary strategies can result in a wider variety of minima in a given energy range. Our findings provide valuable insights into exploring the potential energy surfaces of clusters, mainly medium-sized clusters, which could be the connections between small clusters and nanomaterials.

Bowl-shaped Cluster CuB12 - : A Viable Global Minimum with Twofold Aromaticity.

A low-lying structure is revealed for the CuB12 - cluster, which is bowl-shaped. It consists of a triangular CuB2 base and a B10 rim. Molecular dynamics simulations indicates its structural robustness; at an elevated temperature (600 K), the base rotates reversibly within the B10 perimeter. Chemical bonding analysis detects 2σ- and 3π-delocalized bonds, suggesting double aromaticity. This is also confirmed by two diatropic and concentric ring currents under an external magnetic field.

Planar Elongated B12 Structure in M3B12 Clusters (M = Cu-Au).

Here, it is shown that the M3B12 (M = Cu-Au) clusters' global minima consist of an elongated planar B12 fragment connected by an in-plane linear M3 fragment. This result is striking since this B12 planar structure is not favored in the bare cluster, nor when one or two metals are added. The minimum energy structures were revealed by screening the potential energy surface using genetic algorithms and density functional theory calculations. Chemical bonding analysis shows that the strong electrostatic interactions with the metal compensate for the high energy spent in the M3 and B12 fragment distortion. Furthermore, metals participate in the delocalized π-bonds, which infers an aromatic character to these species.

Planar Hypercoordinate Carbons in Alkali Metal Decorated CE3 2- and CE2 2- Dianions.

After exploring the potential energy surfaces of Mm CE2 p (E=S-Te, M=Li-Cs, m=2, 3 and p=m-2) and Mn CE3 q (E=S-Te, M=Li-Cs, n=1, 2, q=n-2) combinations, we introduce 38 new global minima containing a planar hypercoordinate carbon atom (24 with a planar tetracoordinate carbon and 14 with a planar pentacoordinate carbon). These exotic clusters result from the decoration of V-shaped CE2 2- and Y-shaped CE3 2- dianions, respectively, with alkali counterions. All these 38 systems fulfill the geometrical and electronic criteria to be considered as true planar hypercoordinate carbon systems. Chemical bonding analyses indicate that carbon is covalently bonded to chalcogens and ionically connected to alkali metals.

Li8 Si8 , Li10 Si9 , and Li12 Si10 : Assemblies of Lithium-Silicon Aromatic Units.

We report the global minima structures of Li8 Si8 , Li10 Si9 , and Li12 Si10 systems, in which silicon moieties maintain structural and chemical bonding characteristics similar to those of their building blocks: the aromatic clusters Td -Li4 Si4 and C2v -Li6 Si5 . Electron counting rules, chemical bonding analysis, and magnetic response properties verify the silicon unit's aromaticity persistence. This study demonstrates the feasibility of assembling silicon-based nanostructures from aromatics clusters as building blocks.

Theoretical Prediction of Structures, Vibrational Circular Dichroism, and Infrared Spectra of Chiral Be4B8 Cluster at Different Temperatures.

Lowest-energy structures, the distribution of isomers, and their molecular properties depend significantly on geometry and temperature. Total energy computations using DFT methodology are typically carried out at a temperature of zero K; thereby, entropic contributions to the total energy are neglected, even though functional materials work at finite temperatures. In the present study, the probability of the occurrence of one particular Be4B8 isomer at temperature T is estimated by employing Gibbs free energy computed within the framework of quantum statistical mechanics and nanothermodynamics. To identify a list of all possible low-energy chiral and achiral structures, an exhaustive and efficient exploration of the potential/free energy surfaces is carried out using a multi-level multistep global genetic algorithm search coupled with DFT. In addition, we discuss the energetic ordering of structures computed at the DFT level against single-point energy calculations at the CCSD(T) level of theory. The total VCD/IR spectra as a function of temperature are computed using each isomer's probability of occurrence in a Boltzmann-weighted superposition of each isomer's spectrum. Additionally, we present chemical bonding analysis using the adaptive natural density partitioning method in the chiral putative global minimum. The transition state structures and the enantiomer-enantiomer and enantiomer-achiral activation energies as a function of temperature evidence that a change from an endergonic to an exergonic type of reaction occurs at a temperature of 739 K.

Effects of Temperature on Enantiomerization Energy and Distribution of Isomers in the Chiral Cu13 Cluster.

In this study, we report the lowest energy structure of bare Cu13 nanoclusters as a pair of enantiomers at room temperature. Moreover, we compute the enantiomerization energy for the interconversion from minus to plus structures in the chiral putative global minimum for temperatures ranging from 20 to 1300 K. Additionally, employing nanothermodynamics, we compute the probabilities of occurrence for each particular isomer as a function of temperature. To achieve that, we explore the free energy surface of the Cu13 cluster, employing a genetic algorithm coupled with density functional theory. Moreover, we discuss the energetic ordering of isomers computed with various density functionals. Based on the computed thermal population, our results show that the chiral putative global minimum strongly dominates at room temperature.

Planar Hexacoordinate Carbons: Half Covalent, Half Ionic.

Herein, the first global minima containing a planar hexacoordinate carbon (phC) atom are reported. The fifteen structures belong to the CE3 M3 + (E=S-Te and M=Li-Cs) series and satisfy both geometric and electronic criteria to be considered as a true phC. The design strategy consisted of replacing oxygen in the D3h  CO3 Li3 + structure with heavy and less electronegative chalcogens, inducing a negative charge on the C atom and an attractive electrostatic interaction between C and the alkali-metal cations. The chemical bonding analyses indicate that carbon is covalently bonded to three chalcogens and ionically connected to the three alkali metals.

Embedding a Planar Hypercoordinate Carbon Atom into a [4n+2] π-System.

Through delicate tuning of the electronic structure, we report herein a rational design of seventeen new putative global minimum energy structures containing a planar tetra- or pentacoordinate carbon atom embedded in an aromatic hydrocarbon. These structures are the result of replacing three consecutive hydrogen atoms of an aromatic hydrocarbon by less electronegative groups, forming a multicenter σ-bond with the planar hypercoordinate carbon atom and participating in the π-electron delocalization. This strategy that maximizes both mechanical and electronic effects through aromatic architectures can be extended to several molecular combinations to achieve new and diverse compounds containing planar hypercoordinate carbon centers.

Relativistic effects on the aromaticity of E3M3H3 (E = C-Pb; M = N-Bi) benzene analogues.

The relativistic effects on the aromaticity of a set of benzene analogues, E3M3H3 (E = C-Pb; M = N-Bi) heterocycles, using magnetically induced current density (MICD) and the NICSzz component of the conventional nucleus independent chemical shift (NICS), is hereby examined. The relativistic effects were evaluated by means of four-component relativistic MICD, and two-component NMR relativistic shielding tensor methods. MICD and NICS were also computed in a non-relativistic fashion, to assess the influence of scalar-relativistic and spin-orbit effects. Most of the studied compounds exhibit a net diatropic ring current (aromatic), excluding the nitrogen-containing compounds which are non-aromatic (except for C3N3H3), in agreement with their higher E-N electronegativity difference. Only in the case of bismuth compounds, E3Bi3H3, aromaticity is substantially decreased when relativistic effects are included (mainly due to the spin-orbit contribution). The larger the mass of the system, the larger the magnitude of this change, in line with the expected relativistic effects for heavier elements. The analysis based on the NICSzz computations agrees with that of the MICD, thus supporting both the magnetic behavior and the aromatic character of these compounds.

Aromatic ouroboroi: heterocycles involving a σ-donor-acceptor bond and 4n + 2 π-electrons.

The aromaticity and dynamics of a set of recently proposed neutral 5- and 6-membered heterocycles that are closed by dative (donor-acceptor) or multi-center σ bonds, and have resonance forms with a Hückel number of π-electrons, are examined. The donors and acceptors in the rings include N, O, and F, and B, Be, and Mg, respectively. The planar geometry of the rings, coupled with evidence from different measures of aromaticity, namely the NICSzz, and NICSπzz components of the conventional nucleus independent chemical shifts (NICS), and ring current strengths (RCS), indicate non-trivial degrees of aromaticity in certain cases, including the cyclic C3B2OH6 and C3BOH5 isomers, both with three bonds to the O site in the ring. The former is lower in energy by at least 17.6 kcal mol-1 relative to linear alternatives obtained from molecular dynamics simulations in this work. Some of the other systems examined are best described as non-aromatic. Ring opening, closing, and isomerization are observed in molecular dynamics simulations for some of the systems studied. In a few cases, the ring indeed persists.

(Li6 Si5 )2-5 : The Smallest Cluster-Assembled Materials Based on Aromatic Si5 6- Rings.

Extensive explorations of their potential energy surfaces, combined with high-level quantum chemical computations, strikingly show that the lowest energy structures of the (Li6 Si5 )2-5 systems consist of 2-5 Si5 6- aromatic units, surrounded by Li+ counterions, respectively. These viable gas-phase compounds are the pioneering reported cases of oligomers made by planar aromatic silicon rings. Based on the key evidence that these oligomers are energetically favored, and that their silicon rings aromaticity is thoroughly preserved, the Li6 Si5 cluster is proposed as a potential assembly unit to build silicon-lithium nanostructures, thus opening new paths to design Zintl compounds at the nanoscale level.

Exploring the Potential Energy Surface of Trimetallic Deltahedral Zintl Ions: Lowest-Energy [Sn6Ge2Bi]3- and [(Sn6Ge2Bi)2]4- Structures.

The synthesis and structural characterization of the dimer [(Sn6Ge2Bi)2]4- raise the possibility of obtaining a broad variety of analogous compounds with different Sn/Ge/Bi proportions. Several combinations of nine atoms have been detected by electrospray mass spectrometry as potential assembly units. However, [(Sn6Ge2Bi)2]4- remains as the unique experimentally characterized species in this series. This fact has motivated us to explore its potential energy surface, as well as its monomers' [Sn6Ge2Bi]3-/2-, in an effort to gain insight into the factors that might be privileging the experimental viability of this species. Our results show that the lowest-energy [Sn6Ge2Bi]3- structure remains in its oxidized product [Sn6Ge2Bi]2-, which corresponds to that identified in the dimer [(Sn6Ge2Bi)2]4-. Additionally, local minima, very close in energy to the lowest-energy monomer, are chiral mixtures that dimerize into diverse structures with a probable energetic cost, making them noncompetitive isomers. Finally, the global minimum of the dimer [(Sn6Ge2Bi)2]4- presents the most stable monomers as assembly units. These results show the importance of considering the simultaneity of all of these conditions for the viability of these types of compounds.

Li2 B12 and Li3 B12 : Prediction of the Smallest Tubular and Cage-like Boron Structures.

An intriguing structural transition from the quasi-planar form of B12 cluster upon the interaction with lithium atoms is reported. High-level computations show that the lowest energy structures of LiB12 , Li2 B12 , and Li3 B12 have quasi-planar (Cs ), tubular (D6d ), and cage-like (Cs ) geometries, respectively. The energetic cost of distorting the B12 quasi-planar fragment is overcompensated by an enhanced electrostatic interaction between the Li cations and the tubular or cage-like B12 fragments, which is the main reason of such drastic structural changes, resulting in the smallest tubular (Li2 B12 ) and cage-like (Li3 B12 ) boron structures reported to date.

A Fukui function-guided genetic algorithm. Assessment on structural prediction of Sin (n = 12-20) clusters.

Theoretical studies are essential for the structural characterization of clusters, when it comes to rationalize their unique size-dependent properties and composition. However, the rapid growth of local minima on the potential energy surface (PES), with respect to cluster size, makes the candidate identification a challenging undertaking. In this article, we introduce a hybrid strategy to explore the PES of clusters. This proposal involves the use of a biased initial population of a genetic algorithm procedure. Each individual in this population is built by assembling small fragments, according to the best matching of the Fukui function. The performance of a genetic algorithm procedure. The performance of the method is assessed on the PES exploration of medium-sized Sin clusters (n = 12-20). The most relevant results are: (a) the method converges at almost half of the time used by the canonical version of the GA and, (b) in all the studied cases, with the exception of Si13 and Si16 , the method allowed to identify the global minimum (GM) and other important low-lying structures. Additionally, the apparent deficiency of the proposal to identify the GM was corrected when a Si atom, or other low-lying isomers, were considered to build the clusters. © 2017 Wiley Periodicals, Inc.