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.

Novel algorithms for electronic structure based molecular dynamics with linear system-size scaling

Die vorliegende Arbeit behandelt die Entwicklung und Verbesserung von linear skalierenden Algorithmen für Elektronenstruktur basierte Molekulardynamik. Molekulardynamik ist eine Methode zur Computersimulation des komplexen Zusammenspiels zwischen Atomen und Molekülen bei endlicher Temperatur. Ein entscheidender Vorteil dieser Methode ist ihre hohe Genauigkeit und Vorhersagekraft. Allerdings verhindert der Rechenaufwand, welcher grundsätzlich kubisch mit der Anzahl der Atome skaliert, die Anwendung auf große Systeme und lange Zeitskalen. Ausgehend von einem neuen Formalismus, basierend auf dem großkanonischen Potential und einer Faktorisierung der Dichtematrix, wird die Diagonalisierung der entsprechenden Hamiltonmatrix vermieden. Dieser nutzt aus, dass die Hamilton- und die Dichtematrix aufgrund von Lokalisierung dünn besetzt sind. Das reduziert den Rechenaufwand so, dass er linear mit der Systemgröße skaliert. Um seine Effizienz zu demonstrieren, wird der daraus entstehende Algorithmus auf ein System mit flüssigem Methan angewandt, das extremem Druck (etwa 100 GPa) und extremer Temperatur (2000 - 8000 K) ausgesetzt ist. In der Simulation dissoziiert Methan bei Temperaturen oberhalb von 4000 K. Die Bildung von sp²-gebundenem polymerischen Kohlenstoff wird beobachtet. Die Simulationen liefern keinen Hinweis auf die Entstehung von Diamant und wirken sich daher auf die bisherigen Planetenmodelle von Neptun und Uranus aus. Da das Umgehen der Diagonalisierung der Hamiltonmatrix die Inversion von Matrizen mit sich bringt, wird zusätzlich das Problem behandelt, eine (inverse) p-te Wurzel einer gegebenen Matrix zu berechnen. Dies resultiert in einer neuen Formel für symmetrisch positiv definite Matrizen. Sie verallgemeinert die Newton-Schulz Iteration, Altmans Formel für beschränkte und nicht singuläre Operatoren und Newtons Methode zur Berechnung von Nullstellen von Funktionen. Der Nachweis wird erbracht, dass die Konvergenzordnung immer mindestens quadratisch ist und adaptives Anpassen eines Parameters q in allen Fällen zu besseren Ergebnissen führt.

Symmetry in chemistry from the hydrogen atom to proteins

The last 2 decades have seen discoveries in highly excited states of atoms and molecules of phenomena that are qualitatively different from the “planetary” model of the atom, and the near-rigid model of molecules, characteristic of these systems in their low-energy states. A unified view is emerging in terms of approximate dynamical symmetry principles. Highly excited states of two-electron atoms display “molecular” behavior of a nonrigid linear structure undergoing collective rotation and vibration. Highly excited states of molecules described in the “standard molecular model” display normal mode couplings, which induce bifurcations on the route to molecular chaos. New approaches such as rigid–nonrigid correlation, vibrons, and quantum groups suggest a unified view of collective electronic motion in atoms and nuclear motion in molecules.

Imaging of single molecule diffusion.

In recent years observations at the level of individual atoms and molecules became possible by microscopy and spectroscopy. Imaging of single fluorescence molecules has been achieved but has so far been restricted to molecules in the immobile state. Here we provide methodology for visualization of the motion of individual fluorescent molecules. It is applied to imaging of the diffusional path of single molecules in a phospholipid membrane by using phospholipids carrying one rhodamine dye molecule. For this methodology, fluorescence microscopy was carried to a sensitivity so that single fluorescent molecules illuminated for only 5 ms were resolvable at a signal/noise ratio of 28. Repeated illuminations permitted direct observation of the diffusional motion of individual molecules with a positional accuracy of 30 nm. Such capability has fascinating potentials in bioscience--for example, to correlate biological functions of cell membranes with movements, spatial organization, and stoichiometries of individual components.

The OPIT system part I.

