首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
The thermodynamic properties of surface ceria were investigated through equilibrium isotherms determined by flow titration and coulometric titration measurements on high-surface-area ceria and ceria supported on La-modified alumina (LA). While the surface area of pure ceria was found to be unstable under redox conditions, the extent of reduction at 873 K and a P(O2) of 1.6 × 10−26 atm increased with surface area. Because ceria/LA samples were stable, equilibrium isotherms were determined between 873 and 973 K on a 30-wt% ceria sample. Oxidation enthalpies on ceria/LA were found to vary with the extent of reduction, ranging from −500 kJ/mol O2 at low extents of reduction to near the bulk value of −760 kJ/mol O2 at higher extents. To determine whether +3 dopants could affect the oxidation enthalpies for ceria, isotherms were measured for Sm+3-doped ceria (SDC) and Y+3-doped ceria. These dopants were found to remove the phase transition observed in pure ceria below 973 K but appeared to have minimal effect on the oxidation enthalpies. Implications of these results for catalytic applications of ceria are discussed.  相似文献   

2.
The infrared spectra and stability of CO and H2O sorption over Ag-exchanged ZSM-5 zeolite were investigated by using density function theory (DFT). The changes of NBO charge show that the electron transfers from CO molecule to the Ag+ cation to form an σ-bond, and it accompanies by the back donation of d-electrons from Ag+ cation to the CO (π*) orbital as one and two CO molecules are adsorbed on Ag-ZSM-5. The free energy changes ΔG, −5.55 kcal/mol and 6.52 kcal/mol for one and two CO molecules, illustrate that the Ag+(CO)2 complex is unstable at the room temperature. The vibration frequency of C-O stretching of one CO molecule bonded to Ag+ ion at 2211 cm−1 is in good agreement with the experimental results. The calculated C-O symmetric and antisymmetric stretching frequencies in the Ag+(CO)2 complex shift to 2231 cm−1 and 2205 cm−1 when the second CO molecule is adsorbed. The calculated C-O stretching frequency in CO-Ag-ZSM-5-H2O complex shifts to 2199 cm−1, the symmetric and antisymmetric O-H stretching frequencies are 3390 cm−1 and 3869 cm−1, respectively. The Gibbs free energy change (ΔGH2O) is −6.58 kcal/mol as a H2O molecule is adsorbed on CO-Ag-ZSM-5 complex at 298 K. The results show that CO-Ag-ZSM-5-H2O complex is more stable at room temperature.  相似文献   

3.
Infrared reflection absorption spectroscopy that used buried metal layer substrates (BML-IRRAS) and density functional cluster calculations were employed to investigate the water related oxidation reactions of 2H + H2O/Si(1 0 0)-(2 × 1), 2D + H2O/Si(1 0 0)-(2 × 1), and H2O + H/Si(1 0 0)-(2 × 1). In addition to the oxygen inserted coupled monohydrides, which were previously reported in the former reaction system, we report several other oxidized Si hydride species in our BML-IRRAS experiments. Three new pairs of vibrational bands are identified between 900 and 1000 cm−1. These vibrational frequencies were calculated using Si9 and Si10 cluster models that included all possible structures from zero to five oxygen insertions into the top layer silicon atoms using a B3LYP gradient corrected density functional method with a polarized 6-31G** basis set for all atoms. The three pairs of vibrational modes are assigned to the scissoring modes of adjacent and isolated SiH2 with zero, one, and two oxygen atoms inserted into the Si back bonds. All the other newly observed vibrational peaks related to Si oxidation are also assigned in this study. The Si-O stretching bands observed in the reaction 2D + H2O/Si(1 0 0)-(2 × 1) show an isotope effect, which suggests that in the system 2H + H2O/Si(1 0 0)-(2 × 1) also, hydrogen atom tunneling plays an important role for the insertion of oxygen atoms into Si back bonds that form oxidized adjacent dihydrides.  相似文献   

