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.

Femtosecond laser-based investigations of magnetic circular dichroism in near-threshold photoemission

During the last decades magnetic circular dichroism (MCD) has attracted much interest and evolved into various experimental methods for the investigation of magnetic thin films. For example, synchrotron-based X-ray magnetic circular dichroism (XMCD) displays the absolute values of spin and orbital magnetic moments. It thereby benefits from large asymmetry values of more than 30% due to the excitation of atomic core-levels. Similarly large values are also expected for threshold photoemission magnetic circular dichroism (TPMCD). Using lasers with photon energies in the range of the sample work function this method gives access to the occupied electronic structure close to the Fermi level. However, except for the case of Ni(001) there exist only few studies on TPMCD moreover revealing much smaller asymmetries than XMCD-measurements. Also the basic physical mechanisms of TPMCD are not satisfactorily understood. In this work we therefore investigate TPMCD in one- and two-photon photoemission (1PPE and 2PPE) for ferromagnetic Heusler alloys and ultrathin Co films using ultrashort pulsed laser light. The observed dichroism is explained by a non-conventional photoemission model using spin-resolved band-structure calculations and linear response theory. For the two Heusler alloys Ni2MnGa and Co2FeSi we give first evidence of TPMCD in the regime of two-photon photoemission. Systematic investigations concerning general properties of TPMCD in 1PPE and 2PPE are carried out at ultrathin Co films grown on Pt(111). Here, photon-energy dependent measurements reveal asymmetries of 1.9% in 1PPE and 11.7% in 2PPE. TPMCD measurements at decreased work function even yield larger asymmetries of 6.2% (1PPE) and 17% (2PPE), respectively. This demonstrates that enlarged asymmetries are also attainable for the TPMCD effect on Co(111). Furthermore, we find that the TPMCD asymmetry is bulk-sensitive for 1PPE and 2PPE. This means that the basic mechanism leading to the observed dichroism must be connected to Co bulk properties; surface effects do not play a crucial role. Finally, the enhanced TPMCD asymmetries in 2PPE compared to the 1PPE case are traced back to the dominant influence of the first excitation step and the existence of a real intermediate state. The observed TPMCD asymmetries cannot be interpreted by conventional photoemission theory which only considers direct interband transitions in the direction of observation (Γ-L). For Co(111), these transitions lead to evanescent final states. The excitation to such states, however, is incompatible with the measured bulk-sensitivity of the asymmetry. Therefore, we generalize this model by proposing the TPMCD signal to arise mostly from direct interband transitions in crystallographic directions other than (Γ-L). The necessary additional momentum transfer to the excited electrons is most probably provided by electron-phonon or -magnon scattering processes. Corresponding calculations on the basis of this model are in reasonable agreement with the experimental results so that this approach represents a promising tool for a quantitative description of the TPMCD effect. The present findings encourage an implementation of our experimental technique to time- and spatially-resolved photoemission electron microscopy, thereby enabling a real time imaging of magnetization dynamics of single excited states in a ferromagnetic material on a femtosecond timescale.