Multiscale study on hydrogen mobility in metallic fusion divertor material

For achieving efficient fusion energy production, the plasma-facing wall materials of the fusion reactor should ensure long time operation. In the next step fusion device, ITER, the first wall region facing the highest heat and particle load, i.e. the divertor area, will mainly consist of tiles based on tungsten. During the reactor operation, the tungsten material is slowly but inevitably saturated with tritium. Tritium is the relatively short-lived hydrogen isotope used in the fusion reaction. The amount of tritium retained in the wall materials should be minimized and its recycling back to the plasma must be unrestrained, otherwise it cannot be used for fueling the plasma. A very expensive and thus economically not viable solution is to replace the first walls quite often. A better solution is to heat the walls to temperatures where tritium is released. Unfortunately, the exact mechanisms of hydrogen release in tungsten are not known. In this thesis both experimental and computational methods have been used for studying the release and retention of hydrogen in tungsten. The experimental work consists of hydrogen implantations into pure polycrystalline tungsten, the determination of the hydrogen concentrations using ion beam analyses (IBA) and monitoring the out-diffused hydrogen gas with thermodesorption spectrometry (TDS) as the tungsten samples are heated at elevated temperatures. Combining IBA methods with TDS, the retained amount of hydrogen is obtained as well as the temperatures needed for the hydrogen release. With computational methods the hydrogen-defect interactions and implantation-induced irradiation damage can be examined at the atomic level. The method of multiscale modelling combines the results obtained from computational methodologies applicable at different length and time scales. Electron density functional theory calculations were used for determining the energetics of the elementary processes of hydrogen in tungsten, such as diffusivity and trapping to vacancies and surfaces. Results from the energetics of pure tungsten defects were used in the development of an classical bond-order potential for describing the tungsten defects to be used in molecular dynamics simulations. The developed potential was utilized in determination of the defect clustering and annihilation properties. These results were further employed in binary collision and rate theory calculations to determine the evolution of large defect clusters that trap hydrogen in the course of implantation. The computational results for the defect and trapped hydrogen concentrations were successfully compared with the experimental results. With the aforedescribed multiscale analysis the experimental results within this thesis and found in the literature were explained both quantitatively and qualitatively.

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.

Radiation effects in supported nanoparticles

Nanotechnology applications are entering the market in increasing numbers, nanoparticles being among the main classes of materials used. Particles can be used, e.g., for catalysing chemical reactions, such as is done in car exhaust catalysts today. They can also modify the optical and electronic properties of materials or be used as building blocks for thin film coatings on a variety of surfaces. To develop materials for specific applications, an intricate control of the particle properties, structure, size and shape is required. All these depend on a multitude of factors from methods of synthesis and deposition to post-processing. This thesis addresses the control of nanoparticle structure by low-energy cluster beam deposition and post-synthesis ion irradiation. Cluster deposition in high vacuum offers a method for obtaining precisely controlled cluster-assembled materials with minimal contamination. Due to the clusters small size, however, the cluster-surface interaction may drastically change the cluster properties on deposition. In this thesis, the deposition process of metal and alloy clusters on metallic surfaces is modelled using molecular dynamics simulations, and the mechanisms influencing cluster structure are identified. Two mechanisms, mechanical melting upon deposition and thermally activated dislocation motion, are shown to determine whether a deposited cluster will align epitaxially with its support. The semiconductor industry has used ion irradiation as a tool to modify material properties for decades. Irradiation can be used for doping, patterning surfaces, and inducing chemical ordering in alloys, just to give a few examples. The irradiation response of nanoparticles has, however, remained an almost uncharted territory. Although irradiation effects in nanoparticles embedded inside solid matrices have been studied, almost no work has been done on supported particles. In this thesis, the response of supported nanoparticles is studied systematically for heavy and light ion irradiation. The processes leading to damage production are identified and models are developed for both types of irradiation. In recent experiments, helium irradiation has been shown to induce a phase transformation from multiply twinned to single-crystalline nanoparticles in bimetallic alloys, but the nature of the transition has remained unknown. The alloys for which the effect has been observed are CuAu and FePt. It is shown in this thesis that transient amorphization leads to the observed transition and that while CuAu and FePt do not amorphize upon irradiation in bulk or as thin films, they readily do so as nanoparticles. This is the first time such an effect is demonstrated with supported particles, not embedded in a matrix where mixing is always an issue. An understanding of the above physical processes is essential, if nanoparticles are to be used in applications in an optimal way. This thesis clarifies the mechanisms which control particle morphology, and paves way for the synthesis of nanostructured materials tailored for specific applications.

