Relation between Tolman length and isothermal compressibility for simple liquids

The Tolman length δ 0 of a liquid with a plane surface has attracted increasing theoretical attention in recent years,but the expression of Tolman length in terms of observable quantities is still not very clear.In 2001,Bartell gave a simple expression of Tolman length δ 0 in terms of isothermal compressibility.However,this expression predicts that Tolman length is always negative,which is contrary to the results of molecular dynamics simulations(MDS) for simple liquids.In this paper,this contradiction is analyzed and the reason for the discrepancy in the sign is found.In addition,we introduce a new expression of Tolman length in terms of isothermal compressibility for simple fluids not near the critical points under some weak restrictions.The Tolman length of simple liquids calculated by using this formula is consistent with that obtained using MDS regarding the sign.     

Theoretical study of the permeation of water through TiO2 nanotubes using molecular dynamics simulation

In this study, the theoretical structures of armchair (6, 6) and zigzag (12, 0) TiO₂ nanotubes (TiNTs) were constructed by rolling the (101) layer of an anatase TiO₂ crystal. The (101) layer was made using Materials Studio (MS) by cutting the cleave plane (101) of the anatase TiO₂ crystal. Based on these structures, the basic properties of TiO₂ nanotubes were investigated using MS. Molecular dynamics simulations were performed using the software NAMD to investigate the status and permeation of water through the TiO₂ nanotubes. Structure analysis shows that both the inner and outer walls of the structures were terminated with oxygen atoms. The thicknesses of single tube walls are smaller than that of a perfect triple layer (2.20?Å) in bulk anatase TiO₂. With regard to the bulk Ti–O bond length, the Ti–O bonds in the outer layer are elongated, and are shortened in the inner layer. Molecular dynamics simulation shows that the water molecules in the nanotubes move back and forth, as in one-dimensional Brownian motion. Moreover, the penetration properties of TiNTs are associated with their radii, with the TiNT with larger radii having better penetration properties. Thus, when used in drug delivery or filtration systems, armchair TiNT has a better effect than zigzag TiNT.     

Conformational properties of cyclooctane: a molecular dynamics simulation study

Atomistic molecular dynamics simulations have been used to elucidate the conformational properties of cyclooctane in the gas and bulk liquid phases. Accurate reproduction of the gas phase structure, and of the liquid phase densities and solubility parameters have been used as prerequisites to the prediction of conformational properties. The gas phase results clearly indicate the presence of a conformational mixture consisting of the crown, boat-chair, twist-boat-chair and boat-boat conformers at all temperatures (161, 313 and 400 K) studied. The fraction of the crown family of conformers was found to be relatively insensitive to temperature. However, the relative concentrations of the twist-boat-chair and boat-chair conformations was found to be highly temperature dependent with the boat-chair being favoured at low temperatures. Bulk packing was found to have a profound effect on the conformational properties in the liquid phase. At the temperatures studied (313 and 400 K) the boat-chair family was predominant, with the crown and boat families being essentially absent. The twist-boat-chair conformation was detected in the liquid phase at both temperatures. The pseudorotation pathway for the twist-boat-chair to boat-chair interconversion was prevalent in both gas and liquid phases establishing the conformational flexibility and the relative importance of the twist-boat-chair conformer in comparison to the crown family. The study successfully explains the separate experimental findings in both the gas and liquid phases of cyclooctane.     

Relaxation dynamics and power spectra of liquid water: a molecular dynamics simulation study

We present molecular dynamics simulations of liquid water at normal and supercooled conditions. Autocorrelation functions (ACFs) of several structural quantities and their fourier transforms are obtained and analysed. Structural correlations and relaxation times increase linearly with degree of supercooling. Power spectra of ACFs show increase in librational motion of liquid water with cooling. These modes intensify with supercooling because of structuring and ordering of water molecules. Overall, liquid water structure is homogenous over the temperatures and pressures studied and undergoes fluctuation–dissipation in its local-density variations [English and Tse, Phys. Rev. Lett. 106, 037801 (2011)].     

A molecular dynamics study of fuel droplet evaporation in sub- and supercritical conditions

