动量表象中的氦原子和氢分子

本文推导出动量表象中的基态氦原子和氢分子波函数。对午氦原子,使用了动量变量P_r,P_θ,P_φ,它们分别共轭于球坐标位置变量r,θ,φ。对于氢分子,使用了动量变量P_ξ,P_η,P_φ,它们分别共轭于椭球坐标位置变量ξ,η,φ。     

2,5-二取代噁唑的研究(Ⅲ)——2-联苯基-5-苯基噁唑和其5-对位取代苯基衍生物的合成

合成了2-联苯基-5-苯基噁唑和13种5-对位取代苯基衍生物,测定了它们的熔点,红外紫外吸收光谱,萤光发射光谱,萤光量子产率(φ_i)和激光转换效率(η_i)。     

炭黑填充聚苯乙烯熔体流变行为

研究了炭黑(CB)填充聚苯乙烯(PS)熔体的稳态和动态流变行为. CB/PS复合体系在CB体积分数φ=0.06时发生逾渗转变. 结果表明, 低应变区熔体模量降低主要归因于粒子-粒子及粒子-高分子间作用力的破坏, 高应变下模量的急剧下降则主要与高分子链间解缠结有关. 采用“两相”模型拟合线性动态流变行为, 发现应变放大因子A_f(φ)、填充相模量及松弛指数与温度有关. A_f(φ)~φ关系符合Guth方程和扩散控制的粒子簇聚集模型. “粒子相”形状参数与聚集体分维度均随温度升高而有所降低, 说明CB粒子聚集体因团聚而趋于各向同性, 应变放大效应减弱. “粒子相”特征模量G'_f1(φ)和G"_f0(φ)与φ关系满足标度律. 当φ > 0.06时, G'_f1(φ)和G"_f0(φ)及其标度指数均随温度升高而明显降低, 其G'_f1(φ)变化幅度略大于G"_f0(φ), 说明“粒子相”弹性与黏性组分具有不同的温度依赖性. 随着温度升高, 扩散控制的CB粒子团聚过程加快, 应变放大效应减弱.     

丙烯腈在Cu(100)表面化学吸附的密度泛函理论研究

利用密度泛函方法, 模拟金属铜原子簇Cu₁₄(9,4,1)的(100)表面, 对丙烯腈(CH₂＝CHCN)在Cu(100)面上不同吸附位的吸附状况进行了理论研究. 结果表明: 丙烯腈分子通过端位N原子垂直吸附在金属表面上为弱化学吸附, 部分电荷从丙烯腈分子转移至铜金属簇; 由N原子的孤对电子与金属铜形成弱σ共价键; 顶位是最佳吸附位, 吸附能为40.7391 kJ•mol^－1, N原子与金属表面间的平衡距离为0.2279 nm; 其次为桥位和穴位, 吸附能分别为36.2513和30.2158 kJ•mol^－1, 平衡距离为0.2194和0.2886 nm; 吸附后C≡N键的强度降低, 活化了丙烯腈分子. 化学吸附使体系的熵减小, 是由于丙烯腈分子的平动和转动因吸附而被限制.     

RDF法对非晶态聚醚砜共聚物的近程有序结构研究

用径向分布函数(RDF)研究了聚醚砜共聚物(PES)在r=1.5 nm范围内的近程结构, 结果表明, 分子链间存在近程有序结构, 有序周期约为0.5 nm; 聚醚砜共聚物经退火处理后, 由g(r)图可以看出, 第一个峰变高、变窄; 由G(r)图可以看出, 分子链间第一个有序区内, 分子链间贡献大于链内贡献; 于200 ℃退火处理后有序畴尺寸(rs)略有增大, 原子平均位移(σ)变小, 使得分子链内原子间距分布更加均匀, 相邻分子链间距(rv)变小. 同时透射电镜(TEM)的形态结构分析表明, 经退火处理后, 其聚集态确实形成了某种局部有序的结构.     

采用量子力学方法研究带正电氨基酸分子复合物间非键作用的强度

非键作用是自然界中广泛存在的分子间作用,参与到生命体中最基本的蛋白质与核酸的相互作用等过程。该过程以氨基酸侧链与碱基的相互作用为主,因此准确预测氨基酸侧链与碱基等生物分子的相互作用具有十分重要的科学意义。选用量子力学方法计算了15个含氨基酸侧链的带正电堆积复合物的相互作用能,并与CCSD(T)/CBS方法计算结果进行比较。计算结果表明,SCS-CCSD/CBS方法的均方根偏差为0.06 kcal/mol,线性相关系数为0.999 8;SCS-MI-CCSD/CBS方法的均方根偏差为0.11 kcal/mol,线性相关系数为0.999 9;MP2.5/CBS方法的均方根偏差为0.13 kcal/mol,线性相关系数为0.999 3,验证了上述方法的可靠性。而后进一步使用MP2.5/CBS方法计算精氨酸、组氨酸和赖氨酸侧链与中性分子形成堆积复合物在不同距离下的相互作用能,进而构建并完善了带正电氨基酸分子与中性分子间非键作用强度的标准数据。     

