Thermodynamic concepts in adaptive resolution simulations

This thesis work is devoted to the conceptual and technical development of the Adaptive Resolution Scheme (AdResS), a molecular dynamics method that allows the simulation of a system with different levels of resolution simultaneously. The simulation domain is divided into high and low resolution zones and a transition region that links them, through which molecules can freely diffuse.rnThe first issue of this work regards the thermodynamic consistency of the method, which is tested and verified in a model liquid of tetrahedral molecules. The results allow the introduction of the concept of the Thermodynamic Force, an external field able to correct spurious density fluctuations present in the transition region in usual AdResS simulations.rnThe AdResS is also applied to a system where two different representations with the same degree of resolution are confronted. This simple test extends the method from an Adaptive Resolution Scheme to an Adaptive Representation Scheme, providing a way of coupling different force fields based on thermodynamic consistency arguments. The Thermodynamic Force is successfully applied to the example described in this work as well.rnAn alternative approach of deducing the Thermodynamic Force from pressure consistency considerations allows the interpretation of AdResS as a first step towards a molecular dynamics simulation in the Grand Canonical ensemble. Additionally, such a definition leads to a practical way of determining the Thermodynamic Force, tested in the well studied tetrahedral liquid. The effects of AdResS and this correction on the atomistic domain are analyzed by inspecting the local distribution of velocities, radial distribution functions, pressure and particle number fluctuation. Their comparison with analogous results coming from purely atomistic simulations shows good agreement, which is greatly improved under the effect of the external field.rnA further step in the development of AdResS, necessary for several applications in biophysics and material science, consists of its application to multicomponent systems. To this aim, the high-resolution representation of a model binary mixture is confronted with its coarse-grained representation systematically parametrized. The Thermodynamic Force, whose development requires a more delicate treatment, also gives satisfactory results.rnFinally, AdResS is tested in systems including two-body bonded forces, through the simulation of a model polymer allowed to adaptively change its representation. It is shown that the distribution functions that characterize the polymer structure are in practice not affected by the change of resolution.rnThe technical details of the implementation of AdResS in the ESPResSo package conclude this thesis work.

HCl gas gettering of 3d transition metals for crystalline silicon solar cell concepts

In order to reduce the costs of crystalline silicon solar cells, low-cost silicon materials like upgraded metallurgical grade (UMG) silicon are investigated for the application in the photovoltaic (PV) industry. Conventional high-purity silicon is made by cost-intensive methods, based on the so-called Siemens process, which uses the reaction to form chlorosilanes and subsequent several distillation steps before the deposition of high-purity silicon on slim high-purity silicon rods. UMG silicon in contrast is gained from metallurgical silicon by a rather inexpensive physicochemical purification (e.g., acid leaching and/or segregation). However, this type of silicon usually contains much higher concentrations of impurities, especially 3d transition metals like Ti, Fe, and Cu. These metals are extremely detrimental in the electrically active part of silicon solar cells, as they form recombination centers for charge carriers in the silicon band gap. This is why simple purification techniques like gettering, which can be applied between or during solar cell process steps, will play an important role for such low-cost silicon materials. Gettering in general describes a process, whereby impurities are moved to a place or turned into a state, where they are less detrimental to the solar cell. Hydrogen chloride (HCl) gas gettering in particular is a promising simple and cheap gettering technique, which is based on the reaction of HCl gas with transition metals to form volatile metal chloride species at high temperatures.rnThe aim of this thesis was to find the optimum process parameters for HCl gas gettering of 3d transition metals in low-cost silicon to improve the cell efficiency of solar cells for two different cell concepts, the standard wafer cell concept and the epitaxial wafer equivalent (EpiWE) cell concept. Whereas the former is based on a wafer which is the electrically active part of the solar cell, the latter uses an electrically inactive low-cost silicon substrate with an active layer of epitaxially grown silicon on top. Low-cost silicon materials with different impurity grades were used for HCl gas gettering experiments with the variation of process parameters like the temperature, the gettering time, and the HCl gas concentration. Subsequently, the multicrystalline silicon neighboring wafers with and without gettering were compared by element analysis techniques like neutron activation analysis (NAA). It was demonstrated that HCl gas gettering is an effective purification technique for silicon wafers, which is able to reduce some 3d transition metal concentrations by over 90%. Solar cells were processed for both concepts which could demonstrate a significant increase of the solar cell efficiency by HCl gas gettering. The efficiency of EpiWE cells could be increased by HCl gas gettering by approximately 25% relative to cells without gettering. First process simulations were performed based on a simple model for HCl gas gettering processes, which could be used to make qualitative predictions.

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.

Imiquimod based transcutaneous immunization – insights and novel concepts

This work aims at developing a transcutaneous immunization (TCI) approach in order to activate cytotoxic T-cells. A tumor specific immune response was therefore generated by the TLR7-Agonist imiquimod. Five commercially available creams including the innovators product Aldara® 5% creme were assessed to ascertain their capability to induce an immune response in C57BL/6 mice after dermal administration. Moreover, creams were investigated regarding their imiquimod permeation in a Franz-diffusion cell model. Results obtained from this study were used to develop novel formulation approaches based on dissolved state imiquimod in a submicron scale range. High pressure homogenization ensured emulsification as well as particle size reduction. A freeze dried spreadable solid nanoemulsion based on sucrose fatty acid esters and oil components represented a major formulation approach. Within the scope of this approach the influence of pharmaceutical oils i.e. middle chain triglycerides, avocado oil, jojoba wax, and squalen was assessed towards their TCI performance. Furthermore, an aqueous jojoba wax based emulsion gel was developed. Unlike the innovators product, all formulations demonstrated a distinctly reduced imiquimod permeation across murine skin, a fact particularly evident in case of jojoba wax. Squalen significantly augmented in vivo immune response (p≤0.05 Mann-Whitney-Test). The emulsion gel demonstrated a 10fold decrease of imiquimod permeation. In comparison with the innovators product, the emulsion gel induced an equal immune response with a simultaneously enhanced tumor rejection in a mouse model.