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

Electron localization function at the correlated level

The electron localization function (ELF) has been proven so far a valuable tool to determine the location of electron pairs. Because of that, the ELF has been widely used to understand the nature of the chemical bonding and to discuss the mechanism of chemical reactions. Up to now, most applications of the ELF have been performed with monodeterminantal methods and only few attempts to calculate this function for correlated wave functions have been carried out. Here, a formulation of ELF valid for mono- and multiconfigurational wave functions is given and compared with previous recently reported approaches. The method described does not require the use of the homogeneous electron gas to define the ELF, at variance with the ELF definition given by Becke. The effect of the electron correlation in the ELF, introduced by means of configuration interaction with singles and doubles calculations, is discussed in the light of the results derived from a set of atomic and molecular systems

A chemical Hamiltonian approach study of the basis set superposition error changes on electron densities and one- and two-center energy components

A comparision of the local effects of the basis set superposition error (BSSE) on the electron densities and energy components of three representative H-bonded complexes was carried out. The electron densities were obtained with Hartee-Fock and density functional theory versions of the chemical Hamiltonian approach (CHA) methodology. It was shown that the effects of the BSSE were common for all complexes studied. The electron density difference maps and the chemical energy component analysis (CECA) analysis confirmed that the local effects of the BSSE were different when diffuse functions were present in the calculations

Effect of basis set superposition error on the electron density of molecular complexes

The effect of basis set superposition error (BSSE) on molecular complexes is analyzed. The BSSE causes artificial delocalizations which modify the first order electron density. The mechanism of this effect is assessed for the hydrogen fluoride dimer with several basis sets. The BSSE-corrected first-order electron density is obtained using the chemical Hamiltonian approach versions of the Roothaan and Kohn-Sham equations. The corrected densities are compared to uncorrected densities based on the charge density critical points. Contour difference maps between BSSE-corrected and uncorrected densities on the molecular plane are also plotted to gain insight into the effects of BSSE correction on the electron density

A quantum molecular similarity analysis of changes in molecular electron density caused by basis set flotation and electric field application

Quantum molecular similarity (QMS) techniques are used to assess the response of the electron density of various small molecules to application of a static, uniform electric field. Likewise, QMS is used to analyze the changes in electron density generated by the process of floating a basis set. The results obtained show an interrelation between the floating process, the optimum geometry, and the presence of an external field. Cases involving the Le Chatelier principle are discussed, and an insight on the changes of bond critical point properties, self-similarity values and density differences is performed

Basis set and electron correlation effects on initial convergence for vibrational nonlinear optical properties of conjugated organic molecules

The level of ab initio theory which is necessary to compute reliable values for the static and dynamic (hyper)polarizabilities of three medium size π-conjugated organic nonlinear optical (NLO) molecules is investigated. With the employment of field-induced coordinates in combination with a finite field procedure, the calculations were made possible. It is stated that to obtain reasonable values for the various individual contributions to the (hyper)polarizability, it is necessary to include electron correlation. Based on the results, the convergence of the usual perturbation treatment for vibrational anharmonicity was examined

Estimates of electronic coupling for excess electron transfer in DNA

