首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 281 毫秒
1.
Two new less-symmetric P(OH)3 isomers that are more stable than the C3 structure are found at the density functional theory (B3PW91, B3LYP), MP2, and CCSD(T) levels with the large aug-cc-pvdz/pvtz basis sets. The C1 and C3 structures are qualitatively different from those found for the As(OH)3 molecule. An additional lower lying P(OH)3 structure with Cs symmetry has been obtained. With the largest basis set the Cs isomer is predicted to be the most stable. However, the inclusion of zero-point-energy corrections induces an inversion between the Cs and C1 isomers, with the latter becoming the lowest energy structure at the highest correlated level. Increasing inclusion of electronic correlation effects reduces the energy difference between the C1 and Cs structures while the C1-C3 energy difference and C1-Cs interconversion barrier become larger. In all cases, energy differences and barrier heights are around 1 kcal/mol.  相似文献   

2.
Jacobsen H  Donahue JP 《Inorganic chemistry》2008,47(21):10037-10045
The molecular and electronic structure of monomeric ([Fe(S2C2H2)2](z), [Fe(S2C2(C6H4-p-OCH3)2)2](z)) and dimeric ([{Fe(S2C2H2)2}2](z)) iron bis(dithiolene) complexes, and of their phosphine adducts ([(PH3)Fe(S2C2H2)2](z), [(P(C6H5)3)Fe(S2C2H2)2](z), [(PH3)Fe(S2C2(C6H4-p-OCH3)2)2](z)), carrying various charges (z = 0, 1-, 2-), have been investigated by density functional theory (DFT). Net total spin polarization values S of zero, two, and four have been considered for all neutral model compounds and their dianions, whereas all monoanions have been examined with net total spin polarization values S of one, three, and five. The DFT calculations utilized the pure functional BP86, as well as the hybrid functionals B3LYP and B3LYP*. For the monomers, the calculations reveal the presence of redox non-innocent dithiolene ligands and antiferromagnetic coupling between the ligands and the metal center. For the dimers, complexes with antiferromagnetically coupled iron centers have been found to represent structures of low energy, if not lowest energy structures. The spin-coupling constant of [{Fe(S2C2H2)2}2](2-) is calculated as J = -230 cm(-1). On the basis of the computational results, a model for reversible, electrochemically controlled binding and release of phosphine ligands to iron bis(dithiolene) complexes is proposed. Only BP86 and B3LYP* results, but not those of B3LYP calculations, are in qualitative agreement with experimental findings. BP86 calculations provide the best quantitative match in comparison with the experiment.  相似文献   

