硅-硫二元团簇[(SiS2)nS]-(n=1～4)的结构和稳定性的量子化学研究

用密度泛函(DFT)方法(B3LYP/6-31+G*)研究了硅硫团簇[(SiS2)nS]-(n=1～4)的可能几何构型,得到各稳定构型的电子结构,并计算了相应的振动频率,预测了稳定构型的振动光谱.由其稳定构型的比较可在理论上预测团簇的生长规律,并可初步预测团簇的形成机理.     

LaC4^n（n=—2,—1,0,＋1,＋2）分子簇的结构与稳定性

本文用密度泛函方法研究了ＬａＣ４ｎ（ｎ＝－２,－１,０,＋１,＋２）分子簇的结构与稳定性。振动频率分析表明,在所提出的九个构型中,当ｎ＝－２,０,＋１,＋２时,稀土位于碳环上最稳定,而当ｎ＝－１时,尽管稀土位于碳环上能量最低,但没有找到稳定的构型,我们的结果还指出,稀土元素是分子簇中对外部环境最敏感的部位,即最具有反应活性     

Photoionization and ab initio study of Ba(H2O)n (n = 1-4) clusters

An experimental and theoretical study of the photoionization energies (IE's) of Ba(H(2)O)(n) clusters containing up to n = 4 water molecules has been performed. The clusters were generated by a pick-up source combining laser vaporization with pulsed supersonic expansion, and then photoionized by radiation of 272.5-340 nm. The experimentally determined IE(e)'s for n = 1 to 4 are 4.56 ± 0.05, 4.26 ± 0.05, 3.90 ± 0.05 and 3.71 ± 0.05 eV. This cluster size dependence of IE is reproduced within ±0.06 eV employing the mPW1PW91 density-functional and CCSD(T, Full) quantum-chemical methods combined with the 6-311++G(d,p) basis set for the H and O atoms and three different relativistic effective core potentials for Ba atoms. The calculations indicate that the lowest energy hydration structures represent the most relevant contributions to both the vertical and adiabatic ionization energies. Experimental and theoretical evidence correlates with the progressive surface-delocalization of the electron from the hydration cavity around the Ba atom and suggests that the intra-cluster electron transfer is possible even for small aggregates.     

Study of pi halogen bonds in complexes C2H(4-n)Fn-ClF (n = 0-2)

The pi-halogen bond may be considered, in a broad sense, essentially a pi-hydrogen bond. Using the counterpoise-corrected potential energy surface method (interaction energy optimization), the stationary structures of the C2H(4-n)Fn-ClF (n = 0-2) complexes with all real frequencies have been obtained at the MP2/aug-cc-pVDZ level. For C2H(4-n)Fn-ClF (n = 0-2), the pi-halogen bond has a long distance and is elongated by the F substituent effect. The pi-halogen bond length order is 2.661 A for C2H4-ClF < 2.745 A for C2H3F-ClF < 2.766 A for g-C2H2F2-ClF < 2.8076 A for trans-C2H2F2-ClF < 2.8079 A for cis-C2H2F2-ClF. For three complexes C2H3F-ClF, g-C2H2F2-ClF, and cis-C2H2F2-ClF, the pi-halogen bonds are further shifted and sloped by the F substituent effect. The F substituent effect reduces also the interaction energy of the pi-halogen bond. The interaction energies are respectively -3.7 for C2H4-ClF, -2.8 for C2H3F-ClF, -2.3 for g-C2H2F2-ClF, -1.9 for cis-C2H2F2-ClF, and -1.8 kcal/mol for trans-C2H2F2-ClF, at the CCSD(T)/aug-cc-pVDZ level. The electron correlation contribution of the interaction energy is large for C2H(4-n)Fn-ClF (n = 0-2), which shows that the stabilities of the pi-halogen bond systems results primarily from the dispersion interaction. In the double F substituent systems, the interaction energy of the pi-halogen bond structure with a longer interaction distance is larger than that of the corresponding pi-hydrogen bond structure with a shorter interaction distance. This may be because there are the large electron correlation contributions of the interaction energy, and a secondary interaction between lone pairs of Cl atom and some atoms (H, C) with positive charges in the pi-halogen bond structure.     

Infrared spectra of O2- x (CO2)n clusters (n=1-6): asymmetric docking at the pi* orbital

Isolated superoxide ions solvated by CO2 have been studied by infrared photodissociation spectroscopy and density-functional theory, using CO2 evaporation upon infrared excitation of the O2- x (CO2)n (n=1-6) parent ions. We can assign the observed frequencies to the asymmetric stretch vibration and its combination bands with the symmetric stretch and the overtone of the bending vibration of CO2 in various binding situations. We interpret our findings with the help of density-functional theory. Our data suggest that only one CO2 moiety binds strongly to the O2-, whereas the rest of the CO2 molecules are weakly bound, which is consistent with the experimental spectra. The lobes of the pi* orbital of O2- provide a template for the structure of the microsolvation environment.     

Theoretical study of mixed silicon-lithium clusters Si(n)Li(p)(+) (n=1-6, p=1-2)

Theoretical study on the structures of neutral and singly charged Si(n)Li(p)((+)) (n=1-6, p=1-2) clusters have been carried out in the framework of the density functional theory (DFT) with the B3LYP functional. The structures of the neutral Si(n)Li(p) and cationic Si(n)Li(p)(+) clusters are found to keep the frame of the corresponding Si(n), Li species being adsorbed at the surface. The localization of the lithium cation is not the same one as that of the neutral atom. The Li(+) ion is preferentially located on a Si atom, while the Li atom is preferentially attached at a bridge site. A clear parallelism between the structures of Si(n)Na(p) and those of Si(n)Li(p) appears. The population analysis show that the electronic structure of Si(n)Li(p) can be described as Si(n)(p)(-)+pLi(+) for the small sizes considered. Vertical and adiabatic ionization potentials, adsorption energies, as well as electric dipole moments and static dipolar polarizabilities, are calculated for each considered isomer of neutral species.     

Structures of [(CO2)n(CH3OH)m](-) (n = 1-4, m = 1, 2) cluster anions

The infrared photodissociation spectra of [(CO 2) n (CH 3OH) m ] (-) ( n = 1-4, m = 1, 2) are measured in the 2700-3700 cm (-1) range. The observed spectra consist of an intense broad band characteristic of hydrogen-bonded OH stretching vibrations at approximately 3300 cm (-1) and congested vibrational bands around 2900 cm (-1). No photofragment signal is observed for [(CO 2) 1,2(CH 3OH) 1] (-) in the spectral range studied. Ab initio calculations are performed at the MP2/6-311++G** level to obtain structural information such as optimized structures, stabilization energies, and vibrational frequencies of [(CO 2) n (CH 3OH) m ] (-). Comparison between the experimental and the theoretical results reveals the structural properties of [(CO 2) n (CH 3OH) m ] (-): (1) the incorporated CH 3OH interacts directly with either CO 2 (-) or C 2O 4 (-) core by forming an O-HO linkage; (2) the introduction of CH 3OH promotes charge localization in the clusters via the hydrogen-bond formation, resulting in the predominance of CO 2 (-).(CH 3OH) m (CO 2) n-1 isomeric forms over C 2O 4 (-).(CH 3OH) m (CO 2) n-2 ; (3) the hydroxyl group of CH 3OH provides an additional solvation cite for neutral CO 2 molecules.     

Characteristics of antiaromatic ring pi multi-hydrogen bonds in (H2O)n-C4H4 (n = 1, 2) complexes

By counterpoise-corrected optimization method, the six antiaromatic ring pi multi-hydrogen bond structures with diversiform shapes for (H2O)n-C4H4 (n = 1,2) have been obtained at the MP2/aug-cc-pVDZ level. At the CCSD(T)/aug-cc-pVDZ level, the interaction energy obtained mainly depends on the numbers of H2O and fold numbers of the pi multi-hydrogen bond. The interaction energy order is -2.342 (1a with pi mono-hydrogen) < -2.777 (1b with pi bi-hydrogen) < -4.683 (2a with pi bi-hydrogen) < -4.734 (2b with pi tri-hydrogen) < -4.782 (2c with pi tri-hydrogen) < -5.009 kcal/mol (2d with pi tetra-hydrogen bond). Strangely, why is the interaction energy of the pi bi-hydrogen bond in 1b close to that of the pi mono-hydrogen bond in 1a (their difference is only 15.7%)? The reason is that a pi-type H-bond (as an accompanying interaction) between two lone pairs of the O-atom and a near pair of H-atoms of C4H4 exists shoulder by shoulder in structures 1a, 2a, 2b, and 2c and contributes to the interaction energy. Another accompanying interaction, a repulsive interaction between the pi H-bond (using the H-atom(s) of H2O) and the near pair of H-atoms of C4H4, is also found. For the structures and interaction energies, the pi-type H-bond produces four effects: bending the strong pi H-bond, attracting the pair of H-atoms of C4H4 so that they deviate from the C4 ring plane, showing the interaction energy contribution, and bringing the larger electron correlation contribution. The repulsive interaction also produces four effects: pushing the pair of H-atoms of C4H4 so that they deviate from its ring plane, elongating the distance of the pi H-bond, promoting the formation of pi-type H-bond, and slightly influencing the interaction energy. In the present paper, one C=C bond with two H2O (over and below the ring plane) forms a pi H-bond link in two ways: a strong-weak pi H-bond link and a strong-strong pi H-bond link. The stability contribution of the former is more favorable than the latter. One H2O forms a pi H-bond with C4H4 in two ways. One strong pi H-bond part (over or below the ring plane) always is accompanied by another H-bond part. The accompanying part is either a weak pi H-bond or pi-type H-bond.     

Structures and energetics of Be(n)Si(n) and Be(2n)Si(n) (n = 1-4) clusters

The structures and energies of Be(n)Si(n) and Be(2n)Si(n) (n = 1-4) clusters have been examined in ab initio theoretical electronic structure calculations. Cluster geometries have been established in B3LYP/6-31G(2df) calculations and accurate relative energies determined by the G3XMP2 method. The two atoms readily bond to each other and to other atoms of their own kind. The result is a great variety of low-energy clusters in a variety of structural types.     

Gitonic and distonic alkanonium dications (diprotonated alkane dications C(n)H(2n+4(2+)), n = 1-4)(1)

The structures and stabilities of gitonic and distonic alkanonium dications, i.e., diprotonated alkane dications C(n)H(2n+4)(2+) (n = 1-4), were investigated at the MP4(SDTQ)/6-311G**//MP2/6-31G** level. The global minimum energy structures (2, 4, 7, and 10) of the C(n)H(2n+4)(2+) dications are double C--H protonated alkanes to give structures with two two electron three-center (2e-3c) bonds. Two different dissociation pathways for the dications, viz deprotonation and demethylation, were also computed. Demethylation was found to be the favorable mode of dissociation.     

Structure and properties of the aluminum borates Al(BO2)n and Al(BO2)n(-), (n = 1-4)

The geometrical and electronic structures of Al(BO(2))(n) and Al(BO(2))(n)(-) (n = 1-4) clusters are computed at different levels of theory including density functional theory (DFT), hybrid DFT, double-hybrid DFT, and second-order perturbation theory. All aluminum borates are found to be quite stable toward the BO(2) and BO(2)(-) loss in the neutral and anion series, respectively. Al(BO(2))(4) belongs to the class of hyperhalogens composed of smaller superhalogens, and should possess a large adiabatic electron affinity (EA(ad)) larger than that of its superhalogen building block BO(2). Indeed, the aluminum tetraborate possesses the EA(ad) of 5.6 eV, which, however, is smaller than the EA(ad) of 7.8 eV of the AlF(4) supehalogen despite BO(2) is more electronegative than F. The EA(ad) decrease in Al(BO(2))(4) is due to the higher thermodynamic stability of Al(BO(2))(4) compared to that of AlF(4). Because of its high EA and thermodynamic stability, Al(BO(2))(4) should be capable of forming salts with electropositive counter ions. We optimized KAl(BO(2))(4) as corresponding to a unit cell of a hypothetical KAl(BO(2))(4) salt and found that specific energy and energy density of such a salt are competitive with those of trinitrotoluol (TNT).     

Infrared spectra of the water-nitrogen complexes (H2O)2-(N2)n(n = 1-4) in argon matrices

The infrared spectra of the water-nitrogen complexes trapped in argon matrices have been studied with Fourier transform infrared absorption spectroscopy. The absorption lines of the H20-N2 1:1, 1:2, 1:n, and 2:1 complexes have been confirmed on the basis of the concentration effects. In addition, we have observed a few lines and propose the assignments for the 2:2, 2:3, and 2:4 complexes in the nu1 symmetric stretching and nu2 bending regions of the proton-acceptor molecule, and in the bonded OH stretching region of the proton-donor molecule. The redshifts in the bonded OH stretching mode and blueshifts in the OH bending mode suggest that the hydrogen bonds in the (H2O)2-(N2)n complexes with n = 1-4 are strengthened by the cooperative effects compared to the pure H2O dimer. Two absorption bands due to the 3:n complexes are also observed near the bonded OH stretching region of the H2O trimer.     

Microsolvation of N2H): the nature of interactions in N2H+-(H2)n (n = 1-14) complexes