Evaporation processes of a fuel droplet under sub- and supercritical ambient conditions have been studied using molecular dynamics (MD) simulations. Suspended n-dodecane droplets of various initial diameters evaporating into a nitrogen environment are considered. Both ambient pressure and temperature are varied from sub- to supercritical values, crossing the critical condition of the chosen fuel. Temporal variation in the droplet diameter is obtained and the droplet lifetime is recorded. The time at which supercritical transition happens is determined by calculating the temperature and concentration distributions of the system and comparing with the critical mixing point of the n-dodecane/nitrogen binary system. The dependence of evaporation characteristics on ambient conditions and droplet size is quantified. It is found that the droplet lifetime decreases with increasing ambient pressure and/or temperature. Supercritical transition time decreases with increasing ambient pressure and temperature as well. The droplet heat-up time as well as subcritical to supercritical transition time increases linearly with the initial droplet size d₀, while the droplet lifetime increases linearly with ${\mathrm{d}}_0^2$. A regime diagram is obtained, which indicates the subcritical and supercritical regions as a function of ambient temperature and pressure as well as the initial droplet size.     

有结构壁面上液滴运动特征的耗散粒子动力学模拟

本文采用耗散粒子动力学方法, 研究了恒定外力驱动下液滴在有结构壁面上的运动过程. 通过研究液滴在通过壁面结构时前缘接触点和前进角的变化, 分析了液滴的运动特征. 研究结果表明, 在不同润湿性或不同外力驱动下, 存在使液滴运动最快的“最优”壁面结构, 并对其机理进行了讨论. 本文还探讨了壁面润湿性、热涨落以及外力对液滴运动状态的影响. 关键词： 微流体 有结构壁面 润湿性 耗散粒子动力学(DPD)     

Femtosecond laser ablation of metals:a molecular dynamics simulation study

Using molecular dynamics (MD) methods combining with two-step radiation heating model, the mechanisms of ablation and the thermodynamic states at Ni surface under femtosecond laser irradiation are investigated. Simulation results show that the main mechanisms of ablation are evaporation and tensile stresses generated inside the target. The velocity of stress wave is predicted to be nearly equal to sound velocity. The rates of ablation at different fluences obtained from simulations are in good agreement with experimental data. Superheating phenomenon is also discovered.     

DTO体系分子反应动力学

基于DTO分子( 1A1)的氢同位素效应,得到修正的Born-Oppenheimer(B-O)理论下多体展式分析势能函数.用准经典的Monte-Carlo轨迹法研究了D+ TO(0 ,0)和T+ DO(0 ,0)的分子反应动力学过程.结果表明:在碰撞能量较低时( < 209.2 kJ. mol - 1 ),D+ TO(0 ,0)和T+DO(0,0)反应主要生成O+ DT,并且该反应是无阈能的;有少量的交换反应产物DO (TO)生成,并伴有极少量的络合物产生.碰撞能大于209.2kJ . mol - 1后,逐渐出现分子被完全碰散成D,T ,O 原子的情形.反应 D+TO(0 ,0) →OD + T和 T +DO(0 ,0) →OT + D是无热反应但是有阈能存在.由于D和T原子的同位素效应,置换产物轨线存在非一致性.     

Dynamic characteristics of nanoindentation in Ni:A molecular dynamics simulation study

In this work,three-dimensional molecular dynamics simulation is carried out to elucidate the nanoindentation behaviour of single crystal Ni.The substrate indenter system is modelled using hybrid interatomic potentials including the manybody potential(embedded atom method) and two-body Morse potential.The spherical indenter is chosen,and the simulation is performed for different loading rates from 10 m/s to 200 m/s.Results show that the maximum indentation load and hardness of the system increase with the increase of velocity.The effect of indenter size on the nanoindentation response is also analysed.It is found that the maximum indentation load is higher for the large indenter whereas the hardness is higher for the smaller indenter.Dynamic nanoindentation is carried out to investigate the behaviour of Ni substrate to multiple loading-unloading cycles.It is observed from the results that the increase in the number of loading unloading cycles reduces the maximum load and hardness of the Ni substrate.This is attributed to the decrease in recovery force due to defects and dislocations produced after each indentation cycle.     

Solvation dynamics in methanol: Experimental and molecular dynamics simulation studies

We have investigated the ultrafast dynamics of methanol by time dependent fluorescent shift experiments and molecular dynamic simulations. The experiments were performed with two different probe molecules, 1-aminonaphthalene and coumarin 153. The molecular dynamic simulations employed these probes as well as small atomic and diatomic solutes. We find a previously unobserved fast decay component in the solvation response of methanol. The molecular dynamics results are in good agreement with this experimental result. The origin of this fast response and the linearity of the solvent response are discussed.     

Pt-Ag合金纳米线热力学性质的分子动力学模拟研究