3.
Extensive semiempirical calculations of the hexaanions of IPR (isolated pentagon rule) and non-IPR isomers of C(68)-C(88) and IPR isomers of C(90)-C(98) followed by DFT calculations of the lowest energy structures were performed to find the carbon cages that can provide the most stable isomers of M(3)N@C(2n) clusterfullerenes (M = Sc, Y) with Y as a model for rare earth ions. DFT calculations of isomers of M(3)N@C(2n) (M = Sc, Y; 2n = 68-98) based on the most stable C(2n)(6-) cages were also performed. The lowest energy isomers found by this methodology for Sc(3)N@C(68), Sc(3)N@C(78), Sc(3)N@C(80), Y(3)N@C(78), Y(3)N@C(80), Y(3)N@C(84), Y(3)N@C(86), and Y(3)N@C(88) are those that have been shown to exist by single-crystal X-ray studies as Sc(3)N@C(2n) (2n = 68, 78, 80), Dy(3)N@C(80), and Tb(3)N@C(2n) (2n = 80, 84, 86, 88) clusterfullerenes. Reassignment of the carbon cage of Sc(2)@C(76) to the non-IPR Cs: 17490 isomer is also proposed. The stability of nitride clusterfullerenes was found to correlate well with the stability of the empty 6-fold charged cages. However, the dimensions of the cage in terms of its ability to encapsulate M(3)N clusters were also found to be an important factor, especially for the medium size cages and the large Y(3)N cluster. In some cases the most stable structures are based on the different cage isomers for Sc(3)N and Y(3)N clusters. Up to the cage size of C(84), non-IPR isomers of C(2n)(6-) and M(3)N@C(2n) were found to compete with or to be even more stable than IPR isomers. However, the number of adjacent pentagon pairs in the most stable non-IPR isomers decreases as cage size increases: the most stable M(3)N@C(2n) isomers have three such pairs for 2n = 68-72, two pairs for n = 74-80, and only one pair for n = 82, 84. For C(86) and C(88) the lowest energy IPR isomers are much more stable than any non-IPR isomer. The trends in the stability of the fullerene isomers and the cluster-cage binding energies are discussed, and general rules for stability of clusterfullerenes are established. Finally, the high yield of M(3)N@C(80) (Ih) clusterfullerenes for any metal is explained by the exceptional stability of the C(80)(6-) (Ih: 31924) cage, rationalized by the optimum distribution of the pentagons leading to the minimization of the steric strain, and structural similarities of C(80) (Ih: 31924) with the lowest energy non-IPR isomers of C(760(6-), C(78)(6-), C(82)(6-), and C(84)(6-) pointed out.  相似文献   

4.
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.  相似文献   

5.
Using the first-principles method with the generalized gradient approximation, the authors have studied the structural and electronic properties of Al(12)X(+) (X=C, Si, Ge, Sn, and Pb) clusters in detail. The ground state of Al(12)C(+) is a low symmetry C(s) structure instead of an icosahedron. However, the Si, Ge, Sn, and Pb atom doped cationic clusters favor icosahedral structures. The ground states for Al(12)Si(+) and Al(12)Ge(+) are icosahedra, while the C(5nu) structures optimized from an icosahedron with a vertex capped by a tetravalent atom have the highest binding energy for Al(12)Sn(+) and Al(12)Pb(+) clusters. The I(h) structure and the C(5nu) structure are almost degenerate for Al(12)Ge(+), whose binding energy difference is only 0.03 eV. The electronic properties are altered much by removing an electron from the neutral cluster. The binding strength of a valence electron is enhanced, while the binding energy of the cluster is reduced much. Due to the open electronic shell, the band gaps between the highest occupied molecular orbital and the lowest unoccupied molecular orbital are approximately 0.3 eV for the studied cationic clusters.  相似文献   

6.
We performed a constrained search, combined with density-functional theory optimization, of low-energy geometric structures of silicon clusters Si(39), Si(40), Si(50), Si(60), Si(70), and Si(80). We used fullerene cages as structural motifs to construct initial configurations of endohedral fullerene structures. For Si(39), we examined six endohedral fullerene structures using all six homolog C(34) fullerene isomers as cage motifs. We found that the Si(39) constructed based on the C(34)(C(s):2) cage motif results in a new leading candidate for the lowest-energy structure whose energy is appreciably lower than that of the previously reported leading candidate obtained based on unbiased searches (combined with tight-binding optimization). The C(34)(C(s):2) cage motif also leads to a new candidate for the lowest-energy structure of Si(40) whose energy is notably lower than that of the previously reported leading candidate with outer cage homolog to the C(34)(C(1):1). Low-lying structures of larger silicon clusters Si(50) and Si(60) are also obtained on the basis of preconstructed endohedral fullerene structures. For Si(50), Si(60), and Si(80), the obtained low-energy structures are all notably lower in energy than the lowest-energy silicon structures obtained based on an unbiased search with the empirical Stillinger-Weber potential of silicon. Additionally, we found that the binding energy per atom (or cohesive energy) increases typically >10 meV with addition of every ten Si atoms. This result may be used as an empirical criterion (or the minimal requirement) to identify low-lying silicon clusters with size larger than Si(50).  相似文献   

7.
The electronic and geometric structures, total and binding energies, first and second energy differences, harmonic frequencies, point symmetries, and highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO–LUMO) gaps of small and neutral Bn (n = 2–12) clusters have been investigated using density functional theory (DFT), B3LYP with 6‐311++G(d,p) basis set. Linear, planar, convex, quasi‐planar, three‐dimensional (3D) cage, and open‐cage structures have been found. None of the lowest energy structures and their isomers has an inner atom; i.e., all the atoms are positioned at the surface. Within this size range, the planar and quasi‐planar (convex) structures have the lowest energies. The first and the second energy differences are used to obtain the most stable sizes. A simple growth path is also discussed with the studied sizes and isomers. The results have been compared with previously available theoretical and experimental works. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007  相似文献   

