首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
A theoretical study of the mechanism and kinetics of the OH hydrogen abstraction from isopropylcyclopropane (IPCP) is presented. Optimum geometries, frequencies and gradients have been computed at the BHandHLYP/6-311++G(d,p) level of theory for all stationary points, as well as for additional points along the minimum energy path (MEP). Energies have been improved by single-point calculations at the above geometries using CCSD(T)/6-311++G(d,p) to produce the potential energy surface. The rate coefficients are calculated for the temperature range 260-350 K by using canonical variational theory (CVT) with small-curvature tunneling (SCT) corrections. Our analysis suggests a stepwise mechanism involving the formation of a reactant complex in the entrance channel and a product complex in the exit channel, for all the modeled paths. The reactant complexes are examined in detail, because they exhibit alkene-like structure. The excellent agreement between the overall calculated and experimental rate coefficients at 298 K supports the reliability of the parameters obtained for the temperature dependence and branching ratios of the IPCP + OH reaction, proposed here for the fist time. The expression that best describes the studied reaction is k(overall) = 6.15 x 10(-13)e1747/RT cm3 x molecule(-1) x s(-1). The predicted activation energy is -0.89 kcal/mol.  相似文献   

2.
A theoretical study of the mechanism and kinetics of the OH hydrogen abstraction from hydroxyacetone is presented. Optimum geometries and frequencies have been computed at the BH and HLYP/6-311++G(d,p) level of theory for all stationary points. Energy values have been improved by single-point calculations at the above geometries using CCSD(T)/ 6-311++G(d,p). The rate coefficients are calculated for the temperature range 280-500 K by using conventional transition state theory (TST), including tunneling corrections. Our analysis supports a stepwise mechanism involving the formation of a reactant complex in the entrance channel and a product complex in the exit channel, for all the modeled paths. Four experimental values of the rate constant at 298 K have been previously reported: three of them in great agreement (approximately 3 x 10(-12) cm(3) molecule(-1) s(-1)), and one of them twice larger. The calculations in the present work support the smaller value. A curved Arrhenius plot was found in the studied temperature range; thus the expression that best describes the obtained data is k(280-500)(overall) = 5.29 x 10(-23)T(3.4)e(1623/T) cm(3) molecule(-1) s(-1). The activation energy was found to vary with temperature from -1.33 to +0.15 kcal/mol.  相似文献   

3.
A theoretical study of the mechanism and kinetics of the OH hydrogen abstraction from glyoxal and methylglyoxal is presented. Optimum geometries, frequencies, and gradients have been computed at the BHandHLYP/6-311++G(d,p) level of theory for all the stationary points, as well as for 12 additional points along the minimum energy path (MEP). Energies were obtained by single-point calculations at the above geometries using CCSD(T)/ 6-311++G(d,p) to produce the potential energy surface. The rate coefficients were calculated for the temperature range 200-500 K by using canonical variational theory (CVT) with small-curvature tunneling (SCT) corrections. Our analysis suggests a stepwise mechanism, which involves the formation of a reactant complex. The overall agreement between the calculated and experimental kinetic data is very good. This agreement supports the reliability of the Arrhenius parameters of the glyoxal + OH reaction that are proposed in this work for the first time. The Arrhenius expressions that best describe the studied reactions are k1 = (9.63 +/- 0.23) x l0(-13)exp[(517 +/- 7)/T] and k2 = (3.93 +/- 0.11) x 10(-13)exp[(1060 +/- 8)/T]cm3 molecule(-1)s(-1) for glyoxal and methylglyoxal, respectively.  相似文献   

