Computer simulation of multi-elemental fusion reactor materials

Thermonuclear fusion is a sustainable energy solution, in which energy is produced using similar processes as in the sun. In this technology hydrogen isotopes are fused to gain energy and consequently to produce electricity. In a fusion reactor hydrogen isotopes are confined by magnetic fields as ionized gas, the plasma. Since the core plasma is millions of degrees hot, there are special needs for the plasma-facing materials. Moreover, in the plasma the fusion of hydrogen isotopes leads to the production of high energetic neutrons which sets demanding abilities for the structural materials of the reactor. This thesis investigates the irradiation response of materials to be used in future fusion reactors. Interactions of the plasma with the reactor wall leads to the removal of surface atoms, migration of them, and formation of co-deposited layers such as tungsten carbide. Sputtering of tungsten carbide and deuterium trapping in tungsten carbide was investigated in this thesis. As the second topic the primary interaction of the neutrons in the structural material steel was examined. As model materials for steel iron chromium and iron nickel were used. This study was performed theoretically by the means of computer simulations on the atomic level. In contrast to previous studies in the field, in which simulations were limited to pure elements, in this work more complex materials were used, i.e. they were multi-elemental including two or more atom species. The results of this thesis are in the microscale. One of the results is a catalogue of atom species, which were removed from tungsten carbide by the plasma. Another result is e.g. the atomic distributions of defects in iron chromium caused by the energetic neutrons. These microscopic results are used in data bases for multiscale modelling of fusion reactor materials, which has the aim to explain the macroscopic degradation in the materials. This thesis is therefore a relevant contribution to investigate the connection of microscopic and macroscopic radiation effects, which is one objective in fusion reactor materials research.

Syntheses and Self-Assembling Characteristics of Amphiphilic Star Diblock Copolymers

Uusien polymeeripohjaisten teknologioiden ja materiaalien myötä räätälöityjen polymeerien tarve on kasvanut. Viime vuosituhannen lopussa kehitetyt kontrolloidut polymerointimenetelmät ovat avanneet uusia mahdollisuuksia paitsi monimutkaisten polymeerien synteesiin, myös itsejärjestyvyyteen perustuvien funktionaalisten nanorakenteiden suunnitteluun ja valmistukseen. Nämä voivat jäljitellä luonnossa esiintyviä rakenteita, joita muodostavat esimerkiksi lipidit ja proteiinit. Itsejärjestyvät molekyylit ovat usein amfifiilisiä eli ne koostuvat hydrofiilisistä ja hydrofobisista osista ja polymeereissä nämä osat voivat olla omina lohkoinaan, jolloin puhutaan amfifiilisistä lohko- tai blokkikopolymeereistä. Riippuen järjestyneiden rakenteiden koostumuksesta ja muodosta, amfifiilisiä blokkikopolymeerejä on tutkittu tai jo käytetty nanoteknologiassa, elastomeereissä, voiteluaineissa, pinta-aktiivisina aineina, lääkkeenannostelussa, maaleissa, sekä elektroniikka-, kosmetiikka- ja elintarviketeollisuudessa. Tavallisimmin käytetyt amfifiiliset blokkikopolymeerit ovat olleet lineaarisia, mutta viime aikoina tutkimus on suuntautunut kohti monimutkaisempia rakenteita. Tällaisia ovat esimerkiksi tähtipolymeerit. Tähtimäisissä polymeereissä miselleille tyypillinen ydin-kuori-rakenne säilyy hyvin alhaisissakin polymeerikonsentraatioissa, koska polymeeriketjut ovat kiinni toisissaan yhdessä pisteessä. Siten ne ovat erityisen kiinnostavia tutkimuskohteita erilaisten hydrofobisten orgaanisten yhdisteiden sitomiseksi ja vapauttamiseksi. Tässä työssä on tarkasteltu amfifiilisten tähtipolymeerien itsejärjestymistä vesiliuoksissa sekä kokeellisesti ja tietokonesimulaatioin. Työ koostuu kahdesta osasta: tähtipolymeerien synteesistä makrosyklisillä initiaattoreilla ja amfifiilisten tähtimäisten blokkikopolymeerien ominaisuuksien tutkimisesta.

An approach to palaeoseismicity in the Olkiluoto (sea) area during the early Holocene

