Sticking and erosion at carbon-containing plasma-facing materials in fusion reactors

Controlled nuclear fusion is one of the most promising sources of energy for the future. Before this goal can be achieved, one must be able to control the enormous energy densities which are present in the core plasma in a fusion reactor. In order to be able to predict the evolution and thereby the lifetime of different plasma facing materials under reactor-relevant conditions, the interaction of atoms and molecules with plasma first wall surfaces have to be studied in detail. In this thesis, the fundamental sticking and erosion processes of carbon-based materials, the nature of hydrocarbon species released from plasma-facing surfaces, and the evolution of the components under cumulative bombardment by atoms and molecules have been investigated by means of molecular dynamics simulations using both analytic potentials and a semi-empirical tight-binding method. The sticking cross-section of CH3 radicals at unsaturated carbon sites at diamond (111) surfaces is observed to decrease with increasing angle of incidence, a dependence which can be described by a simple geometrical model. The simulations furthermore show the sticking cross-section of CH3 radicals to be strongly dependent on the local neighborhood of the unsaturated carbon site. The erosion of amorphous hydrogenated carbon surfaces by helium, neon, and argon ions in combination with hydrogen at energies ranging from 2 to 10 eV is studied using both non-cumulative and cumulative bombardment simulations. The results show no significant differences between sputtering yields obtained from bombardment simulations with different noble gas ions. The final simulation cells from the 5 and 10 eV ion bombardment simulations, however, show marked differences in surface morphology. In further simulations the behavior of amorphous hydrogenated carbon surfaces under bombardment with D^+, D^+2, and D^+3 ions in the energy range from 2 to 30 eV has been investigated. The total chemical sputtering yields indicate that molecular projectiles lead to larger sputtering yields than atomic projectiles. Finally, the effect of hydrogen ion bombardment of both crystalline and amorphous tungsten carbide surfaces is studied. Prolonged bombardment is found to lead to the formation of an amorphous tungsten carbide layer, regardless of the initial structure of the sample. In agreement with experiment, preferential sputtering of carbon is observed in both the cumulative and non-cumulative simulations

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.

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.

Morphogenesis of Mammalian Endoplasmic Reticulum and Golgi Apparatus throughout the Cell Cycle

The endoplasmic reticulum (ER) and the Golgi apparatus are organelles that produce, modify and transport proteins and lipids and regulate Ca2+ environment within cells. Structurally they are composed of sheets and tubules. Sheets may take various forms: intact, fenestrated, single or stacked. The ER, including the nuclear envelope, is a single continuous network, while the Golgi shows only some level of connectivity. It is often unclear, how different morphologies correspond to particular functions. Previous studies indicate that the structures of the ER and Golgi are dynamic and regulated by fusion and fission events, cytoskeleton, rate of protein synthesis and secretion, and specific structural proteins. For example, many structural proteins shaping tubular ER have been identified, but sheet formation is much more unclear. In this study, we used light and electron microscopy to study morphological changes of the ER and Golgi in mammalian cells. The proportion, type, location and dynamics of ER sheets and tubules were found to vary in a cell type or cell cycle stage dependent manner. During interphase, ER and Golgi structures were demonstrated to be regulated by p37, a cofactor of the fusion factor p97, and microtubules, which also affected the localization of the organelles. Like previously shown for the Golgi, the ER displayed a tendency for fenestration and tubulation during mitosis. However, this shape change did not result in ER fragmentation as happens to Golgi, but a continuous network was retained. The activity of p97/p37 was found to be important for the reassembly of both organelles after mitosis. In EM images, ER sheet membranes appear rough, since they contain attached ribosomes, whereas tubular membranes appear smooth. Our studies revealed that structural changes of the ER towards fenestrated and tubular direction correlate with loss of ER-bound ribosomes and vice versa. High and low curvature ER membranes have a low and high density of ribosomes, respectively. To conclude, both ER and Golgi architecture depend on fusion activity of p97/p37. ER morphogenesis, particularly of the sheet shape, is intimately linked to the density of membrane bound ribosomes.

