Atomic Velocity Distributions and Diffusing Behaviors in a Nano-Alloy Cluster

This paper studies the velocity distributions and diffusing behaviors of the atoms in a nano-alloy cluster. A series of ternary alloy clusters, Au₅₀Cu₂₅Ni₂₅, and binary alloy clusters, Cu_100−aNi_a are introduced in the molecular dynamics simulations. The velocity distributions of different types of atoms in both static and moving clusters are found to obey the Maxwell’s velocity distribution with the individual mass of atoms and the inner temperature of clusters. Furthermore, the velocity distribution of whole atoms of the cluster is obtained by synthesizing the velocity distributions of compositions in the cluster according to the proportions. The consistency of the atomic motions in a moving cluster is discussed by inspecting the backward velocities of atoms, which are correlated to the translational velocity and the inner temperature of the cluster. The diffusing behaviors of the atoms in a cluster are also investigated from a viewpoint of the interatomic interactions, i.e. Cu atoms enhance the activities of Ni atoms, or Ni atoms reduce that of Cu atoms in the alloy clusters.     

Large-Scale Molecular Dynamics Simulations of Energetic Ni Nanocluster Impact onto the Surface

On the atomic scale, Molecular Dynamics (MD) Simulation of Nano Ni cluster impact on Ni (100) substrate surface have been carried out for energies of E _a = 1–5 eV/atom and total energy of E _T = 195 eV (the total energy of cluster is E _T = nE _a, n is the number of cluster atoms) to understand quantitatively the interaction mechanisms between the cluster atoms and the substrate atoms. The many-body Embedded Atom Method (EAM) was used in this simulation. We investigated the maximum substrate temperature T _max and the time t _max within which this temperature is reached as a function of cluster sizes and the total energy E _T. The temperature T _max is linearly proportional to total cluster energy. For the constant energy per atom and for the cluster size increase, the correlated collisions rapidly transfers energy to the substrate, and the time t _max approached a constant value. For constant total energy the temperature T _max and the time t _max versus different cluster sizes was studied. We showed that the cluster implantation and sputtering atoms from the surface are affected by the cluster size and total kinetic energy of the clusters. Finally time dependence of the number N _dis of disordered atoms in the substrate was observed.     

Cluster cross sections from pickup measurements: are the established methods consistent?

Pickup of several molecules, H(2)O, HBr, and CH(3)OH, and Ar atoms on free Ar(N) clusters has been investigated in a molecular beam experiment. The pickup cross sections of the clusters with known mean sizes, ?≈ 150 and 260 were measured by two independent methods: (i) the cluster beam velocity decrease due to the momentum transfer of the picked up molecules to the clusters, and (ii) Poisson distribution of a selected cluster fragment ion as a function of the pickup pressure. In addition, the pickup cross sections were calculated using molecular dynamics and Monte Carlo simulations. The simulations support the results of the velocity measurements. On the other hand, the Poisson distributions yield significantly smaller cross sections, inconsistent with the known Ar(N) cluster sizes. These results are discussed in terms of: (i) an incomplete coagulation of guest molecules on the argon clusters when two or more molecules are picked up; and (ii) the fragmentation pattern of the embedded molecules and their clusters upon ionization on the Ar cluster. We conclude that the Poisson distribution method has to be cautiously examined, if conclusions should be drawn about the cluster cross section, or the mean cluster size ?, and the number of picked up molecules.     

Methodological approach to study energetic and structural properties of nanostructured cadmium sulfide by using ab-initio molecular dynamics simulations

A single wurtzite phase of cadmium sulfide cluster is investigated by ab-initio molecular dynamics simulations at different temperatures, ranging from 100 K to 600 K. In this study we propose a possible procedure to characterize the CdS quantum dots system by means of molecular dynamics calculations using a standard Car-Parrinello scheme. In order to ensure the accuracy of the numerical approach, preliminary calculations to test pseudopotentials, cutoff and box size on both single atoms systems and Cd–Cd, S–S, Cd–S dimers have been performed. Calculated binding energies and bond lengths are obtained in good agreement with experimental data. Subsequently, an uncapped CdS cluster with size below 2 nm, 48 atoms of cadmium and 48 atoms of sulfur, in a wurtzite geometry was structurally optimized to minimize internal stresses. The CdS cluster has been carefully characterized structurally at several temperatures up to T = 600 K. At the temperature of 340 K atomic diffusion on the surface allows the onset of a new stable atomic configuration.     

