Oberflächenplasmonen-artige Resonanzen auf metallischen Gittern

Die Geometrie einer Metalloberfläche bestimmt die Wechselwirkung zwischen der Oberflächenplasmonenresonanz und anderen Energieformen wie Photonen, anderen Oberflächenplasmonen und molekularen Anregungen. In dieser Arbeit wird der Einfluss dieser Effekte auf die optischen Eigenschaften von metallischen Reliefgittern untersucht.Aufgrund von Modellrechungen werden auf tiefen Gitterstrukturen Resonanzen erwartet, die sich als selbstgekoppelte Oberflächenplasmonen interpretieren lassen. Diese Resonanzen wurden aufgrund der Symmetrie des elektromagnetischen Feldes klassifiziert. Es wurden tiefe Gitter mit unterschiedlichen Profilformen hergestellt, die die experimentelle Beobachtung von drei gekoppelten Resonanzen erlauben. Variationen in der Tiefe und Asymmetrie der Gitter sowie in der experimentellen Geometrie zeigen die theoretisch vorausgesagten Effekte.Fluoreszenzfarbstoffe können mit elektromagnetischen Oberflächenresonanzen Energie austauschen. Die Lokalisierung des elektrischen Feldes von gekoppelten Resonanzen führt zu einer Ortsabhängigkeit der Kopplung, die bei 'freien' Oberflächenplasmonen nicht beobachtet wird. Ein theoretisches Modell ist in der Lage, die experimentellen Befunde weitgehend zu beschreiben. Die Ortsabhägigkeit der photochemischen Zerstörungsrate erlaubt die Anwendung dieses Effektes zur Messung von Diffusionsphänomenen in dünnen Filmen.Des weiteren wurde die Polarisationsabhängigkeit der Anregung von Oberflächenplasmonen wurde in konischer Reflexionsgeometrie und die Rolle der Oberflächenplasmonen in der thermisch induzierten Lichtemission untersucht. Ferner wurde eine vereinfachte Auswertungsroutine zur Anwendung von Gittern in der Untersuchung von dünnen dielektrischen Filmen entwickelt.

Characterisation of the cell-transistor coupling

Die Verbindung von elektrisch aktiven, lebenden Zellen zu extrazellulären Sensorsystemen eröffnet vielfälige Möglichkeiten im Bereich der Biosensorik. Die vorliegende Arbeit leistet einen Beitrag zum tieferen Verständnis der elektrischen Kopplungsmechanismen zwischen den biologischen und elektronischen Teilen solcher Hybridsysteme. Es wurden dazu drei Hauptbereiche bearbeitet:Ein System zur extrazellulären Signalableitung an lebenden Zellen bestehend aus einem Sensorchip, einem Vorverstärkerkopf und einem Hauptverstärker wurde weiterentwickelt.Als Sensoren wurden entweder Metallmikroelektroden-Chips mit 64 Kanälen oder Feldeffekt Transistoren-Chips mit 16 Kanälen (FET) eingesetzt. Es wurden zusätzlich spezielle FET Sensoren mit Rückseitenkontakten hergestellt und eingesetzt.Die elektrische Kopplung von einzelnen Nervenzellen der neuronalen Zell-Linien SH-SY5Y und TR14 oder primär kultivierten Neuronen aus dem Hirnstamm oder dem Hippocampus von embryonalen Ratten mit den extrazellulären Sensoren wurde untersucht. In der 'whole-cell' Patch-Clamp Technik wurden die Beiträge der spannungsgesteuerten Na+- und K+-Ionenkanäle zur extrazellulären Signalform identifiziert. Die Simulation der Signale mit einem Ersatzschaltkreis (Punkt-Kontakt Modell), der in PSPICE implementiert wurde, deutet auf eine starke Abhängigkeit der Signalformen in bezug auf Konzentrationsänderungen von Na+- und K+-Ionen im Volumenbereich zwischen Zelle und den ionensensitiven Transistoren hin. Ein empirisch erweitertes Punkt-Kontakt Modell wurde daraufhin vorgestellt.Im dritten Teil der Arbeit wurden Zellschichten von Kardiomyocyten embryonaler Ratten auf den extrazellulären Sensoren kultiviert. Die Eignung eines solchen Hybridsensors als Modellherz fuer das pharmazeutische Screeing wurde durch Messungen mit Herzstimulanzien und -relaktanzien bestätigt.

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.

