首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
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.  相似文献   

2.
The application of the ab initio genetic algorithm with an embedded gradient has been carried out for the elucidation of global minimum structures of a series of anionic sodium chloride clusters, Na(x)Cl(x+1) (-) (x=1-4), produced in the gas phase using electrospray ionization and studied by photoelectron spectroscopy. These are all superhalogen species with extremely high electron binding energies. The vertical electron detachment energies for Na(x)Cl(x+1) (-) were measured to be 5.6, 6.46, 6.3, and 7.0 eV, for x=1-4, respectively. Our ab initio gradient embedded genetic algorithm program detected the linear global minima for NaCl(2) (-) and Na(2)Cl(3) (-) and three-dimensional structures for the larger species. Na(3)Cl(4) (-) was found to have C(3v) symmetry, which can be viewed as a Na(4)Cl(4) cube missing a corner Na(+) cation, whereas Na(4)Cl(5) (-) was found to have C(4v) symmetry, close to a 3x3 planar structure. Excellent agreement between the theoretically calculated and the experimental spectra was observed, confirming the obtained structures and demonstrating the power of the developed genetic algorithm technique.  相似文献   

3.
The structures of sodium doped water trimers are characterized on the basis of their infrared action spectra in the OH-stretching region and a global optimization approach to identify the lowest energy minima. The most stable structure is an open ring with two contacts of terminal water molecules to the Na atom. This structure explains the dominating feature in the IR depletion spectrum around 3410 cm(-1). Three additional isomer classes were found in an energy window of 12 kJ mol(-1) with vertical ionization energies ranging from ~3.83 eV to ~4.36 eV. These structures show different hydrogen bonding and sodium coordination patterns and are identified by specific spectral features in the IR spectra. The significant abundance of closed rings with an external Na atom, resembling the undoped water trimer, suggests that for larger clusters the picture of the sodium atom being situated on the cluster surface seems adequate.  相似文献   

4.
Coinage and alkali metal mixed clusters, M4Na- (M = Cu, Au) have been investigated experimentally using photoelectron spectroscopy and computationally at correlated ab initio levels. The related Cu4Li-, Ag4Li-, Ag4Na-, and Au4Li- clusters as well as the neutral Cu4Li2 and Cu4Na2 clusters have also been studied computationally. The calculations show that the two lowest isomers of the negatively charged clusters include a pyramidal C4v structure and a planar C2v species. For Cu4Li- and Cu4Na-, the C4v structure is calculated at correlated ab initio level to be 30.9 and 16.9 kJ/mol below the planar C2v isomer, whereas the planar isomers of Au4Li- and Au4Na- are found to be 29.7 and 49.4 kJ/mol below the pyramidal ones. For Ag4Li- and Ag4Na-, the pyramidal isomers are the lowest ones. Comparison of the calculated and measured photoelectron spectra of Cu4Na- and Au4Na- shows that the pyramidal Cu4Na- cluster of C4v symmetry and the planar Au4Na- of C2v symmetry are detected experimentally. Calculations of the magnetically induced current density in Cu4Li- and Cu4Li2 using the Gauge-Including Magnetically Induced Current (GIMIC) method show that strong ring currents are sustained mainly by the highest-occupied molecular orbital primarily derived from the Cu 4s. The GIMIC calculations thus show that the Cu4(2-) ring is -aromatic and that the d orbitals do not play any significant role for the electron delocalization effects. The present study does not support the notion that the square-planar Cu4(2-) is the first example of d-orbital aromatic molecules.  相似文献   

