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.

Nuclear charge radius determination of 7,10 Be and the one-neutron halo nucleus 11 Be

The only nuclear model independent method for the determination of nuclear charge radii of short-lived radioactive isotopes is the measurement of the isotope shift. For light elements (Z < 10) extremely high accuracy in experiment and theory is required and was only reached for He and Li so far. The nuclear charge radii of the lightest elements are of great interest because they have isotopes which exhibit so-called halo nuclei. Those nuclei are characterized by a a very exotic nuclear structure: They have a compact core and an area of less dense nuclear matter that extends far from this core. Examples for halo nuclei are 6^He, 8^He, 11^Li and 11^Be that is investigated in this thesis. Furthermore these isotopes are of interest because up to now only for such systems with a few nucleons the nuclear structure can be calculated ab-initio. In the Institut für Kernchemie at the Johannes Gutenberg-Universität Mainz two approaches with different accuracy were developed. The goal of these approaches was the measurement of the isotope shifts between (7,10,11)^Be^+ and 9^Be^+ in the D1 line. The …first approach is laser spectroscopy on laser cooled Be^+ ions that are trapped in a linear Paul trap. The accessible accuracy should be in the order of some 100 kHz. In this thesis two types of linear Paul traps were developed for this purpose. Moreover, the peripheral experimental setup was simulated and constructed. It allows the efficient deceleration of fast ions with an initial energy of 60 keV down to some eV and an effcient transport into the ion trap. For one of the Paul traps the ion trapping could already be demonstrated, while the optical detection of captured 9^Be^+ ions could not be completed, because the development work was delayed by the second approach. The second approach uses the technique of collinear laser spectroscopy that was already applied in the last 30 years for measuring isotope shifts of plenty of heavier isotopes. For light elements (Z < 10), it was so far not possible to reach the accuracy that is required to extract information about nuclear charge radii. The combination of collinear laser spectroscopy with the most modern methods of frequency metrology …finally permitted the …first-time determination of the nuclear charge radii of (7,10)^Be and the one neutron halo nucleus 11^Be at the COLLAPS experiment at ISOLDE/ CERN. In the course of the work reported in this thesis it was possible to measure the absolute transition frequencies and the isotope shifts in the D1 line for the Be isotopes mentioned above with an accuracy of better than 2 MHz. Combination with the most recent calculations of the mass effect allowed the extraction of the nuclear charge radii of (7,10,11)^Be with an relative accuracy better than 1%. The nuclear charge radius decreases from 7^Be continuously to 10^Be and increases again for 11^Be. This result is compared with predictions of ab-initio nuclear models which reproduce the observed trend. Particularly the "Greens Function Monte Carlo" and the "Fermionic Molecular Dynamic" model show very good agreement.

Solid-state NMR investigation of phosphonic acid based proton conducting materials

In this work, new promising proton conducting fuel cell membrane materials were characterized in terms of their structure and dynamic properties using solid-state nuclear magnetic resonance (NMR) spectroscopy and X-ray diffraction. Structurally different, phosphonic acid (PA) containing materials were systematically evaluated for possible high-temperature operation (e.g. at T>100°C). Notably, 1H, 2H and 31P magic angle spinning (MAS) NMR provided insight into local connectivities and dynamics of the hydrogen bonded network, while packing arrangements were identified by means of heteronuclear dipolar recoupling techniques.rnThe first part of this work introduced rather crystalline, low molecular weight ionomers for proton conducting membranes, where six different geometries such as line, triangle, screw, tetrahedron, square and hexagon, were investigated. The hexagon was identified as the most promising geometry with high-temperature bulk proton conductivities in the range of 10-3 Scm-1 at a relative humidity of 50%. However, 2H NMR and TGA-MS data suggest that the bulk proton transport is mainly due to the presence of crystal water. Single crystal X-ray data revealed that in the tetrahedron phosphonic acids form tetrameric clusters isolating the mobile protons while the phosphonic acids in the hexagon form zigzag-type pathways through the sample.rnThe second part of this work demonstrates how acid-base pairing and the choice of appropriate spacers may influence proton conduction. Different ratios of statistical copolymers of poly (vinylphosphonic acid) and poly (4-vinylpyridine) were measured to derive information about the local structure and chemical changes. Though anhydrous proton conductivities of all statistical copolymers are rather poor, the conductivity increases to 10-2 S cm-1 when exposing the sample to relative humidity of 80%. In contrast to PVPA, anhydride formation of phosphonic acids in the copolymer is not reversible even when exposing the sample to a relative humidity of 100%.rnIn addition, the influence of both spacers and degree of backbone crystallinity on bulk proton conductivity was investigated. Unlike in systems such as poly benzimidazole (PBI), spacers were inserted between the protogenic groups along the backbone. It was found that dilution of the protogenic groups decreases the conductivity, but compared to PVPA, similar apparent activation energies for local motions were obtained from both variable temperature 1H NMR and impedance spectroscopy data. These observations suggest the formation of phosphonic acid clusters with high degrees of local proton motion, where only a fraction of motions contribute to the observable bulk proton conductivity. Additionally, it was shown that gradual changes of the spacer length lead to different morphologies.rnIn summary, applying advanced solid-state NMR and X-ray analysis, structural and dynamic phenomena in proton conducting materials were identified on a molecular level. The results were discussed with respect to different proton conduction mechanisms and may contribute to a more rational design or improvement of proton conducting membranes.rn

