Exploration of Free Energy Surface of the Au10 Nanocluster at Finite Temperature.

The first step in comprehending the properties of Au10 clusters is understanding the lowest energy structure at low and high temperatures. Functional materials operate at finite temperatures; however, energy computations employing density functional theory (DFT) methodology are typically carried out at zero temperature, leaving many properties unexplored. This study explored the potential and free energy surface of the neutral Au10 nanocluster at a finite temperature, employing a genetic algorithm coupled with DFT and nanothermodynamics. Furthermore, we computed the thermal population and infrared Boltzmann spectrum at a finite temperature and compared it with the validated experimental data. Moreover, we performed the chemical bonding analysis using the quantum theory of atoms in molecules (QTAIM) approach and the adaptive natural density partitioning method (AdNDP) to shed light on the bonding of Au atoms in the low-energy structures. In the calculations, we take into consideration the relativistic effects through the zero-order regular approximation (ZORA), the dispersion through Grimme's dispersion with Becke-Johnson damping (D3BJ), and we employed nanothermodynamics to consider temperature contributions. Small Au clusters prefer the planar shape, and the transition from 2D to 3D could take place at atomic clusters consisting of ten atoms, which could be affected by temperature, relativistic effects, and dispersion. We analyzed the energetic ordering of structures calculated using DFT with ZORA and single-point energy calculation employing the DLPNO-CCSD(T) methodology. Our findings indicate that the planar lowest energy structure computed with DFT is not the lowest energy structure computed at the DLPN0-CCSD(T) level of theory. The computed thermal population indicates that the 2D elongated hexagon configuration strongly dominates at a temperature range of 50-800 K. Based on the thermal population, at a temperature of 100 K, the computed IR Boltzmann spectrum agrees with the experimental IR spectrum. The chemical bonding analysis on the lowest energy structure indicates that the cluster bond is due only to the electrons of the 6 s orbital, and the Au d orbitals do not participate in the bonding of this system.

Structure effects of Pt15 clusters for the oxygen reduction reaction: first-principles calculations.

In the present work, the lowest energy structures and electronic properties of Pt15 clusters are investigated using molecular dynamics simulations. The results showed that the most stable configuration is a capped pyramidal structure, which is 0.8 kal mol-1 lower in energy than a layered structure previously reported [V. Kumar and Y. Kawazoe, Evolution of Atomic and Electronic Structure of Pt Clusters: Planar, Layered, Pyramidal, Cage, Cubic, and Octahedral Growth, Phys. Rev. B: Condens. Matter Mater. Phys., 2008, 77, 205418.]. The result is further confirmed by using both the PW91/cc-pVDZ-PP and PBE/PW approaches including the other representative isomers for Pt15. Due to the interesting structure arrangements found, we have investigated the catalytic activities for the oxygen reduction reaction. We found that the most stable Pt15 clusters are plausible catalyts for the ORR according to their interaction with oxygen species, which is consistent with experiments of Pt clusters with atomicity below 20. The results of the structure, electronic, adsorption and vibrational properties of the clusters are provided.

Correction: Structure effects of Pt15 clusters for the oxygen reduction reaction: first-principles calculations.

Correction for 'Structure effects of Pt15 clusters for the oxygen reduction reaction: first-principles calculations' by Peter L. Rodríguez-Kessler et al., Phys. Chem. Chem. Phys., 2023, https://doi.org/10.1039/d2cp05188e.

Hitting the Bull's Eye: Stable HeBeOH+ Complex.

It is now known that the heavier noble gases (Ng=Ar-Rn) show some varying degrees of reactivity with a gradual increase in reactivity along Ar-Rn. However, because of their very small size and very high ionization potential, helium and neon are the hardest targets to crack. Although few neon complexes are isolated at very low temperatures, helium needs very extreme situations like very high pressure. Here, we find that protonated BeO, BeOH+ can bind helium and neon spontaneously at room temperature. Therefore, extreme conditions like very low temperature and/or high pressure will not be required for their experimental isolation. The Ng-Be bond strength is very high for their heavier homologs and the bond strength shows a gradual increase from He to Rn. Moreover, the Ng-Be attractive energy is almost exclusively originated from the orbital interaction which is composed of one Ng(s/pσ )→BeOH+ σ-donation and two weaker Ng(pπ )→BeOH+ π-donations, except for helium. Helium uses its low-lying vacant 2p orbitals to accept π-electron density from BeOH+ . Previously, such electron-accepting ability of helium was used to explain a somewhat stronger helium bond than neon for neutral complexes. However, the present results indicate that such π-back donations are too weak in nature to decide any energetic trend between helium and neon.

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 pentacoordinate carbon in CGa5+ derivatives.

We report a family of systems having a planar pentacoordinate carbon (ppC) based on the next heavier analogue of CAl5+, the ppC system par excellence. Although because of the larger size of Ga, the ppC isomer is not even a local minimum in CGa5+, a single isoelectronic substitution of Ga by smaller sized Be maximizes the bonding in the ppC form. Retaining the 18 valence electron rule, the global minimum structures of CGa4Be, CGa3Be2-, CGa2Be32-, and CGaBe43- clusters and their corresponding lithium salts have a ppC.

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.