5.
Developments in the preparation of new materials for microelectronics are focusing new attention on molecular systems incorporating several arsenic atoms. A systematic investigation of the As2Fn/As2Fn- systems was carried out using Density Functional Theory methods and a DZP++ quality basis set. Global and low-lying local geometric minima and relative energies are discussed and compared. The three types of neutral-anion separations reported in this work are: the adiabatic electron affinity (EAad), the vertical electron affinity (EAvert), and the vertical detachment energy (VDE). Harmonic vibrational frequencies pertaining to the global minimum for each compound are reported. From the first four studied species (As2Fn, n=1-4), all neutral molecules and their anions are shown to be stable with respect to As-As bond breaking. The neutral As2F molecule and its anion are predicted to have Cs symmetry. We find the trans F-As-As-F isomer of C2h symmetry and a pyramidalized vinylidene-like As-As-F2- isomer of Cs symmetry to be the global minima for the As2F2 and As2F2- species, respectively. The lowest lying minima of As2F3 and As2F3- are vinyl radical-like structures F-As-As-F2 of Cs symmetry. The neutral As2F4 global minimum is a trans-bent (like Si2H4) F2-As-As-F2 isomer of C2 symmetry, while its anion is predicted to have an unusual fluorine-bridged (C(1)) structure. The global minima of the neutral As2Fn species, n=5-8, are weakly bound complexes, held together by dipole-dipole interactions. All such structures have the AsFm-AsFn form, where (m,n) is (2,3) for As2F5, (3,3) for As2F6, (4,3) for As2F7), and (5,3) for As2F8. For As2F8 the beautiful pentavalent F4As-AsF4 structure (analogous to the stable AsF5 molecule) lies about 30 kcal/mol above the AsF3 . . . AsF5 complex. The stability of AsF(5) depends crucially on the strong As-F bonds, and replacing one of these with an As-As bond (in F4As-AsF4) has a very negative impact on the molecule's stability. The anions As2Fn-, n=5-8, are shown to be stable with respect to the As-As bond breaking, and we predict that all of them have fluorine-bridged or fluorine-linked structures. The zero-point vibrational energy corrected adiabatic electron affinities are predicted to be 2.28 eV (As2F), 1.95 eV (As2F2), 2.39 eV (As2F3), 1.71 eV (As2F4), 2.72 eV (As2F5), 1.79 eV (As2F6), 5.26 eV (As2F7), and 3.40 eV (As2F8) from the BHLYP method. Vertical detachment energies are rather large, especially for species with fluorine-bridged global minima, having values up to 6.45 eV (As2F7, BHLYP).  相似文献   

6.
The geometrical and electronic properties of the anionic and neutral V2O6 clusters were studied with the spin unrestricted hybrid density functional B3LYP method. The calculated ground states of both clusters are different from the previous theoretical results. The ground state of V2O6- is found to be a doublet with C2v symmetry, while a doublet with D2h symmetry was previously obtained by Vyboishchikov and Sauer. For neutral V2O6, the ground state is an open-shell singlet with D2h symmetry whose energy is very close to that of the triplet state. In contrast, a closed-shell singlet with D2h symmetry was obtained by Vyboishchikov and Sauer, and Calatayud et al. found a triplet ground state with Cs symmetry. Moreover,the calculated adiabatic and vertical detachment energies of the anion cluster are in much better agreement with the experimental results of photoelectron spectroscopy than previous theoretical values.  相似文献   

7.
叠氮二氢硼多聚体结构和性质的理论研究(英文)   总被引:1,自引:0,他引:1  
本文采用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下单体形成二聚体在热力学上是不利的,而形成三聚体和四聚体是有利的.  相似文献   

8.
The equilibrium structure, stability, and electronic properties of the Al(13)X (X=H,Au,Li,Na,K,Rb,Cs) clusters have been studied using a combination of photoelectron spectroscopy experiment and density functional theory. All these clusters constitute 40 electron systems with 39 electrons contributed by the 13 Al atoms and 1 electron contributed by each of the X (X=H,Au,Li,Na,K,Rb,Cs) atom. A systematic study allows us to investigate whether all electrons contributed by the X atoms are alike and whether the structure, stability, and properties of all the magic clusters are similar. Furthermore, quantitative agreement between the calculated and the measured electron affinities and vertical detachment energies enable us to identify the ground state geometries of these clusters both in neutral and anionic configurations.  相似文献   

