Development, implementation and application of electronic structural descriptors to the analysis of the chemical bonding, aromaticity and chemical reactivity

En la literatura sobre mecànica quàntica és freqüent trobar descriptors basats en la densitat de parells o la densitat electrònica, amb un èxit divers segons les aplicacions que atenyin. Per tal de que tingui sentit químic un descriptor ha de donar la definició d'un àtom en una molècula, o ésser capaç d'identificar regions de l'espai molecular associades amb algun concepte químic (com pot ser un parell solitari o zona d'enllaç, entre d'altres). En aquesta línia, s'han proposat diversos esquemes de partició: la teoria d'àtoms en molècules (AIM), la funció de localització electrònica (ELF), les cel·les de Voroni, els àtoms de Hirshfeld, els àtoms difusos, etc. L'objectiu d'aquesta tesi és explorar descriptors de la densitat basats en particions de l'espai molecular del tipus AIM, ELF o àtoms difusos, analitzar els descriptors existents amb diferents nivells de teoria, proposar nous descriptors d'aromaticitat, així com estudiar l'habilitat de totes aquestes eines per discernir entre diferents mecanismes de reacció.

Implementation and application of basis set superposition error-correction schemes to the theoretical modeling of weak intermolecular interactions

This thesis deals with the so-called Basis Set Superposition Error (BSSE) from both a methodological and a practical point of view. The purpose of the present thesis is twofold: (a) to contribute step ahead in the correct characterization of weakly bound complexes and, (b) to shed light the understanding of the actual implications of the basis set extension effects in the ab intio calculations and contribute to the BSSE debate. The existing BSSE-correction procedures are deeply analyzed, compared, validated and, if necessary, improved. A new interpretation of the counterpoise (CP) method is used in order to define counterpoise-corrected descriptions of the molecular complexes. This novel point of view allows for a study of the BSSE-effects not only in the interaction energy but also on the potential energy surface and, in general, in any property derived from the molecular energy and its derivatives A program has been developed for the calculation of CP-corrected geometry optimizations and vibrational frequencies, also using several counterpoise schemes for the case of molecular clusters. The method has also been implemented in Gaussian98 revA10 package. The Chemical Hamiltonian Approach (CHA) methodology has been also implemented at the RHF and UHF levels of theory for an arbitrary number interacting systems using an algorithm based on block-diagonal matrices. Along with the methodological development, the effects of the BSSE on the properties of molecular complexes have been discussed in detail. The CP and CHA methodologies are used for the determination of BSSE-corrected molecular complexes properties related to the Potential Energy Surfaces and molecular wavefunction, respectively. First, the behaviour of both BSSE-correction schemes are systematically compared at different levels of theory and basis sets for a number of hydrogen-bonded complexes. The Complete Basis Set (CBS) limit of both uncorrected and CP-corrected molecular properties like stabilization energies and intermolecular distances has also been determined, showing the capital importance of the BSSE correction. Several controversial topics of the BSSE correction are addressed as well. The application of the counterpoise method is applied to internal rotational barriers. The importance of the nuclear relaxation term is also pointed out. The viability of the CP method for dealing with charged complexes and the BSSE effects on the double-well PES blue-shifted hydrogen bonds is also studied in detail. In the case of the molecular clusters the effect of high-order BSSE effects introduced with the hierarchical counterpoise scheme is also determined. The effect of the BSSE on the electron density-related properties is also addressed. The first-order electron density obtained with the CHA/F and CHA/DFT methodologies was used to assess, both graphically and numerically, the redistribution of the charge density upon BSSE-correction. Several tools like the Atoms in Molecules topologycal analysis, density difference maps, Quantum Molecular Similarity, and Chemical Energy Component Analysis were used to deeply analyze, for the first time, the BSSE effects on the electron density of several hydrogen bonded complexes of increasing size. The indirect effect of the BSSE on intermolecular perturbation theory results is also pointed out It is shown that for a BSSE-free SAPT study of hydrogen fluoride clusters, the use of a counterpoise-corrected PES is essential in order to determine the proper molecular geometry to perform the SAPT analysis.

Adsorption, orientation and thermal decomposition of copper(II) hexafluoroacetylacetonate on rutile TiO2(110)