Image fusion in GRDSS for Land use Mapping

Image fusion techniques are useful to integrate the geometric detail of a high-resolution panchromatic (PAN) image and the spectral information of a low-resolution multispectral (MSS) image, particularly important for understanding land use dynamics at larger scale (1:25000 or lower), which is required by the decision makers to adopt holistic approaches for regional planning. Fused images can extract features from source images and provide more information than one scene of MSS image. High spectral resolution aids in identification of objects more distinctly while high spatial resolution allows locating the objects more clearly. The geoinformatics technologies with an ability to provide high-spatial-spectral-resolution data helps in inventorying, mapping, monitoring and sustainable management of natural resources. Fusion module in GRDSS, taking into consideration the limitations in spatial resolution of MSS data and spectral resolution of PAN data, provide high-spatial-spectral-resolution remote sensing images required for land use mapping on regional scale. GRDSS is a freeware GIS Graphic User Interface (GUI) developed in Tcl/Tk is based on command line arguments of GRASS (Geographic Resources Analysis Support System) with the functionalities for raster analysis, vector analysis, site analysis, image processing, modeling and graphics visualization. It has the capabilities to capture, store, process, analyse, prioritize and display spatial and temporal data.

Drosophila Erect wing (Ewg) controls mitochondrial fusion during muscle growth and maintenance by regulation of the Opa1-like gene

Mitochondrial biogenesis and morphological changes are associated with tissue-specific functional demand, but the factors and pathways that regulate these processes have not been completely identified. A lack of mitochondrial fusion has been implicated in various developmental and pathological defects. The spatiotemporal regulation of mitochondrial fusion in a tissue such as muscle is not well understood. Here, we show in Drosophila indirect flight muscles (IFMs) that the nuclear-encoded mitochondrial inner membrane fusion gene, Opa1-like, is regulated in a spatiotemporal fashion by the transcription factor/co-activator Erect wing (Ewg). In IFMs null for Ewg, mitochondria undergo mitophagy and/or autophagy accompanied by reduced mitochondrial functioning and muscle degeneration. By following the dynamics of mitochondrial growth and shape in IFMs, we found that mitochondria grow extensively and fuse during late pupal development to form the large tubular mitochondria. Our evidence shows that Ewg expression during early IFM development is sufficient to upregulate Opa1-like, which itself is a requisite for both late pupal mitochondrial fusion and muscle maintenance. Concomitantly, by knocking down Opa1-like during early muscle development, we show that it is important for mitochondrial fusion, muscle differentiation and muscle organization. However, knocking down Opa1-like, after the expression window of Ewg did not cause mitochondrial or muscle defects. This study identifies a mechanism by which mitochondrial fusion is regulated spatiotemporally by Ewg through Opa1-like during IFM differentiation and growth.

Speaker verification based on the fusion of speech acoustics and inverted articulatory signals

We propose apractical, feature-level and score-level fusion approach by combining acoustic and estimated articulatory information for both text independent and text dependent speaker verification. From a practical point of view, we study how to improve speaker verification performance by combining dynamic articulatory information with the conventional acoustic features. On text independent speaker verification, we find that concatenating articulatory features obtained from measured speech production data with conventional Mel-frequency cepstral coefficients (MFCCs) improves the performance dramatically. However, since directly measuring articulatory data is not feasible in many real world applications, we also experiment with estimated articulatory features obtained through acoustic-to-articulatory inversion. We explore both feature level and score level fusion methods and find that the overall system performance is significantly enhanced even with estimated articulatory features. Such a performance boost could be due to the inter-speaker variation information embedded in the estimated articulatory features. Since the dynamics of articulation contain important information, we included inverted articulatory trajectories in text dependent speaker verification. We demonstrate that the articulatory constraints introduced by inverted articulatory features help to reject wrong password trials and improve the performance after score level fusion. We evaluate the proposed methods on the X-ray Microbeam database and the RSR 2015 database, respectively, for the aforementioned two tasks. Experimental results show that we achieve more than 15% relative equal error rate reduction for both speaker verification tasks. (C) 2015 Elsevier Ltd. All rights reserved.

