Diffusion of Be, Co and Mn impurities in compound semiconductors and glassy carbon studied by the modified radiotracer technique

The main method of modifying properties of semiconductors is to introduce small amount of impurities inside the material. This is used to control magnetic and optical properties of materials and to realize p- and n-type semiconductors out of intrinsic material in order to manufacture fundamental components such as diodes. As diffusion can be described as random mixing of material due to thermal movement of atoms, it is essential to know the diffusion behavior of the impurities in order to manufacture working components. In modified radiotracer technique diffusion is studied using radioactive isotopes of elements as tracers. The technique is called modified as atoms are deployed inside the material by ion beam implantation. With ion implantation, a distinct distribution of impurities can be deployed inside the sample surface with good con- trol over the amount of implanted atoms. As electromagnetic radiation and other nuclear decay products emitted by radioactive materials can be easily detected, only very low amount of impurities can be used. This makes it possible to study diffusion in pure materials without essentially modifying the initial properties by doping. In this thesis a modified radiotracer technique is used to study the diffusion of beryllium in GaN, ZnO, SiGe and glassy carbon. GaN, ZnO and SiGe are of great interest to the semiconductor industry and beryllium as a small and possibly rapid dopant hasn t been studied previously using the technique. Glassy carbon has been added to demonstrate the feasibility of the technique. In addition, the diffusion of magnetic impurities, Mn and Co, has been studied in GaAs and ZnO (respectively) with spintronic applications in mind.

Interatomic potentials for fusion reactor material simulations

Fusion energy is a clean and safe solution for the intricate question of how to produce non-polluting and sustainable energy for the constantly growing population. The fusion process does not result in any harmful waste or green-house gases, since small amounts of helium is the only bi-product that is produced when using the hydrogen isotopes deuterium and tritium as fuel. Moreover, deuterium is abundant in seawater and tritium can be bred from lithium, a common metal in the Earth's crust, rendering the fuel reservoirs practically bottomless. Due to its enormous mass, the Sun has been able to utilize fusion as its main energy source ever since it was born. But here on Earth, we must find other means to achieve the same. Inertial fusion involving powerful lasers and thermonuclear fusion employing extreme temperatures are examples of successful methods. However, these have yet to produce more energy than they consume. In thermonuclear fusion, the fuel is held inside a tokamak, which is a doughnut-shaped chamber with strong magnets wrapped around it. Once the fuel is heated up, it is controlled with the help of these magnets, since the required temperatures (over 100 million degrees C) will separate the electrons from the nuclei, forming a plasma. Once the fusion reactions occur, excess binding energy is released as energetic neutrons, which are absorbed in water in order to produce steam that runs turbines. Keeping the power losses from the plasma low, thus allowing for a high number of reactions, is a challenge. Another challenge is related to the reactor materials, since the confinement of the plasma particles is not perfect, resulting in particle bombardment of the reactor walls and structures. Material erosion and activation as well as plasma contamination are expected. Adding to this, the high energy neutrons will cause radiation damage in the materials, causing, for instance, swelling and embrittlement. In this thesis, the behaviour of a material situated in a fusion reactor was studied using molecular dynamics simulations. Simulations of processes in the next generation fusion reactor ITER include the reactor materials beryllium, carbon and tungsten as well as the plasma hydrogen isotopes. This means that interaction models, {\it i.e. interatomic potentials}, for this complicated quaternary system are needed. The task of finding such potentials is nonetheless nearly at its end, since models for the beryllium-carbon-hydrogen interactions were constructed in this thesis and as a continuation of that work, a beryllium-tungsten model is under development. These potentials are combinable with the earlier tungsten-carbon-hydrogen ones. The potentials were used to explain the chemical sputtering of beryllium due to deuterium plasma exposure. During experiments, a large fraction of the sputtered beryllium atoms were observed to be released as BeD molecules, and the simulations identified the swift chemical sputtering mechanism, previously not believed to be important in metals, as the underlying mechanism. Radiation damage in the reactor structural materials vanadium, iron and iron chromium, as well as in the wall material tungsten and the mixed alloy tungsten carbide, was also studied in this thesis. Interatomic potentials for vanadium, tungsten and iron were modified to be better suited for simulating collision cascades that are formed during particle irradiation, and the potential features affecting the resulting primary damage were identified. Including the often neglected electronic effects in the simulations was also shown to have an impact on the damage. With proper tuning of the electron-phonon interaction strength, experimentally measured quantities related to ion-beam mixing in iron could be reproduced. The damage in tungsten carbide alloys showed elemental asymmetry, as the major part of the damage consisted of carbon defects. On the other hand, modelling the damage in the iron chromium alloy, essentially representing steel, showed that small additions of chromium do not noticeably affect the primary damage in iron. Since a complete assessment of the response of a material in a future full-scale fusion reactor is not achievable using only experimental techniques, molecular dynamics simulations are of vital help. This thesis has not only provided insight into complicated reactor processes and improved current methods, but also offered tools for further simulations. It is therefore an important step towards making fusion energy more than a future goal.

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.

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.

