Molecular dynamic simulation of lubricant spreading： effect from the substrate and endbead

Molecular dynamic simulations based on a coarse-grained, bead-spring model are adopted to investigate the spreading of both nonfunctional and functional perfluoropolyether (PFPE) on solid substrates. For nonfunctional PFPE, the spreading generally exhibits a smooth profile with a precursor film. The spreading profiles on different substrates are compared, which indicate that the bead-substrate interaction has a significant effect on the spreading behaviour, especially on the formation of the precursor film. For functional PFPE, the spreading generally exhibits a complicated terraced profile. The spreading profiles with different endbeads are compared, which indicate that the endbead-substrate interaction and the endbead--endbead interaction, especially the latter, have a significant effect on the spreading behaviour.     

磁盘润滑膜全氟聚醚的分子动力学模拟研究

纳米润滑膜全氟聚醚(perfluoropolyether,简称PFPE)在固体表面的结构和迁移特性对于计 算机磁盘磁头驱动系统的可靠性具有重要的作用. 采用基于粗粒珠簧模型的分子动力学模拟 方法,考察了不同壁面和端基极性对于PFPE膜静态特性（如分子构型、单体密度分布、端基 密度分布）以及动态特性（如自扩散系数）的影响. 静态特性的模拟结果表明,非极性PFPE 膜在壁面附近具有单层厚度为一个单体直径的层状结构,而极性PFPE膜则具有复杂的层状结 构. 动态特性的模拟结果表明,PFPE膜的扩散能力因壁面作用而增强并因端基极性作用而减 弱.关键词：全氟聚醚膜分子动力学模拟薄膜润滑固液界面结构     

The spreading behaviour of perfluoropolyether droplets on solid surfaces

The spread of perfluoropolyether （PFPE） droplets on solid surfaces has been measured from the top-down view through a microscope system. Effects of substrates, molecular weight and end-group functionality on spreading of the PFPE droplets have been studied experimentally and the results were compared with those by molecular dynamics （MD） simulations. Silicon wafer and diamond-like carbon （DLC） substrates were used to study the effect of substrates on spreading. Two types of PFPE, Z-dol and Z-tetraol, with the same chain structure and various molecular weights （2000 and 4000 g/mol） were employed in experiments. Effect of molecular weight has been investigated through comparing the spreading of Z-dol 2000 and Z-dol 4000, and it is found that the increase of molecular weight will decrease the mobility of PFPE. Comparison between spreading of Z-dol and Z-tetraol of the same molecular weight proved that functional end group plays a significant role on the spreading of PFPE, which confirmed the MD simulation results.     

纳米级润滑膜分子排列取向的拉曼光谱表征技术

利用激光拉曼散射技术,对剪切作用下受限于钢球与石英盘之间的纳米级液晶5 CB的分子排列取向进行研究. 结果表明,在特定的实验条件下,可以得到高信噪比的纳米级润滑膜的拉曼散射信号(20∶1). 同时发现,当激光偏振方向与剪切运动方向平行(垂直)时,所得拉曼信号强度达到最大(小)值,表明纳米级液晶5 CB分子在剪切诱导作用下,沿剪切运动方向趋于定向排列. 另外,当钢球与石英盘之间的剪切速度逐渐增大时,受限的纳米级液晶5 CB的拉曼信号强度也逐渐增大. 最后,利用根据相对光强干涉原理研制的纳米膜厚测量仪对纳米级关键词：薄膜润滑分子排列取向拉曼散射     

Modelling the effect of heterogeneous vaccination on metapopulation epidemic dynamics

Vaccination as an epidemic control strategy has a significant effect on epidemic spreading. In this paper, we propose a novel epidemic spreading model on metapopulation networks to study the impact of heterogeneous vaccination on epidemic dynamics, where nodes represent geographical areas and links connecting nodes correspond to human mobility between areas. Using a mean-field approach, we derive the theoretical spreading threshold revealing a non-trivial dependence on the heterogeneity of vaccination. Extensive Monte Carlo simulations validate the theoretical threshold and also show the complex temporal epidemic behaviours above the threshold.     

TiAl合金薄膜在冷却过程中结构变化的原子尺度计算研究

