分子动力学模拟压水反应堆中氢气对水的影响

通过分子动力学方法模拟了在常温常压下(1 atm, 298 K)和在压水堆环境下(155 atm, 626 K), 水分子数为256, 氢分子数为0, 25, 50, 75和100等不同数目时, 粒子系统的动力学性质和微观结构, 分析了不同氢气对水中溶解氧的影响. 从模拟结果可知, 在常温常压和压水堆环境下, 当氢粒子数分别为0, 25, 50, 75和100时, 粒子系统的均方位移会随氢分子数增加而增加, 并且常温常压下的增长幅度远小于压水堆环境下的增长幅度, 如压水堆环境下氢分子数为75时系统的均方位移约是常温常压下氢分子数为75时系统的均方位移的6.02倍, 比压水堆环境下氢分子数0时系统的均方位移增加了131.88%. 此外, 粒子系统的微观结构, 从径向分布函数看, 在常温常压下随着氢分子数目的增加而小幅度增加, 这与常温常压下因氢气溶解在水中增大了氧离子周围的粒子密度相符合. 而在压水堆环境下, 氢分子数为75, 50, 25与为0时的水比较, 其径向分布均不会有太大的变化, 而分子数为100时会出现明显增加, 与为0时的水比较其径向分布增加了22.00%. 模拟结果表明, 往压水堆中的水加入氢气能明显地抑制水中的溶解氧.     

低温促进表面等离激元共振效应及肌酐的超灵敏表面增强拉曼散射探测

肌酐是肾脏疾病诊断和监测的关键生物标志物,因此,快速、灵敏的肌酐检测非常重要.本文提供了一种通过提高低温下的光子诱导电荷转移效率来促进表面增强拉曼散射光谱(SERS)活性的有效策略.采用种子生长法获得纳米金二十面体(Au₂₀),以此作为SERS活性基底.采用极低温(98 K)SERS检测技术实现对染料分子结晶紫(CV)和生理盐水中的肌酐含量的快速、灵敏检测.实验结果表明:常温296 K下, Au₂₀基底对CV分子的检测限(LOD)低至10^–12 mol/L,并且信号均匀;低温98 K时, CV分子的LOD可达10^–14 mol/L,比296 K时降低2个数量级.最后,使用Au₂₀基底对生理盐水中的肌酐进行无标记检测.结果表明,常温296 K时该SERS基底对肌酐的LOD为10^–6 mol/L, 1619 cm^–1峰的线性相关系数为0.9839.低温98 K时,肌酐的浓度探测极限低至10^–8 mol/L, 1619...     

TiO2/ZnO纳米薄膜界面热导率的分子动力学模拟

采用平衡分子动力学方法及Buckingham势研究了金红石型TiO₂薄膜与闪锌矿型ZnO薄膜构筑的纳米薄膜界面沿晶面[0001](z轴方向)的热导率.通过优化分子模拟初始条件中的截断半径r_c和时间步后,计算并分析了平衡温度、薄膜厚度、薄膜截面大小对热导率的影响.研究表明,薄膜热导率受薄膜温度和厚度的影响很大,当温度由300 K升高600 K时,薄膜的热导率逐渐减小;当薄膜厚度由1.8 nm增大到5 nm时,热导率会逐渐增大;并在此基础 关键词： 热导率 分子动力学 2/ZnO纳米薄膜界面')" href="#">TiO₂/ZnO纳米薄膜界面 数值模拟     

Langmuir-Blodgett膜摩擦分子动力学模拟和机理研究

采用分子动力学方法,模拟了在两块石英基板上由脂肪酸(C₁₅H₃₁COOH)组成的单层Langmuir-Blodgett (LB)膜间的摩擦特性,探究了超薄膜在滑动过程中的摩擦和结构机理.得出对于单层LB膜在滑动过程中,在速度小于60m/s时,随着速度的增大,其剪切压增大;在速度大于60 m/s时,剪切压随速度的增加而减小.其链的倾斜角随着滑动速度的增加而减小.单层膜内的分子之间以氢键方式形成了较大的分子簇,由此导致了剪切压呈现较长的周期性,但在单层膜之间无氢键形 关键词： 分子动力学模拟 朗缪尔布洛杰特膜 纳米摩擦 氢键     

