Magic numbers for classical Lennard‐Jones cluster heat capacities

Solution of disordered microphases in the Bethe approximation

The periodic microphases that self-assemble in systems with competing short-range attractive and long-range repulsive (SALR) interactions are structurally both rich and elegant. Significant theoretical and computational efforts have thus been dedicated to untangling their properties. By contrast, disordered microphases, which are structurally just as rich but nowhere near as elegant, have not been as carefully considered. Part of the difficulty is that simple mean-field descriptions make a homogeneity assumption that washes away all of their structural features. Here, we study disordered microphases by exactly solving a SALR model on the Bethe lattice. By sidestepping the homogenization assumption, this treatment recapitulates many of the key structural regimes of disordered microphases, including particle and void cluster fluids as well as gelation. This analysis also provides physical insight into the relationship between various structural and thermal observables, between criticality and physical percolation, and between glassiness and microphase ordering.

A single-walker approach for studying quasi-nonergodic systems

The jump-walking Monte-Carlo algorithm is revisited and updated to study the equilibrium properties of systems exhibiting quasi-nonergodicity. It is designed for a single processing thread as opposed to currently predominant algorithms for large parallel processing systems. The updated algorithm is tested on the Ising model and applied to the lattice-gas model for sorption in aerogel at low temperatures, when dynamics of the system is critically slowed down. It is demonstrated that the updated jump-walking simulations are able to produce equilibrium isotherms which are typically hidden by the hysteresis effect characteristic of the standard single-flip simulations.

A rare event sampling method for diffusion Monte Carlo using smart darting

We identify a set of multidimensional potential energy surfaces sufficiently complex to cause both the classical parallel tempering and the guided or unguided diffusion Monte Carlo methods to converge too inefficiently for practical applications. The mathematical model is constructed as a linear combination of decoupled Double Wells [(DDW)(n)]. We show that the set (DDW)(n) provides a serious test for new methods aimed at addressing rare event sampling in stochastic simulations. Unlike the typical numerical tests used in these cases, the thermodynamics and the quantum dynamics for (DDW)(n) can be solved deterministically. We use the potential energy set (DDW)(n) to explore and identify methods that can enhance the diffusion Monte Carlo algorithm. We demonstrate that the smart darting method succeeds at reducing quasiergodicity for n ≫ 100 using just 1 × 10(6) moves in classical simulations (DDW)(n). Finally, we prove that smart darting, when incorporated into the regular or the guided diffusion Monte Carlo algorithm, drastically improves its convergence. The new method promises to significantly extend the range of systems computationally tractable by the diffusion Monte Carlo algorithm.

Collapse transitions in a flexible homopolymer chain: application of the Wang-Landau algorithm

The thermodynamic behavior of a continuous homopolymer is described using the Wang-Landau algorithm for chain lengths up to N=561. The coil-globule and liquid-solid transitions are analyzed in detail with traditional thermodynamic and structural quantities. The behavior of the coil-globule transition is well within theoretical and computational predictions for all chain lengths, while the behavior of the liquid-solid transition is much more susceptible to finite-size effects. Certain "magic number" lengths (N=13,55,147,309,561) , whose minimal energy states offer unique icosahedral geometries, are discussed along with chains residing between these special cases. The low temperature behavior near the liquid-solid transition is rich in structural transformations for certain chain lengths, showing many similarities to the behavior of classical clusters with similar interaction potentials. General comments are made on this size dependent behavior and how it affects transition behavior in this model.

Ab Initio Molecular Dynamical Investigation of the Finite Temperature Behavior of the Tetrahedral Au19and Au20Clusters

