Density functional study of polycarbon sulfides CnS (n = 2−9)

TMC-1, the starless core sulfur factory: Discovery of NCS, HCCS, H2CCS, H2CCCS, and C4S and detection of C5S

We report the detection of the sulfur-bearing species NCS, HCCS, H2CCS, H2CCCS, and C4S for the first time in space. These molecules were found towards TMC-1 through the observation of several lines for each species. We also report the detection of C5S for the first time in a cold cloud through the observation of five lines in the 31-50 GHz range. The derived column densities are N(NCS) = (7.8±0.6)×1011 cm-2, N(HCCS) = (6.8±0.6)×1011 cm-2, N(H2CCS) = (7.8±0.8)×1011 cm-2, N(H2CCCS) = (3.7±0.4)×1011 cm-2, N(C4S) = (3.8±0.4)×1010 cm-2, and N(C5S) = (5.0±1.0)×1010 cm-2. The observed abundance ratio between C3S and C4S is 340, that is to say a factor of approximately one hundred larger than the corresponding value for CCS and C3S. The observational results are compared with a state-of-the-art chemical model, which is only partially successful in reproducing the observed abundances. These detections underline the need to improve chemical networks dealing with S-bearing species.

Photoelectron spectroscopy and density functional calculations of C(n)S(m)(-) (n = 2-7; m = 1, 2) clusters

CnSm(-) (n = 2-7; m = 1, 2) clusters were investigated by using photoelectron spectroscopy combined with density functional theory calculations. We found that the vertical detachment energies of both CnS(-) and CnS2(-) (n = 2-7) clusters exhibit a strong odd-even alternation with an increasing number of carbon atoms: the VDEs of even-n clusters are higher than those of adjacent odd-n clusters. The most stable structures of the anionic and neutral CnS (n = 2-7) clusters are linear with the S atom locating at one end of the carbon chain except that the structure of C3S(-) is slightly bent. The ground state isomers of the anionic and neutral CnS2 (n = 2-7) clusters are all linear structures with two S atoms locating at two ends of the carbon chain. The electron affinities of the neutral CnS (n = 2, 4-7) and CnS2 (n = 2-7) clusters are determined based on the experimental adiabatic detachment energies of the corresponding anion species, because the most stable structures of the neutral clusters are similar to those of the corresponding anions.

DFT study of sulfur derivatives of cumulenes and their protonated forms of interstellar interest and calculations of dissociation energies of protonated forms (SC(CH)C(n-2)S)(+) (n = 3-8)

A theoretical study of the sulfur cumulenes SCnS (n = 3-8), CnS ( n = 1-8) and of their protonated forms (SCnS)H(+) and (CnS)H(+) that might exist in the interstellar environment, has been carried out by means of the standard B3LYP/6-311G** method. The geometries and relative energies of singlet and triplet states according to the number of carbons have been computed. Like neutral species, we have found that the ground state of the most stable protonated forms (SC(CH)Cn-2S)(+) and ((HC)Cn-1S)(+), alternates between a triplet state for the even series and a singlet state for the odd series. We provided the data needed to simulate infrared and microwave spectra (vibration frequencies, dipole moments, and rotational constants) for each protonated species (SCnS)H(+) and (CnS)H(+) and for each neutral CnS species. The computing of dissociation energies of the most stable protonated forms (SC(CH)Cn-2S)(+) (n = 3-8) has shown that the lowest values are obtained for the dissociation of compounds with an even number of carbons, in their triplet state, which produce the observed fragments CS and C3S. The dissociation of even protonated forms requires less energy than for the odd protonated forms.

Geometries and stabilities of the carbon clusters with the rhodium impurity: a computational investigation

A density functional study of the RhCn(n = 1-6) clusters with different spin states has been carried out systematically by using the B3LYP/Lan2DZ method. The equilibrium geometries associated with total energies and natural populations of RhCn (n = 1-6) clusters are calculated and presented. Stabilities and electronic properties are discussed in detail. The relative stabilities in term of the calculated fragmentation energies show that the lowest-energy RhCn clusters with rhodium atom being located at terminal of carbon chain are the linear geometries and the ground states of the RhCn clusters alternate between doublet (for n-odd members) and quartet (for n-even members) states. Furthermore, the calculated fragmentation energies of the RhCn show strong even-odd alternations: the RhCn clusters with an odd number of carbon atoms are more stable than those with an even number ones. In addition, we comment on the charge transfer and chemical bonding properties within the clusters.

Electronic spectra of linear isoelectronic clusters C2n+1S and C2n+1Cl+ (n = 0-4): an ab initio study

Structures and stabilities of linear carbon chains C2n+1S and C2n+1Cl+ (n=0-4) in their ground states have been investigated by the CCSD and B3LYP approaches. The CASSCF calculations have been used to determine geometries of selected excited states of both isoelectronic series. Linear C2n+1S cluster has a cumulenic carbon framework, whereas its isoelectronic C2n+1Cl+ has a dominant character of acetylenic structure in the vicinity of terminal Cl. The vertical excitation energies of low-lying excited states have been calculated by the CASPT2 method. Calculations show that the excitation energies have nonlinear size dependence. The 2(1)Sigma+<--X1Sigma+ transition energy in C2n+1S has a limit of 1.78 eV, as the chain size is long enough. The predicted vertical excitation energies for relatively strong 1(1)Pi<--X1Sigma+ and 2(1)Sigma+<--X1Sigma+ transitions are in reasonable agreement with available experimental values. The spin-orbit effect on the spin-forbidden transition in both series is generally small, and the enhancement of the spin-forbidden transition by spin-orbit coupling exhibits geometrical and electronic structural dependence.

