簇合物(p-H2)N-HCCCN的结构和能量

采用两体作用势模型和遗传算法对簇合物(p-H2)NHCCCN的极小能量结构和能量进行了理论研究.结果表明,20个para-H2分子形成HCCCN周围的第一个溶剂层,第一个溶剂层包含三个溶剂环,每个溶剂环都有6个para-H2分子,第19和20个para-H2分子分别聚集在HCCCN分子的N、H原子末端.进一步计算了(p-H2)N-HCCCN的化学势,发现化学势随para-H2分子个数的增加呈震荡变化.     

簇合物(p-H2)N-HCCCN的结构和能量

采用两体作用势模型和遗传算法对簇合物（p-H₂）_N-HCCCN的极小能量结构和能量进行了理论研究.结果表明,20 个para-H₂分子形成HCCCN周围的第一个溶剂层,第一个溶剂层包含三个溶剂环,每个溶剂环都有6 个para-H₂分子,第19和20个para-H₂分子分别聚集在HCCCN分子的N、H原子末端. 进一步计算了（p-H₂）_N-HCCCN的化学势,发现化学势随para-H₂分子个数的增加呈震荡变化.     

成庄无烟煤大分子结构模型及其分子模拟

利用晋城矿区成庄矿煤的工业分析、元素分析、¹³C-NMR、XPS等实验结果,构建了其大分子结构模型.模型中,芳香碳以2、3、4环结构为主,最大环数达五个;脂肪碳以甲、乙基侧链及环烷烃的形式存在.九个氧原子分别以七个羰基(主要为醌基)、一个羟基及一个醚氧的型式存在;两个氮原子以吡咯的形式存在.硫原子含量很低,在模型构建中没有体现.采用分子力学(MM)和分子动力学(MD)方法,对成庄煤结构模型进行能量最小化模拟.结果表明,分子内及分子间芳香层片之间的π-π相互作用,使其以近似平行的方式排列;高煤级煤结构中,短程有序的原因主要是分子间芳香层片的定向排列.分子间的氢键能及范德华能使结构达到最稳构型.     

BH2+与H2O反应机理的量子拓扑研究

采用密度泛函方法B3LYP/6-311+G(d, p)和耦合簇方法CCSD/6-311+G(d, p)研究了BH2+与H2O的气相离子-分子反应机理. 优化得到了反应途径中各驻点的几何构型, 并采用内禀反应坐标法进行追踪. 从量子拓扑学的角度, 讨论了在反应过程中各化学键的变化. 反应(I)经历了一个四元环过渡态, 找到了这个反应的能量过渡态和两个结构过渡态.     

聚集诱导发光材料与纳米晶复合体系动力学新进展北大核心CSCD

近期,中科院大连化物所吴凯丰研究员团队将聚集诱导发光分子(AIEgen)嫁接到纳米晶表面,并研究了这一复合体系的激发态动力学,发现这一复合体系中AIEgen的非辐射分子内运动可以得到有效抑制,这一普适性现象可用于构建各类多功能发光材料。相关工作发表于《物理化学快报》(Journal of Physical Chemistry Letters)上。传统的染料分子由于芳香环的π-π堆积通常表现出聚集诱导猝灭现象(ACQ),阻碍了这些染料分子在不良溶剂和固体发光器件中的广泛应用。AIEgen则刚好相反:在良溶剂中,各种分子内运动可耗散激发态能量,导致发光效率极低;在不良溶剂或者固态薄膜中,分子的团聚有效抑制了上述非辐射分子内运动,表现出强烈的发光行为。     

分子动力学模拟纳米尺寸限制体系下氩溶液中I2的振动能量弛豫

利用平衡态分子动力学方法(EMD)模拟了纳米尺寸限制球壳内I2在Ar溶液中的振动能量转移.计算并讨论了I2振动能量弛豫时间T1随球壳半径、溶剂密度的变化规律.通过分子间相互作用分析,在原子、分子水平上,揭示了随着球壳半径的减小,T1呈逐渐增大趋势的原因.结果表明,球壳的几何限制效应和表面作用对受限溶液密度分布的影响较大,从而导致溶质振动弛豫的显著变化.此外,非限制体系模拟显示,非平衡态分子动力学(NEMD)方法可以得到与平衡态分子动力学方法较一致的振动能量弛豫时间T1.     

分子动力学模拟纳米尺寸限制体系下氩溶液中I2的振动能量弛豫