C60与铈卟啉相互作用的密度泛函理论研究

应用Perdew-Burke-Ernzerhof (PBE)密度泛函理论对5种可能存在的富勒烯与铈卟啉复合物进行了几何结构优化, 通过分子间距离及结合能数值, 确认了C₆₀和铈卟啉可以通过非键相互作用形成超分子主-客体复合物, 且最有可能的作用位点为C₆₀的C5:6键(相邻五元环与六元环共用碳-碳键). 应用扩展过渡态方法对结合能进行分解, 分解结果显示, 静电能对总吸引能的贡献约为60%, 说明静电作用是复合物稳定存在的最主要因素. 最后对复合物中的电子流向进行了研究, 结果表明复合物中电子的转移与主-客体间相互作用有一定联系, 且电子是从主体铈卟啉流向客体C₆₀.     

H2分子在Mg3N2表面吸附的第一性原理研究

采用第一性原理方法,通过计算表面能确定Mg₃N₂（011）为最稳定的吸附表面,分别研究了H₂分子在Mg₃N₂（011）三种终止表面的吸附性质.研究发现H₂分子平行表面放置更有利于吸附,表面能最低的终止表面Model Ⅱ上吸附H₂分子最稳定,主要存在三种化学吸附方式：第一种吸附方式,H₂分子解离成2个H原子分别吸附在N原子上形成双NH基,这是最佳吸附方式;此时H₂分子与Mg₃N₂表面间主要是H原子的1s轨道和N原子的2s、2p轨道发生作用,N-H之间为典型的共价键.第二种吸附方式中H₂分子部分解离,两个H原子吸附在同一个N原子上形成NH₂基.第三种吸附方式中H₂分子解离成两个H原子,一个H原子和表面N原子作用形成NH基,另一个H原子和表面Mg原子作用形成MgH结构.三种吸附方式不存在竞争关系,形成双NH基的吸附方式反应能垒最低,最容易发生.除此之外H₂还能以分子的形式吸附在晶体表面,形成物理吸附.     

O+HCl→OH+Cl反应的立体动力学的准经典轨线研究

运用准经典轨线方法, 基于Peterson从头计算势能面对O+HCl→OH+Cl反应的立体动力学性质进行了研究. 讨论了在31.77和51.04 kJ/mol两种碰撞能情况下极化依赖的微分反应截面(2π/σ)(dσ₀₀/dω_t), (2π/σ)(dσ₂₀/dω_t), (2π/σ)(dσ₂₂₊/dω_t)和(2π/σ)(dσ_21-/dω_t)以及描述k-j′两矢量相关和k-k′-j′三矢量相关的分布函数P(θ_r)和P(φ_r). 计算得到的P(θ_r)分布表明, 产物分子的转动角动量j′具有强烈的取向分布, 并且产物转动角动量的取向效应对散射能的变化比较敏感. 而P(φ_r)的分布表明, 产物分子虽然具有沿着y轴的取向效应, 但是没有明显的定向效应.     

手性有机化合物构效关系的量子拓扑学研究

在拓扑化学理论基础上, 从手性分子中成键原子的结构特征和手性原子所处的局部化学微环境出发, 应用密度泛函理论(DFT), 在B3LYP/6-31＋G(d)水平上计算所得距离代替传统二维距离, 结合分子中各原子的支化度, 应用原子的平衡电负性对分子图进行着色, 并引入手性校正因子进行校正, 得到新型的手性拓扑指数w₁和w₂. 采用偏最小二乘回归法建立了18种手性羟基酸和氨基酸的薄层色谱保留指数R_M的定量结构-保留拓扑模型, 并用这种模型对R_M进行预测, 结果表明预测结果和实验值吻合较好. 同时对拓扑模型采用留一交叉检验法(leave-one-out cross-validation, LOO-CV)和外检验相结合的方法进行测试, 测试结果显示模型具有良好稳定性和较强的预测能力.     

Syntheses and structures of dimesitylbismuth(III) bromide, Mes2BiBr, and bis(diphenyldithiophosphinato)mesitylbismuth(III), MesBi(S2PPh2)2