Experimental studies of the consecutive growth of N2H + (H2)n clusters led to the discovery of an unusual bonding pattern for species with n = 2-4. Theoretical studies revealed that the ligands are located within five well-separated solvation shells that are visible in structures, values of successive enthalpies and entropies of clustering reactions, vibrational motions, the distribution of atomic charges, and interaction energy decomposition components. The pattern of consecutive enthalpy changes for the second shell (n = 2-5) is complicated. This pattern shows anomalous behavior, although its interpretation is not univocal. A large part of consecutive enthalpies for the clustering reactions is a contribution due to the rotational and vibrational properties of clusters which are difficult for adequate modeling in large systems. The structures of clusters are rationalized based on interaction energy contributions of a different nature. Geometries of complexes are determined by prevailing covalent forces.     

Threshold photoionization and density functional theory studies of the niobium carbide clusters Nb3C(n) (n = 1-4) and Nb4C(n) (n = 1-6)

We have used photoionization efficiency spectroscopy to determine ionization potentials (IP) of the niobium-carbide clusters, Nb3C(n) (n = 1-4) and Nb4C(n) (n = 1-6). The Nb3C2 and Nb4C4 clusters exhibit the lowest IPs for the two series, respectively. For clusters containing up to four carbon atoms, excellent agreement is found with relative IPs calculated using density functional theory. The lowest energy isomers are mostly consistent with the development of a 2 x 2 x 2 face-centered cubic structure of Nb4C4. However, for Nb3C4 a low-lying isomer containing a molecular C2 unit is assigned to the experimental IP rather than the depleted 2 x 2 x 2 nanocrystal isomer. For Nb4C5 and Nb4C6, interpretation is less straightforward, but results indicate isomers containing molecular C2 units are the lowest in energy, suggesting that carbon-carbon bonding is preferred when the number of carbon atoms exceeds the number of metal atoms. A double IP onset is observed for Nb4C3, which is attributed to ionization from the both the lowest energy singlet state and a meta-stable triplet state. This work further supports the notion that IPs can be used as a reliable validation for the geometries of metal-carbide clusters calculated by theory.     

