Gas-phase C60Hn+q (n = 0-4, q = 0,1) fullerenes and fulleranes: spectroscopic simulations shed light on cosmic molecular structures.

The discovery of C60, C60+, and C70 in the interstellar medium has ignited a profound interest in the astrochemistry of fullerene and related systems. In particular, the presence of diffuse interstellar bands and their association with C60+ has led to the hypothesis that hydrogenated derivatives, known as fulleranes, may also exist in the interstellar medium and contribute to these bands. In this study, we systematically investigated the structural and spectroscopic properties of C60Hn+q (n = 0-4, q = 0,1) using an automated global minimum search and density functional theory calculations. Our results revealed novel global minimum structures for C60H2 and C60H4, distinct from previous reports. Notably, all hydrogenated fullerenes exhibited lower ionization potentials and higher proton affinities compared to C60. From an astrochemical perspective, our results exposed the challenges in establishing definitive spectroscopic criteria for detecting fulleranes using mid-infrared and UV-Vis spectroscopies. However, we successfully identified distinct electronic transitions in the near-infrared range that serve as distinctive signatures of cationic fulleranes. We strongly advocate for further high-resolution experimental studies to fully explore the potential of these transitions for the interstellar detection of fulleranes.

Revisiting the Spectrum of Co(CN)63-: The Role of Correlation, Solvation, and Vibronic and Spin-Orbit Couplings.

In the present work, we revisit the spectrum of the hexacyanocobaltate(III) ion, [Co(CN)6]3-, which has been considered a prototype complex in the coordination chemistry, with modern quantum chemistry methods. The main features have been describing by revealing the role of different effects, such as vibronic coupling, solvation and spin-orbit coupling. The UV-vis spectrum is composed by two bands (1A1g → 1T1g and 1A1g → 1T2g), characterized by singlet-singlet metal-centered transitions, and a more intense third one, characterized by charge transfer transition. There is also a small band shoulder. The first two are symmetry-forbidden transitions in the Oh group. Their intensity can only be explained by a vibronic coupling mechanism. For the band shoulder, additional to vibronic coupling, spin-orbit coupling is also necessary, since the transition is characterized as singlet to triplet, 1A1g → 3T1g.

Structure, Stability, and Spectroscopic Properties of Small Acetonitrile Cation Clusters.

Ionization and fragmentation pathways induced by ionizing agents are key to understanding the formation of complex molecules in astrophysical environments. Acetonitrile (CH3CN), the simplest organic nitrile, is an important molecule present in the interstellar medium. In this work, DFT and MP2 calculations were performed in order to obtain the low energy structures of the most relevant cations formed from electron-stimulated ion desorption of CH3CN ices. Selected reaction pathways and spectroscopic properties were also calculated. Our results indicate that the most stable acetonitrile cation structure is CH2CNH+ and that hydrogenation can occur successively without isomerization steps until its complete saturation. Moreover, the stability of distinct cluster families formed from the interaction of acetonitrile with small fragments, such as CHn+, C2Hn+, and CHnCNH+, is discussed in terms of their respective binding energies. Some of these molecular clusters are stabilized by hydrogen bonds, leading to species whose infrared features are characterized by a strong redshift of the N-H stretching mode. Finally, the rotational spectra of CH3CN and protonated acetonitrile, CH3CNH+, were simulated using distinct computational protocols based on DFT, MP2, and CCSD(T) considering centrifugal distortion, vibrational-rotational coupling, and vibrational anharmonicity corrections. By adopting an empirical scaling procedure for calculating spectroscopic parameters, we were able to estimate the rotational frequencies of CH3CNH+ with an expected average error below 1 MHz for J values up to 10.

Fragment and cluster ions from gaseous and condensed pyridine produced under electron impact.

We report on direct measurement of all major ion-fragments and cluster-ions formed during high-energy electron impact of 2 keV on gaseous and condensed-phase pyridine. The ion-fragments of the parent pyridine cation are discussed in groups according to the number of atoms from the aromatic ring. The ion yield distributions within these groups show significant shifts towards higher masses for condensed pyridine compared to gaseous pyridine due to hydrogen migration. A wide spectrum of desorbed hydrogenated fragment-ions and ionic clusters with masses up to 320 u are observed for pyridine. The ion yields for the protonated parent molecule (C5H5NH+), the dehydrogenated dimer (C10H9N2+) and the dehydrogenated trimer (C15H12N3+) depend on the mass of the desorbing ionic clusters. The strongest cluster signals are assigned to binding between the parent cation and subunits of the pyridine molecule. Quantum-chemical calculations reveal that the formation of a bond between the pyridine molecules and a carbenium ion is crucial for the stability of selected cluster ions.

Total absorption spectrum of benzene aggregates obtained from two different approaches.