Mode coupling for precise measurements of the electronic g-factor of hydrogen-like ions in Penning traps

The g-factor is a constant which connects the magnetic moment $vec{mu}$ of a charged particle, of charge q and mass m, with its angular momentum $vec{J}$. Thus, the magnetic moment can be writen $ vec{mu}_J=g_Jfrac{q}{2m}vec{J}$. The g-factor for a free particle of spin s=1/2 should take the value g=2. But due to quantum electro-dynamical effects it deviates from this value by a small amount, the so called g-factor anomaly $a_e$, which is of the order of $10^{-3}$ for the free electron. This deviation is even bigger if the electron is exposed to high electric fields. Therefore highly charged ions, where electric field strength gets values on the order of $10^{13}-10^{16}$V/cm at the position of the bound electron, are an interesting field of investigations to test QED-calculations. In previous experiments [H"aff00,Ver04] using a single hydrogen-like ion confined in a Penning trap an accuracy of few parts in $10^{-9}$ was obtained. In the present work a new method for precise measurement of magnetic the electronic g-factor of hydrogen-like ions is discussed. Due to the unavoidable magnetic field inhomogeneity in a Penning trap, a very important contribution to the systematic uncertainty in the previous measurements arose from the elevated energy of the ion required for the measurement of its motional frequencies. Then it was necessary to extrapolate the result to vanishing energies. In the new method the energy in the cyclotron degree of freedom is reduced to the minimum attainable energy. This method consist in measuring the reduced cyclotron frequency $nu_{+}$ indirectly by coupling the axial to the reduced cyclotron motion by irradiation of the radio frequency $nu_{coup}=nu_{+}-nu_{ax}+delta$ where $delta$ is, in principle, an unknown detuning that can be obtained from the knowledge of the coupling process. Then the only unknown parameter is the desired value of $nu_+$. As a test, a measurement with, for simplicity, artificially increased axial energy was performed yielding the result $g_{exp}=2.000~047~020~8(24)(44)$. This is in perfect agreement with both the theoretical result $g_{theo}=2.000~047~020~2(6)$ and the previous experimental result $g_{exp1}=2.000~047~025~4(15)(44).$ In the experimental results the second error-bar is due to the uncertainty in the accepted value for the electron's mass. Thus, with the new method a higher accuracy in the g-factor could lead by comparison to the theoretical value to an improved value of the electron's mass. [H"af00] H. H"affner et al., Phys. Rev. Lett. 85 (2000) 5308 [Ver04] J. Verd'u et al., Phys. Rev. Lett. 92 (2004) 093002-1

Production, radiochemical separation and chemical coupling of radioactive arsenic isotopes to synthesize radiopharmaceuticals for molecular imaging

Noninvasive molecular-imaging technologies are playing a keyrole in drug discovery, development and delivery. Positron Emission Tomography (PET) is such a molecular imaging technology and a powerful tool for the observation of various diseases. However, it is limited by the availability of agents with high selectivity to the target and a physical half-life of the used positron emitting nuclide which matches the biological half-life of the observed process. For the long lasting enrichment of antibodies in tumor tissue few suitable isotopes for PET imaging are currently available. The element arsenic provides a range of isotopes, which could be used for diagnosis and also for endoradiotherapy. This work describes the development of radiochemical separation procedures to separate arsenic isotopes in no-carrier-added (nca) purity from reactor or cyclotron irradiated targets, the development and evaluation of a labeling chemistry to attach these separated arsenic isotopes to monoclonal antibodies, the in vitro and in vivo evaluation of antibodies labeled with radioactive arsenic isotopes and the molecular imaging using small animal PET.

Study of electromagnetic reactions on light nuclei with the Lorentz Integral Transform method

