Fast and accurate numerical solutions in some problems of particle and radiation transport: synthetic acceleration for the method of short characteristics, Doppler-broadened scattering kernel, remote sensing of the cryosphere

The aim of this work is to present various aspects of numerical simulation of particle and radiation transport for industrial and environmental protection applications, to enable the analysis of complex physical processes in a fast, reliable, and efficient way. In the first part we deal with speed-up of numerical simulation of neutron transport for nuclear reactor core analysis. The convergence properties of the source iteration scheme of the Method of Characteristics applied to be heterogeneous structured geometries has been enhanced by means of Boundary Projection Acceleration, enabling the study of 2D and 3D geometries with transport theory without spatial homogenization. The computational performances have been verified with the C5G7 2D and 3D benchmarks, showing a sensible reduction of iterations and CPU time. The second part is devoted to the study of temperature-dependent elastic scattering of neutrons for heavy isotopes near to the thermal zone. A numerical computation of the Doppler convolution of the elastic scattering kernel based on the gas model is presented, for a general energy dependent cross section and scattering law in the center of mass system. The range of integration has been optimized employing a numerical cutoff, allowing a faster numerical evaluation of the convolution integral. Legendre moments of the transfer kernel are subsequently obtained by direct quadrature and a numerical analysis of the convergence is presented. In the third part we focus our attention to remote sensing applications of radiative transfer employed to investigate the Earth's cryosphere. The photon transport equation is applied to simulate reflectivity of glaciers varying the age of the layer of snow or ice, its thickness, the presence or not other underlying layers, the degree of dust included in the snow, creating a framework able to decipher spectral signals collected by orbiting detectors.

