Low-energy Isomers Research Articles

Hydrogen-Bonded Complexes of HPN• and HNP• Radicals with Carbon Monoxide.

Phosphorus mononitride (PN) is a carrier of phosphorus in the interstellar medium. As the simplest derivatives of PN, the radical species HPN• and HNP• have remained elusive. Herein, we report the generation, characterization, and photochemistry of HPN• and HNP• in N2-matrix at 3 K. Specifically, HPN• was formed as a weakly bonded complex with CO in the matrix by 254 nm photolysis of the novel phosphinyl radical HPNCO•. The •NPH-CO complex is extremely unstable, as it undergoes spontaneous isomerization to the lower-energy isomer •PNH-CO through fast quantum mechanical tunneling (QMT) with a half-life of 6.1 min at 3 K. Upon further irradiation at 254 nm, the reverse conversion of •PNH-CO to •NPH-CO along with dehydrogenation to yield PN was observed. The characterization •NPH-CO and •PNH-CO with matrix-isolation IR spectroscopy is supported by D, 15N, and 13C isotope labeling and quantum chemical calculations at the XYGJ-OS/AVTZ level of theory, and the mechanism by hydrogen atom tunneling is consistent with multidimensional instanton theory calculations.

Density functional theory-based study on the structural, electronic and spectral properties of gas-phase PbMg n - (n = 2-12) clusters.

Gas-phase PbMg n- (n = 2-12) cluster structures were globally searched on their potential energy surfaces by means of the CALYPSO prediction software. Structural optimization and calculations of properties such as relative energy and electronic structure were then carried out by density functional theory for each size of low energy isomer. The structural, relative stability, natural charge population, natural electronic configuration and distribution of the strongest peaks of the infrared and Raman spectra of the low energy isomers of PbMg n- (n = 2-12) clusters were systematically investigated in the present work. It was shown that the PbMg7- cluster ground state isomer exhibits the highest stability, for which special electronic excitation and chemical bonding analyses were performed. It is reasonable to believe that this work enriches the structural, spectroscopic and other data of magnesium-based clusters and provides some theoretical basis for possible future experimental syntheses.

Electronic Structures of Small Stoichiometric ZnxOx Clusters.

Stoichiometric ZnxOx clusters in the subnanometer size regime have been the topic of several computational and mass spectrometry studies that showed the particular stability of the stoichiometric species relative to nonstoichiometric species (ZnxOy, x ≠ y). In the current study, we present the angle-resolved anion photoelectron (PE) spectra of stoichiometric ZnxOx-clusters (2 ≤ x ≤ 5), which are interpreted with supporting computational studies that include natural orbital ionization calculations on detachment transition cross sections. All spectra show evidence of Dxh ring structures, which had been predicted to be the most stable structures in previous computational studies. However, a new lowest energy isomer is reported for the Zn2O2 anion and neutral, a bent chain, which is readily reconciled with the most intense feature in the Zn2O2- PE spectrum. The computed PE angular distributions (PADs) associated with the lowest energy cluster structures identified computationally agree with the experimental results, with the exception of Zn5O5-, the experimental PAD of which suggests that strong vibronic coupling may be introducing anomalies. While the lowest lying electronic state of the Zn2O2 chain structure is a triplet state, all neutral ring structures (including Zn2O2, the anion of which also populates the ion beam), favor a singlet electronic state. The computed singlet-triplet splitting of the Dxh structures increases monotonically with x. Overall, we find that the properties of the ring structures evolve smoothly, rather than in the punctuated manner typically seen in the small cluster size regime.

Vanadium oxide clusters in substellar atmospheres. A quantum chemical study

