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.

Characterisation of cellulose- and lignin-based materials using x-ray scattering methods

Wood is an important material for the construction and pulping industries. Using x-ray diffraction the microfibril angle of Sitka spruce wood was studied in the first part of this thesis. Sitka spruce (Picea sitchensis [Bong.] Carr.) is native to the west coast of North America, but due to its fast growth rate, it has also been imported to Europe. So far, its nanometre scale properties have not been systematically characterised. In this thesis the microfibril angle of Sitka spruce was shown to depend significantly on the origin of the tree in the first annual rings near the pith. Wood can be further processed to separate lignin from cellulose and hemicelluloses. Solid cellulose can act as a reducer for metal ions and it is also a porous support for nanoparticles. By chemically reducing nickel or copper in the solid cellulose support it is possible to get small nanoparticles on the surfaces of the cellulose fibres. Cellulose supported metal nanoparticles can potentially be used as environmentally friendly catalysts in organic chemistry reactions. In this thesis the size of the nickel and copper containing nanoparticles were studied using anomalous small-angle x-ray scattering and wide-angle x-ray scattering. The anomalous small-angle x-ray scattering experiments showed that the crystallite size of the copper oxide nanoparticles was the same as the size of the nanoparticles, so the nanoparticles were single crystals. The nickel containing nanoparticles were amorphous, but crystallised upon heating. The size of the nanoparticles was observed to be smaller when the reduction of nickel was done in aqueous ammonium hydrate medium compared to reduction made in aqueous solution. Lignin is typically seen as the side-product of wood industries. Lignin is the second most abundant natural polymer on Earth, and it possesses potential to be a useful material for many purposes in addition to being an energy source for the pulp mills. In this thesis, the morphology of several lignins, which were produced by different separation methods from wood, was studied using small-angle and ultra small-angle x-ray scattering. It was shown that the fractal model previously proposed for the lignin structure does not apply to most of the extracted lignin types. The only lignin to which the fractal model could be applied was kraft lignin. In aqueous solutions the average shape of the low molar mass kraft lignin particles was observed to be elongated and flat. The average shape does not necessarily correspond to the shape of the individual particles because of the polydispersity of the fraction and due to selfassociation of the particles. Lignins, and especially lignosulfonate, have many uses as dispersants, binders and emulsion stabilisers. In this thesis work the selfassociation of low molar mass lignosulfonate macromolecules was observed using small-angle x-ray scattering. By taking into account the polydispersity of the studied lignosulfonate fraction, the shape of the lignosulfonate particles was determined to be flat by fitting an oblate ellipsoidal model to the scattering intensity.

A Formalised Model of the Information and Materials Handling Activities in the Construction Process

A model of the information and material activities that comprise the overall construction process is presented, using the SADT activity modelling methodology. The basic model is further refined into a number of generic information handling activities such as creation of new information, information search and retrieval, information distribution and person-to-person communication. The viewpoint could be described as information logistics. This model is then combined with a more traditional building process model, consisting of phases such as design and construction. The resulting two-dimensional matrix can be used for positioning different types of generic IT-tools or construction specific applications. The model can thus provide a starting point for a discussion of the application of information and communication technology in construction and for measurements of the impacts of IT on the overall process and its related costs.

Magnetically induced currents and magnetic response properties of molecules

The magnetically induced currents in organic monoring and multiring molecules, in Möbius shaped molecules and in inorganic all-metal molecules have been investigated by means of the Gauge-including magnetically induced currents (GIMIC) method. With the GIMIC method, the ring-current strengths and the ring-current density distributions can be calculated. For open-shell molecules, also the spin current can be obtained. The ring-current pathways and ring-current strengths can be used to understand the magnetic resonance properties of the molecules, to indirectly identify the effect of non-bonded interactions on NMR chemical shifts, to design new molecules with tailored properties and to discuss molecular aromaticity. In the thesis, the magnetic criterion for aromaticity has been adopted. According to this, a molecule which has a net diatropic ring current might be aromatic. Similarly, a molecule which has a net paratropic current might be antiaromatic. If the net current is zero, the molecule is nonaromatic. The electronic structure of the investigated molecules has been resolved by quantum chemical methods. The magnetically induced currents have been calculated with the GIMIC method at the density-functional theory (DFT) level, as well as at the self-consistent field Hartree-Fock (SCF-HF), at the Møller-Plesset perturbation theory of the second order (MP2) and at the coupled-cluster singles and doubles (CCSD) levels of theory. For closed-shell molecules, accurate ring-current strengths can be obtained with a reasonable computational cost at the DFT level and with rather small basis sets. For open-shell molecules, it is shown that correlated methods such as MP2 and CCSD might be needed to obtain reliable charge and spin currents. The basis set convergence has to be checked for open-shell molecules by performing calculations with large enough basis sets. The results discussed in the thesis have been published in eight papers. In addition, some previously unpublished results on the ring currents in the endohedral fullerene Sc3C2@C80 and in coronene are presented. It is shown that dynamical effects should be taken into account when modelling magnetic resonance parameters of endohedral metallofullerenes such as Sc3C2@C80. The ring-current strengths in a series of nano-sized hydrocarbon rings are related to static polarizabilities and to H-1 nuclear magnetic resonance (NMR) shieldings. In a case study on the possible aromaticity of a Möbius-shaped [16]annulene we found that, according to the magnetic criterion, the molecule is nonaromatic. The applicability of the GIMIC method to assign the aromatic character of molecules was confirmed in a study on the ring currents in simple monocylic aromatic, homoaromatic, antiaromatic, and nonaromatic hydrocarbons. Case studies on nanorings, hexaphyrins and [n]cycloparaphenylenes show that explicit calculations are needed to unravel the ring-current delocalization pathways in complex multiring molecules. The open-shell implementation of GIMIC was applied in studies on the charge currents and the spin currents in single-ring and bi-ring molecules with open shells. The aromaticity predictions that are made based on the GIMIC results are compared to other aromaticity criteria such as H-1 NMR shieldings and shifts, electric polarizabilities, bond-length alternation, as well as to predictions provided by the traditional Hückel (4n+2) rule and its more recent extensions that account for Möbius twisted molecules and for molecules with open shells.

