Constant pressure molecular dynamics simulations of the 2D r−12 system: Comparison with isochores and isotherms

On morphologies of gold nanoparticles grown from molecular dynamics simulation

The authors use a newly fitted gold embedded atom method potential to simulate the initial nucleation, coalescence, and kinetic growth process of vapor synthesized gold nanoparticles. Overall the population statistics obtained in this work seemed to mirror closely recent experimental HREM observations by Koga and Sugawara [Surf. Sci. 529, 23 (2003)] of inert gas synthesized nanoparticles, in the types of nanoparticles produced and qualitatively in their observance ratio. Our results strongly indicated that early stage coalescence (sintering) events and lower temperatures are the mainly responsible for the occurrence of the Dh and fcc based morphologies, while "ideal" atom by atom growth conditions produced the Ih morphology almost exclusively. These results provide a possible explanation as to why the Dh to Ih occurrence ratio increases as a function of nanoparticle size as observed by Koga and Sugawara.

On the computational calculation of surface free energies for the disordered semihexagonal reconstructed Au(100) surface

Previously we developed a general method for calculating the free energy of any surface constrained to a distinct surface excess number/density. In this paper we show how to combine a range of such surfaces, whose free energies have been calculated, to produce an ad hoc semigrand canonical ensemble of surfaces from which ensemble surface properties can be calculated, including the ensemble surface free energy. We construct such an ensemble for the disordered Au(100) semihexagonal reconstructed surface using a Glue model potential at 1000 K and calculate the ensemble surface free energy to be 0.088 18 eVA(2). The ensemble average surface lateral density was found to be 1.375 (with respect to the bulk), which is in agreement with previous grand canonical Monte Carlo studies.

Application of the constrained fluid λ-integration path to the calculation of high temperature Au(110) surface free energies

Recently a method termed constrained fluid lambda-integration was proposed for calculating the free energy difference between bulk solid and liquid reference states via the construction of a reversible thermodynamic integration path; coupling the two states in question. The present work shows how the application of the constrained fluid lambda-integration concept to solid/liquid slab simulation cells makes possible a generally applicable computer simulation methodology for calculating the free energy of any surface and/or surface defect structure, including surfaces requiring variations in surface atom or density number, such as the (1 x 5) Au(100) or (1 x 2) missing row Au(110) reconstructed surfaces or excess adatom/vacancy/step populated surfaces. We evaluate the methodology by calculating the free energy of various disordered high temperature Au(110) embedded atom method surfaces constrained to differing excess surface atom numbers [including those corresponding to the (1 x 2) missing row reconstructed surface] and obtained the interesting result that at 1000 K (as distinct from lower temperatures) the free energy difference between these surfaces is reduced to zero; a result which is consistent with an expected order-disorder phase transition for the Au(110) surface at such high temperatures.

“Exact” surface free energies of iron surfaces using a modified embedded atom method potential and λ integration

Previously a new universal lambda-integration path and associated methodology was developed for the calculation of "exact" surface and interfacial free energies of solids. Such a method is in principle applicable to any intermolecular potential function, including those based on ab initio methods, but in previous work the method was only tested using a relatively simple embedded atom method iron potential. In this present work we apply the new methodology to the more sophisticated and more accurate modified embedded atom method (MEAM) iron potential, where application of other free- energy methods would be extremely difficult due to the complex many-body nature of the potential. We demonstrate that the new technique simplifies the process of obtaining "exact" surface free energies by calculating the complete set of these properties for the low index surface faces of bcc and fcc solid iron structures. By combining these data with further calculations of liquid surface tensions we obtain the first complete set of exact surface free energies for the solid and liquid phases of a realistic MEAM model system. We compare these predictions to various experimental and theoretical results.

Constrained fluid λ-integration: Constructing a reversible thermodynamic path between the solid and liquid state