As a refractory material, vanadia (solid V$_2$O$_5$) is likely to be found as a condensate in the atmospheres of substellar objects such as exoplanets and brown dwarfs. However, the nature of the nanometer-sized vanadium oxide clusters that partake in the nucleation process is not well understood. We aim to understand the formation of cloud condensation nuclei in oxygen-rich substellar atmospheres by calculating the relevant fundamental properties of the energetically most favorable vanadium oxide molecules and clusters and, investigate how they contribute to the formation of condensation seeds. We applied a hierarchical optimization approach in order to find the most favourable structures for clusters of (VO) N and (VO$_2$) N for N=1-10, and of (V$_2$O$_5$) N for N=1-4, and to calculate their thermodynamical potentials. The candidate geometries are initially optimized by applying classical interatomic potentials; these are then refined at the B3LYP/cc-pVTZ level of theory to obtain accurate zero-point energies and thermochemical quantities. We present previously unreported vanadium oxide cluster structures as the lowest-energy isomers. Moreover, we report revised cluster energies and their thermochemical properties. Chemical equilibrium calculations are used to asses the impact of the updated and newly derived thermodynamic potentials on the gas-phase abundances of vanadium-bearing species. In chemical equilibrium, larger clusters from different stoichiometric families are found to be the most abundant vanadium-bearing species for temperatures below sim 1000 K, while molecular VO is the most abundant between sim 1000 K and sim 2000 K. We determine the nucleation rates of each stoichiometric family for a given (T$_ gas $, p$_ gas $) profile of a brown dwarf using both classical and non-classical nucleation theory. Small differences in the revised Gibbs free energies of the clusters have a large impact on the abundances of vanadium-bearing species in chemical equilibrium at temperatures below sim 1000 K. These abundance changes subsequently have an impact on the nucleation rates of each stoichiometric family. We find that with the revised and more accurate cluster data, non-classical nucleation rates are up to 15 orders of magnitude higher than classical nucleation rates.

Structural prediction of carbon cluster isomers with machine-learning potential

Structural prediction of low-energy isomers of carbon twelve-atom clusters is carried out using the recently developed machine-learning potential GAP-20. The GAP-20 agrees with density-functional theory calculations regarding geometric structures and average C-C bond lengths for most isomers. However, the GAP-20 substantially lowers the energies of cage-like structures, resulting in a wrong ground state. A comparison of the cohesive energies with the density-functional theory points out that the GAP-20 only gives good results for monocyclic rings. Two multicyclic rings appear as new low-energy isomers, which have yet to be discovered in previous research.

An Efficient Growth Pattern Algorithm (GrowPAL) for Cluster Structure Prediction.

Identifying the lowest energy isomers in large clusters is a major challenge. Here, we introduce the Growth Pattern Algorithm (GrowPAL), a new approach that generates initial seeds composed of n+1 atoms from the system with n atoms through an interstitial-type addition (I-type) mechanism. We evaluated the effectiveness of GrowPAL on Lennard-Jones (LJ) clusters with up to n = 80 atoms, verifying the algorithm's ability to find challenging minima such as LJ38 and the partially icosahedral LJ69 with fewer optimizations than existing methods. In addition, we discuss the advantages and limitations of GrowPAL using our deconstruction scheme, which identifies "forebears" structures to study growth pathways. Having evaluated the strengths and weaknesses of GrowPAL, we employed it to explore Sutton-Chen clusters containing 5 to 80 atoms, uncovering three new lowest energy forms. We then applied GrowPAL to boron clusters containing 8 to 15 atoms, successfully identifying all reported minima. Overall, GrowPAL offers a practical solution for efficiently identifying global minima in hierarchical systems, thereby reducing computational costs.

Theoretical prediction of low-energy photoelectron spectra of AlnNi- clusters (n = 1-13).

Mixed-metal clusters have long been studied because of their peculiar properties and how they change with cluster size, composition, and charge state and their potential roles in catalysis. The characterization of these clusters is therefore a very important issue. One of the main experimental tools for characterizing their electronic structure is photoelectron spectroscopy. Theoretical computation completes the task by fully determining the structural properties and matching the theoretical predictions to the measured spectra. We present density functional theory computations of the structural, magnetic, and electronic properties of negatively charged mixed AlnNi- clusters with up to 13 Al atoms. The lowest energy structures of the anionic clusters with up to 7 atoms are also found to be low-energy isomers of the neutral counterparts found in the literature. The 13-atom cluster is found to be a quartet and a perfect icosahedron. The predicted photoelectron spectra are also presented and can be used to interpret future experimental data. We also presented indicators that can be used to determine the potential of these systems for single-atom catalysis. These indicators point to smaller clusters to be more reactive as the gap between the Fermi energy and the center of the d-band increases with cluster size and that Ni occupies an internal site for n = 11-13. We speculate that reactivity can be enhanced if one adds an additional Ni atom. The DFT calculations were performed using the Becke exchange and Perdew-Wang/91 correlation functionals (BPW91), a DFT-optimized all-electron basis set for the aluminum atom, and the Stuttgart small core pseudopotential for the Ni atom. All of the computations used the Gaussian 03 software.