4.
We present a theoretical study of the hydrogen abstraction reactions from CH(3)F and CH(2)F(2) by an ozone molecule. The geometries and harmonic vibrational frequencies of all stationary points are calculated at the MPW1K, BHandHLYP, and MPWB1K levels of theory. The energies of all of the stationary points were refined by using both higher-level (denoted as HL) energy calculations and QCISD(T)/6-311++G(2df,2pd) calculations based on the optimized geometries at the MPW1K/6-31+G(d,p) level of theory. The minimum energy paths (MEPs) were obtained by the MPW1K/6-31+G(d,p) level of theory. Energetic information of the points along the MEPs is further refined by the HL method. The rate constants were evaluated on the basis of the MEPs from the HL level of theory in the temperature range 200-2500 K by using the conventional transition-state theory (TST), the canonical variational transition-state theory (CVT), the microcanonical variational transition-state theory (microVT), the CVT coupled with the small-curvature tunneling (SCT) correction (CVT/SCT), and the microVT coupled with the Eckart tunneling correction (microVT/Eckart) based on the ab initio calculations. A general agreement was found among the TST, CVT, and microVT theories. The fitted three-parameter Arrhenius expressions of the calculated forward CVT/SCT and microVT/Eckart rate constants of the ozonolysis of fluoromethane are k(CVT/SCT)(T) = 2.76 x 10(-34)T(5.81)e((-13975/)(T)) and k(microVT/Eckart)(T) = 1.15 x 10(-34)T(5.97)e((-14530.7/)(T)), respectively. The fitted three-parameter Arrhenius expressions of the calculated forward CVT/SCT and microVT/Eckart rate constants of the ozonolysis of difluoromethane are k(CVT/SCT)(T) = 2.29 x 10(-36)T(6.42)e((-15451.6/)(T)) and k(microVT/Eckart)(T) = 1.31 x 10(-36)T(6.45)e((-15465.8/)(T)), respectively.  相似文献   

5.
The complex potential energy surface and reaction mechanisms for the unimolecular isomerization and decomposition of methyl-nitramine (CH3NHNO2) were theoretically probed at the QCISD(T)/6-311+G*//B3LYP/6-311+G* level of theory. The results demonstrated that there are four low-lying energy channels: (i) the NN bond fission pathway; (ii) a sequence of isomerization reactions via CH3NN(OH)O; (IS2a); (iii) the HONO elimination pathway; (iv) the isomerization and the dissociation reactions via CH3NHONO (IS3). The rate constants of each initial step (rate-determining step) for these channels were calculated using the canonical transition state theory. The Arrhenius expressions of the channels over the temperature range 298-2000 K are k6(T)=1014:8e-46:0=RT , k7(T)=1013:7e-42:1=RT , k8(T)=1013:6e-51:8=RT and k9(T)=1015:6e-54:3=RT s-1, respectively. The calculated overall rate constants is 6.9£10-4 at 543 K, which is in good agreement with the experimental data. Based on the analysis of the rate constants, the dominant pathway is the isomerization reaction to form CH3NN(OH)O at low temperatures, while the NN bond fission and the isomerization reaction to produce CH3NHONO are expected to be competitive with the isomerization reaction to form CH3NN(OH)O at high temperatures.  相似文献   

6.
Diode infrared laser absorption spectroscopy was used to measure the rate constant (k(1)) of the CN + CS(2) reaction for the first time. k(1) was determined to be substantially pressure dependent with a value k(1) = (7.1 ± 0.2 to 41.9 ± 2.9) × 10(-12) cm(3) molecule(-1) s(-1) over 2-40 Torr at 298 K. The potential energy surface (PES) of the reaction was calculated using an ab initio method at B3LYP/6-311++G(d, p)//CCSD(T)/6-311++G(d, p) level of theory. Both experimental and computational results suggest that collision stabilization of the adduct NCSCS may dominate the reaction. The rate constant of the CN + SO(2) reaction was measured to be very slow with an upper limit of k(2) ≤ 3.1 × 10(-14) cm(3) molecule(-1) s(-1), in disagreement with an earlier reported measurement. The PES of this reaction reveals an entrance barrier against formation of the low energy adduct NCOSO, in agreement with the experimental result.  相似文献   

7.
The potential energy surface, including the geometries and frequencies of the stationary points, of the reaction HFCO + OH is calculated using the MP2 method with 6-31+G(d,p) basis set, which shows that the direct hydrogen abstraction route is the most dominating channel with respect to addition and substitution channels. For the hydrogen abstraction reaction, the single-point energies are refined at the QCISD(T) method with 6-311++G(2df,2pd) basis set. The calculated standard reaction enthalpy and barrier height are -17.1 and 4.9 kcal mol(-1), respectively, at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The reaction rate constants within 250-2500 K are calculated by the improved canonical variational transition state theory (ICVT) with small-curvature tunneling (SCT) correction at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The fitted three-parameter formula is k = 2.875 x 10(-13) (T/1000)1.85 exp(-325.0/T) cm(3) molecule(-1) s(-1). The results indicate that the calculated ICVT/SCT rate constant is in agreement with the experimental data, and the tunneling effect in the lower temperature range plays an important role in computing the reaction rate constants.  相似文献   