采用基于嵌入原子势模型的分子动力学方法模拟了熔融TiAl合金薄膜分别在急冷降温和连续降温两个冷却过程中微观结构的变化.通过应用原子平均能量、对分布函数和键对分析技术等对薄膜内原子局域结构进行的分析. 研究表明,在降温过程中,薄膜内局域结构出现分阶段变化的特点.在急冷降温过程中,薄膜结构的变化可以分为3个阶段,而在连续降温过程中,薄膜内局域结构的变化分为2个阶段.关键词：分子动力学TiAl合金薄膜计算机模拟     

Time-dependent thermodynamic properties of the Ising model from damage spreading

The relationship between damage spreading and static thermodynamic properties in the Ising model developed by Coniglioet al. is here extended to include time-dependent thermodynamic quantities. We exploit this new result to measure the time-dependent spin correlation function from damage spreading in the Ising model with heat bath and Glauber dynamics. Until now, only static thermodynamic quantities have been correctly determined from damage spreading, and even then, only with heat bath dynamics. We also show that there are significant differences between the kinetics of damage spreading as found in heat bath and Glauber dynamics.     

Damage spreading in a single-component irreversible reaction process: Dependence of the system's immunity on the Euclidean dimension

The spreading of a globally distributed damage, created in the stationary regime, is studied in a single-component irreversible reaction process, i.e., the BK model [Browne and Kleban,Phys. Rev. A 40, 1615 (1989)]. The BK model describes one variant of the A+AA₂ reaction process on a lattice in contact with a reservoir of A species. The BK model has a single parameter, namely the rate of arrival of A species to the lattice (Y). The model, exhibits an irreversible phase transition between a stationary reactive state with production of A₂ species and a poisoned state with the lattice fully covered by A species. The transition takes place at critical points (Y _C ) which solely depend on the Euclidean dimensiond. It is found that the system is immune ford=1 andd=2, in the sense that even 100% of initial damage is healed within a finite healing period (T _H ). Within the reactive regime,T _H diverges when approachingY _C according toT _H (Y _C –Y)^–, with 1.62 and 1.08 ford=1 andd=2, respectively. Ford=3 a frozen-chaotic transition is found close toY _s 0.4125, i.e., well inside the reactive regime 0YY _C 0.4985. Just atY _S the damageD(t) heals according toD(t) t ^–, with 0.71. For the frozen-chaotic transition atd=3 the order parameter critical exponent 0.997 is determined.     

Simulation of grain boundary effect on characteristics of ZnO thin film transistor by considering the location and orientation of grain boundary

The grain boundaries (GBs) have a strong effect on the electric properties of ZnO thin film transistors (TFTs). A novel grain boundary model was developed to analyse the effect. The model was characterized with different angles between the orientation of the grain boundary and the channel direction. The potential barriers formed by the grain boundaries increase with the increase of the grain boundary angle, so the degradation of the transistor characteristics increases. When a grain boundary is close to the drain edge, the potential barrier height reduces, so the electric properties were improved.     

MgO（001）表面上沉积MgO薄膜过程的分子动力学模拟

采用分子动力学方法模拟了MgO分子连续沉积于MgO（001）表面上的薄膜生长过程,分析了衬底温度和分子入射能对MgO分子在衬底表面上的扩散能力以及对衬底表面覆盖率的影响.模拟结果表明,随着衬底温度的升高,在衬底表面上沉积的MgO分子扩散能力增强,MgO薄膜层中空位缺陷变少.低温下,分子入射能的增大有助于提高衬底表面覆盖率;高温下,表面覆盖率随入射能增大到3.0 eV时达到最大值,入射能继续增大,表面覆盖率减小.关键词：MgO薄膜生长分子动力学计算机模拟表面扩散     

Magnetization reversal process in Fe/Si (001) single-crystalline film investigated by planar Hall effect

A planar Hall effect(PHE) is introduced to investigate the magnetization reversal process in single-crystalline iron film grown on a Si(001) substrate.Owing to the domain structure of iron film and the characteristics of PHE,the magnetization switches sharply in an angular range of the external field for two steps of 90° domain wall displacement and one step of 180°domain wall displacement near the easy axis,respectively.However,the magnetization reversal process near the hard axis is completed by only one step of 90° domain wall displacement and then rotates coherently.The magnetization reversal process mechanism near the hard axis seems to be a combination of coherent rotation and domain wall displacement.Furthermore,the domain wall pinning energy and uniaxial magnetic anisotropy energy can also be derived from the PHE measurement.     

