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.     

Photoelectron spectroscopy and density functional calculations of CuSi(n)- (n = 4-18) clusters

We conducted a combined anion photoelectron spectroscopy and density functional theory study on the structural evolution of copper-doped silicon clusters, CuSi(n)(-) (n = 4-18). Based on the comparison between the experiments and theoretical calculations, CuSi(12)(-) is suggested to be the smallest fully endohedral cluster. The low-lying isomers of CuSi(n)(-) with n ≥ 12 are dominated by endohedral structures, those of CuSi(n)(-) with n < 12 are dominated by exohedral structures. The most stable structure of CuSi(12)(-) is a double-chair endohedral structure with the copper atom sandwiched between two chair-style Si(6) rings or, in another word, encapsulated in a distorted Si(12) hexagonal prism cage. CuSi(14)(-) has an interesting C(3h) symmetry structure, in which the Si(14) cage is composed by three four-membered rings and six five-membered rings.     

On the nature of the low-lying singlet states of 4-(Dimethyl-amino)benzonitrile

4-(N,N-Dimethyl-amino)benzonitrile (DMABN) is a prototype molecule for dual fluorescence. The anomalous emission has been attributed to an intramolecular charge-transfer (ICT) state, and the structure of the latter is still subject to some controversy. We applied a recently developed analytical gradient code for the approximate coupled-cluster singles-and-doubles method CC2 in combination with accurate basis sets to address this problem. Fully optimized excited state structures are presented for the ICT state and the so-called locally excited state, and recent transient IR and Raman measurements on the excited states are interpreted by means of calculated harmonic frequencies. Strong evidence is found for an electronic decoupling of the phenyl and the dimethyl-amino moiety, resulting in a minimum structure for the ICT state with a twisted geometry. In contrast to previous findings, the structure of this state is, at least in the gas phase, not C(2v) symmetric but distorted towards C(s) symmetry. The distortion coordinate is a pyramidalization of the phenyl carbon atom carrying the dimethyl-amino group. The results from the CC2 model are supported by single-point calculations using more elaborate coupled-cluster models (CCSD, CCSDR(3)) and by CASSCF calculations.     

Structure and vibrations of BenOn (n=3-10) clusters

The structure and harmonic vibrations of Be(n)O(n) (n=3-10) clusters have been investigated using density functional theory. All structures are found to be cumulenic D(nh) rings (equal bonds, alternating angles), with one intense out of plane mode and three infrared-active degenerate modes, of which the highest one is extremely intense and asymptotically increases to 1597 cm(-1) for n=10. Comparisons with C(2n) clusters and B(n)N(n) clusters, the structure and bonding type for the Be(n)O(n) clusters are consistent with those of the C(2n) (n=3, 5, 7, ...) clusters and the B(n)N(n) clusters.     

Structure and vibrations of AlnNn (n = 3-9) clusters

The structure and harmonic vibrations of Al(n)N(n) (n = 3-9) clusters have been investigated using the B3LYP (Becke 3-parameter-Lee-Yang-Parr) density functional theory. All structures are found to be cumulenic D(nh) rings (equal bonds, alternating angles), with one intense out-of-plane mode and three infrared-active degenerate modes, of which the highest one is extremely intense and asymptotically increases to 1217 cm(-1) for n = 9. Comparisons with C2n clusters and B(n)N(n) clusters, the structure and bonding type for the Al(n)N(n) clusters are consistent with those of the C2n (n = 3, 5, 7, ...) clusters and the B(n)N(n) clusters.     

The Nature of One‐Dimensional Carbon: Polyynic versus Cumulenic

A question of both fundamental as well as practical importance is the nature of one‐dimensional carbon, in particular whether a one‐dimensional carbon allotrope is polyynic or cumulenic, that is, whether bond‐length alternation occurs or not. By combining the concept of aromaticity and antiaromaticity with the rule of Peierls distortion, the occurrence and magnitude of bond‐length alternation in carbon chains with periodic boundary conditions and corresponding carbon rings as a function of the chain or ring length can be explained. The electronic properties of one‐dimensional carbon depend crucially on the bond‐length alternation. Whereas it is generally accepted that carbon chains in the limit of infinite length have a polyynic structure at the minimum of the potential energy surface with bond‐length alternation, we show here that zero‐point vibrations lead to an effective equalization of all carbon–carbon bond lengths and thus to a cumulenic structure.     

H2O nucleation around Au+