In dieser Arbeit aus dem Bereich der Wenig-Nukleonen-Physik wird die neu entwickelte Methode der Lorentz Integral Transformation (LIT) auf die Untersuchung von Kernphotoabsorption und Elektronenstreuung an leichten Kernen angewendet. Die LIT-Methode ermoeglicht exakte Rechnungen durchzufuehren, ohne explizite Bestimmung der Endzustaende im Kontinuum. Das Problem wird auf die Loesung einer bindungzustandsaehnlichen Gleichung reduziert, bei der die Endzustandswechselwirkung vollstaendig beruecksichtigt wird. Die Loesung der LIT-Gleichung wird mit Hilfe einer Entwicklung nach hypersphaerischen harmonischen Funktionen durchgefuehrt, deren Konvergenz durch Anwendung einer effektiven Wechselwirkung im Rahmem des hypersphaerischen Formalismus (EIHH) beschleunigt wird. In dieser Arbeit wird die erste mikroskopische Berechnung des totalen Wirkungsquerschnittes fuer Photoabsorption unterhalb der Pionproduktionsschwelle an 6Li, 6He und 7Li vorgestellt. Die Rechnungen werden mit zentralen semirealistischen NN-Wechselwirkungen durchgefuehrt, die die Tensor Kraft teilweise simulieren, da die Bindungsenergien von Deuteron und von Drei-Teilchen-Kernen richtig reproduziert werden. Der Wirkungsquerschnitt fur Photoabsorption an 6Li zeigt nur eine Dipol-Riesenresonanz, waehrend 6He zwei unterschiedliche Piks aufweist, die dem Aufbruch vom Halo und vom Alpha-Core entsprechen. Der Vergleich mit experimentellen Daten zeigt, dass die Addition einer P-Wellen-Wechselwirkung die Uebereinstimmung wesentlich verbessert. Bei 7Li wird nur eine Dipol-Riesenresonanz gefunden, die gut mit den verfuegbaren experimentellen Daten uebereinstimmt. Bezueglich der Elektronenstreuung wird die Berechnung der longitudinalen und transversalen Antwortfunktionen von 4He im quasi-elastischen Bereich fuer mittlere Werte des Impulsuebertrages dargestellt. Fuer die Ladungs- und Stromoperatoren wird ein nichtrelativistisches Modell verwendet. Die Rechnungen sind mit semirealistischen Wechselwirkungen durchgefuert und ein eichinvarianter Strom wird durch die Einfuehrung eines Mesonaustauschstroms gewonnen. Die Wirkung des Zweiteilchenstroms auf die transversalen Antwortfunktionen wird untersucht. Vorlaeufige Ergebnisse werden gezeigt und mit den verfuegbaren experimentellen Daten verglichen.

Development of a multi-species method by GC-coupling with inductively coupled plasma isotope dilution mass spectrometry for the simultaneous determination of alkylated lead, mercury and tin compounds

An accurate and sensitive species-specific GC-ICP-IDMS (gas chromatography inductively coupled plasma isotope dilution mass spectrometry) method for the determination of trimethyllead and a multi-species-specific GC-ICP-IDMS method for the simultaneous determination of trimethyllead, methylmercury, and butyltins in biological and environmental samples were developed. They allow the determination of corresponding elemental species down to the low ng g-1 range. The developed synthesis scheme for the formation of isotopically labeled Me3206Pb+ can be used for future production of this spike. The novel extraction technique, stir bar sorptive extraction (SBSE), was applied for the first time in connection with species-specific isotope dilution GC-ICP-MS for the determination of trimethyllead, methylmercury and butyltins. The results were compared with liquid-liquid extraction. The developed methods were validated by the analysis of certified reference materials. The liquid-liquid extraction GC-ICP-IDMS method was applied to seafood samples purchased from a supermarket. The methylated lead fraction in these samples, correlated to total lead, varied in a broad range of 0.01-7.6 %. On the contrary, the fraction of methylmercury is much higher, normally in the range of 80-98 %. The highest methylmercury content of up to 12 µg g-1 has been determined in shark samples, an animal which is at the end of the marine food chain, whereas in other seafood samples a MeHg+ content of less than 0.2 µg g-1 was found. Butyltin species could only be determined in samples, where anthropogenic contaminations must be assumed. This explains the observed broad variation of the butylated tin fraction in the range of <0.3-49 % in different seafood samples. Because all isotope-labelled spike compounds, except trimethyllead, are commercially available, the developed multi-species-specific GC-ICP-IDMS method has a high potential in future for routine analysis.

Beurteilung der Eignung der NPEMFE-Methode (Natural Pulsed Electromagnetic Field of Earth) mit dem "Cereskop" in Rutschungen und in Locker- und Festgesteinen mit Spannungsänderungen im Mittel- und Hochgebirge

