Computer simulations of charged systems in partially periodic geometries

This work presents algorithms for the calculation of the electrostatic interaction in partially periodic systems. The framework for these algorithms is provided by the simulation package ESPResSo, of which the author was one of the main developers. The prominent features of the program are listed and the internal structure is described. In the following, algorithms for the calculation of the Coulomb sum in three dimensionally periodic systems are described. These methods are the foundations for the algorithms for partially periodic systems presented in this work. Starting from the MMM2D method for systems with one non-periodic coordinate, the ELC method for these systems is developed. This method consists of a correction term which allows to use methods for three dimensional periodicity also for the case of two periodic coordinates. The computation time of this correction term is neglible for large numbers of particles. The performance of MMM2D and ELC are demonstrated by results from the implementations contained in ESPResSo. It is also discussed, how different dielectric constants inside and outside of the simulation box can be realized. For systems with one periodic coordinate, the MMM1D method is derived from the MMM2D method. This method is applied to the problem of the attraction of like-charged rods in the presence of counterions, and results of the strong coupling theory for the equilibrium distance of the rods at infinite counterion-coupling are checked against results from computer simulations. The degree of agreement between the simulations at finite coupling and the theory can be characterized by a single parameter gamma_RB. In the special case of T=0, one finds under certain circumstances flat configurations, in which all charges are located in the rod-rod plane. The energetically optimal configuration and its stability are determined analytically, which depends on only one parameter gamma_z, similar to gamma_RB. These findings are in good agreement with results from computer simulations.

Selektiver Nachweis von Uran-236 mittels hochauflösender Resonanzionisations-Massenspektrometrie

Die vorliegende Dissertation beschreibt die Entwicklung, den Aufbau und die Erprobung eines neuartigen lasermassenspektrometrischen Nachweises des Ultraspurenisotops 236U. Das Nuklid 236U wird vorwiegend in Kernreaktoren durch Neutroneneinfang aus 235U gebildet und dient damit als Schlüsselisotop zur Unterscheidung anthropogenen Urans von natürlichem Uran-Vorkommen. Mit seinem Nachweis wird die Untersuchung der Migration von Kernbrennstoff in der Umwelt und die Beantwortung kritischer Fragen in der nuklearen Forensik ermöglicht. Im Rahmen dieser Arbeit wurde das Verfahren der hochauflösenden Resonanzionisations-Massenspektrometrie auf die Anforderungen des selektiven Nachweises von Uran-Isotopen angepasst. Wesentliche Schritte waren hierbei die Untersuchung einer effizienten Atomisation von Uran-Proben, die Identifikation atomarer und autoionisierender Zustände für eine resonante Anregungsleiter, die vollständige Spezifikation der Anregungsleiter für alle Uran-Isotope und schließlich die Umsetzung der Erkenntnisse in ein analytisches Messverfahren. Die optische Selektivität des Verfahrens konnte durch Dichtematrixrechnungen zu ca. 14 Größenordnungen abgeschätzt werden; die Nachweisgrenze des Verfahrens für das Isotopenverhältnis 236U / 238U ist dabei gegenwärtig durch Untergrund begrenzt und beträgt ca. 3 · 10−8. Mit diesen Spezifikationen konnte die Linearität und Präzision des Nachweisverfahrens über einen dynamischen Bereich von vier Größenordnungen nachgewiesen werden.

Strukturbildungsverhalten kettensteifer Polyelektrolyte in wässrigen Lösungen am Beispiel von Poly(para-phenylen)sulfonat

