Structure, Ionization, and Fragmentation of Neutral and Positively Charged Hydrogenated Carbon Clusters: CnHmq+ (n = 1–5, m = 1–4, q = 0–3)

Why an integrated approach between search algorithms and chemical intuition is necessary?

Though search algorithms are appropriate tools for identifying low-energy isomers, fixing several constraints seems to be a fundamental prerequisite to successfully running any structural search program. This causes some potential setbacks as far as identifying all possible isomers, close to the lowest-energy isomer, for any elemental composition. The number of explored candidates, the choice of method, basis set, and availability of CPU time needed to analyze the various initial test structures become necessary restrictions in resolving the issues of structural isomerism reasonably. While one could arrive at new structures through chemical intuition, reproducing or achieving those exact same structures requires increasing the number of variables in any given program, which causes further constraints in exploring the potential energy surface in a reasonable amount of time. Thus, it is emphasized here that an integrated approach between search algorithms and chemical intuition is necessary by taking the C₁₂O₂Mg₂ system as an example. Our initial search through the AUTOMATON program yielded 1450 different geometries. However, through chemical intuition, we found eighteen new geometries within 40.0 kcal mol^-1 at the PBE0-D3/def2-TZVP level. These results indirectly emphasize that an integrated approach between search algorithms and chemical intuition is necessary to further our knowledge in chemical space for any given elemental composition.

From Molecules with a Planar Tetracoordinate Carbon to an Astronomically Known C5H2 Carbene

Ethynylcyclopropenylidene (2), an isomer of C₅H₂, is a known molecule in the laboratory and has recently been identified in Taurus Molecular Cloud-1 (TMC-1). Using high-level coupled-cluster methods up to the CCSDT(Q)/CBS level of theory, it is shown that two isomers of C₅H₂ with a planar tetracoordinate carbon (ptC) atom, (SP-4)-spiro[2.2]pent-1,4-dien-1,4-diyl (11) and (SP-4)-spiro[2.2]pent-1,4-dien-1,5-diyl (13), serve as the reactive intermediates for the formation of 2. Here, a theoretical connection has been established between molecules containing ptC atoms (11 and 13) and a molecule (2) that is present nearly 430 light years away, thus providing evidence for the existence of ptC species in the interstellar medium. The reaction pathways connecting the transition states and the reactants and products have been confirmed by intrinsic reaction coordinate calculations at the CCSDT(Q)/CBS//B3LYP-D3BJ/cc-pVTZ level. While isomer 11 is non-polar (μ = 0), isomers 2 and 13 are polar, with dipole moment values of 3.52 and 5.17 Debye at the CCSD(T)/cc-pVTZ level. Therefore, 13 is also a suitable candidate for both laboratory and radioastronomical studies.

Theoretical Studies of Two Key Low-Lying Carbenes of C5H2 Missing in the Laboratory

The equilibrium geometries and spectroscopic properties of two key singlet carbenes, buta-1,3-diynylcarbene (6) and 2-methylenebicyclo[1.1.0]but-1(3)-en-4-ylidene (9), which have not been experimentally observed to date, are investigated using high-level coupled-cluster (CC) methods. The current theoretical study necessitates new experimental data on C₅H₂ isomers considering the relevance of these molecules to interstellar chemistry. Bent-pentadiynylidene (4) has been missing in the laboratory and the prime focus of our earlier theoretical work. The present theoretical study indicates that isomers 6 and 9 are also viable experimental targets. Apart from ethynylcyclopropenylidene (2), pentatetraenylidene (3), ethynylpropadienylidene (5), and 3-(didehydrovinylidene)cyclopropene (8), which are identified by Fourier transform microwave spectroscopy, the dipole moments of elusive 4, 6, and 9 are also nonzero (μ ≠ 0). The relative energies of these isomers, calculated at the CCSDT(Q)/CBS level of theory, with respect to linear triplet pentadiynylidene (1) reveal that they all lie within 25.1 kcal mol^-1. Therefore, geometric, energetic, aromatic, and spectroscopic parameters are reported here, which may assist the efforts of molecular spectroscopists in the future. Anharmonic vibrational calculations on isomers 6 and 9 indicate that the former is loosely bound and would be challenging to be detected experimentally. Among the undetected carbenes, 9 may be considered as a potential target molecule considering its higher polarity and aromatic nature.

Fully versus constrained statistical fragmentation of carbon clusters and their heteronuclear derivatives

