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.

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.

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.

Nahfeld-induzierte Elektronenemissions-Mikrospektroskopie an stark gekoppelten Plasmonen und metallischen Mikrostrukturen

Plasmonen stellen elektromagnetische Moden in metallischen Strukturen dar, in denen die quasifreien Elektronen im Metall kollektiv oszillieren. Während des letzten Jahrzehnts erfuhr das Gebiet der Plasmonik eine rasante Entwicklung, basierend auf zunehmenden Fortschritten der Nanostrukturierungsmethoden und spektroskopischen Untersuchungsmethoden, die zu der Möglichkeit von systematischen Einzelobjektuntersuchungen wohldefinierter Nanostrukturen führte. Die Anregung von Plasmonen resultiert neben einer radiativen Verstärkung der optischen Streuintensität im Fernfeld in einer nicht-radiativen Überhöhung der Feldstärke in unmittelbarer Umgebung der metallischen Struktur (Nahfeld), die durch die kohärente Ladungsansammlung an der metallischen Oberfläche hervorgerufen wird. Das optische Nahfeld stellt folglich eine bedeutende Größe für das fundamentale Verständnis der Wirkung und Wechselwirkung von Plasmonen sowie für die Optimierung plasmonbasierter Applikationen dar. Die große Herausforderung liegt in der Kompliziertheit des experimentellen Zugangs zum Nahfeld, der die Entwicklung eines grundlegenden Verständisses des Nahfeldes verhinderte.rnIm Rahmen dieser Arbeit wurde Photoemissionselektronenmikroskopie (PEEM) bzw. -mikrospektroskopie genutzt, um ortsaufgelöst die Eigenschaften nahfeld-induzierter Elektronenemission zu bestimmen. Die elektrodynamischen Eigenschaften der untersuchten Systeme wurden zudem mit numerischen, auf der Finiten Integrationsmethode basierenden Berechnungen bestimmt und mit den experimentellen Resultaten verglichen.rnAg-Scheiben mit einem Durchmesser von 1µm und einer Höhe von 50nm wurden mit fs-Laserstrahlung der Wellenlänge 400nm unter verschiedenen Polarisationszuständen angeregt. Die laterale Verteilung der infolge eines 2PPE-Prozesses emittierten Elektronen wurde mit dem PEEM aufgenommen. Aus dem Vergleich mit den numerischen Berechnungen lässt sich folgern, dass sich das Nahfeld an unterschiedlichen Stellen der metallischen Struktur verschiedenartig ausbildet. Insbesondere wird am Rand der Scheibe bei s-polarisierter Anregung (verschwindende Vertikalkomponente des elektrischen Felds) ein Nahfeld mit endlicher z-Komponente induziert, während im Zentrum der Scheibe das Nahfeld stets proportional zum einfallenden elektrischen Feld ist.rnWeiterhin wurde erstmalig das Nahfeld optisch angeregter, stark gekoppelter Plasmonen spektral (750-850nm) untersucht und für identische Nanoobjekte mit den entsprechenden Fernfeldspektren verglichen. Dies erfolgte durch Messung der spektralen Streucharakteristik der Einzelobjekte mit einem Dunkelfeldkonfokalmikroskop. Als Modellsystem stark gekoppelter Plasmonen dienten Au Nanopartikel in sub-Nanometerabstand zu einem Au Film (nanoparticle on plane, NPOP). Mit Hilfe dieser Kombination aus komplementären Untersuchungsmethoden konnte erstmalig die spektrale Trennung von radiativen und nicht-radiativen Moden stark gekoppelter Plasmonen nachgewiesen werden. Dies ist insbesondere für Anwendungen von großer Relevanz, da reine Nahfeldmoden durch den unterdrückten radiativen Zerfall eine große Lebensdauer besitzen, so dass deren Verstärkungswirkung besonders lange nutzbar ist. Ursachen für die Unterschiede im spektralen Verhalten von Fern- und Nahfeld konnten durch numerische Berechnungen identifiziert werden. Sie zeigten, dass das Nahfeld nicht-spärischer NPOPs durch die komplexe Oszillationsbewegung der Elektronen innerhalb des Spaltes zwischen Partikel und Film stark ortsabhängig ist. Zudem reagiert das Nahfeld stark gekoppelter Plasmonen deutlich empfindlicher auf strukturelle Störstellen des Resonators als die Fernfeld-Response. Ferner wurde der Elektronenemissionsmechanismus als optischer Feldemissionsprozess identifiziert. Um den Vorgang beschreiben zu können, wurde die Fowler-Nordheim Theorie der statischen Feldemission für den Fall harmonisch oszillierender Felder modifiziert.