低能Cl原子刻蚀Si(100)表面的分子动力学模拟

使用分子动力学模拟方法研究了不同能量(0.3—10 eV)的Cl原子对表面温度为300 K的Si(100)表面的刻蚀过程.模拟中采用了Tersoff-Brenner势能函数来描述Cl-Si体系的相互作用.模拟结果显示,随着入射Cl原子在表面的吸附达到饱和,Si表面形成一层富Cl反应层.这和实验结果是一致的.反应层厚度随入射能量增加而增加.反应层中主要化合物类型为SiCl,且主要分布于反应层底部.模拟结果发现随初始入射能量的增加,Si的刻蚀率增大.在入射能量为0.3,1和5 eV时,主要的Si刻蚀产物为Si 关键词： 分子动力学 Cl刻蚀Si 分子动力学模拟 微电子机械系统     

CaF2熔化温度的分子动力学模拟

利用壳层模型分子动力学方法,研究了高温高压条件下CaF2的熔化温度,同时计算了温度为300K、压强上升到100GPa时CaF2 的状态方程.研究中考虑了分子动力学模拟的过热熔化,通过晶体的现代熔化理论,对CaF2 的分子动力学模拟熔化温度进行了修正, 获得了高温高压下CaF2的熔化温度.因此,常压下壳层模型分子动力学方法为研究物质熔化提供了一个很好的方法.     

超临界二氧化碳类液-类气区边界线数值分析

超临界二氧化碳(S-CO₂)因在萃取、沉淀、热力循环及化学反应等方面有着十分广阔的应用前景,逐渐成为学术界的重要研究课题.由于在近临界区,可以观察到随温度或压力变化出现大量的物性异变现象,使得各国学者对流体临界点附近区域的研究产生了浓厚兴趣.随着分子动力学模拟技术的快速发展,该技术可辅助传统实验方法用于研究近临界流体的相关物性.为确定S-CO₂在近临界区Widom线范围及类液-类气区的分子结构特征,本文通过分子动力学模拟技术结合聚类分析,研究了温度和压力范围分别在300—350 K和5.5—18.5 MPa下,CO₂密度时间序列变异系数及偏度同Widom线和类液-类气区间的关系.结果表明:S-CO₂在近临界区Widom线的确定可通过连接密度时间序列曲线变异系数极大值点来确定,Widom线沿着临界点开始延伸直到350 K时停止;S-CO₂类液区和类气区的分子分布结构可以用数密度分布的偏度来区分,偏度在类气态时为正值,在类液态时为负值,而在Widom线上达到最大值.     

超临界二氧化碳类液-类气区边界线数值分析

超临界二氧化碳(S-CO₂)因在萃取、沉淀、热力循环及化学反应等方面有着十分广阔的应用前景,逐渐成为学术界的重要研究课题.由于在近临界区,可以观察到随温度或压力变化出现大量的物性异变现象,使得各国学者对流体临界点附近区域的研究产生了浓厚兴趣.随着分子动力学模拟技术的快速发展,该技术可辅助传统实验方法用于研究近临界流体的相关物性.为确定S-CO₂在近临界区Widom线范围及类液-类气区的分子结构特征,本文通过分子动力学模拟技术结合聚类分析,研究了温度和压力范围分别在300—350 K和5.5—18.5 MPa下,CO₂密度时间序列变异系数及偏度同Widom线和类液-类气区间的关系.结果表明:S-CO₂在近临界区Widom线的确定可通过连接密度时间序列曲线变异系数极大值点来确定,Widom线沿着临界点开始延伸直到350 K时停止;S-CO₂类液区和类气区的分子分布结构可以用数密度分布的偏度来区分,偏度在类气态时为正值,在类液态时为负值,而在Widom线上达到最大值.     

