Atomistic simulations of cosolvent/water mixtures

In this thesis, atomistic simulations are performed to investigate hydrophobic solvation and hydrophobic interactions in cosolvent/water binary mixtures. Many cosolvent/water binary mixtures exhibit non-ideal behavior caused by aggregation at the molecular scale level although they are stable and homogenous at the macroscopic scale. Force-field based atomistic simulations provide routes to relate atomistic-scale structure and interactions to thermodynamic solution properties. The predicted solution properties are however sensitive to the parameters used to describe the molecular interactions. In this thesis, a force field for tertiary butanol (TBA) and water mixtures is parameterized by making use of the Kirkwood-Buff theory of solution. The new force field is capable of describing the alcohol-alcohol, water-water and alcohol-water clustering in the solution as well as the solution components’ chemical potential derivatives in agreement with experimental data. With the new force field, the preferential solvation and the solvation thermodynamics of a hydrophobic solute in TBA/water mixtures have been studied. First, methane solvation at various TBA/water concentrations is discussed in terms of solvation free energy-, enthalpy- and entropy- changes, which have been compared to experimental data. We observed that the methane solvation free energy varies smoothly with the alcohol/water composition while the solvation enthalpies and entropies vary nonmonotonically. The latter occurs due to structural solvent reorganization contributions which are not present in the free energy change due to exact enthalpy-entropy compensation. It is therefore concluded that the enthalpy and entropy of solvation provide more detailed information on the reorganization of solvent molecules around the inserted solute. Hydrophobic interactions in binary urea/water mixtures are next discussed. This system is particularly relevant in biology (protein folding/unfolding), however, changes in the hydrophobic interaction induced by urea molecules are not well understood. In this thesis, this interaction has been studied by calculating the free energy (potential of mean force), enthalpy and entropy changes as a function of the solute-solute distance in water and in aqueous urea (6.9 M) solution. In chapter 5, the potential of mean force in both solution systems is analyzed in terms of its enthalpic and entropic contributions. In particular, contributions of solvent reorganization in the enthalpy and entropy changes are studied separately to better understand what are the changes in interactions in the system that contribute to the free energy of association of the nonpolar solutes. We observe that in aqueous urea the association between nonpolar solutes remains thermodynamically favorable (i.e., as it is the case in pure water). This observation contrasts a long-standing belief that clusters of nonpolar molecules dissolve completely in the presence of urea molecules. The consequences of our observations for the stability of proteins in concentrated urea solutions are discussed in the chapter 6 of the thesis.

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.

Sintesi, caratterizzazione e reattività di nuovi Cluster Carbonilici Ni-P

The aim of this work is to explore the chemistry of new heteroatomic molecular compounds never reported in the scientific literature so far: Ni-P carbonyl clusters. First, an attempt was made to illustrate the reasons which brought to the choice of this specific metal-pnictogen couple, such as the interesting properties of Ni-P binary phases, the difficulties related in obtaining structural data for this bulk compounds and the absence of references about nickel-phosphorus molecular compounds (e.g. carbonyl clusters) reported in literature. Then, the general criteria chosen for the reactions between precursors [Ni6(CO)12]2- and PCl3 or POCl3 have been reported. This work has permitted to individuate many new products, of which some have also been isolated and characterised: [Ni11P(CO)18]3-, [Ni23-xP2(CO)30-x]6- (x=0, 1), [HNi31P4(CO)39]5- e [H2Ni31P4(CO)39]4-. Except for the former, a reproducible synthetic path has been refined for all those new Ni-P carbonyl clusters; furthermore some chemical reactivity has been carried out in order to test their characteristics.

Störungstheoretische Behandlung relativistischer Korrekturen zu Gleichgewichtsstrukturen im Rahmen der Coupled-Cluster-Theorie

