首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 218 毫秒
1.
The isotropic 129Xe nuclear magnetic resonance (NMR) chemical shift (CS) in Xe@C60 dissolved in liquid benzene was calculated by piecewise approximation to faithfully simulate the experimental conditions and to evaluate the role of different physical factors influencing the 129Xe NMR CS. The 129Xe shielding constant was obtained by averaging the 129Xe nuclear magnetic shieldings calculated for snapshots obtained from the molecular dynamics trajectory of the Xe@C60 system embedded in a periodic box of benzene molecules. Relativistic corrections were added at the Breit–Pauli perturbation theory (BPPT) level, included the solvent, and were dynamically averaged. It is demonstrated that the contribution of internal dynamics of the Xe@C60 system represents about 8% of the total nonrelativistic NMR CS, whereas the effects of dynamical solvent add another 8%. The dynamically averaged relativistic effects contribute by 9% to the total calculated 129Xe NMR CS. The final theoretical value of 172.7 ppm corresponds well to the experimental 129Xe CS of 179.2 ppm and lies within the estimated errors of the model. The presented computational protocol serves as a prototype for calculations of 129Xe NMR parameters in different Xe atom guest–host systems. © 2013 Wiley Periodicals, Inc.  相似文献   

2.
We calculate the 129Xe chemical shift in endohedral Xe@C60 with systematic inclusion of the contributing physical effects to model the real experimental conditions. These are relativistic effects, electron correlation, the temperature-dependent dynamics, and solvent effects. The ultimate task is to obtain the right result for the right reason and to develop a physically justified methodological model for calculations and simulations of endohedral Xe fullerenes and other confined Xe systems. We use the smaller Xe...C6H6 model to calibrate density functional theory approaches against accurate correlated wave function methods. Relativistic effects as well as the coupling of relativity and electron correlation are evaluated using the leading-order Breit-Pauli perturbation theory. The dynamic effects are treated in two ways. In the first approximation, quantum dynamics of the Xe atom in a rigid cage takes advantage of the centrosymmetric potential for Xe within the thermally accessible distance range from the center of the cage. This reduces the problem of obtaining the solution of a diatomic rovibrational problem. In the second approach, first-principles classical molecular dynamics on the density functional potential energy hypersurface is used to produce the dynamical trajectory for the whole system, including the dynamic cage. Snapshots from the trajectory are used for calculations of the dynamic contribution to the absorption 129Xe chemical shift. The calculated nonrelativistic Xe shift is found to be highly sensitive to the optimized molecular structure and to the choice of the exchange-correlation functional. Relativistic and dynamic effects are significant and represent each about 10% of the nonrelativistic static shift at the minimum structure. While the role of the Xe dynamics inside of the rigid cage is negligible, the cage dynamics turns out to be responsible for most of the dynamical correction to the 129Xe shift. Solvent effects evaluated with a polarized continuum model are found to be very small.  相似文献   

3.
Among rare gases, xenon features an unusually broad nuclear magnetic resonance (NMR) chemical shift range in its compounds and as a non-bonded Xe atom introduced into different environments. In this work we show that (129)Xe NMR chemical shifts in the recently prepared, matrix-isolated xenon compounds appear in new, so far unexplored (129)Xe chemical shift ranges. State-of-the-art theoretical predictions of NMR chemical shifts in compounds of general formula HXeY (Y = H, F, Cl, Br, I, -CN, -NC, -CCH, -CCCCH, -CCCN, -CCXeH, -OXeH, -OH, -SH) as well as in the recently prepared ClXeCN and ClXeNC species are reported. The bonding situation of Xe in the studied compounds is rather different from the previously characterized cases as Xe appears in the electronic state corresponding to a situation with a low formal oxidation state, between I and II in these compounds. Accordingly, the predicted (129)Xe chemical shifts occur in new NMR ranges for this nucleus: ca. 500-1000 ppm (wrt Xe gas) for HXeY species and ca. 1100-1600 ppm for ClXeCN and ClXeNC. These new ranges fall between those corresponding to the weakly-bonded Xe(0) atom in guest-host systems (δ < 300 ppm) and in the hitherto characterized Xe molecules (δ > 2000 ppm). The importance of relativistic effects is discussed. Relativistic effects only slightly modulate the (129)Xe chemical shift that is obtained already at the nonrelativistic CCSD(T) level. In contrast, spin-orbit-induced shielding effects on the (1)H chemical shifts of the H1 atom directly bonded to the Xe center largely overwhelm the nonrelativistic deshielding effects. This leads to an overall negative (1)H chemical shift in the range between -5 and -25 ppm (wrt CH(4)). Thus, the relativistic effects induced by the heavy Xe atom appear considerably more important for the chemical shift of the neighbouring, light hydrogen atom than that of the Xe nucleus itself. The predicted NMR parameters facilitate an unambiguous experimental identification of these novel compounds.  相似文献   

