Geant4 Hadronic Cascade Models and CMS Data Analysis : Computational Challenges in the LHC era

This work belongs to the field of computational high-energy physics (HEP). The key methods used in this thesis work to meet the challenges raised by the Large Hadron Collider (LHC) era experiments are object-orientation with software engineering, Monte Carlo simulation, the computer technology of clusters, and artificial neural networks. The first aspect discussed is the development of hadronic cascade models, used for the accurate simulation of medium-energy hadron-nucleus reactions, up to 10 GeV. These models are typically needed in hadronic calorimeter studies and in the estimation of radiation backgrounds. Various applications outside HEP include the medical field (such as hadron treatment simulations), space science (satellite shielding), and nuclear physics (spallation studies). Validation results are presented for several significant improvements released in Geant4 simulation tool, and the significance of the new models for computing in the Large Hadron Collider era is estimated. In particular, we estimate the ability of the Bertini cascade to simulate Compact Muon Solenoid (CMS) hadron calorimeter HCAL. LHC test beam activity has a tightly coupled cycle of simulation-to-data analysis. Typically, a Geant4 computer experiment is used to understand test beam measurements. Thus an another aspect of this thesis is a description of studies related to developing new CMS H2 test beam data analysis tools and performing data analysis on the basis of CMS Monte Carlo events. These events have been simulated in detail using Geant4 physics models, full CMS detector description, and event reconstruction. Using the ROOT data analysis framework we have developed an offline ANN-based approach to tag b-jets associated with heavy neutral Higgs particles, and we show that this kind of NN methodology can be successfully used to separate the Higgs signal from the background in the CMS experiment.

Patient doses in ct, dental cone beam ct and projection radiography in Finland, with emphasis on paediatric patients

Diagnostic radiology represents the largest man-made contribution to population radiation doses in Europe. To be able to keep the diagnostic benefit versus radiation risk ratio as high as possible, it is important to understand the quantitative relationship between the patient radiation dose and the various factors which affect the dose, such as the scan parameters, scan mode, and patient size. Paediatric patients have a higher probability for late radiation effects, since longer life expectancy is combined with the higher radiation sensitivity of the developing organs. The experience with particular paediatric examinations may be very limited and paediatric acquisition protocols may not be optimised. The purpose of this thesis was to enhance and compare different dosimetric protocols, to promote the establishment of the paediatric diagnostic reference levels (DRLs), and to provide new data on patient doses for optimisation purposes in computed tomography (with new applications for dental imaging) and in paediatric radiography. Large variations in radiation exposure in paediatric skull, sinus, chest, pelvic and abdominal radiography examinations were discovered in patient dose surveys. There were variations between different hospitals and examination rooms, between different sized patients, and between imaging techniques; emphasising the need for harmonisation of the examination protocols. For computed tomography, a correction coefficient, which takes individual patient size into account in patient dosimetry, was created. The presented patient size correction method can be used for both adult and paediatric purposes. Dental cone beam CT scanners provided adequate image quality for dentomaxillofacial examinations while delivering considerably smaller effective doses to patient compared to the multi slice CT. However, large dose differences between cone beam CT scanners were not explained by differences in image quality, which indicated the lack of optimisation. For paediatric radiography, a graphical method was created for setting the diagnostic reference levels in chest examinations, and the DRLs were given as a function of patient projection thickness. Paediatric DRLs were also given for sinus radiography. The detailed information about the patient data, exposure parameters and procedures provided tools for reducing the patient doses in paediatric radiography. The mean tissue doses presented for paediatric radiography enabled future risk assessments to be done. The calculated effective doses can be used for comparing different diagnostic procedures, as well as for comparing the use of similar technologies and procedures in different hospitals and countries.

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.

Modeling of Surface UV Radiation Using Satellite Data