Quantenchemische Untersuchungen von Atomen und Molekülen haben in den letzten Jahren durch die systematische Erweiterung der Methoden und Computerresourcen zunehmend für die Interpretation und Vorhersage experimenteller Ergebnisse an Bedeutung gewonnen. Relativistische Effekte in der Chemie werden zum Beispiel für die gelbe Farbe von Gold und den flüssigen Aggregatzustand von Quecksilber verantwortlich gemacht und müssen daher in quantenchemischen Rechnungen berücksichtigt werden. Relativistische Effekte sind bei leichten Elementen oft so klein, daß sie in vielen quantenchemischen Betrachtungen vernachlässigt werden. Dennoch sind es gerade diese Beiträge, die verbleibende Abweichungen von noch so genauen nichtrelativistischen Rechnungen von ebenso genauen experimentellen Ergebnissen ausmachen können. Relativistische Effekte können auf viele Arten in quantenchemischen Rechnungen berücksichtigt werden. Eine Möglichkeit ist die Störungstheorie. Ein derartiger Ansatz ist die Mass-velocity-Darwin-Näherung, ein anderer die Direkte Störungstheorie. Hier entspricht die relativistische Energiekorrektur erster Ordnung der ersten Ableitung der Energie nach einem relativistischen Störparameter. Für eine Bestimmung der Gleichgewichtsstruktur eines Moleküls müssen die Kräfte auf die Atomkerne bestimmt werden. Diese entsprechen einer ersten Ableitung der Gesamtenergie nach den Kernkoordinaten. Eine Einbeziehung der relativistischen Effekte auf diese Kräfte erfordert daher die gemischte zweite Ableitung der Energie nach dem relativistischen Störparameter und den Kernkoordinaten. Diese relativistischen Korrekturen wurden in dem quantenchemischen Programmpaket ACES2 implementiert. Ein Resultat dieser Arbeit ist, daß nun erstmalig eine Implementierung analytischer Gradienten für die Berechnung relativistischer Korrekturen zu Strukturparametern mit Hilfe der relativistischen Störungstheorie für den Coupled-Cluster-Ansatz bereit steht. Die Coupled-Cluster-Theorie eignet sich besonders gut für die hochgenaue Vorhersage von molekularen Eigenschaften, wie der Gleichgewichtsstruktur. Im Rahmen dieser Arbeit wurde die Basissatzabhängigkeit der relativistischen Beiträge zu Energien, Strukturparametern und harmonischen Schwingungsfrequenzen im Detail untersucht. Für die hier untersuchten Moleküle sind die relativistischen Effekte und Effekte aufgrund der Elektronenkorrelation nicht additiv, so verkürzt die Berücksichtigung relativistischer Effekte bei Hartree-Fock-Rechnungen die Bindung in den Hydrogenhalogeniden, während die Einbeziehung der Elektronenkorrelation durch CCSD(T)-Rechnungen zu einer verlängerten Bindung im Fluorwasserstoff und weniger stark ausgeprägten Korrekturen im Chlor- und Bromwasserstoff führt. Für die anderen hier untersuchten mehratomigen Moleküle findet sich kein einheitlicher Trend; dies unterstreicht die Notwendigkeit expliziter Rechnungen. Damit steht ein leistungsfähiges und vielseitiges Werkzeug für die Berechnung relativistischer Korrekturen auf verschiedenste molekulare Eigenschaften zur Verfügung, das mit modernen, systematisch verbesserbaren quantenchemischen Methoden verknüpft ist. Hiermit ist es möglich, hochgenaue Rechnungen zur Vorhersage und Interpretation von Experimenten durchzuführen.

Studio mediante spettroscopia di assorbimento di raggi x con luce di sincrotrone di cluster di argento su ceria

si vuole studiare la struttura atomica di cluster di argento su ceria, si è utilizzata la spettroscopia di assorbimento di raggi x con luce di sincrotrone. Le misure sono state effettuate all'Esrf di Grenoble e si sono analizzati gli spettri EXAFS di diversi spessori equivalenti di argento (1.5Å , 0.7Å, 0.2Å, 0.1Å). Si tratta il processo di produzione della luce di sincrotrone e le sue proprietà caratteristiche. In particolare si descrive la strumentazione ottica utilizzata: specchi riflettenti per raggi X e monocromatore con geometria a doppio cristallo.

New developments in state-specific multireference coupled-cluster theory

