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.

Systematic studies of protein immobilization by surface plasmon field-enhanced fluorescence spectroscopy

The research interest of this study is to investigate surface immobilization strategies for proteins and other biomolecules by the surface plasmon field-enhanced fluorescence spectroscopy (SPFS) technique. The recrystallization features of the S-layer proteins and the possibility of combining the S-layer lattice arrays with other functional molecules make this protein a prime candidate for supramolecular architectures. The recrystallization behavior on gold or on the secondary cell wall polymer (SCWP) was recorded by SPR. The optical thicknesses and surface densities for different protein layers were calculated. In DNA hybridization tests performed in order to discriminate different mismatches, recombinant S-layer-streptavidin fusion protein matrices showed their potential for new microarrays. Moreover, SCWPs coated gold chips, covered with a controlled and oriented assembly of S-layer fusion proteins, represent an even more sensitive fluorescence testing platform. Additionally, S-layer fusion proteins as the matrix for LHCII immobilization strongly demonstrate superiority over routine approaches, proving the possibility of utilizing them as a new strategy for biomolecular coupling. In the study of the SPFS hCG immunoassay, the biophysical and immunological characteristics of this glycoprotein hormone were presented first. After the investigation of the effect of the biotin thiol dilution on the coupling efficiently, the interfacial binding model including the appropriate binary SAM structure and the versatile streptavidin-biotin interaction was chosen as the basic supramolecular architecture for the fabrication of a SPFS-based immunoassay. Next, the affinity characteristics between different antibodies and hCG were measured via an equilibrium binding analysis, which is the first example for the titration of such a high affinity interaction by SPFS. The results agree very well with the constants derived from the literature. Finally, a sandwich assay and a competitive assay were selected as templates for SPFS-based hCG detection, and an excellent LOD of 0.15 mIU/ml was attained via the “one step” sandwich method. Such high sensitivity not only fulfills clinical requirements, but is also better than most other biosensors. Fully understanding how LHCII complexes transfer the sunlight energy directionally and efficiently to the reaction center is potentially useful for constructing biomimetic devices as solar cells. After the introduction of the structural and the spectroscopic features of LHCII, different surface immobilization strategies of LHCII were summarized next. Among them the strategy based on the His-tag and the immobilized metal (ion) affinity chromatography (IMAC) technique were of great interest and resulted in different kinds of home-fabricated His-tag chelating chips. Their substantial protein coupling capacity, maintenance of high biological activity and a remarkably repeatable binding ability on the same chip after regeneration was demonstrated. Moreover, different parameters related to the stability of surface coupled reconstituted complexes, including sucrose, detergent, lipid, oligomerization, temperature and circulation rate, were evaluated in order to standardize the most effective immobilization conditions. In addition, partial lipid bilayers obtained from LHCII contained proteo-liposomes fusion on the surface were observed by the QCM technique. Finally, the inter-complex energy transfer between neighboring LHCIIs on a gold protected silver surface by excitation with a blue laser (λ = 473nm) was recorded for the first time, and the factors influencing the energy transfer efficiency were evaluated.

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.

Assembly and characterization of supramolecular architectures for biosensor applications