基于原子嵌入势(EAM),采用分子动力学方法,对临界尺寸下的Pt_(0.95)Ag_(0.05)合金纳米线多边形结构的熔化行为进行了计算模拟.结果表明:径向尺寸对Pt_(0.95)Ag_(0.05)合金纳米线的熔点影响较为显著,而长度对其影响较小;引入林德曼因子得到的熔点和用势能-温度变化曲线找到的熔点基本一致;合金纳米线的染色原子由外向内运动;综合分析发现Pt_(0.95)Ag_(0.05)合金纳米线以先外后内的模式进行熔化.     

Pt-Ag合金纳米线热力学性质的分子动力学模拟研究

基于EAM原子嵌入势,采用分子动力学方法,对临界尺寸下的Pt0.95Ag0.05合金纳米线多边形结构的熔化行为进行了计算模拟.结果表明:径向尺寸对Pt0.95Ag0.05合金纳米线的熔点影响较为显著,而长度对其影响较小;引入林德曼因子得到的熔点和用势能-温度变化曲线找到的熔点基本一致;合金纳米线的染色原子由外向内运动;综合分析发现Pt0.95Ag0.05合金纳米线以先外后内的模式进行熔化.     

GPU-accelerated molecular dynamics simulation for study of liquid crystalline flows

We have developed a GPU-based molecular dynamics simulation for the study of flows of fluids with anisotropic molecules such as liquid crystals. An application of the simulation to the study of macroscopic flow (backflow) generation by molecular reorientation in a nematic liquid crystal under the application of an electric field is presented. The computations of intermolecular force and torque are parallelized on the GPU using the cell-list method, and an efficient algorithm to update the cell lists was proposed. Some important issues in the implementation of computations that involve a large number of arithmetic operations and data on the GPU that has limited high-speed memory resources are addressed extensively. Despite the relatively low GPU occupancy in the calculation of intermolecular force and torque, the computation on a recent GPU is about 50 times faster than that on a single core of a recent CPU, thus simulations involving a large number of molecules using a personal computer are possible. The GPU-based simulation should allow an extensive investigation of the molecular-level mechanisms underlying various macroscopic flow phenomena in fluids with anisotropic molecules.     

Structural stability of nano-sized crystals of HMX: A molecular dynamics simulation study

The interaction potential energy and heat of sublimation of nanoparticles of HMX crystal polymorphs are studied by using molecular dynamics methods with a previously developed force field [Bedrov, et al., J. Comput.-Aided Mol. Des. 8 (2001) 77]. Molecular dynamics simulations of nanoparticles with 10, 20, 30, 40, 50, 60, 70, 80, 90, and 100 molecules of HMX are carried out at 300 K. The intermolecular, intramolecular and total interaction energies per mole for the nanoparticles are calculated at 300 K. Then, we have calculated sublimation enthalpy of HMX crystal polymorphs with different sizes. For the all sizes, the β-HMX is found to be the most stable phase, due to having the least total interaction energy. Also, α-HMX is more stable than δ-HMX. An increase in the sublimation enthalpy with the size of the nanoparticle can be seen.     

SiH2体系的分子反应动力学

基于多体展式方法所导出的SiH₂（X¹A₁）分析势能函数,用准经典的Monte Carlo轨线法研究了Si（¹D_g）+H₂（0,0）和H（²S_g）+SiH（0,0）的原子与分子反应动力学过程.研究结果表明：Si（¹D_g）与H₂（0,0）的碰撞在低能时（小于209.20 kJ/mol）生成稳定的络合物SiH₂（X¹A₁）,该反应是无阈能反应;而H（²S_g）与SiH （0,0）碰撞不能生成稳定的络合物,主要发生交换反应H（²S_g）+SiH （0,0）→Si（¹D_g）+H₂（0,0）,该反应也是无阈能反应. 关键词： 2')" href="#">SiH₂ 势能函数 反应截面 轨线     

Solid-liquid phase transition in small water clusters: a molecular dynamics simulation study

Water clusters, (H₂O) _n , of varying sizes (n = 8, 12, 16, 20, 24, 28, 32, 36, and 40) have been studied at different temperatures from 0 to 200 K using molecular dynamics simulations. Transitions between solid and liquid phases were investigated to estimate the melting temperature of the clusters. Although the melting temperatures showed non-monotonic behaviour as a function of cluster size, their general tendency follows the classical relationship T _m ∝ n ^?1/3 to the cluster size n. Moreover, it was observed that the liquid-solid surface tension decreased with the cluster size in a similar way to the liquid-vapour surface tension in bulk water. Upon cooling, ice-like crystals were formed from the smaller clusters with n up to 20, while the larger clusters were transformed to glassy structures. The decrease in the glass transition temperature with the cluster size was observed to be much less than the corresponding melting temperature. The mutual order of the melting and glass-transition temperatures were found to be reversed compared with that observed for bulk water.     