Coupled-cluster theory in its single-reference formulation represents one of the most successful approaches in quantum chemistry for the description of atoms and molecules. To extend the applicability of single-reference coupled-cluster theory to systems with degenerate or near-degenerate electronic configurations, multireference coupled-cluster methods have been suggested. One of the most promising formulations of multireference coupled cluster theory is the state-specific variant suggested by Mukherjee and co-workers (Mk-MRCC). Unlike other multireference coupled-cluster approaches, Mk-MRCC is a size-extensive theory and results obtained so far indicate that it has the potential to develop to a standard tool for high-accuracy quantum-chemical treatments. This work deals with developments to overcome the limitations in the applicability of the Mk-MRCC method. Therefore, an efficient Mk-MRCC algorithm has been implemented in the CFOUR program package to perform energy calculations within the singles and doubles (Mk-MRCCSD) and singles, doubles, and triples (Mk-MRCCSDT) approximations. This implementation exploits the special structure of the Mk-MRCC working equations that allows to adapt existing efficient single-reference coupled-cluster codes. The algorithm has the correct computational scaling of d*N^6 for Mk-MRCCSD and d*N^8 for Mk-MRCCSDT, where N denotes the system size and d the number of reference determinants. For the determination of molecular properties as the equilibrium geometry, the theory of analytic first derivatives of the energy for the Mk-MRCC method has been developed using a Lagrange formalism. The Mk-MRCC gradients within the CCSD and CCSDT approximation have been implemented and their applicability has been demonstrated for various compounds such as 2,6-pyridyne, the 2,6-pyridyne cation, m-benzyne, ozone and cyclobutadiene. The development of analytic gradients for Mk-MRCC offers the possibility of routinely locating minima and transition states on the potential energy surface. It can be considered as a key step towards routine investigation of multireference systems and calculation of their properties. As the full inclusion of triple excitations in Mk-MRCC energy calculations is computational demanding, a parallel implementation is presented in order to circumvent limitations due to the required execution time. The proposed scheme is based on the adaption of a highly efficient serial Mk-MRCCSDT code by parallelizing the time-determining steps. A first application to 2,6-pyridyne is presented to demonstrate the efficiency of the current implementation.

Development and applications of hollow fiber flow field-flow fractionation in the bioanalytical field. Studies of aggregation phenomena in complex protein samples

Recent advances in the fast growing area of therapeutic/diagnostic proteins and antibodies - novel and highly specific drugs - as well as the progress in the field of functional proteomics regarding the correlation between the aggregation of damaged proteins and (immuno) senescence or aging-related pathologies, underline the need for adequate analytical methods for the detection, separation, characterization and quantification of protein aggregates, regardless of the their origin or formation mechanism. Hollow fiber flow field-flow fractionation (HF5), the miniaturized version of FlowFFF and integral part of the Eclipse DUALTEC FFF separation system, was the focus of this research; this flow-based separation technique proved to be uniquely suited for the hydrodynamic size-based separation of proteins and protein aggregates in a very broad size and molecular weight (MW) range, often present at trace levels. HF5 has shown to be (a) highly selective in terms of protein diffusion coefficients, (b) versatile in terms of bio-compatible carrier solution choice, (c) able to preserve the biophysical properties/molecular conformation of the proteins/protein aggregates and (d) able to discriminate between different types of protein aggregates. Thanks to the miniaturization advantages and the online coupling with highly sensitive detection techniques (UV/Vis, intrinsic fluorescence and multi-angle light scattering), HF5 had very low detection/quantification limits for protein aggregates. Compared to size-exclusion chromatography (SEC), HF5 demonstrated superior selectivity and potential as orthogonal analytical method in the extended characterization assays, often required by therapeutic protein formulations. In addition, the developed HF5 methods have proven to be rapid, highly selective, sensitive and repeatable. HF5 was ideally suitable as first dimension of separation of aging-related protein aggregates from whole cell lysates (proteome pre-fractionation method) and, by HF5-(UV)-MALS online coupling, important biophysical information on the fractionated proteins and protein aggregates was gathered: size (rms radius and hydrodynamic radius), absolute MW and conformation.

Cluster carbonilici eterometallici di rodio: Sintesi di nuove specie molecolari e studio della loro reattività

