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.

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.

Effect of phenolic-rich plant materials on protein and lipid oxidation reactions

The antioxidant activity of natural plant materials rich in phenolic compounds is being widely investigated for protection of food products sensitive to oxidative reactions. In this thesis plant materials rich in phenolic compounds were studied as possible antioxidants to prevent protein and lipid oxidation reactions in different food matrixes such as pork meat patties and corn oil-in water emulsions. Loss of anthocyanins was also measured during oxidation in corn oil-in-water emulsions. In addition, the impact of plant phenolics on amino acid level was studied using tryptophan as a model compound to elucidate their role in preventing the formation of tryptophan oxidation products. A high-performance liquid chromatography (HPLC) method with ultraviolet and fluorescence detection (UV-FL) was developed that enabled fast investigation of formation of tryptophan derived oxidation products. Byproducts of oilseed processes such as rapeseed (Brassica rapa L.), camelina (Camelina sativa) and soy meal (Glycine max L.) as well as Scots pine bark (Pinus sylvestris) and several reference compounds were shown to act as antioxidants toward both protein and lipid oxidation in cooked pork meat patties. In meat, the antioxidant activity of camelina, rapeseed and soy meal were more pronounced when used in combination with a commercial rosemary extract (Rosmarinus officinalis). Berry phenolics such as black currant (Ribes nigrum) anthocyanins and raspberry (Rubus idaeus) ellagitannins showed potent antioxidant activity in corn oil-in-water emulsions toward lipid oxidation with and without β-lactoglobulin. The antioxidant effect was more pronounced in the presence of β-lactoglobulin. The berry phenolics also inhibited the oxidation of tryptophan and cysteine side chains of β-lactoglobulin. The results show that the amino acid side chains were oxidized prior the propagation of lipid oxidation, thereby inhibiting fatty acid scission. In addition, the concentration and color of black currant anthocyanins decreased during the oxidation. Oxidation of tryptophan was investigated in two different oxidation models with hydrogen peroxide (H2O2) and hexanal/FeCl2. Oxidation of tryptophan in both models resulted in oxidation products such as 3a-hydroxypyrroloindole-2-carboxylic acid, dioxindolylalanine, 5-hydroxy-tryptophan, kynurenine, N-formylkynurenine and β-oxindolylalanine. However, formation of tryptamine was only observed in tryptophan oxidized in the presence of H2O2. Pine bark phenolics, black currant anthocyanins, camelina meal phenolics as well as cranberry proanthocyanidins (Vaccinium oxycoccus) provided the best antioxidant effect toward tryptophan and its oxidation products when oxidized with H2O2. The tryptophan modifications formed upon hexanal/FeCl2 treatment were efficiently inhibited by camelina meal followed by rapeseed and soy meal. In contrast, phenolics from raspberry, black currant, and rowanberry (Sorbus aucuparia) acted as weak prooxidants. This thesis contributes to elucidating the effects of natural phenolic compounds as potential antioxidants in order to control and prevent protein and lipid oxidation reactions. Understanding the relationship between phenolic compounds and proteins as well as lipids could lead to the development of new, effective, and multifunctional antioxidant strategies that could be used in food, cosmetic and pharmaceutical applications.

Regio- and stereoselectivity of oxidative coupling reactions of phenols : Spirodienones as construction units in lignin

Dimeric phenolic compounds lignans and dilignols form in the so-called oxidative coupling reaction of phenols. Enzymes such as peroxidases and lac-cases catalyze the reaction using hydrogen peroxide or oxygen respectively as oxidant generating phenoxy radicals which couple together according to certain rules. In this thesis, the effects of the structures of starting materials mono-lignols and the effects of reaction conditions such as pH and solvent system on this coupling mechanism and on its regio- and stereoselectivity have been studied. After the primary coupling of two phenoxy radicals a very reactive quinone me-thide intermediate is formed. This intermediate reacts quickly with a suitable nucleophile which can be, for example, an intramolecular hydroxyl group or another nucleophile such as water, methanol, or a phenolic compound in the reaction system. This reaction is catalyzed by acids. After the nucleophilic addi-tion to the quinone methide, other hydrolytic reactions, rearrangements, and elimination reactions occur leading finally to stable dimeric structures called lignans or dilignols. Similar reactions occur also in the so-called lignification process when monolignol (or dilignol) reacts with the growing lignin polymer. New kinds of structures have been observed in this thesis. The dimeric com-pounds with so-called spirodienone structure have been observed to form both in the dehydrodimerization of methyl sinapate and in the beta-1-type cross-coupling reaction of two different monolignols. This beta-1-type dilignol with a spirodienone structure was the first synthetized and published dilignol model compound, and at present, it has been observed to exist as a fundamental construction unit in lignins. The enantioselectivity of the oxidative coupling reaction was also studied for obtaining enantiopure lignans and dilignols. A rather good enantioselectivity was obtained in the oxidative coupling reaction of two monolignols with chiral auxiliary substituents using peroxidase/H2O2 as an oxidation system. This observation was published as one of the first enantioselective oxidative coupling reaction of phenols. Pure enantiomers of lignans were also obtained by using chiral cryogenic chromatography as a chiral resolution technique. This technique was shown to be an alternative route to prepare enantiopure lignans or lignin model compounds in a preparative scale.

