A study of cavitation nucleation in pure water using molecular dynamics simulation

To discover the microscopic mechanism responsible for cavitation nucleation in pure water, nucleation processes in pure water are simulated using the molecular dynamics method. Cavitation nucleation is generated by uniformly stretching the system under isothermal conditions, and the formation and development of cavitation nuclei are simulated and discussed at the molecular level. The processes of energy, pressure, and density are analyzed, and the tensile strength of the pure water and the critical volume of the bubble nuclei are investigated. The results show that critical states exist in the process of cavitation nucleation. In the critical state, the energy, density, and pressure of the system change abruptly, and a stable cavitation nucleus is produced if the energy barrier is broken and the critical volume is exceeded. System pressure and water density are the key factors in the generation of cavitation nuclei. When the critical state is surpassed, the liquid is completely ruptured, and the volume of the cavitation nucleus rapidly increases to larger than 100 nm³; at this point, the surface tension of the bubble dominates the cavitation nucleus, instead of intermolecular forces. The negative critical pressure for bubble nucleation is -198.6 MPa, the corresponding critical volume is 13.84 nm³, and the nucleation rate is 2.42×10³² m^-3·^-1 in pure water at 300 K. Temperature has a significant effect on nucleation: as the temperature rises, nucleation thresholds decrease, and cavitation nucleation occurs earlier.     

石墨烯纳米孔道内受限水介电常数影响因素的分子动力学模拟

受限水的介电特性对于石墨烯电容器的储能效率有着重大的影响.因此本文使用分子动力模方法,研究了不同石墨烯纳米孔道宽度(0.812～10 nm)下内受限水介电常数的分布及其影响因素.结果表明,在石墨烯纳米孔道内部受限水双电层结构可分为空隙层、界面层和体相层三个部分.其整体的介电常数随孔道大小的降低而线性减少.这是由内部体相层的宽度的降低而导致的.此外电极电压大小和石墨烯-水相互作用参数ε_CO的强弱也会显著改变电极表面的双电层(Electric Double Layers, EDLs)结构.其中电压的增大使得介电常数分布的震荡的程度也随之增加,最终导致了整体介电常数的减小.与之类似,亲水态的石墨烯表面(高ε_CO)下受限水分布的震荡程度也显著增加,这导致了整体介电常数的降低.     

低密度Ar熔化及结晶的分子动力学模拟

采用有位移力的ＬｅｎｎａｒｄＪｏｎｅｓ（１２－６）势对微正则系综下低密度Ａｒ系统（简约数密度为ρ＝０．８５）的一级相变过程进行了细致的分子动力学模拟,发现Ａｒ系统的熔化过程是原子的崩塌过程,结晶过程和理想的完整晶体不同,是一活化过程：形核长大过程随温度的降低进行,原来均匀分布在系统中的自由体积呈集中分布,由此系统达到了更稳定的结晶状态。     

Cu固液界面能的分子动力学计算

本文使用分子动力学方法对液相Cu中不同半径晶胚的生长和熔化行为进行了研究. 随着半径的增加, 晶胚生长的临界温度升高. 临界形核过冷度和晶胚半径倒数成正比关系, 这和经典形核理论一致. 由上述关系计算得到Cu的Gibbs-Thomson系数为1.12×10^-7 K·m, Cu的固液界面能为0.146 J/m², Cu的Turnbull系数为0.416, 这些计算值均与实验值一致. 关键词： 分子动力学 经典形核理论 固液界面能     

A molecular dynamics simulation on the convective heat transfer in nanochannels

The understanding of the flow and heat transfer processes for fluid through micro- and nanochannels becomes imperative due to its wide application in micro- and nano-fluidic devices. In this paper, the method to simulate the convective heat transfer process in molecular dynamics is improved based on a previous study. With this method, we simulate a warm dense fluid flowing through a cold parallel-plate nanochannel with constant wall temperature. The characteristics of the velocity and temperature fields are analysed. The temperature difference between the bulk average temperature of fluid and the wall temperature decreases in an exponential form along the flow direction. The Nusselt number for the laminar flow in parallel-plate nanochannel is smaller than its corresponding value at macroscale. It could be attributed to the temperature jump at the fluid–wall interface, which decreases the temperature gradient near the wall. The results also reveal that the heat transfer coefficient is related to the surface wettabilities, which differs from that in the macroscopic condition.     

Influence of water on HFO-1234yf oxidation pyrolysis via ReaxFF molecular dynamics simulation

ABSTRACT