Density functional theory study of the structure, stability, magnetic properties, and (hyper)polarizability of (sub-nm) rare-earth (RE) doped gold clusters: Au5RE with RE = Sc, Y, La-Lu.

Rare-earth doped materials are of immense interest for their potential applications in linear and nonlinear photonics. There is also intense interest in sub-nanometer gold clusters due to their enhanced stability and unique optical, magnetic, and catalytic properties. To leverage their emergent properties, here we report a systematic study of the geometries, stability, electronic, magnetic, and linear and nonlinear optical properties of Au5RE (RE = Sc, Y, La-Lu) clusters using density-functional theory. Several low-energy isomers consisting of planar or non-planar configurations are identified. For most doped clusters, the non-planar configuration is energetically favored. In the case of La-, Pm-, Gd-, and Ho-doped clusters, a competition between planar and non-planar isomers is predicted. A distinct preference for the planar configuration is predicted for Au5Eu, Au5Sm, Au5Tb, Au5Tm, and Au5Yb. The distinction between the planar and non-planar configurations is highlighted by the calculated highest frequencies, with the stretching mode of the non-planar configuration predicted to be stiffer than the planar configuration. The bonding analysis reveals the dominance of the RE-d orbitals in the formation of frontier molecular orbitals, which, in turn, facilitates retaining the magnetic characteristics governed by RE-f orbitals, preventing spin-quenching of rare earths in the doped clusters. In addition, the doped clusters exhibit small energy gaps between frontier orbitals, large dipole moments, and enhanced hyperpolarizability compared to the host cluster.

Accurate ab initio spectroscopic studies of promising interstellar ethanolamine iminic precursors

Context. The detection of NH2CH2CH2OH (ethanolamine) in molecular cloud G+0.693-0.027 adds an additional player to the pre-biotic molecules discovered so far in the interstellar medium (ISM). As this molecule might be formed through condensed-phase hydrogenation steps, detecting one or more of the molecules involved might help to elucidate the chemical pathway leading to its production. Aims. The chemical path involves the formation of four chemical species. In this work, we study the energies of the isomers involved, indicate the best candidates for detection purposes, and provide the distortion constants of the most energetically favoured isomers undetected so far. Methods. We used highly accurate CCSD(T)-F12/cc-pCVTZ-F12 computations to predict the lowest energy isomers as well as their spectroscopic constants, taking corrections for core electron correlation and scalar relativity into account. Results. We studied 14 isomers. We find that the lowest energy isomer proposed in previous studies is not the actual minimum. We provide a set of rotational and distortion constants of the two new most stable isomers together with their fundamental vibrational frequencies in order to guide the search for these important astrochemical precursors of prebiotic molecules in the ISM.

Anion Photoelectron Spectroscopy and Quantum Chemistry Calculations of Gas-Phase TaSi17̅ and TaSi18̅ Clusters: Structural Determination, Bonding Characteristics, and Multiplicity of Structural Forms.

This study explores the structures and chemical bonding properties of TaSi17̅ and TaSi18̅ clusters by employing anion photoelectron spectroscopy and theoretical computations. Utilizing CALYPSO and ABCluster programs for initial structure prediction, B3LYP hybrid functional for optimization, and CCSD(T)/def2-TZVPPD level for energy calculations, the research identifies the most stable isomers of these clusters. Key findings include the identification of two coexisting low-energy isomers for TaSi17̅, exhibiting Ta-endohedral fullerene-like cage structures, and the lowest-energy structures of TaSi17̅ and TaSi18̅ anions can be considered as derived from the TaSi16̅ superatom cluster. The study enhances the understanding of group 14 element chemistry and guides the design of novel inorganic metallic compounds, potentially impacting materials science.

Enumeration, Nomenclature, and Stability Rules of Carbon Nanobelts.