8.
C(2n)N7(-) cluster ions are produced by laser ablating on the K(3)[Fe(CN)6] sample. DFT calculations have been performed for these cluster anions. Various isomeric structures of these clusters are optimized and their energies are compared to find the most stable isomers. The most stable structure for C8N7(-) is similar to that of adenine by theoretical calculation, which is in agreement with the collision-induced dissociation (CID) experimental results. With the increasing even numbers of C atoms from 8 to 16, the N atoms in the double-ring structure are gradually substituted by C atoms from the six-membered ring to the five-membered ring. All these C(2n)N7(-) (n = 3-9) clusters exhibit planar aromatic characters. The energy difference and incremental binding energy analyses show that C(2n)N7(-) (n = 4-8) clusters are more stable than C6N7(-) and C18N7(-), which are consistent with the observed mass spectrum.  相似文献   

9.
We report ab initio calculations of the structures, binding energies, and total spins of the clusters Ni(13), Ni(19), Ni(23), Ni(26), Ni(12)Fe, Ni(11)Fe(2), Ni(18)Fe, Ni(17)Fe(2), Ni(22)Fe, Ni(20)Fe(3), and Ni(25)Fe using a density-functional method that employs linear combination of atomic orbitals as basis sets, nonlocal norm-conserving pseudopotentials, and the generalized gradient approximation for exchange and correlation. Our results show that the Fe-doped Ni clusters, which have icosahedral or polyicosahedral ground-state structures similar to those of the corresponding pure Ni clusters, are most stable with the Fe atoms occupying internal positions, as has also been inferred from experimental results on the adsorption of molecular nitrogen on the cluster surfaces. We also rule out the possibility that the experimentally observed difference between the (nonpolyicosahedral) configurations of N(2)-saturated Ni(26) and N(2)-saturated Ni(25)Fe be due to the influence of the Fe atom on the energy of the underlying metal cluster.  相似文献   

10.
Density functional theory (DFT) calculations of protonated methanol-water mixed clusters, H (+)(MeOH) 1(H 2O) n ( n = 1-8), were extensively carried out to analyze the hydrogen bond structures of the clusters. Various structural isomers were energy optimized, and their relative energies with zero point energy corrections and temperature dependence of the free energies were examined. Coexistence of different morphological isomers was suggested. Infrared spectra were simulated on the basis of the optimized structures. The infrared spectra were also experimentally measured for n = 3-9 in the OH stretching vibrational region. The observed broad bands in the hydrogen-bonded OH stretch region were assigned in comparison with the simulations. From the DFT calculations, the preferential proton location was also investigated. Clear correlations between the excess proton location and the cluster morphology were found.  相似文献   

11.
The geometric and electronic structures of Si(n), Si(n)-, and PSi(n-1) clusters (2 < or = n < or = 13) have been investigated using the ab initio molecular orbital theory formalism. The hybrid exchange-correlation energy functional (B3LYP) and a standard split-valence basis set with polarization functions (6-31+G(d)) were employed to optimize geometrical configurations. The total energies of the lowest energy isomers thus obtained were recalculated at the MP2/aug-cc-pVTZ level of theory. Unlike positively charged clusters, which showed similar structural behavior as that of neutral clusters [Nigam et al., J. Chem. Phys. 121, 7756 (2004)], significant geometrical changes were observed between Si(n) and Si(n)- clusters for n = 6, 8, 11, and 13. However, the geometries of P substituted silicon clusters show similar growth as that of negatively charged Si(n) clusters with small local distortions. The relative stability as a function of cluster size has been verified based on their binding energies, second difference in energy (Delta2 E), and fragmentation behavior. In general, the average binding energy of Si(n)- clusters is found to be higher than that of Si(n) clusters. For isoelectronic PSi(n-1) clusters, it is found that although for small clusters (n < 4) substitution of P atom improves the binding energy of Si(n) clusters, for larger clusters (n > or = 4) the effect is opposite. The fragmentation behavior of these clusters reveals that while small clusters prefer to evaporate monomer, the larger ones dissociate into two stable clusters of smaller size. The adiabatic electron affinities of Si(n) clusters and vertical detachment energies of Si(n)- clusters were calculated and compared with available experimental results. Finally, a good agreement between experimental and our theoretical results suggests good prediction of the lowest energy isomeric structures for all clusters calculated in the present study.  相似文献   

