首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 312 毫秒
1.
Reference completely ab initio 6–3G and nonempirical 3G/MODPOT (ab initio effective core model potential) LCAO -MO -SCF calculations (using the same valence atomic orbital basis) were performed for a series of boron hydrides (B4H10, B5H9, B5H11, and B6H10) and a test 3G/MODPOT + VRDDO (variable retention of diatomic differential overlap for charge conserving integral prescreening) calculation were also performed for B5H9, B6H10, and B10H14. The agreement between the ab initio 6–3G and the corresponding 3G/MODPOT calculations was excellent for valence orbital energies, gross atomic populations, and dipole moments. The results also compared favorably to previous ab initio minimum STO basis results of Lipscomb and coworkers. The 3G/MODPOT + VRDDO calculations verified that for such spatially compact molecules (such as boron hydrides, which are fragments of polyhedra), the VRDDO procedure does not result in a noticeable savings in computer time for molecules of the size and shape of B5H9 and B6H10, in contrast to the savings previously realized for organic molecules of comparable atomic size. However, the agreement in calculational results between the 3G/MODPOT and the 3G/MODPOT +VRDDO results was still as extremely close as it had been for the organic molecules. 3G/MODPOT calculations were also carried out for B8H12, B9H15, B10H14, B10H14?2, 1,2-C2B4H6, and 1,6-C2B4H6 and the results compared to the previous minimum STO basis results. For B10H14, the 3G/MODPOT +VRDDO method led to savings in computer time of 28% over the 3G/MODPOT method itself. The agreement of the 3G/MODPOT results with available experimental photoelectron spectral data for B5H9 and 1,6-C2B4H6 was as good as that of the previous ab initio minimum STO basis calculations.  相似文献   

