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.

Purification, characterization and molecular cloning of a monocot mannose-binding lectin from Remusatia vivipara with nematicidal activity

A mannose-binding lectin (RVL) was purified from the tubers of Remusatia vivipara, a monocot plant by single-step affinity chromatography on asialofetuin-Sepharose 4B. RVL agglutinated only rabbit erythrocytes and was inhibited by mucin, asialomucin, asialofetuin and thyroglobulin. Lectin activity was stable up to 80A degrees C and under wide range of pH (2.0-9.3). SDS-PAGE and gel filtration results showed the lectin is a homotetramer of Mr 49.5 kDa, but MALDI analysis showed two distinct peaks corresponding to subunit mass of 12 kDa and 12.7 kDa. Also the N-terminal sequencing gave two different sequences indicating presence of two polypeptide chains. Cloning of RVL gene indicated posttranslational cleavage of RVL precursor into two mature polypeptides of 116 and 117 amino-acid residues. Dynamic light scattering (DLS) and gel filtration studies together confirmed the homogeneity of the purified lectin and supported RVL as a dimer with Mr 49.5 kDa derived from single polypeptide precursor of 233 amino acids. Purified RVL exerts potent nematicidal activity on Meloidogyne incognita, a root knot nematode. Fluorescent confocal microscopic studies demonstrated the binding of RVL to specific regions of the alimentary-tract and exhibited a potent toxic effect on M. incognita. RVL-mucin complex failed to interact with the gut confirming the receptor mediated lectin interaction. Very high mortality (88%) rate was observed at lectin concentration as low as 30 A mu g/ml, suggesting its potential application in the development of nematode resistant transgenic-crops.

Ion channel formation by zervamicin-IIB. A molecular modelling study.

Zervamicin-IIB (Zrv-IIB) is a 16 residue peptaibol which forms voltage-activated, multiple conductance level channels in planar lipid bilayers. A molecular model of Zrv-IIB channels is presented. The structure of monomeric Zrv-IIB is based upon the crystal structure of Zervamicin-Leu. The helical backbone is kinked by a hydroxyproline residue at position 10. Zrv-IIB channels are modelled as helix bundles of from 4 to 8 parallel helices surrounding a central pore. The monomers are packed with their C-terminal helical segments in close contact, and the bundles are stabilized by hydrogen bonds between glutamine 11 and hydroxyproline 10 of adjacent helices. Interaction energy profiles for movement of three different probes species (K+, Cl- and water) through the central pore are analyzed. The conformations of: (a) the sidechain of glutamine 3; (b) the hydroxyl group of hydroxyproline 10; and (c) the C-terminal hydroxyl group are "optimized" in order to maximize favourable interactions between the channel and the probes, resulting in favourable interaction energy profiles for all three. This suggests that conformational flexibility of polar sidechains enables the channel lining to mimic an aqueous environment.

Charge density based classification of intermolecular interactions in molecular crystals

Evaluation of intermolecular interactions in terms of both experimental and theoretical charge density analyses has produced a unified picture with which to classify strong and weak hydrogen bonds, along with van der Waals interactions, into three regions.

Molecular markers for discriminating Streptococcus pyogenes and S.dysgalactiae subspecies equisimilis

Given the increasing aetiological importance of Streptococcus dysgalactiae subspecies equisimilis in diseases which are primarily attributed to S. pyogenes, molecular markers are essential to distinguish these species and delineate their epidemiology more precisely. Many clinical microbiology laboratories rely on agglutination reactivity and biochemical tests to distinguish them. These methods have limitations which are particularly exacerbated when isolates with mixed properties are encountered. In order to provide additional distinguishing parameters that could be used to unequivocally discriminate these two common pathogens, we assess here three molecular targets: the speB gene, intergenic region upstream of the scpG gene (IRSG) and virPCR. Of these, the former two respectively gave positive and negative results for S. pyogenes, and negative and positive results for S. dysgalactiae subsp. equisimilis. Thus,a concerted use of these nucleic acid-based methods is particularly helpful in epidemiological surveillance to accurately assess the relative contribution of these species to streptococcal infections and diseases.

N-H center dot center dot center dot O and C-H center dot center dotcenter dot O interaction mimicry in the 1:1 molecular complexes of 5,5-diethylbarbituric acid with urea and acetamide

The analogy between N-H center dot center dot center dot O and C-H center dot center dot center dot O intermolecular interactions is studied with variable temperature (180-100 K) single crystal X-ray diffraction analysis.5,5-Diethylbarbituric acid (barbital) forms isostructural molecular complexes (co-crystals) with urea (1) and acetamide (2) that respectively contain these analogous interactions.The behaviour of these two interactions as a function of temperature is very similar. This indicates that the C-H center dot center dot center dot O bond in barbital acetamide plays a similar chemical and structural role as does the N-H center dot center dot center dot O bond in barbital urea. The close relationship between these interactions and their comparable nature is further adduced from the formation of a ternary solid solution (3) of barbital, urea and acetamide. The fact that the C-H center dot center dot center dot O interaction in barbital acetamide is weaker than the N-H center dot center dot center dot O interaction in barbital urea is shown by the fact that acetamide is under expressed and urea is over expressed with respect to the quantities of these substances present in solution prior to crystallization of these ternary crystals.