8.
High-level ab initio calculations have been performed to study the mechanism and kinetics of the reaction of the cyanomethylene radical (HCCN) with the NO. The species involved have been optimized at the B3LYP/6-311++G(3df,2p) level, and their corresponding single-point energies are improved by the CCSD(T)/aug-cc-PVQZ//B3LYP/6-311++G(3df,2p) approach. From the calculated potential energy surface, we have predicted the favorable pathways for the formation of several isomers of a HCCN-NO complex. Barrierless formation of HCN + NCO (P1) is also possible. Formation of HCNO + CN (P3) is endoergic but may become significant at high temperatures. To rationalize the scenario of our calculated results, we also employ the Fukui functions and hard-and-soft acid-and-base (HSAB) theory to seek possible clues. The predicted total rate coefficient, k(total), at He pressure 760 Torr can be represented with the equation k(total) = 1.40 × 10(-7) T(-2.01) exp(3.15 kcal mol(-1)/RT) at T = 298-3000 K in units of cm(3) molecule(-1) s(-1). The predicted total rate coefficients at some available conditions (He pressures of 6, 18, and 30 Torr in the temperature of 298 K) are in reasonable agreement with experimental observation. In addition, the rate constants for key individual product channels are provided in different temperature and pressure conditions.  相似文献   

9.
The mechanisms and kinetics of the reaction of a thiocyanato radical (NCS) with NO were investigated by a high-level ab initio molecular orbital method in conjunction with variational RRKM calculations. The species involved were optimized at the B3LYP/6-311++G(3df,2p) level, and their single-point energies were refined by the CCSD(T)/aug-cc-PVQZ//B3LYP/6-311+G(3df,2p) method. Our calculated results indicate favorable pathways for the formation of several isomers of an NCSNO complex. Formation of OCS + N 2 also is possible, although this pathway involves a substantial energy barrier. The predicted total rate constants, k total, at a 2 torr He pressure can be represented by the following equations: k total = 9.74 x 10 (26) T (-13.88) exp(-6.53 (kcal mol (-1))/ RT) at T = 298-950 K and 1.17 x 10 (-22) T (2.52) exp(-6.86 (kcal mol (-1))/ RT) at T = 960-3000 K, in units of cm (3) molecule (-1) s (-1), and the predicted values are in good agreement with the experimental results in the temperature range of 298-468 K. The calculated results clearly indicate that the branching ratio for R M1 in the temperature range of 298-950 K has the largest value ( R M1 accounts for 0.53-0.39). However, in the higher temperature range (960-3000 K), the formation of OCS + N 2 ( P5) with branching ratio R P5 (0.40-0.79) becomes dominant. The rate constants for key individual product channels are provided for different temperature and pressure conditions.  相似文献   

10.
The kinetics and mechanism of the reaction of the cyanomidyl radical (HNCN) with the hydroxyl radical (OH) have been investigated by ab initio calculations with rate constants prediction. The single and triplet potential energy surfaces of this reaction have been calculated by single-point calculations at the CCSD(T)/6-311+G(3df,2p) level based on geometries optimized at the B3LYP/6-311+G(3df,2p) and CCSD/6-311++G(d,p) levels. The rate constants for various product channels in the temperature range of 300-3000 K are predicted by variational transition-state and Rice-Ramsperger-Kassel-Marcus (RRKM) theories. The predicted total rate constants can be represented by the expressions ktotal=2.66 x 10(+2)xT-4.50 exp(-239/T) in which T=300-1000 K and 1.38x10(-20)xT2.78 exp(1578/T) cm3 molecule(-1) s(-1) where T=1000-3000 K. The branching ratios of primary channels are predicted: k1 for forming singlet HON(H)CN accounts for 0.32-0.28, and k4 for forming singlet HONCNH accounts for 0.68-0.17 in the temperature range of 300-800 K. k2+k7 for producing H2O+NCN accounts for 0.55-0.99 in the high-temperature range of 800-3000 K. The branching ratios of k3 for producing HCN+HNO, k6 for producing H2N+NCO, k8 for forming 3HN(OH)CN, k9 for producing CNOH+3NH, and k5+k10 for producing NH2+NCO are negligible. The rate constants for key individual product channels are provided in a table for different temperature and pressure conditions.  相似文献   