Solar UV radiation is harmful for life on planet Earth, but fortunately the atmospheric oxygen and ozone absorb almost entirely the most energetic UVC radiation photons. However, part of the UVB radiation and much of the UVA radiation reaches the surface of the Earth, and affect human health, environment, materials and drive atmospheric and aquatic photochemical processes. In order to quantify these effects and processes there is a need for ground-based UV measurements and radiative transfer modeling to estimate the amounts of UV radiation reaching the biosphere. Satellite measurements with their near-global spatial coverage and long-term data conti-nuity offer an attractive option for estimation of the surface UV radiation. This work focuses on radiative transfer theory based methods used for estimation of the UV radiation reaching the surface of the Earth. The objectives of the thesis were to implement the surface UV algorithm originally developed at NASA Goddard Space Flight Center for estimation of the surface UV irradiance from the meas-urements of the Dutch-Finnish built Ozone Monitoring Instrument (OMI), to improve the original surface UV algorithm especially in relation with snow cover, to validate the OMI-derived daily surface UV doses against ground-based measurements, and to demonstrate how the satellite-derived surface UV data can be used to study the effects of the UV radiation. The thesis consists of seven original papers and a summary. The summary includes an introduction of the OMI instrument, a review of the methods used for modeling of the surface UV using satellite data as well as the con-clusions of the main results of the original papers. The first two papers describe the algorithm used for estimation of the surface UV amounts from the OMI measurements as well as the unique Very Fast Delivery processing system developed for processing of the OMI data received at the Sodankylä satellite data centre. The third and the fourth papers present algorithm improvements related to the surface UV albedo of the snow-covered land. Fifth paper presents the results of the comparison of the OMI-derived daily erythemal doses with those calculated from the ground-based measurement data. It gives an estimate of the expected accuracy of the OMI-derived sur-face UV doses for various atmospheric and other conditions, and discusses the causes of the differences between the satellite-derived and ground-based data. The last two papers demonstrate the use of the satellite-derived sur-face UV data. Sixth paper presents an assessment of the photochemical decomposition rates in aquatic environment. Seventh paper presents use of satellite-derived daily surface UV doses for planning of the outdoor material weathering tests.

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.

Diffusion of Be, Co and Mn impurities in compound semiconductors and glassy carbon studied by the modified radiotracer technique

The main method of modifying properties of semiconductors is to introduce small amount of impurities inside the material. This is used to control magnetic and optical properties of materials and to realize p- and n-type semiconductors out of intrinsic material in order to manufacture fundamental components such as diodes. As diffusion can be described as random mixing of material due to thermal movement of atoms, it is essential to know the diffusion behavior of the impurities in order to manufacture working components. In modified radiotracer technique diffusion is studied using radioactive isotopes of elements as tracers. The technique is called modified as atoms are deployed inside the material by ion beam implantation. With ion implantation, a distinct distribution of impurities can be deployed inside the sample surface with good con- trol over the amount of implanted atoms. As electromagnetic radiation and other nuclear decay products emitted by radioactive materials can be easily detected, only very low amount of impurities can be used. This makes it possible to study diffusion in pure materials without essentially modifying the initial properties by doping. In this thesis a modified radiotracer technique is used to study the diffusion of beryllium in GaN, ZnO, SiGe and glassy carbon. GaN, ZnO and SiGe are of great interest to the semiconductor industry and beryllium as a small and possibly rapid dopant hasn t been studied previously using the technique. Glassy carbon has been added to demonstrate the feasibility of the technique. In addition, the diffusion of magnetic impurities, Mn and Co, has been studied in GaAs and ZnO (respectively) with spintronic applications in mind.

Oncolytic Adenoviruses with Radiation Therapy for Treatment of Prostate Cancer

Prostate cancer is the most common cancer in males. Although many patients with localized disease can be cured with surgery and radiotherapy, advanced disease and especially castration resistant metastatic disease remains incurable, with a median life expectancy of less than 18 months. Oncolytic adenoviruses (Ads) are a new promising treatment against cancer due to their innate capacity to kill cancer cells. Viral replication in tumor cells leads to oncolysis and production of a multiplicity of new virions that are capable of further destroying cancerous tissue. Oncolytic Ads can be modified for tumor targeted infection and replication and be armed with therapeutic transgenes to maximize the oncolytic effect. Worldwide, clinical trials with oncolytic Ads have demonstrated good safety while the antitumor efficacy remains to be improved. Importantly, the best responses have been reported when oncolytic adenoviruses have been combined with standard cancer treatments, such as chemotherapy and radiation. Further, a challenge in many virotherapy approaches has been the monitoring of virus replication in vivo. Reporter genes have been extensively used as transgenes to evaluate the biodistribution of the virus and activity of specific promoters. However, these techniques are often limited to preclinical evaluation and not amenable to human use. The aim of the thesis was to find and develop new oncolytic Ads with maximum efficacy against metastatic, castration resistant prostate cancer and study them in vitro and in vivo combined to different forms of radiation therapy. Using combination therapy, we were aiming for better antitumor efficacy with reduced side effects. Capsid modified Ads for enhanced transduction were studied. Serotype 3 targeted chimera, Ad5/3, was found to have enhanced infectivity for prostate cancer and was used for developing new viruses for the study. Correlation between Ad-encoded marker peptide secretion and simultaneous viral replication was evaluated and the effects of radiotherapy on viral replication were studied in detail. We found that the repair of double strand breaks caused by ionizing radiation was inhibited by adenoviral proteins and led to autophagic cell death. Both subcutaneous models and intrapulmonary tumor models mimicking metastatic, aggressive disease were used in vivo. Virus efficacy was evaluated by intratumoral injections. Also, intravenous administration was evaluated to study the effectiveness in metastatic disease. Oncolytic adenovirus treatment led to significant tumor growth control and increased the survival rate of the mice. These results were further improved when oncolytic Ads were combined with radiation therapy. Oncolytic Ads expressing human sodium/iodide transporter (hNIS) as a transgene were evaluated for their oncolytic potency and for the functionality of hNIS in vitro and in vivo. Monitoring of viral replication was also assessed using different imaging modalities relative to clinical use. SPECT imaging of tumor-bearing mice was evaluated and combined with simultaneous CT-scanning to obtain important anatomical information on biodistribution, also in a three-dimensional form. It was shown that hNIS-expressing adenoviruses could harbour a bi-functional transgene allowing for localization and imaging of viral replication. Targeted radiotherapy was applied by systemic radioiodide administration and resulted in iodide accumulation into Ad-infected tumor. The combination treatment showed significantly enhanced antitumor efficacy in mice bearing prostate cancer tumors. In summary, the results presented above aim to provide new treatment modalities for castration resistant prostate cancer. Molecular insights were provided for better understanding of the benefits of combined radiation therapy and oncolytic adenoviruses, which will hopefully facilitate the translation of the approach into clinical use for humans.