Durante il mio periodo di Tesi ho lavorato nel campo dei cluster eterometallici di Rh stabilizzati da leganti carbonilici. Inizialmente mi sono occupato della reazione di idrogenazione del cluster [Rh12Sn(CO)27]4- e della specie insatura [Rh12Sn(CO)25]4-, che ho sintetizzato per la prima volta con reattore a microonde. Ho poi testato il prodotto [Rh12Sn(CO)27/25-xHy]4- come catalizzatore nella reazione di idrogenazione di stirene ad etilbenzene. Le prove hanno evidenziato la sua scarsa attività catalitica, anzi il catalizzatore si è rivelato instabile sia in atmosfera di H2, dove si trasforma in [Rh14(CO)25]4-, che di CO/H2, dove genera una specie ancora non ben caratterizzata. Nella seconda parte del mio progetto ho cercato di sintetizzare nuovi cluster carbonilici Rh-Bi. Dopo alcune prove sono riuscito ad isolare e caratterizzare il nuovo cluster icosaedrico [Rh12Bi(CO)27]3-. Ne ho poi studiato la reattività trattandolo con BiCl3, HCl e irraggiandolo con le microonde. Dalla prima reazione ho ottenuto due nuove specie molecolari cioè il dimero [(Rh12Bi(CO)26)2Bi]5-, in cui due unità icosaedriche sono legate insieme da uno ione Bi+, e [Rh14Bi3(CO)27]3-, in cui sono presenti legami incipienti intermolecolari Bi-Bi. Nelle altre due prove sono stati identificati nuovi prodotti non ancora identificati ma tuttora in fase di studio.

Aggregation of hybrid molecules and film formation of colloidal monolayers

A unique characteristic of soft matter is its ability to self-assemble into larger structures. Characterizing these structures is crucial for their applications. In the first part of this work, I investigated DNA-organic hybrid material by means of Fluorescence Correlation Spectroscopy (FCS) and Fluorescence Cross-Correlation Spectroscopy (FCCS). DNA-organic hybrid materials, a novel class of hybrid materials composed of synthetic macromolecules and oligodeoxynucleotide segmenta, are mostly amphiphilic and can self-assemble into supramolecular structures in aqueous solution. A hybrid material of a fluorophore, perylenediimide (PDI), and a DNA segment (DNA-PDI) has been developed in Prof. A. Hermann’s group (University of Groningen). This novel material has the ability to form aggregates through pi-pi stacking between planar PDIs and can be traced in solution due to the fluorescence of PDI. I have determined the diffusion coefficient of DNA-PDI conjugates in aqueous solution by means of FCS. In addition, I investigated whether such DNA-PDIs form aggregates with certain structure, for instance dimers. rnOnce the DNA hybrid material self-assemble into supermolecular structures for instance into micelles, the single molecules do not necessarily stay in one specific micelle. Actually, a single molecule may enter and leave micelles constantly. The average residence time of a single molecule in a certain micelle depends on the nature of the molecule. I have chosen DNA-b-polypropylene oxide (PPO) as model molecules and investigated the residence time of DNA-b-PPO molecules in their according micelles by means of FCCS.rnBesides the DNA hybrid materials, polymeric colloids can also form ordered structures once they are brought to an air/water interface. Here, hexagonally densely packed monolayers can be generated. These monolayers can be deposited onto different surfaces as coating layers. In the second part of this work, I investigated the mechanical properties of such colloidal monolayers using micromechanical cantilevers. When a coating layer is deposited on a cantilever, it can modify the elasticity of the cantilever. This variation can be reflected either by a deflection or by a resonance frequency shift of the cantilever. In turn, detecting these changes provides information about the mechanical properties of the coating layer. rnIn the second part of this work, polymeric colloidal monolayers were coated on a cantilever and homogenous polymer films of a few hundred nanometers in thickness were generated from these colloidal monolayers by thermal annealing or organic vapor annealing. Both the film formation process and the mechanical properties of these resulting homogenous films were investigated by means of cantilever. rnElastic property changes of the coating film, for example upon absorption of organic vapors, induce a deflection of the cantilever. This effect enables a cantilever to detect target molecules, when the cantilever is coated with an active layer with specific affinity to target molecules. In the last part of this thesis, I investigated the applicability of suitably functionalized micromechanical cantilevers as sensors. In particular, glucose sensitive polymer brushes were grafted on a cantilever and the deflection of this cantilever was measured during exposure to glucose solution. rn

Interazione di cluster carbonilici FE-CU con complessi di metalli di transizione e leganti azotati