Olkiluoto Island is situated in the northern Baltic Sea, near the southwestern coast of Finland, and is the proposed location of a spent nuclear fuel repository. This study examined Holocene palaeoseismicity in the Olkiluoto area and in the surrounding sea areas by computer simulations together with acoustic-seismic, sedimentological and dating methods. The most abundant rock type on the island is migmatic mica gneiss, intruded by tonalites, granodiorites and granites. The surrounding Baltic Sea seabed consists of Palaeoproterozoic crystalline bedrock, which is to a great extent covered by younger Mesoproterozoic sedimentary rocks. The area contains several ancient deep-seated fracture zones that divide it into bedrock blocks. The response of bedrock at the Olkiluoto site was modelled considering four future ice-age scenarios. Each scenario produced shear displacements of fractures with different times of occurrence and varying recovery rates. Generally, the larger the maximum ice load, the larger were the permanent shear displacements. For a basic case, the maximum shear displacements were a few centimetres at the proposed nuclear waste repository level, at proximately 500 m b.s.l. High-resolution, low-frequency echo-sounding was used to examine the Holocene submarine sedimentary structures and possible direct and indirect indicators of palaeoseismic activity in the northern Baltic Sea. Echo-sounding profiles of Holocene submarine sediments revealed slides and slumps, normal faults, debris flows and turbidite-type structures. The profiles also showed pockmarks and other structures related to gas or groundwater seepages, which might be related to fracture zone activation. Evidence of postglacial reactivation in the study area was derived from the spatial occurrence of some of the structures, especial the faults and the seepages, in the vicinity of some old bedrock fracture zones. Palaeoseismic event(s) (a single or several events) in the Olkiluoto area were dated and the palaeoenvironment was characterized using palaeomagnetic, biostratigraphical and lithostratigraphical methods, enhancing the reliability of the chronology. Combined lithostratigraphy, biostratigraphy and palaeomagnetic stratigraphy revealed an age estimation of 10 650 to 10 200 cal. years BP for the palaeoseismic event(s). All Holocene sediment faults in the northern Baltic Sea occur at the same stratigraphical level, the age of which is estimated at 10 700 cal. years BP (9500 radiocarbon years BP). Their movement is suggested to have been triggered by palaeoseismic event(s) when the Late Weichselian ice sheet was retreating from the site and bedrock stresses were released along the bedrock fracture zones. Since no younger or repeated traces of seismic events were found, it corroborates the suggestion that the major seismic activity occurred within a short time during and after the last deglaciation. The origin of the gas/groundwater seepages remains unclear. Their reflections in the echo-sounding profiles imply that part of the gas is derived from the organic-bearing Litorina and modern gyttja clays. However, at least some of the gas is derived from the bedrock. Additional information could be gained by pore water analysis from the pockmarks. Information on postglacial fault activation and possible gas and/or fluid discharges under high hydraulic heads has relevance in evaluating the safety assessment of a planned spent nuclear fuel repository in the region.

Affektiivisuuden laskennallinen määritteleminen

One of the most tangled fields of research is the field of defining and modeling affective concepts, i. e. concepts regarding emotions and feelings. The subject can be approached from many disciplines. The main problem is lack of generally approved definitions. However, e.g. linguists have recently started to check the consistency of their theories with the help of computer simulations. Definitions of affective concepts are needed for performing similar simulations in behavioral sciences. In this thesis, preliminary computational definitions of affects for a simple utility-maximizing agent are given. The definitions have been produced by synthetizing ideas from theories from several fields of research. The class of affects is defined as a superclass of emotions and feelings. Affect is defined as a process, in which a change in an agent's expected utility causes a bodily change. If the process is currently under the attention of the agent (i.e. the agent is conscious of it), the process is a feeling. If it is not, but can in principle be taken into attention (i.e. it is preconscious), the process is an emotion. Thus, affects do not presuppose consciousness, but emotions and affects do. Affects directed at unexpected materialized (i.e. past) events are delight and fright. Delight is the consequence of an unexpected positive event and fright is the consequence of an unexpected negative event. Affects directed at expected materialized (i.e. past) events are happiness (expected positive event materialized), disappointment (expected positive event did not materialize), sadness (expected negative event materialized) and relief (expected negative event did not materialize). Affects directed at expected unrealized (i.e. future) events are fear and hope. Some other affects can be defined as directed towards originators of the events. The affect classification has also been implemented as a computer program, the purpose of which is to ensure the coherence of the definitions and also to illustrate the capabilities of the model. The exact content of bodily changes associated with specific affects is not considered relevant from the point of view of the logical structure of affective phenomena. The utility function need also not be defined, since the target of examination is only its dynamics.

Theoretical Studies on Coupled Electron and Proton Transfer in Cytochrome c Oxidase