The kinetics and reactivity of several polyatomic free radicals in reactions with NO2, O2, Cl2, and HCl

In this thesis, the kinetics of several alkyl, halogenated alkyl, and alkenyl free radical reactions with NO2, O2, Cl2, and HCl reactants were studied over a wide temperature range in time resolved conditions. Laser photolysis photoionisation mass spectrometer coupled to a flow reactor was the experimental method employed and this thesis present the first measurements performed with the experimental system constructed. During this thesis a great amount of work was devoted to the designing, building, testing, and improving the experimental apparatus. Carbon-centred free radicals were generated by the pulsed 193 or 248 nm photolysis of suitable precursors along the tubular reactor. The kinetics was studied under pseudo-first-order conditions using either He or N2 buffer gas. The temperature and pressure ranges employed were between 190 and 500 K, and 0.5 45 torr, respectively. The possible role of heterogeneous wall reactions was investigated employing reactor tubes with different sizes, i.e. to significantly vary the surface to volume ratio. In this thesis, significant new contributions to the kinetics of carbon-centred free radical reactions with nitrogen dioxide were obtained. Altogether eight substituted alkyl (CH2Cl, CHCl2, CCl3, CH2I, CH2Br, CHBr2, CHBrCl, and CHBrCH3) and two alkenyl (C2H3, C3H3) free radical reactions with NO2 were investigated as a function of temperature. The bimolecular rate coefficients of all these reactions were observed to possess negative temperature dependencies, while pressure dependencies were not noticed for any of these reactions. Halogen substitution was observed to moderately reduce the reactivity of substituted alkyl radicals in the reaction with NO2, while the resonance stabilisation of the alkenyl radical lowers its reactivity with respect to NO2 only slightly. Two reactions relevant to atmospheric chemistry, CH2Br + O2 and CH2I + O2, were also investigated. It was noticed that while CH2Br + O2 reaction shows pronounced pressure dependence, characteristic of peroxy radical formation, no such dependence was observed for the CH2I + O2 reaction. Observed primary products of the CH2I + O2 reaction were the I-atom and the IO radical. Kinetics of CH3 + HCl, CD3 + HCl, CH3 + DCl, and CD3 + DCl reactions were also studied. While all these reactions possess positive activation energies, in contrast to the other systems investigated in this thesis, the CH3 + HCl and CD3 + HCl reactions show a non-linear temperature dependency on the Arrhenius plot. The reactivity of substituted methyl radicals toward NO2 was observed to increase with decreasing electron affinity of the radical. The same trend was observed for the reactions of substituted methyl radicals with Cl2. It is proposed that interactions of frontier orbitals are responsible to these observations and Frontier Orbital Theory could be used to explain the observed reactivity trends of these highly exothermic reactions having reactant-like transition states.

Lipid assisted membrane entry in specific phospholipid targeted protein systems

The purpose of this thesis project is to study changes in the physical state of cell membranes during cell entry, including how these changes are connected to the presence of ceramide. The role of enzymatical manipulation of lipids in bacterial internalization is also studied. A novel technique, where a single giant vesicle is chosen under the microscope and an enzyme coupled-particle attached to the micromanipulator pipette towards the vesicle, is used. Thus, the enzymatic reaction on the membrane of the giant vesicle can be followed in real-time. The first aim of this study is to develop a system where the localized sphingomyelinase membrane interaction could be observed on the surface of the giant vesicle and the effects could be monitored with microscopy. Domain formation, which resembles acid sphingomyelinase (ASMase), causes CD95 clustering in the cell membrane due to ceramide production (Grassmé et al., 2001a; Grassmé et al., 2001b) and the formation of small vesicles inside the manipulated giant vesicle is observed. Sphingomyelinase activation has also been found to be an important factor in the bacterial and viral invasion process in nonphagocytic cells (Grassmé et al., 1997; Jan et al., 2000). Accordingly, sphingomyelinase reactions in the cell membrane might also give insight into bacterial or viral cellular entry events. We found sphingomyelinase activity in Chlamydia pneumonia elementarybodies (EBs). Interestingly, the bacterium enters host cells by endocytosis but the internalization mechanism of Chlamydia is unknown. The hypothesis is that sphingomyelin is needed for host cell entry in the infection of C. pneumonia. The second project focuses on this subject. The goal of the third project is to study a role of phosphatidylserine as a target for a membrane binding protein. Phosphatidylserine is chosen because of its importance in fusion processes. This will be another example for the importance of lipids in cell targeting, internalization, and externalization.