Growth of Co clusters on thin films Al2O3NiAl(100)

We present a scanning tunnel microscopy study of Co clusters grown through vapor deposition on Al(2)O(3) thin films over NiAl(100) at different coverages and temperatures. Formation of Co clusters was observed at 90, 300, 450, and 570 K. At the three lower temperatures, we find narrow cluster size distributions and the mean sizes (with a diameter of 2.6 nm and a height of 0.7 nm) do not change significantly with the coverage and temperature, until the clusters start to coalesce. Even on 3-4-nm-wide crystalline Al(2)O(3) strips where the deposited Co atoms are confined, the same features sustain. Only at 570 K the normal growth mode where the cluster size increases with the deposition coverage is observed, although the data are less conclusive. A simple modeling of kinetic surface processes on a strip confirms the normal growth mode, but fails to show a favored size unless additional energetic constraints are applied on the cluster sizes. Increasing Co coverages to cluster coalescence, a larger preferable size (mean diameter of 3.5 nm and height of 1.4 nm) appears for growth at 450 K. These two sizes are corroborated by morphology evolution of high Co coverages deposited at 300 K and annealed to 750 K, in which the coalescence is eliminated and the two preferable geometries appear and coexist.     

钯团簇碰撞沉积钯/银基板的分子动力学模拟

运用分子动力学模拟方法研究了常温下较大的钯团簇以不同初始速度撞击不同硬度基板的微观过程,着重分析了沉积形貌的变化、团簇的嵌入深度和原子的扩散程度、基板碰撞接触区域的温度演变以及碰撞过程中团簇与基板间的能量转化,获得了沉积过程中变形形貌、结构特征及能量转化随团簇尺寸、初始速度及基板材质的变化规律.并进一步探究了第二颗团簇以不同时刻沉积时前一团簇的变形形貌及基板接触区域温度变化的特点,发现短时间间隔下第二颗团簇的沉积更有利于团簇与基板的结合.     

Simulations of NaCl Aggregation from Solution: Solvent Determines Topography of Free Energy Landscape

The partition-enabled analysis of cluster histograms (PEACH) method is used to calculate the free energy surface of NaCl aggregation using cluster statistics from MD simulations of small systems (40–90 ions plus solvent) in four solvents. In all cases (pure methanol, pure water, and two methanol/water mixtures) NaCl clusters show a transition from amorphous to rocksalt structure with increasing cluster size. The crossover sizes, and the apparent kinetic barrier to ordering, increase with increasing water content. Implications for the proposed two-step mechanism of NaCl crystal nucleation (in which the ordered structure emerges from a large disordered cluster), and how this mechanism might depend on solvent and on degree of supersaturation, are discussed. In pure water, nonideal crowding effects that promote clustering are identified from systematic concentration-dependent deviations between simulation results and the PEACH model fit. In contrast, the ability of PEACH to fit aggregation statistics in mixed solvents is consistent with negligible interactions between ions in different clusters. © 2018 Wiley Periodicals, Inc.     

Self-assembling of chain molecules in low dimensions: A Monte Carlo study

We studied the self-assembling of linear chain molecules in insoluble monolayers due to attractive interactions. We used lattice Monte Carlo simulations in a two-dimensional system. The molecules consist of segments occupying adjacent lattice sites. The head segments are confined to move along a line whereas the chain segments can arrange in a plane above the heads. Only one interaction parameter is applied. At high densities and small interaction energy the system shows percolation behavior. At moderate and small densities it can be characterized by a monotonous cluster size distribution. Self-assembling occurs at small densities for strong attractive interactions. The corresponding cluster size distributions indicate preferred cluster sizes which depend upon density and interaction strength. With increasing density the clusters grow. The internal cluster structure depends on the cluster size and the interaction parameter. The clusters tend to minimize their total energy. Molecules at cluster margins contribute less to the cluster energy and are mainly disordered. They cause that the cluster properties strongly depend on the cluster size. Large clusters only have minimum energy if the molecules in the cluster are in stretched-out conformation. With decreasing interaction strength the clusters get disordered thereby producing less energy-minimized domain boundaries.     