E3 M3+ (E=C-Pb, M=Li-Cs) Clusters: The Smallest Molecular Stars.

Extensive potential energy surface explorations of twenty-five clusters with the formula E3 M3+ (E=Group 14 element and M=Group 1 element) through density functional theory and high-level ab initio computations reveal that the lowest-energy isomer for all these systems corresponds to a non-classical D3h star-like structure in the singlet state, where three M atoms interact electrostatically with the triangular E3 core, occupying three bridging positions around it. More than 18 200 calculations were done in the search for the minima structures, starting with a first phase at the PBE0/LANL2DZ level and ending with an analysis of the most representative clusters at the CCSD(T)/def2-TZVP//PBE0/def2-TZVP level. The title clusters represent the smallest molecular stars with three planar tetracoordinate E atoms (E=Group 14 element). All these E3 M3+ clusters behave like superalkali cations with small vertical electron affinities (smaller than Cs), large vertical electron detachment energies, and HOMO-LUMO energy gaps. Their energetics, bonding, and electron delocalization are discussed in detail. The high stability of these clusters is reflected from the large dissociation energy needed for different dissociation channels. The electron delocalization is confirmed by the presence of two delocalized π electrons over the E3 core and strong diatropic responses.

Does H4SO5 exist?

The possible existence of H4SO5 in aqueous sulfuric acid is analyzed in detail. For bare H4SO5, the computed free energy barrier for the exergonic transformation of H4SO5 into the H2SO4H2O complex is only 3.8 kcal mol-1. The presence of water or sulfuric acid catalyzes the dehydration to such an extent that it becomes almost a barrierless process. In the gas phase, dehydration of H4SO5 is an autocatalytic reaction as the water molecule produced by the decomposition of one H4SO5 molecule induces further dissociation. Thus, in solution, the surrounding water molecules make the para-sulfuric acid a very vulnerable species to exist. The simulated Raman spectra also corroborate the absence of H4SO5 in solution.

Coaxial Triple-Layered versus Helical Be6 B11- Clusters: Dual Structural Fluxionality and Multifold Aromaticity.

Two low-lying structures are unveiled for the Be6 B11- nanocluster system that are virtually isoenergetic. The first, triple-layered cluster has a peripheral B11 ring as central layer, being sandwiched by two Be3 rings in a coaxial fashion, albeit with no discernible interlayer Be-Be bonding. The B11 ring revolves like a flexible chain even at room temperature, gliding freely around the Be6 prism. At elevated temperatures (1000 K), the Be6 core itself also rotates; that is, two Be3 rings undergo relative rotation or twisting with respect to each other. Bonding analyses suggest four-fold (π and σ) aromaticity, offering a dilute and fluxional electron cloud that lubricates the dynamics. The second, helix-type cluster contains a B11 helical skeleton encompassing a distorted Be6 prism. It is chiral and is the first nanosystem with a boron helix. Molecular dynamics also shows that at high temperature the helix cluster readily converts into the triple-layered one.

How Many Water Molecules Does it Take to Dissociate HCl?

The potential energy surfaces of the HCl(H2O)n (n is the number of water molecules) clusters are systematically explored using density functional theory and high-level ab initio computations. On the basis of electronic energies, the number of water molecules needed for HCl dissociation is four as reported by some experimental groups. However, this number is five owing to the inclusion of entropic factors. Wiberg bond indices are calculated and analyzed, and the results provide a quadratic correlation and classification of clusters according to the nondissociated, partially dissociated, and fully dissociated character of the H-Cl bond. Our computations show that if temperature is not controlled during the experiment, the values obtained for the dipole moment (or for any measurable property) are susceptible to change, providing a different picture of the number of water molecules needed for HCl dissociation in a nanoscopic droplet.

10-π-Electron arenes à la carte: structure and bonding of the [E-(CnHn)-E](n-6) (E = Ca, Sr, Ba; n = 6-8) complexes.

In this paper, we provide solid evidence to show that among an overwhelming structural diversity, alkaline earth metals (Ca, Sr, Ba) have the ability to form inverted sandwich compounds with C6H6, C7H7(+), and C8H8(2+) of Dnh symmetry and general formula [E-(CnHn)-E](n-6) (n = 6-8) with planar 10-π-electron aromatic cores by virtue of transferring two electrons per metal atom to the ring. However, the origin of the orbital interaction between the metals and the carbon ring is quite different; while [E-(C6H6)-E] complexes are dominated by δ-interactions, both π- and δ-interactions are important in [E-(C7H7)-E](+) and [E-(C8H8)-E](2+) complexes.

How strong are the metallocene-metallocene interactions? Cases of ferrocene, ruthenocene, and osmocene.

An exhaustive exploration of the potential energy surfaces of ferrocene, ruthenocene and osmocene dimers has been performed. Our computations involving dispersion show that only four different isomers are present in each metallocene dimer. The collective action of small interaction energies of dispersive nature leads to a dissociation energy of 7.5 kcal mol(-1) for the ferrocene dimer. Dispersion has strong effects on the geometrical parameters, reducing the M···M distances by almost 1 Å. Our results also reveal that inclusion of entropic factors modifies the relative stability of the complexes. The nature of bonding is examined using the energy decomposition analysis and the non-covalent interaction index. Both analyses indicate that dispersion is the major contributing factor in stabilizing a metallocene dimer.

