Cluster dynamics, calculations

The scaling of the relaxation modulus G(t) with time (Eq. 1-1) at the LST was first detected experimentally [5-7]. Subsequently, dynamic scaling based on percolation theory used the relation between diffusion coefficient and longest relaxation time of a single cluster to calculate a relaxation time spectrum for the sum of all clusters [39], This resulted in the same scaling relation for G(t) with an exponent n following Eq. 1-14. 

Interestingly, in the experiments devoted solely to computational chemistry, molecular dynamics calculations had the highest representation (96-98). The method was used in simulations of simple liquids, (96), in simulations of chemical reactions (97), and in studies of molecular clusters (98). One experiment was devoted to the use of Monte Carlo methods to distinguish between first and second-order kinetic rate laws (99). One experiment used DFT theory to study two isomerization reactions (100). 

Both analytic theory and computer simulations are included, and we note that the latter play an especially important role in understanding cluster reactions. Simulations not only provide quantitative results, but they provide insight into the dominant causes of observed behavior, and they can provide likelihood estimates for assessing qualitatively distinct mechanisms that can be used to explain the same experimental data. Simulations can also lead to a greater understanding of dynamical processes occurring in clusters by calculating details which cannot be observed experimentally. 

Perera and Amar (1989) found more detailed support for the structural control of caging in classical dynamics calculations on a model of Br2 in large clusters of Ar and C02. The dissociation channel was found to become closed, as a function of cluster size, between 11 and 12 C02 molecules in the BrJ(C02)M clusters, correlating with the appearance of double-capped minimum energy structures. This correlation was found in the Br2 Ar clusters as well. Collisions between a vibrating diatomic molecule in a cluster and the solvent particles may cause V-T energy transfer and rapid evaporation of the cluster. 

A wide variety of dynamical approximations have been applied to cluster dynamics and kinetics. Most calculations to date are based on simplified potentials and classical mechanics or statistical methods. In the near future, we can expect to see more work with detailed potential energy surfaces (both analytic and implicitly defined by electronic structure calculations) and progress in sorting out quantum effects and treating them more accurately. 

In order to interpret the results of our experiments, optimal-control calculations were performed where a GA controlled 40 independent degrees of freedom in the laser pulses that were used in a molecular dynamics simulation of the laser-cluster interactions for Xejv clusters with sizes ranging from 108 to 5056 atoms/cluster. These calculations, which are reported in detail elsewhere [67], showed optimization of the laser-cluster interactions by a sequence of as many as three laser pulses. Detailed inspection of the simulations revealed that the first pulse in this sequence initiates the cluster ionization and starts the expansion of the cluster, while the second and third pulse optimize two mechanisms that are directly related to the behaviour of the electrons in the cluster. We consistently observe that the second pulse in the three-pulse sequence arrives a time delay where the conditions for enhanced ionization are met. In other words, the second pulse arrives at a time where the ionization of atoms is assisted by the proximity of surrounding ions. The third peak is consistently observed at a delay where the collective oscillation of the quasi-free electrons in the cluster is 7t/2 out of phase with respect to the driving laser field. For a driven and damped oscillator this phase-delay represents an optimum for the energy transfer from the driving force to the oscillator. 

Structure of Br may not be the same as that of the bulk. Some of the molecular dynamics calculations predict that halide anions in water tend to float on the surface of clusters consisting of water molecules rather than within water. This effect may cause a dissimilar solvation structure to that of the bulk. In addition, if the anion is segregated at the surface by surfactants such as large alkylammonium cations, the anion density at the surface should be high and its environment differ from the bulk. This is a preliminary report of the first experimental study of the solution surface by the EXAFS technique. This technique provides us information on the gas/liquid interface, the structure of Langmuir films, and the effect of the interface on chemical reactions. 