The OPIT program is briefly described. OPIT is a basis-set-optimising, self-consistent field, molecular orbital program for calculating properties of closed-shell ground states of atoms and molecules. A file handling technique is then put forward which enables core storage to be used efficiently in large FORTRAN scientific applications programs. Hashing and list processing techniques, of the type frequently used in writing system software and computer operating systems, are here applied to the creation of data files (integral label and value lists etc.). Files consist of a chained series of blocks which may exist in core or on backing store or both. Efficient use of core store is achieved and the processes of file deletion, file re-writing and garbage collection of unused blocks can be easily arranged. The scheme is exemplified with reference to the OPIT program. A subsequent paper will describe a job scheduling scheme for large programs of this sort.

The relevance of the Laplacian of intracule and extracule density distributions for analyzing electron-electron interactions in molecules

A topological analysis of intracule and extracule densities and their Laplacians computed within the Hartree-Fock approximation is presented. The analysis of the density distributions reveals that among all possible electron-electron interactions in atoms and between atoms in molecules only very few are located rigorously as local maxima. In contrast, they are clearly identified as local minima in the topology of Laplacian maps. The conceptually different interpretation of intracule and extracule maps is also discussed in detail. An application example to the C2H2, C2H4, and C2H6 series of molecules is presented

The relevance of the Laplacian of intracule and extracule density distributions for analyzing electron-electron interactions in molecules

A topological analysis of intracule and extracule densities and their Laplacians computed within the Hartree-Fock approximation is presented. The analysis of the density distributions reveals that among all possible electron-electron interactions in atoms and between atoms in molecules only very few are located rigorously as local maxima. In contrast, they are clearly identified as local minima in the topology of Laplacian maps. The conceptually different interpretation of intracule and extracule maps is also discussed in detail. An application example to the C2H2, C2H4, and C2H6 series of molecules is presented

Ab initio study of SO2 molecules interacting with pristine and transition metal covered fullerenes as a possible route for nanofilters

In this work, a systematic study of SO2 molecules interacting with pristine and transition metal (TM) covered C-60 is presented by means of first principles calculations. It is observed that the SO2 molecule interacts weakly with the pristine C-60 fullerene, although the resulting interaction is largely increased when the C-60 structure is covered with Fe, Mn, or Ti atoms and the SO2 Molecules are bounded through the TM atoms. The number of bounded SO2 molecules per TM atoms, in addition to the elevated binding energies per molecules, allows us to conclude that such composites can be used as a template for efficient devices to remove SO2 molecules or, alternatively, as SO2 gas sensor.

Superchemistry: dynamics of coupled atomic and molecular Bose condensates

We analyze the dynamics of a dilute, trapped Bose-condensed atomic gas coupled to a diatomic molecular Bose gas by coherent Raman transitions. This system is shown to result in a new type of “superchemistry,” in which giant collective oscillations between the atomic and the molecular gas can occur. The phenomenon is caused by stimulated emission of bosonic atoms or molecules into their condensate phases.

Using gold nanoparticles in protein crystallography: studies in crystal growth and derivatization

Dissertation for the Master Degree in Structural and Functional Biochemistry

Optical-model potential for electron and positron elastic scattering by atoms

An optical-model potential for systematic calculations of elastic scattering of electrons and positrons by atoms and positive ions is proposed. The electrostatic interaction is determined from the Dirac-Hartree-Fock self-consistent atomic electron density. In the case of electron projectiles, the exchange interaction is described by means of the local-approximation of Furness and McCarthy. The correlation-polarization potential is obtained by combining the correlation potential derived from the local density approximation with a long-range polarization interaction, which is represented by means of a Buckingham potential with an empirical energy-dependent cutoff parameter. The absorption potential is obtained from the local-density approximation, using the Born-Ochkur approximation and the Lindhard dielectric function to describe the binary collisions with a free-electron gas. The strength of the absorption potential is adjusted by means of an empirical parameter, which has been determined by fitting available absolute elastic differential cross-section data for noble gases and mercury. The Dirac partial-wave analysis with this optical-model potential provides a realistic description of elastic scattering of electrons and positrons with energies in the range from ~100 eV up to ~5 keV. At higher energies, correlation-polarization and absorption corrections are small and the usual static-exchange approximation is sufficiently accurate for most practical purposes.

Evaluation of group electronegativities and hardness (softness) of group 14 elements and containing functional groups through density functional theory and correlation with NMR spectra data

