首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 468 毫秒
1.
CO2二聚体分子弱结合作用的DFT计算   总被引:4,自引:0,他引:4  
用密度泛函理论(DFT)的Becke 3LYP方法,在不同基集合(6 31G和6 311G系列)下对平行结构(C 2h)和T形结构(C2v)的CO2二聚体进行ab initio计算.通过计算,得到了CO2二聚体C2h和C2v两种构型的结构参数和离解能,并给出了CO2二聚体相对稳定构型C2h的12个正则振动分析图.结果表明,CO2二聚体的离解能为2 kJ•mol-1,CO2分子之间振动频率很小,从而说明CO2二聚体是弱结合分子.  相似文献   

2.
用量子化学计算方法对CH3CH=·CH与O2气的反应机理进行了理论研究, 在B3LYP/6-311G(d,p) 水平下优化稳定分子结构和寻找过渡态, 并在此构型的基础上, 采用CCSD(T)/6-311G(d,p)方法得到各驻点的高级单点能量. 找到主要路径R(CH3CH=·CH+O2)→m1(trans-CH3CH=CHOO)→m2(形成COO三元环)→m3(C—C键断裂,同时生成CH3CH—O—CHO)→P2(C—O键断裂生成CH3CHO+CHO); 并与C2H3等共轭体系进行了对比.  相似文献   

3.
用ab initio MRSDCI/6-311G(2 df,2 Pd)方法研究了Li_2H分子电子基态的势能面,计算了285个几何构型点的势能值,并采用Simons-Parr-Finlan展开式对这些势能值进行了拟合,得到均方差X~2等于4.64×10~(-6)(hartree~2).Li_2H分子电子基态的平衡几何构型为R_e=0.172nm,相似文献   

4.
三氯醋酸钕与8-羟基喹啉配合物的热化学研究   总被引:8,自引:1,他引:8  
用溶解量热法 ,在自行研制的具有恒定温度环境的新型反应量热计中 ,测定了NdCl3·6H2 O(s)、CCl3COOH(s)和Nd(TCA) 3·3H2 O(s)在1mol/LHCl中的溶解焓。再根据盖斯定律设计了一个热化学循环 ,计算得到了六水合氯化钕与三氯醋酸反应的反应焓ΔrHmθ(2 98.15K) =2 0 1.688kJ/mol,并求出了Nd(TCA) 3·3H2 O(s)的标准生成焓ΔfHmθ[Nd(TCA) 3·3H2 O ,s ,2 98.15K ] =-3 0 5 3 .3kJ/mol。同时测定了Nd(TCA) 3·3H2 O(s) ,C9H7NO(s) ,Nd(TCA) (C9H6 NO) 2 (s)和CCl3COOH(s)在 4mol/LHCl、二甲亚砜和无水乙醇混合溶剂中的溶解焓 ,再根据盖斯定律设计了一个热化学循环 ,计算得到了三氯醋酸钕与 8 羟基喹啉反应的反应焓ΔrHmθ(2 98.15K) =-3 .2 2 6kJ/mol,并求出了Nd(TCA) (C9H6 NO) 2 (s)的标准生成焓ΔfHmθ[Nd(TCA) (C9H6 NO) 2 ,s ,2 98.15K] =-13 5 5 .6kJ/mol。  相似文献   

5.
用 p H滴定法测定了标题化合物在 2 5± 0 .1℃、3 5± 0 .1℃和 4 5± 0 .1℃ KNO3 存在下的离解常数 ;计算了 2 5± 0 .1℃离解反应的焓变Δr Hm和熵变Δr Sm。确定了 2 5± 0 .1℃ KNO3 存在下 ,它与 Co( )、Ni( )及Cu( )形成的配合物的稳定常数  相似文献   

6.
应用规范不变原子轨道法(GIAO)在RHF/6-31G**和B3LYP/6-31G**水平上计算了质子化双氮桥联1,10-菲咯啉大环化合物(H4HAPP2+)C2h和C2h构型的1HNMR,并用TDDFT法计算了H4HAPP2+电子光谱.结果表明,B3LYP/6-31G*优化的C2h构型为较优构型,经谐振频率验证无虚频,C2h构型是H4HAPP2+合理的对称性构型.  相似文献   

