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.

Radiation effects in supported nanoparticles

Nanotechnology applications are entering the market in increasing numbers, nanoparticles being among the main classes of materials used. Particles can be used, e.g., for catalysing chemical reactions, such as is done in car exhaust catalysts today. They can also modify the optical and electronic properties of materials or be used as building blocks for thin film coatings on a variety of surfaces. To develop materials for specific applications, an intricate control of the particle properties, structure, size and shape is required. All these depend on a multitude of factors from methods of synthesis and deposition to post-processing. This thesis addresses the control of nanoparticle structure by low-energy cluster beam deposition and post-synthesis ion irradiation. Cluster deposition in high vacuum offers a method for obtaining precisely controlled cluster-assembled materials with minimal contamination. Due to the clusters small size, however, the cluster-surface interaction may drastically change the cluster properties on deposition. In this thesis, the deposition process of metal and alloy clusters on metallic surfaces is modelled using molecular dynamics simulations, and the mechanisms influencing cluster structure are identified. Two mechanisms, mechanical melting upon deposition and thermally activated dislocation motion, are shown to determine whether a deposited cluster will align epitaxially with its support. The semiconductor industry has used ion irradiation as a tool to modify material properties for decades. Irradiation can be used for doping, patterning surfaces, and inducing chemical ordering in alloys, just to give a few examples. The irradiation response of nanoparticles has, however, remained an almost uncharted territory. Although irradiation effects in nanoparticles embedded inside solid matrices have been studied, almost no work has been done on supported particles. In this thesis, the response of supported nanoparticles is studied systematically for heavy and light ion irradiation. The processes leading to damage production are identified and models are developed for both types of irradiation. In recent experiments, helium irradiation has been shown to induce a phase transformation from multiply twinned to single-crystalline nanoparticles in bimetallic alloys, but the nature of the transition has remained unknown. The alloys for which the effect has been observed are CuAu and FePt. It is shown in this thesis that transient amorphization leads to the observed transition and that while CuAu and FePt do not amorphize upon irradiation in bulk or as thin films, they readily do so as nanoparticles. This is the first time such an effect is demonstrated with supported particles, not embedded in a matrix where mixing is always an issue. An understanding of the above physical processes is essential, if nanoparticles are to be used in applications in an optimal way. This thesis clarifies the mechanisms which control particle morphology, and paves way for the synthesis of nanostructured materials tailored for specific applications.

Temperature dependence of g values in tetragonal ionic Cu++ salts

A temperature dependence has been observed in the spin-Hamiltonian parameters of the Cu++ ion in a tetragonal crystal field and the variation has been interpreted in terms of vibronic effects.

ESR and optical absorption studies in copper diethyldithiocarbamate

ESR investigations at X band and optical-absorption measurements have been reported in single crystals of copper (n) diethyldithiocarbamate Cu[S 2CN(C2H5)2]2 diluted to 0.2% with the corresponding zinc complex. The measurements have been made both at room and liquid-oxygen temperatures. ESR measurements gave the following values for the parameters in spin Hamiltonian g11=2.1085, g=2.023(6), A63= 142.4×10-4 cm-1, A65 = 152.0×10-4 cm-1, B = 22.4×10-4 cm-1, Q~3×10-4 cm-1. Polarized optical absorption study has made possible a proper assignment of the absorption bands to their corresponding transitions. This has led to information regarding the ordering of the MO levels of the complex. The coefficients used in the MO description of the complex have been calculated from the observed parameters. The results show that the metal ligand BIσ bond is purely covalent and that the out-of-plane w bonding is appreciably covalent whereas the in-plane Π bonding is ionic. Further, it is noted that the metal ligand binding is more covalent with sulfur as ligand than with oxygen or nitrogen.

A new class of separators and planarity of chordal graphs

We introduce a new class of clique separators, called base sets, for chordal graphs. Base sets of a chordal graph closely reflect its structure. We show that the notion of base sets leads to structural characterizations of planar k-trees and planar chordal graphs. Using these characterizations, we develop linear time algorithms for recognizing planar k-trees and planar chordal graphs. These algorithms are extensions of the Lexicographic_Breadth_First_Search algorithm for recognizing chordal graphs and are much simpler than the general planarity checking algorithm. Further, we use the notion of base sets to prove the equivalence of hamiltonian 2-trees and maximal outerplanar graphs.

Undamped Oscillations of Collisionless Stellar Systems: Spheres,Spheroids and Discs.

The collisionless Boltzmann equation governing self-gravitating systems such as galaxies has recently been shown to admit exact oscillating solutions with planar and spherical symmetry. The relation of the spherically symmetric solutions to the Virial theorem, as well as generalizations to non-uniform spheres, uniform spheroids and discs form the subject of this paper. These models generalize known families of static solutions. The case of the spheroid is worked out in some detail. Quasiperiodic as well as chaotic time variation of the two axes is demonstrated by studying the surface of section for the associated Hamiltonian system with two degrees of freedom. The relation to earlier work and possible implications for the general problem of collisionless relaxation in self gravitating systems are also discussed.

Intermolecular interactions in pi -conjugated systems: Application to polyenes

A computational scheme has been developed for strongly interacting systems wherein the intermolecular interaction is introduced as a charge-induced-dipole term. Within this approximation, the model Hamiltonian is exactly solved using a valence-bond basis. The validity of the scheme has been checked by use of exact calculations on small model systems. The method has been applied to finite polyenes to study the shifts in the ground-state energies and dipole-allowed excited-state energies in the presence of neighbors. Our calculations show a red shift in the optical gap of the infinite polyene by 0.124 eV, which is rather small compared to the experimental red shift. This is traced to the larger inaccuracy in the calculated shift in the excited state. The calculated shift in the ground-state energies are more accurate and hence the method is better suited for studying the effect of intermolecular interactions on the properties of the ground state.