With recent breakthroughs and advances in synthetic chemistry, carbon nanobelts (CNBs) have become an emerging hot topic in chemistry and materials science. Owing to their unique molecular structures, CNBs have intriguing properties with applications in synthetic materials, host-guest chemistry, optoelectronics, and so on. Although a considerable number of CNBs with diverse forms have been synthesized, no systematic nomenclature is available yet for this important family of macrocycles. Moreover, little is known about the detailed isomerism of CNBs, which, in fact, exhibits greater complexity than that of carbon nanotubes. The copious variety of CNB isomers, along with the underlying structure-property relationships, bears fundamental relevance to the ongoing design and synthesis of novel nanobelts. In this paper, we propose an elegant approach to systematically enumerate, classify, and name all possible isomers of CNBs. Besides the simplest, standard CNBs defined by chiral indices (n, m), the nonstandard CNBs (n, m, l) involve an additional winding index l. Based on extensive quantum chemical calculations, we present a comprehensive study of the relative isomer stability of CNBs containing up to 30 rings. A simple Hückel-based model with a high predictive power reveals that the relative stability of standard CNBs is governed by the π stabilization and the strain destabilization induced by the cylindrical carbon framework, and the former effect prevails over the latter. For nonstandard CNBs, a third stability factor, the H···H repulsion in the benzo[c]phenanthrene-like motifs, is also shown to be important and can be incorporated into the simple quantitative model. In general, lower-energy CNB isomers have a larger HOMO-LUMO gap, suggesting that their thermodynamic stability coincides with kinetic stability. The most stable CNB isomers determined can be considered the optimal targets for future synthesis. These results lay an initial foundation and provide a useful theoretical tool for further research on CNBs and related analogues.

Pd8 Cluster: Too Small to Melt? A BOMD Study.

The question of whether a solid-liquid phase transition occurs in small clusters poses a fundamental challenge. In this study, we attempt to elucidate this phenomenon through a thorough examination of the thermal behavior and structural stability of Pd8 clusters employing ab initio simulations. Initially, a systematic global search is carried out to identify the various isomers of the Pd8 cluster. This is accomplished by employing an ab initio basin-hopping algorithm and using the PBE/SDD scheme integrated in the Gaussian code. The resulting isomers are further refined through reoptimization using the deMon2k package. To ensure the structural firmness of the lowest-energy isomer, we calculated normal modes. The structural stability as a function of temperature is analyzed through the Born-Oppenheimer molecular dynamics (BOMD) approach. Multiple BOMD trajectories at distinct simulated temperatures are examined with data clustering analysis to determine cluster isomers. This analysis establishes a connection between the potential energy landscape and the simulated temperature. To address the question of cluster melting, canonical parallel-tempering BOMD runs are performed and analyzed with the multiple-histogram method. A broad maximum in the heat capacity curve indicates a melting transition between 500 and 600 K. To further examine this transition, the mean-squared displacement and the pair-distance distribution function are calculated. The results of these calculations confirm the existence of a solid-liquid phase transition, as indicated by the heat capacity curve.

Interaction of N2, O2 and H2 Molecules with Superalkalis.

Superalkalis (SAs) are exotic clusters having lower ionization energy than alkali atoms, which makes them strong reducing agents. In the quest for the reduction of diatomic molecules (X2) such as N2, O2, and H2 using Møller-Plesset perturbation theory (MP2), we have studied their interaction with typical superalkalis such as FLi2, OLi3, and NLi4 and calculated various parameters of the resulting SA-X2 complexes. We noticed that the SA-O2 complex and its isomers possess strong ionic interaction, which leads to the reduction of O2 to O2 - anion. On the contrary, there are both ionic and covalent interactions in SA-N2 complexes such that the lowest energy isomers are covalently bonded with no charge transfer from SA. Further, the interaction between SA and H2 leads to weakly bound complexes, which results in the adsorption of H2 molecules. The nature of interaction is found to be closely related to the electron affinity of diatomic molecules. These findings might be useful in the study of the activation, reduction, and adsorption of small molecules, which can be further explored for their possible applications.

C(P)42+: a viable planar tetracoordinate carbon species.

Tetrahedral carbon is a well-accepted structural motif. Theory has paved the way for many rule-breaking planar tetracoordinate carbon structures. Herein, we have explored the possibility of phosphorous supported carbon cluster C(P4)2+ to locate a planar tetracoordinate carbon centre (ptC). Our calculations revealed a penta-atomic planar tetracoordinate carbon (ptC) atom as a local minimum with a significant barrier for interconversion into the lowest energy isomer. The proposed species is kinetically stable and is a probable candidate for experimental detection. The stabilization of the planar structure arises due to delocalization of the p-electrons of the central carbon atom to the surrounding P4 skeleton and significant electrostatic attraction.