Local structure of water studied by Compton scattering

This thesis presents a novel application of x-ray Compton scattering to structural studies of molecular liquids. Systematic Compton-scattering experiments on water have been carried out with unprecedented accuracy at third-generation synchrotron-radiation laboratories. The experiments focused on temperature effects in water, the water-to-ice phase transition, quantum isotope effects, and ion hydration. The experimental data is interpreted by comparison with both model computations and ab initio molecular-dynamics simulations. Accordingly, Compton scattering is found to provide unique intra- and intermolecular structural information. This thesis thus demonstrates the complementarity of the technique to traditional real-space probes for studies on the local structure of water and, more generally, molecular liquids.

The 5-Methyl Group in Thymine Dynamically Influences the Structure of A-Tracts in DNA at the Local and Global Level

DNA sequences containing a stretch of several A:T basepairs without a 5'-TA-3' step are known as A-tracts and have been the subject of extensive investigation because of their unique structural features such as a narrow minor groove and their crucial role in several biological processes. One of the aspects under investigation has been the influence of the 5-methyl group of thymine on the properties of A-tracts. Detailed molecular dynamics simulation studies of the sequences d(CGCAAAUUUGCG) and d(CGCAAATTTGCG) indicate that the presence of the 5-methyl group in thymine increases the frequency of a narrow minor groove conformation, which could facilitate its specific recognition by proteins, and reduce its susceptibility to cleavage by DNase I. The bias toward a wider minor groove in the absence of the thymine 5-methyl group is a static structural feature. Our results also indicate that the presence of the thymine 5-methyl group is necessary for calibrating the backbone conformation and the basepair and dinucleotide step geometry of the core A-tract as well as the flanking CA/TG and the neighboring GC/GC steps, as observed in free and protein-bound DNA. As a consequence, it also fine-tunes the curvature of the longer DNA fragment in which the A-tract is embedded.

The Path to a Porous Semiconductor Multilayer

Thin films are the basis of much of recent technological advance, ranging from coatings with mechanical or optical benefits to platforms for nanoscale electronics. In the latter, semiconductors have been the norm ever since silicon became the main construction material for a multitude of electronical components. The array of characteristics of silicon-based systems can be widened by manipulating the structure of the thin films at the nanoscale - for instance, by making them porous. The different characteristics of different films can then to some extent be combined by simple superposition. Thin films can be manufactured using many different methods. One emerging field is cluster beam deposition, where aggregates of hundreds or thousands of atoms are deposited one by one to form a layer, the characteristics of which depend on the parameters of deposition. One critical parameter is deposition energy, which dictates how porous, if at all, the layer becomes. Other parameters, such as sputtering rate and aggregation conditions, have an effect on the size and consistency of the individual clusters. Understanding nanoscale processes, which cannot be observed experimentally, is fundamental to optimizing experimental techniques and inventing new possibilities for advances at this scale. Atomistic computer simulations offer a window to the world of nanometers and nanoseconds in a way unparalleled by the most accurate of microscopes. Transmission electron microscope image simulations can then bridge this gap by providing a tangible link between the simulated and the experimental. In this thesis, the entire process of cluster beam deposition is explored using molecular dynamics and image simulations. The process begins with the formation of the clusters, which is investigated for Si/Ge in an Ar atmosphere. The structure of the clusters is optimized to bring it as close to the experimental ideal as possible. Then, clusters are deposited, one by one, onto a substrate, until a sufficiently thick layer has been produced. Finally, the concept is expanded by further deposition with different parameters, resulting in multiple superimposed layers of different porosities. This work demonstrates how the aggregation of clusters is not entirely understood within the scope of the approximations used in the simulations; yet, it is also shown how the continued deposition of clusters with a varying deposition energy can lead to a novel kind of nanostructured thin film: a multielemental porous multilayer. According to theory, these new structures have characteristics that can be tailored for a variety of applications, with precision heretofore unseen in conventional multilayer manufacture.