Hole pairing within an extended Anderson impurity model applicable to the high-Tc cuprates

We calculate the binding energy of a hole pair within the extended Anderson Hamiltonian for the high-Tc cuprates including a Cu impurity and an oxygen-derived band. The results indicate that stable hole pairs can be formed for intra-atomic and interatomic Coulomb repulsion strengths larger than 6 and 3.5 eV, respectively. It is also shown that the total hybridization strength between the Cu 3d and oxygen p band should be less than 2.5 eV. The hole pairing takes place primarily within the oxygen-derived p band. The range of parameter values for which hole pairing occurs is also consistent with the earlier photoemission results from these cuprates.

A density matrix renormalization group study of low-lying excitations of polythiophene within a Pariser-Parr-Pople model

Symmetrized density-matrix-renormalization-group calculations have been carried out, within Pariser-Parr-Pople Hamiltonian, to explore the nature of the ground and low-lying excited states of long polythiophene oligomers. We have exploited C-2 symmetry and spin parity of the system to obtain excited states of experimental interest, and studied the lowest dipole allowed excited state and lowest dipole forbidden two photon state, for different oligomer sizes. In the long system limit, the dipole allowed excited state always lies below the lowest dipole forbidden two-photon state which implies, by Kasha rule, that polythiophene fluoresces strongly. The lowest triplet state lies below two-photon state as usual in conjugated polymers. We have doped the system with a hole and an electron and obtained the charge excitation gap and the binding energy of the 1(1)B(u)(-) exciton. We have calculated the charge density of the ground, one-photon and two-photon states for the longer system size of 10 thiophene rings to characterize these states. We have studied bond order in these states to get an idea about the equilibrium excited state geometry of the system. We have also studied the charge density distribution of the singly and doubly doped polarons for longer system size, and observe that polythiophenes do not support bipolarons.

A note on the Berry phase for systems having one degree of freedom

A one-dimensional arbitrary system with quantum Hamiltonian H(q, p) is shown to acquire the 'geometric' phase gamma (C)=(1/2) contour integral c(Podqo-qodpo) under adiabatic transport q to q+q+qo(t) and p to p+po(t) along a closed circuit C in the parameter space (qo(t), po(t)). The non-vanishing nature of this phase, despite only one degree of freedom (q), is due ultimately to the underlying non-Abelian Weyl group. A physical realisation in which this Berry phase results in a line spread is briefly discussed.

Influence of crossed electric and quantizing magnetic fields on the Einstein relation in nonlinear optical, optoelectronic and related materials: Simplified theory, relative comparison and suggestion for experimental determination

An attempt is made to study the Einstein relation for the diffusivity-to-mobility ratio (DMR) under crossed fields' configuration in nonlinear optical materials on the basis of a newly formulated electron dispersion law by incorporating the crystal field in the Hamiltonian and including the anisotropies of the effective electron mass and the spin-orbit splitting constants within the framework of kp formalisms. The corresponding results for III-V, ternary and quaternary compounds form a special case of our generalized analysis. The DMR has also been investigated for II-VI and stressed materials on the basis of various appropriate dispersion relations. We have considered n-CdGeAs2, n-Hg1-xCdxTe, n-In1-xGaxAsyP1-y lattice matched to InP, p-CdS and stressed n-InSb materials as examples. The DMR also increases with increasing electric field and the natures of oscillations are totally band structure dependent with different numerical values. It has been observed that the DMR exhibits oscillatory dependences with inverse quantizing magnetic field and carrier degeneracy due to the Subhnikov-de Haas effect. An experimental method of determining the DMR for degenerate materials in the present case has been suggested. (C) 2010 Elsevier B.V. All rights reserved.

Electron-Paramagnetic-Res study of LiKSO4 in phase IV below -138 degrees C

A detailed single-crystal EPR study of phase IV of lithium potassium sulphate below -138 degrees C has been carried out using NH3+, which substitutes for K+, as the paramagnetic probe. The spin-Hamiltonian parameters have been evaluated at -140 degrees C and yield an isotropic g=2.0034; (AH)XX=(AH)YY=25.3 G and (AH)ZZ=23.8 G; (AN)XX=8.1 G, (AN)YY=21.2 G and (AN)ZZ=25.9 G. In this phase there are 12 magnetically inequivalent K+ sites and their occurrence is ascribed to the loss of a c glide.

Mooij correlation in disordered metals

We point out that the Mooij correlation follows naturally from a dynamically disordered tight-binding Hamiltonian with random modulations of both the diagonal and the off-diagonal matrix elements which are known to act in opposition. The dynamic disorder is treated exactly while the static disorder is incorporated approximately as an effective additional time-dependent disorder affecting the diffusive electron. Such a time translation of static disorder is known to manifest itself in certain limits as a renormalization of the diffusion coefficient. The calculated conductivity exhibits the Mooij correlation at high temperatures, where quantum coherence associated with the static disorder can be ignored.

Poincaré invariance of anomalous gauge theories

We establish the Poincaré invariance of anomalous gauge theories in two dimensions, for both the Abelian and non-Abelian cases, in the canonical Hamiltonian formalism. It is shown that, despite the noncovariant appearance of the constraints of these theories, Poincaré generators can be constructed which obey the correct algebra and yield the correct transformations in the constrained space.