Large-scale coupled-cluster calculations

Coupled-cluster theory provides one of the most successful concepts in electronic-structure theory. This work covers the parallelization of coupled-cluster energies, gradients, and second derivatives and its application to selected large-scale chemical problems, beside the more practical aspects such as the publication and support of the quantum-chemistry package ACES II MAB and the design and development of a computational environment optimized for coupled-cluster calculations. The main objective of this thesis was to extend the range of applicability of coupled-cluster models to larger molecular systems and their properties and therefore to bring large-scale coupled-cluster calculations into day-to-day routine of computational chemistry. A straightforward strategy for the parallelization of CCSD and CCSD(T) energies, gradients, and second derivatives has been outlined and implemented for closed-shell and open-shell references. Starting from the highly efficient serial implementation of the ACES II MAB computer code an adaptation for affordable workstation clusters has been obtained by parallelizing the most time-consuming steps of the algorithms. Benchmark calculations for systems with up to 1300 basis functions and the presented applications show that the resulting algorithm for energies, gradients and second derivatives at the CCSD and CCSD(T) level of theory exhibits good scaling with the number of processors and substantially extends the range of applicability. Within the framework of the ’High accuracy Extrapolated Ab initio Thermochemistry’ (HEAT) protocols effects of increased basis-set size and higher excitations in the coupled- cluster expansion were investigated. The HEAT scheme was generalized for molecules containing second-row atoms in the case of vinyl chloride. This allowed the different experimental reported values to be discriminated. In the case of the benzene molecule it was shown that even for molecules of this size chemical accuracy can be achieved. Near-quantitative agreement with experiment (about 2 ppm deviation) for the prediction of fluorine-19 nuclear magnetic shielding constants can be achieved by employing the CCSD(T) model together with large basis sets at accurate equilibrium geometries if vibrational averaging and temperature corrections via second-order vibrational perturbation theory are considered. Applying a very similar level of theory for the calculation of the carbon-13 NMR chemical shifts of benzene resulted in quantitative agreement with experimental gas-phase data. The NMR chemical shift study for the bridgehead 1-adamantyl cation at the CCSD(T) level resolved earlier discrepancies of lower-level theoretical treatment. The equilibrium structure of diacetylene has been determined based on the combination of experimental rotational constants of thirteen isotopic species and zero-point vibrational corrections calculated at various quantum-chemical levels. These empirical equilibrium structures agree to within 0.1 pm irrespective of the theoretical level employed. High-level quantum-chemical calculations on the hyperfine structure parameters of the cyanopolyynes were found to be in excellent agreement with experiment. Finally, the theoretically most accurate determination of the molecular equilibrium structure of ferrocene to date is presented.

Influence of pollination and seed dispersal on the genetic population structure of two Commiphora species in Madagascar and South Africa

Pollination and seed dispersal are important ecological processes for the regeneration of plant populations and both vectors for gene exchange between plant populations. For my thesis, I studied the pollination ecology of the South African tree Commiphora harveyi (Burseraceae) and compared it with C. guillauminii from Madagascar. Both species have low visitation rates and a low number of pollinating insect species, resulting in a low fruit set. While their pollination ecology is very similar, they differ in their seed dispersal with a low seed dispersal rate in the Malagasy and a high seed dispersal rate in the South African species. This should be reflected in a stronger genetic differentiation among populations in the Malagasy than in the South African species. My results, based on AFLP markers, contradict these expectations, the overall differentiation was lower in the Malagasy (FST = 0.05) than in the South African species (FST = 0.16). However, at a smaller spatial scale (below 3 km), the Malagasy species was genetically more strongly differentiated than the South African species, which was reflected by the high inter-population variance within the sample site (C. guillauminii: 72.2 - 85.5 %; C. harveyi: 8.4 - 14.5 %). This strong differentiation could arise from limited gene flow, which was confirmed by spatial autocorrelation analyses. The shape of the autocorrelogram suggested that gene exchange between individuals occurred only up to 3 km in the Malagasy species, whereas up to 30 km in the South African species. These results on the genetic structure correspond to the expectations based on seed dispersal data. Thus, seed dispersal seems to be a key factor for the genetic structure in plant populations on a local scale.

Stress & structure of thin polymer brush films

Polymer brushes have unique properties with a large variety of possible applications ranging from responsive coatings and drug delivery to lubrication and sensing. For further development a detailed understanding of the properties is needed. Established characterization methods, however, only supply information of the surface. Experimental data about the inner “bulk” structure of polymer brushes is still missing.rnScattering methods under grazing incidence supply structural information of surfaces as well as structures beneath it. Nanomechanical cantilevers supply stress data, which is giving information about the forces acting inside the polymer brush film. In this thesis these two techniques are further developed and used to deepen the understanding of polymer brushes. rnThe experimental work is divided into four chapters. Chapter 2 deals with the preparation of polymer brushes on top of nanomechanical cantilever sensors as well as large area sample by using a “grafting-to” technique. The further development of nanomechanical cantilever readout is subject of chapter 3. In order to simplify cantilever sensing, a method is investigated which allows one to perform multiple bending experiments on top of a single cantilever. To do so, a way to correlate different curvatures is introduced as well as a way to conveniently locate differently coated segments. In chapter 4 the change in structure upon solvent treatment of mixed polymer brushes is investigated by using scattering methods and nanomechanical cantilevers amongst others. This allows one to explain the domain memory effect, which is typically found in such systems. Chapter 5 describes the implementation of a phase shifting interferometer - used for readout of nanomechanical cantilevers - into the µ-focused scattering beamline BW4, allowing simultaneous measurements of stress and structure information. The last experimental chapter 6 deals with the roughness correlation in polymer brushes and its dependence on the chain tethered density.rnIn summary, the thesis deals with utilization of new experimental techniques for the investigation of polymer brushes and further development of the techniques themselves.rn