The research has included the efforts in designing, assembling and structurally and functionally characterizing supramolecular biofunctional architectures for optical biosensing applications. In the first part of the study, a class of interfaces based on the biotin-NeutrAvidin binding matrix for the quantitative control of enzyme surface coverage and activity was developed. Genetically modified ß-lactamase was chosen as a model enzyme and attached to five different types of NeutrAvidin-functionalized chip surfaces through a biotinylated spacer. All matrices are suitable for achieving a controlled enzyme surface density. Data obtained by SPR are in excellent agreement with those derived from optical waveguide measurements. Among the various protein-binding strategies investigated in this study, it was found that stiffness and order between alkanethiol-based SAMs and PEGylated surfaces are very important. Matrix D based on a Nb2O5 coating showed a satisfactory regeneration possibility. The surface-immobilized enzymes were found to be stable and sufficiently active enough for a catalytic activity assay. Many factors, such as the steric crowding effect of surface-attached enzymes, the electrostatic interaction between the negatively charged substrate (Nitrocefin) and the polycationic PLL-g-PEG/PEG-Biotin polymer, mass transport effect, and enzyme orientation, are shown to influence the kinetic parameters of catalytic analysis. Furthermore, a home-built Surface Plasmon Resonance Spectrometer of SPR and a commercial miniature Fiber Optic Absorbance Spectrometer (FOAS), served as a combination set-up for affinity and catalytic biosensor, respectively. The parallel measurements offer the opportunity of on-line activity detection of surface attached enzymes. The immobilized enzyme does not have to be in contact with the catalytic biosensor. The SPR chip can easily be cleaned and used for recycling. Additionally, with regard to the application of FOAS, the integrated SPR technique allows for the quantitative control of the surface density of the enzyme, which is highly relevant for the enzymatic activity. Finally, the miniaturized portable FOAS devices can easily be combined as an add-on device with many other in situ interfacial detection techniques, such as optical waveguide lightmode spectroscopy (OWLS), the quartz crystal microbalance (QCM) measurements, or impedance spectroscopy (IS). Surface plasmon field-enhanced fluorescence spectroscopy (SPFS) allows for an absolute determination of intrinsic rate constants describing the true parameters that control interfacial hybridization. Thus it also allows for a study of the difference of the surface coupling influences between OMCVD gold particles and planar metal films presented in the second part. The multilayer growth process was found to proceed similarly to the way it occurs on planar metal substrates. In contrast to planar bulk metal surfaces, metal colloids exhibit a narrow UV-vis absorption band. This absorption band is observed if the incident photon frequency is resonant with the collective oscillation of the conduction electrons and is known as the localized surface plasmon resonance (LSPR). LSPR excitation results in extremely large molar extinction coefficients, which are due to a combination of both absorption and scattering. When considering metal-enhanced fluorescence we expect the absorption to cause quenching and the scattering to cause enhancement. Our further study will focus on the developing of a detection platform with larger gold particles, which will display a dominant scattering component and enhance the fluorescence signal. Furthermore, the results of sequence-specific detection of DNA hybridization based on OMCVD gold particles provide an excellent application potential for this kind of cheap, simple, and mild preparation protocol applied in this gold fabrication method. In the final chapter, SPFS was used for the in-depth characterizations of the conformational changes of commercial carboxymethyl dextran (CMD) substrate induced by pH and ionic strength variations were studied using surface plasmon resonance spectroscopy. The pH response of CMD is due to the changes in the electrostatics of the system between its protonated and deprotonated forms, while the ionic strength response is attributed from the charge screening effect of the cations that shield the charge of the carboxyl groups and prevent an efficient electrostatic repulsion. Additional studies were performed using SPFS with the aim of fluorophore labeling the carboxymethyl groups. CMD matrices showed typical pH and ionic strength responses, such as high pH and low ionic strength swelling. Furthermore, the effects of the surface charge and the crosslink density of the CMD matrix on the extent of stimuli responses were investigated. The swelling/collapse ratio decreased with decreasing surface concentration of the carboxyl groups and increasing crosslink density. The study of the CMD responses to external and internal variables will provide valuable background information for practical applications.

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.

Grating coupled surface plasmon enhanced fluorescence spectroscopy

