Refine
Year of publication
Document Type
- Doctoral Thesis (25) (remove)
Language
- English (25) (remove)
Has Fulltext
- yes (25)
Is part of the Bibliography
- no (25)
Keywords
- Activation (1)
- Beam dynamics (1)
- CERN PSB (1)
- Chopper (1)
- FAIR (1)
- Fluka (1)
- Gabor lens (1)
- Ion Beam (1)
- LINAC (1)
- Ladder-RFQ (1)
Cancer has become one of the most fatal diseases. The Heidelberg Heavy Ion Cancer Therapy (HICAT) has the potential to become an important and efficient treatment method because of its excellent “Bragg peak” characteristics and on-line irradiation control by the PET diagnostics. The dedicated Heidelberg Heavy Ion Cancer Therapy Project includes two ECR ion sources, a RF linear injector, a synchrotron and three treatment rooms. It will deliver 4*10 high 10 protons, or 1*10 high 10 He, or 1*10 high 9 Carbons, or 5*10 high 8 Oxygens per synchrotron cycle with the beam energy 50-430AMeV for the treatments. The RF linear injector consists of a 400AkeV RFQ and of a very compact 7AMeV IH-DTL accelerator operated at 216.816MHz. The development of the IH-DTL within the HICAT project is a great challenge with respect to the present state of the DTL art because of the following reasons: • The highest operating frequency (216.816MHz) of all IH-DTL cavities; • Extremely large cavity length to diameter ratio of about 11; • IH-DTL with three internal triplets; • The highest effective voltage gain per meter (5.5MV/m); • Very short MEBT design for the beam matching. The following achievements have been reached during the development of the IH-DTL injector for HICAT : The KONUS beam dynamics design with LORASR code fulfills the beam requirement of the HICAT synchrotron at the injection point. The simulations for the IH-DTL injector have been performed not only with a homogeneous input beam, but also with the actual particle distribution from the exit of the HICAT RFQ accelerator as delivered by the PARMTEQ code. The output longitudinal normalized emittance for 95% of all particles is 2.00AkeVns, the emittance growth is less than 24%, while the X-X’ and Y-Y’ normalized emittance are 0.77mmmrad and 0.62mmmrad, respectively. The emittance growth in X-X’ is less than 18%, and the emittance growth in Y-Y’ is less than 5%. Based on the transverse envelopes of the transported particles, the redesign of the buncher drift tubes at the RFQ high energy end has been made to get a higher transit time factor for this novel RFQ internal buncher. An optimized effective buncher gap voltage of 45.4KV has been calculated to deliver a minimized longitudinal beam emittance, while the influence of the effective buncher voltage on the transverse emittance can be neglected. Six different tuning concepts were investigated in detail while tuning the 1:2 scaled HICAT IH model cavity. ‘Volume Tuning’ by a variation of the cavity cross sectional area can compensate the unbalanced capacitance distribution in case of an extreme beta-lambda-variation along an IH cavity. ‘Additional Capacitance Plates’ or copper sheets clamped on drift tube stems are a fast way for checking the tuning sensitivity, but they will be replaced by massive copper blocks mounted on the drift tube girders finally. ‘Lens Coupling’ is an important tuning to stabilize the operation mode and to increase or decrease the coupling between neighboring sections. ‘Tube Tuning’ is the fine tuning concept and also the standard tuning method to reach the needed field distributions as well as the gap voltage distributions. ‘Undercut Tuning’ is a very sensitive tuning for the end sections and with respect to the voltage distribution balance along the structure. The different types of ‘plungers’ in the 3rd and 4th sections have different effects on the resonance frequency and on the field distribution. The different triplet stems and the geometry of the cavity end have been also investigated to reach the design field and voltage distributions. Finally, the needed uniform field distribution along the IH-DTL cavity and the corresponding effective voltage distribution were realized, the remaining maximum gap voltage difference was less than 5% for the model cavity. The several important higher order modes were also measured. The RF tuning of the IH-DTL model cavity delivers the final geometry parameters of the IH-DTL power cavity. A rectangular cavity cross section was adopted for the first time for this IH-DTL cavity. This eases the realization of the volume tuning concept in the 1st and 2nd sections. Lens coupling determines the final distance between the triplet and the girder. The triplets are mounted on the lower cavity half shell. The Microwave Studio simulations have been carried out not only for the HICAT model cavity, but also for the final geometry of the IH-DTL power cavity. The field distribution for the operation mode H110 fits to the model cavity measurement as well as the Higher Order Modes. The simulations prove the IH-DTL geometrical design. On the other hand, the precision of one simulation with 2.3 million mesh points for full cross section area and the CPU time more than 15hours on a DELL PC with Intel Pentium 4 of 2.4GHz and 2.096GRAM were exploited to their limit when calculating the real parameters for the two final machining iterations during production. The shunt impedance of the IH-DTL power cavity is estimated by comparison with the existing tanks to about 195.8MOmega/m, which fits to the simulation result of 200.3MOmega/m with reducing the conductivity to the 5.0*10 high 7 Omega-1m-1. The effective shunt impedance is 153 MOmega/m. The needed RF power is 755kW. The expected quality factor of the IH-DTL cavity is about 15600. The IH-DTL power cavity tuning measurements before cavity copper plating have been performed. The results are within the specifications. There is no doubt that the needed accuracy of the voltage distribution will be reached with the foreseen fine tuning concepts in the last steps.
Studies and measurements of linear coupling and nonlinearities in hadron circular accelerators
(2006)
In this thesis a beam-based method has been developed to measure the strength and the polarity of corrector magnets (skew quadrupoles and sextupoles) in circular accelerators. The algorithm is based on the harmonic analysis (via FFT) of beam position monitor (BPM) data taken turn by turn from an accelerator in operation. It has been shown that, from the differences of the spectral line amplitudes between two consecutive BPMs, both the strength and the polarity of non-linear elements placed in between can be measured. The method has been successfully tested using existing BPM data from the SPS of CERN, since presently the SIS-18 is not equipped with the necessary hardware. The magnet strength of seven SPS extraction sextupoles was measured with a precision of about 10%. The polarities have been unambiguously measured. This method can be used to detect polarity errors and wrong power supply connections during machine commissioning, as well as for a continuous monitoring of the "nonlinearity budget" in superconducting machines. A second beam-based method has been studied for a fast measurement and correction of betatron coupling driven by skew quadrupole field errors and tilted focusing quadrupoles. Traditional methods usually require a time-consuming scan of the corrector magnets in order to minimize the coupling stop band |C|. In this thesis it has been shown how the same correction can be performed in a single machine cycle from the harmonic analysis of multi-BPM data. The method has been successfully applied to RHIC. It has been shown that the stop band |C| (also known in the American literature as Delta-Qmin) measured in a single machine cycle with the new algorithm is compatible with the value obtained by traditional methods. The measurement of the resonance phase Theta defines automatically the best corrector setting, which was found in agreement with the one obtained with a traditional scan. A third theoretical achievement is a new description of the betatron motion close to the difference resonance in presence of linear coupling. Compared to the matrix formalism the motion is parametrized as a function of the resonance driving term f1001 only (which is proven to be an observable), whereas making use of the matrix approach four parameters need to be measured. Formulae describing the exchange of RMS emittances when approaching the resonances have been already derived in the 70s in the smooth approximation. New formulae have been derived here making use of Lie algebra providing a better description of the emittance behavior. The emittance exchange curves are predicted by new formulae with excellent agreement with multi-particle simulations and the counter-intuitive emittance variation along the ring of the emittance is proven to be related to the variation of f1001. A new way to decouple the equations of motion and explicit expressions for the individual single particle invariants have been found. For the first time emittance exchange studies have been carried out in the SIS-18 of GSI. Transverse RMS emittances have been measured during 2005 from rest gas monitor (RGM) data. Crossing the linear coupling resonance, the transverse emittances exchange completely. It has been observed that this effect is reversible. Applications of this manipulation are: emittance equilibration under consideration for future operations of the SIS-18 as booster for the SIS-100; emittance transfer during multi-turn injection to improve the eficiency and to protect the injection septum in high intensity operations, by shifting part of the horizontal emittance into the vertical plane. The emittance exchange curves obtained experimentally have been compared with analytic formulae providing a fast measurement (in few machine cycles only) of the linear coupling stop band |C|. Technical problems prevented the use of the eight skew quadrupoles installed in the SIS-18 to compensate the linear coupling resonance. It has been observed that the emittance exchange curve is highly sensitive to the beam intensity. Multi-particle simulations with 2D PIC space-charge solver have been run to infer heuristic scaling laws able to quantify the observable stop band, to be used for the resonance compensation. The analysis of BPM and RGM data has been performed making use of new software applications developed for this purpose. The bpm2rdt code for the harmonic analysis of BPM data has been written and tested with real data. The software reads the BPM turn-by-turn data and the Twiss parameters. Then it performs the FFT of these data, finds the peaks of the Fourier spectra and infers the RDT fjklm, the strengths ^hjklm and the local terms lambda-jklm. All these observables are printed out together with the corresponding values of the model, computed from the nominal values of strengths and the Twiss parameters. From the FFT of dual-plane BPM data the linear optics (beta functions and phase advances Delta phi) at the corresponding location is also inferred. From the measurement of f1000, the linear coupling coeffcient C (amplitude and phase) is also computed. The code has been tested by using existing SPS data and new RHIC data. For the on-line analysis of RGM data the rgm2emitt code has been written. The application reads in input the raw data files from the RGM and the beam loss monitor (BLM) respectively, the latter created by the RGM on-line software itself. From the RGM data the transverse beam sizes and emittances are inferred and used together with the BLM data to compute the tune shift during the machine cycle.
In the framework of this thesis the intense low energy ion beam transport was investigated. Especially, the beam transport in toroidal magnetic field configurations was discussed, as it may allow the accumulation of high intensive beams in the future. One of the specific tasks is to design an injection system that can be used for the proposed low energy accumulator ring. This thesis regarding beam transport investigations is related to the larger research fields, storage rings used in accelerator physics and non-neutral plasmas. The proposal of building a storage ring with longitudinal guiding magnetic fields was made. Due to natural transversal focussing in magnetic fields it is possible to accumulate very intense charged particle beams, a subject of interest within the physics community. A simulation code (TBT) was written to describe the particle motion in curved segments. Particle in Cell techniques were utilized to simulate a multi particle dynamics. This code allows the user to generate different particle distributions as input parameter. A possibility of reading an external data file was made available so that a measured distribution can be used to compare simulation results with measured ones. A second order cloud in cell method was used to calculate charge density and in turn to solve Poisson’s equation. The circular toroidal coordinate system was used. The drift motion and gyrating motion was proved to be consistent with analytical values. Further simulations were performed to study the self field effects on beam transport. The experiments with single toroidal segments find niche in the work. The experiments were performed to compare the simulation results and gain practical experience. The toroidal segment has similar dimensions (major axis R = 1:3 m, minor axis r = 0:1 m, arc angle 30°) as for a full scale ring design. The main difference lies in the magnetic field strength. The available segments can be operated at room temperature producing 0:6T on axis maximum magnetic field, while for the storage ring design this value is in the range of 5T. The preparatory experiments consisted of building and characterization of the ion source in a first step. Along with the momentum spectrometer and emittance scanner the beam properties were studied. Low mass ion beams He+ and mixed p, H2+, H3+ beams were analyzed. The proton beam consisting of a 48% H+ fraction was extracted regularly and used for further experiments. A moderate beam energy of 10 keV was chosen as operational energy for which 3.08 mA proton beam current was measured. In the second stage, beams were transported through a solenoid and the phase space distribution was measured as a function of the magnetic field for different beam energies. The phase-space as distributions measured in a first stage were simulated backward and then again forward transported through the solenoid. The simulated results were then compared with the measured distribution. The LINTRA transport program was used. The phase-space distribution was further simulated for transport experiments in a toroidal magnetic field. The experiments with a single toroidal segment give basic results necessary to compare the results between transport code (TBT) and measurements. The optical diagnostic provides measurements which can be well compared with the simulated results. A digital camera with a magnetic shield was used to record images in jpeg file format. A subroutine was written to analyze an image file to give the intensity distribution of a given image file. The integrated profile in vertical and horizontal direction was used to calculate the vertical drift and the beam size. The simulated values were in good agreement with the measured ones. The injection system needs most care. The transport program that was used to simulate the beam in the toroid was also used to design the injection system. The injection system with its special field configurations was designed to perform experiments with room temperature segments. The main point to tackle was to smoothly bring the charged particles generated outside the trap into the acceptance of the ring. The designed system consists of two sources, one representing a ring beam and the other one the injection beam. While simulations showed a clear way, how to inject the particle beam via a well positioned solenoid and in combination with a transverse electric field element causing an ExB drift into the main ring acceptance. After construction of these injection elements it will be very important to measure the robustness of such a system with respect to the beam stability- especially of the injection channel.
Design and optimization of the lattice of the superconducting synchrotron SIS300 for slow extraction
(2011)
The superconducting synchrotron SIS300 is planned to be built at the new Facility for Antiproton and Ion Research (FAIR), at GSI-Darmstadt [1]. SIS300 will be a versatile machine, which by means of a low-energy stretcher-mode or a high-energy ramped-mode will provide slowly extracted heavy ion beams towards the experimental areas. To reach the required maximum field of 4.5 T, cos(θ) magnets are necessary. Thus, SIS300 will become the first superconducting synchrotron worldwide with cos(θ) magnets providing resonant slow extraction.
Since SIS300 will be installed in the same tunnel as the SIS100 synchrotron, the dipole layout of SIS300 cannot be freely chosen. Thus, a standard lattice cannot be applied. A redesign of the SIS300 lattice accepting compromises concerning the positions and phase advances between the optical elements has been proposed. Using the analytical model of the slow extraction, firstly proposed by Kobayashi, and the analytical description of the resonance driving modes, a multiobjective optimization algorithm has been developed for the optimization of the lattice under the given boundary conditions. The final goal of the lattice optimization is a higher efficiency of the slow extraction. The results are evaluated by means of tracking simulations performed with the code Elegant.
The field quality in superconducting cos(θ) magnets is determined by the positions of the superconducting cable and the static and time-dependent effects of the current in the cable. Furthermore, the fast ramp rates of 1 T/s in the dipoles, which are fifty times faster than in any other superconducting cos(θ) magnet, together with the fact that the aperture is smaller than in conventional accelerator magnets, makes it extremely difficult to obtain a high-quality magnetic field. The unavoidable field errors affect the beam dynamics and worsen the slow extraction efficiency. Therefore, the field errors in the SIS300 dipoles have been estimated, and their effects have been taken into account in the optimization algorithm. As a result a compensation scheme has been proposed, in which time-dependent gradients in the sextupoles counteract the decay of the sextupole field errors in the dipole magnets during the slow extraction. For the limits where the compensation was no longer possible, tolerances to the magnet field errors have been determined.
Ende der 70ger Jahre, fünf Jahre nach der Einführung des ersten kommerziellen, medizinischen Computertomographen wurde die Tomographie am Los Alamos Scientific Laboratory zum ersten Mal für die Diagnose von Teilchenstrahlen angewendet. Bei der Tomographie wird aus eindimensionalen Projektionen, sogenannten Profilen, welche in möglichst vielen Winkeln um ein Objekt herum aufgenommen werden, ein zweidimensionales Abbild der Dichteverteilung (Slice oder Scheibe) approximiert. Dies ist möglich durch das bereits 1917 von Johann Radon eingeführte Fourier-Scheiben-Theorem. In der Theorie kann die zwei-dimensionale Dichteverteilung exakt ermittelt werden, wenn Projektionen mit einer unendlich feinen Auflösung über unendlich viele Winkel um ein Objekt herum in die Rekonstruktion einbezogen werden. Durch die Rekonstruktion vieler Scheiben kann ein drei-dimensionales Abbild der Dichteverteilung in einem Objekt, in diesem Fall einem Ionenstrahl, berechnet werden, sofern dieses nicht optisch dicht ist.
Die Profile in der nicht-invasiven Strahldiagnose entstehen durch CCD-Kameraaufnahmen von strahlinduzierter Fluoreszenz, welche durch den Einlass von Restgas hervorgerufen wird. Es sind aber auch Profile, welche aus anderen Methoden gewonnen werden (z.B. Gittermessungen) denkbar. An Orten mit hoher Energie ist jedoch eine nicht-invasive Form der Profilaufnahme sowohl für die Qualität des Strahls, wie auch den Schutz der Messgeräte unabdingbar.
In den letzten 40 Jahren wurden im Bereich der Strahltomographie viele wichtige Fortschritte erzielt:
1. Anfangs standen nur sehr wenige Profile zur Verfügung, so dass die Methode der gefilterten Rückprojektion(FBP), welche sich direkt aus dem Fourier-Scheiben-Theorem ableitet und welches auch in der Medizin verwendet wird, nicht angewendet werden kann. Um dieses Problem zu lösen wurden iterative Methoden wie die Algebraische Rekonstruktion (ART) und die Methode der Maximalen Entropie (MEM) für die Strahltomographie erschlossen, so dass auch mit sehr geringer Profilanzahl eine Rücktransformation möglich wurde.
2. Neben der Ortsraumtomographie wurde die Phasenraumtomografie entwickelt, so dass mittlerweile eine Rekonstruktion des sechs-dimensionalen Phasenraumes möglich ist, mit welchem ein Ionenstrahl in seiner Gesamtheit beschrieben werden kann.
3. Die Projektionen wurden lange Zeit durch Aufnahmen von mehreren festen Anschlüssen aus gewonnen (Multi-Port-Technik). Auf diese Weise ist die Anzahl der möglichen Projektionen sehr begrenzt. So entwickelte man später eine Methode welche den Strahl mit Hilfe von Quadrupolen dreht (Quad-Scan-Technik), so dass auf diese Weise von einem Anschluss aus viele Projektionen gemessen werden konnten, so dass sogar die FBP angewendet werden konnte.
4. Die meisten Bestrebungen zielten darauf ab, die Tomographie für eine nicht-invasive Emittanzmessmethode zu nutzen, welches bis heute aufgrund der großen und noch immer zunehmenden Energien in modernen Beschleunigern ein wichtiges Problem ist. Um die Tomographie zur Emittanzmessung zu verwenden, führt man eine Rekonstruktion des Phasenraumes durch. Das Problem ist, dass hierfür das a priori Wissen über die Strahltransportmatrix in die Tomographie mit einfließt, die berechnete Strahltransportmatrix
jedoch nicht mit dem tatsächlichen Strahltransport übereinstimmt, da dieser bei hohen Energien durch auftretende Raumladung nicht-linear verändert wird. Hierzu wurden gute Fortschritte in der Abschätzung der tatsächlichen Transportmatrix gemacht um die Phasenraumtomographie trotzdem mit hinreichend gutem Ergebnis durchführen zu können.
Trotz all dieser Fortschritte und Entwicklungen ist die Tomographie bis heute keine weitverbreitete Methode in der Strahldiagnose. Der Grund ist, dass das Einrichten einer Tomografie eine komplexe Abfolge etlicher Entscheidungen und weitgestreutes Wissen aus vielen unterschiedlichen Bereichen erfordert, dieser nicht zu unterschätzende Mehraufwand jedoch auch durch einen signifikanten Nutzen gerechtfertigt sein muss. Der große Nutzen der Tomographie für die Strahldiagnose und Untersuchung der Strahldynamik ist bis heute allerdings weitgehend unerkannt und weiterhin reduziert auf die Entwicklung einer nicht-invasiven Methode für die Emittanzbestimmung. Ein zweites Hindernis stellte bisher auch die Diskrepanz zwischen Genauigkeit und Platzaufwand dar (hohe Genauigkeit durch viele Projektionen mit Quad-Scan-Technik auf mehreren Metern oder niedrige Genauigkeit durch wenig Projektionen mit Multi-Port-Technik auf weniger als einem Meter). Die Tomografie kann großen Nutzen leisten für die Online-Überwachung wichtiger Maschineneparameter im Strahlbetrieb (Monitoring) als auch für detaillierte Analysen zur Strahldynamik (Modellierung) weit über die Implementierung einer nicht-invasiven Emittanzmessmethode hinaus.
Um dies zu gewährleisten Bedarf es Zweierlei. Zum einen muss die Diskrepanz zwischen Genauigkeit und Platzaufwand aufgehoben werden. Hierzu wurde im Rahmen dieser Arbeit eine rotierbare Vakuumkammer entwickelt die nach dem Vorbild medizinischer Tomographen in mehr als 5000 Winkelschritten um den Strahl herum fahren kann, dabei ein Vakuum von mindestens 10-7mbar aufrecht erhält und einen Platzbedarf von weniger als 400 mm in der Strahlstrecke einnimmt. Zum anderen muss die Implementierung der Tomografie durch eine Angabe von schematischen Schritten und Entscheidungen vereinfacht werden. Eine Strahltomographie muss immer auf ihren jeweiligen Zweck hin implementiert werden, da Einzelelemente der Tomografie wie beispielsweise Messvorrichtung und dadurch die Profilanzahl, zu verwendender Tomographiealgorithmus, zu bestimmende Parameter sich je nach Einsatz unterscheiden können. Jedoch können die dazu nötigen Entscheidungen in ein Schema eingeordnet werden, welches die Implementierung der Tomographie vereinfacht und beschleunigt. Hierzu wurde in dieser Arbeit eine Diagnosepipeline und ein Entscheidungsschema eingeführt, sowie die Implementierung nach diesem Schema am Beispiel einer Strahltomographie für die Frankfurter Neutronenquelle (FRANZ) demonstriert und die entsprechenden Fragen und Entscheidungen diskutiert. Es wird gezeigt, wie sich aus den Messdaten über die Aufbereitung der Daten durch die Tomografie die erforderlichen Standardstrahlparameter für ein Monitoring gewinnen lassen. Zusätzlich wird ein Ebenen-Modell eingeführt, über welches nicht-Standardparameter oder neu modellierte Strahlparameter für detaillierte Analysen der Strahldynamik über die Standardparameter hinaus entwickelt werden können. Diese Arbeit soll ein grundlegendes Konzept für die routinemäßige Implementierung der Tomographie in der Strahldiagnose zur Verfügung stellen. Für die Verwendung zum Monitoring im Strahlbetrieb muss die Bestimmung von Standardparametern noch wesentlich im Zeitaufwand verbessert werden. Die Verwendung der Phasenraumtomographie benötigt noch eine Idee um den arcustangensförmigen Verlauf der berechneten Phasenraumrotationswinkel mit der Forderung der FBP nach äquidistanten Projektionswinkeln verträglicher zu machen.
With the increasing energies and intensities of heavy-ion accelerator facilities, the problem of an excessive activation of the accelerator components caused by beam losses becomes more and more important. Numerical experiments using Monte Carlo transport codes are performed in order to assess the levels of activation. The heavy-ion versions of the codes were released approximately a decade ago, therefore the verification is needed to be sure that they give reasonable results. Present work is focused on obtaining the experimental data on activation of the targets by heavy-ion beams. Several experiments were performed at GSI Helmholtzzentrum für Schwerionenforschung. The interaction of nitrogen, argon and uranium beams with aluminum targets, as well as interaction of nitrogen and argon beams with copper targets was studied. After the irradiation of the targets by different ion beams from the SIS18 synchrotron at GSI, the γ-spectroscopy analysis was done: the γ-spectra of the residual activity were measured, the radioactive nuclides were identified, their amount and depth distribution were detected. The obtained experimental results were compared with the results of the Monte Carlo simulations using FLUKA, MARS and SHIELD. The discrepancies and agreements between experiment and simulations are pointed out. The origin of discrepancies is discussed. Obtained results allow for a better verification of the Monte Carlo transport codes, and also provide information for their further development. The necessity of the activation studies for accelerator applications is discussed. The limits of applicability of the heavy-ion beam-loss criteria were studied using the FLUKA code. FLUKA-simulations were done to determine the most preferable from the radiation protection point of view materials for use in accelerator components.
Studies on the focusing performance of a Gabor lens depending on nonneutral plasma properties
(2013)
The concept of the Gabor lens goes back to an idea by Dennis Gabor, who proposed a magnetron-type trap as an effective diverging lens for electron beams (collecting lens for positive ion beams).
Electrons confined inside the lens volume by orthogonal magnetic and electric fields, create an electric space charge field that causes a radial symmetric focusing force on an ion beam passing through the lens volume.
Since the beginning of the 1990s, a new design of this lens type as well as numerical models to describe the confined plasma cloud have been developed at the Institute for Applied Physics (IAP, Johann Wolfgang Goethe-University Frankfurt).
Thanks to an improved understanding of the plasma confinement as a function of the external fields, two lenses have successfully been tested for low beam currents and remain in operation.
In the scope of this work, the performance of a prototype Gabor lens for the transport of intense, i.e. space charge dominated ion beams, was investigated at the High Current Test Injector (HOSTI) of GSI Helmholtzzentrum für Schwerionenforschung GmbH for the first time.
To ensure an optimal focusing performance of the Gabor lens a homogeneous and stable electron confinement is required. Therefore, new non-interceptive diagnostic methods were developed to investigate the parameters and state of the confined nonneutral plasma column as a function of the external fields.
An essential part of the studies was the time-resolved diagnostic of an occurring plasma instability and the determination of the electron temperature via optical spectroscopy. The latter necessitated the detailed investigation of atomic excitation as well as the measurement of optical-emission cross sections.
A comparison of the results from both experiments i.e. the beam transport measurements at GSI and the diagnostic experiments performed at IAP concerning the plasma state, gave first indications of possible interaction processes between the nonneutral plasma and the ion beam.
This thesis is structured into 7 chapters:
• Chapter 2 gives an overview of the ultrashort high intensity laser interaction with matter. The laser interaction with an induced plasma is described, starting from the kinematics of single electron motion, followed by collective electron effects and the ponderamotive motion in the laser focus and the plasma transparency for the laser beam. The three different mechanisms prepared to accelerate and propagate electrons through matter are discussed. The following indirect acceleration of protons is explained by the Target Normal Sheath Acceleration (TNSA) mechanism. Finally some possible applications of laser accelerated protons are explained briefly.
• Chapter 3 deals with the modeling of geometry and field mapping of magnetic lens. Initial proton and electron distributions, fitted to PHELIX measured data are generated, a brief description of employed codes and used techniques in simulation is given, and the aberrations at the solenoid focal spot is studied.
• Chapter 4 presents a simulation study for suggested corrections to optimize the proton beam as a later beam source. Two tools have been employed in these suggested corrections, an aperture placed at the solenoid focal spot as energy selection tool, and a scattering foil placed in the proton beam to smooth the radial energy beam profile correlation at the focal spot due to chromatic aberrations. Another suggested correction has been investigated, to optimize the beam radius at the focal spot by lens geometry controlling.
• Chapter 5 presents a simulation study for the de-neutralization problem in TNSA caused by the fringing fields of pulsed magnetic solenoid and quadrupole. In this simulation, we followed an electrostatic model, wherethe evolution of both, self and mutual fields through the pulsed magnetic solenoid could be found, which is not the case in the quadrupole and only the growth of self fields could be found. The field mapping of magnetic elements is generated by the Matlab program, while the TraceWin code is employed to study the tracking through magnetic elements.
• Chapter 6 describes the PHELIX laser parameters at GSI with chirp pulse amplification technique (CPA), and Gafchromic Radiochromic film RCF) as a spatial energy resolver film detector. The results of experiments with laser proton acceleration, which were performed in two experimental areas at GSI (Z6 area and PHELIX Laser Hall (PLH)), are presented in section 6.3.
• Chapter 7 includes the main results of this work, conclusions and gives a perspective for future experimental activities.