Novel simulation methods for Coulomb and hydrodynamic interactions

This thesis presents new methods to simulate systems with hydrodynamic and electrostatic interactions. Part 1 is devoted to computer simulations of Brownian particles with hydrodynamic interactions. The main influence of the solvent on the dynamics of Brownian particles is that it mediates hydrodynamic interactions. In the method, this is simulated by numerical solution of the Navier--Stokes equation on a lattice. To this end, the Lattice--Boltzmann method is used, namely its D3Q19 version. This model is capable to simulate compressible flow. It gives us the advantage to treat dense systems, in particular away from thermal equilibrium. The Lattice--Boltzmann equation is coupled to the particles via a friction force. In addition to this force, acting on {it point} particles, we construct another coupling force, which comes from the pressure tensor. The coupling is purely local, i.~e. the algorithm scales linearly with the total number of particles. In order to be able to map the physical properties of the Lattice--Boltzmann fluid onto a Molecular Dynamics (MD) fluid, the case of an almost incompressible flow is considered. The Fluctuation--Dissipation theorem for the hybrid coupling is analyzed, and a geometric interpretation of the friction coefficient in terms of a Stokes radius is given. Part 2 is devoted to the simulation of charged particles. We present a novel method for obtaining Coulomb interactions as the potential of mean force between charges which are dynamically coupled to a local electromagnetic field. This algorithm scales linearly, too. We focus on the Molecular Dynamics version of the method and show that it is intimately related to the Car--Parrinello approach, while being equivalent to solving Maxwell's equations with freely adjustable speed of light. The Lagrangian formulation of the coupled particles--fields system is derived. The quasi--Hamiltonian dynamics of the system is studied in great detail. For implementation on the computer, the equations of motion are discretized with respect to both space and time. The discretization of the electromagnetic fields on a lattice, as well as the interpolation of the particle charges on the lattice is given. The algorithm is as local as possible: Only nearest neighbors sites of the lattice are interacting with a charged particle. Unphysical self--energies arise as a result of the lattice interpolation of charges, and are corrected by a subtraction scheme based on the exact lattice Green's function. The method allows easy parallelization using standard domain decomposition. Some benchmarking results of the algorithm are presented and discussed.

108 Sn studied with intermediate-energy Coulomb excitation

In dieser Arbeit wurde der instabile, Neutronenarme Kern 108Sn mit Hilfe der Coulomb-Anregung bei intermediaeren Energien in inverser Kinematik studiert. Diese Methode wurde bisher zur Untersuchung der ersten angeregten 2+ Zustaende und deren E2 Zerfallsraten in Kernen mit Kernladungszahl Z< 30 angewendet. 108Sn ist somit der Kern mit der groeßten Kernladungszahl, bei dem diese Studien bisher stattfanden. Das Ziel dieses Experiments war die Messung der unbekannten reduzierten Uebergangswahrscheinlichkeit B(E2,0+ -> 2+). Der B(E2)-Wert von 0.230(57) e2b2 wurde relativ zu dem bekannten Wert des Isotops 112Sn bestimmt. Das Experiment wurde an der GSI Darmstadt mit Hilfe des RISING Detektors und des Fragmentseperators (FRS) durchgefuehrt. Sekundaere Strahlen (108Sn, 112Sn) mit einer Energie von ca. 150 MeV pro Nukleon wurden auf ein 386 mg/cm2 dickes 197Au Target geschossen. Die Projektilfragmente wurden mit Hilfe des Fragmentseparators selektiert und identifiziert. Zur Selektion des Reaktionskanals und zur Bestimmung des Winkels der gestreuten Fragmente wurde das Teilchenteleskop CATE, das sich hinter dem Target befand, verwendet. Gammastrahlung, die in Koinzidenz mit den Projektilrestkernen emittiert wurde, wurde in den Germanium-Cluster Detektoren des RISING Detektors nachgewiesen. Der gemessene B(E2,0+ -> 2+)-Wert von 108Sn ist in Uebereinstimmung mit neueren Schalenmodellrechnungen, die auf realistischen effektiven Wechselwirkungen basieren und im Rahmen eines verallgemeinerten Seniorit¨ats-Schemas erklaert werden.

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.

