氯化钾水溶液微观结构的分子动力学模拟研究

氯化钾水溶液不同于氯化钠的,其在低温下只能形成无水盐和冰.现有研究表明在氯化钠水溶液中不仅存在着水分子连续分布的区域与离子团簇,而且这二者的物质组成也与该溶液所形成晶体的分别对应.为了探寻其他溶液中是否存在此类对应关系,并考察溶液的微观结构,本文采用分子动力学方法对氯化钾水溶液进行了研究,表明了该溶液的一些性质.氯化钾水溶液中K⁺-K⁺和Cl^--Cl^-径向分布函数的特征具有一致性,峰的最大值所对应的位置都分别相同,明显不同于氯化钠水溶液的.系列时刻下瞬态图像内O到其最近离子距离中最大值的统计结果表明氯化钾水溶液中存在着一定大小的水分子连续分布的区域,其平均尺寸至少为2.26 nm;瞬态图像中K⁺与其最近邻Cl^-之间的距离主要分布在0.28 nm～0.38 nm之间,占比约为97.4%;溶液中存在着较大和较小两类离子团簇,较大团簇的平均尺寸为1.73nm,平均离子数是25.0,其内部的离子与周围离子之间具有与氯化钾晶体类似的结构;这些结果表明氯化钾水溶液中也...     

氯化钠水溶液结构的研究

利用上海光源(SSRF)的第三代同步辐射光源测定室温下摩尔浓度分别为0.172 mol/L、0.343 mol/L、0.699 mol/L、1.064 mol/L、2.832 mol/L、3.910 mol/L、5.289 mol/L的NaCl水溶液的X射线散射数据。由X射线散射数据可知,随着NaCl水溶液浓度的增大,X射线散射曲线的特征峰由12.6°到13.4°发生偏移。运用Pair Distribution Function(PDF)理论对X射线散射数据进行处理,得到了不同浓度NaCl水溶液及纯水的差值对分布函数,其中的O-O峰随着浓度的增大逐渐分裂为两峰,O-O峰位在0.282 nm处。利用分子动力学模拟研究不同浓度的NaCl水溶液,表明Na+、Cl-的引入对水分子的氢键结构有一定的破坏,当浓度大于15%时,这种效果尤其明显。Na+、Cl-均存在两层水化层,各离子间配位数随浓度的增大而减少。H2O分子的自扩散系数远大于Na+和Cl-的自扩散系数,后两者的值随浓度的增大逐渐减少,Na+和Cl-的水化半径均随浓度增大而降低。     

C36团簇自组装的分子动力学研究

提出了利用C36团簇在气相条件下自组装制备新纳米团簇的设想，并利用分子动力学方法模 拟了包括真实氦气氛作用的碳团簇生长过程，发现环境气体温度是影响最后所生成的团簇结 构的关键因素：C36团簇在1000?K到2000?K的温度范围内，自组装形成保持C36线径特征的 蚕茧状新纳米团簇；在高于2000?K的温度下，最后形成的团簇趋于球状. 关键词： 分子动力学模拟 纳米碳团簇     

金属钛中氦团簇融合的分子动力学模拟

运用分子动力学方法研究了金属钛中氦的扩散聚集行为.在300—800K的温度范围内,模拟了钛基底中氦团簇之间的融合过程.研究发现,温度的升高会加快氦团簇的融合.在300—800K,融合后的氦团簇在所模拟的时间尺度内三维结构保持不变.模拟结果还表明,常温下氦团簇之间的吸引力是导致氦团簇融合的重要因素. 关键词： 氦团簇 团簇融合 分子动力学模拟     

低温下Au959团簇负载于MgO(100)表面后结构变化的分子动力学研究

采用基于Chen-Mbius反演方法,从金属/金属氧化物界面第一性原理计算的粘结能结果中推导出的Au/MgO原子间相互作用势的正则系综(NVT)分子动力学,模拟了在10 K条件下,Au959团簇负载于MgO(100)表面后团簇结构的变化.根据原子对分析技术和对分布函数的分析表明,由于团簇界面处原子间距与载体原子间距相匹配,置于载体上的Au团簇经过一个变形过程后,较其孤立自由表面时的团簇体积变大.     

