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.

The diagnostics of density distribution for inhomogeneous dense DT plasmas using fast protons

The density distribution of inhomogeneous dense deuterium-tritium plasmas in laser fusion is revealed by the energy loss of fast protons going through the plasma. In our simulation of a plasma density diagnostics, the fast protons used for the diagnostics may be generated in the laser-plasma interaction. Dividing a two-dimensional area into grids and knowing the initial and final energies of the protons, we can obtain a large linear and ill-posed equation set. for the densities of all grids, which is solved with the Tikhonov regularization method. We find that the accuracy of the set plan with four proton sources is better than those of the set plans with less than four proton sources. Also we have done the density reconstruction especially. for four proton sources with and without assuming circularly symmetrical density distribution, and find that the accuracy is better for the reconstruction assuming circular symmetry. The error is about 9% when no noise is added to the final energy for the reconstruction of four proton sources assuming circular symmetry. The accuracies for different random noises to final proton energies with four proton sources are also calculated.

Laser-confined fusion

An approach for producing a large quantity of neutrons is proposed. It involves compression of a fuel foil and confinement of the resulting plasma between two intense laser pulses. It is shown that two circularly polarized laser pulses of amplitude a=7 illuminating a deuterium-tritium foil of areal density 3.3 X 10(18) cm(-2) can produce about 4.2 X 10(6) neutrons per joule of the input laser energy.

Neutron source from thin foil confined by two laser pulses

Neutron production from a thin deuterium-tritium (D-T) foil irradiated by two intense femtosecond laser pulses from opposite sides with zero phase difference is studied analytically and numerically. For the interaction of a laser pulse of amplitude a = 7, focal area 100 mu m(2) and areal density 4.4 x 10(18) cm(-2) with a D-T plasma foil, about 1.17 x 10(21) neutron s(-1) can be obtained, much more than from other methods. The profiles of the ion and electron densities are also calculated.

The propagation and deposition of fast muons in dense DT mixture

The propagation of the fast muon population mainly due to collisional effect in a dense deuterium-tritium (DT for short) mixture is investigated and analysed within the framework of the relativistic Fokker-Planck equation. Without the approximation that the muons propagate straightly in the DT mixture, the muon penetration length, the straggling length, and the mean transverse dispersion radius are calculated for different initial energies, and especially for different densities of the densely compressed DT mixture in our suggested muon-driven fast ignition (FI). Unlike laser-driven FI requiring super-high temperature, muons can catalyze DT fusion at lower temperatures and may generate an ignition sparkle before the self-heating fusion follows. Our calculation is important for the feasibility and the experimental study of muon-driven FI.

Bubble regime for ion acceleration in a laser-driven plasma

Proton trapping and acceleration by an electron bubble-channel structure in laser interaction with high-density plasma is investigated by using three-dimensional particle-in-cell simulations. It is shown that protons can be trapped, bunched, and efficiently accelerated for appropriate laser and plasma parameters, and the proton acceleration is enhanced if the plasma consists mainly of heavier ions such as tritium. The observed results are analyzed and discussed in terms of a one-dimensional analytical three-component-plasma wake model.

Observations of hydrogen and helium isotopes in solar cosmic rays

The isotopic composition of hydrogen and helium in solar cosmic rays provides a means of studying solar flare particle acceleration mechanisms since the enhanced relative abundance of rare isotopes, such as ²H, ³H and ³He, is due to their production by inelastic nuclear collisions in the solar atmosphere during the flare. In this work the Caltech Electron/Isotope Spectrometer on the IMP-7 spacecraft has been used to measure this isotopic composition. The response of the dE/dx-E particle telescope is discussed and alpha particle channeling in thin detectors is identified as an important background source affecting measurement of low values of (³He/⁴He).

The following flare-averaged results are obtained for the period, October, 1972 - November, 1973: (²H/¹H) = 7⁺¹⁰_-6 X 10^-6 (1.6 - 8.6 MeV/nuc), (³H/¹H) less than 3.4 x 10^-6 (1.2 - 6.8 MeV/nuc), (³He/⁴He) = (9 ± 4) x 10^-3, (³He/¹H) = (1.7 ± 0.7) x 10^-4 (3.1 - 15.0 MeV/nuc). The deuterium and tritium ratios are significantly lower than the same ratios at higher energies, suggesting that the deuterium and tritium spectra are harder than that of the protons. They are, however, consistent with the same thin target model relativistic path length of ~ 1 g/cm² (or equivalently ~ 0.3 g/cm² at 30 MeV/nuc) which is implied by the higher energy results. The ³He results, consistent with previous observations, would imply a path length at least 3 times as long, but the observations may be contaminated by small ³He rich solar events.

