Molecular dynamics studies of nanoparticle impacts

This thesis concerns the dynamics of nanoparticle impacts on solid surfaces. These impacts occur, for instance, in space, where micro- and nanometeoroids hit surfaces of planets, moons, and spacecraft. On Earth, materials are bombarded with nanoparticles in cluster ion beam devices, in order to clean or smooth their surfaces, or to analyse their elemental composition. In both cases, the result depends on the combined effects of countless single impacts. However, the dynamics of single impacts must be understood before the overall effects of nanoparticle radiation can be modelled. In addition to applications, nanoparticle impacts are also important to basic research in the nanoscience field, because the impacts provide an excellent case to test the applicability of atomic-level interaction models to very dynamic conditions. In this thesis, the stopping of nanoparticles in matter is explored using classical molecular dynamics computer simulations. The materials investigated are gold, silicon, and silica. Impacts on silicon through a native oxide layer and formation of complex craters are also simulated. Nanoparticles up to a diameter of 20 nm (315000 atoms) were used as projectiles. The molecular dynamics method and interatomic potentials for silicon and gold are examined in this thesis. It is shown that the displacement cascade expansionmechanism and crater crown formation are very sensitive to the choice of atomic interaction model. However, the best of the current interatomic models can be utilized in nanoparticle impact simulation, if caution is exercised. The stopping of monatomic ions in matter is understood very well nowadays. However, interactions become very complex when several atoms impact on a surface simultaneously and within a short distance, as happens in a nanoparticle impact. A high energy density is deposited in a relatively small volume, which induces ejection of material and formation of a crater. Very high yields of excavated material are observed experimentally. In addition, the yields scale nonlinearly with the cluster size and impact energy at small cluster sizes, whereas in macroscopic hypervelocity impacts, the scaling 2 is linear. The aim of this thesis is to explore the atomistic mechanisms behind the nonlinear scaling at small cluster sizes. It is shown here that the nonlinear scaling of ejected material yield disappears at large impactor sizes because the stopping mechanism of nanoparticles gradually changes to the same mechanism as in macroscopic hypervelocity impacts. The high yields at small impactor size are due to the early escape of energetic atoms from the hot region. In addition, the sputtering yield is shown to depend very much on the spatial initial energy and momentum distributions that the nanoparticle induces in the material in the first phase of the impact. At the later phases, the ejection of material occurs by several mechanisms. The most important mechanism at high energies or at large cluster sizes is atomic cluster ejection from the transient liquid crown that surrounds the crater. The cluster impact dynamics detected in the simulations are in agreement with several recent experimental results. In addition, it is shown that relatively weak impacts can induce modifications on the surface of an amorphous target over a larger area than was previously expected. This is a probable explanation for the formation of the complex crater shapes observed on these surfaces with atomic force microscopy. Clusters that consist of hundreds of thousands of atoms induce long-range modifications in crystalline gold.

Vibrational dynamics and boson peak in a supercooled polydisperse liquid

Vibrational density of states (VDOS) in a supercooled polydisperse liquid is computed by diagonalizing the Hessian matrix evaluated at the potential energy minima for systems with different values of polydispersity. An increase in polydispersity leads to an increase in the relative population of localized high-frequency modes. At low frequencies, the density of states shows an excess compared to the Debye squared-frequency law, which has been identified with the boson peak. The height of the boson peak increases with polydispersity and shows a rather narrow sensitivity to changes in temperature. While the modes comprising the boson peak appear to be largely delocalized, there is a sharp drop in the participation ratio of the modes that exist just below the boson peak indicative of the quasilocalized nature of the low-frequency vibrations. Study of the difference spectrum at two different polydispersity reveals that the increase in the height of boson peak is due to a population shift from modes with frequencies above the maximum in the VDOS to that below the maximum, indicating an increase in the fraction of the unstable modes in the system. The latter is further supported by the facilitation of the observed dynamics by polydispersity. Since the strength of the liquid increases with polydispersity, the present result provides an evidence that the intensity of boson peak correlates positively with the strength of the liquid, as observed earlier in many experimental systems.

The dynamics of solvation of an electron in the image potential state by a layer of polar adsorbates

Recently, ultrafast two-photon photoemission has been used to study electron solvation at a two-dimensional metal/polar adsorbate interfaces [A. Miller , Science 297, 1163 (2002)]. The electron is bound to the surface by the image interaction. Earlier we have suggested a theoretical description of the states of the electron interacting with a two-dimensional layer of the polar adsorbate [K. L. Sebastian , J. Chem. Phys. 119, 10350 (2003)]. In this paper we have analyzed the dynamics of electron solvation, assuming a trial wave function for the electron and the solvent polarization and then using the Dirac-Frenkel variational method to determine it. The electron is initially photoexcited to a delocalized state, which has a finite but large size, and causes the polar molecules to reorient. This reorientation acts back on the electron and causes its wave function to shrink, which will cause further reorientation of the polar molecules, and the process continues until the electron gets self-trapped. For reasonable values for the parameters, we are able to obtain fair agreement with the experimental observations. (c) 2005 American Institute of Physics.