In search of CS2(H2O)(n=1-4) clusters

Gaussian-3 and MP2/aug-cc-pVnZ methods have been used to calculate geometries and thermochemistry of CS(2)(H2O)n, where n=1-4. An extensive molecular dynamics search followed by optimization using these two methods located two dimers, six trimers, six tetramers, and two pentamers. The MP2/aug-cc-pVDZ structure matched best with the experimental result for the CS(2)(H2O) dimer, showing that diffuse functions are necessary to model the interactions found in this complex. For larger CS(2)(H2O)n clusters, the MP2/aug-cc-pVDZ minima are significantly different from the MP2(full)6-31G* structures, revealing that the G3 model chemistry is not suitable for investigation of sulfur containing van der Waals complexes. Based on the MP2/aug-cc-pVTZ free energies, the concentration of saturated water in the atmosphere and the average amount of CS(2) in the atmosphere, the concentrations of these clusters are predicted to be on the order of 10(5) CS(2)(H2O) clusters.cm(-3) and 10(2) CS(2)(H2O)(2) clusters.cm(-3) at 298.15 K. The MP2/aug-cc-pVDZ scaled harmonic and anharmonic frequencies of the most abundant dimer cluster at 298 K are presented, along with the MP2/aug-cc-pVDZ scaled harmonic frequencies for the CS(2)(H(2)O)(n) structures predicted to be present in a low-temperature molecular beam experiment.     