During 1973 three "³He rich events," containing much more ³He than ²H or ³H were observed on 14 February, 29 June and 5 September. Although the total production cross sections for ²H,³H and ³He are comparable, an upper limit to (²H/³He) and (³H/³He) was 0.053 (2.9-6.8 MeV/nuc), summing over the three events. This upper limit is marginally consistent with Ramaty and Kozlovsky's thick target model which accounts for such events by the nuclear reaction kinematics and directional properties of the flare acceleration process. The 5 September event was particularly significant in that much more ³He was observed than ⁴He and the fluxes of ³He and ¹H were about equal. The range of (³He/⁴He) for such events reported to date is 0.2 to ~ 6 while (³He/¹H) extends from 10^-3 to ~ 1. The role of backscattered and mirroring protons and alphas in accounting for such variations is discussed.

DNA-Mediated Charge Transport Devices for Protein Detection

Detection of biologically relevant targets, including small molecules, proteins, DNA, and RNA, is vital for fundamental research as well as clinical diagnostics. Sensors with biological elements provide a natural foundation for such devices because of the inherent recognition capabilities of biomolecules. Electrochemical DNA platforms are simple, sensitive, and do not require complex target labeling or expensive instrumentation. Sensitivity and specificity are added to DNA electrochemical platforms when the physical properties of DNA are harnessed. The inherent structure of DNA, with its stacked core of aromatic bases, enables DNA to act as a wire via DNA-mediated charge transport (DNA CT). DNA CT is not only robust over long molecular distances of at least 34 nm, but is also especially sensitive to anything that perturbs proper base stacking, including DNA mismatches, lesions, or DNA-binding proteins that distort the π-stack. Electrochemical sensors based on DNA CT have previously been used for single-nucleotide polymorphism detection, hybridization assays, and DNA-binding protein detection. Here, improvements to (i) the structure of DNA monolayers and (ii) the signal amplification with DNA CT platforms for improved sensitivity and detection are described.

First, improvements to the control over DNA monolayer formation are reported through the incorporation of copper-free click chemistry into DNA monolayer assembly. As opposed to conventional film formation involving the self-assembly of thiolated DNA, copper-free click chemistry enables DNA to be tethered to a pre-formed mixed alkylthiol monolayer. The total amount of DNA in the final film is directly related to the amount of azide in the underlying alkylthiol monolayer. DNA monolayers formed with this technique are significantly more homogeneous and lower density, with a larger amount of individual helices exposed to the analyte solution. With these improved monolayers, significantly more sensitive detection of the transcription factor TATA binding protein (TBP) is achieved.

Using low-density DNA monolayers, two-electrode DNA arrays were designed and fabricated to enable the placement of multiple DNA sequences onto a single underlying electrode. To pattern DNA onto the primary electrode surface of these arrays, a copper precatalyst for click chemistry was electrochemically activated at the secondary electrode. The location of the secondary electrode relative to the primary electrode enabled the patterning of up to four sequences of DNA onto a single electrode surface. As opposed to conventional electrochemical readout from the primary, DNA-modified electrode, a secondary microelectrode, coupled with electrocatalytic signal amplification, enables more sensitive detection with spatial resolution on the DNA array electrode surface. Using this two-electrode platform, arrays have been formed that facilitate differentiation between well-matched and mismatched sequences, detection of transcription factors, and sequence-selective DNA hybridization, all with the incorporation of internal controls.

For effective clinical detection, the two working electrode platform was multiplexed to contain two complementary arrays, each with fifteen electrodes. This platform, coupled with low density DNA monolayers and electrocatalysis with readout from a secondary electrode, enabled even more sensitive detection from especially small volumes (4 μL per well). This multiplexed platform has enabled the simultaneous detection of two transcription factors, TBP and CopG, with surface dissociation constants comparable to their solution dissociation constants.