Mechanical properties and structure of gel systems

Understanding the origins of the mechanical properties and its correlation withrnthe microstructure of gel systems is of great scientific and industrial interest. Inrngeneral, colloidal gels can be classified into chemical and physical gels, accordingrnto the life time of the network bonds. The characteristic di↵erences in gelationrndynamics can be observed with rheological measurements.rnAs a model system, a mixture of sodium silicate and low concentration sulfuric acidrnwas used. Nano-sized silica particles grow and aggregate to a system-spanning gelrnnetwork. The influence of the finite solubility of silica at high pH on the gelationrnwas studied with classical and piezo rheometer. The storage modulus of therngel grew logarithmically with time with two distinct growth laws. A relaxationrnat low frequency was observed in the frequency dependent measurements. I attributernthese two behaviors as a sign of structural rearrangements due to the finiternsolubility of silica at high pH. The reaction equilibrium between formation andrndissolution of bonds leads to a finite life time of the bonds and behavior similar tornphysical gel. The frequency dependence was more pronounced for lower water concentrations,rnhigher temperatures and shorter reaction times. With two relaxationrnmodels, I deduced characteristic relaxation times from the experimental data. Besidesrnrheology, the evolution of silica gels at high pH on di↵erent length scales wasrnstudied by NMR and dynamic light scattering. The results revealed that the primaryrnparticles existed already in sodium silicate and aggregated after the mixingrnof reactants due to a chemical reaction. Throughout the aggregation process thernsystem was in its chemical reaction equilibrium. Applying large oscillatory shearrnstrain to the gel allowed for modifying the gel modulus. The e↵ect of shear andrnshear history on the rheological properties of the gel were investigated. The storagernmodulus of the final gel increased with increasing strain. This behavior can be explained with (i) shear-induced aggregate compaction and (ii) combination ofrnbreakage and new formation of bonds.rnIn comparison with the physical gel-like behavior of the silica gel at high pH, typicalrnchemical gel features were exhibited by other gels formed from various chemicalrnreactions. Influences of the chemical structure modification on the gelation wererninvestigated with the piezo-rheometer. The external stimuli can be applied to tunernthe mechanical properties of the gel systems.

Mesozoic fold structure of the Otago and Alpine Schist Belt and its implications on the tectonic evolution of South Island, New Zealand

Erneute Untersuchungen der mesozoischen Faltenstruktur des Otago Schiefergürtels, Südinsel, Neuseeland, zeigen, dass diese aus zwei aufeinander folgenden, ähnlichen, asymmetrischen, offenen bis mäßig engen Großfaltengenerationen im km- Größenbereich besteht anstatt aus den vorher angenommenen Decken- oder Halbfalten. Hauptproblem der Großfaltenstruktur sind Zonen von durchgreifender Boudinage, die in der Nähe der Großfaltenscharniere entstanden sind. Vorherige Bearbeiter deuteten diese Zonen als 'starke Verformungszonen' oder Überschiebungszonen. Diese Arbeit zeigt, dass in diesen Zonen nur durch die asymmetrische Faltung die unteren liegenden Schenkel der Großfalten boudiniert und somit häufig die ansonsten typischen Faltenstrukturen des liegenden Schenkels einer symmetrischen Faltung überprägt wurden. Ein weiteres Problem dieser mesozoischen Großfaltenstruktur ist die Überprägung einer Faltengeneration auf eine frühere. Weil die Verkürzungsrichtung der überprägenden Faltengeneration nicht subparallel zur älteren Faltenachse ist, sondern einen Winkel von rund 30 Grad einschließt, ist ein Wechsel von orthogonalen zu koaxialen Interferenzmustern der Kleinfalten beobachtbar. Folglich ist die Orientierung der Scheitellinie einer überprägenden und überprägten Kleinfalte nicht unbedingt subparallel zur Orientierung der Faltenachse der Großfalte trotz zylindrischer Faltung. Im letzten Teil dieser Arbeit wird die Überprägung der mesozoischen Großfaltenstruktur durch das känozoisch entstandene, transpressionale Alpine Störungssystem, das einen zweiseitigen Falten- und Überschiebungsgürtel im Otago und im Nordwesten anschließenden Alpinen Schiefergürtel bildet, beschrieben.