A novel lambda-integration path is proposed for calculating the Gibbs free energy difference between any arbitrary solid and liquid state needed for the location of melting lines. This technique involves reversibly forcing a liquid state to a solid state across the phase transition along a nonphysical path, thermodynamically coupling the two states directly. The process eliminates the need for coupling to idealized reference states as is presently performed and hence simplifies the location of phase transitions for computer simulation systems. More specifically the path involves a three stage process, whereby, initially a liquid state is transformed to a weakly attractive fluid using linear lambda-integration scaling of the intermolecular potential. In the second stage, the resulting fluid is then constrained to the required solid configurational phase space via the insertion of a periodic lattice of 3D Gaussian wells. The final stage involves reversing to full strength the main intermolecular potential while gradually turning off the constraining 3D Gaussian lattice finally resulting in a stable (or metastable) solid state. Each stage was found to be completely reversible and the resulting change in free energy was thermodynamically integrable. The methodology is demonstrated and validated by calculating solid-liquid coexistence points using the new technique and comparing to those in present literature for the truncated and shifted Lennard-Jones system. The results are found to be in good agreement. The new method is not limited to melting phase transitions and is readily applicable to any simulation methodology, simulation cell size and/or intermolecular potential including ab initio methods.

Diffusion of hard disks and rodlike molecules on surfaces

We study the submonolayer diffusion of hard disks and rodlike molecules on smooth surfaces through numerical simulations and theoretical arguments. We concentrate on the behavior of the various diffusion coefficients as a function of the two-dimensional (2D) number density rho in the case where there are no explicit surface-particle interactions. For the hard disk case, we find that while the tracer diffusion coefficient D(T)(rho) decreases monotonically up to the freezing transition, the collective diffusion coefficient D(C)(rho) is wholly determined by the inverse compressibility which increases rapidly on approaching freezing. We also study memory effects associated with tracer diffusion, and present theoretical estimates of D(T)(rho) from the mode-mode coupling approximation. In the case of rigid rods with short-range repulsion and no orientational ordering, we find behavior very similar to the case of disks with the same repulsive interaction. Both D(T)(rho) and the angular diffusion coefficient D(R)(rho) decrease with rho. Also in this case D(C)(rho) is determined by inverse compressibility and increases rapidly close to freezing. This is in contrast to the case of flexible chainlike molecules in the lattice-gas limit, where D(C)(rho) first increases and then decreases as a function of the density due to the interplay between compressibility and mobility.

Waterlike anomalies for core-softened models of fluids: two-dimensional systems

We use molecular-dynamics simulations in two dimensions to investigate the possibility that a core-softened potential can reproduce static and dynamic anomalies found experimentally in liquid water: (i) the increase in specific volume upon cooling, (ii) the increase in isothermal compressibility upon cooling, and (iii) the increase in the diffusion coefficient with pressure. We relate these anomalies to the shape of the potential. We obtain the phase diagram of the system and identify two solid phases: a square crystal (high-density phase) and a triangular crystal (low-density phase). We also discuss the relation between the anomalies observed and the polymorphism of the solid. Finally, we compare the phase diagram of our model system with experimental data, noting especially the line of temperatures of maximum density, the line of pressures of maximum diffusion constant, and the line of temperatures of minimum isothermal compressibility.

Stability and structure of a supercooled liquid mixture in two dimensions

The structural and thermodynamic properties of a two-dimensional binary mixture of soft discs are reported over a range of temperatures down to large supercoolings using constant NPT molecular dynamics simulations. It is shown that the four orders of magnitude increase in the structural relaxation time is not accompanied by any significant increase in translational or orientational order. The phase diagram, calculated in the temperature/composition plane using thermodynamic integration, exhibits a deep eutectic point that is responsible for stabilizing the amorphous state. Voronoi analysis of the low-temperature ground state reveals a structure characterized by a network of linear arrays of fivefold and sevenfold sites. The heat capacity C(P) exhibits an asymmetric peak with a maximum at T(*)=0.55. It is argued that the initial rapid drop in C(P) for T*<0.55 is an equilibrium result and, hence, the peak in the heat capacity corresponds to the existence of an "enthalpy gap" with a characteristic temperature of T* approximately 0.35. This gap results from a minimum volume change associated with an anharmonic fluctuation.