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.

Thermal Characteristics of SrO-2B(2)O(3) Glasses

Transparent SrO-2B(2)O(3) (SBO) glasses were fabricated via the conventional melt-quenching technique. X-ray diffraction (XRD) and differential thermal analysis (DTA) studies carried out on the as-quenched glasses confirmed their amorphous and glassy nature, respectively. The thermal parameters were evaluated for the as-quenched glass-plates using non-isothermal DTA experiments. The average values of the activation energies for the glass transition and crystallization of these glasses were 800 +/- 10 kJ/mol and 298 +/- 10 kJ/mol respectively. The values of the kinetic parameters that were obtained by different non-isothermal techniques were in close agreement.

Graphene Analogues of BN: Novel Synthesis and Properties

Enthused by the fascinating properties of graphene, we have prepared graphene analogues of BN by a chemical method with a control on the number of layers. The method involves the reaction of boric acid with urea, wherein the relative proportions of the two have been varied over a wide range. Synthesis with a high proportion of urea yields a product with a majority of 1-4 layers. The surface area of BN increases progressively with the decreasing number of layers, and the high surface area BN exhibits high CO, adsorption, but negligible H, adsorption. Few-layer BN has been solubilized by interaction with Lewis bases. We have used first-principles simulations to determine structure, phonon dispersion, and elastic properties of BN with planar honeycomb lattice-based n-layer forms. We find that the mechanical stability of BN with respect to out-of-plane deformation is quite different from that of graphene, as evident in the dispersion of their flexural modes. BN is softer than graphene and exhibits signatures of long-range ionic interactions in its optical phonons. Finally, structures with different stacking sequences of BN have comparable energies, suggesting relative abundance of slip faults, stacking faults, and structural inhomogeneities in multilayer BN.

Shock acceleration in the solar corona

In this thesis acceleration of energetic particles at collisionless shock waves in space plasmas is studied using numerical simulations, with an emphasis on physical conditions applicable to the solar corona. The thesis consists of four research articles and an introductory part that summarises the main findings reached in the articles and discusses them with respect to theory of diffusive shock acceleration and observations. This thesis gives a brief review of observational properties of solar energetic particles and discusses a few open questions that are currently under active research. For example, in a few large gradual solar energetic particle events the heavy ion abundance ratios and average charge states show characteristics at high energies that are typically associated with flare-accelerated particles, i.e. impulsive events. The role of flare-accelerated particles in these and other gradual events has been discussed a lot in the scientific community, and it has been questioned if and how the observed features can be explained in terms of diffusive shock acceleration at shock waves driven by coronal mass ejections. The most extreme solar energetic particle events are the so-called ground level enhancements where particle receive so high energies that they can penetrate all the way through Earth's atmosphere and increase radiation levels at the surface. It is not known what conditions are required for acceleration into GeV/nuc energies, and the presence of both very fast coronal mass ejections and X-class solar flares makes it difficult to determine what is the role of these two accelerators in ground level enhancements. The theory of diffusive shock acceleration is reviewed and its predictions discussed with respect to the observed particle characteristics. We discuss how shock waves can be modeled and describe in detail the numerical model developed by the author. The main part of this thesis consists of the four scientific articles that are based on results of the numerical shock acceleration model developed by the author. The novel feature of this model is that it can handle complex magnetic geometries which are found, for example, near active regions in the solar corona. We show that, according to our simulations, diffusive shock acceleration can explain the observed variations in abundance ratios and average charge states, provided that suitable seed particles and magnetic geometry are available for the acceleration process in the solar corona. We also derive an injection threshold for diffusive shock acceleration that agrees with our simulation results very well, and which is valid under weakly turbulent conditions. Finally, we show that diffusive shock acceleration can produce GeV/nuc energies under suitable coronal conditions, which include the presence of energetic seed particles, a favourable magnetic geometry, and an enhanced level of ambient turbulence.

Computer simulation of H and He effects in fusion reactor materials