The function of NR3B and NR4A orphan nuclear receptors in osteoblasts

Nuclear receptors (NRs) comprise a large family of proteins that mediate the effects of small lipophilic molecules such as steroid hormones. In addition, there are a group of NRs which lack identified natural ligands and are referred as orphan NRs. In this thesis, the function of two such orphan NR families, the NR3B (ERRα, ERRβ and ERRγ) and the NR4A family (NGFI-B, Nurr1 and Nor1), was studied. NR3B and NR4A receptors regulate many biological processes such as energy metabolism and carcinogenesis. In addition, NR3B and NR4A receptors are expressed in bone. Therefore, the signaling and function of NR3B and NR4A orphan nuclear receptors was studied specifically in osteoblasts. NR4A receptors were found to be regulated by NR3B receptors and the Wnt/β-catenin signaling pathway as ERRα, ERRγ and β-catenin repressed the transcriptional activity of NR4A receptors in U2-OS cells. NGFI-B was found to repress the transcriptional activity of ERRγ in HeLa cells. The phytoestrogen equol was identified as a new agonist for ERRγ and ERRβ in PC-3, U2-OS, and SaOS-2 cells. Equol increased the transcriptional activity of ERRγ by increasing ERRγ co-activator binding and by inducing a conformational change in the ligand binding pocket of ERRγ. The growth inhibitory effect of equol on PC-3 prostate cancer cells was decreased by blocking ERRγ expression by siRNA. Therefore, ERRγ could mediate some of the beneficial health effects of equol. The Wnt/β-catenin signaling pathway is important for the differentiation and function of osteoblasts. NR3B and NR4A receptors were found to repress the transcriptional activity mediated by β-catenin in U2-OS cells. The mesenchymal stem cells (MSCs) isolated from ERRα knockout (KO) mice showed diminished proliferation and osteoblastic differentiation compared to the wild-type cells. The overexpression of ERRα in osteoblastic MC3T3-E1 cell line increased their mineralization. Bone sialoprotein (BSP) was shown to be a direct target gene for ERRα and ERRγ as the BSP promoter was activated by ERRα or ERRγ and PGC-1α in HeLa cells. The adipogenic differentiation of ERRα KO MSCs was also decreased and they expressed less adipogenic marker genes. In conclusion, the studies described in this thesis demonstrated that the transcriptional activity of NR3B and NR4A receptors can be regulated by other orphan NRs and signaling pathways in osteoblasts. NR3B receptors can also be regulated by ligands and a new agonist, equol, was identified for ERRβ and ERRγ. New roles for NR3B and NR4A were also identified as they were shown to converge with the Wnt signaling pathway in osteoblasts, ERRγ was shown to mediate the growth inhibitory effect of equol in prostate cancer cells, and ERRα was shown to regulate positively MSC proliferation, osteoblastic differentiation and adipogenesis.

Fast- and drift-ice communities in the Bothnian Bay and the impact of UVA radiation on the Baltic Sea ice ecology