7.
采用密度泛函理论(DFT)的B3LYP泛函, 在6-311G*水平上对B2Cn+(n=1~9)团簇的几何构型和电子结构进行了优化和振动频率计算. 结果表明, 在B2Cn+(n=1~9)团簇的基态构型中, B2C2+、B2C3+为具有D∞h对称性的线形结构, B2C7+为具有Cs对称性的立体环状结构, 其余均为平面构型; 其成键顺序为C—C成键优于B—C 成键, B—C成键优于B—B成键. 进一步得到了B2Cn+(n=1~9)团簇的总能量(ET)、零点能(EZ)、摩尔热容(Cp)、标准熵(S0)以及原子化能(ΔEn+). 其结果显示, 随着n的递增, ET、EZ、Cp、S0和ΔEn+数值均呈现增大趋势, 其中EZ数值呈现近似等梯度的增加趋势. 通过对B2Cn+(n=1~9)团簇基态结构的垂直电子亲合势的研究发现, n为奇数的B2Cn+团簇比n为偶数的稳定.  相似文献   

8.
宁攀  赵建想 《化学研究》2013,(5):493-500
应用密度泛函理论,在B3LYP/6-31G(d)和CAM-B3LYP/6-31G(d)水平上优化了α-联噻吩体系H(C4H2S)n H(n=213)的基态几何构型;与此同时,利用TD-B3LYP方法计算了H(C4H2S)n H的吸收光谱,得到了其垂直激发能和体系大小n的解析表达式;并采用TD-CAM-B3LYP方法研究了其发射光谱.研究表明,H(C4H2S)n H的基态结构呈现三种构型:螺旋上升型、环型和稍有弯曲的带状结构;其中前两种构型是顺式结构,最后一种构型是反式结构;环型结构具有C2对称性,其他两种结构具有C1对称性.此外,其吸收光谱和发射光谱计算值与实验值吻合.  相似文献   

9.
第1题1-1.A:Na2WO4;B:H2WO4或WO3·H2O;C:WO3⑴4FeWO4(s)+8NaOH(l)+O2(g)4Na2WO4(l)+2Fe2O3(s)+4H2O(g)2MnWO4(s)+4NaOH(l)+O2(g)2Na2WO4(l)+2MnO2(s)+2H2O(g)(产物写成Mn2O3或Mn3O4也可)⑵Na2WO4(aq)+2HCl(aq)H2WO4(s)+2NaCl(aq)⑶H2WO4(s)WO3(s)+H2O(g)1-2.[1]2WO3(s)+3C(石墨)2W(s)+3CO2(g)[2]反应在298.15K时的标准焓变、标准熵变和标准自由能变化分别为:ΔH=ΣυΔfH m=[2×0+3×(-393.51)-2×(-842.87)-3×0]kJ/mol=505.21 kJ/molΔS=ΣυS=[2×32.64+3×213.64-2×75.90-3×5.74]J·mol-1·K-1=537.…  相似文献   

10.
姜黄素衍生物电子结构及光谱性质的理论研究   总被引:1,自引:0,他引:1  
在B3LYP/6-31G**和TD B3LYP/6-31G**水平上对标题化合物电子结构和光谱性质进行了研究,为模拟真实条件,计算还考虑了甲醇溶液对电子结构和光谱性质的影响。研究结果表明:标题化合物存在A、B、C三类互变异构体,异构体稳定性为BAC。以ΔU、ΔH作为比较依据,则在气相和甲醇溶液中,A3、B3、C3分别为其异构体的能量最低构象,且B3为最稳定构型。以ΔG作为比较依据,气相中B3为最稳定构型;而甲醇溶液中B4为最稳定构型。异构体A的红外光谱在2503 cm-1处存在特征吸收峰,该吸收峰为其分子内缔合氢键的νO-H;异构体B的红外光谱在1627 cm-1处存在很强的特征吸收峰,该吸收峰为螯合羰基的νC=O。甲醇溶液使分子的ELUMO、EHOMO、ΔEgap均下降,电子光谱λmax发生红移。异构体A、B在气相和甲醇溶液中的λmax均都主要来源于H→L的π→π*跃迁;异构体C在气相中的λmax主要来源于H→L+1、H-1→L的π→π*跃迁,而在甲醇溶液中则主要来源于H-1→L的π→π*跃迁;以上跃迁存在分子内电荷转移现象。  相似文献   