Accelerator Mass Spectrometry and Bayesian Data Analysis

Accelerator mass spectrometry (AMS) is an ultrasensitive technique for measuring the concentration of a single isotope. The electric and magnetic fields of an electrostatic accelerator system are used to filter out other isotopes from the ion beam. The high velocity means that molecules can be destroyed and removed from the measurement background. As a result, concentrations down to one atom in 10^16 atoms are measurable. This thesis describes the construction of the new AMS system in the Accelerator Laboratory of the University of Helsinki. The system is described in detail along with the relevant ion optics. System performance and some of the 14C measurements done with the system are described. In a second part of the thesis, a novel statistical model for the analysis of AMS data is presented. Bayesian methods are used in order to make the best use of the available information. In the new model, instrumental drift is modelled with a continuous first-order autoregressive process. This enables rigorous normalization to standards measured at different times. The Poisson statistical nature of a 14C measurement is also taken into account properly, so that uncertainty estimates are much more stable. It is shown that, overall, the new model improves both the accuracy and the precision of AMS measurements. In particular, the results can be improved for samples with very low 14C concentrations or measured only a few times.

Development of analytical techniques for studies on dispersion of actinides in the environment and characterization of environmental radioactive particles

Radioactive particles from three locations were investigated for elemental composition, oxidation states of matrix elements, and origin. Instrumental techniques applied to the task were scanning electron microscopy, X-ray and gamma-ray spectrometry, secondary ion mass spectrometry, and synchrotron radiation based microanalytical techniques comprising X-ray fluorescence spectrometry, X-ray fluorescence tomography, and X-ray absorption near-edge structure spectroscopy. Uranium-containing low activity particles collected from Irish Sea sediments were characterized in terms of composition and distribution of matrix elements and the oxidation states of uranium. Indications of the origin were obtained from the intensity ratios and the presence of thorium, uranium, and plutonium. Uranium in the particles was found to exist mostly as U(IV). Studies on plutonium particles from Runit Island (Marshall Islands) soil indicated that the samples were weapon fuel fragments originating from two separate detonations: a safety test and a low-yield test. The plutonium in the particles was found to be of similar age. The distribution and oxidation states of uranium and plutonium in the matrix of weapon fuel particles from Thule (Greenland) sediments were investigated. The variations in intensity ratios observed with different techniques indicated more than one origin. Uranium in particle matrixes was mostly U(IV), but plutonium existed in some particles mainly as Pu(IV), and in others mainly as oxidized Pu(VI). The results demonstrated that the various techniques were effectively applied in the characterization of environmental radioactive particles. An on-line method was developed for separating americium from environmental samples. The procedure utilizes extraction chromatography to separate americium from light lanthanides, and cation exchange to concentrate americium before the final separation in an ion chromatography column. The separated radiochemically pure americium fraction is measured by alpha spectrometry. The method was tested with certified sediment and soil samples and found to be applicable for the analysis of environmental samples containing a wide range of Am-241 activity. Proceeding from the on-line method developed for americium, a method was also developed for separating plutonium and americium. Plutonium is reduced to Pu(III), and separated together with Am(III) throughout the procedure. Pu(III) and Am(III) are eluted from the ion chromatography column as anionic dipicolinate and oxalate complexes, respectively, and measured by alpha spectrometry.

A quantum chemical investigation of the metal centres in cytochrome c oxidase

Quantum effects are often of key importance for the function of biological systems at molecular level. Cellular respiration, where energy is extracted from the reduction of molecular oxygen to water, is no exception. In this work, the end station of the electron transport chain in mitochondria, cytochrome c oxidase, is investigated using quantum chemical methodology. Cytochrome c oxidase contains two haems, haem a and haem a3. Haem a3, with its copper companion, CuB, is involved in the final reduction of oxygen into water. This binuclear centre receives the necessary electrons from haem a. Haem a, in turn, receives its electrons from a copper ion pair in the vicinity, called CuA. Density functional theory (DFT) has been used to clarify the charge and spin distributions of haem a, as well as changes in these during redox activity. Upon reduction, the added electron is shown to be evenly distributed over the entire haem structure, important for the accommodation of the prosthetic group within the protein. At the same time, the spin distribution of the open-shell oxidised state is more localised to the central iron. The exact spin density distribution has been disputed in the literature, however, different experiments indicating different distributions of the unpaired electron. The apparent contradiction is shown to be due to the false assumption of a unit amount of unpaired electron density; in fact, the oxidised state has about 1.3 unpaired electrons. The validity of the DFT results have been corroborated by wave function based coupled cluster calculations. Point charges, for use in classical force field based simulations, have been parameterised for the four metal centres, using a newly developed methodology. In the procedure, the subsystem for which point charges are to be obtained, is surrounded by an outer region, with the purpose of stabilising the inner region, both electronically and structurally. Finally, the possibility of vibrational promotion of the electron transfer step between haem a and a3 has been investigated. Calculating the full vibrational spectra, at DFT level, of a combined model of the two haems, revealed several normal modes that do shift electron density between the haems. The magnitude of the shift was found to be moderate, at most. The proposed mechanism could have an assisting role in the electron transfer, which still seems to be dominated by electron tunnelling.