The effect of water molecules on HFO-1234yf oxidation pyrolysis was investigated by ReaxFF-molecular dynamics simulation from 1900 to 4200?K. The initial pyrolysis of HFO-1234yf starts around 2500?K and the water molecules participate in chemical reactions at 2800?K when the reactants pyrolysis reached the highest reaction rate. The primary products including HF, COF₂ and CO₂ are observed at 2600, 2700 and 2900?K, respectively. The influence of water molecules on products is mainly reflected in the promotion activity on the conversion from COF₂ to CO₂ and the generation of HF molecules. Four formation pathways are observed and calculated to further elucidate the procedure of pyrolysis. The main conversion process from H₂O to HF is the ?F?+?H₂O?=?HF+?OH reaction, and the paths from H₂O to ?OH radical and COF₂ to ?CFO radical which are promoted by ?F and ?H radical, respectively, have relatively low energy barriers of 10.44 and 40.29?kJ/mol, and both reaction processes released HF molecules.     

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

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

纳米沟槽表面黏着接触过程的分子动力学模拟研究

采用大规模分子动力学方法研究了刚性球型探头与具有不同纳米沟槽基体表面的黏着接触过程,探讨了表面沟槽结构对载荷-位移曲线、接触引力和拉离力以及材料转移的影响规律.研究结果表明:在相同的压入深度下,与原子级光滑表面的黏着接触过程相比,刚性探头与具有纳米沟槽结构基体表面的接触压力较小,接触加载过程中的引力作用范围较大,并伴随载荷的多次跳跃,且接触引力和拉离力均有减小;当沟槽深度相同时,随着沟槽宽度的增大,接触引力和拉离力逐渐减小,当沟槽宽度逐渐接近探头与光滑表面的接触直径时,接触引力和拉离力又逐渐增大,趋于接近探头与光滑表面的接触过程;当沟槽宽度相同时,随着沟槽深度的增大,接触引力相对减小,拉离力变化不大.     

非绝热分子动力学的量子路径模拟

基于近期发展的经典-量子混合模拟非绝热分子动力学的量子路径方案,本文对5个典型势能面模型进行了模拟,包括单交叉模型、双交叉模型、拓展耦合模型、哑铃模型以及双弓模型.由于难以在严格意义上得到退相干速率,数值模拟中,我们比较了三个不同的退相干速率公式,包括冻结高斯波包近似退相干速率、能量分辨速率以及力分辨速率.在模拟过程中,我们恰当地处理了势能面跳跃时的能量守恒和力的反向问题.通过与全量子动力学模拟的精确结果进行对比发现,对于结构较简单的势能面模型,三种退相干速率都能得到较好的结果;然而对于较复杂的势能面模型,由于复杂量子干涉的原因,与其他混合经典-量子动力学方案类似,量子路径方案仍然难以得到较准确的结果.如何发展更加有效的混合经典-量子模拟方案,是未来研究的重要课题.     

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)].     

Cu-Au体系非对称异质外延行为的分子动力学研究

利用分子动力学模拟方法研究了Cu/Au(001)和Au/Cu(001)异质外延岛的演化行为. 研究结果显示：Cu-Au体系的相互外延行为呈现出明显的非对称性. Cu在Au(001)基体表面可以形成完整的外延结构，而Au在Cu(001)基体表面外延将导致失配位错的出现. 导致非对称外延生长行为的根本原因是外延岛的应变状态的差异和外延岛自身性质的不同. 随着外延岛的长大，Cu外延岛与Au(001)基体的微观失配度由最初的接近宏观失配度的9%左右迅速单调下降，并最终趋于晶格匹配；而Au在Cu(001)基体表面外延的微观失配度则呈现出振荡增加趋势. Cu/Au(001)体系的基体形变主要发生在外延岛的边缘，而Au/Cu(001)体系的基体形变主要发生在外延岛内部所对应的区域.     

Study of mechanical properties of amorphous copper with molecular dynamics simulation

The formation and mechanical properties of amorphous copper are studied using molecular dynamics simulation. The simulations of tension and shearing show that more pronounced plasticity is found under shearing, compared to tension. Apparent strain hardening and strain rate effect are observed. Interestingly, the variations of number density of atoms during deformation indicate free volume creation, especially under higher strain rate. In particular, it is found that shear induced dilatation does appear in the amorphous metal.     

Study of mechanical properties of amorphous copper with molecular dynamics simulation

The formation and mechanical properties of amorphous copper are studied using molecular dynamics simulation. The simulations of tension and shearing show that more pronounced plasticity is found under shearing, compared to tension. Apparent strain hardening and strain rate effect are observed. Interestingly, the variations of number density of atoms during deformation indicate free volume creation, especially under higher strain rate. In particular, it is found that shear induced dilatation does appear in the amorphous metal.     