Solid-state NMR on columnar supramolecular structures

The present work deals with the characterisation of three columnar self-assembled systems, that is, benzene-1,3,5-tricarboxamides, a peripherally thioalkyl-substituted phthalocyanine, and several oligo-(p-phenylenevinylene)s. In order to probe the supramolecular organisation solid-state NMR has been used as the main technique, supported by X-ray measurements, theoretical methods, and thermal analysis. rnrnBenzene-1,3,5-tricarboxamides (BTAs) turned out to be well suited model compounds to study various fundamental supramolecular interactions, such as π-π-interactions, hydrogen bonding, as well as dynamic and steric effects of attached side chains. Six BTAs have been investigated in total, five with a CO-centred amide group bearing different side chains and one with an inverted N-centred amide group. The physical properties of these BTAs have been investigated as a function of temperature. The results indicated that in case of the CO-centred BTAs the stability of the columnar mesophase depends strongly on the nature of the side chains. Further experiments revealed a coplanar orientation of adjacent BTA molecules in the columnar assembly of CO-centred BTAs, whereas the N-centred BTA, showed a deviating not fully coplanar arrangement. These differences were ascribed to distinct hydrogen bonding schemes, involving a parallel alignment of hydrogen bonds in case of CO-centred BTAs and an antiparallel alignment in case of the N-centred counterpart.rnrn The fundamental insights of the supramolecular organisation of BTAs could be partially adapted to an octa-substituted phthalocyanine with thiododecyl moieties. Solid-state NMR in combination with chemical shift calculations determined a tilted herringbone arrangement of phthalocyanine rings in the crystalline phase as well as in the mesophase. Moreover, 1H NMR measurements in the mesophase of this compound suggested an axial rotation of molecules, which is inhibited in the crystalline phase.rnrnAs a third task, the supramolecular assembly of oligo-(p-phenylenevinylene)s of varying length and with different polar head groups have been investigated by a combined X-ray and solid-state NMR study. The results revealed a columnar structure formation of these compounds, being promoted by phase separation of alkyl side chains and aromatic rigid rods. In this system solid-state NMR yielded meaningful insight into the isotropisation process of butoxy and 2-S-methylbutoxy substituted oligo-(p-phenylenevinylene) rods.rn

Interplay of pi–pi stacking and hydrogen bonding in conjugated supramolecular systems studied by solid-state NMR

In this work supramolecular organic systems based on rigid pi-conjugated building blocks and flexible side chains were studied via solid-state NMR spectroscopy. Specifically, these studies focussed on phenylene ethynylene based macrocycles, polymer systems including polythiophenes, and rod-coil copolymers of oligo(p-benzamide) and poly(ethylene glycol). All systems were studied in terms of the local order and mobility. The central topic of this dissertation was to elucidate the role of the flexible side chains in interplay of different non-covalent interactions, like pi-pi-stacking and hydrogen bonding.Combining the results of this work, it can be concluded that the ratio of the rigid block and the attached alkyl side chains can be crucial for the design of an ordered pi-conjugated supramolecular system. Through alkyl side chains, it is also possible to introduce liquid-crystalline phases in the system, which can foster the local order of the system. Moreover in the studied system longer, unbranched alkyl side chains are better suited to stabilize the corresponding aggregation than shorter, branched ones.The combination of non-covalent interactions such as pi-pi-stacking and hydrogen bonding play an important role for structure formation. However, the effect of pi-pi-stacking interaction is much weaker than the effect of hydrogen bonding and is only observed in systems with a suitable local order. Hence, they are often not strong enough to control the local order. In contrast, hydrogen bonds predominantly influence the structural organization and packing. In comparison the size of the alkyl side chains is only of minor importance. The suppression of certain hydrogen bonds can lead to completely different structures and can induce a specific aggregation behavior. Thus, for the design of a supramolecular ordered system the presence of hydrogen bonding efficiently stabilizes the corresponding structure, but the ratio of hydrogen bond forming groups should be kept low to be able to influence the structure selectively.

Structure formation as studied by EPR spectroscopy: From simple solutions to nanoparticles