Membrane Deformability and Membrane Tension of Single Isolated Mitochondria

Mitochondria dynamics is crucial to many biological processes such as mitochondria fusion and fission, which is highly correlated to the mechanics of single mitochondria. However, the mechanobiological coupling of mitochondria has been poorly understood. Here membrane deformability and membrane tension of individual mitochondria isolated from MtDsRed labeled human embryonic T-Rex-293 kidney cells were measured using a micropipette aspiration assay. The results demonstrated that membrane deformation of isolated mitochondria exhibited an elastic transition phase followed by an equilibrium phase, and mitochondrial membrane tension was proportional to the area compressibility. It was also indicated that mitochondrial membrane deformability was significantly affected by physical chemical factors such as osmotic pressure or pH value, and was further correlated to mitochondrial functionality in different respiratory states and Ca2+ regulation. These findings provide a new insight into understanding the mechanical regulation of mitochondrial physiology.

Two overrun phenomena and their effects on fusion yield in the Coulomb explosion of heteronuclear clusters

Two overrun effects in the Coulomb explosion dynamics of heteronuclear clusters have been investigated theoretically by the use of a simplified electrostatic model. When the charge-to-mass ratio of light ions is higher than that of heavy ions, the light ions can overtake the heavy ions inside the cluster and acquire a higher kinetic energy. Further, if the charge density of the heavy ions is twice as high as that of the light ions, i.e. a proposed competitive parameter xi = rho BqB/rho AqA > 2, the inner light ions can overtake those light ions on the surface of the cluster and form a shock shell during the explosion, which might drive the intracluster collision and fusion of the light ions. Different regimes of nuclear fusion are discussed and the corresponding neutron yields are estimated. Our analysis indicates that the probability of intracluster fusion is quite low even if deuterated heteronuclear clusters such as (DI)(n) with large size and high competitive parameter are employed. However, heteronuclear clusters are still a better candidate compared with homonuclear clusters for enhancing the total intercluster fusion yield because both a higher energy region and a higher proportion of deuterons distributing in the energy region can be created in the deuterated heteronuclear clusters.

Battery-Powered RF Pre-Ionization System for the Caltech Magnetohydrodynamically-Driven Jet Experiment: RF Discharge Properties and MHD-Driven Jet Dynamics

This thesis describes investigations of two classes of laboratory plasmas with rather different properties: partially ionized low pressure radiofrequency (RF) discharges, and fully ionized high density magnetohydrodynamically (MHD)-driven jets. An RF pre-ionization system was developed to enable neutral gas breakdown at lower pressures and create hotter, faster jets in the Caltech MHD-Driven Jet Experiment. The RF plasma source used a custom pulsed 3 kW 13.56 MHz RF power amplifier that was powered by AA batteries, allowing it to safely float at 4-6 kV with the cathode of the jet experiment. The argon RF discharge equilibrium and transport properties were analyzed, and novel jet dynamics were observed.

Although the RF plasma source was conceived as a wave-heated helicon source, scaling measurements and numerical modeling showed that inductive coupling was the dominant energy input mechanism. A one-dimensional time-dependent fluid model was developed to quantitatively explain the expansion of the pre-ionized plasma into the jet experiment chamber. The plasma transitioned from an ionizing phase with depressed neutral emission to a recombining phase with enhanced emission during the course of the experiment, causing fast camera images to be a poor indicator of the density distribution. Under certain conditions, the total visible and infrared brightness and the downstream ion density both increased after the RF power was turned off. The time-dependent emission patterns were used for an indirect measurement of the neutral gas pressure.