Quantum Chemical calculations for group 14 elements of Periodic Table (C, Si, Ge, Sn, Pb) and their functional groups have been carried out using Density Functional Theory (DFT) based reactivity descriptors such as group electronegativities, hardness and softness. DFT calculations were performed for a large series of tetracoordinated Sn compounds of the CH3SnRR'X type, where X is a halogen and R and R' are alkyl, halogenated alkyl, alkoxy, or alkyl thio groups. The results were interpreted in terms of calculated electronegativity and hardness of the SnRR'X groups, applying a methodology previously developed by Geerlings and coworkers (J. Phys. Chem. 1993, 97, 1826). These calculations allowed to see the regularities concerning the influence of the nature of organic groups RR' and inorganic group X on electronegativities and hardness of the SnRR'X groups; in this case, it was found a very good correlation between the electronegativity of the fragment and experimental 119Sn chemical shifts, a property that sensitively reflects the change in the valence electronic structure of molecules. This work was complemented with the study of some compounds of the EX and ER types, where E= C, Si, Ge, Sn and R= CH3, H, which was performed to study the influence that the central atom has on the electronegativity and hardness of molecules, or whether these properties are mainly affected for the type of ligand bound to the central atom. All these calculations were performed using the B3PW91 functional together with the 6-311++G** basis set level for H, C, Si, Ge, F, Cl and Br atoms and the 3-21G for Sn and I atoms.

Graphene from Graphite by Chemical and Physical Techniques.

Graphene is a material with extraordinary properties. Its mechanical and electrical properties are unparalleled but the difficulties in its production are hindering its breakthrough in on applications. Graphene is a two-dimensional material made entirely of carbon atoms and it is only a single atom thick. In this work, properties of graphene and graphene based materials are described, together with their common preparation techniques and related challenges. This Thesis concentrates on the topdown techniques, in which natural graphite is used as a precursor for the graphene production. Graphite consists of graphene sheets, which are stacked together tightly. In the top-down techniques various physical or chemical routes are used to overcome the forces keeping the graphene sheets together, and many of them are described in the Thesis. The most common chemical method is the oxidisation of graphite with strong oxidants, which creates a water-soluble graphene oxide. The properties of graphene oxide differ significantly from pristine graphene and, therefore, graphene oxide is often reduced to form materials collectively known as reduced graphene oxide. In the experimental part, the main focus is on the chemical and electrochemical reduction of graphene oxide. A novel chemical route using vanadium is introduced and compared to other common chemical graphene oxide reduction methods. A strong emphasis is placed on electrochemical reduction of graphene oxide in various solvents. Raman and infrared spectroscopy are both used in in situ spectroelectrochemistry to closely monitor the spectral changes during the reduction process. These in situ techniques allow the precise control over the reduction process and even small changes in the material can be detected. Graphene and few layer graphene were also prepared using a physical force to separate these materials from graphite. Special adsorbate molecules in aqueous solutions, together with sonic treatment, produce stable dispersions of graphene and few layer graphene sheets in water. This mechanical exfoliation method damages the graphene sheets considerable less than the chemical methods, although it suffers from a lower yield.

The design, synthesis and characterization of new building blocks for the preparation of molecule-based magnetic materials /