Density functional molecular dynamics simulations have been carried out to understand the finite temperature behavior of Au19 and Au20 clusters. Au20 has been reported to be a unique molecule having tetrahedral geometry, a large HOMO-LUMO energy gap, and an atomic packing similar to that of the bulk gold (Li, J.; et al. Science 2003, 299, 864). Our results show that the geometry of Au19 is exactly identical with that of Au20 with one missing corner atom (called a vacancy). Surprisingly, our calculated heat capacities for this nearly identical pair of gold clusters exhibit dramatic differences. Au20 undergoes a clear and distinct solid-like to liquid-like transition with a sharp peak in the heat capacity curve around 770 K. On the other hand, Au19 has a broad and flat heat capacity curve with continuous melting transition. This continuous melting transition turns out to be a consequence of a process involving a series of atomic rearrangements along the surface to fill in the missing corner atom. This results in a restricted diffusive motion of atoms along the surface of Au19 between 650 to 900 K during which the shape of the ground state geometry is retained. In contrast, the tetrahedral structure of Au20 is destroyed around 800 K, and the cluster is clearly in a liquid-like state above 1000 K. Thus, this work clearly demonstrates that (i) the gold clusters exhibit size sensitive variations in the heat capacity curves and (ii) the broad and continuous melting transition in a cluster, a feature that has so far been attributed to the disorder or absence of symmetry in the system, can also be a consequence of a defect (absence of a cap atom) in the structure.

Accurate modeling of sequential decay in clusters over long time scales: Insights from phase space theory

A general theoretical framework for describing the thermally induced sequential decay in atomic clusters is presented. The scheme relies on a full treatment of individual dissociation steps based on phase space theory (PST), built into a kinetic Monte Carlo (kMC) procedure. This combined PST/kMC approach allows one to follow the evolution of several statistical properties such as the size, the angular momentum, or the temperature of the cluster over arbitrarily long time scales. Quantitative accuracy is achieved by incorporating anharmonicities of the vibrational densities of states, the rigorous conservation of angular momentum via the effective dissociation potential, and a proper calibration of the rate constants. The approach is tested and validated on selected Lennard-Jones clusters in various situations. Several approximations, including a mean-field rate equation treatment, are critically discussed; possible extensions are presented.

Structural transitions in the 309-atom magic number Lennard-Jones cluster

The thermal behavior of the 309-atom Lennard-Jones cluster, whose structure is a complete Mackay icosahedron, has been studied by parallel tempering Monte Carlo simulations. Surprisingly for a magic number cluster, the heat capacity shows a very pronounced peak before melting, which is attributed to several coincident structural transformation processes. The main transformation is somewhat akin to surface roughening and involves a cooperative condensation of vacancies and adatoms that leads to the formation of pits and islands one or two layers thick on the Mackay icosahedron. The second transition in order of importance involves a whole scale transformation of the cluster structure and leads to a diverse set of twinned structures that are assemblies of face-centered-cubic tetrahedra with six atoms along their edges, i.e., one atom more than the edges of the 20 tetrahedra that make up the 309-atom Mackay icosahedron. A surface reconstruction of the icosahedron from a Mackay to an anti-Mackay overlayer is also observed, but with a lower probability.

Equilibrium thermodynamics from basin-sampling

We present a "basin-sampling" approach for calculation of the potential energy density of states for classical statistical models. It combines a Wang-Landau-type uniform sampling of local minima and a novel approach for approximating the relative contributions from local minima in terms of the volumes of basins of attraction. We have employed basin-sampling to study phase changes in atomic clusters modeled by the Lennard-Jones potential and for ionic clusters. The approach proves to be efficient for systems involving broken ergodicity and has allowed us to calculate converged heat capacity curves for systems that could previously only be treated using the harmonic superposition approximation. Benchmarks are also provided by comparison with parallel tempering and Wang-Landau simulations, where these proved feasible.

Parallel tempering simulations of the 13-center Lennard-Jones dipole-dipole cluster (muD=0-->0.5 a.u.)

We investigate the thermodynamic behavior of the thirteen center uniform Lennard-Jones dipole-dipole cluster [(LJDD)(13)] for a wide range of dipole moment strengths. We find a relatively wide range of potential parameters where solid-solid coexistence manifests itself. Using structural characterization methods we determine the shape of the few isomers that contribute to the solid-solid coexistence region. The thermal distributions of the size of the net dipole moment are broad even at the coldest temperatures of the simulation where the (LJDD)(13) cluster is solid.