Spannungsumlagerungen in Mineralen und Gesteinen induzieren in geologisch aktiven Bereichen mikromechanische und seismische Prozesse, wodurch eine schwache natürliche elektromagnetische Strahlung im Niederfrequenzbereich emittiert wird. Die elektromagnetischen Emissionen von nichtleitenden Mineralen sind auf dielektrische Polarisation durch mehrere physikalische Effekte zurückzuführen. Eine gerichtete mechanische Spannung führt zu einer ebenso gerichteten elektromagnetischen Emission. Die Quellen der elektromagnetischen Emissionen sind bekannt, jedoch können sie noch nicht eindeutig den verschiedenen Prozessen in der Natur zugeordnet werden, weshalb im Folgenden von einem seismo-elektromagnetischen Phänomen (SEM) gesprochen wird. Mit der neuentwickelten NPEMFE-Methode (Natural Pulsed Electromagnetic Field of Earth) können die elektromagnetischen Impulse ohne Bodenkontakt registriert werden. Bereiche der Erdkruste mit Spannungsumlagerungen (z.B. tektonisch aktive Störungen, potenzielle Hangrutschungen, Erdfälle, Bergsenkungen, Firstschläge) können als Anomalie erkannt und abgegrenzt werden. Basierend auf dem heutigen Kenntnisstand dieser Prozesse wurden Hangrutschungen und Locker- und Festgesteine, in denen Spannungsumlagerungen stattfinden, mit einem neuentwickelten Messgerät, dem "Cereskop", im Mittelgebirgsraum (Rheinland-Pfalz, Deutschland) und im alpinen Raum (Vorarlberg, Österreich, und Fürstentum Liechtenstein) erkundet und die gewonnenen Messergebnisse mit klassischen Verfahren aus Ingenieurgeologie, Geotechnik und Geophysik in Bezug gesetzt. Unter Feldbedingungen zeigte sich großenteils eine gute Übereinstimmung zwischen den mit dem "Cereskop" erkundeten Anomalien und den mit den konventionellen Verfahren erkundeten Spannungszonen. Auf Grundlage der bisherigen Kenntnis und unter Einbeziehung von Mehrdeutigkeiten werden die Messergebnisse analysiert und kritisch beurteilt.

Numerical coupling of thermal-electric network models and energy-transport equations including optoelectronic semiconductor devices

In this work the numerical coupling of thermal and electric network models with model equations for optoelectronic semiconductor devices is presented. Modified nodal analysis (MNA) is applied to model electric networks. Thermal effects are modeled by an accompanying thermal network. Semiconductor devices are modeled by the energy-transport model, that allows for thermal effects. The energy-transport model is expandend to a model for optoelectronic semiconductor devices. The temperature of the crystal lattice of the semiconductor devices is modeled by the heat flow eqaution. The corresponding heat source term is derived under thermodynamical and phenomenological considerations of energy fluxes. The energy-transport model is coupled directly into the network equations and the heat flow equation for the lattice temperature is coupled directly into the accompanying thermal network. The coupled thermal-electric network-device model results in a system of partial differential-algebraic equations (PDAE). Numerical examples are presented for the coupling of network- and one-dimensional semiconductor equations. Hybridized mixed finite elements are applied for the space discretization of the semiconductor equations. Backward difference formluas are applied for time discretization. Thus, positivity of charge carrier densities and continuity of the current density is guaranteed even for the coupled model.

Sampling methods for low-frequency electromagnetic imaging

For the detection of hidden objects by low-frequency electromagnetic imaging the Linear Sampling Method works remarkably well despite the fact that the rigorous mathematical justification is still incomplete. In this work, we give an explanation for this good performance by showing that in the low-frequency limit the measurement operator fulfills the assumptions for the fully justified variant of the Linear Sampling Method, the so-called Factorization Method. We also show how the method has to be modified in the physically relevant case of electromagnetic imaging with divergence-free currents. We present numerical results to illustrate our findings, and to show that similar performance can be expected for the case of conducting objects and layered backgrounds.

A sampling method for detecting buried objects using electromagnetic scattering

We consider a simple (but fully three-dimensional) mathematical model for the electromagnetic exploration of buried, perfect electrically conducting objects within the soil underground. Moving an electric device parallel to the ground at constant height in order to generate a magnetic field, we measure the induced magnetic field within the device, and factor the underlying mathematics into a product of three operations which correspond to the primary excitation, some kind of reflection on the surface of the buried object(s) and the corresponding secondary excitation, respectively. Using this factorization we are able to give a justification of the so-called sampling method from inverse scattering theory for this particular set-up.