Development of a molecular dynamics-based coalescence model for DSMC simulations of ammonia condensate flows

A coalescence model for homogeneous condensation of ammonia in supersonic expansions to vacuum has been developed using molecular dynamics trajectory calculations. The MD calculations show that the sticking probability increases as the ammonia cluster size increases or the cluster temperature decreases. In addition, the sensitivity of the sticking probability to cluster size decreases as the temperature decreases. Comparison of the Ashgriz-Poo semiempirical coalescence model with MD simulations show that for cluster sizes larger than 100 the former model may be used. To model homogeneous nucleation in an ammonia jet, direct simulation Monte Carlo (DSMC) simulations were performed for different stagnation pressure conditions using the MD simulation outcomes for smaller cluster-cluster collisions and the Ashgriz-Poo model for cluster sizes larger than 100. We found that, by including the combined coalescence model, the average cluster sizes and size distributions predicted by DSMC agree reasonably well with experiment.     

Characterization of large supported metal clusters by optical spectroscopy

Small sodium and silver particles were generated on dielectric substrates like LiF, quartz and sapphire under ultrahigh vacuum conditions. The optical transmission spectra of the clusters were measured as a function of cluster size and shape, for low and high substrate temperatures as well as for s- and p- polarization of the incident light. Excitation of dipolar surface plasmon oscillations in the directions normal and parallel to the substrate surface could be identified. Furthermore, optical spectra for Na and Ag clusters were calculated with the classical Mie theory. The measured spectra vary strongly if the experimental conditions are changed and can be exploited, for example, to characterize the particles with regard to their size and shape. In particular, the axial ratio of the spheroidal clusters could be determined. Its value is considerably different for the two investigated metals and depends on the substrate material. Furthermore, the temperature of the substrate has a pronounced influence on the shape of the particles. At low temperature of T=100 K two-dimensional island growth is predominant. The particles extend only little in the direction perpendicular to the surface and coalesce readily at small coverage of metal atoms. In contrast, the clusters are truly three-dimensional at T=300 K. At this stage, sodium particles still exhibit a rather small axial ratio whereas silver clusters appear almost spherical. Thus, measurements of the optical spectra permit direct in situ monitoring of cluster growth during the nucleation of adsorbed atoms and of temperature induced shape variations. In addition to investigations of the shape of the particles, the quadrupolar surface plasmon mode was observed for Ag clusters.     

Melting and freezing characteristics and structural properties of supported and unsupported gold nanoclusters

Molecular dynamics simulations in conjunction with MEAM potential models have been used to study the melting and freezing behavior and structural properties of both supported and unsupported Au nanoclusters within a size range of 2 to 5 nm. In contrast to results from previous simulations regarding the melting of free Au nanoclusters, we observed a structural transformation from the initial FCC configuration to an icosahedral structure at elevated temperatures followed by a transition to a quasimolten state in the vicinity of the melting point. During the freezing of Au liquid clusters, the quasimolten state reappeared in the vicinity of the freezing point, playing the role of a transitional region between the liquid and solid phases. In essence, the melting and freezing processes involved the same structural changes which may suggest that the formation of icosahedral structures at high temperatures is intrinsic to the thermodynamics of the clusters, rather than reflecting a kinetic phenomenon. When Au nanoclusters were deposited on a silica surface, they transformed into icosahedral structures at high temperatures, slightly deformed due to stress arising from the Au-silica interface. Unlike free Au nanoclusters, an icosahedral solid-liquid coexistence state was found in the vicinity of the melting point, where the cluster consisted of coexisting solid and liquid fractions but retained an icosahedral shape at all times. These results demonstrated that the structural stability in the structures of small Au nanoclusters can be enhanced through interaction with the substrate. Supported Au nanoclusters demonstrated a structural transformation from decahedral to icosahedral motifs during Au island growth, in contrast to the predictions of the minimum-energy growth sequence: icosahedral structures appear first at very small cluster sizes, followed by decahedral structures, and finally FCC structures recovered at very large cluster sizes. The simulations also showed that island shapes are strongly influenced by the substrate, more specifically, the structural characteristic of a Au island is not only a function of size, but also depends on the contact area with the surface, which is controlled by the wetting of the cluster to the substrate.     