On the brink of self-hydration: the water heptadecamer.

For pure, neutral, isolated molecular clusters, (H2O)17 marks the transition from structures with all water molecules on the cluster surface to water self-hydration, i.e., cluster structures around one central water molecule. Getting this right with water model potentials turns out to be challenging. Even the best water potentials currently available, which reproduce collective properties very well, still deliver contradicting results for (H2O)17, when different low-energy isomers from global structure optimizations are examined. Interestingly, ab initio quantum chemistry also struggles with the only seemingly simple question if (H2O)17 is all-surface or water-centered. Hence, although the long history of water potential development may be entering its final phase, it is not quite finished yet.

Exploration of carvacrol aggregation by laser spectroscopy.

Carvacrol is an aromatic monoterpenoid found in thyme oil. Due to its implications for human health, it is important to elucidate its structure and its intramolecular interactions. We have characterised the carvacrol monomer, its complex with water, its dimer, and even its trimer in a supersonic expansion using mass-resolved laser spectroscopy techniques complemented by quantum-chemical computations. The resonance-enhanced multiphoton ionisation spectrum of the monomer features several transitions, which were assigned to the same conformer, confirmed by ion-dip infrared spectroscopy. However, a conclusive assignment of the infrared bands to one of the four conformations of carvacrol remains elusive. The experimental spectra for the monohydrated, the homodimer, and the homotrimer point to the detection of the lowest energy isomer in each case. Their structures are governed by a balance of intramolecular interactions, specifically hydrogen bonding and dispersion forces. Comparison with other similar systems demonstrates that dispersion interactions are key to the stabilisation of the aggregates, being present in all the structures. However, the hydrogen bonding is the dominant force as observed in the lowest-energy conformations.

Magic cluster sizes of cationic and anionic sodium chloride clusters explained by statistical modeling of the complete phase space.

As one of the main components of sea salt aerosols, sodium chloride is involved in numerous atmospheric processes. Gas-phase clusters are ideal models to study fundamental physical and chemical properties of sodium chloride, which are significantly affected by the cluster size. Of particular interest are magic cluster sizes, which exhibit high intensities in mass spectra. In order to understand the origin of these magic cluster sizes, quantum chemical calculations at the CCSD(T)//DFT level are performed, yielding structures and binding energies of neutral (NaCl)x, anionic (NaCl)xCl- and cationic (NaCl)xNa+ clusters up to x = 8. Our calculations show that the clusters can easily isomerize, enabling dissociation into the lowest-energy isomers of the fragments. Energetics can explain the special stability of (NaCl)4Cl-, but (NaCl)4Na+ actually offers low-lying dissociation channels, despite being a magic cluster size. Collision-induced dissociation experiments reveal that the loss of neutral clusters (NaCl)x, x = 2, 4, is in most cases more favorable than the loss of NaCl or the atomic ion, i.e. sodium chloride clusters actually fragment via the cleavage of the entire cluster, not by evaporating small cluster building blocks. This is rationalized by the calculated high stability of even-numbered neutral clusters (NaCl)x, especially x = 2, 4. Analysis of the density of states and rate constants calculated with a modified Rice-Ramsperger-Kassel-Marcus (RRKM) equation called AWATAR - considering all energetically accessible isomers of reactants and fragments - shows that entropic effects are responsible for the magic cluster character of (NaCl)4Na+. In particular, low-lying vibrational modes provide a high density of states of the near-planar cluster. Together with the small contribution of an atomic ion to the sum of states in a loose transition state for dissociation, this leads to a very small unimolecular rate constant for dissociation into (NaCl)4 and Na+, which is the lowest energy fragmentation pathway. Thus, entropic effects may override energetics for certain magic cluster sizes.

Theoretical Studies on the Isomerization Kinetics of Low-Lying Isomers of the SiC4H2 System.