Nel corso di questa Tesi sono state studiate reazioni di cluster carburo-carbonilici Fe/Cu con composti azotati di varia natura: complessi fosforescenti di metalli di transizione, 1-10 fenantrolina, L-amminoacidi e Chinolina. In particolare la reazione tra Fe4C(CO)12Cu2(CH3CN)2 e [Ru(tpy)(bpy)(N4C-C6H4-CN)]+ ha portato alla sintesi dell’addotto Fe4C(CO)12{Cu2Cl[Ru(tpy)(bpy)(N4C-C6H4-CN)]}, sul quale sono state condotte misure di luminescenza (emissione, eccitazione e misura dei tempi di vita degli stati eccitati). Per confronto degli spettri registrati su campioni di adotto in soluzione con quelli del complesso cationico di Ru(II), si è ipotizzato che l’addotto sintetizzato in soluzione dia origine ad un sistema in equilibrio tra le specie legate e dissociate. Le reazioni di Fe4C(CO)12Cu2(CH3CN)2 e [NEt4][Fe5C(CO)14Cu(CH3CN)] con 1-10 fenantrolina hanno permesso di isolare le nuove specie [Fe4C(CO)12(Cuphen)]–, [Fe4C(CO)12(Cuphen)] e [Fe5C(CO)14(Cuphen)]–, sottoforma dei loro sali [Cu(phen)2][Fe4C(CO)12(Cuphen)], [NEt4] [Fe4C(CO)12(Cuphen)], [Fe4C(CO)12(Cuphen)], [NEt4][Fe5C(CO)14(Cuphen)]• CH2Cl2 e [NEt4][Fe5C(CO)14(Cuphen)]•THF. In tali cluster si nota come la natura bidentata di phen e il suo ingombro sterico abbiano causato notevoli riarrangiamenti strutturali rispetto alle specie iniziali contenenti acetonitrile. La sintesi di Fe4C(CO)12(CuQ)2 e [NEt4][Fe5C(CO)14(CuQ)] è avvenuta inaspettatamente a partire dalle reazioni condotte tra Fe4C(CO)12Cu2(CH3CN)2 e [NEt4][Fe5C(CO)14Cu(CH3CN)] con le molecole L-prolina, L-metionina e guanina, a causa della chinolina contenuta come impurezza nei reagenti di partenza. L’esito di questa reazione ha comunque mostrato l’elevata affinità dei cluster per il legante chinolina, sebbene presente in ambiente di reazione in misura sensibilmente inferiore rispetto agli altri reagenti. Tutte le nuove specie sintetizzate sono stati caratterizzate mediante spettroscopia IR e le strutture molecolari sono state determinate mediante diffrazione di raggi X su cristallo singolo.

Exploring Quantum Speed-up Through Cluster-state Computers

The aim of this thesis is to investigate the nature of quantum computation and the question of the quantum speed-up over classical computation by comparing two different quantum computational frameworks, the traditional quantum circuit model and the cluster-state quantum computer. After an introductory survey of the theoretical and epistemological questions concerning quantum computation, the first part of this thesis provides a presentation of cluster-state computation suitable for a philosophical audience. In spite of the computational equivalence between the two frameworks, their differences can be considered as structural. Entanglement is shown to play a fundamental role in both quantum circuits and cluster-state computers; this supports, from a new perspective, the argument that entanglement can reasonably explain the quantum speed-up over classical computation. However, quantum circuits and cluster-state computers diverge with regard to one of the explanations of quantum computation that actually accords a central role to entanglement, i.e. the Everett interpretation. It is argued that, while cluster-state quantum computation does not show an Everettian failure in accounting for the computational processes, it threatens that interpretation of being not-explanatory. This analysis presented here should be integrated in a more general work in order to include also further frameworks of quantum computation, e.g. topological quantum computation. However, what is revealed by this work is that the speed-up question does not capture all that is at stake: both quantum circuits and cluster-state computers achieve the speed-up, but the challenges that they posit go besides that specific question. Then, the existence of alternative equivalent quantum computational models suggests that the ultimate question should be moved from the speed-up to a sort of “representation theorem” for quantum computation, to be meant as the general goal of identifying the physical features underlying these alternative frameworks that allow for labelling those frameworks as “quantum computation”.

Highly accurate quantum chemistry : spin-orbit splittings via multireference coupled-cluster methods and applications in heavy-atom main-group chemistry