铜原子纳米团簇热力学性质的分子动力学模拟研究

本文利用分子动力学模拟方法，研究了CuN（N=80、140、216、312、408、500、628和736）纳米团簇在热化和冷凝过程中结构和热力学性质的变化，模型采用的是Johnson的EAM作用势．模拟结果表明：铜团簇的熔点和凝固。点随其尺寸线性增加，并逐渐向大块晶体靠拢；所有团簇的凝固。羔都低于熔点，出现凝固过程中的滞后现象；在熔点和凝固点附近，团簇都具有负热容特性，负热容是由相变前后团簇内部结构突变引起的。     

纳米团簇熔化过程的分子动力学模拟

本文采用分子动力学结合嵌入原子多体势,模拟了不同半径的Ni纳米团簇的升温熔化过程,研究团簇尺寸对熔点和表面能的影响.模拟结果表明:团簇的熔点显著低于体材料的熔点.团簇熔化的过程首先是在团簇的表面出现预熔,然后向团簇内部扩展,直到整个团簇完全熔为液态.在模拟的纳米尺度范围内,团簇的熔点与团簇尺寸基本成线性关系.团簇的表面能随着团簇尺寸的增大而减小,而且表面能均高于体材料的表面能.     

低温下Cu13团簇负载于Cu(001)表面上结构变化的分子动力学研究

应用分子动力学方法研究温度为10和50 K时具有二十面体结构的Cu13团簇以不同接触条件与Cu(001)表面结合后的结构变化,原子间的相互作用势采用Johnson的嵌入原子方法模型.通过基于原子密度分布函数的分析表明,负载团簇与表面的结合能主要受团簇与载体相接触的最低层原子数及这些原子所具有的不同几何构型影响,同时更高层的原子呈现出不同的几何结构.温度为10 K时,负载团簇的初始位置对团簇几何结构和结合能影响较大.     

金属钛中氦团簇生长行为的分子动力学研究

运用分子动力学方法模拟了常温下金属钛中氦团簇的生长过程.从能量的角度考察了氦团簇的生长机理.研究发现,随着氦团簇的生长,在氦团簇周围逐渐形成位错环缺陷,氦团簇与氦原子的结合能有逐渐下降的趋势,当氦团簇生长到一定尺寸时会通过发射周围缺陷以使得结合能上升,从而增强了进一步吸收氦原子的能力.模拟还发现,随着氦团簇的继续生长,氦团簇的形状由原来的不规则结构逐渐变成了较为规则的棱柱形结构,在随后的生长过程中其生长仅在（001）平面进行,沿［001］轴的厚度几乎不变. 关键词： 氦团簇 缺陷发射 分子动力学模拟     

熔融Cu57团簇在急冷过程中弛豫和局域结构转变的分子动力学研究

采用基于嵌入原子方法的正则系综分子动力学研究熔融Cu₅₇团簇在急冷过程中的弛豫及其局域结构变化.通过对弛豫过程中均方位移、非相干中间散射函数和非Gauss参数三种函数和原子键对随急冷温度不同所发生变化的分析表明,在经过短时间的原子剧烈运动后,急冷温度极大地影响着团簇内原子结构弛豫过程.急冷温度较高时,原子在经历短时间剧烈运动的β弛豫后,进入α弛豫区后以扩散运动为主,随后原子运动表现为非扩散性的原子局域结构重排,团簇内没有出现明显的成核结构.随着温度的降低,原子局域结构的变化在经过短时间原子剧烈运动的β弛豫后,在α弛豫区原子运动表现为扩散性运动,并出现一定数量的不稳定二十面体结构.当急冷温度很低时,在进入α弛豫区后,团簇结构变化逐渐表现为非扩散性原子局域结构重排,形成相当数量的稳定成核二十面体结构.     

