The stable configurations and hydrogen bond nature of (H2O)n clusters (n = 3-6) have been investigated by the B3LYP method at the 6-31++g^** level. Upon calculation, four conclusions have been drawn: (1) In...The stable configurations and hydrogen bond nature of (H2O)n clusters (n = 3-6) have been investigated by the B3LYP method at the 6-31++g^** level. Upon calculation, four conclusions have been drawn: (1) In the (H2O)3-5 clusters, cyclic configurations were confirmed to be the most stable. But in the (H2O)3-4 ones, only cyclic configurations could be observed. From n = 5 ((H2O)5 clusters), three-dimensional configuration could be found: (2) In the (H2O)6 clusters, all configurations are inclined to be three-dimensional except the most stable configuration which is cyclic; (3) The stable order of (H2O)6 clusters indicates that it is the arrangement of hydrogen bond that plays a decisive role in the cluster stabilities, the zero-point energy is also important, and cluster stabilities are independent on the number of hydrogen bonds; (4) There exist strong cooperativity and superadditivity in the (H2O)n clusters.展开更多
The title complex [Cu(L1)(L2)(H2O)]·H2O(1,HL1 = N-(imino(pyridin-2-yl)me-thyl)picolinamidine),HL2 = salicylic acid) has been obtained by volatilization method with L1 prepared from 2,4,6-tripyridyl-1...The title complex [Cu(L1)(L2)(H2O)]·H2O(1,HL1 = N-(imino(pyridin-2-yl)me-thyl)picolinamidine),HL2 = salicylic acid) has been obtained by volatilization method with L1 prepared from 2,4,6-tripyridyl-1,3,5-triazine in situ.1 was fully characterized by single-crystal X-ray diffraction,elemental analysis and FT-IR.This complex exhibits a three-dimensional frame-work constructed through hydrogen bonding and C-H···π stacking interactions.The cyclic voltametric behavior of complex 1 was also investigated.1 belongs to the monoclinic system,space group P21/c with a = 15.112(5),b = 7.115(2),c = 19.899(6) ,β = 112.32°,V = 1979.4(11) 3,Mr = 460.94,Dc = 1.540 g/cm3,F(000) = 948,μ = 1.146 mm-1,Z = 4,the final R = 0.0612 and wR = 0.1813 for 2510 observed reflections with I 2σ(I).展开更多
A novel model for the 1: 2 complex of piperazine with some phenols in solution is established and verified. In CDCl3 solution, one piperazine molecule is tied to two phenol molecules by hydrogen bonds of O-H---N and N...A novel model for the 1: 2 complex of piperazine with some phenols in solution is established and verified. In CDCl3 solution, one piperazine molecule is tied to two phenol molecules by hydrogen bonds of O-H---N and N-H---O. And the protons of >NH and -OH groups exchange quickly and simultaneously between the atoms of phenol oxygen and piperazine nitrogen.展开更多
Three new crystalline compounds 1-3 were successfully obtained by the reactions of 3,3'-dimethoxy-6,6'-dimethyl-2,2'-bipyridine ligand(dmbp) with the corresponding Cu(Ⅰ) salts.Crystal data for 1:orthorhombic ...Three new crystalline compounds 1-3 were successfully obtained by the reactions of 3,3'-dimethoxy-6,6'-dimethyl-2,2'-bipyridine ligand(dmbp) with the corresponding Cu(Ⅰ) salts.Crystal data for 1:orthorhombic Pbca,a = 18.5858(12),b = 8.1821(5),c = 20.6066(13) ,V = 3133.7(3) 3,Z = 8,Dc = 1.843 g/cm3,F(000) = 1696,μ = 3.366 mm-1,the final R = 0.0223 and wR = 0.0542.Crystal data for 2:Orthorhombic Pbca,a = 18.7883(16),b = 8.3249(7),c = 19.0294(17) ,V = 2976.4(4) 3,Z = 8,Dc = 1.731 g/cm3,F(000) = 1552,μ = 4.154 mm-1,the final R = 0.0279 and wR = 0.0680.Crystal data for 3:monoclinic P21/c,a = 13.812(10),b = 9.910(7),c = 23.444(17) ,β = 104.3350(10)°,V = 3090(4) 3,Z = 4,Dc = 1.476 g/cm3,F(000) = 1408,μ = 1.588 mm-1,the final R = 0.0479 and wR = 0.1081.The results of X-ray crystallographic analysis revealed that C14H16ICuN2O2(1) and C14H16BrCuN2O2(2) are isostructural compounds with the dimers connected by C-H···halogen hydrogen bonds to generate a three-dimensional(3D) supramolecular network in 1 and a two-dimensional(2D) sheet structure in 2,respectively,while the mononuclear complex C28H32Cl2Cu2N4O4(3) is ionic.In 3,the [Cu(dmbp)2]+ cations and [ClCuCl]-anions are connected by C-H···Cl hydrogen bonds to form a one-dimensional(1D) chain along the a axis.Therefore,in the three complexes,the C-H···halogen hydrogen bonds dominate their crystal structures.Additionally,The UV luminescent properties of complexes 1-3 were investigated.展开更多
The title compound (C17H16O2, Mr = 252.30) crystallizes in the monoclinic system, space group P21/c with a = 11.5763(14), b = 11.0321(11), c = 11.5094(13) ?, β = 114.581(3)o, V = 1336.7(3) ?3 , Z = 4, Dc = 1.254 g/...The title compound (C17H16O2, Mr = 252.30) crystallizes in the monoclinic system, space group P21/c with a = 11.5763(14), b = 11.0321(11), c = 11.5094(13) ?, β = 114.581(3)o, V = 1336.7(3) ?3 , Z = 4, Dc = 1.254 g/cm3, μ(MoKa) = 0.081 cm-1 and F(000) = 536.00. The final R and wR are 0.0527 and 0.1285, respectively for 3058 observed reflections (I > 2σ(I)). In the title molecule, two phenyl rings are rotated oppositely with respect to the central part C(1)–C(2)=C(3)– C(4) (plane 3) and the dihedral angle between them is 14.8o. The phenone O(1) atom deviates from plane 3 by 0.291 ?. In the crystal the molecules form H-bond chains of R2 (6) and R2 (5) types 1 1 along [001]. The molecule chains interacted through three cooperative C–H…O H-bonds (R3 (11)) 1 in the ‘bay area’ (Fig. 3), extending along [010] and forming layer (100). Between the layers, there exist C–H/π interactions along [101]. Studies on the cooperative C–H…O H-bonds R3 (11) in the 1 similar crystals are also presented.展开更多
The microscopic properties of NaCl-induced phase separation of acetonitrile (ACN)-water mixtures have been studied by proton nuclear magnetic resonance (1H NMR). Acetonitrile-rich phase increases with increasing NaCl ...The microscopic properties of NaCl-induced phase separation of acetonitrile (ACN)-water mixtures have been studied by proton nuclear magnetic resonance (1H NMR). Acetonitrile-rich phase increases with increasing NaCl concentration (cNaCl) at xACN ≈ 0.25. 1H chemical shift of water for acetonitrile-rich phase rapidly decreases with decreasing NaCl mole concentration and that for water-rich phase quickly increases with increasing cNaCl. However, 1H chemical shift of acetonitrile has nothing to do with the molar concentration of NaCl, and it keeps relatively stable for all solutions (±0.002). These results reveal that Na+ and Cl- are rapidly hydrated by water, not by acetonitrile. The change of 1H chemical shift of water has shown that the number of hydrogen bond increases or hydrogen bond strengths with increasing NaCl molarity in mixtures. But hydrogen bond is broken or weaken with the temperature rising. 1H chemical shifts of pure water and the water in acetonitrile-rich phase have been investigated at 293 K, 298 K and 303 K. The hydration number of Na+ (6.05) in water-rich phase is determined by an empirical equation involving 1H chemical shift, temperature and NaCl molarity, which is in good agreement with the literatures.展开更多
By employing the dithiocarbamate salt (K(PPDC), where PPDC = 4?-pyridyl-1-pipe-razine-4-dithiocarbamate) as the functional ligand and di-palladium complexes [(N^N)2Pd2(NO3ˉ)2](NO3ˉ)2 (N^N = 2,2'-bipyri...By employing the dithiocarbamate salt (K(PPDC), where PPDC = 4?-pyridyl-1-pipe-razine-4-dithiocarbamate) as the functional ligand and di-palladium complexes [(N^N)2Pd2(NO3ˉ)2](NO3ˉ)2 (N^N = 2,2'-bipyridine, bpy; 4,4'-dimethylbipyridine, dmbpy) as corner, two novel single metal complexes with Pd(II) centers have been obtained. These organic-metal complexes were characterized by NMR, ESI-MS, elemental analysis, Uv-vis spectra and single-crystal X-ray diffraction analysis. Compound 1?(PF6)2 ([(bpy)Pd(PPDC)]·(PF6)2) crystallizes in triclinic, space group P , a = 8.3968(5), b = 11.5565(7), c = 18.2234(11), α = 97.505(1), β = 91.424(1), γ = 106.146(1)o, C22H24N6S2P2F12Pd, Mr = 832.93, V = 1680.58(18) ?3, Z = 2, Dc = 1.646 Mg/m3, μ(MoKα) = 0.863 mm-1, F(000) = 828, the final R = 0.0455 and wR = 0.1390 for 6981 observed reflections with I 〉 2σ(I). Similarly, compound 2?(PF6)2 ([(dmbpy)Pd(PPDC)]·(PF6)2) also crystallizes in triclinic, space group P , a = 13.9467(3), b = 14.8390(2), c = 17.0632(3) ?, α = 81.8680(10), β = 87.051(2), γ = 83.4590(10)o, C22H25N5S2P2F12Pd, Mr = 819.93, V = 3470.81(11) ?3, Z = 4, Dc = 1.569 Mg/m3, μ(CuKα) = 7.115 mm-1, F(000) = 1632, the final R = 0.0606 and wR = 0.1637 for 12835 observed reflections with I 〉 2σ(I). Crystallography reveals that each metal center coordinates with two N atoms from bpy and two S atoms from PPDC in the square coordination mode. In the crystal structure of complex 1, a weak Pd……Pd interaction can be observed. Interestingly, it was also found that the mononuclear moieties of complex 2 could be packed into a 3-D porous framework via multiple intermolecular C–F……H hydrogen-boding interactions which extended in the a, b, and c axes with PF6ˉ anions frozen inside.展开更多
运用量子化学从头算方法,在MP2/6-311++G(d,p),MP2/6-311++G(2df,2p),MP2/6-311++G(3df,3pd)和QCISD/6-311++G(d,p)水平上,研究了CH3F,CH3Cl和CH3Br作为质子给体与Cl-,Br-作为质子接受体形成的氢键CH3…Y.计算结果表明:6种复合物中C—H...运用量子化学从头算方法,在MP2/6-311++G(d,p),MP2/6-311++G(2df,2p),MP2/6-311++G(3df,3pd)和QCISD/6-311++G(d,p)水平上,研究了CH3F,CH3Cl和CH3Br作为质子给体与Cl-,Br-作为质子接受体形成的氢键CH3…Y.计算结果表明:6种复合物中C—H键收缩,伸缩振动频率增大,形成蓝移氢键.分子中原子(AtomsinMole-cules,AIM)分析表明,这些复合物的电子密度拓扑性质与普通氢键有着本质的不同,在Y…H之间不存在键临界点,而在Y与C之间存在键临界点,因此这些相互作用严格地不能称为氢键.自然键轨道(Natural bond orbital,NBO)分析表明,在这些复合物中弯曲的CH…Y的特殊结构使得分子间超共轭n(Y)→σ*(C—H)减小到可以忽略;质子接受体的电子密度没有转移到σ*(C—H)上,而是转移到了σ*(C—X)(X=F,Cl,Br)上;存在一定程度的重杂化;分子内超共轭相互作用减小使得σ*(C—H)的电子密度减少.这些因素共同导致C—H伸缩振动频率的蓝移.展开更多
采用量子化学计算和分子动力学(MD)模拟,研究了溶剂化物环四亚甲基四硝胺(HMX)/N,N-二甲基甲酰胺(DMF)的结构和分子间相互作用。对浸渍在DMF溶剂中的β-HMX分子的MD模拟表明,HMX的分子构象已经从β转变为α相。在M P2/6-31G*水平上的理...采用量子化学计算和分子动力学(MD)模拟,研究了溶剂化物环四亚甲基四硝胺(HMX)/N,N-二甲基甲酰胺(DMF)的结构和分子间相互作用。对浸渍在DMF溶剂中的β-HMX分子的MD模拟表明,HMX的分子构象已经从β转变为α相。在M P2/6-31G*水平上的理论计算也说明,在D M F溶剂中,α-HM X比β-HM X更稳定。这解释了在HM X/D M F溶剂化物的多晶型中,所有的HMX分子均呈α构象的原因。采用MP2/6-31G*方法,对α-HMX和DMF分子间可能的同型和异型二聚体进行结构优化。结果表明,组分间存在C—H…O氢键相互作用,并且α-HMX/DMF的稳定化能非常接近α-HMX/α-HMX,并远远大于D M F二聚体的。这意味着异型分子间力可和同型分子间力竞争。从热力学的观点来看,共结晶过程可能发生。对α-HM X在D M F中的过饱和溶液的M D模拟表明,分子间相互作用对共结晶有利。这些理论研究对理解在HM X的D M F溶液中为什么发生的是共结晶而不是重结晶提供了有价值的信息。展开更多
基金Project supported by the Natural Science Foundation of Tangshan Teacher’s College (No. 04C06)
文摘The stable configurations and hydrogen bond nature of (H2O)n clusters (n = 3-6) have been investigated by the B3LYP method at the 6-31++g^** level. Upon calculation, four conclusions have been drawn: (1) In the (H2O)3-5 clusters, cyclic configurations were confirmed to be the most stable. But in the (H2O)3-4 ones, only cyclic configurations could be observed. From n = 5 ((H2O)5 clusters), three-dimensional configuration could be found: (2) In the (H2O)6 clusters, all configurations are inclined to be three-dimensional except the most stable configuration which is cyclic; (3) The stable order of (H2O)6 clusters indicates that it is the arrangement of hydrogen bond that plays a decisive role in the cluster stabilities, the zero-point energy is also important, and cluster stabilities are independent on the number of hydrogen bonds; (4) There exist strong cooperativity and superadditivity in the (H2O)n clusters.
基金supported by the University Science Foundation of Anhui Province (No.KJ2009B104)the Applied Chemistry Key Constructing Subject of Anhui Province (No.200802187C)
文摘The title complex [Cu(L1)(L2)(H2O)]·H2O(1,HL1 = N-(imino(pyridin-2-yl)me-thyl)picolinamidine),HL2 = salicylic acid) has been obtained by volatilization method with L1 prepared from 2,4,6-tripyridyl-1,3,5-triazine in situ.1 was fully characterized by single-crystal X-ray diffraction,elemental analysis and FT-IR.This complex exhibits a three-dimensional frame-work constructed through hydrogen bonding and C-H···π stacking interactions.The cyclic voltametric behavior of complex 1 was also investigated.1 belongs to the monoclinic system,space group P21/c with a = 15.112(5),b = 7.115(2),c = 19.899(6) ,β = 112.32°,V = 1979.4(11) 3,Mr = 460.94,Dc = 1.540 g/cm3,F(000) = 948,μ = 1.146 mm-1,Z = 4,the final R = 0.0612 and wR = 0.1813 for 2510 observed reflections with I 2σ(I).
文摘A novel model for the 1: 2 complex of piperazine with some phenols in solution is established and verified. In CDCl3 solution, one piperazine molecule is tied to two phenol molecules by hydrogen bonds of O-H---N and N-H---O. And the protons of >NH and -OH groups exchange quickly and simultaneously between the atoms of phenol oxygen and piperazine nitrogen.
基金supported by the National Natural Science Foundation of China (No.20872057)
文摘Three new crystalline compounds 1-3 were successfully obtained by the reactions of 3,3'-dimethoxy-6,6'-dimethyl-2,2'-bipyridine ligand(dmbp) with the corresponding Cu(Ⅰ) salts.Crystal data for 1:orthorhombic Pbca,a = 18.5858(12),b = 8.1821(5),c = 20.6066(13) ,V = 3133.7(3) 3,Z = 8,Dc = 1.843 g/cm3,F(000) = 1696,μ = 3.366 mm-1,the final R = 0.0223 and wR = 0.0542.Crystal data for 2:Orthorhombic Pbca,a = 18.7883(16),b = 8.3249(7),c = 19.0294(17) ,V = 2976.4(4) 3,Z = 8,Dc = 1.731 g/cm3,F(000) = 1552,μ = 4.154 mm-1,the final R = 0.0279 and wR = 0.0680.Crystal data for 3:monoclinic P21/c,a = 13.812(10),b = 9.910(7),c = 23.444(17) ,β = 104.3350(10)°,V = 3090(4) 3,Z = 4,Dc = 1.476 g/cm3,F(000) = 1408,μ = 1.588 mm-1,the final R = 0.0479 and wR = 0.1081.The results of X-ray crystallographic analysis revealed that C14H16ICuN2O2(1) and C14H16BrCuN2O2(2) are isostructural compounds with the dimers connected by C-H···halogen hydrogen bonds to generate a three-dimensional(3D) supramolecular network in 1 and a two-dimensional(2D) sheet structure in 2,respectively,while the mononuclear complex C28H32Cl2Cu2N4O4(3) is ionic.In 3,the [Cu(dmbp)2]+ cations and [ClCuCl]-anions are connected by C-H···Cl hydrogen bonds to form a one-dimensional(1D) chain along the a axis.Therefore,in the three complexes,the C-H···halogen hydrogen bonds dominate their crystal structures.Additionally,The UV luminescent properties of complexes 1-3 were investigated.
基金This research was supported by the Key Subject Foundation of Jiangsu Province(S1109001)
文摘The title compound (C17H16O2, Mr = 252.30) crystallizes in the monoclinic system, space group P21/c with a = 11.5763(14), b = 11.0321(11), c = 11.5094(13) ?, β = 114.581(3)o, V = 1336.7(3) ?3 , Z = 4, Dc = 1.254 g/cm3, μ(MoKa) = 0.081 cm-1 and F(000) = 536.00. The final R and wR are 0.0527 and 0.1285, respectively for 3058 observed reflections (I > 2σ(I)). In the title molecule, two phenyl rings are rotated oppositely with respect to the central part C(1)–C(2)=C(3)– C(4) (plane 3) and the dihedral angle between them is 14.8o. The phenone O(1) atom deviates from plane 3 by 0.291 ?. In the crystal the molecules form H-bond chains of R2 (6) and R2 (5) types 1 1 along [001]. The molecule chains interacted through three cooperative C–H…O H-bonds (R3 (11)) 1 in the ‘bay area’ (Fig. 3), extending along [010] and forming layer (100). Between the layers, there exist C–H/π interactions along [101]. Studies on the cooperative C–H…O H-bonds R3 (11) in the 1 similar crystals are also presented.
文摘The microscopic properties of NaCl-induced phase separation of acetonitrile (ACN)-water mixtures have been studied by proton nuclear magnetic resonance (1H NMR). Acetonitrile-rich phase increases with increasing NaCl concentration (cNaCl) at xACN ≈ 0.25. 1H chemical shift of water for acetonitrile-rich phase rapidly decreases with decreasing NaCl mole concentration and that for water-rich phase quickly increases with increasing cNaCl. However, 1H chemical shift of acetonitrile has nothing to do with the molar concentration of NaCl, and it keeps relatively stable for all solutions (±0.002). These results reveal that Na+ and Cl- are rapidly hydrated by water, not by acetonitrile. The change of 1H chemical shift of water has shown that the number of hydrogen bond increases or hydrogen bond strengths with increasing NaCl molarity in mixtures. But hydrogen bond is broken or weaken with the temperature rising. 1H chemical shifts of pure water and the water in acetonitrile-rich phase have been investigated at 293 K, 298 K and 303 K. The hydration number of Na+ (6.05) in water-rich phase is determined by an empirical equation involving 1H chemical shift, temperature and NaCl molarity, which is in good agreement with the literatures.
基金supported by the National Natural Science Foundation of China(No.21471011)High-caliber Talents Project of Beijing Municipal Institution
文摘By employing the dithiocarbamate salt (K(PPDC), where PPDC = 4?-pyridyl-1-pipe-razine-4-dithiocarbamate) as the functional ligand and di-palladium complexes [(N^N)2Pd2(NO3ˉ)2](NO3ˉ)2 (N^N = 2,2'-bipyridine, bpy; 4,4'-dimethylbipyridine, dmbpy) as corner, two novel single metal complexes with Pd(II) centers have been obtained. These organic-metal complexes were characterized by NMR, ESI-MS, elemental analysis, Uv-vis spectra and single-crystal X-ray diffraction analysis. Compound 1?(PF6)2 ([(bpy)Pd(PPDC)]·(PF6)2) crystallizes in triclinic, space group P , a = 8.3968(5), b = 11.5565(7), c = 18.2234(11), α = 97.505(1), β = 91.424(1), γ = 106.146(1)o, C22H24N6S2P2F12Pd, Mr = 832.93, V = 1680.58(18) ?3, Z = 2, Dc = 1.646 Mg/m3, μ(MoKα) = 0.863 mm-1, F(000) = 828, the final R = 0.0455 and wR = 0.1390 for 6981 observed reflections with I 〉 2σ(I). Similarly, compound 2?(PF6)2 ([(dmbpy)Pd(PPDC)]·(PF6)2) also crystallizes in triclinic, space group P , a = 13.9467(3), b = 14.8390(2), c = 17.0632(3) ?, α = 81.8680(10), β = 87.051(2), γ = 83.4590(10)o, C22H25N5S2P2F12Pd, Mr = 819.93, V = 3470.81(11) ?3, Z = 4, Dc = 1.569 Mg/m3, μ(CuKα) = 7.115 mm-1, F(000) = 1632, the final R = 0.0606 and wR = 0.1637 for 12835 observed reflections with I 〉 2σ(I). Crystallography reveals that each metal center coordinates with two N atoms from bpy and two S atoms from PPDC in the square coordination mode. In the crystal structure of complex 1, a weak Pd……Pd interaction can be observed. Interestingly, it was also found that the mononuclear moieties of complex 2 could be packed into a 3-D porous framework via multiple intermolecular C–F……H hydrogen-boding interactions which extended in the a, b, and c axes with PF6ˉ anions frozen inside.
文摘运用量子化学从头算方法,在MP2/6-311++G(d,p),MP2/6-311++G(2df,2p),MP2/6-311++G(3df,3pd)和QCISD/6-311++G(d,p)水平上,研究了CH3F,CH3Cl和CH3Br作为质子给体与Cl-,Br-作为质子接受体形成的氢键CH3…Y.计算结果表明:6种复合物中C—H键收缩,伸缩振动频率增大,形成蓝移氢键.分子中原子(AtomsinMole-cules,AIM)分析表明,这些复合物的电子密度拓扑性质与普通氢键有着本质的不同,在Y…H之间不存在键临界点,而在Y与C之间存在键临界点,因此这些相互作用严格地不能称为氢键.自然键轨道(Natural bond orbital,NBO)分析表明,在这些复合物中弯曲的CH…Y的特殊结构使得分子间超共轭n(Y)→σ*(C—H)减小到可以忽略;质子接受体的电子密度没有转移到σ*(C—H)上,而是转移到了σ*(C—X)(X=F,Cl,Br)上;存在一定程度的重杂化;分子内超共轭相互作用减小使得σ*(C—H)的电子密度减少.这些因素共同导致C—H伸缩振动频率的蓝移.
基金National Natural Science Foundation of C hina(NSAF,No.11176029)the Key Research Project of Education Department of SichuanProvince(No.10zd1106)
文摘采用量子化学计算和分子动力学(MD)模拟,研究了溶剂化物环四亚甲基四硝胺(HMX)/N,N-二甲基甲酰胺(DMF)的结构和分子间相互作用。对浸渍在DMF溶剂中的β-HMX分子的MD模拟表明,HMX的分子构象已经从β转变为α相。在M P2/6-31G*水平上的理论计算也说明,在D M F溶剂中,α-HM X比β-HM X更稳定。这解释了在HM X/D M F溶剂化物的多晶型中,所有的HMX分子均呈α构象的原因。采用MP2/6-31G*方法,对α-HMX和DMF分子间可能的同型和异型二聚体进行结构优化。结果表明,组分间存在C—H…O氢键相互作用,并且α-HMX/DMF的稳定化能非常接近α-HMX/α-HMX,并远远大于D M F二聚体的。这意味着异型分子间力可和同型分子间力竞争。从热力学的观点来看,共结晶过程可能发生。对α-HM X在D M F中的过饱和溶液的M D模拟表明,分子间相互作用对共结晶有利。这些理论研究对理解在HM X的D M F溶液中为什么发生的是共结晶而不是重结晶提供了有价值的信息。
基金supported by the National Natural Science Foundation of China(21002006,20452002)Special Program for Key Basic Research of the Ministry of Science and Technology,China(2004-973-36)~~