Methane in carbon nanotube: molecular dynamics simulation

The behaviour of methane molecules inside carbon nanotubes at room temperature is studied using classical molecular dynamics simulations. A methane molecule is represented either by a shapeless super-atom or by a rigid set of five interaction centres localized on atoms. Different loadings of methane molecules ranging from the dense gas density to the liquid density, and the influence of flexibility of the CNT on structural and dynamic properties of confined molecules are considered. The simulation results show the decreases of the diffusion coefficient of methane molecules with density. At higher densities diffusion coefficient values are almost independent of molecular shape, but at low densities one observes faster motion of the super-atom molecule than that for the tetrahedral model of the molecule. For loadings of methane considered here the nanotube flexibility, introduced by the reactive empirical bond order (REBO) potential for interactions between carbon atoms of nanotube, does not have an effect on diffusivity of methane molecules, and its impact on the molecular structure is weak. It is found that methane molecules in the vicinity of the nanotube wall show tripod orientation with respect to the nanotube surface.     

Optical Kerr effect in nematic liquid crystals: a molecular dynamics simulation study

Molecular dynamics simulations of the Gay-Berne model have been undertaken and optical response properties calculated for the isotropic and nematic phases. The components of the optical response were calculated for the coordinate system with respect to the director. Furthermore, the response for the components was separated into the orientational and collision induced contributions in order to analyse the mechanism of the optical response. It was observed, in particular, that one of the depolarized components of the response function does not vanish after long times for the nematic phase, unlike in the isotropic phase. This means that the orientation of the director can be permanently changed by instant irradiation with polarized light. The results give a microscopic picture of the optical Kerr effect in nematic liquid crystals.     

U+CO体系的分子反应动力学研究

基于CUO分子(X3A″)的多体展式分析势能函数,用准经典的MonteCarlo轨迹法研究了U+CO(0,0)的分子反应动力学过程.结果表明:在碰撞能低时(<215kJ.mol-1),可以生成长寿命络合物CUO(X3A″),并且该络合反应是无阈能反应,这一结论与前文用多体项展式理论计算的CUO分子势能曲线结果一致.碰撞能大于4184kJ.mol-1后,逐渐出现置换产物UO和UC,随碰撞能进一步增大,分子将被完全碰散成U,C,O原子,而且反应U+CO(0,0)→UO+C,U+CO(0,0)→UC+O和U+CO(0,0)→U+O+C是有阈能反应基于CUO分子(X3A″)的多体展式分析势能函数,用准经典的MonteCarlo轨迹法研究了U+CO(0,0)的分子反应动力学过程.结果表明:在碰撞能低时(<215kJ.mol-1),可以生成长寿命络合物CUO(X3A″),并且该络合反应是无阈能反应,这一结论与前文用多体项展式理论计算的CUO分子势能曲线结果一致.碰撞能大于4184kJ.mol-1后,逐渐出现置换产物UO和UC,随碰撞能进一步增大,分子将被完全碰散成U,C,O原子,而且反应U+CO(0,0)→UO+C,U+CO(0,0)→UC+O和U+CO(0,0)→U+O+C是有阈能反应 关键词： 玻色-爱因斯坦凝聚体 V型三能级原子 双模压缩态光场 光场正交压缩     

多环芳烃在超临界环己烷中的溶解：分子动力学模拟研究

重油中大量存在的多环芳烃会经过缩合反应会形成稠环芳烃,而溶解扩散的限制会导致稠环芳烃之间的聚合从形成残焦,使得对重油的利用率大大降低,所以多环芳烃热解过程的溶解扩散是重油改质的关键因素之一.因此,基于环己烷对烃类良好的溶解性,本文采用分子动力学模拟的方法研究了多环芳烃及其混合物在超临界环己烷中的溶解行为,结果表明在不同温度与密度下超临界环己烷对于具有NAP均具有良好的溶解性,而温度对于Bghip溶解的影响较小.在多环芳烃混合物的油滴溶解过程中,NAP优先溶解到环己烷相中,而Bghip则集中在粒径减小的油滴中.通过计算体系中多环芳烃的径向分配函数与溶剂化自由能发现环己烷与多环芳烃间之间能相互吸引,环己烷在多环芳烃分子周围形成的溶剂壳能抑制多环芳烃分子间的聚合.温度升高以及密度的将降低会导致多环芳烃的内聚能密度降低,增加了多环芳烃与超临界环己烷之间的相互作用.环己烷密度越小以及温度的升高可以促进多环芳烃或多环芳烃混合物在超临界环己烷中的溶解.