利用平衡态分子动力学方法(EMD)模拟了纳米尺寸限制球壳内I2在Ar溶液中的振动能量转移. 计算并讨论了I2振动能量弛豫时间T1随球壳半径、溶剂密度的变化规律. 通过分子间相互作用分析, 在原子、分子水平上, 揭示了随着球壳半径的减小, T1呈逐渐增大趋势的原因. 结果表明, 球壳的几何限制效应和表面作用对受限溶液密度分布的影响较大, 从而导致溶质振动弛豫的显著变化. 此外, 非限制体系模拟显示, 非平衡态分子动力学(NEMD)方法可以得到与平衡态分子动力学方法较一致的振动能量弛豫时间T1.     

带极性侧链的环[6]芳酰胺的球形自组装

环芳酰胺是一类基于三中心氢键促进,经寡聚前体一步大环合成法得到的刚性大环分子.通过紫外-可见(UV-Vis)光谱、动态光散射(DLS)、扫描电镜(SEM)、透射电镜(TEM)和原子力显微镜(AFM)等实验手段,详细考察了侧链为三甘醇单甲基醚链,由六个苯环单元组成的环[6]芳酰胺的自组装行为.实验结果表明,该大环在1,2-二氯乙烷中发生自组装,其组装聚集体随温度升高产生从聚集体到单分子的解聚变化,至70℃时几乎完全解聚;在由良溶剂(二氯甲烷)和不良溶剂(芳烃类)组成的混合溶剂中,带有三甘醇醚链的环[6]芳酰胺化合物1自组装成微球,结合热稳定性实验和TEM证实是实心微球而非囊泡.进一步发现微球形成和形貌依赖于混合溶剂中不良溶剂的极性和种类,芳烃类溶剂有利于微球形成,而烷烃和极性溶剂则不利,后者更倾向于形成膜的结构.     

溶菌酶蛋白与聚合物防污膜相互作用的分子动力学模拟

采用分子动力学方法研究了溶菌酶蛋白分子(Lysozyme)在2种典型聚合物防污材料(有机硅弹性体聚二甲基硅氧烷PDMS和两性离子类聚磺酸基甜菜碱甲基丙烯酸甲酯SBMA)表面的吸附行为,进一步从微观角度阐释了防污材料的防污机理.通过比较蛋白质与聚合物膜间的作用力和结合能,防污膜表面水化层的动力学性质,以及蛋白质与基底结合位点附近的结构分析表明SBMA有着更优异的防污能力:(1)蛋白质的吸附须要克服两者表面水化层引起的物理障碍和能量势垒,SBMA通过表面氢键、静电作用和笼效应束缚了一层紧密结合的水化层,表面结合水难于脱附,水化层分子的去溶剂化需要克服的能垒高.(2)蛋白质与PDMS的结合能量上更具优势,相比SBMA与蛋白质间的结合更加稳定,不利于蛋白质的脱附.     

RDX分子内振动能量重分配的动力学研究

基于从头算分子动力学(Born-oppenheimer molecular dynamics, BOMD)模拟, 构建了环硝胺六氢-1,3,5-三硝基-1,3,5-三嗪(RDX)单分子不同振动模式之间的耦合矩阵, 并计算了在不同加载能量下从低频振动模式到高频振动模式的最优能量传输路径. 结果表明, RDX单分子中—NNO₂基团更有利于能量局域化, 振动模式v₃和v₄在从低频振动模式到高频振动模式的能量传输过程中扮演着重要角色. 通过对v₃和v₄两个振动模式的进一步分析发现, 加载能量的不同会导致RDX单分子能量传输路径的不同. 当加载能量较低时, RDX单分子倾向于从低频振动模式到中频振动模式再到高频振动模式的能量传输路径; 当加载能量较高时, 能量更倾向于从低频振动模式直接传输到高频振动模式上. 揭示了RDX分子内振动耦合能量转移的微观机制, 为进一步探索RDX将“机械能”转化为“化学能”的微观过程提供了理论基础.     

Computer simulation of the supercritical carbon dioxide fluid (II) internal energy and structure of carbon dioxide system containing one benzene molecule

Using potential models based on ab initio quantum chemical calculations, we study a supercritical CO₂ fluid containing one benzene molecule using Monte Carlo simulations. First, molecular average internal energy is calculated for the whole system and for the first solvation shell of the benzene molecule. This analysis shows that the CO₂ molecules in the first solvation shell have a large energetic stabilization owing to the shape of the solute. In addition to the stabilization, the solute-solvent interactions in the first solvation shell show large fluctuations for both the in-plane and out-of-plane parts. Secondly, an orientational distribution function is defined to investigate the CO₂ fluid structure. This function indicates that the CO₂---CO₂ intermolecular configuration has a large dependence on the temperature of the system for both the whole system, and for the first solvation shell of the solute. Moreover, the benzene molecule is confirmed to control the mutual arrangement between neighboring CO₂ molecules.     