An einer Vielzahl biogener Polyelektrolyte, wie z.B. den Nukleinsäuren DNA und RNA sowie Proteinen, ist die Ausbildung von Strukturhierarchien durch Selbstorganisation von Strukturelementen zu beobachten. Dabei wird das Strukturbildungsverhalten durch eine Kopplung von Wechselwirkungen auf verschiedenen Längenskalen, den kurzreichweitigen ausgeschlossenen Volumen und den langreichweitigen elektrostatischen Wechselwirkungen (Coulomb-Wechselwirkungen), die wiederum durch eine Vielzahl mikroskopische Parameter (z.B. Konformation) beeinflusst werden, bestimmt. Durch diese Komplexität ist es nicht möglich, den für die Strukturbildung hochgeladener Systeme bedeutsamen Beitrag der elektrostatischen Wechselwirkungen isoliert zu betrachten. Aus diesem Grund werden zur Aufklärung von Wechselwirkungs- und Strukturbildungsmechanismen vereinfachte Modell-Systeme herangezogen. Eine Möglichkeit besteht in der Verwendung synthetischer, kettensteifer Polyelektrolyte. Im Rahmen dieser Arbeit wurde das Aggregationsverhalten wässriger Lösungen dodecylsubstituierter Poly(para-phenylen)sulfonate (PPPS) sowie die Beeinflussung der Strukturbildung durch verschiedene Parameter charakterisiert. Als Einflussparameter wurden einerseits die Gegenion-Spezies und andererseits die Temperatur- und Konzentrations-Abhängigkeit untersucht. Hierzu wurden wässrige Lösungen der freien Säuren der PPPS mit Molekulargewichten zwischen MW = 18 kg/mol bis 58 kg/mol mittels Licht-, Röntgen- und Neutronenstreuung sowie durch Licht-, Polarisations-, Transmissionselektronen- und Rasterkraftmikroskopie in einem Konzentrationsbereich von 0,0008 < c < 1,1 g/L untersucht.

The bursting behavior of gas slugs: laboratory and analytical insights into Strombolian volcanic eruptions.

The aim of this thesis is to study how explosive behavior and geophysical signals in a volcanic conduit are related to the development of overpressure in slug-driven eruptions. A first suite of laboratory experiments of gas slugs ascending in analogue conduits was performed. Slugs ascended into a range of analogue liquids and conduit diameters to allow proper scaling to the natural volcanoes. The geometrical variation of the slug in response to the explored variables was parameterised. Volume of gas slug and rheology of the liquid phase revealed the key parameters in controlling slug overpressure at bursting. Founded on these results, a theoretical model to calculate burst overpressure for slug-driven eruptions was developed. The dimensionless approach adopted allowed to apply the model to predict bursting pressure of slugs at Stromboli. Comparison of predicted values with measured data from Stromboli volcano showed that the model can explain the entire spectrum of observed eruptive styles at Stromboli – from low-energy puffing, through normal Strombolian eruptions, up to paroxysmal explosions – as manifestations of a single underlying physical process. Finally, another suite of laboratory experiments was performed to observe oscillatory pressure and forces variations generated during the expansion and bursting of gas slugs ascending in a conduit. Two end-member boundary conditions were imposed at the base of the pipe, simulating slug ascent in closed base (zero magma flux) and open base (constant flux) conduit. At the top of the pipe, a range of boundary conditions that are relevant at a volcanic vent were imposed, going from open to plugged vent. The results obtained illustrate that a change in boundary conditions in the conduit concur to affect the dynamic of slug expansion and burst: an upward flux at the base of the conduit attenuates the magnitude of the pressure transients, while a rheological stiffening in the top-most region of conduit changes dramatically the magnitude of the observed pressure transients, favoring a sudden, and more energetic pressure release into the overlying atmosphere. Finally, a discussion on the implication of changing boundary on the oscillatory processes generated at the volcanic scale is also given.

Search for K isomers in 252,254 No and 260 Sg and investigation of their nuclear structure