Addition of BiBr₃ to Mes₃Bi (Mes = 2, 4, 6-Me₃C₆H₂) in Et₂O gives 86% of Mes₂BiBr (1) as yellow crystals. Reaction of 1 with Ph₂PS₂NH₄ in a 1 : 1 molar ratio gives a quantitative yield of MesBi(S₂PPh₂)₂ (2) rather than the expected dimesitylbismuth compound. The crystal and molecular structures of 1 and 2 were determined at 153 K and 173 K, respectively. They contain Mes₂BiBr molecules with trigonal pyramidal coordination around Bi. The mean Bi---C bond distance is 2.27 Å and the Bi---Br bond distance is 2.690(2) Å. The angles around Bi vary between 89.4 and 106.4°. Intermolecular Bi…Br contacts of 3.795 Å, indicating weak secondary bonding, give rise to zig-zag shaped (Bi---Br)_x chains. In the polymeric chain the coordination geometry around bismuth atoms can be described as pseudo-trigonal bipyramidal. The crystals of 2 consist of discrete monomeric MesBi(S₂PPh₂)₂ molecules with a symmetry plane containing the metal atom and the aromatic ring of the attached mesityl group. The dithiophosphinato ligands exhibit an anisobidentate coordination pattern with long and short phosphorus—sulfur bonds, i.e. P(1)---S(1) 2.051(31) Å and P(1)---S(2) 1.980(3) Å, related to short and long bismuth—sulfur distances, respectively, i.e. Bi---S(1) 2.662(2) Å and Bi---S(2) 3.123(3) Å. This leads to a square-pyramidal geometry around the bismuth atom, with the metal lying 0.33 Å above the basal plane formed by the four sulfur atoms.     

A Bayesian molecular interaction library