Over the last three decades, sensors based on the phenomenon of surface plasmon resonance have proven particularly suitable for real time thin film characterization, gas detection, biomolecular interaction examination and to supplement electrochemical methods. Systems based on prism coupling have been combined with fluorescence detection under the name of surface plasmon fluorescence spectroscopy to increase sensitivity even further. Alternatively, metal gratings can be employed to match photons for plasmon resonance. The real time monitoring of binding reactions not yet been reported in the combination of fluorescence detection and grating coupling. Grating-based systems promise more competitive products, because of reduced operating costs, and offer benefits for device engineering. This thesis is comprised of a comprehensive study of the suitability of grating coupling for fluorescence based analyte detection. Fundamental properties of grating coupled surface plasmon fluorescence spectroscopy are described, as well as issues related to the commercial realization of the method. Several new experimental techniques are introduced and demonstrated in order to optimize performance in certain areas and improve upon capabilities in respect to prism-based systems. Holographically fabricated gratings are characterized by atomic force microscopy and optical methods, aided by simulations and profile parameters responsible for efficient coupling are analyzed. The directional emission of fluorophores immobilized on a grating surface is studied in detail, including the magnitude and geometry of the fluorescence emission pattern for different grating constants and polarizations. Additionally, the separation between the minimum of the reflected intensity and the maximum fluorescence excitation position is examined. One of the key requirements for the commercial feasibility of grating coupling is the cheap and faithful mass production of disposable samples from a given master grating. The replication of gratings is demonstrated by a simple hot embossing method with good reproducibility to address this matter. The in-situ fluorescence detection of analyte immobilization and affinity measurements using grating coupling are described for the first time. The physical factors related to the sensitivity of the technique are assessed and the lower limit of detection of the technique is determined for an exemplary assay. Particular attention is paid to the contribution of bulk fluorophores to the total signal in terms of magnitude and polarization of incident and emitted light. Emission from the bulk can be a limiting factor for experiments with certain assay formats. For that reason, a novel optical method, based on the modulation of both polarization and intensity of the incident beam, is introduced and demonstrated to be capable of eliminating this contribution.

Speciation of Pu(III) in the environmental system humic substances-groundwater-kaolinite

For the safety assessments of nuclear waste repositories, the possible migration of the radiotoxic waste into environment must be considered. Since plutonium is the major contribution at the radiotoxicity of spent nuclear waste, it requires special care with respect to its mobilization into the groundwater. Plutonium has one of the most complicated chemistry of all elements. It can coexist in 4 oxidation states parallel in one solution. In this work is shown that in the presence of humic substances it is reduced to the Pu(III) and Pu(IV). This work has the focus on the interaction of Pu(III) with natural occurring compounds (humic substances and clay minerals bzw. Kaolinite), while Pu(IV) was studied in a parallel doctoral work by Banik (in preparation). As plutonium is expected under extreme low concentrations in the environment, very sensitive methods are needed to monitor its presence and for its speciation. Resonance ionization mass spectrometry (RIMS), was used for determining the concentration of Pu in environmental samples, with a detection limit of 106- 107 atoms. For the speciation of plutonium CE-ICP-MS was routinely used to monitor the behaviour of Pu in the presence of humic substances. In order to reduce the detection limits of the speciation methods, the coupling of CE to RIMS was proposed. The first steps have shown that this can be a powerful tool for studies of pu under environmental conditions. Further, the first steps in the coupling of two parallel working detectors (DAD and ICP_MS ) to CE was performed, for the enabling a precise study of the complexation constants of plutonium with humic substances. The redox stabilization of Pu(III) was studied and it was determined that NH2OHHCl can maintain Pu(III) in the reduced form up to pH 5.5 – 6. The complexation constants of Pu(III) with Aldrich humic acid (AHA) were determined at pH 3 and 4. the logß = 6.2 – 6.8 found for these experiments was comparable with the literature. The sorption of Pu(III) onto kaolinite was studied in batch experiments and it was determine dthat the pH edge was at pH ~ 5.5. The speciation of plutonium on the surface of kaolinite was studied by EXAFS/XANES. It was determined that the sorbed species was Pu(IV). The influence of AHA on the sorption of Pu(III) onto kaolinite was also investigated. It was determined that at pH < 5 the adsorption is enhanced by the presence of AHA (25 mg/L), while at pH > 6 the adsorption is strongly impaired (depending also on the adding sequence of the components), leading to a mobilization of plutonium in solution.

Higher-order calculations in manifestly Lorentz-invariant baryon chiral perturbation theory