Lattice dynamics of hexagonal ice

The lattice dynamics of hexagonal ice is worked out with the force constants deduced from the experimental elastic constants.

Levitation effect in zeolites: Quasielastic neutron scattering and molecular dynamics study of pentane isomers in zeolite NaY

We report the quasielastic neutron scattering (QENS) and molecular dynamics (MD) investigations into diffusion of pentane isomers in zeolite NaY. The molecular cross section perpendicular to the long molecular axis varies for the three isomers while the mass and the isomer-zeolite interaction remains essentially unchanged. Both QENS and MD results show that the branched isomers neopentane and isopentane have higher self-diffusivities as compared with n-pentane at 300 K in NaY zeolite. This result provides direct experimental evidence for the existence of nonmonotonic, anomalous dependence of self-diffusivity on molecular diameter known as the levitation effect. The energetic barrier at the bottleneck derived from MD simulations exists for n-pentane which lies in the linear regime while no such barrier is seen for neopentane which is located clearly in the anomalous regime.Activation energy is in the order E-a(n-pentane)>E-a(isopentane)>E-a(neopentane) consistent with the predictions of the levitation effect. In the liquid phase, it is seen thatD(n pentane)>D(isopentane)>D(neopentane) and E-a(n-pentane)< E-a(isopentane)< E-a(neopentane). Intermediate scattering function for small wavenumbers obtained from MD follows a single exponential decay for neopentane and isopentane. For n-pentane, a single exponential fit provides a poor fit especially at short times. Cage residence time is largest for n-pentane and lowest for neopentane. For neopentane, the width of the self-part of the dynamic structure factor shows a near monotonic decrease with wavenumber. For n-pentane a minimum is seen near k=0.5 A degrees(-1) suggesting a slowing down of motion around the 12-ring window, the bottleneck for diffusion. Finally, the result that the branched isomer has a higher diffusivity as compared with the linear analog is at variation from what is normally seen.

Dislocation dynamics in alpha titanium Dynamique des dislocations dans le titane alphaVersetzungsdynamik in alpha-Titan

The activation area and activation enthalpy are determined as a function of stress and temperature for alpha titanium. The results indicated that plastic flow below about 700°K occurs by a single thermally activated mechanism. Activation area determined by differential-stress creep tests falls in the range 80−8b2 and does not systematically depend on the impurity content. The total activation enthalpy derived from the temperature and strain-rate dependence of flow stress is 1.15 eV. The experimental data support a lattice hardening mechanism as controlling the low-temperature deformation in alpha titanium.

On the generalized continuum model of dipolar solvation dynamics

The continuum model of dipolar solvation dynamics is reviewed. The effects of non-spherical molecular shapes, of non-Debye dielectric relaxation of the polar solvent and of dielectric inhomogeneity of the solvent around the solute dipole are investigated. Several new theoretical results are presented. It is found that our generalized continuum model, which takes into account the dielectric inhomogeneity of the surrounding solvent, provides a description of solvation dynamics consistent with recent experimental results.

Inertial effects in solvation dynamics

The dynamics of solvation of newly created charged species in dense dipolar liquids can proceed at a high speed with time constants often in the subpicosecond domain. The motion of the solvent molecules can be in the inertial limit at such short times. In this paper we present a microscopic study of the effects of inertial motion of solvent molecules on the solvation dynamics of a newly created ion in a model dipolar liquid. Interesting dynamical behavior emerges when the relative contribution of the translational modes in the wave-vector-dependent longitudinal relaxation time is significant. Especially, the theory predicts that the time correlation function of the solvation energy can become oscillatory in some limiting situations. In general, the dynamics becomes faster in the presence of the inertial contribution. We discuss the experimental situations where the inertial effects can be noticeable.

Temperature and concentration dependence of adsorption properties of methane in NaY: a molecular dynamics study

Molecular dynamics calculations on methane sorbed in NaY (Si/Al = 3.0) employing realistic methane-methane and methane-zeolite intermolecular potential functions at different temperatures (50, 150, 220, and 300 K) and concentrations (2, 4, 6, and 8 molecules/cage) are reported. The thermodynamic results are in agreement with the available experimental data. Guest-guest and guest-host radial distribution functions (rdfs), energy distribution functions, distribution of cage occupancy, center-of-cage-center-of-mass (coc-com) rdfs, velocity autocorrelation functions for com and angular motion and the Fourier transformed power spectra, and diffusion coefficients are presented as a function of temperature and concentration. At 50 K, methane is localized near the adsorption site. Site-site migration and essentially free rotational motion are observed at 150 K. Molecules preferentially occupy the region near the inner surface of the alpha-cage. The vibrational frequencies for the com of methane shift toward higher values with decreasing temperature and increasing adsorbate concentration. The observed frequencies for com motion are 36, 53, and 85 cm-1 and for rotational motion at 50 K, 95 and 150 cm-1 in agreement with neutron scattering data. The diffusion coefficients show a type I behavior as a function of loading in agreement with NMR measurements. Cage-to-cage diffusion is found to be always mediated by the surface.