Study of K isomerism in the transfermium region around the deformed shells at N=152, Z=102, and N=162, Z=108 provides important information on the structure of heavy nuclei. Recent calculations suggest that the K-isomerism can enhance the stability of such nuclei against alpha emission and spontaneous fission. Nuclei showing K isomerism have neutron and proton orbitals with large spin projections on the symmetry axis which is due to multi quasiparticle states with aligned spins K. Quasi-particle states are formed by breaking pairs of nucleons and raising one or two nucleons in orbitals near the Fermi surface above the gap, forming high K (multi)quasi-particle states mainly at low excitation energies. Experimental examples are the recently studied two quasi-particle K isomers in 250,256-Fm, 254-No, and 270-Ds. Nuclei in this region, are produced with cross sections ranging from several nb up to µb, which are high enough for a detailed decay study. In this work, K isomerism in Sg and No isotopes was studied at the velocity filter SHIP of GSI, Darmstadt. The data were obtained by using a new data acquisition system which was developed and installed during this work. 252,254-No and 260-Sg were produced in fusion evaporation reactions of 48-Ca and 54-Cr projectiles with 206,208-Pb targets at beam energies close to the Coulomb barrier. A new K isomer was discovered in 252-No at excitation energy of 1.25 MeV, which decays to the ground state rotational band via gamma emission. It has a half-life of about 100 ms. The population of the isomeric state was about 20% of the ground state population. Detailed investigations were performed on 254-No in which two isomeric states (275 ms and 198 µs) were already discovered by R.-D. Herzberg, but due to the higher number of observed gamma decays more detailed information about the decay path of the isomers was obtained in the present work. In 260-Sg, we observed no statistically significant component with a half life different from that of the ground state. A comparison between experimental results and theoretical calculations of the single particle energies shows a fair agreement. The structure of the here studied nuclei is in particular important as single particle levels are involved which are relevant for the next shell closure expected to form the region of the shell stabilized superheavy elements at proton numbers 114, 120, or 126 and neutron number 184. K isomers, in particular, could be an ideal tool for the synthesis and study of these isotopes due to enhanced spontaneous fission life times which could result in higher alpha to spontaneous fission branching ratios and longer half lifes.

Development of empirical potentials for amorphous silica

Development of empirical potentials for amorphous silica Amorphous silica (SiO2) is of great importance in geoscience and mineralogy as well as a raw material in glass industry. Its structure is characterized as a disordered continuous network of SiO4 tetrahedra. Many efforts have been undertaken to understand the microscopic properties of silica by classical molecular dynamics (MD) simulations. In this method the interatomic interactions are modeled by an effective potential that does not take explicitely into account the electronic degrees of freedom. In this work, we propose a new methodology to parameterize such a potential for silica using ab initio simulations, namely Car-Parrinello (CP) method [Phys. Rev. Lett. 55, 2471 (1985)]. The new potential proposed is compared to the BKS potential [Phys. Rev. Lett. 64, 1955 (1990)] that is considered as the benchmark potential for silica. First, CP simulations have been performed on a liquid silica sample at 3600 K. The structural features so obtained have been compared to the ones predicted by the classical BKS potential. Regarding the bond lengths the BKS tends to underestimate the Si-O bond whereas the Si-Si bond is overestimated. The inter-tetrahedral angular distribution functions are also not well described by the BKS potential. The corresponding mean value of theSiOSi angle is found to be ≃ 147◦, while the CP yields to aSiOSi angle centered around 135◦. Our aim is to fit a classical Born-Mayer/Coulomb pair potential using ab initio calculations. To this end, we use the force-matching method proposed by Ercolessi and Adams [Europhys. Lett. 26, 583 (1994)]. The CP configurations and their corresponding interatomic forces have been considered for a least square fitting procedure. The classical MD simulations with the resulting potential have lead to a structure that is very different from the CP one. Therefore, a different fitting criterion based on the CP partial pair correlation functions was applied. Using this approach the resulting potential shows a better agreement with the CP data than the BKS ones: pair correlation functions, angular distribution functions, structure factors, density of states and pressure/density were improved. At low temperature, the diffusion coefficients appear to be three times higher than those predicted by the BKS model, however showing a similar temperature dependence. Calculations have also been carried out on crystalline samples in order to check the transferability of the potential. The equilibrium geometry as well as the elastic constants of α-quartz at 0 K are well described by our new potential although the crystalline phases have not been considered for the parameterization. We have developed a new potential for silica which represents an improvement over the pair potentials class proposed so far. Furthermore, the fitting methodology that has been developed in this work can be applied to other network forming systems such as germania as well as mixtures of SiO2 with other oxides (e.g. Al2O3, K2O, Na2O).

The collapse of linear polyelectrolyte chains in a poor solvent: When does a collapsing polyelectrolyte collect its counter ions?