Molecular dynamic simulations of atom-cluster collision processes

Monomer-cluster collisions of Lennard-Jones argon atoms have been studied using molecular dynamics simulation for target cluster sizes of 2, 3, 4, 5, 10, and 20 atoms. Capture probability of monomers by clusters and the lifetimes of the resulting clusters have been calculated as a function of impact parameter and the total energy of the target cluster. Cluster lifetime is further integrated over all impact parameters to obtain the average lifetime for each cluster size and energy. The average lifetime of the smallest aggregates is shown to be short compared to the collision time between monomers and clusters unless the vapor is highly supersaturated. The formation probability of a new cluster decreases steeply if a minimum lifetime is required for the cluster.     

过渡金属原子簇离子结合能的光离解测量

本实验以激光离解的方式, 较系统地测量了钴、铌及其他一些过渡金属原子簇阳离子的结合能. 实验结果表明: 金属原子簇离子的结合能与其成族原子数关系不大, 而明显地受其组分的影响. 本文还探讨了金属原子簇离子的光离解机理及其几何结构模型。     

The role of dimers in evaporation of small argon clusters

Evaporation of small Lennard-Jones argon clusters has been studied using molecular dynamic simulations. An extensive library of clusters with 4, 5, 6, 11, and 21 atoms has been obtained from an earlier study. Analysis of the evaporation properties of the clusters indicate, that the fraction of dimer evaporations of all evaporation events increases with the total energy of the cluster. The fraction of evaporated dimers from clusters with a constant lifetime is independent of the cluster size for short-lived clusters and increases with cluster size for long-lived clusters. Only a few percent of the clusters which are long lived enough to participate in vapor-liquid nucleation decay by emitting dimers. The mean cluster lifetime as a function of total energy shows the same exponentially decreasing trend for monomer and dimer evaporation channels. The fraction of trimer evaporations is found to be vanishingly small.     

Computational spectroscopy of carbon monoxide isotopomers in helium clusters

The rotational excitation spectrum, including the vibrational shift of the rotational band, of several CO isotopomers solvated in He clusters has been calculated. Reptation quantum Monte Carlo simulations are used in conjunction with an accurate He-CO potential energy surface, which quantitatively describes the rovibrational spectrum of the binary complex. Our simulations, when compared with number-selective infrared spectra taken for different isotopomers, help discriminate among the alternative assignments proposed for cluster sizes around 15 He atoms. The origin of the vibrational band has a red shift that is nearly linear with the cluster size within the first solvation shell and is almost constant up to the largest cluster studied, well beyond completion of the second solvation shell. A blue upturn at even larger sizes would be needed to attain the nanodroplet limit, as recently estimated from the isotopic dependence of the measured R(0) transitions.     

Ion calorimetry: Using mass spectrometry to measure melting points

Calorimetry measurements have been used to probe the melting of aluminum cluster cations with 63 to 83 atoms. Heat capacities were determined as a function of temperature (from 150 to 1050 K) for size-selected cluster ions using an approach based on multicollision-induced dissociation. The experimental method is described in detail and the assumptions are critically evaluated. Most of the aluminum clusters in the size range examined here show a distinct peak in their heat capacities that is attributed to a melting transition (the peak is due to the latent heat). The melting temperatures are below the bulk melting point and show enormous fluctuations as a function of cluster size. Some clusters (for example, n = 64, 68, and 69) do not show peaks in their heat capacities. This behavior is probably due to the clusters having a disordered solid-like phase, so that melting occurs without a latent heat.     