Probing the function of ionotropic and G protein-coupled receptors in surface-confined membranes

This article reports on recent electrical and optical techniques for investigating cellular signaling reactions in artificial and native membranes immobilized on solid supports. The first part describes the formation of planar artificial lipid bilayers on gold electrodes, which reveal giga-ohm electrical resistance and the insertion and characterization of ionotropic receptors therein. These membranes are suited to record a few or even single ion channels by impedance spectroscopy. Such tethered membranes on planar arrays of microelectrodes offer mechanically robust, long-lasting measuring devices to probe the influence of different chemistries on biologically important ionotropic receptors and therefore will have a future impact to probe the function of channel proteins in basic science and in biosensor applications. In a second part, we present complementary approaches to form inside-out native membrane sheets that are immobilized on micrometer-sized beads or across submicrometer-sized holes machined in a planar support. Because the native membrane sheets are plasma membranes detached from live cells, these approaches offer a unique possibility to investigate cellular signaling processes, such as those mediated by ionotropic or G protein-coupled receptors, with original composition of lipids and proteins.

A novel open channel blocker of GABA-A receptors

GABA-A receptors are chloride ion channels composed of five subunits, mediating fast synaptic and tonic inhibition in the mammalian brain. 19 different subunit isoforms have been identified, with the major receptor type in mammalian adult brain consisting of α1, β2, and γ2 subunits. GABA-A receptors are the target of numerous sedating and anxiolytic drugs such as benzodiazepines. The currently known endogenous ligands are GABA, neurosteroids and the endocannabinoid 2- arachidonoyl glycerol (2-AG). The pharmacological properties of this chloride ion channel strictly depend on receptor subunit composition and arrangement. GABA-A receptors bind and are inhibited by epileptogenic agents such as picrotoxin, and cyclodiene insecticides such as dieldrin. We screened aromatic monovalent anions with five-fold symmetry for inhibition of GABA-A receptors. One of the anions, PCCPinhibited currents elicited by GABA with comparable potency as picrotoxin. This inhibition showed all characteristics of an open channel block. The GABA-A receptor ion channel is lined by residues from the M2 membrane-spanning segment. To identify important residues of the pore involved in the interaction with the blocking molecules PCCP-, a mutation scan was performed in combination with subsequent analysis of the expressed mutant proteins using electrophysiological techniques. In a second project we characterised a light-switchable modulator of GABA-A receptors based on propofol. It was my responsibility to investigate the switching kinetics in patch clamp experiments. After its discovery in 1980, propofol has become the most widely used intravenous general anaesthetic. It is commonly accepted that the anaesthesia induced by this unusually lipophilic drug mostly results from potentiation of GABA induced currents. While GABA-A receptors respond to a variety of ligands, they are normally not sensitive towards light. This light sensitivity could be indirectly achieved by using modulators that can be optically switched between an active and an inactive form. We tested an azobenzene derivative of propofol where an aryldiazene unit is directly coupled to the pharmacophore. This molecule was termed azopropofol (AP2). The effect of AP2 on Cl- currents was investigated with electrophysiological techniques using α1β2γ2 GABA-A receptors expressed in Xenopus oocytes and HEK-cells. In the third project we wanted to investigate the functional role of GABA-A receptors in the liver, and their possible involvement in cell proliferation. GABA-A receptors are also found in a wide range of peripheral tissues, including parts of the peripheral nervous system and non-neural tissues such as smooth muscle, the female reproductive system, liver and several cancer tissues. However their precise function in non neuronal or cancerous cells is still unknown. For this purpose we investigated expression, localization and function of the hepatocytes GABA-A receptors in model cell lines and healthy and cancerous hepatocytes.

Nonequilibrium hysteresis and spin relaxation in the mixed-anisotropy dipolar-coupled spin-glass LiHo0.5Er0.5F4