The collapse of linear polyelectrolyte chains in a poor solvent: When does a collapsing polyelectrolyte collect its counter ions? The collapse of polyions in a poor solvent is a complex system and is an active research subject in the theoretical polyelectrolyte community. The complexity is due to the subtle interplay between hydrophobic effects, electrostatic interactions, entropy elasticity, intrinsic excluded volume as well as specific counter-ion and co-ion properties. Long range Coulomb forces can obscure single molecule properties. The here presented approach is to use just a small amount of screening salt in combination with a very high sample dilution in order to screen intermolecular interaction whereas keeping intramolecular interaction as much as possible (polyelectrolyte concentration cp ≤ 12 mg/L, salt concentration; Cs = 10^-5 mol/L). This is so far not described in literature. During collapse, the polyion is subject to a drastic change in size along with strong reduction of free counterions in solution. Therefore light scattering was utilized to obtain the size of the polyion whereas a conductivity setup was developed to monitor the proceeding of counterion collection by the polyion. Partially quaternized PVP’s below and above the Manning limit were investigated and compared to the collapse of their uncharged precursor. The collapses were induced by an isorefractive solvent/non-solvent mixture consisting of 1-propanol and 2-pentanone, with nearly constant dielectric constant. The solvent quality for the uncharged polyion could be quantified which, for the first time, allowed the experimental investigation of the effect of electrostatic interaction prior and during polyion collapse. Given that the Manning parameter M for QPVP4.3 is as low as lB / c = 0.6 (lB the Bjerrum length and c the mean contour distance between two charges), no counterion binding should occur. However the Walden product reduces with first addition of non solvent and accelerates when the structural collapse sets in. Since the dielectric constant of the solvent remains virtually constant during the chain collapse, the counterion binding is entirely caused by the reduction in the polyion chain dimension. The collapse is shifted to lower wns with higher degrees of quaternization as the samples QPVP20 and QPVP35 show (M = 2.8 respectively 4.9). The combination of light scattering and conductivity measurement revealed for the first time that polyion chains already collect their counter ions well above the theta-dimension when the dimensions start to shrink. Due to only small amounts of screening salt, strong electrostatic interactions bias dynamic as well as static light scattering measurements. An extended Zimm formula was derived to account for this interaction and to obtain the real chain dimensions. The effective degree of dissociation g could be obtained semi quantitatively using this extrapolated static in combination with conductivity measurements. One can conclude the expansion factor a and the effective degree of ionization of the polyion to be mutually dependent. In the good solvent regime g of QPVP4.3, QPVP20 and QPVP35 exhibited a decreasing value in the order 1 > g4.3 > g20 > g35. The low values of g for QPVP20 and QPVP35 are assumed to be responsible for the prior collapse of the higher quaternized samples. Collapse theory predicts dipole-dipole attraction to increase accordingly and even predicts a collapse in the good solvent regime. This could be exactly observed for the QPVP35 sample. The experimental results were compared to a newly developed theory of uniform spherical collapse induced by concomitant counterion binding developed by M. Muthukumar and A. Kundagrami. The theory agrees qualitatively with the location of the phase boundary as well as the trend of an increasing expansion with an increase of the degree of quaternization. However experimental determined g for the samples QPVP4.3, QPVP20 and QPVP35 decreases linearly with the degree of quaternization whereas this theory predicts an almost constant value.

Dynamics, rheology and critical properties of colloidal fluid mixtures

Liquids under the influence of external fields exhibit a wide range of intriguing phenomena that can be markedly different from the behaviour of a quiescent system. This work considers two different systems — a glassforming Yukawa system and a colloid-polymer mixture — by Molecular Dynamics (MD) computer simulations coupled to dissipative particle dynamics. The former consists of a 50-50 binary mixture of differently-sized, like-charged colloids interacting via a screened Coulomb (Yukawa) potential. Near the glass transition the influence of an external shear field is studied. In particular, the transition from elastic response to plastic flow is of interest. At first, this model is characterised in equilibrium. Upon decreasing temperature it exhibits the typical dynamics of glassforming liquids, i.e. the structural relaxation time τα grows strongly in a rather small temperature range. This is discussed with respect to the mode-coupling theory of the glass transition (MCT). For the simulation of bulk systems under shear, Lees-Edwards boundary conditions are applied. At constant shear rates γ˙ ≫ 1/τα the relevant time scale is given by 1/γ˙ and the system shows shear thinning behaviour. In order to understand the pronounced differences between a quiescent system and a system under shear, the response to a suddenly commencing or terminating shear flow is studied. After the switch-on of the shear field the shear stress shows an overshoot, marking the transition from elastic to plastic deformation, which is connected to a super-diffusive increase of the mean squared displacement. Since the average static structure only depends on the value of the shear stress, it does not discriminate between those two regimes. The distribution of local stresses, in contrast, becomes broader as soon as the system starts flowing. After a switch-off of the shear field, these additional fluctuations are responsible for the fast decay of stresses, which occurs on a time scale 1/γ˙ . The stress decay after a switch-off in the elastic regime, on the other hand, happens on the much larger time scale of structural relaxation τα. While stresses decrease to zero after a switch-off for temperatures above the glass transition, they decay to a finite value for lower temperatures. The obtained results are important for advancing new theoretical approaches in the framework of mode-coupling theory. Furthermore, they suggest new experimental investigations on colloidal systems. The colloid-polymer mixture is studied in the context of the behaviour near the critical point of phase separation. For the MD simulations a new effective model with soft interaction potentials is introduced and its phase diagram is presented. Here, mainly the equilibrium properties of this model are characterised. While the self-diffusion constants of colloids and polymers do not change strongly when the critical point is approached, critical slowing down of interdiffusion is observed. The order parameter fluctuations can be determined through the long-wavelength limit of static structure factors. For this strongly asymmetric mixture it is shown how the relevant structure factor can be extracted by a diagonalisation of a matrix that contains the partial static structure factors. By presenting first results of this model under shear it is demonstrated that it is suitable for non-equilibrium simulations as well.