This thesis is concerned with calculations in manifestly Lorentz-invariant baryon chiral perturbation theory beyond order D=4. We investigate two different methods. The first approach consists of the inclusion of additional particles besides pions and nucleons as explicit degrees of freedom. This results in the resummation of an infinite number of higher-order terms which contribute to higher-order low-energy constants in the standard formulation. In this thesis the nucleon axial, induced pseudoscalar, and pion-nucleon form factors are investigated. They are first calculated in the standard approach up to order D=4. Next, the inclusion of the axial-vector meson a_1(1260) is considered. We find three diagrams with an axial-vector meson which are relevant to the form factors. Due to the applied renormalization scheme, however, the contributions of the two loop diagrams vanish and only a tree diagram contributes explicitly. The appearing coupling constant is fitted to experimental data of the axial form factor. The inclusion of the axial-vector meson results in an improved description of the axial form factor for higher values of momentum transfer. The contributions to the induced pseudoscalar form factor, however, are negligible for the considered momentum transfer, and the axial-vector meson does not contribute to the pion-nucleon form factor. The second method consists in the explicit calculation of higher-order diagrams. This thesis describes the applied renormalization scheme and shows that all symmetries and the power counting are preserved. As an application we determine the nucleon mass up to order D=6 which includes the evaluation of two-loop diagrams. This is the first complete calculation in manifestly Lorentz-invariant baryon chiral perturbation theory at the two-loop level. The numerical contributions of the terms of order D=5 and D=6 are estimated, and we investigate their pion-mass dependence. Furthermore, the higher-order terms of the nucleon sigma term are determined with the help of the Feynman-Hellmann theorem.

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.

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.

Thin films of polythiophene: Linear and nonlinear optical characterization

Wir haben die linearen und nichtlinearen optischen Eigenschaften von dünnen Schichten und planaren Wellenleitern aus mehreren konjugierten Polymeren (MEH-PPV und P3AT) und Polymeren mit -Elektronen Systemen in der Seitenkette (PVK und PS) untersucht und verglichen. PVK und PS haben relativ kleine Werte des nichtlinearen Brechungsindex n2 bei 532 nm, nämlich (1,2 ± 0,5)10-14 cm2/W und (2,6 ± 0,5) 10-14 cm2/W.rnWir haben die linearen optischen Konstanten von mehreren P3ATs untersucht, insbesondere den Einfluss der Regioregularität und Kettenlänge der Alkylsubstituenten. Wir haben das am besten geeignete Polymere für Wellenleiter Anwendungen identifiziert, welches P3BT-ra genannt ist. Wir haben die linearen optischen Eigenschaften dünner Schichten des P3BT-ra untersucht, die mit Spincoating aus verschiedenen Lösungsmitteln mit unterschiedlichen Siedetemperaturen präparieret wurden. Wir haben festgestellt, dass P3BT-ra Filme aus Toluol-Lösungen die am besten geeigneten Wellenleiter für die intensitätsabhängigen Prismen-Kopplungs Experimente sind, weil diese geringe Wellenleiterdämpfungsverluste bei = 1064 nm haben. rnWir haben die Dispersionen des Wellenleiterdämfungsverlustes gw, des nichtlinearen Brechungsindex n2 und des nichtlinearen Absorptionskoeffizienten 2 von Wellenleitern aus P3BT-ra im Bereich von 700 - 1500 nm gemessen. Wir haben große Werte des nichtlinearen Brechungsindex bis 1,5x10-13 cm2/W bei 1150 nm beobachtet. Wir haben gefunden, dass die Gütenkriterien (“figures of merit“) für rein optische Schalter im Wellenlängebereich 1050 - 1200 nm erfüllt sind. Dieser Bereich entspricht dem niederenergetischen Ausläufer der Zwei-Photonen-Absorption. Die Gütekriterien von P3BT-ra gehören zu den besten der bisher bekannten Werte von konjugierten Polymeren.rnWir haben gefunden, dass P3BT-ra ein vielversprechender Kandidat für integriert-optische Schalter ist, weil es eine gute Kombination aus großer Nichtlinearität dritter Ordnung, geringen Wellenleiterdämpfungverlusten und ausreichender Photostabilität zeigt. rnWir haben einen Vergleich der gemessenen Dispersion von gw, n2 und 2 mit der Theorie durchgeführt. Durch Kurvenanpassung der Dispersion von gw haben wir gefunden, dass Rayleigh-Streuung der dominierende Dämpfungsmechanismus in MEH-PPV und P3BT-ra Wellenleitern ist. Ein quantenmechanischer Ansatz wurde zur Berechnung der nichtlinearen Suszeptibilität dritter Ordnung (3) verwendet, um die gemessenen Spektren von n2 und 2 von P3BT-ra und MEH-PPV zu simulieren. Dies kann erklären, dass sättigbare Absorption und Zwei-Photonen Absorption die hauptsächlichen Effekte sind, welche die Dispersion von n2 und 2 verursachen. 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.