Investigation of the dipole response of nickel isotopes in the presence of a high-frequency electromagnetic field

The electric dipole response of neutron-rich nickel isotopes has been investigated using the LAND setup at GSI in Darmstadt (Germany). Relativistic secondary beams of 56−57Ni and 67−72Ni at approximately 500 AMeV have been generated using projectile fragmentation of stable ions on a 4 g/cm2 Be target and subsequent separation in the magnetic dipole fields of the FRagment Separator (FRS). After reaching the LAND setup in Cave C, the radioactive ions were excited electromagnetically in the electric field of a Pb target. The decay products have been measured in inverse kinematics using various detectors. Neutron-rich 67−69Ni isotopes decay by the emission of neutrons, which are detected in the LAND detector. The present analysis concentrates on the (gamma,n) and (gamma,2n) channels in these nuclei, since the proton and three-neutron thresholds are unlikely to be reached considering the virtual photon spectrum for nickel ions at 500 AMeV. A measurement of the stable 58Ni isotope is used as a benchmark to check the accuracy of the present results with previously published data. The measured (gamma,n) and (gamma,np) channels are compared with an inclusive photoneutron measurement by Fultz and coworkers, which are consistent within the respective errors. The measured excitation energy distributions of 67−69Ni contain a large portion of the Giant Dipole Resonance (GDR) strength predicted by the Thomas-Reiche-Kuhn energy-weighted sum rule, as well as a significant amount of low-lying E1 strength, that cannot be attributed to the GDR alone. The GDR distribution parameters are calculated using well-established semi-empirical systematic models, providing the peak energies and widths. The GDR strength is extracted from the chi-square minimization of the model GDR to the measured data of the (gamma,2n) channel, thereby excluding any influence of eventual low-lying strength. The subtraction of the obtained GDR distribution from the total measured E1 strength provides the low-lying E1 strength distribution, which is attributed to the Pygmy Dipole Resonance (PDR). The extraction of the peak energy, width and strength is performed using a Gaussian function. The minimization of trial Gaussian distributions to the data does not converge towards a sharp minimum. Therefore, the results are presented by a chi-square distribution as a function of all three Gaussian parameters. Various predictions of PDR distributions exist, as well as a recent measurement of the 68Ni pygmy dipole-resonance obtained by virtual photon scattering, to which the present pygmy dipole-resonance distribution is also compared.

Advanced schemes for surface plasmon resonance and plasmon-enhanced fluorescence biosensors

Advanced optical biosensor platforms exploiting long range surface plasmons (LRSPs) and responsive N-isopropylacrylamide (NIPAAm) hydrogel binding matrix for the detection of protein and bacterial pathogen analytes were carried out. LRSPs are optical waves that originate from coupling of surface plasmons on the opposite sites of a thin metallic film embedded between two dielectrics with similar refractive indices. LRSPs exhibit orders of magnitude lower damping and more extended profile of field compared to regular surface plasmons (SPs). Their excitation is accompanied with narrow resonance and provides stronger enhancement of electromagnetic field intensity that can advance the sensitivity of surface plasmon resonance (SPR) and surface plasmon-enhanced fluorescence spectroscopy (SPFS) biosensors. Firstly, we investigated thin gold layers deposited on fluoropolymer surface for the excitation of LRSPs. The study indicates that the morphological, optical and electrical properties of gold film can be changed by the surface energy of fluoropolymer and affect the performance of a SPFS biosensor. A photo-crosslinkable NIPAAm hydrogel was grafted to the sensor surface in order to serve as a binding matrix. It was modified with bio-recognition elements (BREs) via amine coupling chemistry and offered the advantage of large binding capacity, stimuli responsive properties and good biocompatibility. Through experimental observations supported by numerical simulations describing diffusion mass transfer and affinity binding of target molecules in the hydrogel, the hydrogel binding matrix thickness, concentration of BREs and the profile of the probing evanescent field was optimized. Hydrogel with a up to micrometer thickness was shown to support additional hydrogel optical waveguide (HOW) mode which was employed for probing affinity binding events in the gel by means of refractometric and fluorescence measurements. These schemes allow to reach limits of detection (LODs) at picomolar and femtomolar levels, respectively. Besides hydrogel based experiments for detection of molecular analytes, long range surface plasmon-enhanced fluorescence spectroscopy (LRSP-FS) was employed for detection of bacterial pathogens. The influence of capture efficiency of bacteria on surfaces and the profile of the probing field on sensor response were investigated. The potential of LRSP-FS with extended evanescent field is demonstrated for detection of pathogenic E. coli O157:H7 on sandwich immunoassays . LOD as low as 6 cfu mL-1 with a detection time of 40 minutes was achieved.rn