9.
Au(2)H(-) was recognized and confirmed as a minor contamination to typical photoelectron spectra of Au(2) (-), produced by laser vaporization of a pure Au target using an ultrahigh purity helium carrier gas. The hydrogen source was shown to be from trace H impurities present in the bulk gold target. Carefully designed experiments using H(2)- and D(2)-seeded helium carrier gas were used to study the electronic structure of Au(2)H(-) and Au(2)D(-) using photoelectron spectroscopy and density functional calculations. Well-resolved photoelectron spectra with vibrational resolution were obtained for Au(2)H(-) and Au(2)D(-). Two isomers were observed both experimentally and theoretically. The ground state of Au(2)H(-) turned out to be linear with a terminal H atom [Au-Au-H](-) ((1)A(1),C(infinity v)), whereas a linear [Au-H-Au](-) ((1)A(1),D(infinity h)) structure with a bridging H atom was found to be a minor isomer 0.6 eV higher in energy. Calculated electron detachment energies for both isomers agree well with the experimental spectra, confirming their existence in the cluster beam. The observation and confirmation of H impurity in pure gold clusters and the 3.44 A Au-Au distance in the [Au-H-Au](-) isomer presented in the current work provide indirect experimental evidence that the anomalous 3.6 A Au-Au distances observed in gold nanowires is due to an "invisible" hydrogen impurity atom.  相似文献   

10.
The structures of neutral boron clusters, B(11), B(16), and B(17), have been investigated using vibrational spectroscopy and ab initio calculations. Infrared absorption spectra in the wavelength range of 650 to 1550 cm(-1) are obtained for the three neutral boron clusters from the enhancement of their near-threshold ionization efficiency at a fixed UV wavelength of 157 nm (7.87 eV) after resonant absorption of the tunable infrared photons. All three clusters, B(11), B(16), and B(17), are found to possess planar or quasi-planar structures, similar to their corresponding anionic counterparts (B(n) (-)), whose global minima were found previously to be planar, using photoelectron spectroscopy and theoretical calculations. Only minor structural changes are observed between the neutral and the anionic species for these three boron clusters.  相似文献   

11.
The group 14 clusters encapsulated by coinage metals in neutral and anionic states X(10)M(0/-) (X = Ge, Sn, Pb and M = Cu, Ag, Au) are investigated using quantum chemical calculations with the DFT/B3LYP functional and coupled-cluster CCSD(T) theory. Addition of transition metals into the empty cages forms high symmetry endohedral structures, except for Ge(10)Ag(0/-). In agreement with experiments available for X(10)Cu, the D(4d) global minima of the anions are calculated to be magic clusters with large frontier orbital gaps, high vertical and adiabatic detachment energies, and large embedding energies and binding energies as compared to those of the empty cages X(10)(2-). The enhanced stability of these magic clusters can be rationalized by the three-dimensional aromaticity.  相似文献   