4.
The conformational landscape of N-acetylcysteine (NAC) has been investigated by a combined experimental matrix-isolation FT-IR and theoretical methodology. This combination is a powerful tool to study the conformational behavior of relatively small molecules. Geometry optimizations at the HF/3-21 level resulted in 438 different geometries with an energy difference smaller than 22 kJ mol−1. Among these, six conformations were detected with a relative energy difference smaller than 10 kJ mol−1 at the DFT(B3LYP)/6-31++G∗∗ level of theory. These were finally subjected to MP2/6-31++G∗∗ optimizations which resulted in five minima. The vibrational and thermodynamical properties of these conformations were calculated at both the DFT and MP2 methodologies. Experimentally NAC was isolated in an argon matrix at 16 K after being sublimated at 323 K. The most stable MP2 form appeared to be dominant in the experimental spectra but the presence of three other conformations with ΔEMP2 < 10 kJ mol−1 was also demonstrated. The experimentally observed abundance of the H-bond containing conformations appeared to be in good accordance with the predicted MP2 value.  相似文献   

5.
Some initial acid‐catalytic reactions involved in the synthesis of the urea‐formaldehyde resin were theoretically investigated at B3LYP and MP2 levels with solvent effects included. The results suggest that the addition between urea and formaldehyde in neutral condition undergoes with a concerted mechanism represented by a four‐member ring transition state. For this reaction, a notable barrier (above 130 kJ/mol) was identified at all theoretical levels. The reactions between urea and different protonated forms of formaldehyde in acid solution were investigated. The reaction between protonated methanediol with urea can produce the methylol urea cation via an SN2 transition state with a lower barrier of 54.8 kJ/mol. With the mediation of a water molecule, the intra‐molecular proton transfer produces the stable methylol carbonium (NH2CONHCH2+), which plays an important role in the following formation of methylene and methylene ether linkages. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

6.
The adsorption energies for physisorption and the most stable chemisorption of CO2 on the neutral charge of perfect anatase [TiO2] (0 0 1) are −9.03 and −24.66 kcal/mol on the spin-unpolarized and −12.98 and −26.19 kcal/mol on the spin-polarized surface. The small activation barriers of 1.67 kcal/mol on the spin-unpolarized surface and of 6.66 kcal/mol on the spin-unpolarized surface were obtained. The adsorption mechanism of CO2 on the oxygen vacancy defect [TiO2 + VO] surface of anatase TiO2 using density functional theory calculations was investigated. The energetically preferred conversion of CO2 to CO was found either on the spin-unpolarized or spin-polarized surfaces of oxygen vacancy defect surface [TiO2 + VO] as photocatalyst.  相似文献   

7.
Surface segregation studies of Sn and Sb in Cu(1 0 0)-0.14 at.% Sn-0.12 at.% Sb ternary alloy, have been done by making use of Auger Electron Spectroscopy. The method of Linear Temperature Ramp (LTR) was employed, whereby the sample was heated and cooled linearly at a constant rate. The positive heating rate showed both a kinetic segregation profile, as well as a narrow equilibrium segregation region, at higher temperatures. The equilibrium segregation profile was extended by cooling the sample. Sn was first to segregate to the surface due to its higher diffusion coefficient, mainly from a smaller activation energy ESn. Sb, due to its higher segregation energy, eventually replaced Sn from the surface. The modified Darken model was used to simulate the profile yielding the following segregation parameters: Do(Sn) = 6.3 × 10−6 m2/s, Do(Sb) = 2.8 × 10−5 m2/s; ESn = 175.4 kJ/mol, ESb = 186.3 kJ/mol; , ; ΩCu-Sn = 3.4 kJ/mol, ΩCu-Sb = 15.9 kJ/mol and ΩSn-Sb = −5.4 kJ/mol.  相似文献   

8.
Adsorption of parent nitrosamine (NA) on 5T and 10T cluster models of H-ZSM-5 catalyst has been theoretically investigated using quantum chemical B3LYP and MP2 methods. Three stable complexes (A-C) were found on the potential energy surface of interaction between NA and cluster models of H-ZSM-5. NA can interact not only with acidic site of zeolite via the lone electron pair on nitrogen and oxygen atoms (O(N)?HZO) but also with the oxygen atoms of the framework via the hydrogen atoms of NH2 group (NH?O) as well. However, the Lewis acidity of zeolite framework is the dominating factor in the interaction between NA and zeolite. The calculated adsorption enthalpy of NA on 5T and 10T clusters of H-ZSM-5 catalyst at ONIOM(MP2/6-311++G(d,p):HF/6-31+G(d)) level ranges from −19.73 to −40.33 and −63.81 to −73.73 kJ/mol, respectively. Adsorption energy for A-C complexes increases in going from B3LYP method to MP2 one. The results of atoms in molecules (AIM) calculations showed that NH5?O interactions have electrostatic character, whereas O(N)?HZ interactions have partially covalent nature. The results of natural bond orbital (NBO) analysis showed that charge transfer occurs from NA to H-zeolite cluster.  相似文献   