4.
We calibrate the methodology for the calculation of nuclear magnetic resonance (NMR) properties in novel organo-xenon compounds. The available state-of-the-art quantum-chemical approaches are combined and applied to the HXeCCH molecule as the model system. The studied properties are (129)Xe, (1)H, and (13)C chemical shifts and shielding anisotropies, as well as (131)Xe and (2)H nuclear quadrupole coupling constants. The aim is to obtain, as accurately as currently possible, converged results with respect to the basis set, electron correlation, and relativistic effects, including the coupling of relativity and correlation. This is done, on one hand, by nonrelativistic correlated ab initio calculations up to the CCSD(T) level and, on the other hand, for chemical shifts and shielding anisotropies by the leading-order relativistic Breit-Pauli perturbation theory (BPPT) with correlated ab initio and density-functional theory (DFT) reference states. BPPT at the uncorrelated Hartree-Fock level as well as the corresponding fully relativistic Dirac-Hartree-Fock method are found to be inapplicable due to a dramatic overestimation of relativistic effects, implying the influence of triplet instability in this multiply bonded system. In contrast, the fully relativistic second-order Moller-Plesset perturbation theory method can be applied for the quadrupole coupling, which is a ground-state electric property. The performance of DFT with various exchange-correlation functionals is found to be inadequate for the nonrelativistic shifts and shielding anisotropies as compared to the CCSD(T) results. The relativistic BPPT corrections to these quantities can, however, be reasonably predicted by DFT, due to the improved triplet excitation spectrum as compared to the Hartree-Fock method, as well as error cancellation within the five main BPPT contributions. We establish three computationally feasible models with characteristic error margins for future calculations of larger organo-xenon compounds to guide forthcoming experimental NMR efforts. The predicted (129)Xe chemical shift in HXeCCH is in a novel range for this nucleus, between weakly bonded or solvated atomic xenon and xenon in the hitherto characterized molecules.  相似文献   

5.
Four-component relativistic and nonrelativistic molecular orbital calculations were performed for the covalent paramagnetic complex [Ir(CN)5]3−, employing the self-consistent discrete variational method, in the framework of density functional theory. Relativistic effects on the electronic structure and chemical bonding are discussed by comparison of relativistic and nonrelativistic one-electron energy levels, populations, and bond orders. The influence of relativistic effects on calculated absorption energies of the electronic spectrum is briefly assessed. © 1996 John Wiley & Sons, Inc.  相似文献   

6.
7.
An approach for hyperpolarized 129Xe molecular sensors is explored using paramagnetic relaxation agents that can be deactivated upon chemical or enzymatic reaction with an analyte. Cryptophane encapsulated 129Xe within the vicinity of the paramagnetic center experiences fast relaxation that, through chemical exchange of xenon atoms between cage and solvent pool, causes accelerated hyperpolarized 129Xe signal decay in the dissolved phase. In this proof‐of‐concept work, the relaxivity of Gadolinium III‐DOTA on 129Xe in the solvent was increased eightfold through tethering of the paramagnetic molecule to a cryptophane cage. This potent relaxation agent can be ′turned off′ specifically for 129Xe through chemical reactions that spatially separate the GdIII centre from the attached cryptophane cage. Unlike 129Xe chemical shift based sensors, the new concept does not require high spectral resolution and may lead to a new generation of responsive contrast agents for molecular MRI.  相似文献   