Fusion power is an appealing source of clean and abundant energy. The radiation resistance of reactor materials is one of the greatest obstacles on the path towards commercial fusion power. These materials are subject to a harsh radiation environment, and cannot fail mechanically or contaminate the fusion plasma. Moreover, for a power plant to be economically viable, the reactor materials must withstand long operation times, with little maintenance. The fusion reactor materials will contain hydrogen and helium, due to deposition from the plasma and nuclear reactions because of energetic neutron irradiation. The first wall divertor materials, carbon and tungsten in existing and planned test reactors, will be subject to intense bombardment of low energy deuterium and helium, which erodes and modifies the surface. All reactor materials, including the structural steel, will suffer irradiation of high energy neutrons, causing displacement cascade damage. Molecular dynamics simulation is a valuable tool for studying irradiation phenomena, such as surface bombardment and the onset of primary damage due to displacement cascades. The governing mechanisms are on the atomic level, and hence not easily studied experimentally. In order to model materials, interatomic potentials are needed to describe the interaction between the atoms. In this thesis, new interatomic potentials were developed for the tungsten-carbon-hydrogen system and for iron-helium and chromium-helium. Thus, the study of previously inaccessible systems was made possible, in particular the effect of H and He on radiation damage. The potentials were based on experimental and ab initio data from the literature, as well as density-functional theory calculations performed in this work. As a model for ferritic steel, iron-chromium with 10% Cr was studied. The difference between Fe and FeCr was shown to be negligible for threshold displacement energies. The properties of small He and He-vacancy clusters in Fe and FeCr were also investigated. The clusters were found to be more mobile and dissociate more rapidly than previously assumed, and the effect of Cr was small. The primary damage formed by displacement cascades was found to be heavily influenced by the presence of He, both in FeCr and W. Many important issues with fusion reactor materials remain poorly understood, and will require a huge effort by the international community. The development of potential models for new materials and the simulations performed in this thesis reveal many interesting features, but also serve as a platform for further studies.

Dynamo onset as a first-order transition: Lessons from a shell model for magnetohydrodynamics

We carry out systematic and high-resolution studies of dynamo action in a shell model for magnetohydro-dynamic (MHD) turbulence over wide ranges of the magnetic Prandtl number Pr-M and the magnetic Reynolds number Re-M. Our study suggests that it is natural to think of dynamo onset as a nonequilibrium first-order phase transition between two different turbulent, but statistically steady, states. The ratio of the magnetic and kinetic energies is a convenient order parameter for this transition. By using this order parameter, we obtain the stability diagram (or nonequilibrium phase diagram) for dynamo formation in our MHD shell model in the (Pr-M(-1), Re-M) plane. The dynamo boundary, which separates dynamo and no-dynamo regions, appears to have a fractal character. We obtain a hysteretic behavior of the order parameter across this boundary and suggestions of nucleation-type phenomena.

Experimental and Theoretical Charge Density Analysis of Polymorphic Structures: The Case of Coumarin 314 Dye

Experimental charge density distributions in two known conformational polymorphs (orange and yellow) of coumarin 314 dye are analyzed based on multipole modeling of X-ray diffraction data collected at 100 K. The experimental results are compared with the charge densities derived from multipole modeling of theoretical structure factors obtained from periodic quantum calculation with density functional theory (DFT) method and B3LYP/6-31G(d,p) level of theory. The presence of disorder at the carbonyl oxygen atom of ethoxycarbonyl group in the yellow form, which was not identified earlier, is addressed here. The investigationof intermolecular interactions, based on Hirshfeld surface analysis and topological properties via quantum theory of atoms in molecule and total electrostatic interaction energies, revealed significant differences between the polymorphs. The differences of electrostatic nature in these two polymorphic forms were unveiled via construction of three-dimensional deformation electrostatic potential maps plotted over the molecular surfaces. The lattice energies evaluated from ab initio calculations on the two polymorphic forms indicate that the yellow form is likely to be the most favorable thermodynamically. The dipole moments derived from experimental and theoretical charge densities and also from Lorentz tensor approach are compared with the single-molecule dipole moments. In each case, the differences of dipole moments between the polymorphs are identified.

Structural studies of the hen's egg lipovitellins. NMR and fluorescence investigations

The pH dependent reversible association-dissociation reaction of α- and β-lipovitellins from egg yolk has been studied by 1H NMR and fluorescence probe methods. Increased mobility of the choline methyl groups has been demonstrated on dissociation. The lipid methylene resonance of β-lipovitellin shows clear doublet character suggesting that the fatty acid chains exist in distinct environments. The high field component increases with temperature but is suppressed on treatment with pronase, suggesting a significant role for proteins in maintaining the differences in lipid environments. 1-Anilino-8-naphthalene sulfonate has been shown to bind less effectively to the monomeric lipovitellins. This is in agreement with earlier results suggesting that dissociation may be accompanied by increased hydration and conformational changes.

Analysis of the binding of polymyxin B to endotoxic lipid A and core glycolipid using a fluorescent displacement probe