9.
The isomerization of 1-naphthylacetylene diluted in argon was studied behind reflected shock waves in a 2 in i.d. single pulse shock tube over the temperature range 1000-1250 K and overall densities of ∼3 × 10−5 mol/cm3. The only reaction product found in the post shock mixtures was acenaphthylene. The first order rate constant of the isomerization was found to be k = 2.08 × 1012 exp(−54.2 × 103/RT) s−1, where R is expressed in units of cal/K mol. Potential energy surfaces of the cyclization reaction 1-naphthylacetylene → acenaphthylene and the isomerization 1-naphthylacetylene → 2-naphthylacetylene were calculated using the Becke three-parameter hybrid method with Lee-Yang-Parr correlation functional approximation (B3LYP). Structure, energy and frequency calculations were carried out with the Dunning correlation consistent polarized double ζ (cc-pVDZ) basis set. The rate constant (k) for the 1-naphthylacetylene → acenaphthylene cyclization was calculated using transition state theory, the value obtained is k = 3.52 × 1012 exp(−55.9 × 103/RT) s−1, where R is expressed in units of cal/K mol. The agreement between the experiment and the calculations is very good. RRKM calculations were done to transfer k to the pressure of the single pulse shock tube experiments. In view of high temperature and the large molecule involved the deviation from k is very small. The isomerization 1-naphthylacetylene → 2-naphthylacetylene proceeds via the formation of an unstable intermediate 1,2-naphthalenocyclobutene and has a high barrier of ∼83.5 kcal/mol. In view of this high barrier, the isomerization cannot compete with the cyclization that proceeds with a barrier of ∼56 kcal/mol.  相似文献   

10.
High resolution (0.001 cm−1) coherent anti-Stokes Raman spectroscopy (CARS) was used to directly examine the ν1 symmetric stretching mode of the planar symmetric D3h molecules 10BF3 and 11BF3. Simulations of the spectra were done using ν1 rovibrational parameters deduced from published infrared hot-band and difference-band studies and the close similarity to the observed CARS spectra confirms the validity of the infrared constants. No significant perturbations by Fermi resonance or Coriolis interactions with nearby states are observed, in marked contrast to the case of sulfur trioxide, a similar D3h molecule recently studied. In the harmonic approximation, the 10BF3 and 11BF3ν1 Q-branches would be identical since the isotopic substitution is at the center of mass but, interestingly, the ν1 stretching frequency for 11BF3 is found to be 0.198 cm−1higher than for the lighter 10BF3 isotopomer. This counterintuitive result is reproduced almost exactly (0.200 cm−1) by ab initio calculations (B3LYP/cc-pVTZ) that included evaluation of cubic and quartic force constants and xij anharmonicity constants. The ab initio computations also predict to within 1% the ΔB, ΔC changes in the rotational constants in going from the ground state to the v1 = 1 vibrational level. The results illustrate nicely the complementary interplay of modern infrared, Raman, and ab initio methods in obtaining and analyzing rovibrational spectra.  相似文献   

11.
Ethylene adsorption was studied by use of DFT/B3LYP with basis set 6-31G(d,p) in Gaussian’03 software. It was found that ethylene has adsorbed molecularly on all clusters with π adsorption mode. Relative energy values were calculated to be −50.86 kcal/mol, −20.48 kcal/mol, −32.44 kcal/mol and −39.27 kcal/mol for Ni13 nanocluster, Ni10(1 1 1), Ni13(1 0 0) and Ni10(1 1 0) surface cluster models, respectively. Ethylene adsorption energy is inversely proportional to Ni coordination number when Ni10(1 1 1), Ni13(1 0 0) and Ni10(1 1 0) cluster models and Ni13 nanocluster are compared with each other.  相似文献   