12.
A systematic and unbiased structure search based on a genetic algorithm in combination with density functional theory (DFT) procedures has been carried out to locate low-energy isomers of Ga(n) up to n = 25. For the smaller clusters up to n = 8 results are checked by coupled cluster singles and doubles with perturbative triples corrections (CCSD(T)) employing a quadruple zeta type basis set. The CCSD(T) calculations confirm a (3)Π(u) ground state for the dimer. Ga(3) has a doublet ground state 0.2 eV below two quartet states, whereas two isoenergetic triplet states are predicted for Ga(4) with D(4h) and a rhombus structure (D(2h)). Three low-lying isomers with doublet electronic states are found for Ga(5): a W-structure (C(2v)), a planar envelope (C(s)) at 0.015 eV, and a non-planar envelope (C(1)) 0.086 eV above the ground state. A triplet state for a trigonal prism (D(3h)) and a singlet for an open prism (C(2v)) are computed with virtually identical energy for Ga(6). The global minimum for Ga(7) is a capped trigonal prism (C(s)) and that for Ga(8) a distorted cube in D(2h). DFT provides a fair agreement with CCSD(T), deviations in dissociation energies are up to 0.2 eV for n ≤ 8. The structures for Ga(n) are mostly irregular for n ≥ 9, those for Ga(12) to Ga(17) can be derived from the truncated decahedron with D(5h) symmetry though highly distorted by Jahn-Teller effects, for example. For Ga(18) to Ga(23) we find stacks of five- and six-membered rings as global minima, e.g., 5-1-5-1-6 for Ga(18). Ga(24) and Ga(25) consist of layers with packing sequence ABCBA similar to those found for clusters of aluminum. The most important feature of computed cohesive energies is a rapid increase with n: for Ga(25) it reaches 2.46 eV, the experimental bulk value is 2.84 eV. Particularly stable clusters for Ga(n) are seen for n = 7, 14, and 20.  相似文献   

13.
This work presents a study of reactions between neutral and negatively charged Au(n) clusters (n=2,3) and molecular hydrogen. The binding energies of the first and second hydrogen molecule to the gold clusters were determined using density functional theory (DFT), second order perturbation theory (MP2) and coupled cluster (CCSD(T)) methods. It is found that molecular hydrogen easily binds to neutral Au(2) and Au(3) clusters with binding energies of 0.55 eV and 0.71 eV, respectively. The barriers to H(2) dissociation on these clusters with respect to Au(n)H(2) complexes are 1.10 eV and 0.59 eV for n=2 and 3. Although negatively charged Au(n) (-) clusters do not bind molecular hydrogen, H(2) dissociation can occur with energy barriers of 0.93 eV for Au(2) (-) and 1.39 eV for Au(3) (-). The energies of the Au(2)H(2) (-) and Au(3)H(2) (-) complexes with dissociated hydrogen molecules are lower than the energies of Au(2) (-)+H(2) and Au(3) (-)+H(2) by 0.49 eV and 0.96 eV, respectively. There is satisfactory agreement between the DFT and CCSD(T) results for binding energies, but the agreement is not as good for barrier heights.  相似文献   

14.
We perform a systematic investigation on small silicon-doped boron clusters B(n)Si (n=1-7) in both neutral and anionic states using density functional (DFT) and coupled-cluster (CCSD(T)) theories. The global minima of these B(n)Si(0/-) clusters are characterized together with their growth mechanisms. The planar structures are dominant for small B(n)Si clusters with n≤5. The B(6)Si molecule represents a geometrical transition with a quasi-planar geometry, and the first 3D global minimum is found for the B(7)Si cluster. The small neutral B(n)Si clusters can be formed by substituting the single boron atom of B(n+1) by silicon. The Si atom prefers the external position of the skeleton and tends to form bonds with its two neighboring B atoms. The larger B(7)Si cluster is constructed by doping Si-atoms on the symmetry axis of the B(n) host, which leads to the bonding of the silicon to the ring boron atoms through a number of hyper-coordination. Calculations of the thermochemical properties of B(n)Si(0/-) clusters, such as binding energies (BE), heats of formation at 0 K (ΔH(f)(0)) and 298 K (ΔH(f)([298])), adiabatic (ADE) and vertical (VDE) detachment energies, and dissociation energies (D(e)), are performed using the high accuracy G4 and complete basis-set extrapolation (CCSD(T)/CBS) approaches. The differences of heats of formation (at 0 K) between the G4 and CBS approaches for the B(n)Si clusters vary in the range of 0.0-4.6 kcal mol(-1). The largest difference between two approaches for ADE values is 0.15 eV. Our theoretical predictions also indicate that the species B(2)Si, B(4)Si, B(3)Si(-) and B(7)Si(-) are systems with enhanced stability, exhibiting each a double (σ and π) aromaticity. B(5)Si(-) and B(6)Si are doubly antiaromatic (σ and π) with lower stability.  相似文献   