First principles electronic structure calculations have been carried out to investigate the ground state geometry, electronic structure, and the binding energy of [Au(H2O)n]+ clusters containing up to 10 H2O molecules. It is shown that the first coordination shell of Au+ contains two H2O molecules forming a H2O-Au+-H2O structure with C2 symmetry. Subsequent H2O molecules bind to the previous H2O molecules forming stable and fairly rigid rings, each composed of 4 H2O molecules, and leading to a dumbbell structure at [Au(H2O)8]+. The 9th and the 10th H2O molecules occupy locations above the Au+ cation mainly bonded to one H2O from each ring, leading to structures where the side rings are partially distorted and forming structures that resemble droplet formation around the Au+ cation. The investigations highlight quantum effects in nucleation at small sizes and provide a microscopic understanding of the observed incremental binding energy deduced from collision induced dissociation that indicates that [Au(H2O)n]+ clusters with 7-10 H2O molecules have comparable binding energy. The charge on the Au+ is shown to migrate to the outside H2O molecules, suggesting an interesting screening phenomenon.     

一种新型包含平面四配位碳二硼有机化合物的理论研究

采用B3LYP/6-311+G**方法, 研究了一种新型的包含平面四配位碳(ptC)二硼有机化合物C9B2H6的结构、稳定性和振动频率. 计算结果表明, C9B2H6结构的稳定性和两个硼原子的位置有很大关系, 硼原子起给予σ电子和接受π电子的作用. 在C9B2H6的15个异构体中, 最稳定的结构是具有C2v对称性的异构体(1,5), 在异构体(1,5)中, 两个硼原子位于同一个六元环中且与ptC相邻. 而且占据的π轨道说明异构体(1,5)具有10个π电子, 满足4n+2规则. 计算的核独立化学位移(NICS)值显示异构体(1,5)强的芳香性位于C9B2H6的两个三元环而不是两个六元环上.     

Accurate benchmark calculation of the reaction barrier height for hydrogen abstraction by the hydroperoxyl radical from methane. Implications for C(n)H(2n+2) where n = 2 --> 4

The CH4 + HO2(*) reaction is studied by using explicitly correlated coupled-cluster theory with singles and doubles (CCSD-R12) in a large 19s14p8d6f4g3h basis (9s6p4d3f for H) to approach the basis-set limit at the coupled-cluster singles-doubles level. A correction for connected triple excitations is obtained from the conventional CCSD(T) coupled-cluster approach in the correlation-consistent quintuple-zeta basis (cc-pV5Z). The highly accurate results for the methane reaction are used to calibrate the calculations of the hydroperoxyl-radical hydrogen abstraction from other alkanes. For the alkanes C(n)H(2n+2) with n = 2 --> 4, the reactions are investigated at the CCSD(T) level in the correlation-consistent triple-zeta (cc-pVTZ) basis. The results are adjusted to the benchmark methane reaction and compared with those from other approaches that are commonly used in the field such as CBS-QB3, CBS-APNO, and density functional theory. Rate constants are computed in the framework of transition state theory, and the results are compared with previous values available.     

The structure of small gallium nitride clusters

The geometry, vibrational frequencies and stability of the structural isomers of small gallium nitride clusters (n = 2–4) have been investigated using density functional theory. The lowest energy structures are cyclic. The ground electronic state of the cyclic forms for n > 2 is the singlet state. All of the cyclic structures have D_nh symmetry. The caged structures for Ga₄N₄ lie higher in energy than the planar cumulenic monocyclic ring. The Ga‐N bond dominates the structures for many isomers, so that one dissociation channel is loss of a GaN monomer. However, unlike the corresponding boron and aluminum clusters, dissociation into larger fragments is energetically favored. The structural properties of the gallium nitride clusters are similar to those of the analogous AIN (and BN) clusters. © 2000 John Wiley & Sons, Inc. Heteroatom Chem 11:281–286, 2000     

叠氮二氢硼多聚体结构和性质的理论研究(英文)

本文采用DFT-B3LYP方法,以不同基组对叠氮二氢硼多聚体(H2BN3)n (n=1-4)进行计算研究.二聚体(H2BN3)2(C2h对称性)中含B2N2平面四元环结构.船式(Cs对称性)和椅式(C3v对称性)三聚体(H2BN3)3的结合能相近(-122 和 -126 kJ·mol-1),其中均含B3N3六元环结构.拥有B4N4八元环结构的四个四聚体的结合能只有稍微差别.与单体相比,簇合物的结构参数变化较大.由ΔG0T可知,298.2 K下单体形成二聚体在热力学上是不利的,而形成三聚体和四聚体是有利的.     