Selective 3,4-Polymerization Mechanism of Isoprene Catalyzed by Rare Earth Alkyl Complexes

The mechanism of selective 3,4-polymerization reaction of isoprene catalyzed by the rare earth lutecium（Ⅲ） alkyl complexes [2,6-Me2Ph-N-CH2-C（CH2SiMe3）=N-PhMe2-2,6]Lu（CH2SiMe3）2（THF）was investigated by means of the M06/sdd method with solvation effects taken into account.The results show that the structure of the catalyst core remained almost unchanged as the isoprene molecules were alternatively inserted into the complex at two opposite sides.The Gibbs free energies of the coordination complexes,transition state and intermediates indicate that all the isoprene molecules prefer to insert into the complex with the 3,4-polymerization selectivity as catalyzed by the catalyst,which is consistent with the experimental observations.It is found that the insertion reaction of each isoprene is exothermic,which comes mainly from the coordination of the isoprene molecule to the lutecium（Ⅲ） atom.The solvation effects were confirmed important in predicting the Gibbs free energies of the present reaction system.     

Solvation of tetraalkylammonium chlorides in acetonitrile-water mixtures: mass spectrometry and molecular dynamics simulations.

The solvation of tetramethylammonium chloride (Me4NCl) and tetra-n-butylammonium chloride (Bu4NCl) in water-acetonitrile mixtures was investigated by mass spectrometry of clusters isolated from the solution. As far as the positive ions are concerned, clusters composed of alkylammonium ions and acetonitrile molecules only were observed, even for mixtures with high water content. In contrast, for the negative ions, clusters composed of chloride with both water and/or acetonitrile molecules were observed. For the smaller system (Me4NCl) we performed quantum chemical calculations and molecular dynamics simulations. It was found that even though water is present in the solvation shell of Me4N+, only acetonitrile has a strong electrostatic interaction with the cation. Water molecules around Me4N+ form hydrogen bonds with other water molecules, and they interact with Me4N+ mainly via dispersive interactions. These results indicate that Me4N+ behaves like a hydrophobic solute. On the other hand, the interaction of Cl- with water and acetonitrile is of comparable strength and, in both cases, the electrostatic interaction dominates. Herein we demonstrate experimentally and theoretically that positive and negative ions give rise to characteristic solvation structures in mixed solvents: even a relatively small organic cation, such as Me4N+, exhibits a hydrophobic-like solvation shell.     

MnCl_2改性γ-Al_2O_3毛细管填充柱气相色谱分析氢同位素

针对常规气相色谱填充柱分析稳定氢同位素的柱效低、峰宽大、保留时间长等问题,采用MnCl_2改性γ-Al_2O_3填充的石英毛细管柱开展了系统性柱效分析及氢同位素分析技术研究。研究结果表明,使用MnCl_2对γ-Al_2O_3进行改性后,可大大改善单纯的γ-Al_2O_3表面有序度、孔结构和吸附性质,并将正氢(o-H_2)和仲氢(p-H_2)峰洗脱在单一谱峰区域内。制备的长1.0 m、内径0.53 mm的石英毛细填充柱与热导检测器(TCD)级联测试,在体积浓度1至10 m L/L范围内有较好的线性关系,对于低浓度样品检测的相对误差不大于5%。H_2、HD和D_2的保留时间可分别缩短至39、46和60 s,检出限可分别降低至0.046、0.067和0.072 m L/L。毛细管填充柱较常规填充柱具有峰形尖锐、相邻组分分离度高、保留时间短、检出限低等优点,可用于低浓度氢同位素快速测量及氢同位素在线分析。     

5-巯基-1, 3, 4-噻二唑-2-硫酮微溶剂团簇的光谱和理论计算研究

微溶剂作用(即溶剂化过程)广泛存在于所有物理、 化学和生命过程中. 在液相化学反应体系中, 几乎是一切化学反应的基础. 通过傅里叶变换拉曼光谱(FT-Raman)并结合密度泛函理论(DFT), 表征了固态5-巯基- 1, 3, 4-噻二唑-2-硫酮(MTT)的结构, 并进一步确认了MTT在乙腈、 甲醇和水中微溶剂团簇的大小和氢键位点. 通过探究MTT在不同溶剂及pH条件下的紫外-可见吸收光谱(UV-Vis), 结合含时密度泛函理论(TD-DFT)计算, 揭示了溶剂和pH对MTT电子跃迁带的影响, 进一步解释了其光谱位移. 结合能量计算可以得出, MTT分别与1个乙腈、 2个甲醇和2个水分子形成MTT(CH₃CN), MTT(CH₃OH)₂ 和MTT(H₂O)₂团簇.     

Structure and stability of thiourea with water, DFT and MP2 calculations