We present a study of the model spin-glass LiHo0.5Er0.5F4 using simultaneous ac susceptibility, magnetization, and magnetocaloric effect measurements along with small angle neutron scattering (SANS) at sub-Kelvin temperatures. All measured bulk quantities reveal hysteretic behavior when the field is applied along the crystallographic c axis. Furthermore, avalanchelike relaxation is observed in a static field after ramping from the zero-field-cooled state up to 200–300 Oe. SANS measurements are employed to track the microscopic spin reconfiguration throughout both the hysteresis loop and the related relaxation. Comparing the SANS data to inhomogeneous mean-field calculations performed on a box of one million unit cells provides a real-space picture of the spin configuration. We discover that the avalanche is being driven by released Zeeman energy, which heats the sample and creates positive feedback, continuing the avalanche. The combination of SANS and mean-field simulations reveal that the conventional distribution of cluster sizes is replaced by one with a depletion of intermediate cluster sizes for much of the hysteresis loop.

Accurate rotational constant and bond lengths of hexafluorobenzene by femtosecond rotational Raman coherence spectroscopy and ab initio calculations

The gas-phase rotational motion of hexafluorobenzene has been measured in real time using femtosecond (fs) time-resolved rotational Raman coherence spectroscopy (RR-RCS) at T = 100 and 295 K. This four-wave mixing method allows to probe the rotation of non-polar gas-phase molecules with fs time resolution over times up to ∼5 ns. The ground state rotational constant of hexafluorobenzene is determined as B 0 = 1029.740(28) MHz (2σ uncertainty) from RR-RCS transients measured in a pulsed seeded supersonic jet, where essentially only the v = 0 state is populated. Using this B 0 value, RR-RCS measurements in a room temperature gas cell give the rotational constants B v of the five lowest-lying thermally populated vibrationally excited states ν7/8, ν9, ν11/12, ν13, and ν14/15. Their B v constants differ from B 0 by between −1.02 MHz and +2.23 MHz. Combining the B 0 with the results of all-electron coupled-cluster CCSD(T) calculations of Demaison et al. [Mol. Phys.111, 1539 (2013)] and of our own allow to determine the C-C and C-F semi-experimental equilibrium bond lengths r e(C-C) = 1.3866(3) Å and r e(C-F) = 1.3244(4) Å. These agree with the CCSD(T)/wCVQZ r e bond lengths calculated by Demaison et al. within ±0.0005 Å. We also calculate the semi-experimental thermally averaged bond lengths r g(C-C)=1.3907(3) Å and r g(C-F)=1.3250(4) Å. These are at least ten times more accurate than two sets of experimental gas-phase electron diffraction r g bond lengths measured in the 1960s.

Fiducial and differential cross sections of Higgs boson production measured in the four-lepton decay channel in pp collisions at √s=8 TeV with the ATLAS detector

Measurements of fiducial and differential cross sections of Higgs boson production in the H →ZZ* → 4ℓ decay channel are presented. The cross sections are determined within a fiducial phase space and corrected for detection efficiency and resolution effects. They are based on 20.3 fb−1 of pp collision data, produced at √s = 8 TeV centre-of-mass energy at the LHC and recorded by the ATLAS detector. The differential measurements are performed in bins of transverse momentum and rapidity of the four-lepton system, the invariant mass of the subleading lepton pair and the decay angle of the leading lepton pair with respect to the beam line in the four-lepton rest frame, as well as the number of jets and the transverse momentum of the leading jet. The measured cross sections are compared to selected theoretical calculations of the Standard Model expectations. No significant deviation from any of the tested predictions is found. c

Chemical analysis and argon geochronology from dive SAR_1 on Cornacya seamount of the Sardinia Channel