Structures of apo and GTP-bound molybdenum cofactor biosynthesis protein MoaC from Thermus thermophilus HB8

The first step in the molybdenum cofactor (Moco) biosynthesis pathway involves the conversion of guanosine triphosphate (GTP) to precursor Z by two proteins (MoaA and MoaC). MoaA belongs to the S-adenosylmethioninedependent radical enzyme superfamily and is believed to generate protein and/or substrate radicals by reductive cleavage of S-adenosylmethionine using an Fe-S cluster. MoaC has been suggested to catalyze the release of pyrophosphate and the formation of the cyclic phosphate of precursor Z. However, structural evidence showing the binding of a substrate-like molecule to MoaC is not available. Here, apo and GTP-bound crystal structures of MoaC from Thermus thermophilus HB8 are reported. Furthermore, isothermal titration calorimetry experiments have been carried out in order to obtain thermodynamic parameters for the protein-ligand interactions. In addition, molecular-dynamics (MD) simulations have been carried out on the protein-ligand complex of known structure and on models of relevant complexes for which X-ray structures are not available. The biophysical, structural and MD results reveal the residues that are involved in substrate binding and help in speculating upon a possible mechanism.

Single-File Diffusion of Confined Water Inside SWNTs: An NMR Study

We report a nuclear magnetic resonance (NMR) study of confined water inside similar to 1.4 nm diameter single-walled carbon nanotubes (SWNTs). We show that the confined water does not freeze even up to 223 K. A pulse field gradient (PFG) NMR method is used to determine the mean squared displacement (MSD) of the water molecules inside the nanotubes at temperatures below 273 K, where the bulk water outside the nanotubes freezes and hence does not contribute to the proton NMR signal. We show that the mean squared displacement varies as the square root of time, predicted for single-file diffusion in a one-dimensional channel. We propose a qualitative understanding of our results based on available molecular dynamics simulations.

Single-File Diffusion of Water Inside Narrow Carbon Nanorings

We use atomistic molecular dynamics (MD) simulations to study the diffusion of water molecules confined inside narrow (6,6) carbon nanorings. The water molecules form two oppositely polarized chains. It is shown that the effective interaction between these two chains is repulsive in nature. The computed mean-squared displacement (MSD) clearly shows a scaling with time similar to t(1/2), which is consistent with single-file diffusion (SFD). The time up to which the water molecules undergo SFD is shown to be the lifetime of the water molecules inside these chains. Simulations of "uncharged" water molecules inside the nanoring show the formation of several water chains and yield SFD. These observations conclusively prove that the diffusion is Fickian when there is a single chain of water and SFD is observed only when two or more chains are present.

Hydration of ions under confinement

Molecular dynamics simulations are used to examine the changes in water density and hydration characteristics of NaCl solutions confined in slit-shaped graphitic pores. Using a structural signature, we define the hydration limit as the salt concentration at which a sharp drop in the hydration number is observed. At small pores (H = 8.0-10 angstrom), confined water does not possess bulk-like features and remains in a layered arrangement between two surfaces. Despite this high degree of confinement, ions are able to form a quasi-2D hydration shell between two surfaces. Our results indicate the strong propensity of ions to form the first hydration shell, even under extremely confined aqueous environments. The hydration of ions is seen to weakly perturb the oxygen density distributions between two surfaces. The hydration number of Na+ reduces to about 4.15 at a pore width of H = 0.8 nm, when compared with the bulk hydration number of 6.25. At larger pore widths, above H = 16 angstrom, where bulk-like water densities are observed in the central regions of the pore, the hydration number is above 6.

Investigation of the effect of nonlocal scale on ultrasonic wave dispersion characteristics of a monolayer graphene