Theoretical study of Al(n) and Al(n)O (n = 2-10) clusters

The stable structures, energies, and electronic properties of neutral, cationic, and anionic clusters of Al(n) (n = 2-10) are studied systematically at the B3LYP/6-311G(2d) level. We find that our optimized structures of Al5(+), Al9(+), Al9(-), Al10, Al10(+), and Al10(-) clusters are more stable than the corresponding ones proposed in previous literature reports. For the studied neutral aluminum clusters, our results show that the stability has an odd/even alternation phenomenon. We also find that the Al3, Al7, Al7(+), and Al7(-) structures are more stable than their neighbors according to their binding energies. For Al7(+) with a special stability, the nucleus-independent chemical shifts and resonance energies are calculated to evaluate its aromaticity. In addition, we present results on hardness, ionization potential, and electron detachment energy. On the basis of the stable structures of the neutral Al(n) (n = 2-10) clusters, the Al(n)O (n = 2-10) clusters are further investigated at the B3LYP/6-311G(2d), and the lowest-energy structures are searched. The structures show that oxygen tends to either be absorbed at the surface of the aluminum clusters or be inserted between Al atoms to form an Al(n-1)OAl motif, of which the Al(n-1) part retains the stable structure of pure aluminum clusters.     

Vibrational spectroscopy and structures of Ni+(C2H2)(n) (n =1-4) complexes

