Water-ethanol mixtures by molecular dynamics and x-ray Compton scattering

Water-ethanol mixtures are commonly used in industry and house holds. However, quite surprisingly their molecular-level structure is still not completely understood. In particular, there is evidence that the local intermolecular geometries depend significantly on the concentration. The aim of this study was to gain information on the molecular-level structures of water-ethanol mixtures by two computational methods. The methods are classical molecular dynamics (MD), where the movement of molecules can be studied, and x-ray Compton scattering, in which the scattering cross section is sensitive to the electron momentum density. Firstly, the water-ethanol mixtures were studied with MD simulations, with the mixture concentration ranging from 0 to 100%. For the simulations well-established force fields were used for the water and ethanol molecules (TIP4P and OPLS-AA, respectively). Moreover, two models were used for ethanol, rigid and non-rigid. In the rigid model the intramolecular bond lengths are fixed, whereas in the non-rigid model the lengths are determined by harmonic potentials. Secondly, mixtures with three different concentrations employing both ethanol models were studied by calculating the experimentally observable x-ray quantity, the Compton profile. In the MD simulations a slight underestimation in the density was observed as compared to experiment. Furthermore, a positive excess of hydrogen bonding with water molecules and a negative one with ethanol was quantified. Also, the mixture was found more structured when the ethanol concentration was higher. Negligible differences in the results were found between the two ethanol models. In contrast, in the Compton scattering results a notable difference between the ethanol models was observed. For the rigid model the Compton profiles were similar for all the concentrations, but for the non-rigid model they were distinct. This leads to two possibilities of how the mixing occurs. Either the mixing is similar in all concentrations (as suggested by the rigid model) or the mixing changes for different concentrations (as suggested by the non-rigid model). Either way, this study shows that the choice of the force field is essential in the microscopic structure formation in the MD simulations. When the sources of uncertainty in the calculated Compton profiles were analyzed, it was found that more statistics needs to be collected to reduce the statistical uncertainty in the final results. The obtained Compton scattering results can be considered somewhat preliminary, but clearly indicative of the behaviour of the water-ethanol mixtures when the force field is modified. The next step is to collect more statistics and compare the results with experimental data to decide which ethanol model describes the mixture better. This way, valuable information on the microscopic structure of water-ethanol mixtures can be found. In addition, information on the force fields in the MD simulations and on the ability of the MD simulations to reproduce the microscopic structure of binary liquids is obtained.

Spectroscopy and photochemistry of hydrogen peroxide and its complexes in solid rare gases

Pdf

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.

Hydrogen soil deposition and atmospheric variations in the boreal zone

This research has been prompted by an interest in the atmospheric processes of hydrogen. The sources and sinks of hydrogen are important to know, particularly if hydrogen becomes more common as a replacement for fossil fuel in combustion. Hydrogen deposition velocities (vd) were estimated by applying chamber measurements, a radon tracer method and a two-dimensional model. These three approaches were compared with each other to discover the factors affecting the soil uptake rate. A static-closed chamber technique was introduced to determine the hydrogen deposition velocity values in an urban park in Helsinki, and at a rural site at Loppi. A three-day chamber campaign to carry out soil uptake estimation was held at a remote site at Pallas in 2007 and 2008. The atmospheric mixing ratio of molecular hydrogen has also been measured by a continuous method in Helsinki in 2007 - 2008 and at Pallas from 2006 onwards. The mean vd values measured in the chamber experiments in Helsinki and Loppi were between 0.0 and 0.7 mm s-1. The ranges of the results with the radon tracer method and the two-dimensional model were 0.13 - 0.93 mm s-1 and 0.12 - 0.61 mm s-1, respectively, in Helsinki. The vd values in the three-day campaign at Pallas were 0.06 - 0.52 mm s-1 (chamber) and 0.18 - 0.52 mm s-1 (radon tracer method and two-dimensional model). At Kumpula, the radon tracer method and the chamber measurements produced higher vd values than the two-dimensional model. The results of all three methods were close to each other between November and April, except for the chamber results from January to March, while the soil was frozen. The hydrogen deposition velocity values of all three methods were compared with one-week cumulative rain sums. Precipitation increases the soil moisture, which decreases the soil uptake rate. The measurements made in snow seasons showed that a thick snow layer also hindered gas diffusion, lowering the vd values. The H2 vd values were compared to the snow depth. A decaying exponential fit was obtained as a result. During a prolonged drought in summer 2006, soil moisture values were lower than in other summer months between 2005 and 2008. Such conditions were prevailing in summer 2006 when high chamber vd values were measured. The mixing ratio of molecular hydrogen has a seasonal variation. The lowest atmospheric mixing ratios were found in the late autumn when high deposition velocity values were still being measured. The carbon monoxide (CO) mixing ratio was also measured. Hydrogen and carbon monoxide are highly correlated in an urban environment, due to the emissions originating from traffic. After correction for the soil deposition of H2, the slope was 0.49±0.07 ppb (H2) / ppb (CO). Using the corrected hydrogen-to-carbon-monoxide ratio, the total hydrogen load emitted by Helsinki traffic in 2007 was 261 t (H2) a-1. Hydrogen, methane and carbon monoxide are connected with each other through the atmospheric methane oxidation process, in which formaldehyde is produced as an important intermediate. The photochemical degradation of formaldehyde produces hydrogen and carbon monoxide as end products. Examination of back-trajectories revealed long-range transportation of carbon monoxide and methane. The trajectories can be grouped by applying cluster and source analysis methods. Thus natural and anthropogenic emission sources can be separated by analyzing trajectory clusters.

Lewis Acid-Lewis Base Mediated Metal-Free Hydrogen Activation and Catalytic Hydrogenation

Organocatalysis, the use of organic molecules as catalysts, is attracting increasing attention as one of the most modern and rapidly growing areas of organic chemistry, with countless research groups in both academia and the pharmaceutical industry around the world working on this subject. The literature review of this thesis mainly focuses on metal-free systems for hydrogen activation and organocatalytic reduction. Since these research topics are relatively new, the literature review also highlights the basic principles of the use of Lewis acid-Lewis base pairs, which do not react irreversibly with each other, as a trap for small molecules. The experimental section progresses from the first observation of the facile heterolytical cleavage of hydrogen gas by amines and B(C6F5)3 to highly active non-metal catalysts for both enantioselective and racemic hydrogenation of unsaturated nitrogen-containing compounds. Moreover, detailed studies of structure-reactivity relationships of these systems by X-ray, neutron diffraction, NMR methods and quantum chemical calculations were performed to gain further insight into the mechanism of hydrogen activation and hydrogenation by boron-nitrogen compounds.

Long-Term Excess Mortality after Aneurysmal Subarachnoid Hemorrhage

Objective There is high case-fatality rate and loss of productive life-years related to aneurysmal subarachnoid hemorrhage (aSAH) but little data on long-term survival of SAH patients. We aim to evaluate long-term excess mortality and related risk factors after aSAH. Methods One year survivors (n=3080) after aSAH from Department of Neurosurgery in Helsinki between 1980 and 2007 were reviewed for this retrospective follow-up study. Follow-up started one year after SAH and continued until death or the end of 2008 (36 960 patient-years). Mortality and relative survival ratio (RSR) were compared with matched general population. Results After 20 years, survivors of aSAH showed 18% excess mortality compared to general population. Risk factors included: old age; poor preoperative clinical condition; conservative aneurysm treatment; multiple aneurysms; and unfavourable clinical outcome at 3 months. Conclusion Even after initially favourable recovery, patients with aSAH experience excess mortality in the long run. Cardiovascular and cerebrovascular diseases are prominent in this population.