In this article, an ultrasonic wave propagation in graphene sheet is studied using nonlocal elasticity theory incorporating small scale effects. The graphene sheet is modeled as an isotropic plate of one-atom thick. For this model, the nonlocal governing differential equations of motion are derived from the minimization of the total potential energy of the entire system. An ultrasonic type of wave propagation model is also derived for the graphene sheet. The nonlocal scale parameter introduces certain band gap region in in-plane and flexural wave modes where no wave propagation occurs. This is manifested in the wavenumber plots as the region where the wavenumber tends to infinite or wave speed tends to zero. The frequency at which this phenomenon occurs is called the escape frequency. The explicit expressions for cutoff frequencies and escape frequencies are derived. The escape frequencies are mainly introduced because of the nonlocal elasticity. Obviously these frequencies are function of nonlocal scaling parameter. It has also been obtained that these frequencies are independent of y-directional wavenumber. It means that for any type of nanostructure, the escape frequencies are purely a function of nonlocal scaling parameter only. It is also independent of the geometry of the structure. It has been found that the cutoff frequencies are function of nonlocal scaling parameter (e(0)a) and the y-directional wavenumber (k(y)). For a given nanostructure, nonlocal small scale coefficient can be obtained by matching the results from molecular dynamics (MD) simulations and the nonlocal elasticity calculations. At that value of the nonlocal scale coefficient, the waves will propagate in the nanostructure at that cut-off frequency. In the present paper, different values of e(o)a are used. One can get the exact e(0)a for a given graphene sheet by matching the MD simulation results of graphene with the results presented in this paper. (C) 2010 Elsevier B.V. All rights reserved.

Irradiation effects in graphene and related materials

Nanomaterials with a hexagonally ordered atomic structure, e.g., graphene, carbon and boron nitride nanotubes, and white graphene (a monolayer of hexagonal boron nitride) possess many impressive properties. For example, the mechanical stiffness and strength of these materials are unprecedented. Also, the extraordinary electronic properties of graphene and carbon nanotubes suggest that these materials may serve as building blocks of next generation electronics. However, the properties of pristine materials are not always what is needed in applications, but careful manipulation of their atomic structure, e.g., via particle irradiation can be used to tailor the properties. On the other hand, inadvertently introduced defects can deteriorate the useful properties of these materials in radiation hostile environments, such as outer space. In this thesis, defect production via energetic particle bombardment in the aforementioned materials is investigated. The effects of ion irradiation on multi-walled carbon and boron nitride nanotubes are studied experimentally by first conducting controlled irradiation treatments of the samples using an ion accelerator and subsequently characterizing the induced changes by transmission electron microscopy and Raman spectroscopy. The usefulness of the characterization methods is critically evaluated and a damage grading scale is proposed, based on transmission electron microscopy images. Theoretical predictions are made on defect production in graphene and white graphene under particle bombardment. A stochastic model based on first-principles molecular dynamics simulations is used together with electron irradiation experiments for understanding the formation of peculiar triangular defect structures in white graphene. An extensive set of classical molecular dynamics simulations is conducted, in order to study defect production under ion irradiation in graphene and white graphene. In the experimental studies the response of carbon and boron nitride multi-walled nanotubes to irradiation with a wide range of ion types, energies and fluences is explored. The stabilities of these structures under ion irradiation are investigated, as well as the issue of how the mechanism of energy transfer affects the irradiation-induced damage. An irradiation fluence of 5.5x10^15 ions/cm^2 with 40 keV Ar+ ions is established to be sufficient to amorphize a multi-walled nanotube. In the case of 350 keV He+ ion irradiation, where most of the energy transfer happens through inelastic collisions between the ion and the target electrons, an irradiation fluence of 1.4x10^17 ions/cm^2 heavily damages carbon nanotubes, whereas a larger irradiation fluence of 1.2x10^18 ions/cm^2 leaves a boron nitride nanotube in much better condition, indicating that carbon nanotubes might be more susceptible to damage via electronic excitations than their boron nitride counterparts. An elevated temperature was discovered to considerably reduce the accumulated damage created by energetic ions in both carbon and boron nitride nanotubes, attributed to enhanced defect mobility and efficient recombination at high temperatures. Additionally, cobalt nanorods encapsulated inside multi-walled carbon nanotubes were observed to transform into spherical nanoparticles after ion irradiation at an elevated temperature, which can be explained by the inverse Ostwald ripening effect. The simulation studies on ion irradiation of the hexagonal monolayers yielded quantitative estimates on types and abundances of defects produced within a large range of irradiation parameters. He, Ne, Ar, Kr, Xe, and Ga ions were considered in the simulations with kinetic energies ranging from 35 eV to 10 MeV, and the role of the angle of incidence of the ions was studied in detail. A stochastic model was developed for utilizing the large amount of data produced by the molecular dynamics simulations. It was discovered that a high degree of selectivity over the types and abundances of defects can be achieved by carefully selecting the irradiation parameters, which can be of great use when precise pattering of graphene or white graphene using focused ion beams is planned.