模拟盒参数对汽液体系模拟过程的影响

为探讨模拟盒参数对分子动力学模拟影响这个长期被忽略的问题,采用Verlet-List搜寻法和Leapfrog差分算法,对分子数目为1000的汽液共存体系进行分子动力学模拟,得到了系统温度为100K时,汽液两相密度、表面张力随系统切片数、汽相空间尺寸、液膜厚度以及步长等模拟盒参数的变化情况.研究结果表明,系统切片数不会影响表面张力及两相密度;体系z轴无量纲长度为60时,模拟值与实验值吻合较好;液膜无量纲厚度达到13以后,液相密度将趋于稳定,两侧界面分子不会相互影响;分子动力学步长不影响平衡后的系统密度,但对表面张力的计算影响较明显.     

A molecular dynamics simulation of nanoscale convective heat transfer with the effect of axial heat conduction

Effective heat dissipation from nano-fluidic devices is sometimes necessary to ensure their performance and lifespan. In the molecular dynamics simulation of nanoscale convective heat transfer, thermostats cannot be directly applied to the fluid because of the non-uniform temperature distribution. Periodic boundary is typically utilised, but unrealistic axial heat conduction exists when there is a temperature difference between the outlet and images of inlet atoms. In this paper, the effect of axial conduction caused by periodic boundary is investigated through the Péclet number (Pe). Taking viscous dissipation into consideration, the magnitude of outlet thermal diffusion is observed to decrease with increasing Pe. The local average temperature of fluid changes in an exponential form except in the region close to the outlet. Results show that the contribution of outlet axial conduction to the local average temperature is less than 2.0% when Pe > 10. The main reason is that the magnitude of fluid velocity and viscous heat dissipation in nanochannels is much larger than that in macro-channels at the same Péclet number.     

正则系综条件下空化空泡形成的分子动力学模拟

液体中空化现象的研究对减少空蚀破坏, 提高空化空泡的有效利用具有重要意义.本文采用分子动力学模拟的方法, 对正则系综条件下系统中空化的发生特性进行研究, 分析空化发生的机理, 讨论温度、数密度等因素对空化发生的影响, 并与格子Boltzmann方法进行了比较.模拟结果表明: 温度和数密度的变化, 都对系统中的空化产生显著影响.其中, 温度升高, 使系统中空化空泡的形成由稳定变得不稳定, 最终难以形成.数密度降低, 则会促进空化空泡的形成.随着数密度的减小, 温度对空化空泡形成的影响程度下降.     

碳纳米锥力学特性的分子动力学研究

结合原子间短程作用势(Brenner势)和长程作用势(Lennard-Jones势),利用分子动力学方法对各种锥角的碳纳米锥进行拉伸和压缩实验,获得其载荷-应变关系曲线、受拉/压载荷极限、应变极限和构形演变等力学特性,并与等量原子组成的碳纳米管进行比较研究.研究结果表明,等量碳原子组成的碳纳米锥的受拉/压载荷极限随着锥角的增大先是增大后减小,受拉/压应变极限则随着锥角的增大而增大.与碳纳米锥相比,等量碳原子组成的碳纳米管的受拉/压载荷极限和应变极限显得既不突出也不逊色.在受压构形演化方面,与碳纳米管丰富的径向屈曲/扭转/侧向屈曲组合形变不同,112.88°和83.62°锥角的碳纳米锥受压沿轴向完美内陷,而60.0°和38.94°锥角的碳纳米锥受压发生侧向屈曲.     

分子动力学模拟方法在非线性光学中的应用

有机分子非线性光学材料在频率转换、电光调制和双光子吸收等方面具有重要的应用.介绍了近年来分子动力学模拟方法在有机分子非线性光学性质理论研究中的主要应用,包括电场极化效应、局域场因子、非线性极化率和双光子吸收等.此外,结合最新的科研工作,介绍了分子动力学模拟方法在溶剂效应和聚集效应研究中发挥的重要作用. 关键词： 分子动力学模拟 溶剂效应 聚集效应 双光子吸收     

Temperature dependence of the thermal conductivity of water: a molecular dynamics simulation study using the SPC/E model

In this study, molecular dynamics simulations of SPC/E (extended simple point charge) water model have been carried out in the canonical (NVT fixed) ensemble over the range of temperatures 300–550 K with Ewald summation. The evaluated thermal conductivity for SPC/E water overestimates the experimental data at 300–550 K. In accordance with experimental data, SPC/E predicts a maximum in the thermal conductivity at 400 K. The temperature dependence of thermal conductivity of SPC/E water was discussed.