8.
9.
To elucidate the physical origin of relativistic changes of molecular properties, exact theorems, perturbation theory, and Hartree-Fock-Slater-Pauli calculations are exploited. The relativistic molecular virial theorem offers insight into the relativistic and nonrelativistic, kinetic, and potential energy contributions to the bond energy. In general, there exist two contributions to the relativistic correction of a molecular property: the relativistic change at the nonrelativistic equilibrium geometry and the change of the nonrelativistic property due to the relativistic change of the equilibrium geometry. Sometimes the first and sometimes the second contribution is the dominant one. Accurate numerical results for H+2-like systems are obtained using direct relativistic double perturbation theory. In some cases, near-degenerate perturbation theory is mandatory. Relativistic changes of chemical bond energies are often proportional to the density change in the K-shell when the bond is formed. Relativistic corrections to many properties (and also to the 1s2-correlation energy) are often proportional to Z2α2. © 1996 John Wiley & Sons, Inc.  相似文献   

10.
Ab initio accurate all-electron relativistic molecular orbital Dirac–Fock self-consistent field calculations are reported for the linear symmetric XeF2 molecule at various internuclear distances with our recently developed relativistic universal Gaussian basis set. The nonrelativistic limit Hartree–Fock calculations were also performed for XeF2 at various internuclear distances. The relativistic correction to the electronic energy of XeF2 was calculated as ~ ?215 hartrees (?5850 eV) by using the Dirac–Fock method. The dominant magnetic part of the Breit interaction correction to the nonrelativistic interelectron Coulomb repulsion was included in our calculations by both the Dirac–Fock–Breit self-consistent field and perturbation methods. The calculated Breit correction is ~6.5 hartrees (177 eV) for XeF2. The relativistic Dirac–Fock as well as the nonrelativistic HF wave functions predict XeF2 to be unbound, due to neglect of electron correlation effects. These effects were incorporated for XeF2 by using various ab initio post Hartree–Fock methods. The calculated dissociation energy obtained using the MP 2(full) method with our extensive basis set of 313 primitive Gaussians that included d and f polarization functions on Xe and F is 2.77 eV, whereas the experimental dissociation energy is 2.78 eV. The calculated correlation energy is ~ ?2 hartrees (?54 eV) at the predicted internuclear distance of 1.986 Å, which is in excellent agreement with the experimental Xe—F distance of 1.979 Å in XeF2. In summary, electron correlation effects must be included in accurate ab initio calculations since it has been shown here that their inclusion is crucial for obtaining theoretical dissociation energy (De) close to experimental value for XeF2. Furthermore, relativistic effects have been shown to make an extremely significant contribution to the total energy and orbital binding energies of XeF2. © 1995 John Wiley & Sons, Inc.  相似文献   

11.
Chemical shift tensors obtained from solid-state NMR spectroscopy are very sensitive reporters of structure and dynamics in proteins. While accurate 13C and 15N chemical shift tensors are accessible by magic angle spinning (MAS) NMR, their quantum mechanical calculations remain challenging, particularly for 15N atoms. Here we compare experimentally determined backbone 13Cα and 15NH chemical shift tensors by MAS NMR with hybrid quantum mechanics/molecular mechanics/molecular dynamics (MD-QM/MM) calculations for the carbohydrate-binding domain of galectin-3. Excellent agreement between experimental and computed 15NH chemical shift anisotropy values was obtained using the Amber ff15ipq force field when solvent dynamics was taken into account in the calculation. Our results establish important benchmark conditions for improving the accuracy of chemical shift calculations in proteins and may aid in the validation of protein structure models derived by MAS NMR.  相似文献   

12.
Advancement of hyperpolarized 129Xe MRI technology toward clinical settings demonstrates the considerable interest in this modality for diagnostic imaging. The number of contrast agents, termed biosensors, for 129Xe MRI that respond to specific biological targets, has grown and diversified. Directly functionalized xenon-carrying macrocycles, such as the large family of cryptophane-based biosensors, are good for localization-based imaging and provide contrast before and after binding events occur. Noncovalently functionalized constructs, such as cucurbituril- and cyclodextrin-based biosensors, benefit from commercial availability and optimal exchange dynamics for CEST imaging. In this work, we report the first directly functionalized cucurbituril used as a xenon biosensor. Biotinylated cucurbit[7]uril (btCB7) gives rise to a 129Xe hyperCEST response at the unusual shift of δ=28 ppm when bound to its protein target with substantial CEST contrast. We posit that the observed chemical shift is due to the deformation of btCB7 upon binding to avidin, caused by proximity to the protein surface. Conformational searches and molecular dynamics (MD) simulations support this hypothesis. This construct combines the strengths of both families of biosensors, enables a multitude of biological targets through avidin conjugation, and demonstrates the advantages of functionalized cucurbituril-based biosensors.  相似文献   

