Refine
Document Type
- Doctoral Thesis (5)
- Master's Thesis (1)
- Working Paper (1)
Language
- English (7)
Has Fulltext
- yes (7)
Is part of the Bibliography
- no (7)
Keywords
- Monte-Carlo-Simulation (7) (remove)
Institute
In the work presented herein the microscopic transport model BAMPS (Boltzmann Approach to Multi-Parton Scatterings) is applied to simulate the time evolution of the hot partonic medium that is created in Au+Au collisions at the Relativistic Heavy Ion Collider (RHIC) and in Pb+Pb collisions at the recently started Large Hadron Collider (LHC). The study is especially focused on the investigation of the nuclear modification factor R_{AA}, that quantifies the suppression of particle yields at large transverse momentum with respect to a scaled proton+proton reference, and the simultaneous description of the collective properties of the medium in terms of the elliptic flow v_{2} within a common framework.
In the present work we applied the Optically read out PArticle track Chamber, OPAC, for the measurement of radial dose distributions, d(r), around tracks of heavy ions passing through the gas-filled sensitive volume of the chamber. The measured data were compared with d(r) functions derived from data calculated with the Monte Carlo particle transport code, TRAX – which is used for the heavy ion therapy planning at GSI. To measure this quantity we have used here an optically read out time projection chamber (OPAC) with a parallel-drift field and one or several electron and light amplification stages. The two dimensional projection of the three dimensional ionization pattern caused by the ionizing particle passing through the chamber is captured by an image intensified CCD camera. The work is motivated by the role the radial dose distribution plays in the estimation of the relative biological effectiveness (RBE) of heavy ions, e.g. in radiation therapy and in radiation protection. The most successful model for high-dose irradiation with ions (applicable e.g. for heavy ion therapy) is found to be the local effect model (LEM). The present work intends to deliver measured data for one of the basic physical parameters which serve as input for the application of the local effect model: the radial dose distribution, d(r). The first goal of our measurement program was the measurement of d(r) distributions around carbon ions of different energies from 400 MeV/u down to the Bragg peak regions. We found an excellent agreement between the measured and simulated distributions at all carbon energies for the r–range in which the measurements deliver useful results. The lower limit of this range is about 100 nm and the upper limit is 6000 nm at a resolution of down to 33 nm - if scaled to water density. Despite the simplifications in the TRAX code (e.g. binary encounter theory for the emission ionization electrons), the discrepancies between the simulated and measured d(r) distributions are found to be lower than the measurement uncertainties at most measured carbon ion energies in almost the whole observed r-range. Hence, within the limitations of our measurements we can conclude that the precision of TRAX is sufficient to simulate the d(r) distributions around carbon ions to serve as input parameter for therapy planning. However, this conclusion is only valid for larger radial distances (r >100 nm). For smaller radial distances the measured data are dominated by the diffusion. Apart from carbon ion tracks, tracks of very heavy ions (40Ar, 84Kr and 238U) were also measured with OPAC. The simulated d(r) values were typically slightly or significantly higher than the measured data in the 100 nm < r < 5000 nm region. The experience has shown: the heavier or the faster the ion, the higher the discrepancies. On the one hand, we found a surprisingly good agreement between measurements and simulations if the ions had energies of around 50 MeV/u (i.e. relatively low energy). On the other hand, at higher energies, simulated data underestimate the measured ones by up to a factor of two in the region of 100 nm < r < 1000 nm for 84Kr (E = 650 MeV/u) or in the region of 100 nm < r < 6000 nm for 238U (E = 1 GeV/u). A possible reason for these discrepancies is that the BEA model, used in TRAX for the production ionization electrons, is not adequate for very heavy projectiles. The energy values of the very heavy ions were selected with the aim of comparing the track structures - and namely the d(r) distributions - of ions with largely different atomic mass but similar LET values. From the Z-dependency of the stopping power we know that for heavier ions a higher specific ion energy (expressed in MeV/u) is required to provide the same LET. For example the common LET of 315 keV/micro-m was achieved at largely different specific energy levels of 4,4 MeV/u for 12C, 65 MeV/u for 40Ar and 650 MeV/u for 84Kr ions. The difference in the track structures was expected mainly due to the different ion velocities and thus e.g. different ranges of d-electrons. This expectation could be confirmed by the measurements. The reason why - in line with the simulations - no strong differences could be observed in the d(r) distributions of the argon and krypton ions is the relatively small difference in the velocities of the both ion types in conjunction with the limited range in r, where the data can be compared. In contrary, the d(r) function of the carbon ion shows a qualitatively different behavior than the heavier ions inside the observable radius-range - in agreement with the simulations.
Lattice simulation of a center symmetric three dimensional effective theory for SU(2) Yang-Mills
(2010)
We present lattice simulations of a center symmetric dimensionally reduced effective field theory for SU(2) Yang Mills which employ thermal Wilson lines and three-dimensional magnetic fields as fundamental degrees of freedom. The action is composed of a gauge invariant kinetic term, spatial gauge fields and a potential for the Wilson line which includes a "fuzzy" bag term to generate non-perturbative fluctuations between Z(2) degenerate ground states. The model is studied in the limit where the gauge fields are set to zero as well as the full model with gauge fields. We confirm that, at moderately weak coupling, the "fuzzy" bag term leads to eigenvalue repulsion in a finite region above the deconfining phase transition which shrinks in the extreme weak-coupling limit. A non-trivial Z(N) symmetric vacuum arises in the confined phase. The effective potential for the Polyakov loop in the theory with gauge fields is extracted from the simulations including all modes of the loop as well as for cooled configurations where the hard modes have been averaged out. The former is found to exhibit a non-analytic contribution while the latter can be described by a mean-field like ansatz with quadratic and quartic terms, plus a Vandermonde potential which depends upon the location within the phase diagram. Other results include the exact location of the phase boundary in the plane spanned by the coupling parameters, correlation lengths of several operators in the magnetic and electric sectors and the spatial string tension. We also present results from simulations of the full 4D Yang-Mills theory and attempt to make a qualitative comparison to the 3D effective theory.
This paper uses the co-incidence of extreme shocks to banks’ risk to examine within country and across country contagion among large EU banks. Banks’ risk is measured by the first difference of weekly distances to default and abnormal returns. Using Monte Carlo simulations, the paper examines whether the observed frequency of large shocks experienced by two or more banks simultaneously is consistent with the assumption of a multivariate normal or a student t distribution. Further, the paper proposes a simple metric, which is used to identify contagion from one bank to another and identify “systemically important” banks in the EU.
In this work the preparation of organic donor-acceptor thin films was studied. A chamber for organic molecular beam deposition was designed and integrated into an existing deposition system for metallic thin films. Furthermore, the deposition system was extended by a load-lock with integrated bake-out function, a chamber for the deposition of metallic contacts via stencil mask technique and a sputtering chamber. For the sublimation of the organic compounds several effusion cells were designed. The evaporation characteristic and the temperature profile within the cells was studied. Additionally, a simulation program was developed, which calculates the evaporation characteristics of different cell types. The following processes were integrated: evaporation of particles, migration on the cell walls and collisions in the gas phase. It is also possible to consider a temperature gradient within the cell. All processes can be studied separately and their relative strength can be varied. To verify the simulation results several evaporation experiments with different cell types were employed. The thickness profile of the prepared thin films was measured position-dependently. The results are in good agreement with the simulation. Furthermore, the simulation program was extended to the field of electron beam induced deposition (EBID). The second part of this work deals with the preparation and characterization of organic thin films. The focus hereby lies on the charge transfer salt (BEDT-TTF)(TCNQ), which has three known structure variants. Thin films were prepared by different methods of co-evaporation and were studied with optical microscopy, X-ray diffraction and energy dispersive X-ray spectroscopy (EDX).The formation of the monoclinic phase of (BEDT-TTF)(TCNQ) could be shown. As a last part tunnel structures were prepared as first thin film devices and measured in a He4 cryostat.
Statistical physics of power flows on networks with a high share of fluctuating renewable generation
(2010)
Renewable energy sources will play an important role in future generation of electrical energy. This is due to the fact that fossil fuel reserves are limited and because of the waste caused by conventional electricity generation. The most important sources of renewable energy, wind and solar irradiation, exhibit strong temporal fluctuations. This poses new problems for the security of supply. Further, the power flows become a stochastic character so that new methods are required to predict flows within an electrical grid. The main focus of this work is the description of power flows in a electrical transmission network with a high share of renewable generation of electrical energy. To define an appropriate model, it is important to understand the general set-up of a stable system with fluctuating generation. Therefore, generation time series of solar and wind power are compared to load time series for whole Europe and the required balancing or storage capacities analyzed. With these insights, a simple model is proposed to study the power flows. An approximation to the full power flow equations is used and evaluated with Monte-Carlo simulations. Further, approximations to the distributions of power flows along the links are analytically derived. Finally, the results are compared to the power flows calculated from the generation and load data.
As a part of this thesis, a Monte Carlo-based code has been developed capable of simulating the transition of proton beam properties to neutron beam properties as it occurs in the Li-7(p, n)Be-7 reaction. It is able to reproduce not only the angle-integrated energy distributions but it is also capable of predicting the angle-dependent neutron spectra as measured at Forschungszentrum Karlsruhe (Karlsruhe, Germany) and Physikalisch-Technische Bundesanstalt (Braunschweig, Germany). Since the code retains all three spatial dimensions as well as all three velocity dimensions, it provides very detailed information on the neutron beam. The resulting data can aid in many different aspects, for example it can be used in shielding construction, or for lithium target design. In this work, the code is used to predict the neutron beam properties expected at the Frankfurt Neutron Source at Stern-Gerlach-Zentrum (FRANZ) facility. For different proton beam energies, the neutron distribution in x/p_x, y/p_y, and z/p_z is shown as well as a Mollweide projection, which illustrates the kinematic collimation effect that limits the neutron cone opening angle to less than 180 degree.