2.
Ab initio MODPOT/VRDDO/MERGE calculations were carried out on carcinogenic 3-methylcholanthrene (3-MCA) and its metabolites. The results for 3-MCA were compared to our earlier similar calculations for carcinogenic benzo(a)pyrene (BP). Both compounds 3-MCA and BP are carcinogenic and are metabolically activated by similar mechanisms but in different positions. Both the calculated wave functions for 3-MCA and BP and the electrostatic molecular potential contour maps generated from these wave functions correctly reflect the similarity of mechanisms of metabolic activation and the differences in position. Our calculated results both for BP and for 3-MCA reflect accurately their experimentally observed behavior. Thus this combination of theoretical techniques can be used with confidence to describe the behavior of the polycyclic aromatic hydrocarbons (PAH's) and their metabolites. The ab initio MODPOT/VRDDO method incorporates two very desirable options into our fast ab initio Gaussian programs: MODPOT –ab initio effective core model potentials—and a charge-conserving integral prescreening approximation which we named VRDDO (variable retention of diatomic differential overlap). For orbital energies and population analysis the MODPOT/VRDDO results agree to essentially three decimal places with completely ab initio calculations using the same valence atomic basis set. For this series of very closely related congeners our recent MERGE technique which allows reuse of integrals from a common skeletal fragment was used. The ab initio MODPOT/VRDDO/MERGE calculations were carried out for 3-MCA, 3-MCA oxides, 3-MCA dihydrodiols, and 3-MCA dihydrodiolepoxides. The metabolites investigated were 3-MCA 9,10-oxide; 3-MCA 7,8-oxide; 3-MCA 9,10-dihydrodiol [trans(axial, axial); trans(equatorial, equatorial); cis(axial, equatorial); cis(equatorial, axial)]; and 3-MCA 9,10-dihydrodiol–7,8-epoxide [for both conformations A and B of the dihydrodiol and for all stereoisomers of the dihydrodiolepoxides relative to below and above the plane: ααα, and ααβ αβα αββ βαα βαβ ββα and βββ (most stable)]. Calculations were also carried out for opening of the C7? O? C8 epoxide ring both towards C7 and C8 for the most stable isomer Aβββ (above the ring). Opening the epoxide ring between C7 and O leads to a more stable intermediate than opening the epoxide ring between C8 and O. Again, however, as with opening the epoxide ring in BP 7,8-dihydrodiol–9,10-epoxide there is no buildup of positive charge on C7 in the 3-MCA metabolites as postulated by some cancer researchers, but rather the C7 becomes slightly more negative. Nor is there a buildup of negative charge on the O atom, but rather it becomes slightly more positive. As the epoxide ring is opened further than 90° for the O? C7? C8 or O? C8? C7 angles, there appears to be a possible mixing of configurations that is being investigated further.  相似文献   

3.
Ab initio MODPOT /VRDDO calculations have been carried out on carcinogenic benzo(a)pyrene and its metabolites. The MODPOT /VRDDO method incorporates two very desirable options into our fast ab initio Gaussian programs: MODPOT —ab initio effective core model potentials—and a charge-conserving integral prescreening approximation which we named VRDDO (variable retention of diatomic differential overlap). For orbital energies and population analyses the MODPOT /VRDDO results agree to essentially three decimal places with completely ab initio calculations using the same valence atomic basis set. For this series of very closely related congeners a new MERGE technique was implemented that allows reuse of integrals of a common skeletal fragment. Since our program computes integrals efficiently by blocks, reusing information common to the block, it was more difficult to implement a MERGE technique than for integral programs which calculate the integrals one-byone. The MODPOT /VRDDO calculations were performed for benzo(a)pyrene (BP), BP oxides, BP dihydrodiols, and BP dihydrodiol epoxides. The metabolites investigated were BP-7,8-oxide, BP-4,5-oxide, BP-7,8-dihydrodiol [cis(e, a), cis(a, e), trans(e, e), and trans(a, a)], and BP-7,8-dihydrodiol-9,10-epoxide [β,β,β (the most stable), β,β,α; α,α,β, and α,α,α all derived from cis-BP-7,8-dihydrodiol and β,α,β; α,β,β and α,β,β derived from trans-BP-7,8-dihydrodiol]. Several different conformations were calculated for each of the BP dihydrodiols and BP dihydrodiol epoxides. Calculations were carried out for the opening of the C9—O—C10 epoxide ring both toward C9 and C10 for the, most stable β,β,β isomer of BP-7,8-dihydrodiol-9,10-epoxide. Opening the epoxide ring between C10 and O leads to a more stable intermediate than opening the epoxide ring between C9 and C10. However, there is no buildup of positive charge in C10 as has been postulated by some cancer researchers, but rather the C10 becomes slightly more negative. Nor is there a buildup of negative charge on the O atom. rather it becomes slightly less negative. As the epoxide ring is opened further than 90° for the O—C9—C10 or O—C10—C9 angles, there appears to be a possible mixing of configurations that is being investigated further.  相似文献   

4.
Ab initio MODPOT /VRDDO /MERGE calculations have been carried out for all the different position isomers of nitrocubane from mononitrocubane through octanitrocubane for a perfect symmetrical cubic cubane skeleton and for mononitrocubane through septanitrocubane for the almost cubic experimentally determined cubane skeleton. These calculations were carried out with our own rapid efficient ab initio programs which also incorporate a number of desirable computational strategies for calculations on large molecules. The skeletal total overlap population of the cubane skeleton (a theoretical index we showed years ago to be sensitive and predictive of stability of energetic molecular frameworks) indicates that successive nitration seems to increase the stability of the cubane skeleton. Successive nitration also seems to increase the total overlap population of the C? NO2 bond. There are subtle differences depending on the exact positional isomer for a constant number of nitro groups—but the overall trend is definite. We have also generated electrostatic molecular potential contour (EMPC ) maps around these nitrocubanes. These maps are indicative of preferred positions of electrophilic and nucleophilic attack as a function of the number of nitro groups or their positions. These EMPC maps can also indicate, to a first approximation, a limit on how close these molecules may be able to approach each other in a crystal.  相似文献   

5.
The Structure of an unusual Tetramere of Lithium Phenoxide: [C6H5OLi · C4H8O]4 · C6H5OH Single crystals of lithium phenoxide have been obtained from THF. In the structure (P 21/n, Z = 4, a = 11.69 Å, b = 21.15 Å, c = 18.55 Å, β = 91.11°) four lithium atoms and four phenoxide oxygen atoms are cubically arranged. Further, each lithium atom coordinates the oxygen atom of a tetrahydrofuran molecule. The ideal cubeform structure is disturbed by one phenol molecule which is coordinated in addition to four phenoxide and four THF molecules. Hence, one edge of the cube (Li4? O4) is substituted by the coordination of the phenol oxygen atom O5 with Li4 and hydrogen bonding between O4 and the hydroxy group of phenol. Van der Waals forces are the only interaction between these complexes.  相似文献   

6.
The [C6H6O]+˙ ion from C6H5OC4H9 is structurally similar to the molecular ion from phenol as demonstrated by 13C labeling.  相似文献   

7.
Quantum chemical ab initio MODPOT /VRDDO calculations have been carried out on the following aminonitrobenzenes for which crystal structures had been determined experimentally: 4-nitroaniline; N,N-dimethyl-p-nitroaniline; 2,4,6-trinitroaniline; 1,3-diamino-2,4,6-trinitrobenzene (DATB—Form I); 1,3,5-triamino-2,4,6-trinitrobenzene (TATB); 2,3,4,6-tetranitroaniline; N-methyl-N,2,4,6-tetranitroaniline (Tetryl); and N-(β,β,β-trifluoroethyl)-N,2,4,6-tetranitroaniline. These quantum chemical calculations were performed on the molecules in their conformations as found in their crystal structures. The calculations were carried out with our own ab initio programs which also incorporate as options several desirable features for calculations on large molecules: ab initio effective core model potentials (MODPOT) which enable calculations of valence electrons only explicitly, yet accurately, and a charge conserving integral prescreening evaluation (which we named VRDDO-variable retention of diatomic differential overlap) especially effective for spatially extended molecules. Aminonitrobenzenes are especially interesting since there are inherent intramolecular ring distortions and deviations from planarity and intramolecular hydrogen bonds as well as intermolecular hydrogen bonds causing further deviations from planarity. The theoretical indices resulting from the quantum chemical calculations are relevant to a number of properties and behavioral characteristics of these molecules, both intramolecular and intermolecular. The charges on the atoms [from the gross atomic populations (GAP 's)] are needed for calculation of the atomic multipole–atomic multipole electrostatic contributions (a dominant factor) to the intermolecular interaction energies. These electrostatic interaction energies are part of the input necessary for calculations on the crystal packing and densities of these molecules. These GAP 's are also of value in interpreting the experimental photoelectron and ESCA spectra of these molecules. The total overlap populations (TOP 's) between atoms are related to the inherent bond strengths and can serve as a quantitative replacement for the old empirical bond length-bond order-bond energy relationship still used by explosives chemists to identify the “target bonds” (the weakest bonds). The TOP 's are of considerable value in predicting and tracing initiation and subsequent steps of explosive phenomena. The molecular orbital energies of the lowest unoccupied orbitals are of interest since nitroexplosives have been implicated in testicular toxicity and the initial metabolic activation appears to proceed through a one-electron reduction of the nitroexplosive.  相似文献   

8.
Fragmentation of the pyridine ring upon K‐shell excitation/ionization has been studied with gaseous 2‐, 3‐ and 4‐methylpyridine by the electron‐impact method. Ab initio molecular orbital (MO) calculations were also carried out to explore electronic states correlating with specific fragments. Some specific fragmentation channels were identified from the ionic fragments enhanced characteristically at the N 1s edge. Yields of the C2HN+ and C5H5+/C5H6+ ions show that the fission of the N? C2 and C4? C5/C5? C6 bonds of the ring is likely to occur after the N 1s excitation and ionization. Ab initio MO calculations for the 2‐methylpyridine molecule indicate that the dissociation channels to produce these ions are only accessible through the excited states of the parent molecular dication, which can be formed by Auger decays after the N 1s ionization. Fragment ions via hydrogen rearrangement are produced as well, but the rearrangement is not a phenomenon specific to the K‐shell excitation/ionization. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

9.
We report the synthesis of [n]manganoarenophanes (n=1, 2) featuring boron, silicon, germanium, and tin as ansa‐bridging elements. Their preparation was achieved by salt‐elimination reactions of the dilithiated precursor [Mn(η5‐C5H4Li)(η6‐C6H5Li)]?pmdta (pmdta=N,N,N′,N′,N′′‐pentamethyldiethylenetriamine) with corresponding element dichlorides. Besides characterization by multinuclear NMR spectroscopy and elemental analysis, the identity of two single‐atom‐bridged derivatives, [Mn(η5‐C5H4)(η6‐C6H5)SntBu2] and [Mn(η5‐C5H4)(η6‐C6H5)SiPh2], could also be determined by X‐ray structural analysis. We investigated for the first time the reactivity of these ansa‐cyclopentadienyl–benzene manganese compounds. The reaction of the distannyl‐bridged complex [Mn(η5‐C5H4)(η6‐C6H5)Sn2tBu4] with elemental sulfur was shown to proceed through the expected oxidative addition of the Sn?Sn bond to give a triatomic ansa‐bridge. The investigation of the ring‐opening polymerization (ROP) capability of [Mn(η5‐C5H4)(η6‐C6H5)SntBu2] with [Pt(PEt3)3] showed that an unexpected, unselective insertion into the Cipso?Sn bonds of [Mn(η5‐C5H4)(η6‐C6H5)SntBu2] had occurred.  相似文献   

10.
The electrostatic molecular potential contour maps were calculated for carcinogenic 3-methylcholanthrene (3-MCA) and a number of its metabolites {3-MCA 7,8-oxide and 3-MCA 9,10-oxide; 3-MCA 7,8-dihydrodiols[several stereoisomers: A trans(equatorial, equatorial) and A cis(equatorial, axial)]; 3-MCA 9,10-dihydrodiol–7,8-epoxide A βββ and 3-MCA 9,10-dihydrodiol–7,8-epoxide A αβα}. The maps were generated from our ab initio MODPOT/VRDDO/MERGE wave functions calculated for these species. The results of these maps for 3-MCA [similarly to our results for the maps we generated for benzo(a)pyrene (BP)] show that these electrostatic molecular potential contour maps can be used to indicate favored positions of attack for electrophilic species, such as “electrophilic” oxygen to form an epoxide as well as for positive ion attack. The 3-MCA maps indicate the favored site for attack and the pathways. The maps around 3-MCA 9,10-oxide and around 3-MCA 9,10-dihydrodiol-7,8-epoxide indicate the directional preferences for proton assisted epoxide ring opening. The maps around the 3-MCA dihydrodiols indicate that while for certain stereoisomers the “electrophilic” oxygen will prefer to attack from below, for other isomers it will prefer to attack from above. This gives great insight into the stereochemical preference for formation of different 3-MCA 9,10-dihydrodiol–7,8-epoxides.  相似文献   

11.
Preparation, 11B, 13C, 1H NMR and Vibrational Spectra of Monoethoxyhydro-closo-dodecaborate(2–), and the Crystal Structure of [(C5H5N)2CH2][B12H11(OC2H5)] By treatment of Na2[B12H12] with dry HF in ethanol Na2[B12H11(OC2H5)] is formed which has been separated by ion exchange chromatography on diethylaminoethyl(DEAE) cellulose from the starting compound and by-products. The X-ray structure determination of [(C5H5N)2CH2][B12H11(OC2H5)] (monoclinic, space group P21/m, a = 9.1906(3), b = 12.6612(8), c = 9.3640(12) Å, β = 112.947(6)°, Z = 2) reveals the complete ordering of the anion sublattice. The 11B nmr spectrum exhibits the characteristic feature (1:5:5:1) of a mono substituted B12 cage with a strong down-field shift of ipso-B at +6.5 ppm. In the 13C nmr spectrum a triplet at 67.9 ppm of the methylene group and a quartet at 19.5 ppm of the methyl group is observed. Correspondingly, the 1H nmr spectrum shows two multiplets at 3.7 and 1.3 as expected for an ethoxy substituent, and a multiplet at 2.1 ppm due to the protons of the boron cluster. The i.r. and Raman spectra exhibit strong CH stretching vibrations between 2 963 and 2 863 cm?1, and in the i.r. spectrum the CO and BO stretching frequencies of the B? O? C bridge are observed at 1 175 and 1 140 cm?1.  相似文献   

12.
We report the first example of aryl hydrogen scrambling occurring in a molecular anion prior to or accompanying fragmentation, i.e. for the reaction [M]?· → [M ? H2NO2.]? from o-NO2? C6H4? X? C6H5 (X = O or S). Proximity effects occur in these spectra when X = CO, NH, O, or S, and certain of these have been substantiated by 2H and 18O labelling.  相似文献   

13.
In the complex salt [η6‐1‐chloro‐2‐(pyrrolidin‐1‐yl)benzene](η5‐cyclopentadienyl)iron(II) hexafluoridophosphate, [Fe(C5H5)(C10H12ClN)]PF6, (I), the complexed cyclopentadienyl and benzene rings are almost parallel, with a dihedral angle between their planes of 2.3 (3)°. In a related complex salt, (η5‐cyclopentadienyl){2‐[η6‐2‐(pyrrolidin‐1‐yl)phenyl]phenol}iron(II) hexafluoridophosphate, [Fe(C5H5)(C16H17NO)]PF6, (II), the analogous angle is 5.4 (1)°. In both complexes, the aromatic C atom bound to the pyrrolidine N atom is located out of the plane defined by the remaining five ring C atoms. The dihedral angles between the plane of these five ring atoms and a plane defined by the N‐bound aromatic C atom and two neighboring C atoms are 9.7 (8) and 5.6 (2)° for (I) and (II), respectively.  相似文献   

14.
The mol­ecules of the title compound, C7H5BrO2, form zigzag chains running along the b axis and are stacked in layers perpendicular to the a axis. Intermolecular bonding occurs through hydrogen bonds linking the hydroxyl and carbonyl groups, with an O?O distance of 2.804 (4) Å. The Br atom deviates significantly from the plane of the ring and the aldehyde group is twisted by 7.1 (5)° around the Csp2—Caryl bond. The geometry of the mol­ecule in the crystal is compared to that given by ab initio quantum mechanical calculations for the isolated mol­ecule, using a molecular orbital Hartree–Fock method and density functional theory.  相似文献   

15.
We have explored two areas of approximately rigorous calculations for computing nonempirical wave functions for heavy and/or large molecules orders of magnitude faster than with conventional ab-initio methods but with the same chemical accuracy. First, we have developed and used a series of programs (starting from our new fast sets of ab-initio Gaussian SCF and SCF -CI programs) incorporating ab-initio effective core model potentials (MOD -POT ) which allow one to treat only the valence electrons explicitly, plus a charge conserving integral prescreening, which cuts down significantly on the number of integrals that have to be calculated, stored, or processed for a large molecule. We have named this latter procedure VRDDO (variable retention of diatomic differential overlap). With these MODPOT and MODPOT /VRDDO methods we have explored a variety of small, medium, and large systems ranging from electron affinities of atoms through to molecules of biological interest and large boron hydrides. The results compared to ab-initio SCF or SCF /CI calcuations are very good, usually within 0.001 to 0.002 a.u. for orbital energies and gross atomic populations (GAPS ) and even better along potential energy curves. Secondly, we have explored the use of the MS -Xα method for less conventional molecules and properties than those for which it is customarily employed.  相似文献   

16.
The Structures of two Hydrates of Sodium Phenoxide: C6H5ONa · H2O and C6H5ONa · 3 H2O In the monohydrate of sodium phenoxide sodium is coordinated by 4 oxygen atoms having an average distance Na? O of about 2.631 Å being arranged in form of a distorted tetrahedron. The oxygen atoms of water and phenoxid serve as bridging ligands. Hence, the structure can be considered as a network with a general formula [Na[4]O]. Moreover, the oxygen atoms are linked via hydrogen bonding. In the trihydrate of sodium phenoxide sodium is surrounded with 5 oxygen atoms with an average distance of 2.39 Å forming a tetragonal pyramide. The oxygen of the phenoxide, however, does not participate in the coordination of the sodium ion. The coordination polyhedrons are connected by sharing edges and verteces. The resulting layer can be described by the general formula [Na[5]O2[2]O[2]O[1]]. Via hydrogen bonding the phenoxide ions are attached to this layer.  相似文献   

17.
The absolute bimolecular rate constants for the reactions of C6H5 with 2‐methylpropane, 2,3‐dimethylbutane and 2,3,4‐trimethylpentane have been measured by cavity ringdown spectrometry at temperatures between 290 and 500 K. For 2‐methylpropane, additional measurements were performed with the pulsed laser photolysis/mass spectrometry, extending the temperature range to 972 K. The reactions were found to be dominated by the abstraction of a tertiary C H bond from the molecular reactant, resulting in the production of a tertiary alkyl radical: C6H5 + CH(CH3)3 → C6H6 + t‐C4H9 (1) (1) C6H5 + (CH3)2CHCH(CH3)2 → C6H6 + t‐C6H13 (2) (2) C6H5 + (CH3)2CHCH(CH3)CH(CH3)2 → C6H6 + t‐C8H17 (3) (3) with the following rate constants given in units of cm3 mol−1 s−1: k1 = 10(11.45 ± 0.18) e−(1512 ± 44)/T k2 = 10(11.72 ± 0.15) e−(1007 ± 124)/T k3 = 10(11.83 ± 0.13) e−(428 ± 108)/T © 1999 John Wiley & Sons, Inc. Int J Chem Kinet 31: 645–653, 1999  相似文献   

18.
The synthesis of ansa complexes has been studied intensively owing to their importance as homogeneous catalysts and as precursors of metal‐containing polymers. However, paramagnetic non‐metallocene derivatives are rare and have been limited to examples with vanadium and titanium. Herein, we report an efficient procedure for the selective dilithiation of paramagnetic sandwich complex [Cr(η5‐C5H5)(η6‐C6H6)], which allows the preparation of a series of [n]chromoarenophanes (n=1, 2, 3) that feature silicon, germanium, and tin atoms at the bridging positions. The electronic and structural properties of these complexes were probed by X‐ray diffraction analysis, cyclic voltammetry, and by UV/Vis and EPR spectroscopy. The spectroscopic parameters for the strained and less strained complexes (i.e., with multiple‐atom linkers) indicate that the unpaired electron resides primarily in a d orbital on chromium(I); this result was also supported by density functional theory (DFT) calculations. We did not observe a correlation between the experimental UV/Vis and EPR data and the degree of molecular distortion in these ansa complexes. The treatment of tin‐bridged complex [Cr(η5‐C5H4)(η6‐C6H5)SntBu2] with [Pt(PEt3)3] results in the non‐regioselective insertion of the low‐valent Pt0 fragment into the Cipso? Sn bonds in both the five‐ and six‐membered rings, thereby furnishing a bimetallic complex. This observed reactivity suggests that ansa complexes of this type are promising starting materials for the synthesis of bimetallic complexes in general and also underline their potential to undergo ring‐opening processes to yield new metal‐containing polymers.  相似文献   

19.
Quantum chemical results will be presented on drugs, carcinogens, teratogens, and endogenous biomolecules using our new nonempirical ab initio MODPOT /VRDDO method, which incorporates as options to our ab initio LCAO -MO -SCF /CI programs ab initio effective core model potentials (MODPOT ) permitting one to calculate only the valence electrons explicitly yet accurately and an integral prescreening technique (VRDDO , variable retention of diatomic differential overlap) especially effective for spatially extended molecules. For molecules of the size of those of interest the MODPOT /VRDDO calculations run an order-of-magnitude faster than with our own fast ab initio programs and still retain accuracy to the third decimal place for the valence electron properties. We have also just implemented a new efficient MERGE technique which allows us to reuse integrals from a common skeletal fragment and only to have to recalculate those for a new atom or group or a change in its position. Examples will be presented of the use of this technique on a carcinogenic polycyclic aromatic hydrocarbon and its metabolites. The pKa's, oil-water partition, and drug distribution coefficients as a sensitive function of pH have been measured for a number of drugs as well as for relevant endogenous biomolecules. The pH dependence of the lipophilicities of such molecules has profound implication on appropriate use of such data in QSAR studies.  相似文献   

20.
The title compound, raloxifene hydro­chloride, C28H28NO4S+·Cl?, belongs to the benzo­thio­phene class of antiosteoporotic drugs. In the molecular cation, the 2‐phenol ring sustains a dihedral angle of 45.3 (1)° relative to the benzo­[b]­thio­phene system. The benzo­[b]­thio­phene and phenyl ring planes are twisted with respect to the carbonyl plane, with the smallest twist component occurring between the phenyl and carbonyl planes. The N atom bears the positive charge in the molecular cation and the piperidine ring adopts an almost perfect chair conformation. The Cl? anion is involved in the formation of N—H?Cl and O—H?Cl intermolecular hydrogen bonds, which lead to the formation of a layer of molecular cations.  相似文献   

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

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