12.
The heats of adsorption at several coverages of the linear and bridged CO species (denoted L and B, respectively) adsorbed on the Pt0 sites of the 2.9 wt% Pt/10% K/Al2O3 catalyst are determined using the Adsorption Equilibrium Infrared spectroscopy method. The addition of K on 2.9% Pt/Al2O3 modifies significantly the adsorption of CO on the Pt particles: (a) the ratio L/B is decreased from 8.4 to 1, (b) a new adsorbed CO species is detected with an IR band at 1763 cm−1, (c) the heats of adsorption of L and B CO species are significantly altered and the positions of their IR bands are shifted. The heats of adsorption of L CO species are decreased: i.e. 206 and 105 kJ/mol at low coverages on Pt/Al2O3 and Pt/K/Al2O3 respectively. Two B CO species denoted B1 and B2, with different heats of adsorption are observed on Pt/K/Al2O3. The heats of adsorption of B2 CO species (major B CO species) are significantly larger than those measured in the absence of K: i.e. 94 and 160 kJ/mol at low coverages on Pt/Al2O3 and Pt/K/Al2O3 respectively, whereas those of B1 CO species (minor species) are similar: 90 kJ/mol at low coverages. These values are consistent with the qualitative High Resolution Electron Energy Loss Spectrometry literature data on Pt(1 1 1) modified by potassium.  相似文献   

13.
Middle infrared absorption, Raman scattering and proton magnetic resonance relaxation measurements were performed for [Zn(NH3)4](BF4) in order to establish relationship between the observed phase transitions and reorientational motions of the NH3 ligands and BF4 anions. The temperature dependence of spin-lattice relaxation time (T1(1H)) and of the full width at half maximum (FWHM) of the bands connected with ρr(NH3), ν2(BF4) and ν4(BF4) modes in the infrared and in the Raman spectra have shown that in the high temperature phase of [Zn(NH3)4](BF4)2 all molecular groups perform the following stochastic reorientational motions: fast (τR≈10−12 s) 120° flips of NH3 ligands about three-fold axis, fast isotropic reorientation of BF4 anions and slow (τR≈10−4 s) isotropic reorientation (“tumbling”) of the whole [Zn(NH3)4]2+ cation. Mean values of the activation energies for uniaxial reorientation of NH3 and isotropic reorientation of BF4 at phases I and II are ca. 3 kJ mol−1 and ca. 5 kJ mol−1, respectively. At phases III and IV the activation energies values for uniaxial reorientation of both NH3 and of BF4 equal to ca. 7 kJ mol−1. Nearly the same values of the activation energies, as well as of the reorientational correlation times, at phases III and IV well explain existence of the coupling between reorientational motions of NH3 and BF4. Splitting some of the infrared bands at TC2=117 K suggests reducing of crystal symmetry at this phase transition. Sudden narrowing of the bands connected with ν2(BF4), ν4(BF4) and ρr(NH3) modes at TC3=101 K implies slowing down (τR?10−10 s) of the fast uniaxial reorientational motions of the BF4 anions and NH3 ligands at this phase transition.  相似文献   