Sviluppo di una metodologia per la modellazione geologica-geomeccanica ed il monitoraggio di reservoir adibiti a siti di stoccaggio di gas naturale

La valutazione dei rischi associati all’operatività dei sistemi di stoccaggio, quali la sismicità indotta e la subsidenza, è requisito basilare per una loro corretta gestione e progettazione, e passa attraverso la definizione dell’influenza sullo stato tensionale delle variazioni di pressione di poro nel sottosuolo. Principale scopo di questo progetto è lo sviluppo di una metodologia in grado di quantificare le deformazioni dei reservoir in funzione della pressione di poro, di tarare i modelli utilizzati con casi studio che presentino dati di monitoraggio reali, tali da consentire un confronto con le previsioni di modello. In questa tesi, la teoria delle inomogeneità è stata utilizzata, tramite un approccio semianalitico, per definire le variazioni dei campi elastici derivanti dalle operazioni di prelievo e immissione di fluidi in serbatoi geologici. Estensione, forma e magnitudo delle variazioni di stress indotte sono state valutate tramite il concetto di variazione dello sforzo critico secondo il criterio di rottura di Coulomb, tramite un’analisi numerica agli elementi finiti. La metodologia sviluppata è stata applicata e tarata su due reservoir sfruttati e riconvertiti a sistemi di stoccaggio che presentano dataset, geologia, petrofisica, e condizioni operative differenti. Sono state calcolate le variazioni dei campi elastici e la subsidenza; è stata mappata la variazione di sforzo critico di Coulomb per entrambi i casi. I risultati ottenuti mostrano buon accordo con le osservazioni dei monitoraggi, suggerendo la bontà della metodologia e indicando la scarsa probabilità di sismicità indotta. Questo progetto ha consentito la creazione di una piattaforma metodologica di rapido ed efficace utilizzo, per stimare l’influenza dei sistemi di stoccaggio di gas sullo stato tensionale della crosta terrestre; in fase di stoccaggio, permette di monitorare le deformazioni e gli sforzi indotti; in fase di progettazione, consente di valutare le strategie operative per monitorare e mitigare i rischi geologici associati a questi sistemi.

Hydrodynamical simulations of early-type galaxies: effects of galaxy shape and stellar dynamics on hot coronae