13.
Ab initio all-electron fully relativistic Dirac–Fock self-consistent field and Dirac–Fock–Breit calculations are reported for the XeF4 molecule at various internuclear distances assuming the experimental D4h geometry with our recently developed relativistic universal Gaussian basis set. The nonrelativistic limit Hartree–Fock calculations were also performed for XeF4 at various internuclear distances. The calculated relativistic correction to the total energy of molecule at the Dirac–Fock level is ~ ?5856 eV, whereas the magnetic part of the Breit correction to the electron-electron interaction is calculated as ~ 177 eV. The electron correlation effects were included in the nonrelativistic Hartree–Fock calculations using the second-order Møller-Plesset (MP 2) theory, and the calculated correlation energy for XeF4 is ?71 eV. The basis-set superposition error (BSSE ) was estimated by using the counterpoise method for Xe and F. The inclusion of both the relativistic and electron correlation effects in the calculated total energies of F, Xe, and XeF4 predicts the Xe—F bond length and dissociation energy of XeF4 as 1.952 Å and 5.59 eV, respectively, which are in excellent agreement with the experimental values of 1.953 Å and 5.69 eV, respectively, for XeF4. The contribution of the electron correlation and relativistic effects to the dissociation energy of XeF4 is 8.11 and 0.05 eV, respectively. The Breit interaction, however, contributes only 0.02 eV to the dissociation energy of XeF4. Electron correlation is most significant for the prediction of an accurate value of dissociation energy, whereas relativistic effects are very important for the prediction of spin-orbital splitting as well as the energies of the orbitals, especially the inner orbitals of XeF4. © 1995 John Wiley & Sons, Inc.  相似文献   

14.
The principal relativistic heavy-atom effects on the nuclear magnetic resonance (NMR) shielding tensor of the heavy atom itself (HAHA effects) are calculated using ab initio methods at the level of the Breit-Pauli Hamiltonian. This is the first systematic study of the main HAHA effects on nuclear shielding and chemical shift by perturbational relativistic approach. The dependence of the HAHA effects on the chemical environment of the heavy atom is investigated for the closed-shell X(2+), X(4+), XH(2), and XH(3) (-) (X=Si-Pb) as well as X(3+), XH(3), and XF(3) (X=P-Bi) systems. Fully relativistic Dirac-Hartree-Fock calculations are carried out for comparison. It is necessary in the Breit-Pauli approach to include the second-order magnetic-field-dependent spin-orbit (SO) shielding contribution as it is the larger SO term in XH(3) (-), XH(3), and XF(3), and is equally large in XH(2) as the conventional, third-order field-independent spin-orbit contribution. Considering the chemical shift, the third-order SO mechanism contributes two-thirds of the difference of approximately 1500 ppm between BiH(3) and BiF(3). The second-order SO mechanism and the numerically largest relativistic effect, which arises from the cross-term contribution of the Fermi contact hyperfine interaction and the relativistically modified spin-Zeeman interaction (FC/SZ-KE), are isotropic and practically independent of electron correlation effects as well as the chemical environment of the heavy atom. The third-order SO terms depend on these factors and contribute both to heavy-atom shielding anisotropy and NMR chemical shifts. While a qualitative picture of heavy-atom chemical shifts is already obtained at the nonrelativistic level of theory, reliable shifts may be expected after including the third-order SO contributions only, especially when calculations are carried out at correlated level. The FC/SZ-KE contribution to shielding is almost completely produced in the s orbitals of the heavy atom, with values diminishing with the principal quantum number. The relative contributions converge to universal fractions for the core and subvalence ns shells. The valence shell contribution is negligible, which explains the HAHA characteristics of the FC/SZ-KE term. Although the nonrelativistic theory gives correct chemical shift trends in present systems, the third-order SO-I terms are necessary for more reliable predictions. All of the presently considered relativistic corrections provide significant HAHA contributions to absolute shielding in heavy atoms.  相似文献   