Diese Dissertation demonstriert und verbessert die Vorhersagekraft der Coupled-Cluster-Theorie im Hinblick auf die hochgenaue Berechnung von Moleküleigenschaften. Die Demonstration erfolgt mittels Extrapolations- und Additivitätstechniken in der Single-Referenz-Coupled-Cluster-Theorie, mit deren Hilfe die Existenz und Struktur von bisher unbekannten Molekülen mit schweren Hauptgruppenelementen vorhergesagt wird. Vor allem am Beispiel von cyclischem SiS_2, einem dreiatomigen Molekül mit 16 Valenzelektronen, wird deutlich, dass die Vorhersagekraft der Theorie sich heutzutage auf Augenhöhe mit dem Experiment befindet: Theoretische Überlegungen initiierten eine experimentelle Suche nach diesem Molekül, was schließlich zu dessen Detektion und Charakterisierung mittels Rotationsspektroskopie führte. Die Vorhersagekraft der Coupled-Cluster-Theorie wird verbessert, indem eine Multireferenz-Coupled-Cluster-Methode für die Berechnung von Spin-Bahn-Aufspaltungen erster Ordnung in 2^Pi-Zuständen entwickelt wird. Der Fokus hierbei liegt auf Mukherjee's Variante der Multireferenz-Coupled-Cluster-Theorie, aber prinzipiell ist das vorgeschlagene Berechnungsschema auf alle Varianten anwendbar. Die erwünschte Genauigkeit beträgt 10 cm^-1. Sie wird mit der neuen Methode erreicht, wenn Ein- und Zweielektroneneffekte und bei schweren Elementen auch skalarrelativistische Effekte berücksichtigt werden. Die Methode eignet sich daher in Kombination mit Coupled-Cluster-basierten Extrapolations-und Additivitätsschemata dafür, hochgenaue thermochemische Daten zu berechnen.

Internally contracted multireference coupled-cluster methods

This thesis details the development of quantum chemical methods for the accurate theoretical description of molecular systems with a complicated electronic structure. In simple cases, a single Slater determinant, in which the electrons occupy a number of energetically lowest molecular orbitals, offers a qualitatively correct model. The widely used coupled-cluster method CCSD(T) efficiently includes electron correlation effects starting from this determinant and provides reaction energies in error by only a few kJ/mol. However, the method often fails when several electronic configurations are important, as, for instance, in the course of many chemical reactions or in transition metal compounds. Internally contracted multireference coupled-cluster methods (ic-MRCC methods) cure this deficiency by using a linear combination of determinants as a reference function. Despite their theoretical elegance, the ic-MRCC equations involve thousands of terms and are therefore derived by the computer. Calculations of energy surfaces of BeH2, HF, LiF, H2O, N2 and Be3 unveil the theory's high accuracy compared to other approaches and the quality of various hierarchies of approximations. New theoretical advances include size-extensive techniques for removing linear dependencies in the ic-MRCC equations and a multireference analog of CCSD(T). Applications of the latter method to O3, Ni2O2, benzynes, C6H7NO and Cr2 underscore its potential to become a new standard method in quantum chemistry.

Higher-order molecular properties and excitation energies in single-reference and multireference coupled-cluster theory

Coupled-cluster (CC) theory is one of the most successful approaches in high-accuracy quantum chemistry. The present thesis makes a number of contributions to the determination of molecular properties and excitation energies within the CC framework. The multireference CC (MRCC) method proposed by Mukherjee and coworkers (Mk-MRCC) has been benchmarked within the singles and doubles approximation (Mk-MRCCSD) for molecular equilibrium structures. It is demonstrated that Mk-MRCCSD yields reliable results for multireference cases where single-reference CC methods fail. At the same time, the present work also illustrates that Mk-MRCC still suffers from a number of theoretical problems and sometimes gives rise to results of unsatisfactory accuracy. To determine polarizability tensors and excitation spectra in the MRCC framework, the Mk-MRCC linear-response function has been derived together with the corresponding linear-response equations. Pilot applications show that Mk-MRCC linear-response theory suffers from a severe problem when applied to the calculation of dynamic properties and excitation energies: The Mk-MRCC sufficiency conditions give rise to a redundancy in the Mk-MRCC Jacobian matrix, which entails an artificial splitting of certain excited states. This finding has established a new paradigm in MRCC theory, namely that a convincing method should not only yield accurate energies, but ought to allow for the reliable calculation of dynamic properties as well. In the context of single-reference CC theory, an analytic expression for the dipole Hessian matrix, a third-order quantity relevant to infrared spectroscopy, has been derived and implemented within the CC singles and doubles approximation. The advantages of analytic derivatives over numerical differentiation schemes are demonstrated in some pilot applications.