Effect of water vapour stabilization of raw materials on the quality of the inhalation product

Generation of raw materials for dry powder inhalers by different size reduction methods can be expected to influence physical and chemical properties of the powders. This can cause differences in particle size, size distribution, shape, crystalline properties, surface texture and energy. These physical properties of powders influence the behaviour of particles before and after inhalation. Materials with an amorphous surface have different surface energy compared to materials with crystalline surface. This can affect the adhesion and cohesion of particles. Changes in the surface nature of the drug particles results in a change in product performance. By stabilization of the raw materials the amorphous surfaces are converted into crystalline surfaces. The primary aim of the study was to investigate the influence of the surface properties of the inhalation particles on the quality of the product. The quality of the inhalation product is evaluated by measuring the fine particle dose (FPD). FDP is the total dose of particles with aerodynamic diameters smaller than 5,0 μm. The secondary aim of this study was to achieve the target level of the FPD and the stability of the FPD. This study was also used to evaluate the importance of the stabilization of the inhalation powders. The study included manufacturing and analysing drug substance 200 μg/dose inhalation powder batches using non-stabilized or stabilized raw materials. The inhaler formulation consisted of micronized drug substance, lactose <100μm and micronized lactose <10μm. The inhaler device was Easyhaler®. Stabilization of the raw materials was done in different relative humidity, temperature and time. Surface properties of the raw materials were studied by dynamic vapour sorption, scanning electron microscopy and three-point nitrogen adsorption technique. Particle size was studied by laser diffraction particle size analyzer. Aerodynamic particle size distribution from inhalers was measured by new generation impactor. Stabilization of all three raw materials was successful. A clear difference between nonstabilized and stabilized raw materials was achieved for drug substance and lactose <10μm. However for lactose <100μm the difference wasn’t as clear as wanted. The surface of the non-stabilized drug substance was more irregular and the particles had more roughness on the surface compared to the stabilized drug substances particles surface. The surface of the stabilized drug particles was more regular and smoother than non-stabilized. Even though a good difference between stabilized and non-stabilized raw materials was achieved, a clear evidence of the effect of the surface properties of the inhalation particles on the quality of the product was not observed. Stabilization of the raw materials didn’t lead to a higher FPD. Possible explanations for the unexpected result might be too rough conditions in the stabilization of the drug substance or smaller than wanted difference in the degree of stabilization of the main component of the product <100μm. Despite positive effects on the quality of the product were not seen there appears to be some evidence that stabilized drug substance results in smaller particle size of dry powder inhalers.

Nanoscale studies of materials using x-ray synchrotron radiation and mesoscopic modelling

X-ray synchrotron radiation was used to study the nanostructure of cellulose in Norway spruce stem wood and powders of cobalt nanoparticles in cellulose support. Furthermore, the growth of metallic clusters was modelled and simulated in the mesoscopic size scale. Norway spruce was characterized with x-ray microanalysis at beamline ID18F of the European Synchrotron Radiation Facility in Grenoble. The average dimensions and the orientation of cellulose crystallites was determined using x-ray microdiffraction. In addition, the nutrient element content was determined using x-ray fluorescence spectroscopy. Diffraction patterns and fluorescence spectra were simultaneously acquired. Cobalt nanoparticles in cellulose support were characterized with x-ray absorption spectroscopy at beamline X1 of the Deutsches Elektronen-Synchrotron in Hamburg, complemented by home lab experiments including x-ray diffraction, electron microscopy and measurement of magnetic properties with a vibrating sample magnetometer. Extended x-ray absorption fine structure spectroscopy (EXAFS) and x-ray diffraction were used to solve the atomic arrangement of the cobalt nanoparticles. Scanning- and transmission electron microscopy were used to image the surfaces of the cellulose fibrils, where the growth of nanoparticles takes place. The EXAFS experiment was complemented by computational coordination number calculations on ideal spherical nanocrystals. The growth process of metallic nanoclusters on cellulose matrix is assumed to be rather complicated, affected not only by the properties of the clusters themselves, but essentially depending on the cluster-fiber interfaces as well as the morphology of the fiber surfaces. The final favored average size for nanoclusters, if such exists, is most probably a consequence of these two competing tendencies towards size selection, one governed by pore sizes, the other by the cluster properties. In this thesis, a mesoscopic model for the growth of metallic nanoclusters on porous cellulose fiber (or inorganic) surfaces is developed. The first step in modelling was to evaluate the special case of how the growth proceeds on flat or wedged surfaces.

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.