14.
New types of hybrid material have been synthesized by using four different methods of immobilization of humic acid (HA) on chitin. The most stable hybrid material toward the change of medium acidity was then utilized as sorbent for Cr(III).The HA was extracted from peat soil of Gambut District, South Kalimantan, Indonesia, using the recommended procedure of International Humic Substances Society (IHSS), while the chitin was isolated from crab shell waste through deproteination using 3.5% (w/v) NaOH and followed by removal of inorganic impurities using 1 M HCl. The four methods of immobilization of HA on chitin were (i) Method A: chitin powder (4 g) was gently poured into the stirred solution of 0.4 g HA in 40 mL of 0.01 M NaOH. After overnight stirring, the solid was separated, washed with water, and dried in oven at 70 °C. (ii) Method B: gelatinous chitin (40 g) in 250 mL of 0.5 M HCl was reacted with HA (4 g) in 500 mL of 0.5 M NaOH and aged for 24 h. The product was washed with water and dried. (iii) Method C: HA powder (0.5 g) was mixed with the stirred gel of chitin (2.5 g) in 60 mL of CaCl2 saturated methanol and the mixture was then washed with the mixed solution of 25 mL of 2 M sodium citrate and ethylene glycol 1:1. The solid was separated, washed with water, and dried. (iv) Method D: the solution of HA (0.056 g) in 10 mL of 0.01 M NaOH was reacted with the gel of chitin (0.2 g) in 10 mL of CaCl2 saturated methanol. After 24 h stirring, the solid was separated from the reaction medium, washed with the mixed solution of 2 M sodium citrate and ethylene glycol 1:1, and followed by washing with water and drying. Parameters investigated in this study consisted of the stability test of the immobilized HA, as well as the rate constant (k1), capacity (b), and energy (E) of sorption as well as the rate constant of desorption (k−1). The k1 and k−1 were determined according to a kinetic model of first order sorption reaching equilibrium, while the b and E were determined according to the Langmuir isotherm model.Compared to HA, Methods, A, C, and D; Method B produced the most stable immobilization of HA on chitin. The hybrid material (Chitin-HA) synthesized through Method B was stable in the acidity range that equivalent to pH 2.0-11.0. At the acidity giving maximum sorption, i.e. pH 5, the presence of immobilized HA on the Chitin-HA enhanced more than three times the k1 and k−1, i.e. from 0.057 min−1 and 8.51 × 10−4 (min−1) (mol/L) for chitin to 0.183 min−1 and 3.27 × 10−3 (min−1) (mol/L) for the Chitin-HA. On the contrary, the presence of HA on Chitin-HA only gave small increase on b and small decrease on E. The values of b and E for Cr(III) on chitin were 1.45 × 10−2 mol/g and 23.12 kJ/mol, respectively, while those on Chitin-HA were 1.78 × 10−2 mol/g and 19.95 kJ/mol, respectively.  相似文献   

15.
在B3LYP/6-31+G*水平上研究BF4-与甘氨酸间氢键作用特征,并在B3LYP/6-311++G**水平上计算单点能. 对二聚物几何结构、能量、氢键成键特征进行分析. 分子中的原子拓扑分析表明氢键成键原子间存在(3,-1)关键点,电子密度和Laplacian量落在氢键范围内. 进一步对氢键形成导致H原子净电荷、偶极矩、能量及体积的改变进行系统分析.  相似文献   

16.
The adsorption of 1,3-cyclohexadiene, 1,4-cyclohexadiene, cyclohexene and cyclohexane on Pt(1 1 1) was studied using ab initio density functional theory. For 1,3-cyclohexadiene three adsorption modes were distinguished: bridge 1,2-di-σ/3,4-π, hollow 1,4-di-σ/2,3-π and bridge 1,4-di-σ/2,3-π with adsorption energies of −155, −147 and −75 kJ/mol, respectively. Three stable adsorption modes were also identified for 1,4-cyclohexadiene: bridge quadra-σ, hollow di-σ/π and bridge di-π with adsorption energies of −146 kJ/mol, −142 kJ/mol and −88 kJ/mol, respectively. Cyclohexene was found to adsorb in six modes: 4 di-σ and 2 π-adsorption modes. The preferred configuration was found to be boat di-σ with an adsorption energy of −81 kJ/mol. The three other di-σ adsorption modes have comparable adsorption energies, ranging from −64 to −69 kJ/mol. Molecular strain and CPt bonding energies are used to elucidate stability trends. Cyclohexane is found to adsorb only at the hollow site whereby the axial hydrogen atoms are positioned over surface Pt-atoms with an adsorption energy of −37 kJ/mol. The calculations correctly predict the weakening of the axial CH bonds and provide a possible explanation for the large shift in the vibrational frequencies.  相似文献   

17.
High-resolution (0.0015-0.0035 cm−1) infrared spectra of isotopically enriched 11BF3 have been examined in detail. The analysis of the combination and overtone states within the region of study, from 1650 to 4600 cm−1, led to the assignment of over 25,000 transitions. The major perturbations were due to the Fermi resonances between states possessing one quantum of v3 and three quanta of v4. With corrections through the quadratic rotational terms, the equilibrium Be and Ce values have been determined; 0.3462679(7) cm−1 and 0.1731151(6) cm−1, respectively. An improved set of equilibrium rotational constants for 10BF3, consistent with this analysis of 11BF3 are also given. The averaged equilibrium values for both isotopomers lead to a B-F bond distance of re = 130.704 ± 0.005 pm. All of the quadratic anharmonic constants, with the exception of x33 were independently determined from experiment. For the first time for BF3, a normal force field analysis was performed that utilized the experimentally determined, fundamental harmonic vibrational frequencies.  相似文献   