The relative stabilities of thiourea in water are investigated computationally by considering thiourea–water complexes containing up to 1–6 water molecules (CS(NH₂)₂(H₂O)_n=1–6) using density functional theory and MP2 ab initio molecular orbital theory. The results show that the thiourea complex is stable and has an unusually high affinity for incoming water molecules. The clusters are progressively stabilized by the addition of water molecules, as indicated by the increasing of the binding energy. The binding energy of the cluster to each H₂O molecule is about 33 kJ mol⁻¹ for n=1–5.The C–S bond, N–C bond distance, Mulliken populations and binding energy keep approximately constant as the clusters increase in size with an increasing number of H₂O molecules. As the solvation progresses, the C–S distance increases monotonically while the Mulliken populations on the C–S bond reduces monotonically with the addition of each H₂O molecule, indicating that the C–S bond of the thiourea unit in the clusters is de-stabilized with an increasing number of H₂O molecules. Charge transfers for the clusters are mainly found at N, S atoms of the thiourea.     

Tautomeric properties, conformations and structure of N-(2-hydroxy-5-chlorophenyl) salicylaldimine

The crystal structure of N-(2-hydroxy-5-chlorophenyl) salicylaldimine (C₁₃H₁₀NO₂Cl) was determined by X-ray analysis. It crystallizes orthorhombic space group P2₁2₁2₁ with a=12.967(2) Å, b=14.438(3) Å, c=6.231(3) Å, V=1166.5(6) Å³, Z=4, D_c=1.41 g cm⁻³ and μ(MoK)=0.315 mm⁻¹. The title compound is thermochromic and the molecule is nearly planar. Both tautomeric forms (keto and enol forms in 68(3) and 32(3)%, respectively) are present in the solid state. The molecules contain strong intramolecular hydrogen bonds, N1–H1O1/O2 (2.515(1) and 2.581(2) Å) for the keto form and O1–H01N1 for the enol one. There is also strong intermolecular O2–HO1 hydrogen bonding (2.599(2) Å) between neighbouring molecules. Minimum energy conformations AM1 were calculated as a function of the three torsion angles, θ₁(N1–C7–C6–C5), θ₂(C8–N1–C7–C6) and θ₃(C9–C8–N1–C7), varied every 10°. Although the molecule is nearly planar, the AM1 optimized geometry of the title compound is not planar. The non-planar conformation of the title compound corresponding to the optimized X-ray structure is the most stable conformation in all calculations.     

DNA碱基(对)及其Zn2+复合物对CO2,N2,H2小分子的吸附

用M062X/6-31+G^*方法探讨了腺嘌呤(A)、 胸腺嘧啶(T)、 鸟嘌呤(G)、 胞嘧啶(C)及其碱基对(AT, GC)以及Zn²⁺复合物(AAA-Zn²⁺, AAT-Zn²⁺和GGC-Zn²⁺)对混合小分子H₂, N₂, CO₂的吸附情况, 系统研究了其相互作用模式及吸附强度, 预测了常见混合气体分子与碱基(对)及复合物的吸附位置. 研究表明, CO₂倾向于以氢键的形式结合到碱基(对)的氨基氢或亚氨基氢上, 而N₂和H₂分子则倾向于结合到这些碱基(对)的平面π电子上, 以堆垛的形式存在. 根据吸附强度大小, 预测了由这些碱基为骨架合成的金属有机骨架(MOF)吸附材料对小分子的选择性吸附顺序为H₂22. 研究表明, 以AT对结合金属Zn²⁺为节点的纯天然碱基对构成的MOF要比实验合成的AA碱基对与Zn²⁺结合的MOF具备更好的吸附和分离性能.     

钇对石墨烯储氢性能的影响

应用基于密度泛函理论的第一性原理方法研究过渡金属钇（Y）修饰对石墨烯储氢性能的影响。考虑Y原子在石墨烯上易形成团簇,采用B原子掺杂有效阻止了团簇形成。通过模拟计算得到的改性体系稳定、储氢性能优异,可吸附6个H₂分子,平均吸附能范围为-0.539到-0.655 eV （per H₂）,理论上满足理想的氢吸附能范围。经Bader电荷初步计算和基于Y/B/graphene （G）体系吸附H₂分子的电子态密度及电荷差分密度图分析得,Y原子与石墨烯间通过电荷转移产生结合,与H₂分子则发生典型的Kubas型相互作用。Y原子改变了H₂分子与石墨烯基的电荷分布,成为连接两者电子云的桥梁,从而增强了H₂分子的吸附能。改性石墨烯体系吸附的均为氢分子,有利于在环境温度和压力条件下进行循环控制,是具有良好发展前景的储氢材料之一。