15.
Xenon fluoride radicals were generated by solid-state chemical reactions of mobile fluorine atoms with xenon atoms trapped in Ar matrix. Highly resolved electron spin resonance spectra of XeF* were obtained in the temperature range of 5-25 K and the anisotropic hyperfine parameters were determined for magnetic nuclei 19F, 129Xe, and 131Xe using naturally occurring and isotopically enriched xenon. Signs of parallel and perpendicular hyperfine components were established from analysis of temperature changes in the spectra and from numerical solutions of the spin Hamiltonian for two nonequivalent magnetic nuclei. Thus, the complete set of components of hyperfine- and g-factor tensors of XeF* were obtained: 19F (Aiso=435, Adip=1249 MHz) and 129Xe (Aiso=-1340, Adip=-485 MHz); g(parallel)=1.9822 and g(perpendicular)=2.0570. Comparison of the measured hyperfine parameters with those predicted by density-functional theory (DFT) calculations indicates, that relativistic DFT gives true electron spin distribution in the 2Sigma+ ground-state, whereas nonrelativistic theory underestimates dramatically the electron-nuclear contact Fermi interaction (Aiso) on the Xe atom. Analysis of the obtained magnetic-dipole interaction constants (Adip) shows that fluorine 2p and xenon 5p atomic orbitals make a major contribution to the spin density distribution in XeF*. Both relativistic and nonrelativistic calculations give close magnetic-dipole interaction constants, which are in agreement with the measured values. The other relativistic feature is considerable anisotropy of g-tensor, which results from spin-orbit interaction. The orbital contribution appears due to mixing of the ionic 2Pi states with the 2Sigma+ ground state, and the spin-orbit interaction plays a significant role in the chemical bonding of XeF*.  相似文献   

16.
《Microporous Materials》1994,2(2):127-136
The adsorption isotherms and 129Xe nuclear magnetic resonance (NMR) chemical shifts of xenon and the adsorption isotherms of carbon monoxide of Cu(II)- and Cu(I)-exchanged zeolites NaY were measured. The former zeolites of 53, 75, and 95% exchange degrees were investigated after various pretreatment steps comprising dehydration, oxidation and reduction with CO at 420°C as well as long-term CO reduction at 470°C. The Cu(I)Y zeolite of 70% exchange degree was prepared via a solid-state exchange procedure with CuCl and subjected to dehydration at 420°C. In all cases, except the dehydrated zeolites, almost linear xenon adsorption isotherms and linear 129Xe NMR chemical shift versus xenon concentration curves running parallel to each other are obtained. In contrast, the chemical shift curves for the dehydrated zeolites are non-linear at low xenon concentrations turning towards negative chemical shift values at very low pressures. The whole body of the experimental xenon data can be explained quantitatively with a unifying approach on the basis of a site adsorption model where the sites are (i) two types of cuprous ions of much different adsorption strength and 129Xe chemical shift, (ii) Na+ cations, (iii) Lewis acid sites generated through autoreduction and reduction of Cu2+ by CO, and (iv) framework sites free of cations. These five types of sites are each characterized by Langmuir adsorption isotherm constants and local 129Xe NMR chemical shifts. The adsorption site concentrations in the various zeolites are evaluated. The supercage Cu(I) concentration values are in nice agreement with the results deduced from the CO adsorption isotherm measurements.  相似文献   

17.
We present perturbational ab initio calculations of the nuclear-spin-dependent relativistic corrections to the nuclear magnetic resonance shielding tensors that constitute, together with the other relativistic terms reported by us earlier, the full leading-order perturbational set of results for the one-electron relativistic contributions to this observable, based on the (Breit-)Pauli Hamiltonian. These contributions are considered for the H(2)X (X = O,S,Se,Te,Po) and HX (X = F,Cl,Br,I,At) molecules, as well as the noble gas (Ne, Ar, Kr, Xe, Rn) atoms. The corrections are evaluated using the relativistic and magnetic operators as perturbations on an equal footing, calculated using analytical linear and quadratic response theory applied on top of a nonrelativistic reference state provided by self-consistent field calculations. The (1)H and heavy-atom nuclear magnetic shielding tensors are compared with four component, nearly basis-set-limit Dirac-Hartree-Fock calculations that include positronic excitations, as well as available literature data. Besides the easy interpretability of the different contributions in terms of familiar nonrelativistic concepts, the accuracy of the present perturbational scheme is striking for the isotropic part of the shielding tensor, for systems including elements up to Xe.  相似文献   