Hydration of the chloride ion in concentrated aqueous solutions using neutron scattering and molecular dynamics

Neutron scattering experiments were performed on 6 m LiCl solutions in order to obtain the solvation structure around the chloride ion. Molecular dynamics simulations on systems mirroring the concentrated electrolyte conditions of the experiment were carried out with a variety of chloride force-fields. In each case the simulations were run with both full ionic charges and employing the electronic continuum correction (implemented through charge scaling) to account effectively for electronic polarisation. The experimental data were then used to assess the successes and shortcomings of the investigated force-fields. We found that due to the very good signal-to-noise ratio in the experimental data, they provide a very narrow window for the position of the first hydration shell of the chloride ion. This allowed us to establish the importance of effectively accounting for electronic polarisation, as well as adjusting the ionic size, for obtaining a force-field which compares quantitatively to the experimental data. The present results emphasise the utility of performing neutron diffraction with isotopic substitution as a powerful tool in gaining insight and examining the validity of force-fields in concentrated electrolyte solutions.     

金属Ni纳米线凝固行为的分子动力学模拟

采用EAM镶嵌原子作用势，通过经典的分子动力学模拟方法研究了不同冷却速度下的金属Ni纳米线的凝固行为，并给出了纳米线在凝固区域的结构演变过程.利用键对分析技术研究了在不同冷却速度下体系中的原子团簇在降温过程中的变化情况.研究表明，纳米线的凝固起始于表面原子，并且随着冷却速度的降低，Ni纳米线的微观结构从非晶态过渡到多壳螺旋结构，最终达到稳定的面心立方结构.多壳螺旋结构同时具有确定的结晶温度和长程无序、短程有序的非晶结构的特征. 关键词： 纳米线 凝固行为 分子动力学 键对分析     

氯化钠晶体在高温高压下热物理参数的分子动力学计算

用分子动力学方法和优化的Tosi-Fumi有效两体势函数计算了NaCl晶体温度从298?K到773?K的等温压缩线和体积模量、热膨胀系数和格临爱森系数等热物理参量.计算结果表明，NaCl晶体在温度从298?K到1073?K、压力从0到80?GPa范围内的格临爱森系数γ=γ ₀(V/V₀)^q的近似公式成立，且q≈1078. 关键词： 分子动力学计算 格临爱森系数 氯化钠晶体     

Molecular dynamic studies of the solubility of sodium chloride: fast calculations using seed crystalline cluster probe

Theoretical predictions of solubility, typically accomplished by comparing the chemical potential of pure solid and solution, currently suffer from a lack of accuracy. We suggest an alternative method for predicting solubility based on molecular dynamics simulations of the behaviour of a small seed crystalline cluster probe in solutions of varying concentrations. The size dynamics of a properly chosen seed cluster that dissolves in unsaturated solutions and grows in size in supersaturated solutions is indicative of the saturation point. This approach is illustrated by its application to NaCl in water.     

Simulation of plasma doping process by using the localized molecular dynamics method

Plasma doping is the candidate for semiconductor doping. Accurate simulation of the doping technology is needed for the advanced integrated circuit manufacturing. In this paper, the plasma doping process simulation is performed by using the localized molecular dynamics method. Models that involve the statistics of the implanted compositions, angles and energies are developed. The effect of the model on simulation results is studied. The simulation results about the doping concentration profile are supported by experimental data.     

The structure,dynamics and solvation mechanisms of ions in water from long time molecular dynamics simulations: a case study of CaCl2 (aq) aqueous solutions