Nickel cation-acetylene complexes of the form Ni(+)(C(2)H(2))(n), Ni(+)(C(2)H(2))Ne, and Ni(+)(C(2)H(2))(n)Ar(m) (n = 1-4) are produced in a molecular beam by pulsed laser vaporization. These ions are size-selected and studied in a time-of-flight mass spectrometer by infrared laser photodissociation spectroscopy in the C-H stretch region. The fragmentation patterns indicate that the coordination number is 4 for this system. The n = 1-4 complexes with and without rare gas atoms are also investigated with density functional theory. The combined IR spectra and theory show that pi-complexes are formed for the n = 1-4 species, causing the C-H stretches in the acetylene ligands to shift to lower frequencies. Theory reveals that there are low-lying excited states nearly degenerate with the ground state for all the Ni(+)(C(2)H(2))(n) complexes. Although isomeric structures are identified for rare gas atom binding at different sites, the attachment of rare gas atoms results in only minor perturbations on the structures and spectra for all complexes. Experiment and theory agree that multiple acetylene binding takes place to form low-symmetry structures, presumably due to Jahn-Teller distortion and/or ligand steric effects. The fully coordinated Ni(+)(C(2)H(2))(4) complex has a near-tetrahedral structure.     

Global optimization of ionic Mg(n)F(2n) (n=1-30) clusters