纳米摩擦中极性有机分子超薄膜的结构、对称性及能量机理

采用分子动力学方法, 模拟了由脂肪酸C_nH_2n+1COOH和C₁₇H₃₅COOH (n=12,13,14,15,16,17) 按1:1比例组成的7种混合单层Langmuir-Blodgett (LB)膜和由C₁₆H₃₃COOH 分子组成的单层膜的摩擦性质. 结果显示: 1) 随着混合单层膜内的不同分子链长差的减小, 其剪切压随之减小, 摩擦力主要来自单层膜间的库伦作用; 2) 混合膜内的两种不同分子的尾基排列对其摩擦性能影响较大, 当混合LB膜中所有分子尾基全同排列时剪切压较小. 当分子链长差为1 个C-C键长时, 分子尾基排列对膜的摩擦性质影响较大. 3) 同种分子尾基全同排列组成的单层膜, 当上下两单层膜的尾基呈镜面对称时, 其剪切压随着分子链长的增加而减小, 摩擦力主要来自膜间的库伦作用; 当上下两单层膜的尾基呈中心对称时, 膜间摩擦力主要来自膜间的范德华 (VDW) 作用. 关键词： 分子动力学模拟 纳米摩擦 薄膜 库伦能     

Ni熔体凝固过程中临界晶核和亚临界晶核的分子动力学模拟

本文用分子动力学模拟研究了Ni熔体以不同冷速凝固后微观结构的演变规律, 并通过理论计算确定出了Ni熔体凝固后获得理想非晶的临界条件. 模拟结果发现冷速小于10¹¹ K/s时, Ni 熔体凝固后形成晶态组织; 冷速在10¹¹ K/s到10^14.5 K/s之间时, Ni熔体凝固后形成由晶态结构与非晶态结构组成的混合组织. 冷速小于10¹⁰ K/s, Ni 熔体凝固后形成的晶态组织具有fcc结构; 冷速在10¹⁰ K/s到10^14.5 K/s之间时, Ni熔体凝固后组织中的晶态由fcc和hcp结构层状镶嵌排列构成. 通过分析模拟结果和计算结果, 确定出了Ni熔体凝固后形成理想非晶的临界冷速为10^14.5 K/s. 并发现Ni熔体中临界晶核(冷速等于10^14.5 K/s)和亚临界晶核(冷速大于10^14.5 K/s) 均由fcc和hcp组成层状偏聚结构, 这表明Ni熔体中生长的晶体、临界晶核和晶胚的结构是相同的. 关键词： 分子动力学模拟 晶体团簇 临界冷速 结构     

Clusterization of water molecules as deduced from statistical mechanical approach

Using the methods of statistical mechanics we have shown that a homogeneous water network is unstable and spontaneously disintegrates to the nonhomogeneous state (i.e. peculiar clusters), which can be treated as an ordinary state of liquid water. The major peculiarity of the concept is that it separates the paired potential into two independent components—the attractive potential and the repulsive one, which in turn should feature a very different dependence on the distance from the particle (a water molecule in the present case). We choose the interaction potential as a combination of the ionic crystal potential and the vibratory potential associated with the elastic properties of the water system as a whole. The number ℵ of water molecules that enters a cluster is calculated as a function of several parameters, such as the dielectric constant, the mass of a water molecule, the distance between nearest molecules, and the vibrations of nearest molecules in their nodes. The number of H₂O molecules that comprise a cluster is estimated as about ℵ ≈ 900, which agrees with the available experimental data. Presented at the 2nd International Conference “Physics of Liquid Matter: Modern Problems” (September 2003, Kyiv, Ukraine)     

Interaction potential and selection rules in the rotational spectrum of water molecule during its collision with atoms

