Feasibility study of symmetric solution of molecular argon flow inside microscale nozzles

The computational cost of numerical methods in microscopic-scales such as molecular dynamics (MD) is a deterrent factor that limits simulations with a large number of particles. Hence, it is desirable to decrease the computational cost and run time of simulations, especially for problems with a symmetrical domain. However, in microscopic-scales, implementation of symmetric boundary conditions is not straight-forward. Previously, the present authors have successfully used a symmetry boundary condition to solve molecular flows in constant-area channels. The results obtained with this approach agree well with the benchmark cases. Therefore, it has provided us with a sound ground to further explore feasibility of applying symmetric solutions of micro-fluid flows in other geometries such as variable-area ducts. Molecular flows are solved for the whole domain with and without the symmetric boundary condition. Good agreement has been reached between the results of the symmetric solution and the whole domain solution. To investigate robustness of the proposed method, simulations are conducted for different values of affecting parameters including an external force, a flow density, and a domain length. The results indicate that the symmetric solution is also applicable to variable-area ducts such as micro-nozzles.     

An ab initio-based approach to phase diagram calculations for GaN(0 0 0 1) surfaces

Surface phase diagrams of GaN(0 0 0 1)-(2 × 2) and pseudo-(1 × 1) surfaces are systematically investigated by using our ab initio-based approach. The phase diagrams are obtained as functions of temperature T and Ga beam equivalent pressure p_Ga by comparing chemical potentials of Ga atom in the vapor phase with that on the surface. The calculated results imply that the (2 × 2) surface is stable in the temperature range of 700-1000 K at 10⁻⁸ Torr and 900-1400 K at 10⁻² Torr. This is consistent with experimental stable temperature range for the (2 × 2). On the other hand, the pseudo-(1 × 1) phase is stable in the temperature range less than 700 K at 10⁻⁸ Torr and less than 1000 K at 10⁻² Torr. Furthermore, the stable region of the pseudo-(1 × 1) phase almost coincides with that of the (2 × 2) with excess Ga adatom. This suggests that Ga adsorption or desorption during GaN MBE growth can easily change the pseudo-(1 × 1) to the (2 × 2) with Ga adatom and vice versa.     

不可恢复道路堵塞路径选择问题及其算法

针对不可恢复道路堵塞的路径选择问题，分析了堵塞发生的不同起始时间对通过被堵塞边的费用(时间)产生的影响，这种影响引起的后效性使得这个问题不适合用动态规划中逆序推算方法求解。本设计了一种算法(MDA)，对经典的Dijkstra算法中的权值计算进行了修正，证明了算法的复杂性为O(n^2)，并用该算法对模型进行了求解。最后通过一个算例分析，对算法进行了验证。     

复相介质及界面结构温度场跨尺度计算

本文针对复相介质及界面结构的跨层次热问题,提出了计算界面结构导热系数的分子动力学方法,复相介质宏观有效导热系数和温度函数的均匀化与多尺度渐近展开方法以及界面结构局域温度场的量子修正。建立了微观与宏观相关联的跨尺度计算模式,初步实现了从量子力学到连续模型的跨越。     

The Dynamical Structure Function of the One‐Component Plasma Revisited

We present high precision Molecular Dynamics simulations for the dynamical structure function, S (k, ω), of the classical Coulomb One Component Plasma (OCP), for a wide range of the coupling parameter Γ (from 0.05 to 10 000). We follow the positive‐to‐negative transition of the slope of the ω (k) dispersion curve at small wave numbers caused by the onset of correlations with increasing coupling. The high signal‐to‐noise ratio of the data over several orders of magnitude allows examination of a wide dynamical range of S (k, ω), including extreme values of ω and k, and the identification of waves at higher harmonics of the plasma frequency. Sum rules are found to be accurately satisfied. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Influences of the Structure of Lipids on Thermal Stability of Lipid Membranes

The binding free energy (BFE) of lipid to lipid bilayer is a critical factor to determine the thermal or mechanical stability of the bilayer. Although the molecular structure of lipids has significant impacts on BFE of the lipid, there lacks a systematic study on this issue. In this paper we use coarse-grained molecular dynamics simulation to investigate this problem for several typical phospholipids. We find that both the tail length and tail unsaturation can significantly affect the BFE of lipids but in opposite way, namely, BFE decreases linearly with increasing length, but increases linearly with addition of unsaturated bonds. Inspired by the specific structure of cholesterol which is a crucial component of biomembrane, we also find that introduction of carbo-ring-like structures to the lipid tail or to the bilayer may greatly enhance the stability of the bilayer. Our simulation also shows that temperature can influence the bilayer stability and this effect can be significant when the bilayer undergoes phase transition. These results may be helpful to the design of liposome or other self-assembled lipid systems.     

Laser ablation of metals by femtosecond pulses: Theoretical and experimental study

We have investigated ultrashort laser micromachining of metals, both from the point of view of the basic physical processes, and the technological implications. The process of hole drilling of Ni with ≈300 fs SHG (λ = 527 nm) Nd-glass and Al samples with 100 fs Ti:sapphire (λ = 800 nm) laser pulses, respectively, has been experimentally addressed by using time-gated optical emission spectroscopy of the ablated material and SEM analysis of the targets. The ablation process has also been analyzed by classical, molecular dynamics (MD) simulations, by using a Morse potential to describe the interaction between the atoms, and taking into account the electron heat diffusion contribution. The dependence of the ablation depth on laser fluence, as measured by SEM analysis, is in good agreement with the numerical simulations and is also well correlated with the optical emission yield of the expanding plume.     

SiC晶界薄膜的变电荷分子动力学模拟

基于迭代变电荷方法,用分子动力学模拟了SiC中的晶界薄膜．从原子尺度上模拟了不同的晶界薄膜的结构．观察到了晶粒与晶界薄膜间的电荷转移并且晶界薄膜的厚度与电荷转移有关．该结果提供了晶界存在空间电荷的直接证据,并证明静电作用与晶界薄膜的平衡厚度密切相关． 关键词： 分子动力学 变电荷 晶界薄膜     

团簇异构体生长概率的理论预测

由于实验表征团簇结构的困难,理论预测团簇在生长过程中各种异构体的生成概率便十分重要.考虑到分子动力学可模拟高温时各异构体的相对丰度,把问题归结为高温时所形成的各种异构体在室温环境中如何转化,并建立了统计力学模型,得到了异构体相互演化所需时间与环境温度的解析关系.以C₁₂团簇为例,在2500 K时所形成的最可几异构体在室温环境中需要超过10¹²年才能转化为最低势能异构体. 关键词： 团簇 异构体 统计物理 分子动力学