Structure and vibrations of Ga(n)N(n) (n = 3-10) clusters

The structure and harmonic vibrations of Ga(n)N(n) (n = 3-10) clusters have been investigated using the B3LYP (Becke 3-parameter-Lee-Yang-Parr) density functional theory. All structures are found to be cumulenic D(nh) rings (equal bonds, alternating angles), with one intense out of plane mode and three infrared-active degenerate modes, of which the highest one is extremely intense and asymptotically increases to 1029 cm(-1) for n = 10. Comparisons with C2n, B(n)N(n), and Al(n)N(n) clusters, the structure and bonding type for the Ga(n)N(n) (n=3-10) clusters are consistent with those of the C2n (n = 3, 5, 7, ...) clusters, the B(n)N(n) (n = 3-10), and Al(n)N(n) (n = 3-9) clusters.     

Carbon chains and the (5,5) single-walled nanotube: structure and energetics versus length

Reliable thermochemistry is computed for infinite stretches of pure-carbon materials including acetylenic and cumulenic carbon chains, graphene sheet, and single-walled carbon nanotubes (SWCNTs) by connection to the properties of finite size molecules that grow into the infinitely long systems. Using ab initio G3 theory, the infinite cumulenic chain (:C[double bond]C[double bond]C[double bond]C:) is found to be 1.9+/-0.4 kcal/mol per carbon less stable in free energy at room temperature than the acetylenic chain (.C[triple bond]C-C[triple bond]C.) which is 24.0 kcal/mol less stable than graphite. The difference between carbon-carbon triple, double, and single bond lengths (1.257, 1.279, and 1.333 A, respectively) in infinite chains is evident but much less than with small hydrocarbon molecules. These results are used to evaluate the efficacy of similar calculations with the less rigorous PM3 semiempirical method on the (5,5) SWCNT, which is too large to be studied with high-level ab initio methods. The equilibrium electronic energy change for C(g)-->C[infinite (5,5) SWCNT] is -166.7 kcal/mol, while the corresponding free energy change at room temperature is -153.3 kcal/mol (6.7 kcal/mol less stable than graphite). A threefold alternation (6.866, 6.866, and 6.823 A) in the ring diameter of the equilibrium structure of infinitely long (5,5) SWCNT is apparent, although the stability of this structure over the constant diameter structure is small compared to the zero point energy of the nanotube. In general, different (n,m) SWCNTs have different infinite tube energetics, as well as very different energetic trends that vary significantly with length, diameter, and capping.     

Density-functional study of structural and electronic properties of Si(n)C(n) (n=1-10) clusters

Density-functional theory with generalized gradient approximation for the exchange-correlation potential has been used to calculate the structural and electronic structure of Si(n)C(n) (n=1-10) clusters. The geometries are found to undergo a structural change from two dimensional to three dimensional when the cluster size n equals 4. Cagelike structures are favored as the cluster size increases. A distinct segregation between the silicon and carbon atoms is observed for these clusters. It is found that the C atoms favor to form five-membered rings as the cluster size n increases. However, the growth motif for Si atoms is not observed. The Si(n)C(n) clusters at n=2, 6, and 9 are found to possess relatively higher stability. On the basis of the lowest-energy geometries obtained, the size dependence of cluster properties such as binding energy, HOMO-LUMO gap, Mulliken charge, vibrational spectrum, and ionization potential has been computed and analyzed. The bonding characteristics of the clusters are discussed.     

Novel pentagonal silicon rings and nanowheels stabilized by flat pentacoordinate carbon(s)

It is predicted by accurate density functional and coupled-cluster theory that planar [Si(5)C](2-) and [Si(5)C](1-) rings can be stabilized by flat pentacoordinate carbon-silicon bonds. The energy difference of the [Si(5)C](2-) dianion from the lowest energy three-dimensional isomer is about 12.2 kcal∕mol at the level of the density functional theory using the Becke 3-parameter (exchange), Lee, Yang and Parr functional, and the triple-ζ doubly polarized basis sets. Stable composite [Si(5)C](2) structures are formed either as nanowheels with axial C-C bonds of 1.51 A? or as isoenergetic pentagonal graphiticlike layers with double C-C distance (3.02 A?) and almost double aromaticity index, based on nucleus independent chemical shifts. Both of these structures are at least 12 kcal∕mol lower in energy than the lowest energy Si(10)C(2) structure reported in the literature, but about 5 kcal∕mol higher than the lowest energy structure found here.     