15.
The electronic structure and chemical bonding in the Ta 3 (-) cluster are investigated using photoelectron spectroscopy and density functional theory calculations. Photoelectron spectra are obtained for Ta 3 (-) at four photon energies: 532, 355, 266, and 193 nm. While congested spectra are observed at high electron binding energies, several low-lying electronic transitions are well resolved and compared with the theoretical calculations. The electron affinity of Ta 3 is determined to be 1.35 +/- 0.03 eV. Extensive density functional calculations are performed at the B3LYP/Stuttgart +2f1g level to locate the ground-state and low-lying isomers for Ta 3 and Ta 3 (-). The ground-state for the Ta 3 (-) anion is shown to be a quintet ( (5)A 1') with D 3 h symmetry, whereas two nearly isoenergetic states, C 2 v ( (4)A 1) and D 3 h ( (6)A 1'), are found to compete for the ground-state for neutral Ta 3. A detailed molecular orbital analysis is performed to elucidate the chemical boding in Ta 3 (-), which is found to possess multiple d-orbital aromaticity, commensurate with its highly symmetric D 3 h structure.  相似文献   

16.
Small boron clusters have been shown to be planar from a series of combined photoelectron spectroscopy and theoretical studies. However, a number of boron clusters are quasiplanar, such as B(7)(-) and B(12)(-). To elucidate the nature of the nonplanarity in these clusters, we have investigated the electronic structure and chemical bonding of two isoelectronic Al-doped boron clusters, AlB(6)(-) and AlB(11)(-). Vibrationally resolved photoelectron spectra were obtained for AlB(6)(-), resulting in an accurate electron affinity (EA) for AlB(6) of 2.49 ± 0.03 eV. The photoelectron spectra of AlB(11)(-) revealed the presence of two isomers with EAs of 2.16 ± 0.03 and 2.33 ± 0.03 eV, respectively. Global minimum structures of both AlB(6)(-) and AlB(11)(-) were established from unbiased searches and comparison with the experimental data. The global minimum of AlB(6)(-) is nearly planar with a central B atom and an AlB(5) six membered ring, in contrast to that of B(7)(-), which possesses a C(2v) structure with a large distortion from planarity. Two nearly degenerate structures were found for AlB(11)(-) competing for the global minimum, in agreement with the experimental observation. One of these isomers with the lower EA can be viewed as substituting a peripheral B atom by Al in B(12)(-), which has a bowl shape with a B(9) outer ring and an out-of-plane inner B(3) triangle. The second isomer of AlB(11)(-) can be viewed as an Al atom interacting with a B(11)(-) cluster. Both isomers of AlB(11)(-) are perfectly planar. It is shown that Al substitution of a peripheral B atom in B(7)(-) and B(12)(-) induces planarization by slightly expanding the outer ring due to the larger size of Al.  相似文献   

17.
During photoelectron spectroscopy experiments, the spectra of B(11)O(-) and B(10)Au(-) clusters are found to exhibit similar patterns except for a systematic spectral shift of ~0.5 eV, hinting that they possess similar geometric structures. The electron affinities are measured to be 4.02 ± 0.04 eV for B(11)O and 3.55 ± 0.02 eV for B(10)Au. DFT calculations at the B3LYP level show that B(11)O(-) and B(10)Au(-) adopt similar C(1) ((1)A) ground states, which are based on the quasiplanar B(10) cluster interacting with a BO unit and Au, respectively. The B(11)O(-) and B(10)Au(-) clusters are thus valent isoelectronic because both BO and Au can be viewed as monovalent units, forming highly covalent B-BO and B-Au bonds analogous to the B-H bond in B(10)H(-). For B(10)Au(-), we also find a highly symmetric D(10h) ((1)A(1g)) planar molecular wheel as a minimum on the potential energy surface. However, it is 45 kcal/mol above the ground state at the B3LYP level and not viable for experimental observation. Natural bond orbital analyses reveal interesting covalent versus ionic B-Au bonding in the C(1) B(10)Au(-) and D(10h) B(10)Au(-) structures, respectively, providing insight for the design of D(nh) MB(n) molecular wheels.  相似文献   