12.
Density functional theory calculations are performed on small cationic transition metal doped silver clusters, Ag5X+ (X = Sc, Ti, V, Cr, Mn, Fe, Co, and Ni) using the B3LYP and BP86 functionals. Several two-dimensional and three-dimensional isomers with the dopant at a high coordinated site are found to be close in energy. The relative energy of the isomers is checked with CCSD(T). The interaction between the dopant 3d electrons and the host is discussed by considering the density of states and the shape of the molecular orbitals. A large local spin magnetic moment on the dopant atom is predicted.  相似文献   

13.
The geometric structures, isomeric stabilities, and potential energy profiles of various isomers and transition states in Si(3)H(2) neutral, cation and anion are investigated at the coupled-cluster singles, doubles (triples) level of theory. For the geometrical survey, the basis sets used are of the Dunning's correlation consistent basis sets of triple-zeta quality (cc-pVTZ) for the neutral and cation and the Dunning's correlation consistent basis sets of double-zeta quality with diffuse functions (aug-cc-pVDZ) for the anion. For the final energy calculations, the aug-cc-pVTZ: Dunning's correlation consistent basis sets of triple-zeta quality with diffuse functions and cc-pVQZ: Dunning's correlation consistent basis sets of quadruple-zeta quality basis sets are used for the neutral and the aug-cc-pVTZ ones for the cation and anion. The global minimum neutral (I-1: (1)A(1)) has the same framework as that (cyclopropenylidene) of the C(3)H(2) molecule. Other low-lying three isomers (I-2, I-3, and I-4) are also predicted to be within 20 kJ/mol. Five transition states are optimized and their energy relationships with the isomers are clarified. The geometric structure of the global minimum cation (C-1: (2)A(1)) has the same framework as that of the neutral, but that of the anion (A-1: (2)A(')) differs very much from those of the neutral and cation. The calculated vertical and adiabatic ionization potentials from the global minimum neutral (I-1) are 7.85 and 7.77 eV, respectively. The adiabatic electron affinity of the neutral I-1 and the electron detachment energy of the global minimum anion (A-1) are predicted to be 1.21 and 1.92 eV, respectively. The two-electron three-centered bond is widely observed in the present Si(3)H(2) neutral, cation, and anion. The contour plots of their localized molecular orbitals clearly show the existence of such nonclassical chemical bonds.  相似文献   

14.
Equilibrium structures of H(4-n)Ti(OH)n (n = 2-4) molecules and the Ti(OH)4 dimer and trimers were optimized at the B3LYP level of theory. Theoretical vibrational frequencies of TiO stretching modes obtained with several basis sets were compared with the existing experimental frequencies of these vibrations, and the 6-31+G(d) set was chosen for cluster calculations. Only one energy minimum was found for the [Ti(OH)4](2) dimer, but two isomers without symmetry elements stabilized by internal hydrogen bonds and two isomers, belonging to C(s) and C(i) point groups, with free OH groups were found as minima at the [Ti(OH)4](3) potential energy surface. The structure with the linear arrangement of hexacoordinated titanium atoms in the Ti3O12 skeleton may be proposed for trimeric species observed in liquid titanium alkoxides as the only structure satisfying experimental spectroscopic evidence about the presence of center of inversion in these species. Frequency changes of TiO4 modes which accompany the oligomer formation are predicted and discussed.  相似文献   

15.
从团簇角度对TiNi形状记忆合金进行了量子化学从头算研究。设计并优化了等原子比(TiNi)x(x=2~4)簇的多种可能几何结构,并对较稳定构型进行电子结构的分析。结果表明,等原子比的(TiNi)n团簇以TiNi成键为主要分子骨架,小团簇有较多能量接近的异构体,TiTi成键对能量降低有较大贡献。  相似文献   

16.
Important structural isomers of NH(4) (+)(H(2)O)(n=4,6) have been studied by using density functional theory, Moller-Plesset second order perturbation theory, and coupled-cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. The zero-point energy (ZPE) correction to the complete basis set limit of the CCSD(T) binding energies and free energies is necessary to identify the low energy structures for NH(4) (+)(H(2)O)(n=4,6) because otherwise wrong structures could be assigned for the most probable structures. For NH(4) (+)(H(2)O)(6), the cage-type structure, which is more stable than the previously reported open structure before the ZPE correction, turns out to be less stable after the ZPE correction. In first principles Car-Parrinello molecular dynamics simulations around 100 K, the combined power spectrum of three lowest energy isomers of NH(4) (+)(H(2)O)(4) and two lowest energy isomers of NH(4) (+)(H(2)O)(6) explains each experimental IR spectrum.  相似文献   

17.
A combination of experiment and density functional theory was used to investigate the energetics of CO adsorption onto several small M(x)S(y)(+) (M = Mo, W; x/y = 2/6, 3/7, 5/7, 6/8) clusters as a probe of their atomic and electronic structure. Experimentally, tandem mass spectrometry was used to measure the relative yields of M(x)S(y)(+)(CO)(n) cluster adducts formed by collisions between a beam of mass-selected M(x)S(y)(+) cluster ions and CO molecules in a high-pressure collision cell (hexapole ion guide). The most probable M(x)S(y)(+)(CO)(n) adducts observed are those with n < or = x, that is, only one CO molecule bound to each metal site. The notable exception is the M(5)S(7)(+) cluster, for which the n = 6 adduct is found to have nearly the same intensity as the n = x = 5 adduct. Density functional calculations were used to search for the lowest energy structures of the bare M(x)S(y)(+) clusters and to obtain their relative stability for sequential CO binding. The calculated trends in CO binding energies were then compared to the experimental adduct distributions for assigning the ground-state structures. In this way, it was possible to distinguish between two nearly isoenergetic ground-state isomers for the M(2)S(6)(+) and M(3)S(7)(+) clusters, as only one isomer gave a calculated CO stabilization energy trend that was consistent with the experimental data. Similar comparisons of predicted and observed CO adsorption trends also provide evidence for assigning the ground-state structures of the M(5)S(7)(+) and M(6)S(8)(+) clusters. The latter contain metallic cores with most of the sulfur atoms bonded along the edges or in the faces of the metal core structure. The n = 6 and 7 adducts of M(5)S(7)(+) are predicted to be more stable than the n = x = 5 adduct, but only the n = 6 adduct is observed experimentally. The DFT calculations show that the n = 7 adduct undergoes internal bond breaking whereas the n = 6 framework is stable, albeit highly distorted. For the M(6)S(8)(+) cluster, the calculations predict that the two lowest energy isomers can bind more than six CO molecules without fragmentation; however, the apparent binding energy drops significantly for adducts with n > 6. In general, the ability of these small M(x)S(y)(+) clusters to bind more CO molecules than the number of metal atoms is a balance between the gain in CO adsorption energy versus the strain introduced into the cluster structure caused by CO crowding, the consequences of which can be fragmentation of the M(x)S(y)(+)(CO)(n) cluster adduct (n > x).  相似文献   

18.
The molecular structures of low-lying isomers of anionic and neutral sodium auride clusters have been studied computationally at the second-order M?ller-Plesset perturbation theory level using quadruple-ζ basis sets augmented with a double set of polarization functions. The first vertical detachment energies were calculated at the M?ller-Plesset level as the energy difference between the cluster anion and the corresponding neutral cluster. The photodetachment energies of higher-lying ionization channels were calculated by adding electronic excitation energies of the neutral clusters to the first vertical detachment energy. The excitation energies were calculated at the linear response approximate coupled-cluster singles and doubles level using the anionic cluster structures. The obtained ionization energies for NaAu(-), NaAu(2)(-), NaAu(3)(-), NaAu(4)(-), Na(2)Au(2)(-), Na(2)Au(3)(-), Na(3)Au(3)(-), and Na(2)Au(4)(-) were compared to values deduced from experimental photoelectron spectra. Comparison of the calculated photoelectron spectra for a few energetically low-lying isomers shows that the energetically lowest cluster structures obtained in the calculations do not always correspond to the clusters produced experimentally. Spin-component-scaled second-order M?ller-Plesset perturbation theory calculations shift the order of the isomers such that the observed clusters more often correspond to the energetically lowest structure, whereas the spin-component-scaled approach does not improve the photodetachment energies of the sodium aurides. The potential energy surface of the sodium aurides is very soft, with several low-lying isomers requiring an accurate electron correlation treatment. The calculations show that merely the energetic criterion is not a reliable means to identify the structures of the observed sodium auride clusters; other experimental information is needed to ensure a correct assignment of the cluster structures. The cluster structures of nonstoichiometric anionic sodium aurides have been determined by comparing calculated ionization energies for low-lying structures of the anionic clusters with experimental data.  相似文献   

19.
The electronic structures of the cluster compound alpha-Nb(3)Cl(8) and the intercalated phase beta'-NaNb(3)Cl(8) have been studied by core level and valence band X-ray and ultraviolet photoelectron spectroscopy (XPS and UPS), diffuse reflectance spectroscopy, and charge-self-consistent molecular orbital (CSC-EH) and band structure (CSC-EH-TB) calculations. The crystal structures of the two compounds consist of layers of interconnected Nb(3)Cl(13) units. XP and UP valence band spectra as well as the band structure calculations show well separated sets of Cl 3p levels at lower energy (higher binding energy) and Nb 4d levels at higher energy (lower binding energy), indicative of mainly ionic Nb-Cl bonding. The UP spectra of alpha-Nb(3)Cl(8) reveal a triple-peak structure for the Nb 4d levels, corresponding to the 1a(1), 1e, and 2a(1) metal-metal bonding orbitals of a seven-electron Nb(3) cluster as suggested by theory. The valence band shapes are in good agreement with the theoretical density-of-states curve. The relative intensities in the XP valence band spectra of beta'-NaNb(3)Cl(8) evidence an additional electron in the Nb 4d orbitals. The better resolved UP spectra, however, show a broad pattern for the Nb 4d levels which is neither expected from a single cluster model nor given by the theoretical density-of-states curve. Possible origins for this discrepancy between experiment and theory are discussed. Electrostatic interactions between Na(+) and Cl(-) ions in beta'-NaNb(3)Cl(8) lead to a narrowing of the Cl 2p core level and Cl 3p valence band signals, the former being shifted to higher binding energy in comparison to alpha-Nb(3)Cl(8). Two rather narrow absorptions in the optical spectra of alpha-Nb(3)Cl(8) are assigned to the transitions from the 1a(1) and 1e levels into the singly occupied 2a(1) orbital. These absorptions are missing in the optical spectra of beta'-NaNb(3)Cl(8), in agreement with the theoretical expectations for an eight electron Nb(3) cluster.  相似文献   

20.
Density functional theory B3PW91/6-31+G* calculations on BenCm (n=1-10; m=1, 2, ..., to 11-n) clusters have been carried out to examine the effect of cluster size, relative composition, binding energy per atom, HOMO-LUMO gap, vertical ionization potential, and electron affinity on their relative stabilities. The most stable planar cyclic conformations of these clusters always show at least a set of two carbon atoms between two beryllium atoms, while structures where beryllium atoms cluster together, or allow the intercalation of one carbon atom between two of them, generally seem to be the least stable ones. Clusters containing 1, 2, and 3 beryllium atoms (Be2C8, Be3C6, Be2C6, BeC6, Be2C4, BeC4, Be2C2, and BeC2) are identified as clusters of "magic numbers" in terms of their high binding energy per atom, high HOMO-LUMO gap, vertical ionization potential, and second difference in energy per beryllium atom.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号