Communication pathways from the anti-codon region to the aminoacylation site in Methionyl tRNA Synthetase: Molecular dynamics simulations and the structure network analysis

The enzymes of the family of tRNA synthetases perform their functions with high precision by synchronously recognizing the anticodon region and the aminoacylation region, which are separated by ?70 in space. This precision in function is brought about by establishing good communication paths between the two regions. We have modeled the structure of the complex consisting of Escherichia coli methionyl-tRNA synthetase (MetRS), tRNA, and the activated methionine. Molecular dynamics simulations have been performed on the modeled structure to obtain the equilibrated structure of the complex and the cross-correlations between the residues in MetRS have been evaluated. Furthermore, the network analysis on these simulated structures has been carried out to elucidate the paths of communication between the activation site and the anticodon recognition site. This study has provided the detailed paths of communication, which are consistent with experimental results. Similar studies also have been carried out on the complexes (MetRS + activated methonine) and (MetRS + tRNA) along with ligand-free native enzyme. A comparison of the paths derived from the four simulations clearly has shown that the communication path is strongly correlated and unique to the enzyme complex, which is bound to both the tRNA and the activated methionine. The details of the method of our investigation and the biological implications of the results are presented in this article. The method developed here also could be used to investigate any protein system where the function takes place through long-distance communication.

A molecular dynamics study of a PbO.SiO2 glass and melt

The structures of a PbO.SiO2 glass and melt have been studied using molecular dynamics simulation employing Born-Mayer-Huggins pair potentials. Various pair distribution functions are presented and discussed. Pb-Pb correlations persist in the melt, in agreement with experimental observations. The calculated and experimental radial distribution functions are compared.

Dynamics of activationless reactions in solution

Recent picosecond and subpicosecond laser spectroscopy experiments have revealed several chemically and biologically important reactions in solution in which the reaction potential surface does not present a barrier to the motion along the reaction coordinate.The dynamics of these reactions display diverse and interesting behavior. They include the dependence of relaxation rate on the solvent viscosity, the solvent polarity, the temperature, and the wavelength of the exciting light. In this article we review the recent developments in the theoretical description of activationless processes in solution and compare them with the available experimental results

Sequence Independent Recombination Triple Helices: A Molecular Dynamics Study

Recent experimental studies have shown that the Rec-A mediated homologous recombination reaction involves a triple helical intermediate, in which the third strand base forms hydrogen bonds with both the bases in the major groove of the Watson-Crick duplex. Such 'mixed' hydrogen bonds allow formation of sequence independent triplexes. DNA triple helices involving 'mixed' hydrogen bonds have been studied, using model building, molecular mechanics (MM) and molecular dynamics (MD). Models were built for a tripler comprising all four possible triplets viz., G.C*C, C.G*G, A.T*T and T.A*A. To check the stability of all the 'mixed' hydrogen bonds in such triplexes and the conformational preferences of such tripler structures, MD studies were carried out starting from two structures with 30 degrees and 36 degrees twist between the basepairs. It was observed that though the two triplexes converged towards a similar structure, the various hydrogen bonds between the WC duplex and the third strand showed differential stabilities. An MD simulation with restrained hydrogen bonds showed that the resulting structure was stable and remained close to the starting structure. These studies help us in defining stable hydrogen bond geometries involving the third strand and the WC duplex. It was observed that in the C.G*G triplets the N7 atom of the second strand is always involved in hydrogen bonding. In the G.C*C triplets, either N3 or O2 in the third strand cytosine can interchangeably act as a hydrogen bond acceptor.

Molecular theory of ion solvation dynamics in water, acetonitrile and methanol: A unified microscopic description of collective dynamics in dipolar liquids

A recently developed microscopic theory of solvation dynamics in real dipolar liquids is used to calculate, for the first time, the solvation time correlation function in liquid acetonitrile, water and methanol. The calculated results are in excellent agreement with known experimental and computer simulation studies.

Solvation dynamics in liquid water. A novel interplay between librational and diffusive modes

A microscopic calculation of the solvation dynamics of an ion in liquid water is presented. The calculated solvation time correlation function shows an ultrafast Gaussian decay which carries about 70%–90% of the strength followed by a biexponential decay with time constants equal to 250 fs and 1 ps. These results are in excellent agreement with the computer simulations of Maroncelli and Fleming and also with the experimental findings of Barbara and Jarzeba. In addition, we find that both the rotational librations and the intermolecular translational vibrational modes of water contribute significantly to the initial Gaussian decay.