If water molecules are in an atmosphere of buffer monatomic gas, so that the collision frequency with atoms is much higher than the collision frequency between molecules, rotational levels of the molecule with the same angular momentum quantum number j are split into two sets with even and odd quantum numbers m′ (projection on the Oz′ axis of the molecular coordinate system). Such separation arises due to the interaction of the dipole moment of water molecule with the induced dipole on an atom during particle collision. The symmetry of the molecule-atom interaction potential allows transitions between levels of the same parity, but forbids transitions between sets. Experiments on saturation of rotational transitions in a rarefied gas mixture consisting of heavy water and argon vapors are interpreted. These experiments confirm the assumption that rotational levels of the water molecule are indeed separated into two sets due to the interaction of the molecule with the environment.     

绝缘油在水和酸共同作用下的分子动力学模拟

为探讨水分和温度对变压器油中酸扩散行为的影响,在分子模拟软件内建立了水分子、矿物油分子、甲酸分子的复合模型,水的体积分数分别为1%,3%,5%,甲酸的质量分数为3%.通过分子动力学方法模拟计算了甲酸分子在不同含水量油模型中的扩散系数、径向分布函数以及复合模型的总动能;此外,选取了水体积分数为1%的复合模型为研究对象,以温度为变量进行动力学模拟.结果表明:随着水含量的增多,甲酸的扩散系数并没有逐渐增大,甲酸分子与水分子之间的氢键数目呈逐渐减少的趋势.而随着温度的升高,油模型具有的动能越大,酸分子的扩散逐渐增强.油模型中,含水量越多,水和酸在热运动下碰撞油分子的能量越大;温度越高,相互之间的碰撞越激烈.模拟结果说明水分增加或温度增高,均对油分子的结构稳定性造成一定程度的破坏.     

Ab initio molecular dynamics investigations of structural, electronic and dynamical properties of water in supercritical carbon dioxide

Ab initio molecular dynamics simulations of a solitary perdeuterated water molecule solvated in supercritical carbon dioxide have been performed along an isotherm at three different densities. Electron donor-acceptor interactions between the oxygen atom of water and the carbon atom of CO₂ as well as hydrogen bonded interactions between the two molecules have been shown to play a dominant role in the solvation. The mean dipole moment of the water molecule increases with the density of the solution, from a value of 1.85 D at low density to around 2.15 D at the highest density. The increase in the solvent density causes the water molecule to exhibit a range of behavior, from a free molecule to one that interacts strongly with CO₂. A blue shift in the bending mode of water has been observed with increasing solvent density. The carbon dioxide molecules which are present in the first neighbor shell of water are found to exhibit larger propensity to deviate from a linear geometry in their instantaneous configurations.      

Experimental analysis of the electrode products emitted byhigh-current vacuum arcs

The mass distribution of particles produced in the high-current vacuum arc was investigated. The experiments were concentrated on evaluating the spatial mass distribution emitted in the radial as well as in the azimuthal directions calculated from the mass deposition profile on collectors surrounding the arc discharge. The experiments were carried out in a vacuum chamber evacuated to an ambient pressure <10^-4 Pa. High-current arcs in the range from 2-7.5 kA were drawn between butt contacts of 31 and 55 mm in diameter (anode and cathode, respectively) both of a copper-chromium alloy (CuCr25). The surface mass deposited along multi-segment collectors was measured by a micro densitometer, and an inductively coupled plasma (ICP) spectrometer. Two angular mass deposit distributions were determined: the azimuthal distribution on the plane parallel to the cathode surface, and the radial distribution as a function of the angle with respect to the cathode plane. Both distributions were anisotropic and the structure of the deposition layer depended on the angle of incidence of the particles onto the substrate, the density of the particle flux and other factors. The mass deposited on the collectors consisted mostly of chromium molecules (approx. 80% of Cr and 20% Cu) for CuCr25 or CuCr40 electrodes     