Full configuration interaction potential energy curves for the X (1)Sigma(g) (+), B (1)Delta(g), and B(') (1)Sigma(g) (+) states of C(2): a challenge for approximate methods

The C(2) molecule exhibits unusual bonding and several low-lying excited electronic states, making the prediction of its potential energy curves a challenging test for quantum chemical methods. We report full configuration interaction results for the X (1)Sigma(g) (+), B (1)Delta(g), and B(') (1)Sigma(g) (+) states of C(2), which exactly solve the electronic Schrodinger equation within the space spanned by a 6-31G( *) basis set. Within the D(2h) subgroup used by most electronic structure programs, these states all have the same symmetry ((1)A(g)), and all three states become energetically close for interatomic distances beyond 1.5 A. The quality of several single-reference ab initio methods is assessed by comparison to the benchmark results. Unfortunately, even coupled-cluster theory through perturbative triples using an unrestricted Hartree-Fock reference exhibits large nonparallelity errors (>20 kcal mol(-1)) for the ground state. The excited states are not accurately modeled by any commonly used single-reference method, nor by configuration interaction including full quadruple substitutions. The present benchmarks will be helpful in assessing theoretical methods designed to break bonds in ground and excited electronic states.     

Isomer stability of N6C6H6 cages

Recent theoretical studies have identified carbon-nitrogen cages that are potentially stable high energy density materials (HEDM). One such molecule is an N(6)C(6)H(6) cage in which a six-membered ring of nitrogen is bonded to C(3)H(3) triangles on both sides. This molecule is based on the structure of the most stable N(12) cage, with six carbon atoms substituted into the structure. In the current study, several N(6)C(6)H(6) isomers (including the previously studied cage) are examined by theoretical calculations to determine which is actually the most stable. Stability will be evaluated from two points of view: (1) thermodynamic stability of one isomer versus another and (2) kinetic stability of each isomer as determined by the energetics of bond breaking. Density functional theory (B3LYP), perturbation theory (MP2 and MP4), and coupled-cluster theory (CCSD(T)) are used in this study, along with the correlation-consistent basis sets of Dunning. Trends in thermodynamic and kinetic stability are discussed.     

Singlet-triplet energy separation of cyclobutylidene

Ab initio (MP2, CCSD(T)) and density functional theory (BLYP, B3LYP) calculations provide insight concerning novel aspects of structure and bonding in cyclobutylidene (1). Singlet cyclobutylidene ((1)1) adopts a bicyclobutane-like structure (C(s) symmetry) that includes a weak, transannular bonding interaction between the carbene carbon and the opposing CH(2) group. Conformational ring inversion in (1)1 occurs through a transition state of C(2)(v)() symmetry (TS(1)1) with an enthalpy barrier of approximately 3 kcal/mol. Stabilization afforded the singlet state by the transannular interaction appears to be largely offset by a loss of hyperconjugative stabilization from the adjacent C-H bonds. Triplet cyclobutylidene ((3)1) exhibits a C(2)(v)() structure and conventional bonding. The triplet state lies 5.9 kcal/mol above the singlet ground state at the CCSD(T)/TZP//CCSD(T)/DZP level of theory. The singlet-triplet energy gap of cyclobutylidene (-5.9 kcal/mol) lies between that of an acyclic analogue, dimethylcarbene (-1.6 kcal/mol), and a highly strained analogue, cyclopropylidene (-13.8 kcal/mol). The magnitude of the energy gap suggests that triplet cyclobutylidene ((3)1) will be thermally accessible under a variety of experimental conditions.     

Structure and stability of B5C and C5B clusters

Geometry optimizations and vibrational frequencies of B5C and C5B clusters were calculated with the Becke-3LYP method using the 6-311+G(d) basis set and some stable configurations of B5C and C5B clusters have been found. The most stable structure of B5C is a planar six-membered ring. However, for C5B clusters, the most stable structure is linear with a boron atom in position 3. Various configurations of B5C clusters containing three-membered boron rings have predominance in energy, whereas various configurations of C5B clusters containing three-membered carbon rings are disadvantageous in energy. In B5C clusters, isomer2 can be converted into isomer1 by surmounting an energy barrier of 43.83 kJ.mol(-1). In C5B clusters, the conversions of isomer5 into isomer2 and isomer7 into isomer2 have energy barriers of 19.66 and 20.57 kJ.mol(-1), respectively.