18.
Over 8000 line positions and intensities of phosphine (PH3) at 3 μm have been measured at 0.0115 cm−1 resolution with the McMath-Pierce Fourier Transform spectrometer at Kitt Peak. The observed line intensities ranged from 4.13 × 10−6 to 4.69 × 10−2 cm−2 atm−1 at 296 K, for line positions between 2724.477 and 3601.652 cm−1. This region spans eight interacting vibrational states: 3ν2 (2940.8 cm−1), 2ν2 + ν4 (3085.6 cm−1), ν2 + 2ν4 (3214.9 cm−1), ν1 + ν2 (3307.6 cm−1), ν2 + ν3 (3310.5 cm−1), 3ν4 (∼3345 cm−1), ν1 + ν4 (3426.9 cm−1), and ν3 + ν4 (3432.9 cm−1). Assignments have been determined for all the bands except 3ν4 (a weak band in a highly congested area) for a total of 4232 transitions. The total integrated intensity for this region is 5.70 cm−2 atm−1 near 296 K, and assigned lines account for 79% of the observed absorption. The two strongest bands in the region are ν1 + ν4 and ν3 + ν4 with band strengths at 296 K of 1.61 and 2.01 cm−2 atm−1, respectively. An empirical database of PH3 line parameters (positions, intensities, and assignments) is now available. Lower state energies (corresponding to assignments from this study) and line widths from the literature are included; default values are used for unassigned features.  相似文献   

19.
The X2Πg, 2Σg+ and 2Δg states of AgBr2 have been studied through benchmark ab initio CASSCF + Averaged Coupled Pair Functional (ACPF) and DFT calculations using especially developed valence basis sets to study the transition energies, geometries, vibrational frequencies, Mulliken charges and spin densities. The spin-orbit (SO) effects were included through the effective hamiltonian formalism using the |ΛSΣ〉 ACPF energies as diagonal elements. At the ACPF level, the ground state is 2Πg, in contradiction with ligand-field theory and Hartree-Fock results. The ACPF adiabatic excitation energies of the 2Σg+ and 2Δg states are 3825 and 20 152 cm−1, respectively. The inclusion of the SO effects leads to a pure Ω = 3/2 (2Πg) ground state, a Ω = 1/2 (97% 2Πg + 3% 2Σg+) A state, a Ω = 1/2 (3% 2Πg + 97% 2Σg+) B state, a Ω = 5/2 (2Δg) C state and a Ω = 3/2 (99% 2Δg) D state. The B97, B3LYP and PBE0 functionals, which were shown to yield accurate transition energies for CuCl2, overestimate the X2Πg-2Σg+ Te by around 25% but provide a qualitative energetic ordering in agreement with CASSCF and ACPF results. The nature of the bonding in the X2Πg ground state is different from that of AgCl2 since the Mulliken charge on the metal is 0.95 while the spin density is only 0.39. DFT strongly delocalizes the spin density providing even smaller values of around 0.13 on Ag not only for the ground state, but also for the 2Σg+ state.  相似文献   

20.
To model the adsorption of Au+ cation in aqueous solution on the semiconductor surface, the interactions of Au+ and hydrated Au+ cations with clean Si(1 1 1) surface were investigated by using hybrid density functional theory (B3LYP) and Møller-Plesset second-order perturbation (MP2) methods. Si(1 1 1) surface was described with Si7H11, Si11H17 and Si22H21 clusters. The effect of the basis set superposition error (BSSE) was taken into account by applying the counterpoise correction. The calculated results indicated that the binding energies between hydrated Au+ cations and clean Si(1 1 1) surface are large, suggesting a strong interaction between hydrated Au+ cations and the semiconductor surface. The bonding nature of the chemical adsorption of Au+ to Si surface can be classified as partial covalent as well as ionic bonding. As the number of water molecules increases, the water molecules form hydrogen bond network with one another and only one water molecule binds directly to the Au+ cation. The Au+ cation in aqueous solution will safely attach to the clean Si(1 1 1) surface.  相似文献   

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

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