We have investigated the adsorption and thermal decomposition of copper hexafluoroacetylacetonate (Cu-11(hfaC)(2)) on single crystal rutile TiO2(110). Low energy electron diffraction shows that room temperature saturation coverage of the Cu-II(hfac)(2) adsorbate forms an ordered (2 x 1) over-layer. X-ray and ultra-violet photoemission spectroscopy of the saturated surface were recorded as the sample was annealed in a sequential manner to reveal decomposition pathways. The results show that the molecule dissociatively adsorbs by detachment of one of the two ligands to form hfac and Cu-1(hfac) which chemisorb to the substrate at 298 K. These ligands only begin to decompose once the surface temperature exceeds 473 K where Cu core level shifts indicate metallisation. This reduction from Cu(I) to Cu(0) takes place in the absence of an external reducing agent and without disproportionation and is accompanied by the onset of decomposition of the hfac ligands. Finally, C K-edge near edge X-ray absorption fine structure experiments indicate that both the ligands adsorb aligned in the < 001 > direction and we propose a model in which the hfac ligands adsorb on the 5-fold coordinated Ti atoms and the Cu-1(hfac) moiety attaches to the bridging O atoms in a square planar geometry. The calculated tilt angle for these combined geometries is approximately 10 degrees to the surface normal.

Electronic structure of point defects in controlled self-doping of the TiO2 (110) surface: Combined photoemission spectroscopy and density functional theory study

Point defects in metal oxides such as TiO2 are key to their applications in numerous technologies. The investigation of thermally induced nonstoichiometry in TiO2 is complicated by the difficulties in preparing and determining a desired degree of nonstoichiometry. We study controlled self-doping of TiO2 by adsorption of 1/8 and 1/16 monolayer Ti at the (110) surface using a combination of experimental and computational approaches to unravel the details of the adsorption process and the oxidation state of Ti. Upon adsorption of Ti, x-ray and ultraviolet photoemission spectroscopy (XPS and UPS) show formation of reduced Ti. Comparison of pure density functional theory (DFT) with experiment shows that pure DFT provides an inconsistent description of the electronic structure. To surmount this difficulty, we apply DFT corrected for on-site Coulomb interaction (DFT+U) to describe reduced Ti ions. The optimal value of U is 3 eV, determined from comparison of the computed Ti 3d electronic density of states with the UPS data. DFT+U and UPS show the appearance of a Ti 3d adsorbate-induced state at 1.3 eV above the valence band and 1.0 eV below the conduction band. The computations show that the adsorbed Ti atom is oxidized to Ti2+ and a fivefold coordinated surface Ti atom is reduced to Ti3+, while the remaining electron is distributed among other surface Ti atoms. The UPS data are best fitted with reduced Ti2+ and Ti3+ ions. These results demonstrate that the complexity of doped metal oxides is best understood with a combination of experiment and appropriate computations.

Comparison of subcellular partitioning, distribution, and internal speciation of Cu between Cu-tolerant and naive populations of Dendrodrilus rubidus Savigny

When considering contaminated site ecology and ecological risk assessment a key question is whether organisms that appear unaffected by accumulation of contaminants are tolerant or resistant to those contaminants. A population of Dendrodrilus rubidus Savigny earthworms from the Coniston Copper Mines, an area of former Cu mining, exhibit increased tolerance and accumulation of Cu relative to a nearby non-Cu exposed population. Distribution of total Cu between different body parts (posterior, anterior, body wall) of the two populations was determined after a 14 day exposure to 250 mg Cu kg(-1) in Cu-amended soil. Cu concentrations were greater in Coniston earthworms but relative proportions of Cu in different body parts were the same between populations. Cu speciation was determined using extended X-ray absorption fine structure spectroscopy (EXAFS). Cu was coordinated to 0 atoms in the exposure soil but to S atoms in the earthworms. There was no difference in this speciation between the different earthworm populations. In another experiment earthworms were exposed to a range of Cu concentrations (200-700 mg Cu kg(-1)). Subcellular partitioning of accumulated Cu was determined. Coniston earthworms accumulated more Cu but relative proportions of Cu in the different fractions (cytosol > granular > tissue fragments, cell membranes, and intact cells) were the same between populations. Results suggest that Coniston D. rubidus are able to survive in the Cu-rich Coniston Copper Mines soil through enlargement of the same Cu storage reservoirs that exist in a nearby non-Cu exposed population.

Model catalyst studies of the strong metal-support interaction: Surface structure identified by STM on Pd nanoparticles on TiO2(110)