We describe a library of molecular fragments designed to model and predict non-bonded interactions between atoms. We apply the Bayesian approach, whereby prior knowledge and uncertainty of the mathematical model are incorporated into the estimated model and its parameters. The molecular interaction data are strengthened by narrowing the atom classification to 14 atom types, focusing on independent molecular contacts that lie within a short cutoff distance, and symmetrizing the interaction data for the molecular fragments. Furthermore, the location of atoms in contact with a molecular fragment are modeled by Gaussian mixture densities whose maximum a posteriori estimates are obtained by applying a version of the expectation-maximization algorithm that incorporates hyperparameters for the components of the Gaussian mixtures. A routine is introduced providing the hyperparameters and the initial values of the parameters of the Gaussian mixture densities. A model selection criterion, based on the concept of a `minimum message length' is used to automatically select the optimal complexity of a mixture model and the most suitable orientation of a reference frame for a fragment in a coordinate system. The type of atom interacting with a molecular fragment is predicted by values of the posterior probability function and the accuracy of these predictions is evaluated by comparing the predicted atom type with the actual atom type seen in crystal structures. The fact that an atom will simultaneously interact with several molecular fragments forming a cohesive network of interactions is exploited by introducing two strategies that combine the predictions of atom types given by multiple fragments. The accuracy of these combined predictions is compared with those based on an individual fragment. Exhaustive validation analyses and qualitative examples (e.g., the ligand-binding domain of glutamate receptors) demonstrate that these improvements lead to effective modeling and prediction of molecular interactions.     

Similarity searching in databases of flexible 3D structures using smoothed bounded distance matrices

This paper describes a method for calculating the similarity between pairs of chemical structures represented by 3D molecular graphs. The method is based on a graph matching procedure that accommodates conformational flexibility by using distance ranges between pairs of atoms, rather than fixing the atom pair distances. These distance ranges are generated using triangle and tetrangle bound smoothing techniques from distance geometry. The effectiveness of the proposed method in retrieving other compounds of like biological activity is evaluated, and the results are compared with those obtained from other, 2D-based methods for similarity searching.     

单占据片断轨道的组合式(Kost+Perkin)定域化

为了使开壳层片断轨道具有正确的集居数,正则轨道必须定域化.但是,在确保目标轨道Φ_sdan单占据性的同时,Kost定域化也破坏了Φ_s特定的对称性.实际计算表明,Perkin定域化可以弥补Kost程序的缺陷,将单占据的Φ_s转化成高度定域的、对称的片断轨道.在片断分子中,C-H_R键长r的选择和Kost定域化的方式对Perkin定域化的成败具有重大的影响.当Gaussian基组为STO-3G3-21G和4-31G时,r应为0.1nm;但在6-31G水平下,必须r=0.09nm.     

Lattice dynamics of solid hydrogen chloride

Lattice dynamics of solid hydrogen chloride is studied, assuming a rigid molecule, in the harmonic and pair potential (Lennard-Jones interaction) approximation between atoms with the inclusion of electrostatic interactions between point dipoles placed on atomic centres. The potential parameters for each of the different non-bonded atom pairs were obtained by means of an optimization routine to give a least-squares fit to observed zone centre (k = 0) frequencies, equilibrium conditions and lattice energy of the lattice.     

Rigid-CLL: avoiding constant-distance computations in cell linked-lists algorithms

Many of the existing molecular simulation tools require the efficient identification of the set of nonbonded interacting atoms. This is necessary, for instance, to compute the energy values or the steric contacts between atoms. Cell linked-lists can be used to determine the pairs of atoms closer than a given cutoff distance in asymptotically optimal time. Despite this long-term optimality, many spurious distances are anyway computed with this method. Therefore, several improvements have been proposed, most of them aiming to refine the volume of influence for each atom. Here, we suggest a different improvement strategy based on avoiding to fill cells with those atoms that are always at a constant distance of a given atom. This technique is particularly effective when large groups of the particles in the simulation behave as rigid bodies as it is the case in simplified models considering only few of the degrees of freedom of the molecule. In these cases, the proposed technique can reduce the number of distance computations by more than one order of magnitude, as compared with the standard cell linked-list technique. The benefits of this technique are obtained without incurring in additional computation costs, because it carries out the same operations as the standard cell linked-list algorithm, although in a different order. Since the focus of the technique is the order of the operations, it might be combined with existing improvements based on bounding the volume of influence for each atom.     

Technetium complexes with thioether ligands—III. Synthesis and structural characterization of cationic nitridotechnetium(V) complexes with thiacrown ethers

The first representative of a new class of TcN complexes with thiacrown ethers have been prepared by ligand exchange reaction of NBu₄[TcNCl₄] with 1,4,8,11-tetrathiacyclotetradecane (14S4), 1,5,9,13-tetrathiacyclohexadecane (16S4), 1,5,9,13-tetrathiacyclohexadecane-3,11-diole (16S4-(OH)₂) and 1,4,7,10,13,16-hexathiacyclooctadecane (18S6). The crystal structure of [TcNCl(14S4)]TcNCl₄ 1) consists of couples of independent cations with the metal in the oxidation state + 5 and hexavalent TcNCl₄⁻ anions. In the complex cation the metal is six-coordinated in a rather distorted octahedral geometry, being directly bound to four sulphur atoms from the macrocyclic ligand in the equatorial plane and to the nitrido atom and to one chlorine atom in the axial positions. The strong trans influence of the nitrido ligand causes an extreme lengthening of the Tc---Cl bond distance to 2.718 Å. The octahedral molecular structure of [TcNCl(18S6)]TcNCl₄ (3) is comparable with that of 1, but only four sulphur atoms of the thiacrown ether form the equatorial plane, two sulphur atoms remain non-coordinated, and the nitrido and Cl⁻ ligands are in axial positions. The most interesting feature in the structure of [TcNCl(16S4-(OH)₂)]Cl (5) is the observation of an exceptionally long Tc---N distance of 1.95 Å.     

定域片断轨道基组的对称化

片断的UHF运算不能保证每个片断轨道具有确切的电子占据数,故Kost定域化是必需的.当片断产生于多键断裂时,在确保目标轨道单占据性的同时,Kost定域破坏了轨道基组的对称性.为此,在Kost定域化后,必须对单占据轨道作2×2对称化旋转后,再作有条件的RHF运算.以乙烯基片断CHCH(波二烯分子中的一个片断)为例,详述了对称化的方法、原理和计算程序.以C-H片断为例,细述三单键片断轨道基组对称化的特殊性.介绍C-HR参考键长选择的判据,探讨键长与选择Gaussian基组大小的关系.     

聚丙烯腈PAN分子链的第一性原理研究

The equilibrium geometries of the polyacrylonitrile (PAN) chain was theoretically studied using the Hartree-Fock method at the STO-3G levels. As for the optimized structures, the average distance of the C atom couple in the main chain is 155.6 pm; the average distance of the C atom couple in the branch chain is 149.7 pm; the average distance between a C atom bonding with N atom is 115.5 pm. For the charge distribution, because of the influence of a N atom with its comparatively larger negative charge, the C atoms in the main chain are different in their charge distribution. Finally, the vibration models of the chain have been analyzed to clarify the reaction sequence of dehydrogenation and cyclization during pre oxidation and carbonization of the polyacrylonitrile.