The global optimization basin-hopping (BH) method has been used to locate the global minima (GM) of Mg(n)F(2n) (n=1-30) clusters using a Born-Mayer-type potential. Some of the GM were particularly difficult to find, requiring more than 1.5 x 10(4) BH steps. We have found that both the binding energy per MgF2 unit and the effective volume of the GM isomers increase almost linearly with n, and that cluster symmetry decreases with cluster size. The data derived from the BH runs reveal a growing density of local minima just above the GM as n increases. Despite this, the attraction basin around each GM is relatively large, since after all their atomic coordinates are randomly displaced by values as high as 2.0 bohrs, the perturbed structures, upon reoptimization, relax back to the GM in more than 50% of the cases (except for n=10 and 11). The relative stabilities derived from energy second differences suggest that n=8,10,13,15, and 20 are probably the magic numbers for these systems. Mass spectrum experiments would be very useful to clarify this issue.     

Structures and vibrational spectra of the sulfur-rich oxides SnO (n = 4-9): the importance of pi*-pi* interactions

The structures of a large number of isomers of the sulfur oxides S(n)O with n = 4-9 have been calculated at the G3X(MP2) level of theory. In most cases, homocyclic molecules with exocyclic oxygen atoms in an axial position are the global minimum structures. Perfect agreement is obtained with experimentally determined structures of S(7)O and S(8)O. The most stable S(4)O isomer as well as some less stable isomers of S(5)O and S(6)O are characterized by a strong pi*-pi* interaction between S==O and S==S groups, which results in relatively long S--S bonds with internuclear distances of 244-262 pm. Heterocyclic isomers are less stable than the global minimum structures, and this energy difference approximately increases with the ring size: 17 (S(4)O), 40 (S(5)O), 32 (S(6)O), 28 (S(7)O), 45 (S(8)O), and 54 kJ mol(-1) (S(9)O). Owing to a favorable pi*-pi* interaction, preference for an axial (or endo) conformation is calculated for the global energy minima of S(7)O, S(8)O, and S(9)O. Vapor-phase decomposition of S(n)O molecules to SO(2) and S(8) is strongly exothermic, whereas the formation of S(2)O and S(8) is exothermic if n<7, but slightly endothermic for S(7)O, S(8)O, and S(9)O. The calculated vibrational spectra of the most stable isomers of S(6)O, S(7)O, and S(8)O are in excellent agreement with the observed data.     

Molecular structures and energetics of the (TiO2)n (n = 1-4) clusters and their anions

The (TiO2)n clusters and their anions for n = 1-4 have been studied with coupled cluster theory [CCSD(T)] and density functional theory (DFT). For n > 1, numerous conformations are located for both the neutral and anionic clusters, and their relative energies are calculated at both the DFT and CCSD(T) levels. The CCSD(T) energies are extrapolated to the complete basis set limit for the monomer and dimer and calculated up to the triple-zeta level for the trimer and tetramer. The adiabatic and vertical electron detachment energies of the anionic clusters to the ground and first excited states of the neutral clusters are calculated at both levels and compared with the experimental results. The comparison allows for the definitive assignment of the ground-state structures of the anionic clusters. Anions of the dimer and tetramer are found to have very closely lying conformations within 2 kcal/mol at the CCSD(T) level, whereas that of the trimer does not. In addition, accurate clustering energies and heats of formation are calculated for the neutral clusters and compared with the available experimental data. Estimates of the titanium-oxygen bond energies show that they are stronger than the group VIB transition metal-oxygen bonds except for tungsten. The atomization energies of these clusters display much stronger basis set dependence than the clustering energies. This allows the calculation of more accurate heats of formation for larger clusters on the basis of calculated clustering energies.