Model catalysts of Pd nanoparticles and films on TiO2 (I 10) were fabricated by metal vapour deposition (MVD). Molecular beam measurements show that the particles are active for CO adsorption, with a global sticking probability of 0.25, but that they are deactivated by annealing above 600 K, an effect indicative of SMSI. The Pd nanoparticles are single crystals oriented with their (I 11) plane parallel to the surface plane of the titania. Analysis of the surface by atomic resolution STM shows that new structures have formed at the surface of the Pd nanoparticles and films after annealing above 800 K. There are only two structures, a zigzag arrangement and a much more complex "pinwheel" structure. The former has a unit cell containing 7 atoms, and the latter is a bigger unit cell containing 25 atoms. These new structures are due to an overlayer of titania that has appeared on the surface of the Pd nanoparticles after annealing, and it is proposed that the surface layer that causes the SMSI effect is a mixed alloy of Pd and Ti, with only two discrete ratios of atoms: Pd/Ti of 1: 1 (pinwheel) and 1:2 (zigzag). We propose that it is these structures that cause the SMSI effect. (c) 2005 Elsevier Inc. All rights reserved.

Theory of enhanced magnetic anisotropy induced by Pt adatoms on two-dimensional Co clusters supported on a Pt(111) substrate

Calculations are reported of the magnetic anisotropy energy of two-dimensional (2D) Co nanostructures on a Pt(111) substrate. The perpendicular magnetic anisotropy (PMA) of the 2D Co clusters strongly depends on their size and shape, and rapidly decreases with increasing cluster size. The PMA calculated is in reasonable agreement with experimental results. The sensitivity of the results to the Co-Pt spacing at the interface is also investigated and, in particular, for a complete Co monolayer we note that the value of the spacing at the interface determines whether PMA or in-plane anisotropy occurs. We find that the PMA can be greatly enhanced by the addition of Pt adatoms to the top surface of the 2D Co clusters. A single Pt atom can induce in excess of 5 meV to the anisotropy energy of a cluster. In the absence of the Pt adatoms the PMA of the Co clusters falls below 1 meV/Co atom for clusters of about 10 atoms whereas, with Pt atoms added to the surface of the clusters, a PMA of 1 meV/Co atom can be maintained for clusters as large as about 40 atoms. The effect of placing Os atoms on the top of the Co clusters is also considered. The addition of 5d atoms and clusters on the top of ferromagnetic nanoparticles may provide an approach to tune the magnetic anisotropy and moment separately.

Ab initio and tight-binding calculations of noncollinear magnetism in manganese clusters

Combining ab initio and tight-binding calculations, we have studied the noncollinear magnetism in manganese clusters. The oscillations in the per-atom moments observed experimentally are reproduced theoretically. The tendency of antiferromagnetic coupling between near neighbors leads to noncollinear coupling between atoms within the clusters. For clusters containing 12, 13, 15, 19, and 23 atoms, the geometrical structures were optimized from ab initio calculations with collinear coupled spin moments among different atomic sites. For larger clusters such as Mn-36 and Mn-55, the geometries are taken as portions of an fcc structure. Although the local atomic moments have high values close to 4 mu(B), the net moments lie in the range of 0.4-1.2 mu(B)/atom. Taking the noncollinear coupling into account brings the calculated magnetic moments much closer to the experimental results.

On defining base units in terms of fundamental constants

The definitions of the base units of the international system of units have been revised many times since the idea of such an international system was first conceived at the time of the French revolution. The objective today is to define all our units in terms of 'invariants of nature', i.e. by referencing our units to the fundamental constants of physics, or the properties of atoms, rather than the characteristics of our planet or of artefacts. This situation is reviewed, particularly in regard to finding a new definition of the kilogram to replace its present definition in terms of a prototype material artefact.

Microwave spectrum of 2-aminopyridine

The microwave spectra of 2-aminopyridine-NH2, -ND2, and of both of the two possible -NHD molecules have been observed and assigned in the 0+ vibrational state of the amino group inversion vibration; the assignment for three of the molecules in the 0− state is also made. From intensity measurements the 0+-0− splitting is estimated to be 135 ± 25 cm−1 for the -NH2 molecule and 95 ± 30 cm−1 for the -ND2 molecule. The rotational constants are interpreted in terms of a structure in which the amino group is bent about 32° out of the molecular plane, the c coordinates of the two amino H atoms being 0.21 and 0.28 Å. Stark effect measurements give a dipole moment of about 0.9 D which is almost entirely in the b axis, and which changes quite significantly between the 0+ and 0− states.

A potential energy surface for the ground state of formaldehyde, H2CO(1A1)