The low-lying isomers of SiC4H2 are investigated to understand the kinetics of isomerization pathways using density functional theory. In our earlier work, we studied the various possible isomers (J. Phys. Chem. A, 2020, 124, 987-1002) and the chemical bonding of low-lying isomers of SiC4H2 (J. Phys. Chem. A, 2022, 126, 9366-9374). Among them, four isomers, 1-ethynyl-3-silacycloprop-1-en-3-ylidene (1), 3-silapent-1,4-diyn-3-ylidene (2), 1-silapent-1,2,3,4-tetraen-1-ylidene (4), and 1-silapent-2,4-diyn-1-ylidene (5) have already been identified in the laboratory. The previously known theoretical isomer 2-methylene-1-silabicyclo[1.1.0]but-1(3)-en-4-ylidene (3) and the newly identified unknown isomer through the present kinetic studies 5-silabicyclo[2.1.0]pent-1(4),2-dien-5-ylidene (N6) remain elusive in the laboratory to date. The isomerization pathways of the low-lying isomers of SiC4H2 are predicted through the transition state structures. Intrinsic reaction coordinate analysis identifies the minimum energy reaction pathways connecting the transition state from one isomer to another of the investigated system. The present kinetic data reveal the isomerization of global minimum energy isomer 1 to thermodynamically stable low-lying isomers, 2 and 5. Interestingly, isomer 3 interconverts to the experimentally known low-energy isomer 4, the second most thermodynamically stable isomer among them. The thermodynamic and kinetic parameters of the low-lying isomers of SiC4H2 are also documented in this work. The rate coefficient and equilibrium constant for isomerization reactions are calculated using the Rice-Ramsperger-Kassel-Marcus theory. The equilibrium constant delineates the difficulties in forming N6 and 3 through the isomerization pathways. Furthermore, ab initio molecular dynamics studies dictate the stability of low-lying isomers of SiC4H2 within the time scale of the simulation.

Diffusion quantum Monte Carlo study on magnesium clusters as large as nanoparticles

Nanoscale magnesium clusters are important potential hydrogen storage materials, and density functional theory (DFT) is mainly used for their theoretical investigation. The results of the coupled-cluster theory at the singles and doubles level with a perturbative treatment of triples [CCSD(T)] were employed previously to choose proper exchange–correlation (XC) functionals in DFT calculations for magnesium clusters, but it is too expensive to be applied to Mgn with n > 7. The diffusion Monte Carlo (DMC) method is employed in this work to study magnesium clusters up to nanosize. The error of atomization energies with DMC using single-determinant-Jastrow (SDJ) trial wavefunctions has been shown to be somewhat larger than that of CCSD(T) for many molecules. However, cohesive energies with DMC using SDJ for Mgn with n ≤ 7 are in excellent agreement with those of CCSD(T) using the aug-cc-pVQZ basis set, with a difference of less than 1 kcal/mol. DMC results are employed to investigate the performance of different XC functionals on magnesium clusters. Our results indicate that the PBE0 functional is the best XC functional for determining the lowest-energy isomer when compared with DMC results, while the RPBE functional is the best XC functional for calculating cohesive energies per atom of these magnesium clusters with a mean absolute error of 0.5 kcal/mol. These XC functionals are expected to provide reasonable results for even larger magnesium clusters.

Effect of unpaired electron number elements (Al, Cr, Mn) doping in Fe2O3 on ortho to para hydrogen conversion at 77 K

Para hydrogen (p-H2), as a low-energy spin isomer of hydrogen molecule, has an important significance in the storage and transportation of liquid hydrogen. The ortho hydrogen (o-H2) to p-H2 conversion catalyst accelerates the intrinsically slow o-H2 to p-H2 conversion and becomes an essential component of the hydrogen liquefaction process. In this paper, doped Fe2O3 with different hetero elements (X-FO, X = Al, Cr, Mn) are prepared by a simple co-precipitation method to explore the effect of different unpaired electron numbers in doping elements on catalytic performance. It is shown that the disorder in the internal crystal structure of X-FO caused by doping of hetero elements is the main reason for the increased performance of X-FO, as this leads directly to an increase in the internal magnetic disorder. The catalytic performance is enhanced with the increase in the number of unpaired electrons in doping elements. Mn-FO with high unpaired electron number exhibits optimal catalytic performance. When the hydrogen flow rate is 100 mL min−1, p-H2 content is 49.45 % after catalytic conversion at 77 K. In this paper, a simple method is adopted, which has advantages in terms of production process and cost, and provides experimental support for the large-scale production of catalysts.