11.
CH3(2A′)自由基与臭氧反应机理的量子化学研究   总被引:2,自引:0,他引:2  
用量子化学UMP2方法,在6-311++G**基组水平上研究了CH3(2A′)自由基与臭氧反应机理,全参数优化了反应过程中反应物、中间体、过渡态和产物的几何构型,在UQCISD(T)/6-311++G**水平上计算了它们的能量;并对它们进行了振动分析,以确定中间体和过渡态的真实性;同时应用经典过渡态理论计算了反应的速率常数,并与实验值进行了比较, CH3自由基与臭氧反应速率常数的理论计算结果为: 4.73×10-14 cm3•molecule-1•s-1,与实验报导的结果(k=2.52×10-14 cm3•molecule-1•s-1)很接近,同时发现CH3(2A′)自由基与O3的反应是强放热反应.  相似文献   

12.
The mechanisms of the SH (SD) radicals with Cl2 (R1), Br2 (R2), and BrCl (R3) are investigated theoretically, and the rate constants are calculated using a dual-level direct dynamics method. The optimized geometries and frequencies of the stationary points are calculated at the MP2/6-311G(d,p) and MPW1K/6-311G(d,p) levels. Higher-level energies are obtained at the approximate QCISD(T)/6-311++G(3df, 2pd) level using the MP2 geometries as well as by the multicoefficient correlation method based on QCISD (MC-QCISD) using the MPW1K geometries. Complexes with energies less than those of the reactants or products are located at the entrance or the exit channels of these reactions, which indicate that the reactions may proceed via an indirect mechanism. The enthalpies of formation for the species XSH/XSD (X = Cl and Br) are evaluated using hydrogenation working reactions method. By canonical variational transition-state theory (CVT), the rate constants of SH and SD radicals with Cl2, Br2, and BrCl are calculated over a wide temperature range of 200-2000 K at the a-QCISD(T)/6-311++G(3df, 2pd)//MP2/6-311G(d, p) level. Good agreement between the calculated and experimental rate constants is obtained in the measured temperature range. Our calculations show that for SH (SD) + BrCl reaction bromine abstraction (R3a or R3a') leading to the formation of BrSH (BrSD) + Cl in a barrierless process dominants the reaction with the branching ratios for channels 3a and 3a' of 99% at 298 K, which is quite different from the experimental result of k3a'/k3' = 54 +/- 10%. Negative activation energies are found at the higher level for the SH + Br2 and SH + BrCl (Br-abstraction) reactions; as a result, the rate constants show a slightly negative temperature dependence, which is consistent with the determination in the literature. The kinetic isotope effects for the three reactions are "inverse". The values of kH/kD are 0.88, 0.91, and 0.69 at room temperature, respectively, and they increase as the temperature increases.  相似文献   

13.
The hydrogen abstract reactions of OH radicals with HOF (R1), HOCl (R2), and HOBr (R3) have been studied systematically by a dual-level direct-dynamics method. The geometries and frequencies of all the stationary points are optimized at the MP2/6-311+G(2d, 2p) level of theory. A hydrogen-bonded complex is located at the product channel for the OH + HOBr reaction. To improve the energetics information along the minimum energy path (MEP), single-point energy calculations are carried out at the CCSD(T)/6-311++G(3df, 3pd) level of theory. Interpolated single-point energy (ISPE) method is employed to correct the energy profiles for the three reactions. It is found that neither the barrier heights (DeltaE) nor the H-O bond dissociation energies [D(H-O)] exhibit any clear-cut linear correlations with the halogen electronegative. The decrease of DeltaE and D(H-O) for the three reactions are in order of HOF > HOBr > HOCl. Rate constants for each reaction are calculated by canonical variational transition-state theory (CVT) with a small-curvature tunneling correction (SCT) within 200-2000 K. The agreement of the rate constants with available experimental values for reactions R2 and R3 at 298 K is good. Our results show that the variational effect is small while the tunneling correction has an important contribution in the calculation of rate constants in the low-temperature range. Due to the lack of the kinetic data of these reactions, the present theoretical results are expected to be useful and reasonable to estimate the dynamical properties of these reactions over a wide temperature range where no experimental value is available.  相似文献   