Early-Type galaxies (ETGs) are embedded in hot (10^6-10^7 K), X-ray emitting gaseous haloes, produced mainly by stellar winds and heated by Type Ia supernovae explosions, by the thermalization of stellar motions and occasionally by the central super-massive black hole (SMBH). In particular, the thermalization of the stellar motions is due to the interaction between the stellar and the SNIa ejecta and the hot interstellar medium (ISM) already residing in the ETG. A number of different astrophysical phenomena determine the X-ray properties of the hot ISM, such as stellar population formation and evolution, galaxy structure and internal kinematics, Active Galactic Nuclei (AGN) presence, and environmental effects. With the aid of high-resolution hydrodynamical simulations performed on state-of-the-art galaxy models, in this Thesis we focus on the effects of galaxy shape, stellar kinematics and star formation on the evolution of the X-ray coronae of ETGs. Numerical simulations show that the relative importance of flattening and rotation are functions of the galaxy mass: at low galaxy masses, adding flattening and rotation induces a galactic wind, thus lowering the X-ray luminosity; at high galaxy masses the angular momentum conservation keeps the central regions of rotating galaxies at low density, whereas in non-rotating models a denser and brighter atmosphere is formed. The same dependence from the galaxy mass is present in the effects of star formation (SF): in light galaxies SF contributes to increase the spread in Lx, while at high galaxy masses the halo X-ray properties are marginally sensitive to SF effects. In every case, the star formation rate at the present epoch quite agrees with observations, and the massive, cold gaseous discs are partially or completely consumed by SF on a time-scale of few Gyr, excluding the presence of young stellar discs at the present epoch.

Theoretical study of correlated systems using hybrid functionals

Diese Arbeit unterstreicht das Potential von Hybridfunktionalen (B3LYP) für die Untersuchung einer großen Bandbreite von Systemen. Durch die Einbeziehung der exakten Hartree-Fock Austauschenergie kann B3LYP für molekulare und kristalline Systeme eingesetzt werden. Zum Beispiel können stark korrelierte Systeme mit B3LYP erfolgreich erforscht werden. Die elektronische Struktur von PAHs wurde mit B3LYP Hybriddichtefunktionalen untersucht. Mit der ∆SCF-Methode wurden Elektronenbindungsenergien bestimmt, welche die mit UPS gewonnenen experimentellen Resultate bestätigen und ergänzen. Symmetrieeigenschaften der molekularen Orbitale wurden analysiert, um eine Zuordnung und Einschätzung der zugehörigen Signalstärke zu ermöglichen. Während σ-artige Orbitale nur schwer durch UPS-Messungen an dünnen Filmen detektiert werden können, bieten Rechnungen eine detaillierte Einsicht in die verborgenen Teile der Spektren.rnWeiterhin wurden π−π-Komplexe untersucht, welche von verschiedenen Donor- und Akzeptor-Molekülen gebildet werden. Die Moleküle basieren auf polyzyklischen, aromatischen Kohlenwasserstoffen. Für Ladungstransferkomplexe finden DFT Rechnungen ein Minimum in der Oberfläche der potentiellen Energie. Diese attraktive Wechselwirkung wird durch Coulombanziehung verursacht. Allerdings ist die Coulombanziehung nicht die stärkste Wechselwirkung in Ladungstransferkomplexen. Die Einbeziehung von van der Waals-Korrekturen verbessert den intermolekularen Abstand und die Bindungsenergie.rnEine Verkleinerung der intermolekularen Abstände führt zu einer großen Verschiebung der HOMO- und LUMO-Energie.rnAus der Klasse der kristallinen korrelierten Systeme wurden Rb4O6 und FeSe untersucht. Im Falle von Rb4O6 führen Ladungsordnung und Korrelationen zu einem isolierenden Grundzustand. Das hypothetische druckabhängige Phasendiagramm wurde untersucht. Eine Erhöhung des Drucks führt zu einer vergrößerten Bandlücke. Bei etwa 75 GPa wird die Bandbreite W größer als der Bandabstand U und das System nimmt einen homogen gemischt valenten Zustand mit teilweise besetzten π−π-Orbitalen an. Für Drücke ab 160 GPa wird W sehr viel größer als U und das System wird metallisch.rnIm Fall von FeSe finden wir eine korrelierte und isolierende Phase bei hohen Drücken, während das System bei niedrigen Drücken supraleitendes Verhalten zeigt. Die Berechnungen der Elektronenstruktur mit dem Hybridfunktional B3LYP führt zum korrekten halbleitenden Grundzustand in der NiAs- und MnP-Struktur von FeSe. Die Rolle der Korrelationen, der Stöchiometrie und der Nähe zum Magnetismus wird besprochen. Im Speziellen wird gezeigt, dass die Phase mit NiAs-Struktur starke lokale Korrelationen aufweist, was zu einem halbleitenden Zustand in einem weiten Druckbereich führt.

Application of Submerged Grouted Anchors in Sheet Pile Quay Walls