Sarcya 1 dive explored a previously unknown 12 My old submerged volcano, labelled Cornacya. A well developed fracturation is characterised by the following directions: N 170 to N-S, N 20 to N 40, N 90 to N 120, N 50 to N 70, which corresponds to the fracturation pattern of the Sardinian margin. The sampled lavas exhibit features of shoshonitic suites of intermediate composition and include amphibole-and mica-bearing lamprophyric xenoliths which are geochemically similar to Ti-poor lamproites. Mica compositions reflect chemical exchanges between the lamprophyre and its shoshonitic host rock suggesting their simultaneous emplacement. Nd compositions of the Cornacya K-rich suite indicate that continental crust was largely involved in the genesis of these rocks. The spatial association of the lamprophyre with the shoshonitic rocks is geochemically similar to K-rich and TiO2-poor igneous suites, emplaced in post-collisional settings. Among shoshonitic rocks, sample SAR 1-01 has been dated at 12.6±0.3 My using the 40Ar/39Ar method with a laser microprobe on single grains. The age of the Cornacya shoshonitic suite is similar to that of the Sisco lamprophyre from Corsica, which similarly is located on the western margin of the Tyrrhenian Sea. Thus, the Cornacya shoshonitic rocks and their lamprophyric xenolith and the Sisco lamprophyre could represent post-collisional suites emplaced during the lithospheric extension of the Corsica-Sardinia block, just after its rotation and before the Tyrrhenian sea opening. Drilling on the Sardinia margin (ODP Leg 107) shows that the upper levels of the present day margin (Hole 654) suffered tectonic subsidence before the lower part (Hole 652). The structure of this lower part is interpreted as the result of an eastward migration of the extension during Late Miocene and Early Pliocene times. Data of Cornacya volcano are in good agreement with this model and provide good chronological constraints for the beginning of the phenomenon.

Chronostratigraphic framework at sites GeoB14403 and GeoB14414 from the Gela Basin, Sicily Channel

A multi-proxy chronological framework along with sequence-stratigraphic interpretations unveils composite Milankovitch cyclicity in the sedimentary records of the Last GlacialeInterglacial cycle at NE Gela Basin on the Sicilian continental margin. Chronostratigraphic data (including foraminifera-based eco-biostratigraphy and d18O records, tephrochronological markers and 14C AMS radiometric datings) was derived from the shallow-shelf drill sites GeoB14403 (54.6 m recovery) and GeoB14414 (27.5 m), collected with both gravity and drilled MeBo cores in 193 m and 146 m water depth, respectively. The recovered intervals record Marine Isotope Stages and Substages (MIS) from MIS 5 to MIS 1, thus comprising major stratigraphic parts of the progradational deposits that form the last 100-ka depositional sequence. Calibration of shelf sedimentary units with borehole stratigraphies indicates the impact of higher-frequency (20-ka) sea level cycles punctuating this 100-ka cycle. This becomes most evident in the alternation of thick interstadial highstand (HST) wedges and thinner glacial forced-regression (FSST) units mirroring seaward shifts in coastal progradation. Albeit their relatively short-lived depositional phase, these subordinate HST units form the bulk of the 100-ka depositional sequence. Two mechanisms are proposed that likely account for enhanced sediment accumulation ratios (SAR) of up to 200 cm/ka during these intervals: (1) intensified activity of deep and intermediate Levantine Intermediate Water (LIW) associated to the drowning of Mediterranean shelves, and (2) amplified sediment flux along the flooded shelf in response to hyperpycnal plumes that generate through extreme precipitation events during overall arid conditions. Equally, the latter mechanism is thought to be at the origin of undulated features resolved in the acoustic records of MIS 5 Interstadials, which bear a striking resemblance to modern equivalents forming on late-Holocene prodeltas of other Mediterranean shallow-shelf settings.

Paleo-sea surface temperature calculations in the equatorial east Atlantic

We present two ~270 kyr paleo-sea surface temperature (SST) records from the Equatorial Divergence and the South Equatorial Current derived from Mg/Ca ratios in the planktic foraminifer Globigerinoides sacculifer. The present study suggests that the magnesium signature of G. sacculifer provides a seasonal SST estimate from the upper ~50 m of the water column generated during upwelling in austral low-latitude fall/winter. Common to both down-core records is a glacial-interglacial amplitude of ~3°-3.5°C for the last climatic changes and lower Holocene and glacial oxygen isotope stage 2 temperatures compared with interglacial stage 5.5 and glacial stage 6 temperatures, respectively. The comparison to published SST estimates from alkenones, oxygen isotopes, and foraminiferal transfer function from the same core material pinpoints discrepancies and conformities between methods.

Importance of Nuclear Data Uncertainties in Criticality Calculations

To study the propagation of the uncertainty from basic data across different scale and physics phenomena -> through complex coupled multi-physics and multi-scale simulations