Coarse-graining and quantum-classical adaptive coupling in soft matter

Die vorliegende Arbeit untersucht den Zusammenhang zwischen Skalen in Systemen weicher Materie, der für Multiskalen-Simulationen eine wichtige Rolle spielt. Zu diesem Zweck wurde eine Methode entwickelt, die die Approximation der Separierbarkeit von Variablen für die Molekulardynamik und ähnliche Anwendungen bewertet. Der zweite und größere Teil dieser Arbeit beschäftigt sich mit der konzeptionellen und technischen Erweiterung des Adaptive Resolution Scheme'' (AdResS), einer Methode zur gleichzeitigen Simulation von Systemen mit mehreren Auflösungsebenen. Diese Methode wurde auf Systeme erweitert, in denen klassische und quantenmechanische Effekte eine Rolle spielen.rnrnDie oben genannte erste Methode benötigt nur die analytische Form der Potentiale, wie sie die meisten Molekulardynamik-Programme zur Verfügung stellen. Die Anwendung der Methode auf ein spezielles Problem gibt bei erfolgreichem Ausgang einen numerischen Hinweis auf die Gültigkeit der Variablenseparation. Bei nicht erfolgreichem Ausgang garantiert sie, dass keine Separation der Variablen möglich ist. Die Methode wird exemplarisch auf ein zweiatomiges Molekül auf einer Oberfläche und für die zweidimensionale Version des Rotational Isomer State (RIS) Modells einer Polymerkette angewandt.rnrnDer zweite Teil der Arbeit behandelt die Entwicklung eines Algorithmus zur adaptiven Simulation von Systemen, in denen Quanteneffekte berücksichtigt werden. Die Quantennatur von Atomen wird dabei in der Pfadintegral-Methode durch einen klassischen Polymerring repräsentiert. Die adaptive Pfadintegral-Methode wird zunächst für einatomige Flüssigkeiten und tetraedrische Moleküle unter normalen thermodynamischen Bedingungen getestet. Schließlich wird die Stabilität der Methode durch ihre Anwendung auf flüssigen para-Wasserstoff bei niedrigen Temperaturen geprüft.

Mode-coupling theory: generalizations, high dimensions and microscopic dynamics

This work contains several applications of the mode-coupling theory (MCT) and is separated into three parts. In the first part we investigate the liquid-glass transition of hard spheres for dimensions d→∞ analytically and numerically up to d=800 in the framework of MCT. We find that the critical packing fraction ϕc(d) scales as d²2^(-d), which is larger than the Kauzmann packing fraction ϕK(d) found by a small-cage expansion by Parisi and Zamponi [J. Stat. Mech.: Theory Exp. 2006, P03017 (2006)]. The scaling of the critical packing fraction is different from the relation ϕc(d)∼d2^(-d) found earlier by Kirkpatrick and Wolynes [Phys. Rev. A 35, 3072 (1987)]. This is due to the fact that the k dependence of the critical collective and self nonergodicity parameters fc(k;d) and fcs(k;d) was assumed to be Gaussian in the previous theories. We show that in MCT this is not the case. Instead fc(k;d) and fcs(k;d), which become identical in the limit d→∞, converge to a non-Gaussian master function on the scale k∼d^(3/2). We find that the numerically determined value for the exponent parameter λ and therefore also the critical exponents a and b depend on the dimension d, even at the largest evaluated dimension d=800. In the second part we compare the results of a molecular-dynamics simulation of liquid Lennard-Jones argon far away from the glass transition [D. Levesque, L. Verlet, and J. Kurkijärvi, Phys. Rev. A 7, 1690 (1973)] with MCT. We show that the agreement between theory and computer simulation can be improved by taking binary collisions into account [L. Sjögren, Phys. Rev. A 22, 2866 (1980)]. We find that an empiric prefactor of the memory function of the original MCT equations leads to similar results. In the third part we derive the equations for a mode-coupling theory for the spherical components of the stress tensor. Unfortunately it turns out that they are too complex to be solved numerically.