Einzelmolekülspektroskopische und quantenchemische Untersuchungen zum elektronischen Energietransfer

In der vorliegenden Arbeit wurden Untersuchungen zum Mechanismus, der Dynamik und der Kontrolle des elektronischen Energietransfers in multichromophoren Modellsystemen durchgeführt. Als Untersuchungsmethoden wurden hauptsächlich die konfokale Einzelmolekülspektroskopie und die Quantenchemie eingesetzt. Der Aufbau des Einzelmolekülmikroskops wurde bezüglich der Anregungs- und Detektionskomponenten variiert, um die unterschiedlichen Experimente durchzuführen. Die quantenchemischen Rechnungen wurden auf Dichtefunktional- und Coupled-Cluster-Niveau durchgeführt. Die aus den Rechnungen erhaltenen zusätzlichen Informationen über experimentell zum Teil schwer zugängliche Eigenschaften der Farbstoffe unterstützten die Interpretation der experimentellen Befunde. rnIn früheren Untersuchungen der AG Basché wurden die Energietransfer-Raten von Donor-Akzeptor-Systemen gemessen, die erhebliche Abweichungen von nach der Förster-Theorie berechneten Raten zeigten. Daher war ein Ziel der vorliegenden Arbeit, diese Abweichungen zu erklären. Zu diesem Zweck wurde die Geometrie der Diaden experimentell untersucht, sowie die elektronische Kopplung zwischen den Chromophoren quantenchemisch berechnet. Die relative Orientierung der Chromophore in den Diaden wurde in einem Einzelmolekül-Experiment mit rotierender Anregungspolarisation abgefragt. Die erhaltenen Winkelverteilungen konnten schließlich eindeutig auf die Flexibilität der die Chromophore verbrückenden Oligophenyl-Einheiten zurückgeführt werden. Die Unterschiede der gemessenen Energietransfer-Raten zu den nach der Förster-Theorie ermittelten Werten konnten jedoch nicht über die molekulare Flexibilität der Systeme erklärt werden. Aufklärung über die Diskrepanzen zur Förster-Theorie ergaben die quantenchemischen Rechnungen. In Rahmen dieser Arbeit wurde zum ersten Mal die Coupled-Cluster-Theorie zur Berechnung der elektronischen Kopplung eingesetzt. Die Betrachtung der isolierten Chromophore reichte aber nicht aus, um die gemessenen Abweichungen von der Förster-Theorie zu erklären. Erst über die Berücksichtigung der molekularen Brücke konnten die gefunden Abweichungen erklärt werden. Die deutliche Verstärkung der elektronischen Kopplung ist auf die Polarisierbarkeit der Brücke zurückzuführen.rnNach diesen Betrachtungen stand die Kontrolle des Energietransfers im Fokus der weiteren Untersuchungen. In den durchgeführten Einzelmolekülexperimenten wurden die Chromophore der Donor-Akzeptor-Systeme selektiv mit zwei Laserpulsen unterschiedlicher Wellenlänge angeregt. Beim gleichzeitigen Anregen beider Chromophore wurde Singulett-Singulett-Annihilation (SSA) induziert, ein Energietransferprozess, bei dem die Anregungsenergie vom vorigen Akzeptor zum vorigen Donor übertragen wird. Da über SSA Fluoreszenzphotonen gelöscht wurden, konnte über den Abstand der Laserpulse die Fluoreszenzintensität des einzelnen Moleküls moduliert werden. Konzeptionell verwandte Einzelmolekülexperimente wurden an einem weiteren molekularen System durchgeführt, das aus einem Kern und einer Peripherie bestand. Fluoreszenzauszeiten des Gesamtsystems bei selektiver Anregung des Kerns wurden auf die Population eines Triplett-Zustandes zurückgeführt, der die Fluoreszenz der Peripherie löschte. rnAbschließend wurde der SSA-Prozess zwischen zwei gleichartigen Chromophoren untersucht. Es wurde eine Methode entwickelt, die es zum ersten Mal erlaubte, die SSA-Zeitkonstante individueller Moleküle zu bestimmen. Hierfür wurden die Daten der gemessenen Photonen-Koinzidenzhistogramme mittels eines im Rahmen dieser Arbeit hergeleiteten analytischen Zusammenhangs ausgewertet, der über Monte-Carlo-Simulationen bestätigt wurde.