With the sensitivity and selectivity obtained from the multiplexed, two working electrode array, an electrochemical signal-on assay for activity of the human methyltransferase DNMT1 was incorporated. DNMT1 is the most abundant human methyltransferase, and its aberrant methylation has been linked to the development of cancer. However, current methods to monitor methyltransferase activity are either ineffective with crude samples or are impractical to develop for clinical applications due to a reliance on radioactivity. Electrochemical detection of methyltransferase activity, in contrast, circumvents these issues. The signal-on detection assay translates methylation events into electrochemical signals via a methylation-specific restriction enzyme. Using the two working electrode platform combined with this assay, DNMT1 activity from tumor and healthy adjacent tissue lysate were evaluated. Our electrochemical measurements revealed significant differences in methyltransferase activity between tumor tissue and healthy adjacent tissue.

As differential activity was observed between colorectal tumor tissue and healthy adjacent tissue, ten tumor sets were subsequently analyzed for DNMT1 activity both electrochemically and by tritium incorporation. These results were compared to expression levels of DNMT1, measured by qPCR, and total DNMT1 protein content, measured by Western blot. The only trend detected was that hyperactivity was observed in the tumor samples as compared to the healthy adjacent tissue when measured electrochemically. These advances in DNA CT-based platforms have propelled this class of sensors from the purely academic realm into the realm of clinically relevant detection.

CYCLIC PEPTIDE HEPATOTOXINS FROM FRESH-WATER CYANOBACTERIAL (BLUE-GREEN-ALGAE) WATERBLOOMS COLLECTED IN CENTRAL CHINA

Toxic cyanobacteria (blue-green algae) waterblooms have been found in several Chinese water bodies since studies began there in 1984. Waterbloom samples for this study contained Anabaena circinalis, Microcystis aeruginosa and Oscillatoria sp. Only those waterblooms dominated by Microcystis aeruginosa were toxic by the intraperitoneal (i.p.) mouse bioassay. Signs of poisoning were the same as with known hepatotoxic cyclic peptide microcystins. One toxic fraction was isolated from each Microcystis aeruginosa sample. Two hepatotoxic peptides were purified from each of the fractions by high-performance liquid chromatography and identified by amino acid analysis followed by low and high resolution fast-atom bombardment mass spectrometry (FAB-MS). LD50 i.p. mouse values for the two toxins were 245-mu-g/kg (Toxin A) and 53-mu-g/g (Toxin B). Toxin content in the cells was 0.03 to 3.95 mg/g (Toxin A) and 0.18 to 3.33 mg/kg (Toxin B). The amino acid composition of Toxin A was alanine [1], arginine [2], glutamic acid [1] and beta-methylaspartic acid [1]; for Toxin B it was the same, except one of the arginines was replaced with a leucine. Low- and high-resolution FAB-MS showed that the molecular weights were 1,037 m/z (Toxin A) and 994 m/z (Toxin B), with formulas of C49H76O12N13 (Toxin A) and C49H75O12N10 (Toxin B). It was concluded that Toxin A is microcystin-RR and Toxin B is microcystin-LR, both known cyclic heptapeptide hepatotoxins isolated from cyanobacteria in other parts of the world. Sodium borohydride reduction of microcystin-RR yielded dihydro-microcystin-RR (m/z = 1,039), an important intermediate in the preparation of tritium-labeled toxin for metabolism and fate studies.

Environmental Isotopes as Indicators for Ground Water Recharge to Fractured Granite

To assess the contribution of accumulated winter precipitation and glacial meltwater to the recharge of deep ground water flow systems in fracture crystalline rocks, measurements of environmental isotope ratios, hydrochemical composition, and in situ parameters of ground water were performed in a deep tunnel. The measurements demonstrate the significance of these ground water recharge components for deep ground water flow systems in fractured granites of a high alpine catchment in the Central Alps, Switzerland. Hydrochemical and in situ parameters, as well as d18O in ground water samples collected in the tunnel, show only small temporal variations. The precipitation record of d18O shows seasonal variations of ~14‰ and a decrease of 0.23‰ ± 0.03‰ per 100 m elevation gain. d2H and d18O in precipitation are well correlated and plot close to the meteoric water line, as well as d2H and d18O in ground water samples, reflecting the meteoric origin of the latter. The depletion of 18O in ground water compared to 18O content in precipitation during the ground water recharge period indicates significant contributions from accumulated depleted winter precipitation to ground water recharge. The hydrochemical composition of the encountered ground water, Na-Ca-HCO3-SO4(-F), reflects an evolution of the ground water along the flowpath through the granite body. Observed tritium concentrations in ground water range from 2.6 to 16.6 TU, with the lowest values associated with a local negative temperature anomaly and anomalous depleted 18O in ground water. This demonstrates the effect of local ground water recharge from meltwater of submodern glacial ice. Such localized recharge from glaciated areas occurs along preferential flowpaths within the granite body that are mainly controlled by observed hydraulic active shear fractures and cataclastic faults.