CONTEXT: The study of molecular aggregation effects on the electronic spectrum is essential for the development of optoelectronic devices. However, investigating the entire valence absorption spectrum of aggregates using quantum mechanical methods is a challenging task. In this work, we perform systematic simulations of the absorption spectrum of benzene molecular clusters up to 35 eV applying two approaches based on time-dependent density functional theory. The results show that depending on the dimer packing, different energy shifts occur for the symmetry allowed [Formula: see text] transition, in comparison to the monomer. The transition intensity increases for the band around 6 eV for larger aggregates from the monomer to dimers and tetramer, indicating the occurrence of the symmetry forbidden (in [Formula: see text] point group) [Formula: see text] [Formula: see text] transition. The benzene crystal exhibits a large redshift following the experimental spectrum. Also, the continuum regions of all spectra show a good agreement with the experiments both in gas and solid phases. METHODS: Geometry optimization of the monomer was carried out with Gaussian 09 software using the PBE0/def2-TZVP level of theory. We used dimers and tetramer molecular geometries extracted from the experimental crystal structure. The absorption spectra were directly obtained by the Liouville-Lanczos TDDFT approach with plane waves basis set or indirectly by TDDFT pseudo-spectra calculated in a [Formula: see text] basis followed by analytic continuation procedure to obtain complex polarizability. The former is available at Quantum ESPRESSO, and the latter was calculated using Gaussian 09 with the post-processing performed with a code previously developed in our group.

Comparison among several vibronic coupling methods.

A comparison of four approaches to account the vibronic coupling in photoabsorption is performed. The methods considered are nuclear ensemble (NE), direct vibronic coupling (DVC), adiabatic Hessian (AH), and vertical gradient (VG). The case study is the symmetry-forbidden [Formula: see text] [Formula: see text]A[Formula: see text] [Formula: see text] [Formula: see text] [Formula: see text]A[Formula: see text] (n [Formula: see text] [Formula: see text]) transition in formaldehyde. Being forbidden in the equilibrium geometry, this transition is entirely induced by vibronic coupling and constitutes an appropriate case to study the performance of different methods. From DVC, it is found that mode 1 (C=O out-of-plane bending) is the most inducing, followed by mode 6 (in-plane C-H asymmetric stretching) and finally by mode 2 (in-plane C-H asymmetric bending). We were able to correlate 17 out of 20 structures obtained from NE with these modes, showing that these two methods, although different in principle, give comparable results. The simulated spectra were obtained for all methods and compared, and each one has its own advantage. In what concerns the transition studied, NE gives the best description of the spectrum, DVC is the only one that easily gives an absolute value for OOS, and AH and VG are the computationally less expensive methods. From the latter two, VG is the less demanding on computational grounds, since it does not require the excited state Hessian.

Potential curves for inner-shell states of CO calculated at multiconfigurational self-consistent field level.

A general strategy to calculate potential curves at multiconfigurational self-consistent field (MCSCF) level for inner-shell states is reported in this paper. Convergence is commonly very tough for inner-shell states, especially at this level of calculation, due to the problem of variational collapse of the inner-shell wave function to the ground or to a low-lying excited state. The present method allows to avoid this drawback by a sequence of constrained optimization in the orbital mixing step. The specific states studied are that resulting from transitions X (1)Σ(+) → (C 1s(-1) π(∗)) (1,3)Π of CO. Accurate values are achieved for transition energies and vibrational splittings. A comparison is made with other approach, i.e., inner-shell CI based on a MCSCF wave function optimized for ground or low-lying excited states. This last approach is shown to fail in describing the whole potential curve.

The problem of hole localization in inner-shell states of N2 and CO2 revisited with complete active space self-consistent field approach.

Potential energy curves for inner-shell states of nitrogen and carbon dioxide molecules are calculated by inner-shell complete active space self-consistent field (CASSCF) method, which is a protocol, recently proposed, to obtain specifically converged inner-shell states at multiconfigurational level. This is possible since the collapse of the wave function to a low-lying state is avoided by a sequence of constrained optimization in the orbital mixing step. The problem of localization of K-shell states is revisited by calculating their energies at CASSCF level based on both localized and delocalized orbitals. The localized basis presents the best results at this level of calculation. Transition energies are also calculated by perturbation theory, by taking the above mentioned MCSCF function as zeroth order wave function. Values for transition energy are in fairly good agreement with experimental ones. Bond dissociation energies for N(2) are considerably high, which means that these states are strongly bound. Potential curves along ground state normal modes of CO(2) indicate the occurrence of Renner-Teller effect in inner-shell states.

Theoretical investigation on the stability of negatively charged formic acid clusters.