11.
Ozone water reaction including a complex was studied at the MP2/6-311++G(d,p) and CCSD/6-311++G(2df,2p)//MP2/6-311++G(d,p) levels of theory. The interaction between water oxygen and central oxygen of ozone produces stable H2O-O3 complex with no barrier. With decomposition of this complex through H-abstraction by O3 and O-abstraction by H2O, three possible product channels were found. Intrinsic reaction coordinate, topological analyses of atom in molecule, and vibrational frequency calculation have been used to confirm the preferred mechanism. Thermodynamic data at T = 298.15 K and atmospheric pressure have been calculated. The results show that the production of hydrogen peroxide is the main reaction channel with ΔG = ?21.112 kJ mol-1.  相似文献   

12.
Y M Choi  J Park  M C Lin 《Chemphyschem》2004,5(5):661-668
The kinetics and mechanism of the reaction of C6H5 with CH3CHO have been investigated experimentally and theoretically. The total rate constant for the reaction has been measured by means of the cavity ring-down spectrometry (CRDS) in the temperature range 299-501 K at pressures covering 20-75 Torr. The overall bimolecular rate constant can be represented by the expression k = (2.8 +/- 0.2) x 10(11) exp[-(700 +/- 30)/T] cm3 mol-1 s-1, which is slightly faster than for the analogous C6H5 + CH2O reaction determined with the same method in the same temperature range. The reaction mechanism for the C6H5 + CH3CHO reaction was also explored with quantum-chemical calculations at various hybrid density functional theories (DFTs) and using ab initio high-level composite methods. The theories predict that the reaction may occur by two hydrogen-abstraction and two addition channels with the aldehydic hydrogen-abstraction reaction being dominant. The rate constant calculated by the transition state theory for the aldehydic hydrogen-abstraction reaction is in good agreement with the experimental result after a very small adjustment of the predicted reaction barrier (+0.3 kcal mol-1). Contributions from other product channels are negligible under our experimental conditions. For combustion applications, we have calculated the rate constants for key product channels in the temperature range of 298-2500 K under atmospheric-pressure conditions; they can be represented by the following expressions in units of cm 3mol-1 s-1: k1,cho = 8.8 x 10(3)T2.6 exp(-90/T), k2,ch3 = 6.0 x 10(1)T3.3 exp(-950/T), k3a(C6H5COCH3 + H) = 4.2 x 10(5)T0.6 exp(-410/T) and k3b(C6H5CHO + CH3) = 6.6 x 10(9)T-0.5 exp(-310/T).  相似文献   

13.
The kinetics and equilibrium of the allyl radical reaction with molecular oxygen have been studied in direct measurements using temperature-controlled tubular flow reactor coupled to a laser photolysis/photoionization mass spectrometer. In low-temperature experiments (T = 201-298 K), association kinetics were observed, and the measured time-resolved C(3)H(5) radical signals decayed exponentially to the signal background. In this range, the determined rate coefficients exhibited a negative temperature dependence and were observed to depend on the carrier-gas (He) pressure {p = 0.4-36 Torr, [He] = (1.7-118.0) × 10(16) cm(-3)}. The bimolecular rate coefficients obtained vary in the range (0.88-11.6) × 10(-13) cm(3) s(-1). In higher-temperature experiments (T = 320-420 K), the C(3)H(5) radical signal did not decay to the signal background, indicating equilibration of the reaction. By measuring the radical decay rate under these conditions as a function of temperature and following typical second- and third-law procedures, plotting the resulting ln K(p) values versus 1/T in a modified van't Hoff plot, the thermochemical parameters of the reaction were extracted. The second-law treatment resulted in values of ΔH(298)° = -78.3 ± 1.1 kJ mol(-1) and ΔS(298)° = -129.9 ± 3.1 J mol(-1) K(-1), with the uncertainties given as one standard error. When results from a previous investigation were taken into account and the third-law method was applied, the reaction enthalpy was determined as ΔH(298)° = -75.6 ± 2.3 kJ mol(-1).  相似文献   