14.
A direct dynamics method is employed to study the kinetics of the multiple channel reaction CH(3)OCl + Cl. The potential energy surface (PES) information is explored from ab initio calculations. Two reaction channels, Cl- and H-abstractions, have been identified. The optimized geometries and frequencies of the stationary points and the minimum-energy paths (MEPs) are calculated at the MP2 level of theory using the 6-311G(d, p) and cc-pVTZ basis sets, respectively. The single-point energies along the MEPs are further refined at the G3(MP2)//MP2/6-311G(d, p), G3//MP2/6-311G(d, p), as well as by the multicoefficient correlation method based on QCISD (MC-QCISD) using the MP2/cc-pVTZ geometries. The enthalpies of formation for the species CH(3)OCl and CH(2)OCl are calculated via isodesmic reactions. The rate constants of the two reaction channels are evaluated by using the variational transition-state theory over a wide range of temperature, 220-2200 K. The calculated rate constants exhibit the slightly negative temperature dependence and show good agreement with the available experimental data at room temperature at the G3(MP2)//MP2/6-311G(d, p) level. The present calculations indicate that the two channels are competitive at low temperatures while H-abstraction plays a more important role with the increase of temperature. The calculated k(1a)/k(1) ratio of 0.5 at 298 K is in general agreement with the experimental one, 0.8 +/- 0.2. The high rate constant for CH(3)OCl + Cl shows that removal by reaction with Cl atom is a potentially important loss process for CH(3)OCl in the polar stratosphere.  相似文献   

15.
The kinetics of the NCCO + NO(2) reaction was studied by transient infrared laser absorption spectroscopy. The total rate constant of the reaction was measured to be k = (2.1 ± 0.1) × 10(-11) cm(3) molecule(-1) s(-1) at 298 K. Detection of products and consideration of possible secondary chemistry shows that CO(2) + NO + CN is the primary product channel. The rate constants of the NCCO + CH(4) and NCCO + C(2)H(4) reactions were also measured, obtaining upper limits of k (NCCO + CH(4)) ≤ 7.0 × 10(-14) cm(3) molecule(-1) s(-1) and k (NCCO + C(2)H(4)) ≤ 5.0 × 10(-15) cm(3) molecule(-1) s(-1). Ab initio calculations on the singlet and triplet potential energy surfaces at B3LYP/6-311++G**//CCSD(T)/6-311++G** levels of theory show that the most favorable reaction pathway occurs on the singlet surface, leading to CO(2) + NO + CN products, in agreement with experiment.  相似文献   

16.
The multiple-channel reactions OH + CH3NHC(O)OCH3 --> products are investigated by direct dynamics method. The optimized geometries, frequencies, and minimum energy path are all obtained at the MP2/6-311+G(d,p) level, and energetic information is further refined by the BMC-CCSD (single-point) method. The rate constants for every reaction channels, R1, R2, R3, and R4, are calculated by canonical variational transition state theory with small-curvature tunneling correction over the temperature range 200-1000 K. The total rate constants are in good agreement with the available experimental data and the two-parameter expression k(T) = 3.95 x 10(-12) exp(15.41/T) cm3 molecule(-1) s(-1) over the temperature range 200-1000 K is given. Our calculations indicate that hydrogen abstraction channels R1 and R2 are the major channels due to the smaller barrier height among four channels considered, and the other two channels to yield CH3NC(O)OCH3 + H2O and CH3NHC(O)(OH)OCH3 + H2O are minor channels over the whole temperature range.  相似文献   