Misonidazole binding in murine liver tissue: a marker for cellular hypoxia in vivo

The hepatic microcirculation is believed to cause variable cellular oxygenation within the organ. In this study a marker of cellular hypoxia was used to demonstrate liver oxygen tension gradients in vivo. Covalent binding of misonidazole adducts to cellular macromolecules is enhanced by hypoxia. Autoradiographs of liver from mice treated with radiolabeled misonidazole demonstrated enhanced binding of adducts within hepatocytes surrounding hepatic veins. Livers from both hypoxic and normal mice had characteristic autoradiographic grain patterns reflecting regional oxygen tension variation in vivo. Differential binding of misonidazole adducts formed in hypoxic cells could have an application in studies of liver physiology and biochemistry.

Evolução hidrogeoquímica das águas sulfúreas de Entre-os-Rios: avaliação preliminar

O principal objectivo desta dissertação foi avaliar a evolução hidrogeoquímica das águas minerais de Entre‐os‐Rios, para uma melhor compreensão do modelo hidrogeológico conceptual deste sistema hidromineral. Desta forma, foram coligidos diversos dados hidroquímicos, quer das nascentes clássicas (Torre, Curveira, Ardias, Arcos Esquerda e Arcos Direita), quer do furo Barbeitos. Foram compiladas e analisadas oitenta análises hidroquímicas no período 1938‐2012, incluindo características organolépticas (cheiro, cor e turbidez), diversas propriedades físico-químicas (temperatura, pH, condutividade eléctrica, sulfuração, etc), os principais catiões e aniões (bicarbonato, fluoreto, sódio, lítio, etc) e os elementos vestigiários (chumbo, tungsténio, boro, etc). Além disso, foram integrados os dados históricos disponíveis de finais do século XIX e inícios do século XX. Foram igualmente reunidos e discutidos alguns dados isotópicos (oxigénio‐18, deutério e trítio). O recurso hidromineral de Entre‐os‐Rios está condicionado pela litologia e pelas condições tectónicas. As análises químicas revelaram que as águas minerais de Entre‐os‐Rios apresentam uma estabilidade química nos últimos 100 anos. Estas águas são orto‐ a hipertermais, fracamente mineralizadas, de reacção alcalina, sulfídricas, bicarbonatadas sódicas, carbonatadas e muito fluoretadas. Estas características são claramente distintas das águas normais da região. As águas de Entre‐os‐Rios são muito semelhantes às águas minerais de S. Vicente e, em diversos parâmetros, bastante diferentes das águas minerais das Caldas da Saúde. Os dados isotópicos permitiram concluir que as águas de Entre‐os‐Rios têm uma origem meteórica, com um tempo de residência longo no sistema aquífero, e que são, muito provavelmente, submodernas, com uma recarga anterior a 1952. Na região de Entre‐os‐Rios coexistem três sistemas aquíferos, um sistema granítico superficial, livre e um sistema livre a semi‐confinado, ambos com circulação de águas normais, e um sistema aquífero granítico, profundo, confinado, com circulação de água mineral.

Effects of Soil Compaction and Organic Carbon Content on Preferential Flow in Loamy Field Soils

Abstract: Preferential flow and transport through macropores affect plant water use efficiency and enhance leaching of agrochemicals and the transport of colloids, thereby increasing the risk for contamination of groundwater resources. The effects of soil compaction, expressed in terms of bulk density (BD), and organic carbon (OC) content on preferential flow and transport were investigated using 150 undisturbed soil cores sampled from 15 × 15–m grids on two field sites. Both fields had loamy textures, but one site had significantly higher OC content. Leaching experiments were conducted in each core by applying a constant irrigation rate of 10 mm h−1 with a pulse application of tritium tracer. Five percent tritium mass arrival times and apparent dispersivities were derived from each of the tracer breakthrough curves and correlated with texture, OC content, and BD to assess the spatial distribution of preferential flow and transport across the investigated fields. Soils from both fields showed strong positive correlations between BD and preferential flow. Interestingly, the relationships between BD and tracer transport characteristics were markedly different for the two fields, although the relationship between BD and macroporosity was nearly identical. The difference was likely caused by the higher contents of fines and OC at one of the fields leading to stronger aggregation, smaller matrix permeability, and a more pronounced pipe-like pore system with well-aligned macropores.