Redox-linked proton transfer by cytochrome c oxidase

The respiratory chain is found in the inner mitochondrial membrane of higher organisms and in the plasma membrane of many bacteria. It consists of several membrane-spanning enzymes, which conserve the energy that is liberated from the degradation of food molecules as an electrochemical proton gradient across the membrane. The proton gradient can later be utilized by the cell for different energy requiring processes, e.g. ATP production, cellular motion or active transport of ions. The difference in proton concentration between the two sides of the membrane is a result of the translocation of protons by the enzymes of the respiratory chain, from the negatively charged (N-side) to the positively charged side (P-side) of the lipid bilayer, against the proton concentration gradient. The endergonic proton transfer is driven by the flow of electrons through the enzymes of the respiratory chain, from low redox-potential electron donors to acceptors of higher potential, and ultimately to oxygen. Cytochrome c oxidase is the last enzyme in the respiratory chain and catalyzes the reduction of dioxygen to water. The redox reaction is coupled to proton transport across the membrane by a yet unresolved mechanism. Cytochrome c oxidase has two proton-conducting pathways through which protons are taken up to the interior part of the enzyme from the N-side of the membrane. The K-pathway transfers merely substrate protons, which are consumed in the process of water formation at the catalytic site. The D-pathway transfers both substrate protons and protons that are pumped to the P-side of the membrane. This thesis focuses on the role of two conserved amino acids in proton translocation by cytochrome c oxidase, glutamate 278 and tryptophan 164. Glu278 is located at the end of the D-pathway and is thought to constitute the branching point for substrate and pumped protons. In this work, it was shown that although Glu278 has an important role in the proton transfer mechanism, its presence is not an obligatory requirement. Alternative structural solutions in the area around Glu278, much like the ones present in some distantly related heme-copper oxidases, could in the absence of Glu278 support the formation of a long hydrogen-bonded water chain through which proton transfer from the D-pathway to the catalytic site is possible. The other studied amino acid, Trp164, is hydrogen bonded to the ∆-propionate of heme a3 of the catalytic site. Mutation of this amino acid showed that it may be involved in regulation of proton access to a proton acceptor, a pump site, from which the proton later is expelled to the P-side of the membrane. The ion pair that is formed by the ∆-propionate of heme a3 and arginine 473 is likely to form a gate-like structure, which regulates proton mobility to the P-side of the membrane. The same gate may also be part of an exit path through which water molecules produced at the catalytically active site are removed towards the external side of the membrane. Time-resolved optical and electrometrical experiments with the Trp164 to phenylalanine mutant revealed a so far undetected step in the proton pumping mechanism. During the A to PR transition of the catalytic cycle, a proton is transferred from Glu278 to the pump site, located somewhere in the vicinity of the ∆-propionate of heme a3. A mechanism for proton pumping by cytochrome c oxidase is proposed on the basis of the presented results and the mechanism is discussed in relation to some relevant experimental data. A common proton pumping mechanism for all members of the heme-copper oxidase family is moreover considered.

Development of atmospheric pressure ionization ion mobility spectrometry and ion mobility spectrometry mass spectrometry

This study is focused on the development and evaluation of ion mobility instrumentation with various atmospheric pressure ionization techniques and includes the following work. First, a high-resolution drift tube ion mobility spectrometer (IMS), coupled with a commercial triple quadrupole mass spectrometer (MS), was developed. This drift tube IMS is compatible with the front-end of commercial Sciex mass spectrometers (e.g., Sciex API-300, 365, and 3000) and also allows easy (only minor modifications are needed) installation between the original atmospheric pressure ion source and the triple quadrupole mass spectrometer. Performance haracteristics (e.g.,resolving power, detection limit, transmission efficiency of ions) of this IMS-MS instrument were evaluated. Development of the IMS-MS instrument also led to a study where a proposal was made that tetraalkylammonium ions can be used as chemical standards for ESI-IMS. Second, the same drift tube design was also used to build a standalone ion mobility spectrometer equipped with a Faraday plate detector. For this highresolution (resolving power about 100 shown) IMS device, a multi-ion source platform was built, which allows the use of a range of atmospheric pressure ionization methods, such as: corona discharge chemical ionization (CD-APCI), atmospheric pressure photoionization (APPI), and radioactive atmospheric pressure chemical ionization (R-APCI). The multi-ion source platform provides easy switching between ionization methods and both positive and negative ionization modes can be used. Third, a simple desorpion/ionization on silicon (DIOS) ion source set-up for use with the developed IMS and IMS-MS instruments was built and its operation demonstrated. Fourth, a prototype of a commercial aspiration-type ion mobility spectrometer was mounted in front of a commercial triple quadrupole mass spectrometer. The set-up, which is simple, easy to install, and requires no major modifications to the MS, provides the possibility of gathering fundamental information about aspiration mobility spectrometry.