Systematic long time (5–20 ns) molecular dynamics (MD) simulations have been carried out to study the structural and dynamical properties of CaCl₂ aqueous solutions over a wide range of concentrations (≤9.26 m) in this study. Our simulations reveal totally different structural characteristics of those yielded from short time (≤1 ns) MD simulations [A.A. Chialvo and J.M. Simonson, J. Chem. Phys. 119, 8052 (2003); T. Megyes, I. Bako, S. Balint, T. Grosz, and T. Radnai, J. Mol. Liq. 129, 63 (2006)]. An apparent discontinuity was found at 4–5 m of CaCl₂ in various properties including ion–water coordination number and self-diffusion coefficient of ions, which were first noticed by Phutela and Pitzer in their thermodynamic modelling [R.C. Phutela and K.S. Pitzer, J. Sol. Chem. 12, 201 (1983)]. In this study, residence time was first taken into consideration in the study of Ca²⁺–Cl^? ion pairing, and it was found that contact ion pair and solvent-sharing ion pair start to form at the CaCl₂(aq) concentrations of about 4.5 and 4 m, respectively, which may be responsible for the apparent discontinuity. In addition, the residence time of water molecules around Ca²⁺ or Cl^? showed that the hydration structures of Ca²⁺ and Cl^? are flexible with short residence time (<1 ns). It needs to be pointed out that it takes much longer simulation time for the CaCl₂–H₂O system to reach equilibrium than what was assumed in previous studies.     

Studies on the structures and interactions of glutathione in aqueous solution by molecular dynamics simulations and NMR spectroscopy

All-atom molecular simulations and temperature-dependent NMR have been used to investigate the conformations and hydrogen bonds of glutathione (GSH) in aqueous solution. The simulations start from three different initial conformations. The properties are characterized by intramolecular distances, radius of gyration, root-mean-square deviation, and solvent-accessible surface. GSH is highly flexible in aqueous solutions in the simulations. Moreover, conformations can covert between “extended” and “folded” states. Interestingly, the two different hydrogen atoms in cysteine (H_N2) and glycin (H_N3) show different capabilities in forming NH?OW hydrogen bonds. The temperature-dependent NMR results of the different amide hydrogen atoms also show agreements with the MD simulations. Competing formation of GSH hydrogen-bonding interactions in aqueous solutions leads to hydrogen-bonding networks and the distribution of conformations. These changes will affect the activity of GSH under physiological conditions.     

Cu87团簇在升温与急冷过程中结构变化的原子尺度模拟

本文采用基于嵌入原子法的正则系综分子动力学方法在原子尺度上计算了包含87个原子的Cu87金属团簇在连续升温和急冷降温时的结构演化过程。根据原子平均势能、对分布函数、原子堆积结构和主要原子键对数目随温度的变化表明,温度的不同极大地影响团簇内的原子堆积结构。在升温过程中,随着温度的升高,团簇内原子堆积结构出现由密排六方、二十面体直到无序堆积的变化。在急冷降温过程中,随着急冷温度的降低,团簇内由出现的一定数量的二十面体和面心立方的局域结构、数量不一的HCP,FCC和二十面体局域结构,直到急冷温度较低时的一定数量的二十面体局域结构。     

Mechanism of microweld formation and breakage during Cu-Cu wire bonding investigated by molecular dynamics simulation

Currently,wire bonding is the most popular first-level interconnection technology used between the die and package terminals,but even with its long-term and excessive usage,the mechanism of wire bonding has not been completely evaluated.Therefore,fundamental research is still needed.In this study,the mechanism of microweld formation and breakage during Cu-Cu wire bonding was investigated by using molecular dynamics simulation.The contact model for the nanoindentation process between the wire and substrate was developed to simulate the contact process of the Cu wire and Cu substrate.Elastic contact and plastic instability were investigated through the loading and unloading processes.Moreover,the evolution of the indentation morphology and distributions of the atomic stress were also investigated.It was shown that the loading and unloading curves do not coincide,and the unloading curve exhibited hysteresis.For the substrate,in the loading process,the main force changed from attractive to repulsive.The maximum von Mises stress increased and shifted from the center toward the edge of the contact area.During the unloading process,the main force changed from repulsive to attractive.The Mises stress reduced first and then increased.Stress concentration occurs around dislocations in the middle area of the Cu wire.