Recent experimental results on negatively charged formic acid clusters generated by the impact of (252)Cf fission fragments on icy formic acid target are compared to quantum mechanical calculations. Structures for the clusters series, (HCOOH)nOH(-), where 2 < or = n < or = 4, are proposed based on ab initio electronic structure methods. The results show that cluster growth does not have a regular pattern of nucleation. A stability analysis was performed considering the commonly defined stability function. Temporal behavior of the clusters was evaluated by Born-Oppenheimer molecular dynamics to check the mechanism that provides cluster stability. The evaluated temporal profiles indicate the importance of hydrogen atom migration between the formic acid moieties in maintaining the stability of the structures and the water formation due to hydrogen abstraction by the hydroxyl approach.

Rotational spectrum simulations of asymmetric tops in an astrochemical context.

Rotational spectroscopy plays a major role in the field of observational astrochemistry, enabling the detection of more than 200 species including a plethora of complex organic molecules in different space environments. Those line detections allow correctly determining the sources and physical properties, as well as exploring their morphology, evolutionary stage, and chemical evolution pathways. In this context, quantum chemistry is a powerful tool to the investigation of the molecular inventory of astrophysical environments, guiding laboratory experiments and assisting in both line assignments and extrapolation of the experimental data to unexplored frequency ranges. In the present work, we start by briefly reviewing the rotational model Hamiltonian for asymmetric tops beyond the rigid-rotor approximation, including rotational-vibrational, centrifugal, and anharmonic effects. Then, aiming at further contributing to the recording and analysis of laboratory microwave spectroscopy by means of accessible, less demanding quantum chemical methods, we performed density functional theory (DFT) calculations of the spectroscopic parameters of astrochemically relevant species, followed by their rotational spectrum simulations. Furthermore, dispersion-correction effects combined with different functionals were also investigated. Case studies are the asymmetric tops H2CO, H2CS, c-HCOOH, t-HCOOH, and HNCO. Spectroscopic parameter predictions were overall very close to experiment, with mean percentage errors smaller than 1% for zeroth order and [Formula: see text] for first-order constants. We discuss the implications and impacts of those constants on spectrum simulations, and compare line-frequency predictions at millimeter wavelengths. Moreover, theoretical spectroscopic parameters of c-HCOOH and HNCO are introduced for the first time in this work.

Acrylic acid hydrodeoxygenation reaction mechanism over molybdenum carbide studied by DFT calculations.

Platinum- and palladium-based catalysts are commonly used in hydrogenation reactions, but they present a great disadvantage of being quite expensive. In most cases, they can be substituted by cheaper alternative catalysts formed by transition metal carbides, such as molybdenum carbide (Mo2C). Among the reactions that can be catalyzed by Mo2C, hydrodeoxygenation (HDO) presents a great technological interest, especially in biofuel production. Nonetheless, the selectivity of carbides in HDO reactions of fatty acids is not well understood yet. In the present work, the reaction mechanism of the acrylic acid HDO over Mo2C, a fatty acid model molecule, was studied by density functional theory (DFT), with Perdew-Burke-Ernzerhof (PBE) functional and periodic boundary conditions. A global mechanism is proposed, divided in four steps, from acrylic acid to propane. In the first reaction step, decomposition by C-OH bond cleavage, with 24 kcal mol- 1 of activation energy, dominates over C=C and C=O hydrogenation. This result is in line with the absence of propanoic acid among the products and the formation of acrolein, as shown in an experimental work previously published. The proposed global mechanism is in fair agreement with the experimental findings. The main product is propane, which has the same number of carbon atoms of the reactant. This mechanism can be viewed as a model for HDO of any fatty acid catalyzed by Mo2C, since acrylic acid has the minimal structural features of fatty acids, i.e., a carboxyl group and a C=C double bond. Graphical Abstract HDO over Mo2C provides a product with same carbon atoms number of the reactant.

Characterization by 27Al NMR, X-ray absorption spectroscopy, and density functional theory techniques of the species responsible for benzene hydrogenation in Y zeolite-supported carburized molybdenum catalysts.

Carburized molybdenum catalysts supported on a dealuminated NaH-Y zeolite were prepared by carburization under a 20% methane in hydrogen flow of two precursors obtained by adsorption of molybdenum hexacarbonyl, one containing 5 wt % and the other 10 wt % Mo, and a third one was prepared by impregnation with aqueous ammonium heptamolybdate, containing 5 wt % Mo. The three catalysts displayed very distinct behaviors in the benzene hydrogenation reaction at atmospheric pressure and 363 K. By using XANES spectroscopy at the molybdenum L edge, EXAFS and XANES spectroscopy at the molybdenum K edge, and 27Al solid-state NMR spectroscopy, it was shown that different carburized molybdenum species exist in each sample. In the catalyst containing 10 wt % Mo, formation of molybdenum carbide nanoparticles was observed, with an estimated diameter of 1.8 nm. In the catalyst containing 5 wt % Mo and prepared by carburization of adsorbed molybdenum hexacarbonyl, formation of molybdenum oxycarbide dimers is proposed. In the latter case, density functional theory calculations have led to a dimer structure which is compatible with EXAFS results. In the catalyst prepared by impregnation with ammonium heptamolybdate solution followed by carburization, the molybdenum seems to interact with extraframework alumina to produce highly disordered mixed molybdenum-aluminum oxycarbides.