The aim of this thesis was to study ecology of Baltic Sea ice from two perspectives. In the first two studies, sea-ice ecology from riverine-influenced fast ice to drift ice in the Bothnian Bay was investigated, whereas the last two studies focus on the sensitivity of sea-ice bacteria and algae to UVA examined in situ. The seasonal sea ice cover is one of the main characteristics of the Baltic Sea, and despite the brackish parental water, the ice structure is similar to polar ice with saline brine inclusions, the sea ice habitat. The decreasing seawater salinity from the northern Baltic Sea to the Bothnian Bay translates to decreasing brine volumes along the gradient, governing the size and community structure of the food webs in ice. However, the drift and fast ice in the Bothnian Bay may differ greatly in this sense, as drift ice may have been formed at more southern locations. Rafting and the formation of snow ice are common processes in the ice field of the Bothnian Bay. As evidenced in this thesis, rafting altered the vertical distribution of organisms and snow-ice formation provided habitable space in the better-illuminated, nitrogen-rich surface layer. The divergence between fast and drift ice became apparent at the more advanced stages, and chlorophyte biomass decreased from fast to drift ice, while the opposite held true for protozoan and metazoan biomass. The brine volumes affected the communities somewhat, and a higher percentage of flagellate species was generally linked to lower brine volumes, whereas chain-forming diatoms were mostly concentrated in layers with larger brine volumes. These results add to knowledge of the ecological significance of the ice cover lasting up to 7 months per year in this area. Sea-ice food webs are generally light-limited, but while increasing light irradiances typically enhance the primary production and further, the secondary production in sea ice, any increase in solar radiation also includes an increase in harmful UVA radiation. The Baltic Sea ice microbial communities were clearly sensitive to UVA and the responses were strongly linked to the earlier light history, as well as to the solar irradiances they were exposed to. The increased biomass of chlorophytes and pennate diatoms, when UVA was excluded, indicates that their normally minor contribution to the biomass in the upper layers of sea ice might be partly dictated by UVA. The effects of UVA on bacterial production in Baltic Sea ice mostly followed the responses in algal growth, but occasionally the exposure to UVA even enhanced the bacterial production. The dominant bacterial class, Flavobacteria, seemed to be UVA-tolerant, whereas all the Alpha-, Beta- and Gammaproteobacteria present in the surface layer showed UVA sensitivity. These results indicate that changes in the light field of ice may alter the community structure and affect the functioning of ice food webs, and are of importance when the effects of thinning of the ice cover are assessed.

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.

Hyperfine effects in the nuclear magnetic resonance of paramagnetic molecules

Paramagnetic, or open-shell, systems are often encountered in the context of metalloproteins, and they are also an essential part of molecular magnets. Nuclear magnetic resonance (NMR) spectroscopy is a powerful tool for chemical structure elucidation, but for paramagnetic molecules it is substantially more complicated than in the diamagnetic case. Before the present work, the theory of NMR of paramagnetic molecules was limited to spin-1/2 systems and it did not include relativistic corrections to the hyperfine effects. It also was not systematically expandable. --- The theory was first expanded by including hyperfine contributions up to the fourth power in the fine structure constant α. It was then reformulated and its scope widened to allow any spin state in any spatial symmetry. This involved including zero-field splitting effects. In both stages the theory was implemented into a separate analysis program. The different levels of theory were tested by demonstrative density functional calculations on molecules selected to showcase the relative strength of new NMR shielding terms. The theory was also tested in a joint experimental and computational effort to confirm assignment of 11 B signals. The new terms were found to be significant and comparable with the terms in the earlier levels of theory. The leading-order magnetic-field dependence of shielding in paramagnetic systems was formulated. The theory is now systematically expandable, allowing for higher-order field dependence and relativistic contributions. The prevailing experimental view of pseudocontact shift was found to be significantly incomplete, as it only includes specific geometric dependence, which is not present in most of the new terms introduced here. The computational uncertainty in density functional calculations of the Fermi contact hyperfine constant and zero-field splitting tensor sets a limit for quantitative prediction of paramagnetic shielding for now.

The Radiation, snow characteristics and albedo at summit (RASCALS) expedition report

The RASCALS expedition spent over three weeks at the Summit camp research station near the top of the Greenland Ice Sheet during polar summer 2010. During this time, detailed measurements of the physical and optical properties of Arctic perennial snow were carried out concurrently with snow albedo and reflectance measurements. Favorable weather conditions during the campaign enabled the collection of a large dataset on Arctic snow albedo and associated quantities for use in developing and validating remote sensing algorithms for snow albedo using satellites. This report provides a description of the measurements and conditions during the campaign. RASCALS-retkikunnan tehtävä oli tutkia Grönlannin mannerjäätikön lumen fysikaalisia ja optisia ominaisuuksia sekä Auringon valon vuorovaikutusta lumen kanssa. Retikunta vietti hieman yli kolme viikkoa mannerjäätikön keskellä sijaitsevalla Summit Camp-tutkimusasemalla tehden mittauksia. Sääolot suosivat kampanjaa, jonka seurauksena onnistuttiin keräämään laaja ja monipuolinen tietoaineisto mannerjäätikön lumen pintakerroksesta ja eritoten lumen heijastavuuden (albedon)käyttäytymisestä. Aineisto on hyödyllinen kehitettäessä ja varmennettaessa lumen albedon kaukokartoitusmenetelmiä satelliiteilla.