Two new families of building blocks have been prepared and fully characterized and their coordination chemistry exploited for the preparation of molecule-based magnetic materials. The first class of compounds were prepared by exploiting the chemistry of 3,3'-diamino-2,2'-bipyridine together with 2-pyridine carbonyl chloride or 2-pyridine aldehyde. Two new ligands, 2,2'-bipyridine-3,3'-[2-pyridinecarboxamide] (Li, 2.3) and N'-6/s(2-pyridylmethyl) [2,2'bipyridine]-3,3'-diimine (L2, 2.7), were prepared and characterized. For ligand L4, two copper(II) coordination compounds were isolated with stoichiometrics [Cu2(Li)(hfac)2] (2.4) and [Cu(Li)Cl2] (2.5). The molecular structures of both complexes were determined by X-ray crystallography. In both complexes the ligand is in the dianionic form and coordinates the divalent Cu(II) ions via one amido and two pyridine nitrogen donor atoms. In (2.4), the coordination geometry around both Cu11 ions is best described as distorted trigonal bipyramidal where the remaining two coordination sites are satisfied by hfac counterions. In (2.5), both Cu(II) ions adopt a (4+1) distorted square pyramidal geometry. One copper forms a longer apical bond to an adjacent carbonyl oxygen atom, whereas the second copper is chelated to a neighboring Cu-Cl chloride ion to afford chloride bridged linear [Cu2(Li)Cl2]2 tetramers that run along the c-axis of the unit cell. The magnetic susceptibility data for (2.4) reveal the occurrence of weak antiferromagnetic interactions between the copper(II) ions. In contrast, variable temperature magnetic susceptibility measurements for (2.5) reveal more complex magnetic properties with the presence of ferromagnetic exchange between the central dimeric pair of copper atoms and weak antiferromagnetic exchange between the outer pairs of copper atoms. The Schiff-base bis-imine ligand (L2, 2.7) was found to be highly reactive; single crystals grown from dry methanol afforded compound (2.14) for which two methanol molecules had added across the imine double bond. The susceptibility of this ligand to nucleophilic attack at its imine functionality assisted via chelation to Lewis acidic metal ions adds an interesting dimension to its coordination chemistry. In this respect, a Co(II) quaterpyridine-type complex was prepared via a one-pot transformation of ligand L2 in the presence of a Lewis acidic metal salt. The rearranged complex was characterized by X-ray crystallography and a reaction mechanism for its formation has been proposed. Three additional rearranged complexes (2.13), (2.17) and (2.19) were also isolated when ligand (L2, 2.7) was reacted with transition metal ions. The molecular structures of all three complexes have been determined by X-ray crystallography. The second class of compounds that are reported in this thesis, are the two diacetyl pyridine derivatives, 4-pyridyl-2,6-diacetylpyridine (5.5) and 2,2'-6,6'-tetraacetyl-4,4'-bipyridine (5.15). Both of these compounds have been designed as intermediates for the metal templated assembly of a Schiff-base N3O2 macrocycle. From compound (5.15), a covalently tethered dimeric Mn(II) macrocyclic compound of general formula {[Mn^C^XJCl-FkO^Cl-lO.SFbO (5.16) was prepared and characterized. The X-ray analysis of (5.16) reveals that the two manganese ions assume a pentagonal-bipyramidal geometry with the macrocycle occupying the pentagonal plane and the axial positions being filled by a halide ion and a H2O molecule. Magnetic susceptibility data reveal the occurrence of antiferromagnetic interactions between covalently tethered Mn(II)-Mn(II) dimeric units. Following this methodology a Co(II) analogue (5.17) has also been prepared which is isostructural with (5.16).

Investigation of molecular polarizabilities and derivatives in halomethanes

The research undertaken was to obtain absolute Raman intensities for the symmetric stretching vibrations of the methyl halides, CH3X with (X=F, CI, Br), by experiment and theory. The intensities were experimentally measured using the Ar+ ion gas laser as excitation source, a Spex 14018 double monochromator and a RCA C-31034 photomultiplier tube as detector. These intensities arise from changes in the derivative of the polarizability (8 a'), with respect to vibration along a normal coordinate (8qi). It was intended that these derivatives obtained with respect to normal coordinates would be converted to derivatives with respect to internal coordinates, for a quantitative comparison with theory. Theoretical numerical polarizability derivatives for the stretching vibrations are obtained using the following procedure. A vibration was simulated in the molecule by increasi.ng and decreasing the respective bond by the amount ±o.oosA for the C-H bonds and ±o.oIA for the C-X (X=F, CI, Br) bond. The derivative was obtained by taking the difference in the polarizability for the equilibrium geometry and the geometry when a particular bond is changed. This difference, when divided by the amount of change in each bond and the number of bonds present results in the derivative of the polarizability with respect to internal coordinate i.e., !1u/!1r. These derivatives were obtained by two methods: I} ab initio molecular orbital calculation and 2} theory of atoms in molecules (AIM) analysis. Due to errors in the experimental setup only a qualitative analysis of the results was undertaken relative to the theory. Theoretically it is predicted that the symmetric carbonhalogen stretch vibrations are more intense than the respective carbon-hydrogen stretch, but only for the methyl chloride and bromide. The carbon fluorine stretch is less intense than the carbon-hydrogen stretch, a fact which is attributed to the small size and high electronegativity of the fluorine atom. The experimental observations are seen to agree qualitatively with the theory results. It is hoped that when the experiment is repeated, a quantitative comparison can be made. The analysis by the theory of atoms in molecules, along with providing polarizabilities and polarizability derivatives, gives additional information outlined below. The theory provides a pictorial description of the main factors contributing to the molecular polarizability and polarizability derivative. These contributions are from the charge transfer and atomic dipole terms i.e., transfer of charge from one atom to another and the reorganization of atomic electronic charge distribution due to presence of an electric field. The linear relationship between polarizability and molecular volume was also observed.