The effect of deposition temperature on the intermixing and microstructure of Fe/Ni thin film

The physical vapour deposition of Ni atoms on α-Fe(001) surface under different deposition temperatures were simulated by molecular dynamics to study the intermixing and microstructure of the interfacial region. The results indicate that Ni atoms hardly penetrate into Fe substrate while Fe atoms easily diffuse into Ni deposition layers. The thickness of the intermixing region is temperature-dependent, with high temperatures yielding larger thicknesses. The deposited layers are mainly composed of amorphous phase due to the abnormal deposition behaviour of Ni and Fe. In the deposited Ni-rich phase, the relatively stable metallic compound B₂ structured FeNi is found under high deposition temperature conditions.     

Interface mixing behaviour of Lennard–Jones FCC (100) thin film

A modified Monte Carlo method combined with quenched molecular dynamics simulation is used to determine mixing energetics and concentration profiles at interface for systems containing mono-and bilayers of adatoms adsorbed on FCC (100) crystal surface. The systems under consideration are constructed via Lennard–Jones potential at temperatures near 0 K. For systems with monolayer of adatoms, intermixing at the interface becomes preferable with increasing magnitude of the potential well-depth ratio of adatom to substrate atom. The increasing tendency of intermixing is linearly enhanced when the adatom becomes smaller than the substrate atom, otherwise, the intermixing trend is non-linear and weaker. For systems with bilayers of adatoms, complex development of concentration profile is observed along with increasing magnitude of the potential well-depth ratio and atomic size difference between adatom and substrate atom. This behaviour is related to the interplay between contributions of asymmetric bond interaction and relaxation to minimise the total energy of the system.     

Molecular dynamics study of the mechanical characteristics of Ni/Cu bilayer using nanoindentation

<正>In the present work,a three-dimensional molecular dynamics simulation is carried out to perform the nanoindentation experiment on Ni single crystal.The substrate indenter system is modeled using hybrid interatomic potentials including the many-body potential embedded atom method(EAM),and two-body morse potential.To simulate the indentation process,a spherical indenter(diameter = 80 A,1 A=0.1 nm) is chosen.The results show that the mechanical behaviour of a monolithic Ni is not affected by crystalline orientation.To elucidate the effect of a heterogeneous interface, three bilayer interface systems are constructed,namely Ni(100)/Cu(111),Ni(110)/Cu(111),and Ni(111)/Cu(111).The simulations along these systems clearly describe that mechanical behaviour directly depends on the lattice mismatch. The interface with the smaller mismatch between the specified crystal planes is proved to be harder and vice versa.To describe the relationship between film thickness and interface effect,we choose various values of film thickness ranging from 20 A to 50 A to perform the nanoindentation experiment.It is observed that the interface is significant only for the relatively small thickness of film and the separation between interface and the indenter tip.It is shown that with the increase in film thickness,the mechanical behaviour of the film shifts more toward that of monolithic material.     

Theoretical study of the hydrogen abstraction reaction of CH3OH with NCO

A direct dynamics method is employed to study the mechanism and kinetics of the hydrogen abstraction reaction of CH₃OH with NCO. The optimized geometries and frequencies of the stationary points and the minimum-energy paths (MEPs) are obtained at the MP2/6-311G(d,p) level. In order to obtain more accurate potential energy surface (PES) information and provide more credible energy data for kinetic calculation, the single-point energies along the MEPs are further computed at QCISD(T)/6-311+G(d,p) and G3MP2 levels. The rate constants for two channels, the methyl-H abstraction channel and hydroxyl-H abstraction channel, are calculated by canonical variational transition state theory (CVT) with small-curvature tunneling (SCT) contributions over the wide temperature region 220–1500?K. The theoretical overall rate constants are in good agreement with the available experimental data. For the title reaction, the methyl-H abstraction channel is dominant, while the hydroxyl-H abstraction channel is negligible over the whole temperature region.     