Dansylcadaverine, a cationic fluorescent probe binds to bacterial lipopolysaccharide and lipid A, and is displaced competitively by other compounds which possess affinity toward endotoxins. The binding parameters of dansylcadaverine for lipid A were determined by Scatchard analysis to be two apparently equivalent sites with apparent dissociation constants (Kd) ranging between 16 μM to 26 μM, while that obtained for core glycolipid from Salmonella minnesota Re595 yielded a Kd of 22 μM to 28 μM with three binding sites. The Kd of polymyxin B for lipid A was computed from dansylcadaverine displacement by the method of Horovitz and Levitzki (Horovitz, A., and Levitzki, A. (1987) Proc. Natl. Acad. Sci. USA 84, 6654–6658). The applicability of this method for analyzing fluorescence data was validated by comparing the Kds of melittin for lipid A obtained by direct Scatchard analysis, and by the Horovitz-Levitzki method. The displacement of dansylcadaverine from lipid A by polymyxin B was distinctly biphasic with Kds for polymyxin B-lipid A interactions corresponding to 0.4 μM and 1.5 μM, probably resulting as a consequence of lipid A being a mixture of mono- and di-phosphoryl species. This was not observed with core glycolipid, for which the Kd for polymyxin was estimated to range from 1.1 μM to 5.8 μM. The use of dansylcadaverine as a displacement probe offers a novel and convenient method of quantitating the interactions of a wide variety of substances with lipid A.

Sequencing of T-superfamily conotoxins from Conus virgo: Pyroglutamic acid identification and disulfide arrangement by MALDI mass spectrometry

De novo mass spectrometric sequencing of two Conus peptides, Vi1359 and Vi1361, from the vermivorous cone snail Conus virgo, found off the southern Indian coast, is presented. The peptides, whose masses differ only by 2 Da, possess two disulfide bonds and an amidated C-terminus. Simple chemical modifications and enzymatic cleavage coupled with matrix assisted laser desorption ionization (MALDI) mass spectrometric analysis aided in establishing the sequences of Vi1359, ZCCITIPECCRI-NH2, and Vi1361, ZCCPTMPECCRI-NH2, Which differ only at residues 4 and 6 (Z = pyroglutamic acid). The presence of the pyroglutamyl residue at the N-terminus was unambiguously identified by chemical hydrolysis of the cyclic amide, followed by esterification. The presence of Ile residues in both the peptides was confirmed from high-energy collision induced dissociation (CID) studies, using the observation Of W-n- and d(n)-ions as a diagnostic. Differential cysteine labeling, in conjunction with MALDI-MS/MS, permitted establishment of disulfide connectivity in both peptides as Cys2-Cys9 and Cys3-Cys10. The cysteine pattern clearly reveals that the peptides belong to the class of T-superfamily conotoxins, in particular the T-1 superfamily.

Characterization of the interaction of lipid A and lipopolysaccharide with human serum albumin: implications for an endotoxin carrier function for albumin

The interactions of lipid A and lipopolysaccharide (LPS) with human serum albumin (HSA) were examined using fluorescence methods. Lipid A binds HSA with a stoichiometry of 2:1 with dissociation constants of 1.0 µM and 6.0 µM for the high- and low-affinity interactions, respectively. Lipid A displaces HSA-bound dansylsarcosine competitively, but not HSA-bound warfarin, suggesting that domain III-A, and not domain 11-A, is a lipid A binding site. Domain I does not contribute a site for lipid A. Based on these data, and the structural similarity between subdomains III-A and III-B, it is proposed that these two regions of HSA represent the high- and low-affinity sites of interaction of lipid A. Whole LPS also binds HSA, displacing dansylsarcosine, and its lipid A moiety appears to be the interaction site. However, there are differences between LPS and free lipid A. Polymyxin B forms ternary complexes with LPS bound to HSA, suggesting that the regions on LPS recognized by HSA and polymyxin B are different. The observed affinity of lipid A for HSA and mass action effects due to its abundance in the circulation would imply a major LPS carrier function for HSA.

Unfolding of Plasmodium falciparum Triosephosphate Isomerase in Urea and Guanidinium Chloride: Evidence for a Novel Disulfide Exchange Reaction in a Covalently Cross-Linked Mutant†