18.
Quantum chemical calculations of the nuclear shielding tensor, the nuclear quadrupole coupling tensor, and the spin-rotation tensor are reported for the Xe dimer using ab initio quantum chemical methods. The binary chemical shift delta, the anisotropy of the shielding tensor Delta sigma, the nuclear quadrupole coupling tensor component along the internuclear axis chi( parallel ), and the spin-rotation constant C( perpendicular ) are presented as a function of internuclear distance. The basis set superposition error is approximately corrected for by using the counterpoise correction (CP) method. Electron correlation effects are systematically studied via the Hartree-Fock, complete active space self-consistent field, second-order M?ller-Plesset many-body perturbation, and coupled-cluster singles and doubles (CCSD) theories, the last one without and with noniterative triples, at the nonrelativistic all-electron level. We also report a high-quality theoretical interatomic potential for the Xe dimer, gained using the relativistic effective potential/core polarization potential scheme. These calculations used valence basis set of cc-pVQZ quality supplemented with a set of midbond functions. The second virial coefficient of Xe nuclear shielding, which is probably the experimentally best-characterized intermolecular interaction effect in nuclear magnetic resonance spectroscopy, is computed as a function of temperature, and compared to experiment and earlier theoretical results. The best results for the second virial coefficient, obtained using the CCSD(CP) binary chemical shift curve and either our best theoretical potential or the empirical potentials from the literature, are in good agreement with experiment. Zero-point vibrational corrections of delta, Delta sigma, chi (parallel), and C (perpendicular) in the nu=0, J=0 rovibrational ground state of the xenon dimer are also reported.  相似文献   

19.
The 129Xe NMR line shapes of xenon adsorbed in the nanochannels of the (+/-)-[Co(en)3]Cl3 ionic crystal have been calculated by grand canonical Monte Carlo (GCMC) simulations. The results of our GCMC simulations illustrate their utility in predicting 129Xe NMR chemical shifts in systems containing a transition metal. In particular, the nanochannels of (+/-)-[Co(en)3]Cl3 provide a simple, yet interesting, model system that serves as a building block toward understanding xenon chemical shifts in more complex porous materials containing transition metals. Using only the Xe-C and Xe-H potentials and shielding response functions derived from the Xe@CH4 van der Waals complex to model the interior of the channel, the GCMC simulations correctly predict the 129Xe NMR line shapes observed experimentally (Ueda, T.; Eguchi, T.; Nakamura, N.; Wasylishen, R. E. J. Phys. Chem. B 2003, 107, 180-185). At low xenon loading, the simulated 129Xe NMR line shape is axially symmetric with chemical-shift tensor components delta(parallel) = 379 ppm and delta(perpendicular) = 274 ppm. Although the simulated isotropic chemical shift, delta(iso) = 309 ppm, is overestimated, the anisotropy of the chemical-shift tensor is correctly predicted. The simulations provide an explanation for the observed trend in the 129Xe NMR line shapes as a function of the overhead xenon pressure: delta(perpendicular) increased from 274 to 292 ppm, while delta(parallel) changed by only 3 ppm over the entire xenon loading range. The overestimation of the isotropic chemical shifts is explained based upon the results of quantum mechanical 129Xe shielding calculations of xenon interacting with an isolated (+/-)-[Co(en)3]Cl3 molecule. The xenon chemical shift is shown to be reduced by about 12% going from the Xe@[Co(en)3]Cl3 van der Waals complex to the Xe@C2H6 fragment.  相似文献   

20.
The miscibility of two-component polymer blends has been investigated using xenon-129 (129Xe) nuclear magnetic resonance (NMR) to probe the phase morphology. The chemical shift of 129Xe dissolved in a given polymer is unique, thus heterogeneous blends with large domain sizes exhibit two 129Xe NMR lines. When a single resonance is obtained, the data are consistent with miscibility, yielding an upper bound on the domain size. The temperature dependence of the relative solubilities and chemical shifts of 129Xe dissolved in the pure components may allow a determination of the phase morphology in blends exhibiting a single resonance. The method is used to demonstrate that polychloroprene and 25% epoxidized 1,4-polyisoprene form a miscible blend.  相似文献   

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

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