Messung des elektrischen Neutron-Formfaktors in der Reaktion 3 He(e,e'n)pp bei Q 2 = 1.58 (GeV/c) 2

Die elektromagnetischen Nukleon-Formfaktoren sind fundamentale Größen, welche eng mit der elektromagnetischen Struktur der Nukleonen zusammenhängen. Der Verlauf der elektrischen und magnetischen Sachs-Formfaktoren G_E und G_M gegen Q^2, das negative Quadrat des Viererimpulsübertrags im elektromagnetischen Streuprozess, steht über die Fouriertransformation in direkter Beziehung zu der räumlichen Ladungs- und Strom-Verteilung in den Nukleonen. Präzise Messungen der Formfaktoren über einen weiten Q^2-Bereich werden daher für ein quantitatives Verständnis der Nukleonstruktur benötigt.rnrnDa es keine freien Neutrontargets gibt, gestaltet sich die Messung der Neutron-Formfaktoren schwierig im Vergleich zu der Messung am Proton. Konsequenz daraus ist, dass die Genauigkeit der vorhandenen Daten von Neutron-Formfaktoren deutlich geringer ist als die von Formfaktoren des Protons; auch der vermessene Q^2-Bereich ist kleiner. Insbesondere der elektrische Sachs-Formfaktor des Neutrons G_E^n ist schwierig zu messen, da er aufgrund der verschwindenden Nettoladung des Neutrons im Verhältnis zu den übrigen Nukleon-Formfaktoren sehr klein ist. G_E^n charakterisiert die Ladungsverteilung des elektrisch neutralen Neutrons und ist damit besonders sensitiv auf die innere Struktur des Neutrons.rnrnIn der hier vorgestellten Arbeit wurde G_E^n aus Strahlhelizitätsasymmetrien in der quasielastischen Streuung vec{3He}(vec{e}, e'n)pp bei einem Impulsübertrag von Q^2 = 1.58 (GeV/c)^2 bestimmt. Die Messung fand in Mainz an der Elektronbeschleunigeranlage Mainzer Mikrotron innerhalb der A1-Kollaboration im Sommer 2008 statt. rnrnLongitudinal polarisierte Elektronen mit einer Energie von 1.508 GeV wurden an einem polarisierten ^3He-Gastarget, das als effektives, polarisiertes Neutrontarget diente, gestreut. Die gestreuten Elektronen wurden in Koinzidenz mit den herausgeschlagenen Neutronen detektiert; die Elektronen wurden in einem magnetischen Spektrometer nachgewiesen, durch den Nachweis der Neutronen in einer Matrix aus Plastikszintillatoren wurde der Beitrag der quasielastischen Streuung am Proton unterdrückt.rnrnAsymmetrien des Wirkungsquerschnitts bezüglich der Elektronhelizität sind bei Orientierung der Targetpolarisation in der Streuebene und senkrecht zum Impulsübertrag sensitiv auf G_E^n / G_M^n; mittels deren Messung kann G_E^n bestimmt werden, da der magnetische Formfaktor G_M^n mit vergleichsweise hoher Präzision bekannt ist. Zusätzliche Messungen der Asymmetrie bei einer Polarisationsorientierung parallel zum Impulsübertrag wurden genutzt, um systematische Fehler zu reduzieren.rnrnFür die Messung inklusive statistischem (stat) und systematischem (sys) Fehler ergab sich G_E^n = 0.0244 +/- 0.0057_stat +/- 0.0016_sys.

Emissione di compton inverso e applicazioni astrofisiche

Nella tesi vengono trattati i fenomeni di urto tra fotoni e elettroni, nella fattispecie si riconoscono tre urti differenti: lo scattering Thomson, Compton e Compton Inverso (IC). L'attenzione viene focalizzata sull'effetto IC. Una volta riportati i risultati sui tre processi sopracitati, il testo prosegue con la trattazione del problema di interesse astrofisico, ovvero la Comptonizzazione. La Comptonizzazione è il caso più generale in cui IC e effetto Compton competono, e stabilisce come lo spettro fotonico viene modificato a causa delle interazioni col plasma. Vengono in seguito trattati due interessanti applicazioni astrofisiche della Comptonizzazione: l' effetto Sunyaev-Zel'Dovich e il Synchrotron-self-Compton.

Processi di scattering in astrofisica

Lo scattering è un processo fisico che descrive le interazioni tra particelle e radiazione elettromagnetica. Esistono diversi tipi di scattering, dove la natura di questi processi dipende dall'energia della particella e della radiazione. Quando l'elettrone fermo subisce un urto elastico con un fotone di bassa energia si ha lo \emph{scattering Thomson}. In questo processo l'onda elettromagnetica accelera la particella; successivamente dopo l'urto l'onda viene diffusa con la stessa frequenza iniziale. Questo processo accade a basse energie poichè ad energie più alte prevale l'effetto dovuto allo \emph{scattering Compton}. In questo meccanismo un fotone interagisce con un elettrone cedendogli energia e in particolare si ha una variazione di frequenza per il fotone. Quando è l'elettrone a trasferire al fotone parte della sua energia si ha \emph{Inverse Compton}. Gli ultimi due processi in realtà sono situazioni limite di un caso più generale chiamato \emph{Comptonizzazione}. Considerando un plasma rarefatto di fotoni ed elettroni, la Comptonizzazione stabilisce come l'evoluzione dello spettro viene modificato a causa dell'interazione tra fotoni ed elettroni. Nel secondo capitolo di questo elaborato verrà esaminata l'emissione in banda X della radiazione elettromagnetica in seguito a fenomeni di accrescimento intorno a buchi neri; nello specifico si analizzerà l'emissione soft e hard per buchi neri galattici.

Parity violating asymmetry in the electron deuteron quasielastic scattering, strange vector and axial vector form factors and beam normal spin asymmetries

In dieser Arbeit wurde die paritätsverletzende Asymmetrie in derrnquasielastischen Elektron-Deuteron-Streuung bei Q^2=0.23 (GeV/c)^2 mitrneinem longitudinal polarisierten Elektronstrahl bei einer Energie von 315rnMeV bestimmt. Die Messung erfolgte unter Rückwärtswinkeln. Der Detektor überdeckte einen polaren Streuwinkelbereichrnzwischen 140 und 150 deg. Das Target bestand aus flüssigemrnDeuterium in einer Targetzelle mit einer Länge von 23.4 cm. Dierngemessene paritätsverletzende Asymmetrie beträgt A_{PV}^d = (-20.11 pm 0.87_{stat} pm 1.03_{syst}), wobei der erste Fehler den statistischenrnFehlereitrag und der zweite den systematischen Fehlerbeitrag beschreibt. Ausrnder Kombination dieser Messung mit Messungen der paritätsverletzendenrnAsymmetrie in der elastischen Elektron-Proton-Streuung bei gleichem Q^2rnsowohl bei Vorwärts- als auch bei Rückwärtsmessungen können diernVektor-Strange-Formfaktoren sowie der effektive isovektorielle und isoskalarernVektorstrom des Protons, der die elektroschwachen radiativen Anapolkorrekturenrnenthält, bestimmt werden. Diese Arbeit umfasst ausserdem die Bestimmungrnder Asymmetrien bei einem transversal polarisierten Elektronstrahl sowohl beirneinem Proton- als auch einem Deuterontarget unter Rückwärtswinkeln beirnImpulsüberträgen von Q^2=0.10 (GeV/c)^2, Q^2=0.23 (GeV/c)^2rnund Q^2=0.35 (GeV/c)^2. Die im Experiment beobachteten Asymmetrien werdenrnmit theoretischen Berechnungen verglichen, welche den Imaginärteil der Zweiphoton-Austauschamplitude beinhalten.

Emissione di Compton inverso e applicazioni astrofisiche

Con questa Tesi si vuole illustrare L'effetto Compton inverso e alcune delle sue applicazioni astrofisiche. Poichè l'effetto Compton inverso è un processo di scattering si è ritenuto importante spiegare, brevemente, la diffusione Thomson e la diffusione Compton per mostrare le differenze tra i tre processi. Successivamente si è studiata l'emissione per Compton inverso dovuta a una singola interazione, la potenza e lo spettro del caso generale. Tra le applicazioni astrofisiche si sono riportate la Comptonizzazione e il Synchrotron self-Compton, per poi approfondire ognuna con una esempio specifico. Nel caso della Comptonizzazione si è discusso l'effetto Sunyaev-zeldovich mentre i BL Lacertae rappresentano l'approfondimento relatico al Synchrotron self-Compton.

Emissione di Compton inverso e applicazioni astrofisiche

Descrizione dell'effetto Compton inverso (IC) e sue implicazioni nell'ambito astrofisico. Vengono anche trattati altri effetti, tutti correlati con lo scattering IC, quali Comptonizzazione, Synchrotron-self-Compton e effetto Sunyaev-Zeldovich.

The neutron background of the XENON100 dark matter experiment

TheXENON100 experiment, installed underground at the LaboratoriNazionali del Gran Sasso, aims to directly detect dark matter in the form of weakly interacting massive particles (WIMPs) via their elastic scattering off xenon nuclei. This paper presents a study on the nuclear recoil background of the experiment, taking into account neutron backgrounds from (alpha, n) reactions and spontaneous fission due to natural radioactivity in the detector and shield materials, as well as muon-induced neutrons. Based on MonteCarlo simulations and using measured radioactive contaminations of all detector components, we predict the nuclear recoil backgrounds for the WIMP search results published by theXENON100 experiment in 2011 and 2012, 0.11(-0.04)(+0.08) events and 0.17(-0.07)(+0.12) events, respectively, and conclude that they do not limit the sensitivity of the experiment.

Accurate evaluation of hadronic uncertainties in spin-independent WIMP-nucleon scattering: Disentangling two- and three-flavor effects

We show how to avoid unnecessary and uncontrolled assumptions usually made in the literature about soft SU(3) flavor symmetry breaking in determining the two-flavor nucleon matrix elements relevant for direct detection of weakly interacting massive particles (WIMPs). Based on SU(2) chiral perturbation theory, we provide expressions for the proton and neutron scalar couplings fp,nu and fp,nd with the pion-nucleon σ term as the only free parameter, which should be used in the analysis of direct detection experiments. This approach for the first time allows for an accurate assessment of hadronic uncertainties in spin-independent WIMP-nucleon scattering and for a reliable calculation of isospin-violating effects. We find that the traditional determinations of Vfpu−fnu and fpd−fnd are off by a factor of 2.

Crystalline, Mixed-Valence Manganese Analogue of Prussian Blue: Magnetic, Spectroscopic, X-ray and Neutron Diffraction Studies

The compound of stoichiometry Mn(II)3[Mn(III)(CN)6]2·zH2O (z = 12−16) (1) forms air-stable, transparent red crystals. Low-temperature single crystal optical spectroscopy and single crystal X-ray diffraction provide compelling evidence for N-bonded high-spin manganese(II), and C-bonded low-spin manganese(III) ions arranged in a disordered, face-centered cubic lattice analogous to that of Prussian Blue. X-ray and neutron diffraction show structured diffuse scattering indicative of partially correlated (rather than random) substitutions of [Mn(III)(CN)6] ions by (H2O)6 clusters. Magnetic susceptibility measurements and elastic neutron scattering experiments indicate a ferrimagnetic structure below the critical temperature Tc = 35.5 K.

Molecular-Based Magnetism in Bimetallic Two-Dimensional Oxalate-Bridged Networks. An X-ray and Neutron Diffraction Study

Bimetallic, oxalate-bridged compounds with bi- and trivalent transition metals comprise a class of layered materials which express a large variety in their molecular-based magnetic behavior. Because of this, the availability of the corresponding single-crystal structural data is essential to the successful interpretation of the experimental magnetic results. We report in this paper the crystal structure and magnetic properties of the ferromagnetic compound {[N(n-C3H7)4][MnIICrIII(C2O4)3]}n (1), the crystal structure of the antiferromagnetic compound {[N(n-C4H9)4][MnIIFeIII(C2O4)3]}n (2), and the results of a neutron diffraction study of a polycrystalline sample of the ferromagnetic compound {[P(C6D5)4][MnIICrIII(C2O4)3]}n (3). Crystal data:  1, rhombohedral, R3c, a = 9.363(3) Å, c = 49.207(27) Å, Z = 6; 2, hexagonal, P63, a = 9.482(2) Å, c = 17.827(8) Å, Z = 2. The structures consist of anionic, two-dimensional, honeycomb networks formed by the oxalate-bridged metal ions, interleaved by the templating cations. Single-crystal field dependent magnetization measurements as well as elastic neutron scattering experiments on the manganese(II)−chromium(III) samples show the existence of long-range ferromagnetic ordering behavior below Tc = 6 K. The magnetic structure corresponds to an alignment of the spins perpendicular to the network layers. In contrast, the manganese(II)−iron(III) compound expresses a two-dimensional antiferromagnetic ordering.

Preliminary approach to neutron instrument selecction at ESSB-Bilbao based on experience at ISIS molecular spectroscopy group

Collaborative efforts between the Neutronics and Target Design Group at the Instituto de Fusión Nuclear and the Molecular Spectroscopy Group at the ISIS Pulsed Neutron and Muon Source date back to 2012 in the context of the ESS-Bilbao project. The rationale for these joint activities was twofold, namely: to assess the realm of applicability of the low-energy neutron source proposed by ESS-Bilbao - for details; and to explore instrument capabilities for pulsed-neutron techniques in the range 0.05-3 ms, a time range where ESS-Bilbao and ISIS could offer a significant degree of synergy and complementarity. As part of this collaboration, J.P. de Vicente has spent a three-month period within the ISIS Molecular Spectroscopy Group, to gain hands-on experience on the practical aspects of neutron-instrument design and the requisite neutron-transport simulations. To date, these activities have resulted in a joint MEng thesis as well as a number of publications and contributions to national and international conferences. Building upon these previous works, the primary aim of this report is to provide a self-contained discussion of general criteria for instrument selection at ESS-Bilbao, the first accelerator-driven, low-energy neutron source designed in Spain. To this end, Chapter 1 provides a brief overview of the current design parameters of the accelerator and target station. Neutron moderation is covered in Chapter 2, where we take a closer look at two possible target-moderator-reflector configurations and pay special attention to the spectral and temporal characteristics of the resulting neutron pulses. This discussion provides a necessary starting point to assess the operation of ESSB in short- and long-pulse modes. These considerations are further explored in Chapter 3, dealing with the primary characteristics of ESS-Bilbao as a short- or long-pulse facility in terms of accessible dynamic range and spectral resolution. Other practical aspects including background suppression and the use of fast choppers are also discussed. The guiding principles introduced in the first three chapters are put to use in Chapter 4 where we analyse in some detail the capabilities of a small-angle scattering instrument, as well as how specific scientific requirements can be mapped onto the optimal use of ESS-Bilbao for condensed-matter research. Part 2 of the report contains additional supporting documentation, including a description of the ESSB McStas component, a detailed characterisation of moderator response and neutron pulses, and estimates ofparameters associated with the design and operation of neutron choppers. In closing this brief foreword, we wish to thank both ESS-Bilbao and ISIS for their continuing encouragement and support along the way.

Force-free twisted magnetospheres of neutron stars

Context. The X-ray spectra observed in the persistent emission of magnetars are evidence for the existence of a magnetosphere. The high-energy part of the spectra is explained by resonant cyclotron upscattering of soft thermal photons in a twisted magnetosphere, which has motivated an increasing number of efforts to improve and generalize existing magnetosphere models. Aims. We want to build more general configurations of twisted, force-free magnetospheres as a first step to understanding the role played by the magnetic field geometry in the observed spectra. Methods. First we reviewed and extended previous analytical works to assess the viability and limitations of semi-analytical approaches. Second, we built a numerical code able to relax an initial configuration of a nonrotating magnetosphere to a force-free geometry, provided any arbitrary form of the magnetic field at the star surface. The numerical code is based on a finite-difference time-domain, divergence-free, and conservative scheme, based of the magneto-frictional method used in other scenarios. Results. We obtain new numerical configurations of twisted magnetospheres, with distributions of twist and currents that differ from previous analytical solutions. The range of global twist of the new family of solutions is similar to the existing semi-analytical models (up to some radians), but the achieved geometry may be quite different. Conclusions. The geometry of twisted, force-free magnetospheres shows a wider variety of possibilities than previously considered. This has implications for the observed spectra and opens the possibility of implementing alternative models in simulations of radiative transfer aiming at providing spectra to be compared with observations.

On the absolute calibration of bench-top small-angle X-ray scattering instruments: a comparison of different standard methods

Absolute calibration relates the measured (arbitrary) intensity to the differential scattering cross section of the sample, which contains all of the quantitative information specific to the material. The importance of absolute calibration in small-angle scattering experiments has long been recognized. This work details the absolute calibration procedure of a small-angle X-ray scattering instrument from Bruker AXS. The absolute calibration presented here was achieved by using a number of different types of primary and secondary standards. The samples were: a glassy carbon specimen, which had been independently calibrated from neutron radiation; a range of pure liquids, which can be used as primary standards as their differential scattering cross section is directly related to their isothermal compressibility; and a suspension of monodisperse silica particles for which the differential scattering cross section is obtained from Porod's law. Good agreement was obtained between the different standard samples, provided that care was taken to obtain significant signal averaging and all sources of background scattering were accounted for. The specimen best suited for routine calibration was the glassy carbon sample, due to its relatively intense scattering and stability over time; however, initial calibration from a primary source is necessary. Pure liquids can be used as primary calibration standards, but the measurements take significantly longer and are, therefore, less suited for frequent use.

Elastic and inelastic scattering of fast neutrons in fusion reactor materials

In this work, the angular distributions for elastic and. inelastic scattering of fast neutrons in fusion .reactor materials have been studied. Lithium and lead material are likely to be common components of fusion reactor wall configuration design. The measurements were performed using an associated particle time-of- flight technique. The 14 and 14.44 Mev neutrons were produced by the T(d,n} 4He reaction with deuterons being accelerated in a 150kev SAMES type J accelerator at ASTON and in.the 3. Mev DYNAMITRON at the Joint Radiation Centre, Birmingham respectively. The associated alpha-particles and fast. neutrons were detected.by means of a plastic scintillator mounted on a fast focused photomultiplier tube. The samples used were extended flat plates of thicknesses up to 0.9 mean-free-path for Lithium and 1.562 mean-free-path for Lead. The differential elastic scattering cross-sections were measured for 14 Mev neutrons for various thicknesses of Lithium and Lead in the angular range from zero to; 90º. In addition, the angular distributions of elastically scattered 14,.44 Mev .neutrons from Lithium samples were studied in the same angular range. Inelastic scattering to the 4.63 Mev state in 7Li and the 2.6 Mev state, and 4.1 Mev state in 208Pb have:been :measured.The results are compared to ENDF/B-IV data files and to previous measurements. For the Lead samples the differential neutron scattering:cross-sections for discrete 3 Mev ranges and the angular distributions were measured. The increase in effective cross-section due to multiple scattering effects,as the sample thickness increased:was found to be predicted by the empirical .relation ....... A good fit to the exoerimental data was obtained using the universal constant............ The differential elastic scattering cross-section data for thin samples of Lithium and Lead were analyzed in terms of optical model calculations using the. computer code. RAROMP. Parameter search procedures produced good fits to the·cross-sections. For the case of thick samples of Lithium and Lead, the measured angular distributions of :the scattered neutrons were compared to the predictions of the continuous slowing down model.