14.
朱维晃  吴丰昌  黄廷林 《色谱》2008,26(5):550-553
通过高效液相色谱法研究了3-(2-吡啶基)-5,6-二苯基-1,2,4-三嗪(PDT)和Fe(Ⅱ)的配合物[Fe(PDT)3]2+的面式和经式两种几何异构体之间的动力学平衡过程。结果表明:不同温度(30,35,40,45 ℃)下,两种几何异构体含量(x)之间的相互转变均符合动力学一级反应,其xeln[(xe-x0)/(xe-x)]值和反应时间t(min)之间的关系分别为:xeln[(xe-x0)/(xe-x)]=0.082t+0.729 (r2=0.9911,T=45 ℃),xeln[(xe-x0)/(xe-x)]=0.049t+0.598 (r2=0.9987,T=40 ℃),xeln[(xe-x0)/(xe-x)]=0.022t+0.586 (r2=0.9987,T=35 ℃),xeln[(xe-x0)/(xe-x)]=0.012t+0.591(r2=0.9988,T=30 ℃)。两种异构体之间的动力学相互转变过程中的活化焓(ΔH)、活化熵(ΔS)和活化能(ΔEa)分别为:ΔH=103.84 kJ·mol-1,ΔS=271.93 J·mol-1·K-1,ΔEa=86.74 kJ·mol-1 (面式异构体向经式异构体转变);ΔH=106.47 kJ·mol-1,ΔS=257.65 J·mol-1·K-1,ΔEa=94.43 kJ·mol-1 (经式异构体向面式异构体转变)。  相似文献   