The low-mass jets formed with the aid of the pre-ionization system were extremely narrow and collimated near the electrodes, with peak density exceeding that of jets created without pre-ionization. The initial neutral gas distribution prior to plasma breakdown was found to be critical in determining the ultimate jet structure. The visible radius of the dense central jet column was several times narrower than the axial current channel radius, suggesting that the outer portion of the jet must have been force free, with the current parallel to the magnetic field. The studies of non-equilibrium flows and plasma self-organization being carried out at Caltech are relevant to astrophysical jets and fusion energy research.

Duplication-degeneration as a mechanism of gene fission and the origin of new genes in Drosophila species

Gene fission and fusion, the processes by which a single gene is split into two separate genes and two adjacent genes are fused into a single gene, respectively, are among the primary processes that generate new genes(1-4). Despite their seeming reversibi

Characterize dynamic conformational space of human CCR5 extracellular domain by molecular modeling and molecular dynamics simulation

The chemokine receptor CCR5 is the receptor for several chemokines and major coreceptor for R5 human immunodeficiency virus type-1 strains entry into cell. Three-dimensional models of CCR5 were built by using homology modeling approach and 1 ns molecular dynamics (MD) simulation, because studies of site-directed mutagenesis and chimeric receptors have indicated that the N-terminus (Nt) and extracellular loops (ECLs) of CCR5 are important for ligands binding and viral fusion and entry, special attention was focused on disulfide bond function, conformational flexibility, hydrogen bonding, electrostatic interactions, and solvent-accessible surface area of Nt and ECLs of this protein part. We found that the extracellular segments of CCR5 formed a well-packet globular domain with complex interactions occurred between them in a majority of time of MID simulation, but Nt region could protrude from this domain sometimes. The disulfide bond Cys20-Cys269 is essential in controlling specific orientation of Nt region and maintaining conformational integrity of extracellular domain. RMS comparison analysis between conformers revealed the ECL1 of CCR5 stays relative rigid, whereas the ECL2 and Nt are rather flexible. Solvent-accessible surface area calculations indicated that the charged residues within Nt and ECL2 are often exposed to solvent. Integrating these results with available experimental data, a two-step gp120-CCR5 binding mechanism was proposed. The dynamic interaction of CCR5 extracellular domain with gp120 was emphasized. (C) 2004 Elsevier B.V. All rights reserved.

Controlling electronic structures by irradiation in single-walled SiC nanotubes: a first-principles molecular dynamics study

Using first-principles molecular dynamics simulations, the displacement threshold energy and defect configurations are determined in SiC nanotubes. The simulation results reveal that a rich variety of defect structures (vacancies, Stone-Wales defects and antisite defects) are formed with threshold energies from 11 to 64 eV. The threshold energy shows an anisotropic behavior and exhibits a dramatic decrease with decreasing tube diameter. The electronic structure can be altered by the defects formed by irradiation, which suggests that the electron irradiation may be a way to use defect engineering to tailor electronic properties of SiC nanotubes.

Bubble Dynamics and Heat Transfer in Microgravity Pool Boiling

Boiling is an extremely complicated and illusive process. Microgravity experiments offer a unique opportunity to study the complex interactions without external forces, such as buoyancy, which can affect the bubble dynamics and the related heat transfer. Furthermore, they can also provide a means to study the actual influence of gravity on the boiling. Two research projects on pool boiling in microgravity have been conducted aboard the Chinese recoverable satellites. Ground-based experiments both in normal gravity and in short-term microgravity in the Drop Tower Beijing and numerical simulations have also been performed. Steady boiling of R113 on thin platinum wires was studied with a temperature-controlled heating method, while quasi-steady boiling of FC-72 on a plane plate was investigated with an exponentially increasing heating voltage. It was found that the bubble dynamics in microgravity has a distinct difference from that in normal gravity, and that the heat transfer characteristic is depended upon the bubble dynamics. Lateral motions of bubbles on the heaters were observed before their departure in microgravity. The surface oscillation of the merged bubbles due to lateral coalescence between adjacent bubbles drove it to detach from the heaters. Slight enhancement of heat transfer on wires is observed in microgravity, while diminution is evident for high heat flux in the plate case.