Oxidation behavior with quantum dots formation from amorphous GaAs thin films

ABSTRACT

We investigated the oxidation behaviour of an amorphous GaAs thin film deposited onto a micro/nanotextured Si surface by an electron beam. After the deposited film was exposed to air, microcrystallites were formed with octahedral cubic and monoclinic structures of arsenic oxides. Short time exposure after thin film deposition showed the formation of cubic arsenolite while long time exposure showed the formation of monoclinic claudetites as well as cubic arsenolites. These oxide microcrystallites at the GaAs thin film surface disappeared after the sample annealing process. However, the amorphous GaAs thin film included high-density GaAs nanodots. From UV and inverse photoemission spectroscopies, the thin film showed n-type band structure with an energy gap of 2.73?eV. Photoluminescence measurement showed an emission peak at (450–513)?nm with the energy of (2.41–2.75)?eV corresponding to dot size of (4.1–4.5)?nm.     

Annealing effect of NiO/Co90Fe10 thin films: From bilayer to nanocomposite

Exchange-biased bilayers are widely used in the pinned layers of spintronic devices. While magnetic field annealing (MFA) was routinely engaged during the fabrication of these devices, the annealing effect of NiO/CoFe bilayers is not yet reported. In this paper, the transition from NiO/Co₉₀Fe₁₀ bilayer to nanocomposite single layer was observed through rapid thermal annealing at different temperatures under magnetic field. The as-deposited and low-temperature (<623 K) annealed samples had rock salt (NiO) and face center cubic (Co₉₀Fe₁₀) structures. On the other hand, annealing at 623 K and 673 K resulted in nanocomposite single layers composed of oxides (matrix) and alloys (precipitate), due to grain boundary oxidization and strong interdiffusion in the NiO/CoFe and CoFe/SiO₂ interfaces. The structural transition was accompanied by the reduction of grain sizes, re-ordering of crystallites, incensement of roughness, and reduction of Ni²⁺. When measured at room temperature, the bilayers exhibited soft magnetism with small room-temperature coercivity. The nanocomposite layers exhibited an enhanced coercivity due to the changes in the magnetization reversal mechanism by pinning from the oxides. At 10 K, the increased antiferromagnetic anisotropy in the NiO resulted in enhanced coercivity and exchange bias in the bilayers. The nanocomposites exhibited weaker exchange bias compared with the bilayers due to frustrated interfacial spins. This investigation on how the magnetic properties of exchange-biased bilayers are influenced by magnetic RTA provides insights into controlling the magnetization reversal properties of thin films.     

Molecular dynamics study on water self-diffusion in aqueous mixtures of methanol,ethylene glycol and glycerol: investigations from the point of view of hydrogen bonding

Molecular dynamics simulations have been performed to investigate the aqueous binary mixtures of alcohols, including methanol, ethylene glycol (EG) and glycerol of molalities ranging from 1 to 5 m at the temperatures of 273, 288 and 298 K, respectively. The primary purpose of this paper is to investigate the mechanism of water self-diffusion in water-alcohol mixtures from the point of view of hydrogen bonding. The effects of temperature and concentration on water self-diffusion coefficient are evaluated quantitatively in this work. Temperature and concentration to some extent affect the hydrogen bonding statistics and dynamics of the binary mixtures. It is shown that the self-diffusion coefficient of water molecules decreases as the concentration increases or the temperature decreases. Moreover, calculations of mean square displacements of water molecules initially with different number n of H-bonds indicate that the water self-diffusion coefficient decreases as n increases. We also studied the aggregation of alcohol molecules by the hydrophobic alkyl groups. The largest cluster size of the alkyl groups clearly increases as the concentration increases, implying the emergence of a closely connected network of water and alcohols. The clusters of water and alcohol that interacted could block the movement of water molecules in binary mixtures. These findings provide insight into the mechanisms of water self-diffusion in aqueous binary mixtures of methanol, EG and glycerol.     

非磁性材料磁电阻效应的研究进展

探讨了非磁性材料与磁性材料磁电阻效应的异同，分析了几种有代表性的非磁性材料的巨大磁电阻效应，概述了目前非磁材料磁电阻效应的研究和应用状况。