Proton translocation coupled to electron transfer reactions in terminal oxidases

Terminal oxidases are the final proteins of the respiratory chain in eukaryotes and some bacteria. They catalyze most of the biological oxygen consumption on Earth done by aerobic organisms. During the catalytic reaction terminal oxidases reduce dioxygen to water and use the energy released in this process to maintain the electrochemical proton gradient by functioning as a redox-driven proton pump. This membrane gradient of protons is extremely important for cells as it is used for many cellular processes, such as transportation of substrates and ATP synthesis. Even though the structures of several terminal oxidases are known, they are not sufficient in themselves to explain the molecular mechanism of proton pumping. In this work we have applied a complex approach using a variety of different techniques to address the properties and the mechanism of proton translocation by the terminal oxidases. The combination of direct measurements of pH changes during catalytic turnover, time-resolved potentiometric electrometry and optical spectroscopy, made it possible to obtain valuable information about various aspects of oxidase functioning. We compared oxygen binding properties of terminal oxidases from the distinct heme-copper (CcO) and cytochrome bd families and found that cytochrome bd has a high affinity for oxygen, which is 3 orders of magnitude higher than that of CcO. Interestingly, the difference between CcO and cytochrome bd is not only in higher affinity of the latter to oxygen, but also in the way that each of these enzymes traps oxygen during catalysis. CcO traps oxygen kinetically - the molecule of bound dioxygen is rapidly reduced before it can dissociate. Alternatively, cytochrome bd employs an alternative mechanism of oxygen trapping - part of the redox energy is invested into tight oxygen binding, and the price paid for this is the lack of proton pumping. A single cycle of oxygen reduction to water is characterized by translocation of four protons across the membrane. Our results make it possible to assign the pumping steps to discrete transitions of the catalytic cycle and indicate that during in vivo turnover of the oxidase these four protons are transferred, one at a time, during the P→F, F→OH, Oh→Eh, and Eh→R transitions. At the same time, each individual proton translocation step in the catalytic cycle is not just a single reaction catalyzed by CcO, but rather a complicated sequence of interdependent electron and proton transfers. We assume that each single proton translocation cycle of CcO is assured by internal proton transfer from the conserved Glu-278 to an as yet unidentified pump site above the hemes. Delivery of a proton to the pump site serves as a driving reaction that forces the proton translocation cycle to continue.

Transforming growth factor -beta signaling through Smad3 in allergy : Studies in the mechanisms of asthma, atopic dermatitis and allergic contact reactions

Transforming growth factor β signalling through Smad3 in allergy Allergic diseases, such as atopic dermatitis, asthma, and contact dermatitis are complex diseases influenced by both genetic and environmental factors. It is still unclear why allergy and subsequent allergic disease occur in some individuals but not in others. Transforming growth factor (TGF)-β is an important immunomodulatory and fibrogenic factor that regulates cellular processes in injured and inflamed skin. TGF-β has a significant role in the regulation of the allergen-induced immune response participating in the development of allergic and asthmatic inflammation. TGF-β is known to be an immunomodulatory factor in the progression of delayed type hypersensitivity reactions and allergic contact dermatitis. TGF-β is crucial in regulating the cellular responses involved in allergy, such as differentiation, proliferation and migration. TGF-β signals are delivered from the cytoplasm to the nucleus by TGF-β signal transducers called Smads. Smad3 is a major signal transducer in TGF-β -signalling that controls the expression of target genes in the nucleus in a cell-type specific manner. The role of TGF-β-Smad3 -signalling in the immunoregulation and pathophysiology of allergic disorders is still poorly understood. In this thesis, the role of TGF-β-Smad -signalling pathway using Smad3 -deficient knock out mice in the murine models of allergic diseases; atopic dermatitis, asthma and allergic contact reactions, was examined. Smad3-pathway regulates allergen induced skin inflammation and systemic IgE antibody production in a murine model atopic dermatitis. The defect in Smad3 -signalling decreased Th2 cytokine (IL-13 and IL-5) mRNA expression in the lung, modulated allergen induced specific IgG1 response, and affected mucus production in the lung in a murine model of asthma. TGF-β / Smad3 -signalling contributed to inflammatory hypersensitivity reactions and disease progression via modulation of chemokine and cytokine expression and inflammatory cell recruitment, cell proliferation and regulation of the specific antibody response in a murine model of contact hypersensitivity. TGF-β modulates inflammatory responses - at least partly through the Smad3 pathway - but also through other compensatory, non-Smad-dependent pathways. Understanding the effects of the TGF-β signalling pathway in the immune system and in disease models can help in elucidating the multilevel effects of TGF-β. Unravelling the mechanisms of Smad3 may open new possibilities for treating and preventing allergic responses, which may lead to severe illness and loss of work ability. In the future the Smad3 signalling pathway might be a potential target in the therapy of allergic diseases.