A model potential energy function for the ground state of H2CO has been derived which covers the whole space of the six internal coordinates. This potential reproduces the experimental energy, geometry and quadratic force field of formaldehyde, and dissociates correctly to all possible atom, diatom and triatom fragments. Thus there are good reasons for believing it to be close to the true potential energy surface except in regions where both hydrogen atoms are close to the oxygen. It leads to the prediction that there should be a metastable singlet hydroxycarbene HCOH which has a planar trans structure and an energy of 2•31 eV above that of equilibrium formaldehyde. The reaction path for dissociation into H2 + CO is predicted to pass through a low symmetry transition state with an activation energy of 4•8 eV. Both of these predictions are in good agreement with recently published ab initio calculations.

Intensities in infrared bands V: Bond polar properties in dimethylacetylene

The absolute intensities of all except one of the infra-red fundamental vibration bands of dimethyl acetylene have been determined, and the results have been used to compute polar properties of the C—H and C—C bonds. It has been found that if the very probable assumption is made that the acetylenic carbon atoms carry a residual negative charge, the hydrogen atoms in the C—H bonds must carry a residual positive charge. The probable value of the C—H dipole is about 04 Debye, and that of the C—C bond about 1 Debye. Comparisons have been made with the results of similar work with related molecules.

ASYM20: a program for force constant and normal coordinate calculations, with a critical review of the theory involved

The theory of harmonic force constant refinement calculations is reviewed, and a general-purpose program for force constant and normal coordinate calculations is described. The program, called ASYM20. is available through Quantum Chemistry Program Exchange. It will work on molecules of any symmetry containing up to 20 atoms and will produce results on a series of isotopomers as desired. The vibrational secular equations are solved in either nonredundant valence internal coordinates or symmetry coordinates. As well as calculating the (harmonic) vibrational wavenumbers and normal coordinates, the program will calculate centrifugal distortion constants, Coriolis zeta constants, harmonic contributions to the α′s. root-mean-square amplitudes of vibration, and other quantities related to gas electron-diffraction studies and thermodynamic properties. The program will work in either a predict mode, in which it calculates results from an input force field, or in a refine mode, in which it refines an input force field by least squares to fit observed data on the quantities mentioned above. Predicate values of the force constants may be included in the data set for a least-squares refinement. The program is written in FORTRAN for use on a PC or a mainframe computer. Operation is mainly controlled by steering indices in the input data file, but some interactive control is also implemented.

Characterization of a gamma-aminobutyric acid transport system of rhizobium leguminosarum bv. viciae 3841

Spontaneous mutants of Rhizobium leguminosarum bv. viciae 3841 were isolated that grow faster than the wild type on gamma-aminobutyric acid (GABA) as the sole carbon and nitrogen source. These strains (RU1736 and RU1816) have frameshift mutations (gtsR101 and gtsR102, respectively) in a GntR-type regulator (GtsR) that result in a high rate of constitutive GABA transport. Tn5 mutagenesis and quantitative reverse transcription-PCR showed that GstR regulates expression of a large operon (pRL100242 to pRL100252) on the Sym plasmid that is required for GABA uptake. An ABC transport system, GtsABCD (for GABA transport system) (pRL100248-51), of the spermidine/putrescine family is part of this operon. GtsA is a periplasmic binding protein, GtsB and GtsC are integral membrane proteins, and GtsD is an ATP-binding subunit. Expression of gtsABCD from a lacZ promoter confirmed that it alone is responsible for high rates of GABA transport, enabling rapid growth of strain 3841 on GABA. Gts transports open-chain compounds with four or five carbon atoms with carboxyl and amino groups at, or close to, opposite termini. However, aromatic compounds with similar spacing between carboxyl and amino groups are excellent inhibitors of GABA uptake so they may also be transported. In addition to the ABC transporter, the operon contains two putative mono-oxygenases, a putative hydrolase, a putative aldehyde dehydrogenase, and a succinate semialdehyde dehydrogenase. This suggests the operon may be involved in the transport and breakdown of a more complex precursor to GABA. Gts is not expressed in pea bacteroids, and gtsB mutants are unaltered in their symbiotic phenotype, suggesting that Bra is the only GABA transport system available for amino acid cycling.

Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favorability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB)versus Thr(HB)-Asn(nHB). Significantly (p < or = 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi1 rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi1 conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi1 and chi2 conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favorability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http://www.rubic.rdg.ac.uk/betapairprefsparallel/).