The complexity of life is based on an effective energy transduction machinery, which has evolved during the last 3.5 billion years. In aerobic life, the utilization of the high oxidizing potential of molecular oxygen powers this machinery. Oxygen is safely reduced by a membrane bound enzyme, cytochrome c oxidase (CcO), to produce an electrochemical proton gradient over the mitochondrial or bacterial membrane. This gradient is used for energy-requiring reactions such as synthesis of ATP by F0F1-ATPase and active transport. In this thesis, the molecular mechanism by which CcO couples the oxygen reduction chemistry to proton-pumping has been studied by theoretical computer simulations. By building both classical and quantum mechanical model systems based on the X-ray structure of CcO from Bos taurus, the dynamics and energetics of the system were studied in different intermediate states of the enzyme. As a result of this work, a mechanism was suggested by which CcO can prevent protons from leaking backwards in proton-pumping. The use and activation of two proton conducting channels were also enlightened together with a mechanism by which CcO sorts the chemical protons from pumped protons. The latter problem is referred to as the gating mechanism of CcO, and has remained a challenge in the bioenergetics field for more than three decades. Furthermore, a new method for deriving charge parameters for classical simulations of complex metalloenzymes was developed.

Mechanisms and molecular regulation of mammalian tooth replacement

In most non-mammalian vertebrates, such as fish and reptiles, teeth are replaced continuously. However, tooth replacement in most mammals, including human, takes place only once and further renewal is apparently inhibited. It is not known how tooth replacement is genetically regulated, and little is known on the physiological mechanism and evolutionary reduction of tooth replacement in mammals. In this study I have attempted to address these questions. In a rare human condition cleidocranial dysplasia, caused by a mutation in a Runt domain transcription factor Runx2, tooth replacement is continued. Runx2 mutant mice were used to investigate the molecular mechanisms of Runx2 function. Microarray analysis from dissected embryonic day 14 Runx2 mutant and wild type dental mesenchymes revealed many downstream targets of Runx2, which were validated using in situ hybridization and tissue culture methods. Wnt signaling inhibitor Dkk1 was identified as a candidate target, and in tissue culture conditions it was shown that Dkk1 is induced by FGF4 and this induction is Runx2 dependent. These experiments demonstrated a connection between Runx2, FGF and Wnt signaling in tooth development and possibly also in tooth replacement. The role of Wnt signaling in tooth replacement was further investigated by using a transgenic mouse model where Wnt signaling mediator β-catenin is continuously stabilized in dental epithelium. This stabilization led to activated Wnt signaling and to the formation of multiple enamel knots. In vitro and transplantation experiments were performed to examine the process of extra tooth formation. We showed that new teeth were continuously generated and that new teeth form from pre-existing teeth. A morphodynamic activator-inhibitor model was used to simulate enamel knot formation. By increasing the intrinsic production rate of the activator (β-catenin), the multiple enamel knot phenotype was reproduced by computer simulations. It was thus concluded that β-catenin acts as an upstream activator of enamel knots, closely linking Wnt signaling to the regulation of tooth renewal. As mice do not normally replace teeth, we used other model animals to investigate the physiological and genetic mechanisms of tooth replacement. Sorex araneus, the common shrew was earlier reported to have non-functional tooth replacement in all antemolar tooth positions. We showed by histological and gene expression studies that there is tooth replacement only in one position, the premolar 4 and that the deciduous tooth is diminished in size and disappears during embryogenesis without becoming functional. The growth rates of deciduous and permanent premolar 4 were measured and it was shown by competence inference that the early initiation of the replacement tooth in relation to the developmental stage of the deciduous tooth led to the inhibition of deciduous tooth morphogenesis. It was concluded that the evolutionary loss of deciduous teeth may involve the early activation of replacement teeth, which in turn suppress their predecessors. Mustela putorius furo, the ferret, has a dentition that resembles that of the human as ferrets have teeth that belong to all four tooth families, and all the antemolar teeth are replaced once. To investigate the replacement mechanism, histological serial sections from different embryonic stages were analyzed. It was noticed that tooth replacement is a process which involves the growth and detachment of the dental lamina from the lingual cervical loop of the deciduous tooth. Detachment of the deciduous tooth leads to a free successional dental lamina, which grows deeper into the mesenchyme, and later buds the replacement tooth. A careful 3D analysis of serial histological sections was performed and it was shown that replacement teeth are initiated from the successional dental lamina and not from the epithelium of the deciduous tooth. The molecular regulation of tooth replacement was studied and it was shown by examination of expression patterns of candidate regulatory genes that BMP/Wnt inhibitor Sostdc1 was strongly expressed in the buccal aspect of the dental lamina, and in the intersection between the detaching deciduous tooth and the successional dental lamina, suggesting a role for Sostdc1 in the process of detachment. Shh was expressed in the enamel knot and in the inner enamel epithelium in both generations of teeth supporting the view that the morphogenesis of both generations of teeth is regulated by similar mechanisms. In summary, histological and molecular studies on different model animals and transgenic mouse models were used to investigate tooth replacement. This thesis work has significantly contributed to the knowledge on the physiological mechanisms and molecular regulation of tooth replacement and its evolutionary suppression in mammals.

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.