Specificity for the Exchange of Phospholipids Through Polymyxin B Mediated Intermembrane Molecular Contacts

Structural specificity for the direct vesicle−vesicle exchange of phospholipids through stable molecular contacts formed by the antibiotic polymyxin B (PxB) is characterized by kinetic and spectroscopic methods. As shown elsewhere [Cajal, Y., Rogers, J., Berg, O. G., & Jain, M. K. (1996) Biochemistry 35, 299−308], intermembrane molecular contacts between anionic vesicles are formed by a small number of PxB molecules, which suggests that a stoichiometric complex may be responsible for the exchange of phospholipids. Larger clusters containing several vesicles are formed where each vesicle can make multiple contacts if sterically allowed. In this paper we show that the overall process can be dissected into three functional steps: binding of PxB to vesicles, formation of stable vesicle−vesicle contacts, and exchange of phospholipids. Polycationic PxB binds to anionic vesicles. Formation of molecular contacts and exchange of monoanionic phospholipids through PxB contacts does not depend on the chain length of the phospholipid. Only monoanionic phospholipids (with methanol, serine, glycol, butanol, or phosphatidylglycerol as the second phosphodiester substituent in the head group) exchange through these contacts, whereas dianionic phosphatidic acid does not. Selectivity for the exchange was also determined with covesicles of phosphatidylmethanol and other phospholipids. PxB does not bind to vesicles of zwitterionic phosphatidylcholine, and its exchange in covesicles is not mediated by PxB. Vesicles of dianionic phospholipids, like phosphatidic acid, bind PxB; however, this phospholipid does not exchange. The structural features of the contacts are characterized by the spectroscopic and chemical properties of PxB at the interface. PxB in intermembrane contacts is readily accessible from the aqueous phase to quenchers and reagents that modify amino groups. Results show that PxB at the interface can exist in two forms depending on the lipid/PxB ratio. Additional studies show that stable PxB-mediated vesicle−vesicle contacts may be structurally and functionally distinct from “stalks”, the putative transient intermediate for membrane fusion. The phenomenon of selective exchange of phospholipids through peptide-mediated contacts could serve as a prototype for intermembrane targeting and sorting of phospholipids during their biosynthesis and trafficking in different compartments of a cell. The protocols and results described here also extend the syllogistic foundations of interfacial equilibria and catalysis.

Study of molecular reorientation and quantum rotational tunneling in tetramethylammonium selenate by H-1 NMR

H-1 NMR spin-lattice relaxation time measurements have been carried out in [(CH3)(4)N](2)SeO4 in the temperature range 389-6.6K to understand the possible phase transitions, internal motions and quantum rotational tunneling. A broad T, minimum observed around 280K is attributed to the simultaneous motions of CH3 and (CH3)(4)N groups. Magnetization recovery is found to be stretched exponential below 72 K with varying stretched exponent. Low-temperature T-1 behavior is interpreted in terms of methyl groups undergoing quantum rotational tunneling. (c) 2007 Elsevier Inc. All rights reserved.

Orientational melting and reorientational motion in a cubane molecular crystal: A molecular simulation study

Detailed molecular simulations are carried out to investigate the effect of temperature on orientational order in cubane molecular crystal. We report a transition from an orientationally ordered to an orientationally disordered plastic crystalline phase in the temperature range 425-450 K. This is similar to the experimentally reported transition at 395 K. The nature of this transition is first order and is associated with a 4.8% increase in unit Cell volume that is comparable to the experimentally reported unit cell volume change of 5.4% (Phys. Rev. Lett. 1997, 78, 4938). An orientational order parameter, eta(T), has been defined in terms of average angle of libration of a molecular 3-fold axis and the orientational melting has been characterized by using eta(T). The orientational melting is associated with an anomaly in specific heat at constant pressure (C-p) and compressibility (kappa). The enthalpy of transition and entropy of transition associated with this orientational melting are 20.8 J mol(-1) and 0.046 J mol(-1) K-1, respectively. The structure of crystalline as well as plastic crystalline phases is characterized by using various radial distribution functions and orientational distribution functions. The coefficient of thermal expansion of the plastic crystalline phase is more than twice that of the crystalline phase.

Large molecular motions are tolerated in crystals of diamine double salt of trans-chlorocinnamic acids with trans-1,2-diaminocyclohexane

Contrary to the general assumption that photoreactions in crystals may not proceed with large molecular motions, a pedal-like motion prompted by electronic excitation is believed to be involved during the β-dimer formation from the crystals of the diamine double salt of trans-2,4-dichlorocinnamic acid and trans-1,2-diaminocyclohexane.

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.