Sheet pile walls are one of the oldest earth retention systems utilized in civil engineering projects. They are used for various purposes; such as excavation support system, cofferdams, cut-off walls under dams, slope stabilization, waterfront structures, and flood walls. Sheet pile walls are one of the most common types of quay walls used in port construction. The worldwide increases in utilization of large ships for transportation have created an urgent need of deepening the seabed within port areas and consequently the rehabilitation of its wharfs. Several methods can be used to increase the load-carrying capacity of sheet-piling walls. The use of additional anchored tie rods grouted into the backfill soil and arranged along the exposed wall height is one of the most practical and appropriate solutions adopted for stabilization and rehabilitation of the existing quay wall. The Ravenna Port Authority initiated a project to deepen the harbor bottom at selected wharves. An extensive parametric study through the finite element program, PLAXIS 2D, version 2012 was carried out to investigate the enhancement of using submerged grouted anchors technique on the load response of sheet-piling quay wall. The influence of grout-ties area, length of grouted body, anchor inclination and anchor location were considered and evaluated due to the effect of different system parameters. Also a comparative study was conducted by Plaxis 2D and 3D program to investigate the behavior of these sheet pile quay walls in terms of horizontal displacements induced along the sheet pile wall and ground surface settlements as well as the anchor force and calculated factor of safety. Finally, a comprehensive study was carried out by using different constitutive models to simulate the mechanical behavior of the soil to investigate the effect of these two models (Mohr-Coulomb and Hardening Soil) on the behavior of these sheet pile quay walls.

Grundlegende Untersuchungen zum Verhalten ausgewählter Seltenerdelemente bei Explosionen und zur Anwendbarkeit spezieller Methoden zu deren Analyse

Bisher ist bei forensischen Untersuchungen von Explosionen die Rückverfolgung der verwendeten Sprengstoffe begrenzt, da das Material in aller Regel bei der Explosion zerstört wird. Die Rückverfolgung von Sprengstoffen soll mit Hilfe von Identifikations-Markierungssubstanzen erleichtert werden. Diese stellen einen einzigartigen Code dar, der auch nach einer Sprengung wiedergefunden und identifiziert werden kann. Die dem Code zugeordneten, eindeutigen Informationen können somit ausgelesen werden und liefern der Polizei bei der Aufklärung weitere Ansätze.rnZiel der vorliegenden Arbeit ist es, das Verhalten von ausgewählten Seltenerdelementen (SEE) bei Explosion zu untersuchen. Ein auf Lanthanoidphosphaten basierender Identifikations-Markierungsstoff bietet die Möglichkeit, verschiedene Lanthanoide innerhalb eines einzelnen Partikels zu kombinieren, wodurch eine Vielzahl von Codes generiert werden kann. Somit kann eine Veränderung der Ausgangszusammensetzung des Codes auch nach einer Explosion durch die Analyse eines einzelnen Partikels sehr gut nachvollzogen und somit die Eignung des Markierungsstoffes untersucht werden. Eine weitere Zielsetzung ist die Überprüfung der Anwendbarkeit der Massenspektrometrie mit induktiv gekoppeltem Plasma (ICP-MS) und Partikelanalyse mittels Rasterelektronenmikroskopie (REM) für die Analyse der versprengten Identifikations-Markierungssubstanzen. rnDie Ergebnisbetrachtungen der ICP-MS-Analyse und REM-Partikelanalyse deuten zusammenfassend auf eine Fraktionierung der untersuchten Lanthanoide oder deren Umsetzungsprodukte nach Explosion in Abhängigkeit ihrer thermischen Belastbarkeit. Die Befunde zeigen eine Anreicherung der Lanthanoide mit höherer Temperaturbeständigkeit in größeren Partikeln, was eine Anreicherung von Lanthanoiden mit niedrigerer Temperaturbeständigkeit in kleineren Partikeln impliziert. Dies lässt sich in Ansätzen durch einen Fraktionierungsprozess in Abhängigkeit der Temperaturstabilität der Lanthanoide oder deren Umsetzungsprodukten erklären. Die der Fraktionierung zugrunde liegenden Mechanismen und deren gegenseitige Beeinflussung bei einer Explosion konnten im Rahmen dieser Arbeit nicht abschließend geklärt werden.rnDie generelle Anwendbarkeit und unter Umständen notwendige, komplementäre Verwendung der zwei Methoden ICP-MS und REM-Partikelanalyse wird in dieser Arbeit gezeigt. Die ICP-MS stellt mit großer untersuchter Probenfläche und hoher Genauigkeit eine gute Methode zur Charakterisierung der Konzentrationsverhältnisse der untersuchten Lanthanoide dar. Die REM-Partikelanalyse hingegen ermöglicht im Falle von Kontamination der Proben mit anderen Lanthanoid-haltigen Partikeln eine eindeutige Differenzierung der Elementvergesellschaftung pro Partikel. Sie kann somit im Gegensatz zur ICP-MS Aufschluss über die Art und Zusammensetzung der Kontamination geben. rnInnerhalb der vorgenommenen Untersuchungen stellte die bei der ICP-MS angewandte Probennahmetechnik eine ideale Art der Probennahme dar. Bei anderen Oberflächen könnte diese jedoch in Folge der in verschiedenen Partikelgrößen resultierenden Fraktionierung zu systematisch verfälschten Ergebnissen führen. Um die generelle Anwendbarkeit der ICP-MS im Hinblick auf die Analyse versprengter Lanthanoide zu gewährleisten, sollte eine Durchführung weiterer Sprengungen auf unterschiedlichen Probenoberflächen erfolgen und gegebenenfalls weitere Probennahme-, Aufschluss- und Anreicherungsverfahren evaluiert werden.rn