17.
Reactions of hydroxyl (OH) radicals with 1-butene (k(1)), trans-2-butene (k(2)), and cis-2-butene (k(3)) were studied behind reflected shock waves over the temperature range 880-1341 K and at pressures near 2.2 atm. OH radicals were produced by shock-heating tert-butyl hydroperoxide, (CH(3))(3)-CO-OH, and monitored by narrow-line width ring dye laser absorption of the well-characterized R(1)(5) line of the OH A-X (0, 0) band near 306.7 nm. OH time histories were modeled using a comprehensive C(5) oxidation mechanism, and rate constants for the reaction of OH with butene isomers were extracted by matching modeled and measured OH concentration time histories. We present the first high-temperature measurement of OH + cis-2-butene and extend the temperature range of the only previous high-temperature study for both 1-butene and trans-2-butene. With the potential energy surface calculated using CCSD(T)/6-311++G(d,p)//QCISD/6-31G(d), the rate constants and branching fractions for the H-abstraction channels of the reaction of OH with 1-butene were calculated in the temperature range 300-1500 K. Corrections for variational and tunneling effects as well as hindered-rotation treatments were included. The calculations are in good agreement with current and previous experimental data and with a recent theoretical study.  相似文献   

18.
冀永强  冯文林  徐振峰  雷鸣  郝茂荣 《化学学报》2001,59(12):2099-2104
采用DFT(B3LYP)方法,分别在6-311g(d,p),6-311++g(d,p)和自洽相关基组cc-pVIZ水平上优化了基态硝基甲烷和自由基H,OH,CH3,CH2[^3B1]以及O[^3P]等发生吸氢反应时的过渡态结构,并计算了反应的位垒。研究表明,对同一反应,不同基组下优化得到的过渡态几何结构基本一致;反应位垒数值的大小也基本接近,经校正,硝基甲烷同自由基反应位垒的理论计算值同实验结果基本吻合。  相似文献   

19.
The reaction C(2)H(5) + HBr --> C(2)H(6) + Br has been theoretically studied over the temperature range from 200 to 1400 K. The electronic structure information is calculated at the BHLYP/6-311+G(d,p) and QCISD/6-31+G(d) levels. With the aid of intrinsic reaction coordinate theory, the minimum energy paths (MEPs) are obtained at the both levels, and the energies along the MEP are further refined by performing the single-point calculations at the PMP4(SDTQ)/6-311+G(3df,2p)//BHLYP and QCISD(T)/6-311++G(2df,2pd)//QCISD levels. The calculated ICVT/SCT rate constants are in good agreement with available experimental values, and the calculate results further indicate that the C(2)H(5) + HBr reaction has negative temperature dependence at T < 850 K, but clearly shows positive temperature dependence at T > 850 K. The current work predicts that the kinetic isotope effect for the title reaction is inverse in the temperature range from 200 to 482 K, i.e., k(HBr)/k(DBr) < 1.  相似文献   

20.
Arrhenius rate expressions were determined for beta-scission of phenoxyl radical from 1-phenyl-2-phenoxyethanol-1-yl, PhC*(OH)CH2OPh (V). Ketyl radical V was competitively trapped by thiophenol to yield PhCH(OH)CH2OPh in competition with beta-scission to yield phenoxyl radical and acetophenone. A basis rate expression for hydrogen atom abstraction by sec-phenethyl alcohol, PhC*(OH)CH3, from thiophenol, log(k(abs)/M(-1) s(-1)) = (8.88 +/- 0.24) - (6.07 +/- 0.34)/theta, theta = 2.303RT, was determined by competing hydrogen atom abstraction with radical self-termination. Self-termination rates for PhC*(OH)CH3 were calculated using the Smoluchowski equation employing experimental diffusion coefficients of the parent alcohol, PhCH(OH)CH3, as a model for the radical. The hydrogen abstraction basis reaction was employed to determine the activation barrier for the beta-scission of phenoxyl from 1-phenyl-2-phenoxyethanol-1-yl (V): log(k beta)/s(-1)) = (12.85 +/- 0.22) - (15.06 +/- 0.38)/theta, k beta (298 K) ca. (64.0 s(-1) in benzene), and log(k beta /s(-1)) = (12.50 +/- 0.18) - (14.46 +/- 0.30)/theta, k beta (298 K) = 78.7 s(-1) in benzene containing 0.8 M 2-propanol. B3LYP/cc-PVTZ electronic structure calculations predict that intramolecular hydrogen bonding between the alpha-OH and the -OPh leaving group of ketyl radical (V) stabilizes both ground- and transition-state structures. The computed activation barrier, 14.9 kcal/mol, is in good agreement with the experimental activation barrier.  相似文献   

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

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