In this thesis, three nitroxide based ionic systems were used to investigate structure and dynamics of their respective solutions in mixed solvents by means of electron paramagnetic resonance (EPR) and electron nuclear double resonance (ENDOR) spectroscopy at X- and W-band (9.5 and 94.5 GHz, respectively). rnFirst, the solvation of the inorganic radical Fremy’s salt (K2ON(SO3)2) in isotope substituted binary solvent mixtures (methanol/water) was investigated by means of high-field (W-band) pulse ENDOR spectroscopy and molecular dynamics (MD) simulations. From the analysis of orientation-selective 1H and 2H ENDOR spectra the principal components of the hyperfine coupling (hfc) tensor for chemically different protons (alcoholic methyl vs. exchangeable protons) were obtained. The methyl protons of the organic solvent approach with a mean distance of 3.5 Å perpendicular to the approximate plane spanned by ON(S)2 of the probe molecule. Exchangeable protons were found to be distributed isotropically, approaching closest to Fremy’s salt from the hydrogen-bonded network around the sulfonate groups. The distribution of exchangeable and methyl protons as found in MD simulations is in full agreement with the ENDOR results. The solvation was found to be similar for the studied solvent ratios between 1:2.3 and 2.3:1 and dominated by an interplay of H-bond (electrostatic) interactions and steric considerations with the NO group merely involved into H-bonds.rnFurther, the conformation of spin labeled poly(diallyldimethylammonium chloride) (PDADMAC) solutions in aqueous alcohol (methanol, ethanol, n-propanol, ethylene glycol, glycerol) mixtures in dependence of divalent sodium sulfate was investigated with double electron-electron resonance (DEER) spectroscopy. The DEER data was analyzed using the worm-like chain model which suggests that in organic-water solvent mixtures the polymer backbones are preferentially solvated by the organic solvent. We found a less serve impact on conformational changes due to salt than usually predicted in polyelectrolyte theory which stresses the importance of a delicate balance of hydrophobic and electrostatic interactions, in particular in the presence of organic solvents.rnFinally, the structure and dynamics of miniemulsions and polymerdispersions prepared with anionic surfactants, that were partially replaced by a spin labeled fatty acid in presence and absence of a lanthanide beta-diketonate complex was characterized by CW EPR spectroscopy. Such miniemulsions form multilayers with the surfactant head group bound to the lanthanide ion. Beta-diketonates were formerly used as NMR shift reagents and nowadays find application as luminescent materials in OLEDs and LCDs and as contrast agent in MRT. The embedding of the complex into a polymer matrix results in an easy processable material. It was found that the structure formation takes place in miniemulsion and is preserved during polymerization. For surfactants with carboxyl-head group a higher order of the alkyl chains and less lateral diffusion is found than for sulfat-head groups, suggesting a more uniform and stronger coordination to the metal ion. The stability of these bilayers depends on the temperature and the used surfactant which should be considered for the used polymerization temperature if a maximum output of the structured regions is wished.

Tailoring molecular weight of the conjugated model polymer MEH-PPV and its structure-property relations

In der vorliegenden Arbeit wurden Struktur-Eigenschaftsbeziehungen des konjugierten Modell-Polymers MEH-PPV untersucht. Dazu wurde Fällungs-fraktionierung eingesetzt, um MEH-PPV mit unterschiedlichem Molekulargewicht (Mw) zu erhalten, insbesondere MEH-PPV mit niedrigem Mw, da dieses für optische Wellenleiterbauelemente optimal geeignet ist Wir konnten feststellen, dass die Präparation einer ausreichenden Menge von MEH-PPV mit niedrigem Mw und geringer Mw-Verteilung wesentlich von der geeigneten Wahl des Lösungsmittels und der Temperatur während der Zugabe des Fällungsmittels abhängt. Alternativ dazu wurden UV-induzierte Kettenspaltungseffekte untersucht. Wir folgern aus dem Vergleich beider Vorgehensweisen, dass die Fällungsfraktionierung verglichen mit der UV-Behandlung besser geeignet ist zur Herstellung von MEH-PPV mit spezifischem Mw, da das UV-Licht Kettendefekte längs des Polymerrückgrats erzeugt. 1H NMR and FTIR Spektroskopie wurden zur Untersuchung dieser Kettendefekte herangezogen. Wir konnten außerdem beobachten, dass die Wellenlängen der Absorptionsmaxima der MEH-PPV Fraktionen mit der Kettenlänge zunehmen bis die Zahl der Wiederholeinheiten n  110 erreicht ist. Dieser Wert ist signifikant größer als früher berichtet. rnOptische Eigenschaften von MEH-PPV Wellenleitern wurden untersucht und es konnte gezeigt werden, dass sich die optischen Konstanten ausgezeichnet reproduzieren lassen. Wir haben die Einflüsse der Lösungsmittel und Temperatur beim Spincoaten auf Schichtdicke, Oberflächenrauigkeit, Brechungsindex, Doppelbrechung und Wellenleiter-Dämpfungsverlust untersucht. Wir fanden, dass mit der Erhöhung der Siedetemperatur der Lösungsmittel die Schichtdicke und die Rauigkeit kleiner werden, während Brechungsindex, Doppelbrechung sowie Wellenleiter-Dämpfungsverluste zunahmen. Wir schließen daraus, dass hohe Siedetemperaturen der Lösungsmittel niedrige Verdampfungsraten erzeugen, was die Aggregatbildung während des Spincoatings begünstigt. Hingegen bewirkt eine erhöhte Temperatur während der Schichtpräparation eine Erhöhung von Schichtdicke und Rauhigkeit. Jedoch nehmen Brechungsindex und der Doppelbrechung dabei ab.rn Für die Schichtpräparation auf Glassubstraten und Quarzglas-Fasern kam das Dip-Coating Verfahren zum Einsatz. Die Schichtdicke der Filme hängt ab von Konzentration der Lösung, Transfergeschwindigkeit und Immersionszeit. Mit Tauchbeschichtung haben wir Schichten von MEH-PPV auf Flaschen-Mikroresonatoren aufgebracht zur Untersuchung von rein-optischen Schaltprozessen. Dieses Verfahren erweist sich insbesondere für MEH-PPV mit niedrigem Mw als vielversprechend für die rein-optische Signalverarbeitung mit großer Bandbreite.rn Zusätzlich wurde auch die Morphologie dünner Schichten aus anderen PPV-Derivaten mit Hilfe von FTIR Spektroskopie untersucht. Wir konnten herausfinden, dass der Alkyl-Substitutionsgrad einen starken Einfluss auf die mittlere Orientierung der Polymerrückgrate in dünnen Filmen hat.rn