Analytische Ableitungen der Energie für den "Equation-of-Motion-Coupled-Cluster"-Ansatz zur Berechnung von Ionisationspotentialen

In der vorliegenden Arbeit wird die Theorie der analytischen zweiten Ableitungen für die EOMIP-CCSD-Methode formuliert sowie die durchgeführte Implementierung im Quantenchemieprogramm CFOUR beschrieben. Diese Ableitungen sind von Bedeutung bei der Bestimmung statischer Polarisierbarkeiten und harmonischer Schwingungsfrequenzen und in dieser Arbeit wird die Genauigkeit des EOMIP-CCSD-Ansatzes bei der Berechnung dieser Eigenschaften für verschiedene radikalische Systeme untersucht. Des Weiteren können mit Hilfe der ersten und zweiten Ableitungen vibronische Kopplungsparameter berechnet werden, welche zur Simulation von Molekülspektren in Kombination mit dem Köppel-Domcke-Cederbaum (KDC)-Modell - in der Arbeit am Beispiel des Formyloxyl (HCO2)-Radikals demonstriert - benötigt werden.rnrnDer konzeptionell einfache EOMIP-CC-Ansatz wurde gewählt, da hier die Wellenfunktion eines Radikalsystems ausgehend von einem stabilen geschlossenschaligen Zustand durch die Entfernung eines Elektrons gebildet wird und somit die Problematik der Symmetriebrechung umgangen werden kann. Im Rahmen der Implementierung wurden neue Programmteile zur Lösung der erforderlichen Gleichungen für die gestörten EOMIP-CC-Amplituden und die gestörten Lagrange-Multiplikatoren zeta zum Quantenchemieprogramm CFOUR hinzugefügt. Die unter Verwendung des Programms bestimmten Eigenschaften werden hinsichtlich ihrer Leistungsfähigkeit im Vergleich zu etablierten Methoden wie z.B. CCSD(T) untersucht. Bei der Berechnung von Polarisierbarkeiten und harmonischen Schwingungsfrequenzen liefert die EOMIP-CCSD-Theorie meist gute Resultate, welche nur wenig von den CCSD(T)-Ergebnissen abweichen. Einzig bei der Betrachtung von Radikalen, für die die entsprechenden Anionen nicht stabil sind (z.B. NH2⁻ und CH3⁻), liefert der EOMIP-CCSD-Ansatz aufgrund methodischer Nachteile keine aussagekräftige Beschreibung. rnrnDie Ableitungen der EOMIP-CCSD-Energie lassen sich auch zur Simulation vibronischer Kopplungen innerhalb des KDC-Modells einsetzen.rnZur Kopplung verschiedener radikalischer Zustände in einem solchen Modellpotential spielen vor allem die Ableitungen von Übergangsmatrixelementen eine wichtige Rolle. Diese sogenannten Kopplungskonstanten können in der EOMIP-CC-Theorie besonders leicht definiert und berechnet werden. Bei der Betrachtung des Photoelektronenspektrums von HCO2⁻ werden zwei Alternativen untersucht: Die vertikale Bestimmung an der Gleichgewichtsgeometrie des HCO2⁻-Anions und die Ermittlung adiabatischer Kraftkonstanten an den Gleichgewichtsgeometrien des Radikals. Lediglich das adiabatische Modell liefert bei Beschränkung auf harmonische Kraftkonstanten eine qualitativ sinnvolle Beschreibung des Spektrums. Erweitert man beide Modelle um kubische und quartische Kraftkonstanten, so nähern sich diese einander an und ermöglichen eine vollständige Zuordnung des gemessenen Spektrums innerhalb der ersten 1500 cm⁻¹. Die adiabatische Darstellung erreicht dabei nahezu quantitative Genauigkeit.