Spinpolarisierte Rastertunnelmikroskopie magnetischer Nanostrukturen am Beispiel von bcc-Co/Fe(110), Fe/Mo(110) und Kupfer-Phthalocyanin/Fe(110)

In dieser Arbeit werden, nach einer Einführung in die spinpolarisierte Rastertunnelmikroskopie und -spektroskopie als experimentelle Methode zur Untersuchung magnetischer Nanostrukturen, Ergebnisse zur spinpolarisierten elektronischen Struktur in Abhängigkeit von der Kristallstruktur am Beispiel ultradünner Co-Schichten sowie in Abhängigkeit von der Magnetisierungsrichtung für ultradünne Fe-Schichten vorgestellt. Hochaufgelöste Messungen zeigen die ortsabhängige Spinpolarisation auf einem einzelnen Kupfer-Phthalocyanin Molekül. rnrnKobalt wurde durch pseudomorphes Wachstum auf den (110)-Oberflächen der kubisch raumzentrierten Metalle Chrom und Eisen deponiert. Im Unterschied zu früheren Berichten in der Literatur lassen sich nur zwei Lagen Co in der kubisch raumzentrierten (bcc) Ordnung stabilisieren. Die bcc-Co Schichten auf der Fe(110)-Oberfläche zeigen keine Anzeichen von epitaktischen Verzerrungen. rnDickere Schichten rekonstruieren in eine dicht gepackte Struktur (hcp/fcc). Durch die bcc Ordnung wird die Spinpolarisation von Kobalt auf P=62% erhöht (hcp-Co: P=45%). rnrnDie temperaturabhängige Spinreorientierung (SRT) ultradünner Filme Fe/Mo(110) wurde mit spinpolarisierter Spektroskopie untersucht. Eine Neuausrichtung der Magnetisierung aus der senkrechten [110]-Achse in die in der Ebene liegenden [001]-Achse wird bei T=(13,2+-0,5)K festgestellt, wobei es sich um einen diskontinuierlichen Reorientierungsübergang handelt, d.h. die freie Energie weist innerhalb eines bestimmten Temperaturbereichs gleichzeitig zwei Minima auf. Weiterhin wird in der Mono- und Doppellage Fe/Mo(110 eine Abhängigkeit der elektronischen Struktur von der Ausrichtung der magnetisch leichten Achse und von der Magnetisierung beobachtet. rnrnDie Untersuchung des spinpolarisierten Ladungstransports durch ein Kupfer-Phthalocyanin-Molekül auf der Fe/Mo(110) Oberfläche liefert einen wesentlichen Beitrag zum Verständnis des Spintransports an der Grenzfläche zwischen Metall und organischem Molekül. Die HOMO-LUMO-Energielücke des freien Moleküls wird durch die Wechselwirkung mit der Metalloberfläche mit Grenzflächenzuständen gefüllt. Diese Zustände reduzieren die Spinpolarisation des durch das Molekül fließenden Tunnelstroms durch einen zusätzlichen unpolarisierten Strombeitrag um einen Faktor zwei. Spinpolarisierte hybridisierte Grenzflächenzustände mit größerem Abstand zur Fermi-Energie führen in Abhängigkeit von der Position auf dem Molekül zu weiteren Beiträgen zur effektiven Spinpolarisation. Diese Untersuchungen belegen die Möglichkeit einer effektiven Spininjektion in organische Halbleiter und damit das Potential dieser Materialien für die weitere Entwicklung von Spintronik-Bauteilen.

Computer simulations of two-dimensional colloidal crystals under confinement and shear

In this thesis we are presenting a broadly based computer simulation study of two-dimensional colloidal crystals under different external conditions. In order to fully understand the phenomena which occur when the system is being compressed or when the walls are being sheared, it proved necessary to study also the basic motion of the particles and the diffusion processes which occur in the case without these external forces. In the first part of this thesis we investigate the structural transition in the number of rows which occurs when the crystal is being compressed by placing the structured walls closer together. Previous attempts to locate this transition were impeded by huge hysteresis effects. We were able to determine the transition point with higher precision by applying both the Schmid-Schilling thermodynamic integration method and the phase switch Monte Carlo method in order to determine the free energies. These simulations showed not only that the phase switch method can successfully be applied to systems with a few thousand particles and a soft crystalline structure with a superimposed pattern of defects, but also that this method is way more efficient than a thermodynamic integration when free energy differences are to be calculated. Additionally, the phase switch method enabled us to distinguish between several energetically very similar structures and to determine which one of them was actually stable. Another aspect considered in the first result chapter of this thesis is the ensemble inequivalence which can be observed when the structural transition is studied in the NpT and in the NVT ensemble. The second part of this work deals with the basic motion occurring in colloidal crystals confined by structured walls. Several cases are compared where the walls are placed in different positions, thereby introducing an incommensurability into the crystalline structure. Also the movement of the solitons, which are created in the course of the structural transition, is investigated. Furthermore, we will present results showing that not only the well-known mechanism of vacancies and interstitial particles leads to diffusion in our model system, but that also cooperative ring rotation phenomena occur. In this part and the following we applied Langevin dynamics simulations. In the last chapter of this work we will present results on the effect of shear on the colloidal crystal. The shear was implemented by moving the walls with constant velocity. We have observed shear banding and, depending on the shear velocity, that the inner part of the crystal breaks into several domains with different orientations. At very high shear velocities holes are created in the structure, which originate close to the walls, but also diffuse into the inner part of the crystal.

Nucleation studies on graphics processing units

Ein System in einem metastabilen Zustand muss eine bestimmte Barriere in derrnfreien Energie überwinden um einen Tropfen der stabilen Phase zu formen.rnHerkömmliche Untersuchungen nehmen hierbei kugelförmige Tropfen an. Inrnanisotropen Systemen (wie z.B. Kristallen) ist diese Annahme aber nicht ange-rnbracht. Bei tiefen Temperaturen wirkt sich die Anisotropie des Systems starkrnauf die freie Energie ihrer Oberfläche aus. Diese Wirkung wird oberhalb derrnAufrauungstemperatur T R schwächer. Das Ising-Modell ist ein einfaches Mo-rndell, welches eine solche Anisotropie aufweist. Wir führen großangelegte Sim-rnulationen durch, um die Effekte, die mit einer endlichen Simulationsbox ein-rnhergehen, sowie statistische Ungenauigkeiten möglichst klein zu halten. DasrnAusmaß der Simulationen die benötigt werden um sinnvolle Ergebnisse zu pro-rnduzieren, erfordert die Entwicklung eines skalierbaren Simulationsprogrammsrnfür das Ising-Modell, welcher auf verschiedenen parallelen Architekturen (z.B.rnGrafikkarten) verwendet werden kann. Plattformunabhängigkeit wird durch ab-rnstrakte Schnittstellen erreicht, welche plattformspezifische Implementierungs-rndetails verstecken. Wir benutzen eine Systemgeometrie die es erlaubt eine Ober-rnfläche mit einem variablen Winkel zur Kristallebene zu untersuchen. Die Ober-rnfläche ist in Kontakt mit einer harten Wand, wobei der Kontaktwinkel Θ durchrnein Oberflächenfeld eingestellt werden kann. Wir leiten eine Differenzialglei-rnchung ab, welche das Verhalten der freien Energie der Oberfläche in einemrnanisotropen System beschreibt. Kombiniert mit thermodynamischer Integrationrnkann die Gleichung benutzt werden, um die anisotrope Oberflächenspannungrnüber einen großen Winkelbereich zu integrieren. Vergleiche mit früheren Mes-rnsungen in anderen Geometrien und anderen Methoden zeigen hohe Überein-rnstimung und Genauigkeit, welche vor allem durch die im Vergleich zu früherenrnMessungen wesentlich größeren Simulationsdomänen erreicht wird. Die Temper-rnaturabhängigkeit der Oberflächensteifheit κ wird oberhalb von T R durch diernKrümmung der freien Energie der Oberfläche für kleine Winkel gemessen. DiesernMessung lässt sich mit Simulationsergebnissen in der Literatur vergleichen undrnhat bessere Übereinstimmung mit theoretischen Voraussagen über das Skalen-rnverhalten von κ. Darüber hinaus entwickeln wir ein Tieftemperatur-Modell fürrndas Verhalten um Θ = 90 Grad weit unterhalb von T R. Der Winkel bleibt bis zu einemrnkritischen Feld H C quasi null; oberhalb des kritischen Feldes steigt der Winkelrnrapide an. H C wird mit der freien Energie einer Stufe in Verbindung gebracht,rnwas es ermöglicht, das kritische Verhalten dieser Größe zu analysieren. Die harternWand muss in die Analyse einbezogen werden. Durch den Vergleich freier En-rnergien bei geschickt gewählten Systemgrößen ist es möglich, den Beitrag derrnKontaktlinie zur freien Energie in Abhängigkeit von Θ zu messen. Diese Anal-rnyse wird bei verschiedenen Temperaturen durchgeführt. Im letzten Kapitel wirdrneine 2D Fluiddynamik Simulation für Grafikkarten parallelisiert, welche u. a.rnbenutzt werden kann um die Dynamik der Atmosphäre zu simulieren. Wir im-rnplementieren einen parallelen Evolution Galerkin Operator und erreichen

Numerical simulation of some viscoelastic fluids

Numerical simulation of the Oldroyd-B type viscoelastic fluids is a very challenging problem. rnThe well-known High Weissenberg Number Problem" has haunted the mathematicians, computer scientists, and rnengineers for more than 40 years. rnWhen the Weissenberg number, which represents the ratio of elasticity to viscosity, rnexceeds some limits, simulations done by standard methods break down exponentially fast in time. rnHowever, some approaches, such as the logarithm transformation technique can significantly improve rnthe limits of the Weissenberg number until which the simulations stay stable. rnrnWe should point out that the global existence of weak solutions for the Oldroyd-B model is still open. rnLet us note that in the evolution equation of the elastic stress tensor the terms describing diffusive rneffects are typically neglected in the modelling due to their smallness. However, when keeping rnthese diffusive terms in the constitutive law the global existence of weak solutions in two-space dimension rncan been shown. rnrnThis main part of the thesis is devoted to the stability study of the Oldroyd-B viscoelastic model. rnFirstly, we show that the free energy of the diffusive Oldroyd-B model as well as its rnlogarithm transformation are dissipative in time. rnFurther, we have developed free energy dissipative schemes based on the characteristic finite element and finite difference framework. rnIn addition, the global linear stability analysis of the diffusive Oldroyd-B model has also be discussed. rnThe next part of the thesis deals with the error estimates of the combined finite element rnand finite volume discretization of a special Oldroyd-B model which covers the limiting rncase of Weissenberg number going to infinity. Theoretical results are confirmed by a series of numerical rnexperiments, which are presented in the thesis, too.

Monte Carlo simulations of nucleation of colloidal crystals

A crystal nucleus in a finite volume may exhibit phase coexistence with a surrounding fluid. The thermodynamic properties of the coexisting fluid (pressure and chemical potential) are enhanced relative to their coexistence values. This enhancement is uniquely related to the surface excess free energy. rnA model for weakly attractive soft colloidal particles is investigated, the so called Asakura-Oosawa model. In simulations, this model allows for the calculation of the pressure in the liquid using the virial formula directly. The phase coexistence pressure in the thermodynamic limit is obtained from the interface velocity method. We introduce a method by which the chemical potential in dense liquids can be measured. There is neither a need to locate the interface nor to compute the anisotropic interfacial tension to obtain nucleation barriers. Therefore, our analysis is appropriate for nuclei of arbitrary shape. Monte Carlo simulations over a wide range of nucleus volumes yield to nucleation barriers independent from the total system volume. The interfacial tension is determined via the ensemble-switch method, hence a detailed test of classical nucleation theory is possible. The anisotropy of the interfacial tension and the resulting non-spherical shape has only a minor effect on the barrier for the Asakura-Oosawa model.

Propionsäureabbau in NawaRo-Biogasanlagen

Die Energiewende ist begleitet von dem Ausbau erneuerbarer Energien. Dabei spielt die Energiegewinnung aus Biomasse eine wichtige Rolle. Der optimale Betrieb einer Biogasanlage erfordert eine stabile Methanproduktion, welche jedoch durch die Akkumulation von Propionsäure nachhaltig gestört werden kann. Aus diesem Grund ist der mikrobielle Abbau dieser Substanz von besonderem Interesse. Die Thermodynamik des anaeroben bakteriellen Abbaus von Propionsäure erfordert die syntrophe Verwertung des entstehenden Wasserstoffs durch Wasserstoff-verbrauchende Mikroorganismen, beispielsweise methanogene Archaea.rnMit dem Ziel, die Erkenntnislage der Propionat-Verwertung in NawaRo-Biogasanlagen zu erweitern, sollten Propionat-verwertende Anreicherungskulturen aus NawaRo-Biogasanlagen etabliert, charakterisiert und molekularbiologisch analysiert werden.rnAus landwirtschaftlichen Biogasanlagen wurden reproduzierbar Propionat-verwertende Anreicherungskulturen mittels anaerober Kultivierungstechniken etabliert. Die anaerob Propionat-verwertende Aktivität der Kulturen blieb über Jahre erhalten und konnte unter verschiedenen Bedingungen charakterisiert werden. Die Analyse der sukzessiven Diversität von vier Anreicherungskulturen ermöglichte einen Einblick in die sich während der Propionat-Verwertung sukzessiv verändernde mikrobielle Diversität. Dabei wurden die aus der 16S rDNA-Analyse resultierenden Sequenzcluster MP-1 (Cryptanaerobacter sp./ Pelotomaculum sp.), MP-6 und MP-15 (beide ''Candidatus Cloacamonas sp. ''), sowie MP-9 (Syntrophobacter sulfatireducens) als potentiell Propionat-verwertende Schlüsselspezies identifiziert. Mit S. sulfatireducens wurde eine bekannte syntroph Propionat-verwertende Spezies gefunden. Die Sequenzen von MP-1 waren nahe verwandt mit Pelotomaculum schinkii, ebenfalls eine beschriebene syntroph Propionat-verwertende Spezies. Bei dem nächsten Verwandten der Cluster MP-6 und MP-15 handelte es sich um ''Candidatus Cloacamonas acidaminovorans'', eine bisher unkultivierbare Spezies, dessen Genom für den gesamten Abbauweg der syntrophen Propionat-Oxidation codiert. Syntrophobacter sulfatireducens kam zusammen mit Vertretern der methanogenen Gattungen Methanoculleus, Methanosaeta und Methanomethylovorans vor. Als methanogener Partner von Cryptanaerobacter sp./ Pelotomaculum sp. dominierte die Gattung Methanosarcina. Aufgrund der starken Präsenz der syntroph Acetat oxidierenden Spezies Tepidanaerobacter acetatoxydans (Sequenz-Cluster MP-3), sowie potentiell homoacetogener Arten, wurde zudem ein theoretischer Zusammenhang der Propionat-Verwertung mit der syntrophen Acetat-Oxidation und der autotrophen Homoacetogenese vorgeschlagen.rn

Computer simulation of rod-coil block copolymers and tetrapod,polymer mixtures

In this thesis I present a new coarse-grained model suitable to investigate the phase behavior of rod-coil block copolymers on mesoscopic length scales. In this model the rods are represented by hard spherocylinders, whereas the coil block consists of interconnected beads. The interactions between the constituents are based on local densities. This facilitates an efficient Monte-Carlo sampling of the phase space. I verify the applicability of the model and the simulation approach by means of several examples. I treat pure rod systems and mixtures of rod and coil polymers. Then I append coils to the rods and investigate the role of the different model parameters. Furthermore, I compare different implementations of the model. I prove the capability of the rod-coil block copolymers in our model to exhibit typical micro-phase separated configurations as well as extraordinary phases, such as the wavy lamellar state, percolating structuresrnand clusters. Additionally, I demonstrate the metastability of the observed zigzag phase in our model. A central point of this thesis is the examination of the phase behavior of the rod-coil block copolymers in dependence of different chain lengths and interaction strengths between rods and coil. The observations of these studies are summarized in a phase diagram for rod-coil block copolymers. Furthermore, I validate a stabilization of the smectic phase with increasing coil fraction.rnIn the second part of this work I present a side project in which I derive a model permitting the simulation of tetrapods with and without grafted semiconducting block copolymers. The effect of these polymers is added in an implicit manner by effective interactions between the tetrapods. While the depletion interaction is described in an approximate manner within the Asakura-Oosawa model, the free energy penalty for the brush compression is calculated within the Alexander-de Gennes model. Recent experiments with CdSe tetrapods show that grafted tetrapods are clearly much better dispersed in the polymer matrix than bare tetrapods. My simulations confirm that bare tetrapods tend to aggregate in the matrix of excess polymers, while clustering is significantly reduced after grafting polymer chains to the tetrapods. Finally, I propose a possible extension enabling the simulation of a system with fluctuating volume and demonstrate its basic functionality. This study is originated in a cooperation with an experimental group with the goal to analyze the morphology of these systems in order to find the ideal morphology for hybrid solar cells.

Theoretical Calculations of Acid Dissociation Constants: A Review Article

Acid dissociation constants, or pKa values, are essential for understanding many fundamental reactions in chemistry. These values reveal the deprotonation state of a molecule in a particular solvent. There is great interest in using theoretical methods to calculate the pKa values for many different types of molecules. These include molecules that have not been synthesized, those for which experimental pKa determinations are difficult, and for larger molecules where the local environment changes the usual pKa values, such as for certain amino acids that are part of a larger polypeptide chain. Chemical accuracy in pKa calculations is difficult to achieve, because an error of 1.36 kcal/mol in the change of free energy for deprotonation in solvent results in an error of 1 pKa unit. In this review the most valuable methods for determining accurate pKa values in aqueous solution are presented for educators interested in explaining or using these methods for their students.

Absolute pKa Determinations for Substituted Phenols

The CBS-QB3 method was used to calculate the gas-phase free energy difference between 20 phenols and their respective anions, and the CPCM continuum solvation method was applied to calculate the free energy differences of solvation for the phenols and their anions. The CPCM solvation calculations were performed on both gas-phase and solvent-phase optimized structures. Absolute pKa calculations with solvated phase optimized structures for the CPCM calculations yielded standard deviations and root-mean-square errors of less than 0.4 pKa unit. This study is the most accurate absolute determination of the pKa values of phenols, and is among the most accurate of any such calculations for any group of compounds. The ability to make accurate predictions of pKa values using a coherent, well-defined approach, without external approximations or fitting to experimental data, is of general importance to the chemical community. The solvated phase optimized structures of the anions are absolutely critical to obtain this level of accuracy, and yield a more realistic charge separation between the negatively charged oxygen and the ring system of the phenoxide anions.

Experimentation with different thermodynamic cycles used for pKa calculations on carboxylic acids using Complete Basis Set and Gaussian-n Models combined with CPCM Continuum Solvation Methods

Complete basis set and Gaussian-n methods were combined with Barone and Cossi's implementation of the polarizable conductor model (CPCM) continuum solvation methods to calculate pKa values for six carboxylic acids. Four different thermodynamic cycles were considered in this work. An experimental value of −264.61 kcal/mol for the free energy of solvation of H+, ΔGs(H+), was combined with a value for Ggas(H+) of −6.28 kcal/mol, to calculate pKa values with cycle 1. The complete basis set gas-phase methods used to calculate gas-phase free energies are very accurate, with mean unsigned errors of 0.3 kcal/mol and standard deviations of 0.4 kcal/mol. The CPCM solvation calculations used to calculate condensed-phase free energies are slightly less accurate than the gas-phase models, and the best method has a mean unsigned error and standard deviation of 0.4 and 0.5 kcal/mol, respectively. Thermodynamic cycles that include an explicit water in the cycle are not accurate when the free energy of solvation of a water molecule is used, but appear to become accurate when the experimental free energy of vaporization of water is used. This apparent improvement is an artifact of the standard state used in the calculation. Geometry relaxation in solution does not improve the results when using these later cycles. The use of cycle 1 and the complete basis set models combined with the CPCM solvation methods yielded pKa values accurate to less than half a pKa unit. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001

Accurate pKa Calculations for Carboxylic Acids Using Complete Basis Set and Gaussian-n Models Combined with CPCM Continuum Solvation Methods”

Complete Basis Set and Gaussian-n methods were combined with CPCM continuum solvation methods to calculate pKa values for six carboxylic acids. An experimental value of −264.61 kcal/mol for the free energy of solvation of H+, ΔGs(H+), was combined with a value for Ggas(H+) of −6.28 kcal/mol to calculate pKa values with Cycle 1. The Complete Basis Set gas-phase methods used to calculate gas-phase free energies are very accurate, with mean unsigned errors of 0.3 kcal/mol and standard deviations of 0.4 kcal/mol. The CPCM solvation calculations used to calculate condensed-phase free energies are slightly less accurate than the gas-phase models, and the best method has a mean unsigned error and standard deviation of 0.4 and 0.5 kcal/mol, respectively. The use of Cycle 1 and the Complete Basis Set models combined with the CPCM solvation methods yielded pKa values accurate to less than half a pKa unit.

Further Quantum Mechanical Evidence that Difluorotoluene does not Hydrogen Bond

Calculations were run on the methylated DNA base pairs adenine:thymine and adenine:difluorotoluene to further investigate the hydrogen-bonding properties of difluorotoluene (F). Geometries were optimized using hybrid density functional theory. Single-point calculations at the MP2(full) level were performed to obtain more rigorous energies. The functional counterpoise method was used to correct for the basis set superposition error (BSSE), and the interaction energies were also corrected for fragment relaxation. These corrections brought the B3LYP and MP2 interaction energies into excellent agreement. In the gas phase, the Gibbs free energies calculated at the B3LYP and MP2 levels of theory predict that A and T will spontaneously form an A:T pair while A:F spontaneously dissociates into A and F. Solvation effects on the pairing of the bases were explored using implicit solvent models for water and chloroform. In aqueous solution, both A:T and A:F are predicted to dissociate into their component monomers. Semiempirical calculations were performed on small sections of B-form DNA containing the two pairs, and the results provide support for the concept that base stacking is more important than hydrogen bonding for the stability of the A:F pair within a DNA helix.

Accurate relative pKa calculations for carboxylic acids using complete basis set and Gaussian-n models combined with continuum solvation methods

The complete basis set methods CBS-4, CBS-QB3, and CBS-APNO, and the Gaussian methods G2 and G3 were used to calculate the gas phase energy differences between six different carboxylic acids and their respective anions. Two different continuum methods, SM5.42R and CPCM, were used to calculate the free energy differences of solvation for the acids and their anions. Relative pKa values were calculated for each acid using one of the acids as a reference point. The CBS-QB3 and CBS-APNO gas phase calculations, combined with the CPCM/HF/6-31+G(d)//HF/6-31G(d) or CPCM/HF/6-31+G(d)//HF/6-31+G(d) continuum solvation calculations on the lowest energy gas phase conformer, and with the conformationally averaged values, give results accurate to ½ pKa unit. © 2001 American Institute of Physics.

First-principle studies of intermolecular and intramolecular catalysis of protonated cocaine

We have performed a series of first-principles electronic structure calculations to examine the reaction pathways and the corresponding free energy barriers for the ester hydrolysis of protonated cocaine in its chair and boat conformations. The calculated free energy barriers for the benzoyl ester hydrolysis of protonated chair cocaine are close to the corresponding barriers calculated for the benzoyl ester hydrolysis of neutral cocaine. However, the free energy barrier calculated for the methyl ester hydrolysis of protonated cocaine in its chair conformation is significantly lower than for the methyl ester hydrolysis of neutral cocaine and for the dominant pathway of the benzoyl ester hydrolysis of protonated cocaine. The significant decrease of the free energy barrier, ∼4 kcal/mol, is attributed to the intramolecular acid catalysis of the methyl ester hydrolysis of protonated cocaine, because the transition state structure is stabilized by the strong hydrogen bond between the carbonyl oxygen of the methyl ester moiety and the protonated tropane N. The relative magnitudes of the free energy barriers calculated for different pathways of the ester hydrolysis of protonated chair cocaine are consistent with the experimental kinetic data for cocaine hydrolysis under physiologic conditions. Similar intramolecular acid catalysis also occurs for the benzoyl ester hydrolysis of (protonated) boat cocaine in the physiologic condition, although the contribution of the intramolecular hydrogen bonding to transition state stabilization is negligible. Nonetheless, the predictability of the intramolecular hydrogen bonding could be useful in generating antibody-based catalysts that recruit cocaine to the boat conformation and an analog that elicited antibodies to approximate the protonated tropane N and the benzoyl O more closely than the natural boat conformer might increase the contribution from hydrogen bonding. Such a stable analog of the transition state for intramolecular catalysis of cocaine benzoyl-ester hydrolysis was synthesized and used to successfully elicit a number of anticocaine catalytic antibodies.

Pople's Gaussian-3 model chemistry applied to an investigation of (H2O)8 water clusters

The Gaussian-3 method developed by Pople and coworkers has been used to calculate the free energy of neutral octamer clusters of water, (H2O)8. The most energetically stable structures are in excellent agreement with those determined from experiment and those predicted from previous high-level calculations. Cubic structures are favored over noncubic structures over all temperature ranges studied. The D2d cubic structure is the lowest free energy structure and dominates the potential energy and free energy hypersurfaces from 0 K to 298 K.