Determination of effective charge of flexible polyelectrolytes by light scattering

The most important property controlling the physicochemical behaviour of polyelectrolytes and their applicability in different fields is the charge density on the macromolecular chain. A polyelectrolyte molecule in solution may have an effective charge density which is smaller than the actual charge density determined from its chemical structure. In the present work an attempt has been made to quantitatively determine this effective charge density of a model polyelectrolyte by using light scattering techniques. Flexible linear polyelectrolytes with a Poly(2-Vinylpyridine) (2-PVP) backbone are used in the present study. The polyelectrolytes are synthesized by quaternizing the pyridine groups of 2-PVP by ethyl bromide to different quaternization degrees. The effect of the molar mass, degree of quaternization and solvent polarity on the effective charge is studied. The results show that the effective charge does not vary much with the polymer molar mass or the degree of quaternization. But a significant increase in the effective charge is observed when the solvent polarity is increased. The results do not obey the counterion condensation theory proposed by Manning. Based on the very low effective charges determined in this study, a new mechanism for the counterion condensation phenomena from a specific polyelectrolyte-counterion interaction is proposed

Hydrogels for DNA isolation from blood

This dissertation describes the synthesis of surface attached hydrogel biomaterials, characterization of their properties, evaluation of structuring concepts and the investigation of these materials in the isolation of DNA from human whole blood. Photosensitive hydrogel precursor materials on the basis of hydroxyethylmethacrylate (HEMA) were synthesized by free radical polymerization. In order to obtain surface bound hydrogel films, the precursors were deposited on a suitable substrate and subsequently irradatiated with UV - light to accomplish the formation of crosslinks in the film and create surface attachment. The composition of the polymerization precursor materials was determined by comprehensive NMR and GPC studies, revealing the copolymerizationrnbehaviour of the used monomers - HEMA derivatives and the photocrosslinkerrnMABP - and their respective distribution in the hydrogel precursors. The degree of crosslinking of the hydrogels was characterized with UV/vis spectroscopy. Stress-strain measurements were conducted in order to investigate the mechanical properties of the biomaterials. Moreover, the swelling process and biomolecule adsorption properties of the hydrogels were investigated with SPR/OW spectroscopy. For this, the deposition and binding of the hydrogels on gold or SiO2 surfaces was facilitated with photocrosslinkable adhesion promotors. The produced hydrogels were mechanically rigid and stablernunder the conditions of PCR and blood lysis. Furthermore, strategies towards the increase of hydrogel surface structure and porosity with porosigens, 2D laser interference lithography and photocleavable blockcopolymers were investigated. At last, a combinatorial strategy was used for the determination of the usefulness of hydrogels for the isolation from DNA from blood. A series of functionalized hydrogel precursors were synthesized, transferred to the surface inside a PCR tube and subsequently screened in regard to DNA adsorption properties with Taqman quantitative PCR. This approach yielded a promising candidate for a functional PCR tube coating that would allow the entire DNA isolation procedure being carried out in a single reaction container.rnThereforce, the practical application of such macromolecular architectures can be envisioned to improve industrial DNA diagnostic processes.

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.