An effective method for accurate prediction of the first hyperpolarizability of alkalides

The proper theoretical calculation method for nonlinear optical (NLO) properties is a key factor to design the excellent NLO materials. Yet it is a difficult task to obatin the accurate NLO property of large scale molecule. In present work, an effective intelligent computing method, as called extreme learning machine-neural network (ELM-NN), is proposed to predict accurately the first hyperpolarizability (β(0)) of alkalides from low-accuracy first hyperpolarizability. Compared with neural network (NN) and genetic algorithm neural network (GANN), the root-mean-square deviations of the predicted values obtained by ELM-NN, GANN, and NN with their MP2 counterpart are 0.02, 0.08, and 0.17 a.u., respectively. It suggests that the predicted values obtained by ELM-NN are more accurate than those calculated by NN and GANN methods. Another excellent point of ELM-NN is the ability to obtain the high accuracy level calculated values with less computing cost. Experimental results show that the computing time of MP2 is 2.4-4 times of the computing time of ELM-NN. Thus, the proposed method is a potentially powerful tool in computational chemistry, and it may predict β(0) of the large scale molecules, which is difficult to obtain by high-accuracy theoretical method due to dramatic increasing computational cost.     

6,12-二乙炔基茚并[1,2-b]芴系列衍生物的非线性光学性质

采用密度泛函理论(DFT) CAM-B3LYP方法对6,12-二乙炔基茚并[1,2-b]芴系列衍生物的极化率(α_s)和第二超极化率(γ_s)进行研究. 结果表明, 此类分子具有较大的γ_s值. 用乙炔基硅烷基和氧原子取代茚并[1,2-b]芴分子6,12位的氢原子后, 分子的几何构型发生改变, 进而影响其非线性光学(NLO)性质. 连接乙炔基硅烷基的分子α_s值和γ_s值均增大, 而连接氧原子的分子α_s值和γ_s值均减小. 茚并[1,2-b]芴环2,8位取代基R(R=H, F, CH₃)的不同, 对分子的γ_s值也有一定的影响, R为CH₃时分子的α_s值和γ_s值均较大. 由含时密度泛函理论(TD-DFT)方法计算的吸收光谱分析可知, 与茚并[1,2-b]芴系列分子相比, 引入乙炔基硅烷基的分子共轭性增强, 最大吸收波长红移; 引入氧原子的分子几何结构扭曲, 共轭性降低, 最大吸收波长蓝移.     

表面活性剂在溶液中聚集形态的动力学模拟

用耗散颗粒动力学模拟方法(DPD)展示了表面活性剂分子在溶液中的聚集形态，用扩散程度表征了表面活性剂溶液中的自组装情况。结果发现：这种分子动力学模拟方法能够直观地得到表面活性剂的聚集形态；随着表面活性剂的浓度增加，聚集形态依次从球状胶束、棒状或虫状胶束，六角状相，向层状相变化。     

Quantum Chemistry Studies on the Fe-Cu Interactions and 31p NMR in Fe（CO）3（Ph2Ppy）2（CuXn） （Xn = Cl2^2-, Cl-, Br-）

In order to study the Fe-Cu interactions and their effects on 31p NMR, the structures of mononuclear complex Fe（CO）3fPhzPpy）a 1 and binuclear complexes Fe（CO）3（PhEPpy）z（CuXn） （2： Xn = Cl2^2-, 3： Xn = Cl-, 4： Xn = Br-） are calculated by density functional theory （DFT） PBE0 method. For complexes 1, 3 and 4, the 31p NMR chemical shifts calculated by PBE0-GIAO method are in good agreement with experimental results. The 31p chemical shift is 82.10 ppm in the designed complex 2. The Fe-Cu interactions （including Fe→Cu and Fe←Cu charge transfer） mainly exhibit the indirect interactions. Moreover, the Fe-Cu（I） interactions （mostly acting as σFe-p→4Scu and aFe-C→4Scu charge transfer） in complexes 3 and 4 are stronger than Fe-Cu（Ⅱ） interactions （mostly acting as σFe-p→4Scu and σFe-p←4Sc,） in complex 2. In complex 2, the stronger Fe←Cu interac- tions, acting as σFe-p←44SCu charge transfer, increase the electron density on P nucleus, which causes the upfield 31p chemical shift compared with mononuclear complex 1. For 3 and 4, although a little deshielding for P nucleus is derived from the delocalization of σFe-p→4Scu due to the Fe→Cu interactions, the stronger σFe-c→np charge-transfer finally increases the electron density on P nucleus. As a result, an upfield 31p chemical shift is observed compared with 1. The stability follows the order of 2〉3=4, indicating that Fe（CO）3（PhzPpy）2（CuCl2） is stable and could be synthesized experimentally. The N-Cu（Ⅱ） interaction plays an important role in the stability of 2. Because the delocalization of σFe-p→4SCu and σFe-c→πc-o weakens the a bonds of Fe-C and ~r bonds of CO, it is favorable for increasing the catalytic activity of binuclear complexes. Complexes 3 and 4 are expected to show higher catalytic activity compared to 2.     

8-羟基喹啉铝光电性质的Ab initio和DFT研究