Scale bridging concepts in molecular simulation : coarse-graining and thermodynamic coupling

In dieser Arbeit wurden Simulation von Flüssigkeiten auf molekularer Ebene durchgeführt, wobei unterschiedliche Multi-Skalen Techniken verwendet wurden. Diese erlauben eine effektive Beschreibung der Flüssigkeit, die weniger Rechenzeit im Computer benötigt und somit Phänomene auf längeren Zeit- und Längenskalen beschreiben kann.rnrnEin wesentlicher Aspekt ist dabei ein vereinfachtes (“coarse-grained”) Modell, welches in einem systematischen Verfahren aus Simulationen des detaillierten Modells gewonnen wird. Dabei werden ausgewählte Eigenschaften des detaillierten Modells (z.B. Paar-Korrelationsfunktion, Druck, etc) reproduziert.rnrnEs wurden Algorithmen untersucht, die eine gleichzeitige Kopplung von detaillierten und vereinfachten Modell erlauben (“Adaptive Resolution Scheme”, AdResS). Dabei wird das detaillierte Modell in einem vordefinierten Teilvolumen der Flüssigkeit (z.B. nahe einer Oberfläche) verwendet, während der Rest mithilfe des vereinfachten Modells beschrieben wird.rnrnHierzu wurde eine Methode (“Thermodynamische Kraft”) entwickelt um die Kopplung auch dann zu ermöglichen, wenn die Modelle in verschiedenen thermodynamischen Zuständen befinden. Zudem wurde ein neuartiger Algorithmus der Kopplung beschrieben (H-AdResS) der die Kopplung mittels einer Hamilton-Funktion beschreibt. In diesem Algorithmus ist eine zur Thermodynamischen Kraft analoge Korrektur mit weniger Rechenaufwand möglich.rnrnAls Anwendung dieser grundlegenden Techniken wurden Pfadintegral Molekulardynamik (MD) Simulationen von Wasser untersucht. Mithilfe dieser Methode ist es möglich, quantenmechanische Effekte der Kerne (Delokalisation, Nullpunktsenergie) in die Simulation einzubeziehen. Hierbei wurde zuerst eine Multi-Skalen Technik (“Force-matching”) verwendet um eine effektive Wechselwirkung aus einer detaillierten Simulation auf Basis der Dichtefunktionaltheorie zu extrahieren. Die Pfadintegral MD Simulation verbessert die Beschreibung der intra-molekularen Struktur im Vergleich mit experimentellen Daten. Das Modell eignet sich auch zur gleichzeitigen Kopplung in einer Simulation, wobei ein Wassermolekül (beschrieben durch 48 Punktteilchen im Pfadintegral-MD Modell) mit einem vereinfachten Modell (ein Punktteilchen) gekoppelt wird. Auf diese Weise konnte eine Wasser-Vakuum Grenzfläche simuliert werden, wobei nur die Oberfläche im Pfadintegral Modell und der Rest im vereinfachten Modell beschrieben wird.