Sterol-binding proteins in late endosomes : Regulation of endosome motility and lipid metabolism

Despite its bad reputation in the mass media, cholesterol is an indispensable constituent of cellular membranes and vertebrate life. It is, however, also potentially lethal as it may accumulate in the arterial intima causing atherosclerosis or elsewhere in the body due to inherited conditions. Studying cholesterol in cells, and research on how the cell biology of cholesterol affects on system level is essential for a better understanding of the disease states associated with cholesterol and for the development of new therapies for these conditions. On its way to the cell, exogenous cholesterol traverses through endosomes, transport vesicles involved in internalizing material to cells, and needs to be transported out of this compartment. This endosomal pool of cholesterol is important for understanding both the common disorders of metabolism and the more rare hereditary disorders of cholesterol metabolism. The study of cholesterol in cells has been hampered by the lack of bright fluorescent sterol analogs that would resemble cholesterol enough to be used in cellular studies. In the first study of my thesis, we present a new sterol analog, Boron-Dipyrromethene (BODIPY)-cholesterol for visualizing sterols in living cells and organism. This fluorescent cholesterol derivative is shown to behave similarly to cholesterol both by atomic scale computer simulations and biochemical experiments. We characterize its localization inside different types of living cells and show that it can be used to study sterol trafficking in living organisms. Two sterol binding proteins associated with the endosomal membrane; the Niemann-Pick type C disease protein 1 (NPC1) and the Oxysterol Binding Protein Related Protein 1 (ORP1) are the subjects of the rest of this study. Sensing cholesterol on endosomes, transporting lipids away from this compartment and the effects these lipids play on cellular metabolism are considered. In the second study we characterize how the NPC1 protein affects lipid metabolism. We show that this cholesterol binding protein affects synthesis of triglycerides and that genetic polymorphisms or a genetic defect in the NPC1 gene affect triglyceride on the whole body level. These effects take place via regulation of carbon fluxes to different lipid classes in cells. In the third part we characterize the effects of another endosomal sterol binding protein, ORP1L on the function and motility of endosomes. Specifically we elucidate how a mutation in the ability of ORP1L to bind sterols affects its behavior in cells, and how a change in ORP1L levels in cells affects the localization, degradative capacity and motility of endosomes. In addition we show that ORP1L manipulations affect cholesterol balance also in macrophages, a cell type important for the development of atherosclerosis.

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.

Towards Understanding Powder Behavior Via Simulation

The aim of this study was to investigate powder and tablet behavior at the level of mechanical interactions between single particles. Various aspects of powder packing, mixing, compression, and bond formation were examined with the aid of computer simulations. The packing and mixing simulations were based on spring forces interacting between particles. Packing and breakage simulations included systems in which permanent bonds were formed and broken between particles, based on their interaction strengths. During the process, a new simulation environment based on Newtonian mechanics and elementary interactions between the particles was created, and a new method for evaluating mixing was developed. Powder behavior is a complicated process, and many of its aspects are still unclear. Powders as a whole exhibit some aspects of solids and others of liquids. Therefore, their physics is far from clear. However, using relatively simple models based on particle-particle interaction, many powder properties could be replicated during this work. Simulated packing densities were similar to values reported in the literature. The method developed for describing powder mixing correlated well with previous methods. The new method can be applied to determine mixing in completely homogeneous materials, without dividing them into different components. As such, it can describe the efficiency of the mixing method, regardless of the powder's initial setup. The mixing efficiency at different vibrations was examined, and we found that certain combinations of amplitude, direction, and frequencies resulted in better mixing while using less energy. Simulations using exponential force potentials between particles were able to explain the elementary compression behavior of tablets, and create force distributions that were similar to the pressure distributions reported in the literature. Tablet-breaking simulations resulted in breaking strengths that were similar to measured tablet breaking strengths. In general, many aspects of powder behavior can be explained with mechanical interactions at the particle level, and single particle properties can be reliably linked to powder behavior with accurate simulations.

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.