Molecular dynamic simulation of core–shell structure: study of the interaction between modified surface of nano-SiO2 and PAMAA in vacuum and aqueous solution

The interaction between acrylamide acrylicacid copolymer (PAMAA) and the modified surface of nano-SiO₂ is investigated using the molecular dynamic (MD) simulation. The binding energies (E_binding) of interface, the concentration profiles of PAMAA and functional groups (carboxyl and acylamino) of corresponding model, the mean square displacements (MSD) and diffusion coefficients (D) of PAMAA in four systems with different modifiers are all calculated at 325 K in vacuum. Vinyl trimethoxy silane (VTEOS) shows best modification effect in the systems mentioned above. Furthermore, the effects of temperature on the interaction between VTEOS modified surface of nano-SiO₂ and PAMAA are studied at 300, 325, 350, 375 and 400 K in aqueous solution. Interesting results show that, water molecular layer reduces with the increase of temperature, and then improves the interaction between PAMAA and VTEOS modified surface of nano-SiO₂. The corresponding E_binding of interface, the radial distribution functions (RDF) of carbon atoms on the surface and oxygen atoms of water molecules, the concentration profiles of PAMAA on the surface of nano-SiO₂, the MSD and D of PAMAA are all studied seriously to find the reason of this counterintuitive phenomenon.     

大型先进压水堆中次锕系核素嬗变特性

随着核电事业的快速发展,核电厂卸载的乏燃料越来越多。如何处置核电站乏燃料中的次锕系核素（MA）既是核燃料再利用的重要过程,又是闭式循环中的关键步骤。如果处置得当,不仅可以提高燃料的利用率,而且可以将MA变成同位素燃料电池、中子源等有用的核素。国际上认可的处置方法是分离-嬗变,但是嬗变MA的难点是嬗变堆型的选取和如何提高嬗变率。压水堆(PWR)是国内外最成熟的堆型和商业运行的主要堆型,也是现阶段最具有可能进行MA嬗变的堆型。于是,本文利用MCNP程序研究了压水堆嬗变MA的特性,通过研究MA嬗变棒的设计、添加位置和添加量等对压水堆堆芯有效增殖因子的影响,初步探索出最佳的压水堆嬗变MA的设计方案,为我国现阶段进行压水堆嬗变MA奠定了理论基础。     

Structure correlations of water molecules in a concentrated aqueous solution of ethanol

Monte Carlo calculation results for water–ethanol mixtures calculated in the isothermal and isobaric ensemble at T=300 K are presented. The analysis of the radial distribution functions for hydrogen bonding between ethanol molecule and water molecules was done. The concentrations of ethanol 0.75 show characteristics features of hydrogen-bonded liquids. The influence of the polar solvent molecule on water local structure and the hydrogen bond nets ware analyzed. Radial distribution functions for water–water interactions calculated in pure water and water–ethanol system are compared.     

Role of oxygen functional groups for structure and dynamics of interfacial water on low rank coal surface: a molecular dynamics simulation

Molecular dynamics simulations were employed to study the effects of oxygen functional groups for structure and dynamics properties of interfacial water molecules on the subbituminous coal surface. Because of complex composition and structure, the graphite surface modified by hydroxyl, carboxyl and carbonyl groups was used to represent the surface model of subbituminous coal according to XPS results, and the composing proportion for hydroxyl, carbonyl and carboxyl is 25:3:5. The hydration energy with ?386.28 kJ/mol means that the adsorption process between water and coal surface is spontaneous. Density profiles for oxygen atoms and hydrogen atoms indicate that the coal surface properties affect the structural and dynamic characteristics of the interfacial water molecules. The interfacial water exhibits much more ordering than bulk water. The results of radial distribution functions, mean square displacement and local self-diffusion coefficient for water molecule related to three oxygen moieties confirmed that the water molecules prefer to absorb with carboxylic groups, and adsorption of water molecules at the hydroxyl and carbonyl is similar.