15.
Two new isostructural two-dimensional (2D) coordination polymers exhibiting spin crossover (SCO) behavior of formulation [Fe(4,4'-bipy)(2)(NCX)(2)]·4CHCl(3) (4,4'-bipy = 4,4'-bipyridine; X = S [1·4CHCl(3)], Se [2·4CHCl(3)]) have been synthesized and characterized, and both undergo cooperative spin transitions (ST). For 1·4CHCl(3) the ST takes place in two steps with critical temperatures of T(c1)(down) = 143.1 K, T(c2)(down) = 91.2 K, T(c1)(up) = 150.7 K, and T(c2)(up) = 112.2 K. 2·4CHCl(3) displays half ST characterized by T(c)(down) = 161.7 K and T(c)(up) = 168.3 K. The average enthalpy and entropy variations and cooperativity parameters associated with the ST have been estimated to be ΔH(1)(av) = 5.18 kJ mol(-1), ΔS(1)(av) = 35 J K(-1) mol(-1), and Γ(1) = 2.8 kJ mol(-1) and ΔH(2)(av) = 3.55 kJ mol(-1), ΔS(2)(av) = 35 J K(-1) mol(-1), and Γ(2) = 2.6 kJ mol(-1) for 1·4CHCl(3), and ΔH(av) = 6.25 kJ mol(-1), ΔS(av) = 38.1 J K(-1) mol(-1), and Γ = 3.2 kJ mol(-1) for 2·4CHCl(3). At T > [T(c1) (1·4CHCl(3)); T(c) (2·4CHCl(3))], both compounds are in the space group P2/c while at T < [T(c1) (1·4CHCl(3)); T(c) (2·4CHCl(3))] they change to the C2/c space group and display an ordered checkerboard-like arrangement of iron(II) sites where the high- and low-spin states coexist at 50%.  相似文献   

16.
在(298.15 ±0.01) K下用转动弹热量计测定了离子液体硫酸乙酯-1-甲基-3-乙基咪唑(EMIES)及合成它的原料1-甲基咪唑的恒容燃烧热,通过计算得到它们的标准燃烧焓 分别为(-2671±2) 和(-286.3±0.5) kJ·mol-1;标准生成焓 分别为(-3060±3) kJ·mol-1和(-2145±4) kJ·mol-1.结合文献上硫酸二乙酯的标准生成焓数据,得到了合成离子液体EMIES的反应热(-102.3±1.0) kJ·mol-1,与合成实验中观察到的强烈放热现象是一致的.根据离子液体EMIES的热容数据,计算了不同温度下EMIES的标准生成焓.  相似文献   

17.
o-Xylene could be a good candidate to represent the family of aromatic hydrocarbons in a surrogate fuel. This study uses computational chemistry to calculate standard enthalpies of formation at 298 K, Δ(f)H°(298 K), standard entropies at 298 K, S°(298 K), and standard heat capacities C(p)°(T) over the temperature range 300 K to 1500 K for ten target species present in the low-temperature oxidation mechanism of o-xylene: o-xylene (1), 2-methylbenzyl radical (2), 2-methylbenzylperoxy radical (3), 2-methylbenzyl hydroperoxide (4), 2-(hydroperoxymethyl)benzyl radical (5), 2-(hydroperoxymethyl)benzaldehyde (6), 1-ethyl-2-methylbenzene (7), 2,3-dimethylphenol (8), 2-hydroxybenzaldehyde (9), and 3-hydroxybenzaldehyde (10). Δ(f)H°(298 K) values are weighted averages across the values calculated using five isodesmic reactions and five composite calculation methods: CBS-QB3, G3B3, G3MP2, G3, and G4. The uncertainty in Δ(f)H°(298 K) is also evaluated. S°(298 K) and C(p)°(T) values are calculated at B3LYP/6-311G(d,p) level of theory from molecular properties and statistical thermodynamics through evaluation of translational, rotational, vibrational, and electronic partition functions. S°(298 K) and C(p)°(300 K) values are evaluated using the rigid-rotor-harmonic-oscillator model. C(p)°(T) values at T ≥ 400 K are calculated by treating separately internal rotation contributions and translational, external rotational, vibrational, and electronic contributions. The thermochemical properties of six target species are used to develop six new additivity groups taking into account the interaction between two substituents in ortho (ORT/CH2OOH/ME, ORT/ET/ME, ORT/CHO/OH, ORT/CHO/CH2OOH) or meta (MET/CHO/OH) positions, and the interaction between three substituents (ME/ME/OH123) located one beside the other (positions numbered 1, 2, 3) for two- or three-substituted benzenic species. Two other additivity groups are also developed using the thermochemical properties of benzenic species taken from the literature: the C/CB/H2/OO and the CB/CO groups. These groups extend the capacities of the group additivity method to deal with substituted benzenic species.  相似文献   

18.
It has been proposed that the (13)C-octanoic acid breath test (OBT) provides a safe, non-radioactive means of measuring gastric emptying. However, deuterated octanoic acid provides a better marker when compared with scintigraphy, as the kinetics are less complex than those of the (13)C label. The appearance of (2)H in saliva is modelled as a two-compartment body water system, using an asymmetric triangular gastric emptying function. This study compared the (2)H-octanoic acid saliva test (OST) with the OBT in measuring altered states of gastric emptying in the nutritional context of diet manipulation. Gastric emptying was measured using the OST and OBT in a three-way crossover study involving 12 healthy male and female subjects (mean BMI = 23.4 kg/m(2), aged 24-57 years). Following an overnight fast, subjects were given an egg meal, labelled with 10 microL/kg body weight (2)H-octanoic acid and 100 microL (13)C-octanoic acid. The meal was nutritionally manipulated to provide a 1 MJ, 2 MJ or 3 MJ meal. Breath and saliva samples were collected at regular intervals for 6 h, with further saliva samples being collected over four subsequent days. (2)H isotopic enrichment in saliva and (13)C isotopic enrichment in breath were analysed using isotope ratio mass spectrometry and the data fitted to the respective gastric emptying models. The half excretion time (T(1/2) (D)), time to maximum emptying rate (T(1) (D)) and time when emptying is complete (T(2) (D)) were calculated from the (2)H saliva test data, and the lag time (T(lag) (C)), half excretion time (T(1/2) (C)), latency phase (T(lat) (C)) and ascension time (T(asc) (C)) were calculated from the (13)C breath test data. Overall, the OBT correlated well with the OST, with a significant relationship between T(1/2) (C) and T(1/2) (D), a significant relationship between T(lat) (C) and T(1) (D) and finally a significant relationship between T(asc) (C) and T(2) (D). Gastric emptying measured using the OST was significantly faster with the 1 MJ meal (DeltaT(1/2) (D) = -0.77 h vs. 2 MJ, p = 0.004). Increases were also seen when the meal size was increased from 2 MJ to 3 MJ (DeltaT(1/2) (D) = +0.44 h vs. 2 MJ), but these were not significant. These trends were mirrored in the OBT data, with significant differences between 1 MJ and 2 MJ (DeltaT(1/2) (C) = -0.63 h vs. 2 MJ, p = 0.013) and non-significant increases with the larger 3 MJ meal (DeltaT(1/2) (C) = +0.10 h vs. 2 MJ). Total meal calorie content was shown to have an effect on gastric emptying using both the OBT and the OST. The deuterium method allows the direct calculation of the gastric emptying function and could be used as an alternative to gamma scintigraphy, allowing further validation of the (13)C-octanoic acid breath test.  相似文献   

19.
Structures, enthalpy (Δ(f)H°(298)), entropy (S°(T)), and heat capacity (C(p)(T)) are determined for a series of nitrocarbonyls, nitroolefins, corresponding nitrites, and their carbon centered radicals using the density functional B3LYP and composite CBS-QB3 calculations. Enthalpies of formation (Δ(f)H°(298)) are determined at the B3LYP/6-31G(d,p), B3LYP/6-31+G(2d,2p), and composite CBS-QB3 levels using several work reactions for each species. Entropy (S) and heat capacity (C(p)(T)) values from vibration, translational, and external rotational contributions are calculated using the rigid-rotor-harmonic-oscillator approximation based on the vibration frequencies and structures obtained from the density functional studies. Contribution to Δ(f)H(T), S, and C(p)(T) from the analysis on the internal rotors is included. Recommended values for enthalpies of formation of the most stable conformers of nitroacetone cc(═o)cno2, acetonitrite cc(═o)ono, nitroacetate cc(═o)no2, and acetyl nitrite cc(═o)ono are -51.6 kcal mol(-1), -51.3 kcal mol(-1), -45.4 kcal mol(-1), and -58.2 kcal mol(-1), respectively. The calculated Δ(f)H°(298) for nitroethylene c═cno2 is 7.6 kcal mol(-1) and for vinyl nitrite c═cono is 7.2 kcal mol(-1). We also found an unusual phenomena: an intramolecular transfer reaction (isomerization) with a low barrier (3.6 kcal mol(-1)) in the acetyl nitrite. The NO of the nitrite (R-ONO) in CH(3)C(═O')ONO moves to the C═O' oxygen in a motion of a stretching frequency and then a shift to the carbonyl oxygen (marked as O' for illustration purposes).  相似文献   

20.
This paper reports the first quantitative ab initio prediction of the disproportionation/combination ratio of alkyl+alkyl reactions using CH3+C2H5 as an example. The reaction has been investigated by the modified Gaussian-2 method with variational transition state or Rice-Ramsperger-Kassel-Marcus calculations for several channels producing (1) CH4+CH2CH2, (2) C3H8, (3) CH4CH3CH, (4) H2+CH3CHCH2, (5) H2+CH3CCH3, and (6) C2H6+CH2 by H-abstraction and association/decomposition mechanisms through singlet and triplet potential energy paths. Significantly, the disproportionation reaction (1) producing CH4+C2H4 was found to occur primarily by the lowest energy path via a loose hydrogen-bonding singlet molecular complex, H3CHC2H4, with a 3.5 kcal/mol binding energy and a small decomposition barrier (1.9 kcal/mol), instead of a direct H-abstraction process. Bimolecular reaction rate constants for the formation of the above products have been calculated in the temperature range 300-3000 K. At 1 atm, formation of C3H8 is dominant below 1200 K. Over 1200 K, the disproportionation reaction becomes competitive. The sum of products (3)-(6) accounts for less than 0.3% below 1500 K and it reaches around 1%-4% above 2000 K. The predicted rate constant for the disproportionation reaction with multiple reflections above the complex well, k1=5.04 x T(0.41) exp(429/T) at 200-600 K and k1=1.96 x 10(-20) T(2.45) exp(1470/T) cm3 molecule(-1) s(-1) at 600-3000 K, agrees closely with experimental values. Similarly, the predicted high-pressure rate constants for the combination reaction forming C3H8 and its reverse dissociation reaction in the temperature range 300-3000 K, k2(infinity)=2.41 x 10(-10) T(-0.34) exp(259/T) cm3 molecule(-1) s(-1) and k(-2)(infinity)=8.89 x 10(22) T(-1.67)exp(-46 037/T) s(-1), respectively, are also in good agreement with available experimental data.  相似文献   

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

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