Phase changes in Lennard-Jones mixed clusters with composition ArnXe6−n (n=0,1,2)

We have carried out parallel tempering Monte Carlo calculations on the binary six-atom mixed Lennard-Jones clusters, Ar(n)Xe(6-n) (n=0,1,2). We have looked at the classical configurational heat capacity C(V)(T) as a probe of phase behavior. All three clusters show a feature in the heat capacity in the region of 15-20 K. The Ar(2)Xe(4) cluster exhibits a further peak in the heat capacity near 7 K. We have also investigated dynamical properties of the Ar(2)Xe(4) cluster as a function of temperature using molecular dynamics. We report the interbasin isomerization rate and the bond fluctuation parameter obtained from these calculations. At 7 K, the isomerization rate is on the order of 0.01 ns(-1); at 20 K, the isomerization rate is greater than 10 ns(-1). Furthermore, at 7 K, the bond fluctuation parameter is less than 3%; at 20 K, it is in the range of 10-15% (depending on the sampling time used). Using this information, together with Monte Carlo quenching data, we assign the 15-20 K feature in the heat capacity to a solid-liquid phase change and the 7-K peak to a solid-solid phase change. We believe this is the smallest Lennard-Jones cluster system yet shown to exhibit solid-solid phase change behavior.

Size-temperature phase diagram for small Lennard-Jones clusters

The size-temperature "phase diagram" for small Lennard-Jones clusters LJn (n< or =147) is constructed based on the analysis of the heat capacities Cv(T) , computed by the parallel tempering Monte Carlo method. Two types of "phase transitions" are often observed: the higher-temperature solid-liquid transition, due to melting of the cluster core, and the lower-temperature melting that occurs within the surface layer. The latter transition can only be observed for clusters with Mackay packing of the overlayer.

Parameter space minimization methods: Applications to Lennard-Jones–dipole-dipole clusters

The morphology of the uniform Lennard-Jones-dipole-dipole cluster with 13 centers (LJDD)13 is investigated over a relatively wide range of values of the dipole moment. We introduce and compare several necessary modifications of the basin-hopping algorithm for global optimization to improve its efficiency. We develop a general algorithm for T=0 Brownian dynamics in curved spaces, and a graph theoretical approach necessary for the elimination of dissociated states. We find that the (LJDD)13 cluster has icosahedral symmetry for small to moderate values of the dipole moment. As the dipole moment increases, however, its morphology shifts to an hexagonal antiprism, and eventually to a ring.

Effect of potential energy distribution on the melting of clusters

We find that the potential energy distribution of atoms in clusters can consistently explain many important phenomena related to the phase changes of clusters, such as the nonmonotonic variation of melting temperature with size, the dependence of melting, boiling, and sublimation temperatures on the interatomic potentials, the existence of a surface-melted phase, and the absence of a premelting peak in heat capacity curves. We also find a new type of premelting mechanism in double icosahedral Pd19 clusters, where one of the two internal atoms escapes to the surface at the premelting temperature.

Approach to ergodicity in monte carlo simulations

The approach to the ergodic limit in Monte Carlo simulations is studied using both analytic and numerical methods. With the help of a stochastic model, a metric is defined that enables the examination of a simulation in both the ergodic and nonergodic regimes. In the nonergodic regime, the model implies how the simulation is expected to approach ergodic behavior analytically, and the analytically inferred decay law of the metric allows the monitoring of the onset of ergodic behavior. The metric is related to previously defined measures developed for molecular dynamics simulations, and the metric enables the comparison of the relative efficiencies of different Monte Carlo schemes. Applications to Lennard-Jones 13-particle clusters are shown to match the model for Metropolis, J-walking, and parallel tempering based approaches. The relative efficiencies of these three Monte Carlo approaches are compared, and the decay law is shown to be useful in determining needed high temperature parameters in parallel tempering and J-walking studies of atomic clusters.