Spin-orbit splitting for inner-shell 2p states.

A strategy to calculate spin-orbit splitting for inner-shell transitions at an ab initio level is presented. The initial wave function is calculated for a spinless Hamiltonian at a multiconfigurational level, with just a few configurations, followed by multireference configuration interaction in order to establish a set of singlet and triplet states at 2p excitation edge. Then, the full Breit-Pauli Hamiltonian is formed and diagonalized on this state basis. The spin-orbit splitting is determined by a graphical procedure depending on the intensity of the transition from ground state. The specific states studied are those originating from 2p transitions in argon, HCl, H2S, and PH3.

Transition energy and potential energy curves for ionized inner-shell states of CO, O2 and N 2 calculated by several inner-shell multiconfigurational approaches.

Potential energy curves and inner-shell ionization energies of carbon monoxide, oxygen and nitrogen molecules were calculated using several forms of the inner-shell multiconfigurational self-consistent field (IS-MCSCF) method-a recently proposed protocol to obtain specifically converged inner-shell states at this level. The particular forms of the IS-MCSCF method designated IS-GVB-PP, IS-FVBL and IS-CASSCF stand for perfect pairing generalized valence bond, full valence bond-like MCSCF and complete active space self consistent field, respectively. A comparison of these different versions of the IS-MCSCF method was carried out for the first time. The results indicate that inner-shell states are described accurately even for the simplest version of the method (IS-GVB-PP). Dynamic correlation was recovered by multireference configuration interaction or multireference perturbation theory. For molecules not having equivalent atoms, all methods led to comparable and accurate transition energies. For molecules with equivalent atoms, the most accurate results were obtained by multireference perturbation theory. Scalar relativistic effects were accounted for using the Douglas-Kroll-Hess Hamiltonian.

Ab initio simulation of changes in geometry, electronic structure, and Gibbs free energy caused by dehydration of hydrotalcites containing Cl⁻ and CO3²â» counteranions.

This ab initio study was performed to better understand the correlation between intercalated water molecules and layered double hydroxides (LDH), as well as the changes that occur by the dehydration process of Zn-Al hydrotalcite-like compounds containing Cl⁻ and CO3²â» counterions. We have verified that the strong interaction among intercalated water molecules, cointercalated anions, and OH groups from hydroxyl layers is reflected in the thermal stability of these compounds. The Zn(2/3)Al(1/3)(OH)2Cl(1/3)·2/3H2O hydrotalcite loses all the intercalated water molecules around 125 °C, while the Zn(2/3)Al(1/3)(OH)2(CO3)(1/6)·4/6H2O compound dehydrates at about 175 °C. These values are in good agreement with experimental data. The interlayer interactions were discussed on the basis of electron density difference analyses. Our calculation shows that the electron density in the interlayer region decreases during the dehydration process, inducing the migration of the Cl⁻ anion and the displacement of the hydroxyl layer from adjacent layers. Changes in these compound structures occur to recover part of the hydrogen bonds broken due to the removal of water molecules. It was observed that the chloride ion had initially a lower Löwdin charge (Cl(-0.43)), which has increased its absolute value (Cl(-0.58)) after the water molecules removal, while the charges on carbonate ions remain invariant, leading to the conclusion that the Cl⁻ anion can be more influenced by the amount of water molecules in the interlayer space than the CO3²â» anion in hydrotalcite-like compounds.

Intensity of d-d symmetry-forbidden electronic transition in Cr(CO)6.

Absolute absorption intensities (oscillator strengths) are calculated for the d-d symmetry-forbidden transition in hexacarbonyl chromium. The vibronic coupling mechanism is taken into account in a way that represents an alternative to the traditional perturbative approach of Herzberg and Teller. In the so-called direct method, the electronic transition moment is directly expanded in a power series of the vibrational normal coordinates of suitable symmetry. In the present case, i.e., d-d ligand field transitions, or more specifically (1)A(1g) --> (1)T(1g) and (1)A(1g) --> (1)T(2g) transitions, the dipole selection rule is broken by vibronic interaction induced by normal modes that transform like T(1u) and T(2u) representations of the O(h) group. An analysis of the relative importance of normal modes in promoting electronic transitions is carried out.