18.
A home-made magnetic-bottle time-of-flight anion photoelectron spectrometer(PES)for the investigation of binary metal cluster geometry and electron structure is described. The photoelectron spectrometer is installed near the first space focus of home-made reflectron time of flight mass spectrometer(RTOFMS),coupled with laser ablation,pulse supersonic molecular carrier gas cluster source. The magnetic-bottle photoelectron spectrometer's resolution is about 0. 1 eV for 1 eV photoelectrons. The adiabatic electron affinity energies of neutral clusters and some features relative to their excited states can be obtained from the spectra,i. e. ,from the anion's spectra,not only the features of the anion but also the neutral clusters' features can be investigated. The detailed design,construction,and operation of the new apparatus are presented. And studied PbM-(M = Cu,Ag,Au)binary metal cluster anions,the results give clear diagram about their structures and the bond interactions. The adiabatic electron affinity energies obtained by the photoelectron spectrometer agree well with the calculated results using relativistic density functional theory(DFT)method. It show that this anion photoelectron spectrometer can be well used in studying binary metal cluster anions in the experiment condition.  相似文献   

19.
The structures and the electronic properties of two aluminum-doped boron clusters, AlB(7)(-) and AlB(8)(-), were investigated using photoelectron spectroscopy and ab initio calculations. The photoelectron spectra of AlB(7)(-) and AlB(8)(-) are both broad, suggesting significant geometry changes between the ground states of the anions and the neutrals. Unbiased global minimum searches were carried out and the calculated vertical electron detachment energies were used to compare with the experimental data. We found that the Al atom does not simply replace a B atom in the parent B(8)(-) and B(9)(-) planar clusters in AlB(7)(-) and AlB(8)(-). Instead, the global minima of the two doped-clusters are of umbrella shapes, featuring an Al atom interacting ionically with a hexagonal and heptagonal pyramidal B(7) (C(6v)) and B(8) (C(7v)) fragment, respectively. These unique umbrella-type structures are understood on the basis of the special stability of the quasi-planar B(7)(3-) and planar B(8)(2-) molecular wheels derived from double aromaticity.  相似文献   

20.
The electronic structure of Sn(n) (-) clusters (n=4-45) was examined using photoelectron spectroscopy at photon energies of 6.424 eV (193 nm) and 4.661 eV (266 nm) to probe the semiconductor-to-metal transition. Well resolved photoelectron spectra were obtained for small Sn(n) (-) clusters (n< or =25), whereas more congested spectra were observed with increasing cluster size. A distinct energy gap was observed in the photoelectron spectra of Sn(n) (-) clusters with n< or =41, suggesting the semiconductor nature of small neutral tin clusters. For Sn(n) (-) clusters with n> or =42, the photoelectron spectra became continuous and no well-defined energy gap was observed, indicating the onset of metallic behavior for the large Sn(n) clusters. The photoelectron spectra thus revealed a distinct semiconductor-to-metal transition for Sn(n) clusters at n=42. The spectra of small Sn(n) (-) clusters (n< or =13) were also compared with those of the corresponding Si(n) (-) and Ge(n) (-) clusters, and similarities were found between the spectra of Sn(n) (-) and those of Ge(n) (-) in this size range, except for Sn(12) (-), which led to the discovery of stannaspherene (the icosahedral Sn(12) (2-)) previously [L. F. Cui et al., J. Am. Chem. Soc. 128, 8391 (2006)].  相似文献   

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

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