利用abinitioHF和密度泛函理论B3LYP等方法,对金属有机配合物8-羟基喹啉铝(AlQ₃)进行几何结构优化,探索分子内部电子跃迁的机理.结果表明,电子从基态跃迁到低激发态时主要为π-π^*跃迁;电荷从含氧的苯酚环转移至含氮的吡啶环上,包括两环之间C→C转移和O→N转移,与金属离子关系不大.考虑到配体对发光性质的贡献,进一步设计了3种AlQ₃的衍生物.     

聚合物PVP与表面活性剂AOT相互作用的介观模拟

用耗散颗粒动力学模拟(DPD)方法研究了聚乙烯吡咯烷酮(PVP)与2-乙基己基琥珀酸酯磺酸钠(AOT)之间的相互作用.在三维模拟格子中,聚合物链均方末端距〈r²〉随着表面活性剂浓度的增加呈现一种首先减小,接着增加,然后又减小的趋势.构型和结构分析表明,AOT的加入能够引起聚合物链的二面角分布发生改变,这意味着AOT与PVP产生了相互作用.同时表面活性剂/聚合物体系的聚集形态也可以在DPD三维模拟格子中直观显现出来.     

咪唑氧自由基分子NLO性质的密度泛函理论研究

采用密度泛函理论(DFT) UB3LYP方法, 在6-31g(d)水平上对2,2’-(1,2-乙炔基-4,1-亚苯基)双[4,4’,5,5’-四氢]咪唑氧自由基分子及其异构体的自旋耦合性质进行分析, 并结合有限场(FF)方法计算它们的非线性光学(NLO)系数, 以探讨咪唑氧环在共轭链不同位置时体系的自旋耦合规律和NLO系数. 结果表明, 所有体系基态自旋符合自旋极化规则, 它们的极化率随自旋多重度的增加而减小; 一阶超极化率因受分子对称性影响, 对称性不同其一阶超极化率的变化也不同; 二阶超极化率呈现随着自旋多重度的增加而增加的趋势. 从理论上探讨这些自由基分子自旋耦合规律与NLO活性的关系, 为有机自由基NLO材料的分子设计与实验研究提供一定的理论依据.     

气液界面上阴离子表面活性剂单层膜的分子动力学模拟

用分子动力学方法研究了阴离子表面活性剂十二烷基硫酸钠(SDS)在气液界面上的结构和动力学性质. 选择单分子占有面积分别为0.45和0.68 nm²的两个模拟体系, 通过径向分布函数表征了单层膜的厚度, 并根据疏水链中碳原子与极性头中硫原子之间组成的矢量分布和取向函数, 对比了不同界面单层膜的有序排列情况. 结果表明在分子占有面积较小达到饱和吸附的情况下, 界面上的SDS具有较好的有序性. 通过计算气液界面附近水分子的扩散系数发现: 由于氢键和静电作用的影响, 界面区域内的水分子较本体溶液中的水分子有较弱的迁移能力.     

HF分子电子基态和激发态的MRCI计算

应用舍Davidson修正的多参考组态相互作用（MRCI）方法,在aug-cc-pVTZ基组水平上对HF基态及最低的多个单重和三重电子激发态进行了势能扫描计算.结合群论原理及分子的离解极限,分析了电子态势能曲线的特征,得出激发态B¹S⁺对应的离解极限为H⁺+F^-（¹S）.基于势能曲线,数值求解核运动的径向Schrodinge方程,得到J=0时束缚电子态X¹S⁺,B¹S⁺,C¹P和D¹S⁺的振动能级和转动常数,继而进行数据拟合得到电子态的光谱常数,基态X¹S⁺:ω_e=4146.94 cm^-1,ω_ez_e=88.08 cm^-1,B_e=21.22 cm^-1,a=0.785 cm^-1;B¹S⁺态:ω_e=1131.37 cm^-1,ω_ex_e=17.28 cm^-1,B_e=3.96 cm^-1,a_e=0.0215 cm^-1,C¹P态:ω_e=2696.37 cm^-1,ω_ex_e=73.43 cm^-1,B_e=15.91 cm^-1,a_e=0.776 cm^-1,D¹S⁺态:ω_e=3104.22 cm^-1,ω_ex_e=118.92 cm^-1,B_e=17.25 cm^-1,a_e=0.992 cm^-1,拟合结果与实验值吻合的较好.     

NaH和AlH分子基态(X~1Σ~+)的结构和解析势能函数

根据群论及原子分子反应静力学的有关原理,推导出NaH和AlH分子基态(X~1Σ~+)的合理离解极限.使用密度泛函方法中的B3LYP、B3PW91和MPW1PW91方法,在6-311++G,6-311++G(3df,3pd),cc-pVQZ和aug-cc-pVQZ基组下对NaH和AlH分子的基态进行结构优化计算,使用优选的B3PW91/6-311++G(3df,3pd)对基态单点能扫描计算,然后用最小二乘法拟合Murrell-Sorbie函数,得到对应的势能函数参数及光谱常数.结果表明,采用Murrell-Sorbie函数计算所得的光谱常数与实验结果符合的很好,能精确地描述NaH和AlH分子基态的势能函数.