受限于纳米狭缝中的四缔合Lennard-Lones流体的相平衡和界面张力研究

分别用改进的基础测量理论和平均球近似理论表达短程作用和长程作用对四缔合Lennard-Jones流体的过剩自由能的贡献. 在密度函泛理论的框架下, 研究了平均密度等温线, 密度分布, 未缔合分子在平衡汽相和液相中的分布, 相平衡以及平衡时的界面张力等热力学性质. 分析了缔合能量, 流体-固体作用和孔宽对受限于纳米狭缝中的四缔合Lennard-Jones流体相行为的影响.     

Asymptotic analysis of microscopic impenetrability constraints for atomistic systems

We analyze systems of atomistic interactions on a triangular lattice allowing for fracture under a geometric condition on the triangles corresponding to a microscopic impenetrability constraint. Such systems can be thought as a computational simulation of materials undergoing brittle fracture. We show that in the small-deformation regime such approximation can be validated analytically in the framework of variational models of fracture. Conversely, in a finite-deformation regime various pathologies show that the continuum approximation of such a system differs from the usual variational representations of fracture and either needs new types of formulations on the continuum, or a proper interpretation of the atomistic constraints limiting their range and adapting them to a dynamical framework.     

Impact of different potentials on the structures and energies of clusters

Cluster studies have attracted much interest in the past decades because of their extraordinary properties. To describe the interaction between atoms or molecules and predict the energies and structures, potential functions are developed. However, different potentials generally produce different structures and energies for a cluster. To study the effect of potentials on the structure of a cluster, He clusters in the size range of 13-140 are investigated by Lennard-Jones (LJ), Pirani, and Hartree-Fock-dispersion individual damping (HFD-ID) potential with dynamic lattice searching (DLS) method. Potential function curves, cluster structures, bonds, and energies of the global minima are compared. The results show that cluster energies decrease with the values of the potential functions, the differences between structures depend upon the disagreements of the potentials, and the preferable motif of a cluster changes from icosahedron to decahedron with the increase of the derivative of the short-range part of the potentials.     

Crossover from normal diffusion to single-file diffusion of particles in a one-dimensional channel: LJ particles in zeolite zsm-22

The crossover from normal Fickian diffusion, where mean-squared displacement goes as t to single-file diffusion, where <z²> t?^0.5 is studied as function of particle size confined in zeolite zsm-22 using molecular dynamics simulations. The simulation results indicate that the crossover is smooth as the particle size increases and the diffusion reaches single file through a series of subdiffusion processes. A small number of mutual passage of particles can destroy the correlated single-file diffusion. The density dependence on single-file mobility obtained from molecular dynamics simulations are in very good agreement with theoretical predictions.     

A dynamic lattice searching method with interior operation for unbiased optimization of large Lennard-Jones clusters

For improving the efficiency of dynamic lattice searching (DLS) method for unbiased optimization of large Lennard-Jones (LJ) clusters, a variant of the interior operation (IO) proposed by Takeuchi was combined with DLS. The method is named as DLS-IO. In the method, the IO moves outer atoms with higher energy toward the coordinates center, i.e., (0, 0, 0), of a cluster and a local minimization (LM) follows each IO. This makes the interior atoms more compact and the outer atoms more uniformly distributed with lower potential energy. Therefore, the starting structure for DLS operations is closer to the global optimum compared with the randomly generated structures. On the other hand, a method to identify the central atom is proposed for the central vacancy problem. Optimizations of LJ(500), LJ(561), LJ(660), LJ(665), and LJ(670) were investigated with the DLS-IO, and the structural transition during the optimization was analyzed. It was found that the method is efficient and unbiased for optimization of large LJ clusters, and it may be a promising approach to be universally used for structural optimizations.     

Multibaric-multithermal ensemble molecular dynamics simulations

We present new generalized-ensemble molecular dynamics simulation algorithms, which we refer to as the multibaric-multithermal molecular dynamics. We describe three algorithms based on (1) the Nosé thermostat and the Andersen barostat, (2) the Nosé-Poincaré thermostat and the Andersen barostat, and (3) the Gaussian thermostat and the Andersen barostat. The multibaric-multithermal simulations perform random walks widely both in the potential-energy space and in the volume space. Therefore, one can calculate isobaric-isothermal ensemble averages in wide ranges of temperature and pressure from only one simulation run. We test the effectiveness of the multibaric-multithermal algorithm by applying it to a Lennard-Jones 12-6 potential system.     

钯团簇形成和增长机理的Monte Carlo研究

利用Monte Carlo（MC）方法和Lennard-Jones加Axilord-Teller （LJ＋AT）势能函数,研究了气相中钯团簇的形成过程和增长机理.发现具有二十面体结构的Pd13团簇可以在气相中自发形成,较大的团簇通过在Pd13二十面体结构的表面添加原子组成四面体的方式形成.通过分析团簇结构和能量之间的关系,发现除了Pd13和Pd55以外,Pd19和Pd39团簇也具有五次对称性,都是比较稳定的结构.     

Vapour–liquid phase equilibrium and surface tension of fully flexible Lennard–Jones chains

We report measurements of the vapour–liquid coexistence densities and surface tension of fully flexible Lennard–Jones chain molecules ranging in length from 4 to 60 beads. We demonstrate that the surface tension for all chain lengths collapses to a single master curve when plotted according to the universal parachor correlation. We find a universal parachor exponent 3.79 ± 0.05 for conditions close to the critical point, with a deviation observed for the longest chains far below the critical point.