Elucidation of the conformational free energy landscape in H.pylori LuxS and its implications to catalysis

Background: One of the major challenges in understanding enzyme catalysis is to identify the different conformations and their populations at detailed molecular level in response to ligand binding/environment. A detail description of the ligand induced conformational changes provides meaningful insights into the mechanism of action of enzymes and thus its function. Results: In this study, we have explored the ligand induced conformational changes in H. pylori LuxS and the associated mechanistic features. LuxS, a dimeric protein, produces the precursor (4,5-dihydroxy-2,3-pentanedione) for autoinducer-2 production which is a signalling molecule for bacterial quorum sensing. We have performed molecular dynamics simulations on H. pylori LuxS in its various ligand bound forms and analyzed the simulation trajectories using various techniques including the structure network analysis, free energy evaluation and water dynamics at the active site. The results bring out the mechanistic details such as co operativity and asymmetry between the two subunits, subtle changes in the conformation as a response to the binding of active and inactive forms of ligands and the population distribution of different conformations in equilibrium. These investigations have enabled us to probe the free energy landscape and identify the corresponding conformations in terms of network parameters. In addition, we have also elucidated the variations in the dynamics of water co-ordination to the Zn2+ ion in LuxS and its relation to the rigidity at the active sites. Conclusions: In this article, we provide details of a novel method for the identification of conformational changes in the different ligand bound states of the protein, evaluation of ligand-induced free energy changes and the biological relevance of our results in the context of LuxS structure-function. The methodology outlined here is highly generalized to illuminate the linkage between structure and function in any protein of known structure.

Tetrahydrofuran Clathrate Hydrate Formation

We report on the formation of tetrahydrofuran clathrate hydrate studied by x-ray Raman scattering measurements at the oxygen K edge. A comparison of x-ray Raman spectra measured from water-tetrahydrofuran mixtures and tetrahydrofuran hydrate at different temperatures supports stochastic hydrate formation models rather than models assuming hydrate precursors. This is confirmed by molecular dynamics simulations and density functional theory calculations of x-ray Raman spectra. In addition, changes in the spectra of tetrahydrofuran hydrate with temperatures close to the hydrate's dissociation temperature were observed and may be connected to changes in hydrate's local structure due to the formation of hydrogen bonds between guest and water molecules.

Ions in water: Role of attractive interactions in size dependent diffusivity maximum

A molecular dynamics study of model ions in water is reported. The van der Waals diameter of both the cations and anions is varied. We have carried out two sets of simulations-with and without dispersion interaction-between the ion and water. Self-diffusivity of the ions exhibits an anomalous maximum as a function of the van der Waals diameter for both these sets. This existence of a maximum in self-diffusivity when there is no dispersion interaction between the ion and the water is attributed to the attractive term from electrostatic interactions. Detailed analysis of this effect shows that the solvent shell is more strongly defined in the presence of dispersion interactions. A smaller ion exhibits biexponential decay while a single exponential decay is seen for the ion with maximum diffusivity in the self-part of the intermediate scattering function. The solvent structure around the ion appears to determine much of the dynamics of the ion. Interesting trends are seen in the activation energies and these can be understood in terms of the levitation effect. (C) 2010 American Institute of Physics. doi:10.1063/1.3481656]