The conformational stability of Plasmodium falciparum triosephosphate isomerase (TIMWT) enzyme has been investigated in urea and guanidinium chloride (GdmCl) solutions using circular dichroism, fluorescence, and size-exclusion chromatography. The dimeric enzyme is remarkably stable in urea solutions. It retains considerable secondary, tertiary, and quaternary structure even in 8 M urea. In contrast, the unfolding transition is complete by 2.4 M GdmCl. Although the secondary as well as the tertiary interactions melt before the perturbation of the quaternary structure, these studies imply that the dissociation of the dimer into monomers ultimately leads to the collapse of the structure, suggesting that the interfacial interactions play a major role in determining multimeric protein stability. The C-m(urea)/C-m(GdmCl) ratio (where C-m is the concentration of the denaturant required at the transition midpoint) is unusually high for triosephosphate isomerase as compared to other monomeric and dimeric proteins. A disulfide crosslinked mutant protein (Y74C) engineered to form two disulfide cross-links across the interface (13-74') and (13'-74) is dramatically destablized in urea. The unfolding transition is complete by 6 M urea and involves a novel mechanism of dimer dissociation through intramolecular thiol-disulfide exchange.

Studies on the ionic transport and structural investigations of La0.5Li0.5TiO3 perovskite synthesized by wet chemical methods and the effect of Ce, Zr substitution at Ti site

La0.5Li0.5TiO3 perovskite was synthesized by various wet chemical methods. By adopting low temperature methods of preparation lithium loss from the material is prevented. La0.5Li0.5TiO3 (LLTO) was formed with cubic symmetry at 1473 K. LLTO was formed at relatively lower temperature by using hydrothermal preparation method. PVA gel-decomposition route yield tetragonal LLTO on annealing the dried gel at 1473 K. By using gel-carbonate route LiTi2O4 minor phase was found to remain even after heat-treatment at 1473 K. The hydroxylation of LLTO was done in deionized water as well as in dilute acetic acid medium. By hydroxylation process incorporation of hydroxyls and leaching out of Li+ was observed from the material. The Li+ concentration of these compositions was examined by AAS. The electrical conductivities of these compositions were measured by dc and ac impedance techniques at elevated temperatures. The activation energies of electrical conduction for these compositions were estimated from the experimental results. The measured activation energy of Li+ conduction is 0.34 eV. Unhydroxylated samples exhibit only Li+ conduction, whereas, the hydroxylated LLTO show proton conductivity at 298-550 K in addition to Li+ conductivity. The effect of Zr or Ce substitution in place of Ti were attempted. La0.5Li0.5ZrO3 Perovskite was not formed; instead pyrochlore phase (La2Zr2O7) along with monoclinic ZrO2 phases was observed above 1173 K; below 1173 K cubic ZrO2 is stable. (La0.5Li0.5)(2)CeO4 solid solution was formed in the case of Ce substitution at Ti sublattice on heat-treatment up to 1673 K. (c) 2005 Springer Science + Business Media, Inc.

Design of mechanism-based inhibitors of transthyretin amyloidosis: Studies with biphenyl ethers and new structural templates

Transthyretin (TTR), a tetrameric thyroxine (T4) carrier protein, is associated with a variety of amyloid diseases. In this study, we explore the potential of biphenyl ethers (BPE), which are shown to interact with a high affinity to its T4 binding site thereby preventing its aggregation and fibrillogenesis. They prevent fibrillogenesis by stabilizing the tetrameric ground state of transthyretin. Additionally, we identify two new structural templates (2-(5-mercapto-[1,3,4]oxadiazol-2-yl)-phenol and 2,3,6-trichloro-N-(4H-[1,2,4]triazol-3-yl) represented as compounds 11 and 12, respectively, throughout the manuscript) exhibiting the ability to arrest TTR amyloidosis. The dissociation constants for the binding of BPEs and compound 11 and 12 to TTR correlate with their efficacies of inhibiting amyloidosis. They also have the ability to inhibit the elongation of intermediate fibrils as well as show nearly complete (> 90%) disruption of the preformed fibrils. The present study thus establishes biphenyl ethers and compounds 11 and 12 as very potent inhibitors of TTR fibrillization and inducible cytotoxicity.

Photoswitchable multivalent sugar ligands: Synthesis, isomerization, and lectin binding studies of azobenzene-glycopyranoside derivatives

Coating of azobenzene chromophore with multivalent sugar ligands has been accomplished. Such sugar coating allows the study of the isomerization properties of this chromophore in aqueous solutions. The predominantly cis-isomer-containing photostationary state (PS) mixture of these azobenzene derivatives is found to be stable for hours. The rate constants for their isomerization, as well as the Arrhenius activation energies, are determined experimentally. An assessment of the lectin binding properties of the lactoside bearing isomeric azobenzene derivatives, by isothermal calorimetric methods, reveals the existence of an unusual cooperativity in their binding to lectin peanut agglutinin. Thermodynamic parameters evaluated for the trans and the PS mixture are discussed, in detail, for the lactoside bearing bivalent azobenzene derivative.