Ionic liquids in a multiscale study

Trotz des hohen Interesse an Ionischen Flüssigkeiten wird das zielgerichtete Design und die Anwendung Ionischer Flüssigkeiten durch fehlendes grundlegendes Verständnis erschwert. Deshalb wurde die Balance der molekularen Wechselwirkungen in Ionischen Flüssigkeiten studiert, um die Eigenschaften dieser zu verstehen und die Kraftfeldentwicklung im Rahmen des Multiskalenansatzes zu systematisieren. Es wurden reine Imidazolium-basierte Ionische Flüssigkeiten, Mischungen mit kleinen Molekülen und eine protische Ionische Flüssigkeit mit ab-initio-Methoden, hauptsächlich Car-Parrinello-Molekulardynamik, untersucht. Weiterhin wurden Eigenschaften der Flüssigphase mit denen von Ionenpaaren verglichen.rnIm Fokus standen die molekularen elektrostatischen Eigenschaften und es wurde gezeigt, dass Coulomb-Wechselwirkungen zu einzigartigen Charakteristika führten. So waren die Ionen-Nettoladungen stets reduziert, die molekularen Dipolmomentverteilungen sehr breit, elektronische Polarisation war entscheidend. Die elektrostatischen Eigenschaften waren allgemein lokal auf molekularen Größen- und Zeitskalen und hingen stark von Phasenzustand und Zusammensetzung ab. Für andere molekulare Eigenschaften, wie der Neigung zu dispersiven Kontakten oder Wasserstoffbrücken, wurde gezeigt, dass sie einen entscheidenden Einfluss auf die Feinstruktur Ionischer Flüssigkeiten hatten. Das Gleichgewicht der Wechselwirkungen zeigte sich auch in Leistungsspektren, die sich aus den ab-initio-Molekulardynamiksimulationen ergaben. Diese boten einen neuen Weg für den Vergleich zum Experiment und für einen Einblick in die schnelle Dynamik Ionischer Flüssigkeiten.

Setup of a Penning trap for precision laser spectroscopy at HITRAP

Ion traps have been established as a powerful tool for ion cooling and laser spectroscopy experiments since a long time ago. SpecTrap, one of the precision experiments associated to the HITRAP facility at GSI, is implementing a Penning trap for studies of large bunches of externally produced highly charged ions. The extremely strong electric and magnetic fields that exist around the nuclei of heavy elements drastically change their electronic properties, such as energy level spacings and radiative lifetimes. The electrons can therefore serve as sensitive probes for nuclear properties such as size, magnetic moment and spatial distribution of charge and magnetization. The energies of forbidden fine and hyperfine structure transitions in such ions strongly depend on the nuclear charge and shift from the microwave domain into the optical domain. Thus, they become accessible for laser spectroscopy and its potentially high accuracy. A number of such measurements has been performed in storage rings and electron beam ion traps and yielded results with relative accuracies in the 10