Electronic coupling Vda is one of the key parameters that determine the rate of charge transfer through DNA. While there have been several computational studies of Vda for hole transfer, estimates of electronic couplings for excess electron transfer (ET) in DNA remain unavailable. In the paper, an efficient strategy is established for calculating the ET matrix elements between base pairs in a π stack. Two approaches are considered. First, we employ the diabatic-state (DS) method in which donor and acceptor are represented with radical anions of the canonical base pairs adenine-thymine (AT) and guanine-cytosine (GC). In this approach, similar values of Vda are obtained with the standard 6-31 G* and extended 6-31++ G* basis sets. Second, the electronic couplings are derived from lowest unoccupied molecular orbitals (LUMOs) of neutral systems by using the generalized Mulliken-Hush or fragment charge methods. Because the radical-anion states of AT and GC are well reproduced by LUMOs of the neutral base pairs calculated without diffuse functions, the estimated values of Vda are in good agreement with the couplings obtained for radical-anion states using the DS method. However, when the calculation of a neutral stack is carried out with diffuse functions, LUMOs of the system exhibit the dipole-bound character and cannot be used for estimating electronic couplings. Our calculations suggest that the ET matrix elements Vda for models containing intrastrand thymine and cytosine bases are essentially larger than the couplings in complexes with interstrand pyrimidine bases. The matrix elements for excess electron transfer are found to be considerably smaller than the corresponding values for hole transfer and to be very responsive to structural changes in a DNA stack

Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps

Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps Darbellay Bastien Institut de Minéralogie et Géochimie Résumé grand public Comme toute matière, la roche est sensible à son environnement et cherche à s'adapter pour acquérir un état stable (état d'équilibre). Les changements des conditions physiques (température et pression) vont ainsi impliquer des modifications dans la roche. Le métamorphisme est l'étude de ces changements. Les minéraux qui constituent la roche peuvent modifier, leur structure, leur chimie ou être remplacer par d'autres minéraux plus stables. Il est ainsi crucial de déterminer les processus responsables et limitant de la croissance minérale. Trois processus permettent la croissance ; (1) la dissolution des éléments du réactant, (2) le transport de ces éléments vers le site de croissance, (3) l'incorporation de ces éléments dans la nouvelle structure. Cette thèse se focalise sur les structures des minéraux de haute pression (forme, zonation chimique, structure interne) pour essayer de déterminer les facteurs importants à l'origine de leur état final. Les zones d'étude se situent dans la zone de Sésia. La première partie traite de la problématique liée à l'incorporation d'un élément dans une structure minérale. A l'image de la croissance humaine, les irrégularités minéralogiques permettent de mettre en lumière un dysfonctionnement de la croissance due à un excès ou à une carence d'un élément. Bien dosé, cet élément est cependant essentiel à la croissance. Les zoisites (épidotes) des métabasites de la région de Cima di Bonze montrent une zonation chimique en sablier. Dans cette zonation la teneur en fer excède la capacité maximum que peut contenir la structure orthorhombique de la zoisite. Des défauts de structure permettent l'accommodation de cet excès. La zoisite peut ainsi adapter sa structure pour permettre l'incorporation d'une relativement grande quantité de fer. Les études précédentes montraient, pour des conditions similaires, la formation de deux épidotes distinctes. La deuxième partie se penche sur la compétition entre le minéral qui fait sa croissance et les minéraux (réactants) qui l'entourent. Les métapélites de la région du Monte Mucrone contiennent des grenats atollaires. Des études détaillées de la texture et de la zonation chimique du grenat ainsi qu'une modélisation thermodynamique ont permis de mieux cerner les facteurs importants responsables de la forme atollaire. Cette structure est obtenue par un changement du comportement de la croissance du grenat le long d'un chemin P-T hercynien. Dans un premier stade, le grenat croît rapidement et consume peu le quartz de la matrice. La croissance se fait ainsi le long des jointures des grains de quartz. Dans un second temps, les changements de conditions PT donnent une croissance lente du grenat et une forte consommation du quartz. Le grenat peut ainsi développer sa forme dodécaédrale classique. La troisième partie s'intéresse aux distances de transport par diffusion d'un élément (ici l'argon) durant la haute pression. Pour ce faire, un profile d'âges 40Ar/39Ar sur biotite a été mesuré depuis un veine de haute pression riche en argon jusque dans son encaissant (granitoïd du Monte Mucrone). Le profile montre une répartition des âges suivant une courbe de diffusion. Le transport se fait sur une longueur de deux centimètre avec l'aide d'un fluide. Il est réduit à une échelle millimétrique quand la phase fluide disparaît. Cette étude montre ainsi les difficultés de transport des éléments durant la haute pression ne permettant pas un rééquilibrage de la roche à grande échelle. Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps Darbellay Bastien Institut de Minéralogie et Géochimie Résumé de thèse Les processus de croissance (diffusion des éléments et les réactions d'interface) et les conditions dans lesquelles les minéraux grandissent (température, pression, fluide, composition chimique de la roche), déterminent la texture ainsi que la zonation des minéraux. Cette thèse se focalise, par le biais de textures peu communes, sur trois différents processus impliqués dans la croissance minérale à haute pression (Zone de Sésia, Alpes de l'Ouest, Italie). L'incorporation d'un élément dans une structure minérale ne peut se faire que dans des sites en accord avec la taille et la charge ionique de l'élément. De plus, la balance de charge doit être maintenue dans le minéral. La régularité de la structure cristalline fixe ainsi une limite maximum de concentration d'un élément donné. Les zoisites provenant des métabasites de la région de Cima di Bonze montrent des zonations en sablier caractérisées par une concentration anormale en fer. La zonation se marque par une différente teinte de biréfringence et par un plus grand angle d'extinction que le reste de la zoisite. Une inter-croissance de clinozoisite à l'intérieur de la structure orthorhombique de la zoisite peut ainsi être suspectée. Les analyses XRD (diffraction des rayons x) ainsi que les analyses Raman ne confirment pas cette suspicion. Seules les analyses TEM (microscope à électrons transmis) montrent des défauts de structure pouvant être interprétés comme des modules de clinozoisite. Ils ne peuvent cependant pas être considérés comme une phase thermodynamique. Un nouveau trou d'immiscibilité entre deux zoisite (X ep= 0.1 and Xep = 0.15) a ainsi pu être établi. Dans les métapélites la région du Monte Mucrone, des grenats fortement zonés montrent une évolution texturale singulière. Ils présentent une forme initiale de `champignon' qui se développe pour former une structure atollaire finale. L'étude conjuguée de la structure 3D et des zonations, ainsi que l'établissement d'un model thermodynamique, indiquent que ces structures proviennent de deux épisodes de croissances : (1) La croissance du grenat durant un chemin prograde hercynien (de 525 °C et 6.2 kbar à 640 °C et 9 kbar) permet la formation des textures atollaires. Elles sont le résultat d'une croissance poecilitique initiale suivie d'une croissance idiomorphique du grenat. (2) La structure est rendue plus complexe par la cristallisation d'un grenat homogène tout autour ainsi qu'à l'intérieur du grenat hercynien durant la haute pression alpine (550 °C and 20 kbar). L'arrivée de l'eau durant la haute pression facilite le transport d'éléments et permet une cristallisation rapide du grenat. La diffusion peut être un facteur limitant de la croissance minéralogique. Elle a aussi une grande importance pour la géochronologie. Une veine de haute pression à l'intérieur du granitoïde du Monte Mucrone a été étudiée dans le but de déterminer la distance de diffusion de l'argon. Le profile d'âges 40Ar/39Ar sur biotites, établi de la veine vers le métagranitoïde, suit une courbe de diffusion. Les âges sont élevés proche de la veine (800 Ma) puis décroissent jusqu'à des âges homogènes (170-150 Ma) à deux centimètres de la veine. La présence de fluide, marqué par de hautes concentrations en chlore, permet une diffusion centimétrique. Cependant, la distance est réduite à une échelle millimétrique quand le fluide est absent. Les très faibles distances de diffusion préservent les âges pré-alpins et impliquent un événement géologique pour les âges de 170-150 Ma. Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps Darbellay Bastien Institut de Minéralogie et Géochimie Thesis abstract Rock textures and zonings are the consequence of growth processes (element diffusion and interface reaction) steered by the environment in which they grew (pressure, differential stress, temperature, fluid and rock composition). The thesis presented here focuses on three different topics, each of it dealing with aspects of mineral growth processes during subduction, in a high-pressure environment. All studies were conducted in the Sesia Zone of the Western European Alps, Italy. The first study addresses the crystallography and geochemistry of element incorporation in zoisite, one of the major hydrous minerals found in subduction zone rocks. Elements can be incorporated into a mineral structure only on crystallographic sites that offer enough space for the ion and the overall charge balance has to be maintained. Element concentrations are hence limited. Incorporation of some elements produces complex zoning, including hourglass like patterns, which are the focus of the first contribution. Zoisites from Cima di Bonze (Sesia Zone) show spectacular hourglass zoning defined by Fe-content variations. The hourglass zones have a distinct birefringence and a different extinction angle than the regular part of the zoisite. We show by detailed XRD (X-ray diffraction) and confocal Raman analyses that the high Fe-zones are nevertheless zoisite, and not clinozoisite as one might expect. High resolution TEM (transmission electron microscopy) analyses show planar defects on (100) that can be interpreted as small-scale clinozoisite modules. However, these clinozoisites cannot be interpreted as a distinctive thermodynamic phase and the entire mineral has to be considered as zoisite. The miscibility gap between two zoisites (Xep = 0.1 and Xep = 0.15) can be then definite at 550 ± 50°C and 14 to 20 Kbar. Strongly zoned garnets in quartz rich metapelite from the Monte Mucrone area (Sesia Zone) show evolution form 3D mushroom to atoll structure. The second contribution presents textural investigations, garnet zoning and thermodynamic modeling that demonstrate that atoll garnets are the result of two distinctive growth events. (1) Garnet atoll structure is already formed during a prograde Hercynian path from 525 °C and 6.2 kbar to 640 °C and 9 kbar. It results in an initial poikilitic growth followed by a final idiomorphic growth event. (2) Alpine HP garnet are homogenous (550 °C and 20 kbar) and grew around and also inside the Hercynian garnet. Lack of prograde Alpine garnet and fast growth of the HP garnet is explained by the absence of water during much of the prograde path. Water saturation was only observed towards the end, close towards the peak metamorphic conditions. Diffusion could be a limiting factor for crystal growth. It has also a great importance in geochronology. HP vein inside the metagranitoide of the Monte Mucrone (~300 Ma) was investigated to determine argon diffusion scales during high-pressure metamorphism. 40Ar/39Ar biotite ages profile from the vein toward the metagranodiorite show a diffusion curve: old ages (800 Ma) located close to the vein decrease until homogenous 170-150 Ma ages are obtained, two centimeter away from the vein. Centimeter-scale diffusion occurs with help of a fluid phase marked by high chlorine concentrations. Argon diffusion is reduced to a millimeter scale when free fluid is absent. Very short diffusion distance permits to preserve pre-Alpine ages. The 170-150 Ma ages are considered to be geologic meaningful, probably resulting from the extensional tectonics linked to opening of the Tethian ocean.

Correlative light and electron microscopy in parasite research.

The interaction of a parasite and a host cell is a complex process, which involves several steps: (1) attachment to the plasma membrane, (2) entry inside the host cell, and (3) hijacking of the metabolism of the host. In biochemical experiments, only an event averaged over the whole cell population can be analyzed. The power of microscopy, however, is to investigate individual events in individual cells. Therefore, parasitologists frequently perform experiments with fluorescence microscopy using different dyes to label structures of the parasite or the host cell. Though the resolution of light microscopy has greatly improved, it is not sufficient to reveal interactions at the ultrastructural level. Furthermore, only specifically labeled structures can be seen and related to each other. Here, we want to demonstrate the additional value of electron microscopy in this area of research. Investigation of the different steps of parasite-host cell interaction by electron microscopy, however, is often hampered by the fact that there are only a few cells infected, and therefore it is difficult to find enough cells to study. A solution is to profit from low magnification, hence large overview, and specific location of the players by fluorescence labels in a light microscope with the high power resolution and structural information provided by an electron microscope, in short by correlative light and electron microscopy.

Analysis of acute brain slices by electron microscopy: A correlative light-electron microscopy workflow based on Tokuyasu cryo-sectioning.

Acute brain slices are slices of brain tissue that are kept vital in vitro for further recordings and analyses. This tool is of major importance in neurobiology and allows the study of brain cells such as microglia, astrocytes, neurons and their inter/intracellular communications via ion channels or transporters. In combination with light/fluorescence microscopies, acute brain slices enable the ex vivo analysis of specific cells or groups of cells inside the slice, e.g. astrocytes. To bridge ex vivo knowledge of a cell with its ultrastructure, we developed a correlative microscopy approach for acute brain slices. The workflow begins with sampling of the tissue and precise trimming of a region of interest, which contains GFP-tagged astrocytes that can be visualised by fluorescence microscopy of ultrathin sections. The astrocytes and their surroundings are then analysed by high resolution scanning transmission electron microscopy (STEM). An important aspect of this workflow is the modification of a commercial cryo-ultramicrotome to observe the fluorescent GFP signal during the trimming process. It ensured that sections contained at least one GFP astrocyte. After cryo-sectioning, a map of the GFP-expressing astrocytes is established and transferred to correlation software installed on a focused ion beam scanning electron microscope equipped with a STEM detector. Next, the areas displaying fluorescence are selected for high resolution STEM imaging. An overview area (e.g. a whole mesh of the grid) is imaged with an automated tiling and stitching process. In the final stitched image, the local organisation of the brain tissue can be surveyed or areas of interest can be magnified to observe fine details, e.g. vesicles or gold labels on specific proteins. The robustness of this workflow is contingent on the quality of sample preparation, based on Tokuyasu's protocol. This method results in a reasonable compromise between preservation of morphology and maintenance of antigenicity. Finally, an important feature of this approach is that the fluorescence of the GFP signal is preserved throughout the entire preparation process until the last step before electron microscopy.

Electron microscopy of RecA-DNA complexes: Two different states, their functional significance and relation to the solved crystal structure

Seven different electron microscopy techniques habe been employed to study the RecA protein of E. coli. This review provides a summary of the conclusions that have been drawn from these studies, and attempts to relate these observations to models for the role of RecA protein in homologous recombination.

Synapses on motoneuron dendrites in the brachial section of the frog spinal cord: a computer-aided electron microscopic study of cobalt-filled cells.

Cobalt-labelled motoneuron dendrites of the frog spinal cord at the level of the second spinal nerve were photographed in the electron microscope from long series of ultrathin sections. Three-dimensional computer reconstructions of 120 dendrite segments were analysed. The samples were taken from two locations: proximal to cell body and distal, as defined in a transverse plane of the spinal cord. The dendrites showed highly irregular outlines with many 1-2 microns-long 'thorns' (on average 8.5 thorns per 100 microns 2 of dendritic area). Taken together, the reconstructed dendrite segments from the proximal sites had a total length of about 250 microns; those from the distal locations, 180 microns. On all segments together there were 699 synapses. Nine percent of the synapses were on thorns, and many more close to their base on the dendritic shaft. The synapses were classified in four groups. One third of the synapses were asymmetric with spherical vesicles; one half were symmetric with spherical vesicles; and one tenth were symmetric with flattened vesicles. A fourth, small class of asymmetric synapses had dense-core vesicles. The area of the active zones was large for the asymmetric synapses (median value 0.20 microns 2), and small for the symmetric ones (median value 0.10 microns 2), and the difference was significant. On average, the areas of the active zones of the synapses on thin dendrites were larger than those of synapses on large calibre dendrites. About every 4 microns 2 of dendritic area received one contact. There was a significant difference between the areas of the active zones of the synapses at the two locations. Moreover, the number per unit dendritic length was correlated with dendrite calibre. On average, the active zones covered more than 4% of the dendritic area; this value for thin dendrites was about twice as large as that of large calibre dendrites. We suggest that the larger active zones and the larger synaptic coverage of the thin dendrites compensate for the longer electrotonic distance of these synapses from the soma.