Planar pentacoordinate carbons in CBe5(4-) derivatives.

The potential energy surfaces of a series of clusters with formula CBe5Lin(n-4) (n = 1 to 5) have been systematically explored. Our computations show that the lithium cations preserve the CBe5(4-) pentagon, such that the global minimum structure for these series of clusters has a planar pentacoordinate carbon (ppC) atom. The systems are primarily connected via a network of multicenter σ-bonds, in which the C atom acts as σ-acceptor and this acceptance of charge is balanced by the donation of the 2pz electrons to the π-cloud. The induced magnetic field analysis suggests that the clusters with formula CBe5Lin(n-4) (n = 1 to 5) are fully delocalized. The fact that these ppC-containing clusters are the lowest-energy forms on the corresponding potential energy surfaces raises expectations that these species can be prepared experimentally in the gas phase.

Planar tetracoordinate carbons with a double bond in CAl3E clusters.

The potential energy surfaces of a series of clusters with the formula CAl3E (E = P, As, Sb, Bi) are systematically explored using density functional theory and high level ab initio calculations. The global minimum structure of these clusters contains a planar tetracoordinate carbon atom. The presence of a C=E double bond is supported by the Wiberg bond indices, the adaptive natural density partitioning analysis, and the magnetic response. Our results show that these planar tetracoordinate carbon clusters are both thermodynamically and kinetically viable species in the gas phase.

Intramolecular structure and dynamics of mequinol and guaiacol in the gas phase: Rotationally resolved electronic spectra of their S1 states.

The molecular structures of guaiacol (2-methoxyphenol) and mequinol (4-methoxyphenol) have been studied using high resolution electronic spectroscopy in a molecular beam and contrasted with ab initio computations. Mequinol exhibits two low frequency bands that have been assigned to electronic origins of two possible conformers of the molecule, trans and cis. Guaiacol also shows low frequency bands, but in this case, the bands have been assigned to the electronic origin and vibrational modes of a single conformer of the isolated molecule. A detailed study of these bands indicates that guaiacol has a vibrationally averaged planar structure in the ground state, but it is distorted along both in-plane and out-of-plane coordinates in the first electronically excited state. An intramolecular hydrogen bond involving the adjacent -OH and -OCH3 groups plays a major role in these dynamics.

Photoinduced C-C reactions on insulators toward photolithography of graphene nanoarchitectures.

On-surface chemistry for atomically precise sp(2) macromolecules requires top-down lithographic methods on insulating surfaces in order to pattern the long-range complex architectures needed by the semiconductor industry. Here, we fabricate sp(2)-carbon nanometer-thin films on insulators and under ultrahigh vacuum (UHV) conditions from photocoupled brominated precursors. We reveal that covalent coupling is initiated by C-Br bond cleavage through photon energies exceeding 4.4 eV, as monitored by laser desorption ionization (LDI) mass spectrometry (MS) and X-ray photoelectron spectroscopy (XPS). Density functional theory (DFT) gives insight into the mechanisms of C-Br scission and C-C coupling processes. Further, unreacted material can be sublimed and the coupled sp(2)-carbon precursors can be graphitized by e-beam treatment at 500 °C, demonstrating promising applications in photolithography of graphene nanoarchitectures. Our results present UV-induced reactions on insulators for the formation of all sp(2)-carbon architectures, thereby converging top-down lithography and bottom-up on-surface chemistry into technology.

Structural evolution of small gold clusters doped by one and two boron atoms.

The potential energy surfaces (PES) of a series of gold-boron clusters with formula Aun B (n = 1-8) and Aum B2 (m = 1-7) have been explored using a modified stochastic search algorithm. Despite the complexity of the PES of these clusters, there are well-defined growth patterns. The bonding of these clusters is analyzed using the adaptive natural density partitioning and the natural bonding orbital analyses. Reactivity is studied in terms of the molecular electrostatic potential.

Exploring the potential energy surface of E2P4 clusters (E=Group 13 element): the quest for inverse carbon-free sandwiches.

Inverse carbon-free sandwich structures with formula E2P4 (E=Al, Ga, In, Tl) have been proposed as a promising new target in main-group chemistry. Our computational exploration of their corresponding potential-energy surfaces at the S12h/TZ2P level shows that indeed stable carbon-free inverse-sandwiches can be obtained if one chooses an appropriate Group 13 element for E. The boron analogue B2P4 does not form the D(4h)-symmetric inverse-sandwich structure, but instead prefers a D(2d) structure of two perpendicular BP2 units with the formation of a double B-B bond. For the other elements of Group 13, Al-Tl, the most favorable isomer is the D(4h) inverse-sandwich structure. The preference for the D(2d) isomer for B2P4 and D(4h) for their heavier analogues has been rationalized in terms of an isomerization-energy decomposition analysis, and further corroborated by determination of aromaticity of these species.