Interacting Bose-Fermi mixtures in optical lattices

In this thesis, we investigate mixtures of quantum degenerate Bose and Fermi gases of neutral atoms in threedimensional optical lattices. Feshbach resonances allow to control interspecies interactions in these systems precisely, by preparing suitable combinations of internal atomic states and applying external magnetic fields. This way, the system behaviour can be tuned continuously from mutual transparency to strongly interacting correlated phases, up to the stability boundary.rnThe starting point for these investigations is the spin-polarized fermionic band insulator. The properties of this non-interacting system are fully determined by the Pauli exclusion principle for the occupation of states in the lattice. A striking demonstration of the latter can be found in the antibunching of the density-density correlation of atoms released from the lattice. If bosonic atoms are added to this system, isolated heteronuclear molecules can be formed on the lattice sites via radio-frequency stimulation. The efficiency of this process hints at a modification of the atom number distribution over the lattice caused by interspecies interaction.rnIn the following, we investigate systems with tunable interspecies interaction. To this end, a method is developed which allows to assess the various contributions to the system Hamiltonian both qualitatively and quantitatively by following the quantum phase diffusion of the bosonic matter wave.rnBesides a modification of occupation number statistics, these measurements show a significant renormalization of the bosonic Hubbard parameters. The final part of the thesis considers the implications of this renormalization effect on the many particle physics in the mixture. Here, we demonstrate how the quantum phase transition from a bosonic superfluid to a Mott insulator state is shifted towards considerably shallower lattices due to renormalization.

Temperatur- und ortsabhängige Zustandsdichte des organischen Supraleiters Kappa-(BEDT-TTF) 2 Cu[N(CN) 2]Br

Im Rahmen dieser Arbeit wurde die temperatur- und ortsabhängige Zustandsdichte des organischen Supraleiters kappa-(BEDT-TTF)2Cu[N(CN)2]Br mit Rastertunnelspektroskopie bei tiefen Temperaturen untersucht.rnZusätzlich zur bereits bekannten supraleitenden Energielücke wird dabei eine logarithmische Unterdrückung der Zustandsdichte an der Fermikante beobachtet, die auch oberhalb der kritischen Temperatur erhalten bleibt. In der vorliegenden Arbeit wird gezeigt, dass sich dieses Verhalten durch ein für ungeordnete elektronische Systeme entwickeltes Modell unter Berücksichtigung von Coulomb-Wechselwirkungen beschreiben lässt. Die daraus resultierenden Fluktuationen der elektronischen Struktur führen zu einer Verbreiterung der gemessenen supraleitenden Energielücke, die sich durch sehr kleine Kohärenzmaxima im entsprechenden Quasiteilchenanregungsspektrum äußert. Dieses Verhalten wurde bereits beobachtet, konnte jedoch bisher nicht erklärt werden. Die theoretische Beschreibung der logarithmischen Unterdrückung trägt somit zusätzlich zum Verständnis des supraleitenden Beitrags bei, sodass die gesamte Zustandsdichte vollständig beschrieben werden kann. Die Analyse der gemessenen supraleitenden Energielücke wurde für verschiedene Symmetrien des Ordnungsparameters durchgeführt, wobei die beste Übereinstimmung für die Annahme einer d-wellenartigen Symmetrie mit zwei unterschiedlich stark ausgeprägten Energielücken gefunden wurde. Der Paarbildungsmechanismus, der zur Bindung zweier Elektronen zu einem Cooper-Paar führt, kann mit einer $d$-wellenartigen Symmetrie nicht durch die in konventionellen Supraleitern gefundene Elektron-Phonon-Kopplung erklärt werden. Stattdessen wird in Analogie zur Hochtemperatur-Supraleitung eine durch antiferromagnetische Spin-Wechselwirkungen induzierte Kopplung der Elektronen vermutet. Dies wird zum einen durch die oberhalb der kritischen Temperatur auftretende, zweite Energielücke und zum anderen durch die zwischen 4,66 und 5,28 liegende Kopplungsstärke 2Delta/(kB Tc) unterstützt, die deutlich größer als für konventionelle Supraleiter mit Elektron-Phonon-Kopplung ist.