Electronic spectra of heteroatom-containing isoelectronic carbon chains C2nS and C2nCl+ (n=1–5)

Structures and stabilities of carbon chains C(2n)S and C2(n)Cl+ (n=1-5) in their ground states have been investigated by the density functional theory and the coupled cluster approach using single and double substitutions. The complete active space self-consistent-field method has been used for geometry optimization of selected excited states in both series. Calculations show that both C(2n)S (n=1-5) and C2(n)Cl+ (n=3-5) have linear structures in the triplet ground state 3Sigma-, while C2Cl+ and C4Cl+ have nonlinear structures in the ground state 3A". The vertical transition energies and emission energies by the multiconfigurational second-order perturbation theory in linear clusters C(2n)S and C2(n)Cl+ exhibit similar size dependences. In comparison with the available experimental observations, the predicted excitation energies for the allowed 2 3Sigma- <--X 3Sigma- transitions have an accuracy of no more than 0.24 eV. Spin-orbit coupling configuration interaction calculations indicate that the spin-forbidden 2 1Sigma+<--X 3Sigma- transition in these species has an oscillator strength with the magnitude of 10(-4)-10(-5), and they may be observable experimentally.

Small ScCn Cyclic Clusters: A Density Functional Study of Their Structure and Stability

A theoretical study of the ScCn, ScCn+, and ScCn- (n = 1-10) cyclic clusters has been carried out employing the B3LYP density functional method. Predictions for several molecular properties that could help in their possible experimental characterization, such as equilibrium geometries, electronic structures, dipole moments, and vibrational frequencies, are reported. All ScCn cyclic clusters are predicted to have doublet ground states. For cationic clusters the ground state is alternate between singlets (n-even species) and triplets (n-odd members). In the case of anionic clusters the singlet-triplet separation is relatively small, with the singlets being favored in most cases. In general, even-odd parity effects are also observed for different properties, such as incremental binding energies, ionization energies, and electron affinities. For all neutral, cationic, and anionic clusters it is found that cyclic species are more stable than their open-chain counterparts. Therefore, cyclic structures are the most interesting possible targets for an experimental search of scandium-doped carbon clusters.

Theoretical study of the electronic structure of CnS (n=1–6) thiocumulenes

Linear sulfur-carbon chains C(n)S (n=1-6) of astronomical interest were examined by means of several theoretical methods. The three smallest compounds of the series were chosen to evaluate the performance of several computational models, including Hartree-Fock theory, density functional theory with the Becke's three parameter exchange functional and the correlation functional of Lee, Yang, and Parr (B3LYP), and electron-correlated methods (second-order Moller-Plesset perturbation method (MP2), configuration interaction method including single and double excitations (CISD), and quadratic configuration interaction method including single and double excitations (QCISD) in combination with a large variety of basis sets. The systematic comparison between the experiment and theory indicates that the B3LYP/6-311G** method can be considered suitable for the study of the electronic structures of the C(n)S compounds. The electronic ground states of the C(n)S molecules alternate between 1Sigma and 3Sigma for odd and even values of n, respectively. The B3LYP/6-311G** wave functions for these electronic ground states were analyzed by means of the atoms in molecules (AIM) and natural bond orbital (NBO) methods. Both approaches suggest that the electronic structures for the singlet and triplet compounds must be considered separately. According to the NBO method, singlet compounds can be properly represented by acetylenic structures with alternating single and triple bonds (S[triple bond]C-C[triple bond]C...). However, triplet compounds are better described by means of double bond-double bond cumulenic structures (S=C=C=C=C...) as a consequence of the average between different alpha and beta electronic densities. AIM delocalization indexes and NBO interactions between localized orbitals also indicate that these structures are strongly pi delocalized. Finally, the different singlet and triplet structures proposed provide a consistent explanation for the geometries, dipole moments, and spin-density values of the C(n)S compounds studied.

On the Competition between Linear and Cyclic Isomers in Second-Row Dicarbides

Second-row dicarbides C(2)X (X = Na-Cl) are investigated with quantum mechanical techniques. The cyclic-linear competition in these systems is studied, and the bonding scheme for these compounds is discussed in terms of the topological analysis of the electronic density. C(2)Na, C(2)Mg, C(2)Al, and C(2)Si are found to prefer a C(2)(v)-symmetric arrangement corresponding to a T-shape structure. On the other hand, for C(2)P, C(2)S, and C(2)Cl the linear isomer is predicted to be the ground state. A detailed analysis of the variation of the electronic energy and orbital energies with the geometry has been carried out. A simple theoretical model, taking into account the main interactions between the valence orbitals of both fragments, the X atom and the C(2) molecule, allows an interpretation of the main features of these compounds.