TD-DFT) they calculated the transition energies and dipole moments for NMA both in vacuum and in an aqueous solution. Moreover, in the treatment of the solvent they compared two different approaches, i.e., a polarizable-continuum method (COSMO) and a supermolecule approach. For the latter, the authors performed molecular-dynamics calculations using a force-field model and, subsequently, extracted a cluster containing the solute and 3 water molecules that form hydrogen bonds to the solute. Averages over 90 such configurations were ultimately determined. 

Two methods are in common use for simulating molecular liquids the Monte Carlo method (MC) and molecular dynamics calculations (MD). Both depend on the availability of reasonably accurate potential energy surfaces and both are based on statistical classical mechanics, taking no account of quantum effects. In the past 10-15 years quantum Monte Carlo methods (QMC) have been developed that allow intramolecular degrees of freedom to be studied, but because of the computational complexity of this approach results have only been reported for water clusters. 

K. Uehara, M. Ishitobi, T. Oda, and Y. Hiwatari, First-principle molecular dynamics calculation of selenium clusters. Mol. Simul., 18 (1997), 385-394. 

However, these molecular dynamics calculations suffer some limitations the empirical nature of the potential (especially for the metal-support interaction) and the arbitrary separation between the metal-metal and metal-support interactions (the metal-metal potential is probably perturbed near the interface). Indeed, according to the type of potential used, very different results are obtained. In the case of Pd/MgO, a mean dilatation [91] or contraction [92] is observed. For finite-temperature molecular dynamics, the calculations are limited to very short times and it is not sure that the equilibrium shape is reached. As we have seen in the last section the cluster shape can be blocked for a long time on facetted metastable shapes. 

Nauchitel and Pertsin have studied the melting properties of 13-, 19-, and 55-particle Lennard-Jones clusters.Questioning the validity of results obtained from free-volume simulations of such systems, they have used hard-sphere boundaries to constrain their clusters to finite volumes. The results of Nauchitel and Pertsin are most interesting for the 55-particle cluster. For certain ranges of temperature and mean density, structural evidence for surface melting was obtained projections of the cluster s coordinates, and radial density distribution functions, like those given in Fig. 17, characterize the cluster as a 13-particle icosahedral core surrounded by a fluidlike shell. However, dynamic calculations like those described for other clusters in the previous section have yet to be obtained to determine how fluidlike these outer atoms really are. 

As an alternative to the normal-mode method, Monte Carlo and molecular dynamics calculations have been performed on small clusters. Monte Carlo and molecular dynamics methods have the virtue of being exact, within calculable error bars, subject to the constraint of the approximate intermo-lecular interactions that are used. Prior to about six years ago both methods were restricted to systems projjerly described by classical mechanics. This restriction implied that systems for which tunneling or low-temjjerature vibrations were important at best could be treated approximately. 

In a much more recent study, Sieck et al.91 used molecular-dynamics calculations in studying some few, selected Srv clusters, i.e., for N = 25, 29, 35, 71, and 239, in connection with a parameterized density-functional method. The fact that they found several isomers for clusters of these sizes should not surprise, but is rather a confirmation of the finding for Lennard-Jones clusters... 

Although there exist accurate ab-initio calculations on gold cluster, these calculations are quite expensive and time consuming. On the other hand molecular-dynamics (MD) simulations provide an easy tool to obtain the structural stability of clusters. One of the authors (S.E.) has performed molecular dynamics simulations for gold clusters [14] using empirical potentials based on purely experimental data. 

Molecular dynamics calculations have been made on the motions of the Na+ ion in mordenite zeolites.9 Selenium clusters doped with Na (i.e. Na2Sen) show a Raman band in the range 165-225 cm 1 due to Na-Se motions.10 Raman microspectra have been reported for caesium oxides, e.g. an alg mode of Cs20 was seen at 103 cm 1.11... 

Narita, I. and Oku, T. (2002). Molecular dynamics calculation ofH2 gas storage in C50 and BjgNjg clusters. Diam. Relat. Mater., 11, 945-8. 

---