Molecular Dynamics Studies of the Kinetics of Phase Changes in Clusters IV： Crystal Nucleation from Molten （NaCl）256 and （NaCl）500 Clusters

Molecular dynamics computer simulation based on the Born-Mayer-Huggins potential function has been carried out to study the effects of duster size and temperature on the nucleation rate of sodium chloride dusters in the temperature range of 580 K to 630 K. Clusters with 256 and 500 NaCl molecules have been studied and the results have been compared with those obtained from 108 molecule dusters. The melting point (MP) of the clusters were observed to increase with the size of the clusters and can be well described by a linear equation MP =1107(37)-1229(23)N^-1/3(N is the number of molecules in the duster).The nucleation rate was found to decrease with increasing the duster size or temperature. Various nucleation theories have been used to interpret the nucleation rates obtained from this molecular dynamics simulation. It is possible to use a constant diffuse interface thickness to interpret the nucleation rate from the diffuse interface theory in the temperature range of this study. However, the interfacinl free energy estimated from classical nucleation theory and diffuse interface theory increases too fast with increasing the temperature while that from Gran-Gunton theory does not change with changing temperatures.The sizes of critical nuclei estimated from all the theories are smaller than those estimated from our simulations.     

Chirally directed formation of nanometer-scale proline clusters

Ion mobility measurements, combined with molecular mechanics simulations, are used to study enantiopure and racemic proline clusters formed by electrospray ionization. Broad distributions of cluster sizes and charge states are observed, ranging from clusters containing only a few proline units to clusters that contain more than 100 proline units (i.e., protonated clusters of the form [xPro + nH](n+) with x = 1 to >100 and n = 1-7). As the sizes of clusters increase, there is direct evidence for nanometer scale, chirally induced organization into specific structures. For n = 4 and 5, enantiopure clusters of approximately 50 to 100 prolines assemble into structures that are more elongated than the most compact structure that is observed from the racemic proline clusters. A molecular analogue, cis-4-hydroxy-proline, displays significantly different behavior, indicating that in addition to the rigidity of the side chain ring, intermolecular interactions are important in the formation of chirally directed clusters. This is the first case in which assemblies of chirally selective elongated structures are observed in this size range of amino acid clusters. Relationships between enantiopurity, cluster shape, and overall energetics are discussed.     

Argon cluster-ion sputter yield: Molecular dynamics simulations on silicon and equation for estimating total sputter yield

Argon Gas Cluster-Ion Beam sources have become widely-used on X-ray photoelectron spectroscopy (XPS) and secondary ion mass spectrometry (SIMS) instruments in recent years, but there is little reference data on sputter yields in the literature as yet. Total sputter yield reference data is needed in order to calibrate the depth scale, of XPS or SIMS depth profiles. We previously published a semi-empirical ‘Threshold’ equation for estimating cluster total sputter yield from the energy-per-atom of the cluster and the effective monatomic sputter threshold of the material. This has been shown to agree extremely well with sputter yield measurements on a range or organic and inorganic materials for clusters of around a thousand atoms. Here we use the molecular dynamics (MD) approach to explore a wider range of energy and cluster size than is easy to do experimentally to high precision. We performed MD simulations using the ‘Large-scale Atomic/Molecular Massively Parallel Simulator’ (LAMMPS) parallel MD code on high-performance computer (HPC) systems. We performed 1150 simulations of individual collisions with a silicon (100) surface as an archetypal inorganic substrate, for cluster sizes between 30 and 3000 argon atoms and energies in the range 5 to 40 eV per atom. This corresponds to the most important regime for experimental cluster depth-profiling in SIMS and XPS. Our MD results show a dependence on cluster size as well as energy-per-atom. Using the exponent previously suggested by Paruch et al., we modified the Threshold model equation published previously to take this into account. The modified Threshold equation fits all our MD results extremely well, building on its success in fitting experimental sputter yield measurements. This work is submitted to the volume dedicated to Dr. Martin P Seah, MBE, who was a great influence on the early career of one of the authors (PJC) and who himself made many valuable contributions to the literature on sputtering as it relates to surface and interface analysis.