The Microcanonical Metropolis Monte Carlo (MMMC) method has been shown to describe reasonably well fragmentation of clusters composed of identical atomic species. However, this is not so clear in the case of heteronuclear clusters as some regions of phase space might be inaccessible due to the different mobility of the different atomic species, the existence of large isomerization barriers, or the quite different chemical nature of the possible intermediate species. In this paper, we introduce a constrained statistical model that extends the range of applicability of the MMMC method to such mixed clusters. The method is applied to describe fragmentation of isolated clusters with high, moderate, and no heteronuclear character, namely, C_nH_m, C_nN, and C_n clusters for which experimental fragmentation branching ratios are available in the literature. We show that the constrained statistical model describes fairly well fragmentation of C_nH_m clusters in contrast with the poor description provided by the fully statistical model. The latter model, however, works pretty well for both C_n and C_nN clusters, thus showing that the ultimate reason for this discrepancy is the inability of the MMMC method to selectively explore the whole phase space. This conclusion has driven us to predict the fragmentation patterns of the C₄N cluster for which experiments are not yet available.

Disentangling the complex spectrum of the ethynyl cation

The ethynyl cation, C2H+, is of great importance in astrophysical media and in combustion. It is involved in the formation of larger organic compounds and in their decomposition mechanisms. Here, we investigate the low-lying electronic states of this cation using pure ab initio methodologies. The evolution of its potential energy surfaces along the stretching and bending coordinates reveals a high density of electronic states that favours mutual interactions and the mixing of wavefunctions. The ground state is of 3Π space symmetry and the lowest singlet state (1Π) is found to be a quasi-linear-quasi-linear Renner-Teller system. Our work suggests that the (spin-)rovibronic spectrum of such a molecular system is complicated, because of the contributions of multiple couplings, including Renner-Teller, vibronic and spin-orbit. We also deduced the adiabatic ionization energy of the ethynyl radical, in good agreement with recent measurements. In summary, our work shows that the ethynyl cation, in spite of its small size, still represents a challenging molecular problem to be solved.

Ion-induced molecular growth in clusters of small hydrocarbon chains

We report on studies of collisions between 3 keV Ar⁺ projectile ions and neutral targets of isolated 1,3-butadiene (C₄H₆) molecules and cold, loosely bound clusters of these molecules. We identify molecular growth processes within the molecular clusters that appears to be driven by knockout processes and that could result in the formation of (aromatic) ring structures. These types of reactions are not unique to specific projectile ions and target molecules, but will occur whenever atoms or ions with suitable masses and kinetic energies collide with aggregates of matter, such as carbonaceous grains in the interstellar medium or aerosol nanoparticles in the atmosphere.

M3C: A Computational Approach To Describe Statistical Fragmentation of Excited Molecules and Clusters

The Microcanonical Metropolis Monte Carlo method, based on a random sampling of the density of states, is revisited for the study of molecular fragmentation in the gas phase (isolated molecules, atomic and molecular clusters, complex biomolecules, etc.). A random walk or uniform random sampling in the configurational space (atomic positions) and a uniform random sampling of the relative orientation, vibrational energy, and chemical composition of the fragments is used to estimate the density of states of the system, which is continuously updated as the random sampling populates individual states. The validity and usefulness of the method is demonstrated by applying it to evaluate the caloric curve of a weakly bound rare gas cluster (Ar₁₃), to interpret the fragmentation of highly excited small neutral and singly positively charged carbon clusters (C_n, n = 5,7,9 and C_n⁺, n = 4,5) and to simulate the mass spectrum of the acetylene molecule (C₂H₂).

Communication: On the first ionization threshold of the C2H radical

The slow photoelectron spectrum of the ethynyl radical has been recorded for the first time by using the DESIRS beamline of the SOLEIL synchrotron facility. Ethynyl was generated using a microwave discharge flow tube. The observation of the X⁺Π3←XΣ+2 transition allowed the first direct measurement of the adiabatic ionization threshold of this radical (E_I = 11.641(5) eV). The experimental results are supported by ab initio calculations. Our preliminary investigation of the cationic ground state potential energy surfaces predicts a non-negligible Renner-Teller effect which has not been discussed previously.

On the absolute photoionization cross section and dissociative photoionization of cyclopropenylidene

We report the determination of the absolute photoionization cross section of cyclopropenylidene, c-C₃H₂, and the heat of formation of the C₃H radical and ion derived by the dissociative ionization of the carbene.