Controlling the spin of metal atoms adsorbed on oxide surfaces: Ni on regular and defective sites of the MgO(001) surface

We have analyzed the relative energy of nonmagnetic and magnetic low-lying electronic states of Ni atoms adsorbed on regular and defective sites of the MgO(001) surface. To this end cluster and periodic surface models are used within density functional theory. For Ni atoms adsorbed on oxygen vacancies at low coverage, the interaction energy between the metal and the support is much larger than on regular sites. Strong bonding results in a diamagnetic adsorbed species and the energy required to reach the high-spin state increases. Moreover, a correlation appears between the low-spin to high-spin energy difference and the interaction energy hypothesizing that it is possible to prepare the surface to tune the high-spin to low-spin energy difference. Magnetic properties of adsorbed thin films obtained upon increasing coverage are more difficult to interpret. This is because the metallic bond is readily formed and dominates over the effect of the atoms directly bound to the vacancy.

Analytical Dirac-Hartree-Fock-Slater screening function for atoms (Z=1-92)

An analytical approximation, depending on five parameters, for the atomic screening function is proposed. The corresponding electrostatic potential takes a simple analytical form (superposition of three Yukawa potentials) well suited to most practical applications. Parameters in the screening function, determined by an analytical fitting procedure to Dirac-Hartree-Fock-Slater (DHFS) self-consistent data, are given for Z=1¿92. The reliability of this analytical approach is demonstrated by showing that (a) Born cross sections for elastic scattering of fast charged particles by the present analytical field and by the DHFS field practically coincide and (b) one-electron binding energies computed from the independent-particle model with our analytical field (corrected for exchange and electrostatic self-interaction) agree closely with the DHFS energy eigenvalues.

Dinàmica de sistemes mesoscòpics amb correlacions quàntiques mitjançant la formulació de de Broglie-Bohm

Per a determinar la dinàmica espai-temporal completa d’un sistema quàntic tridimensional de N partícules cal integrar l’equació d’Schrödinger en 3N dimensions. La capacitat dels ordinadors actuals permet fer-ho com a molt en 3 dimensions. Amb l’objectiu de disminuir el temps de càlcul necessari per a integrar l’equació d’Schrödinger multidimensional, es realitzen usualment una sèrie d’aproximacions, com l’aproximació de Born–Oppenheimer o la de camp mig. En general, el preu que es paga en realitzar aquestes aproximacions és la pèrdua de les correlacions quàntiques (o entrellaçament). Per tant, és necessari desenvolupar mètodes numèrics que permetin integrar i estudiar la dinàmica de sistemes mesoscòpics (sistemes d’entre tres i unes deu partícules) i en els que es tinguin en compte, encara que sigui de forma aproximada, les correlacions quàntiques entre partícules. Recentment, en el context de la propagació d’electrons per efecte túnel en materials semiconductors, X. Oriols ha desenvolupat un nou mètode [Phys. Rev. Lett. 98, 066803 (2007)] per al tractament de les correlacions quàntiques en sistemes mesoscòpics. Aquesta nova proposta es fonamenta en la formulació de la mecànica quàntica de de Broglie– Bohm. Així, volem fer notar que l’enfoc del problema que realitza X. Oriols i que pretenem aquí seguir no es realitza a fi de comptar amb una eina interpretativa, sinó per a obtenir una eina de càlcul numèric amb la que integrar de manera més eficient l’equació d’Schrödinger corresponent a sistemes quàntics de poques partícules. En el marc del present projecte de tesi doctoral es pretén estendre els algorismes desenvolupats per X. Oriols a sistemes quàntics constituïts tant per fermions com per bosons, i aplicar aquests algorismes a diferents sistemes quàntics mesoscòpics on les correlacions quàntiques juguen un paper important. De forma específica, els problemes a estudiar són els següents: (i) Fotoionització de l’àtom d’heli i de l’àtom de liti mitjançant un làser intens. (ii) Estudi de la relació entre la formulació de X. Oriols amb la aproximació de Born–Oppenheimer. (iii) Estudi de les correlacions quàntiques en sistemes bi- i tripartits en l’espai de configuració de les partícules mitjançant la formulació de de Broglie–Bohm.

Eina per a la visualització de molècules

Aquest TFC consisteix a construir una eina que permet visualitzar molècules en aquest format de visualització. L'eina haurà de permetre navegar al voltant de la molècula i emprar tècniques de selecció que permetin identificar àtoms, calcular distàncies i angles de torsió. A més, s'haurà de poder definir un eix i sobre aquest eix generar una rotació i gravar-la en una pel·lícula (en el format que sigui més senzill).

Relacions estructura-propietat per un conjunt d'hidrocarburs a partir de nous descriptors tridimensionals derivats de la semblança molecular quàntica

En aquest article es defineixen uns nous índexs tridimensionals per a la descripció de les molècules a partir de paràmetres derivats de la Teoria de la Semblança Molecular i de les distàncies euclidianes entre els àtoms i les càrregues atòmiques efectives. Aquests indexs,anomenats 3D, s'han aplicat a l'estudi de les relacions estructura-propietat d'una família d'hidrocarburs, i han demostrat una capacitat de descripció de tres propietats de la família (temperatura d'ebullició, temperatura de fusió i densitat) molt més acurada que quan s'utilitzen els indexs 2D clàssics

Growth and characterization of CoO ultra thin films

In this report we present the growth process of the cobalt oxide system using reactive electron beam deposition. In that technique, a target of metallic cobalt is evaporated and its atoms are in-flight oxidized in an oxygen rich reactive atmosphere before reaching the surface of the substrate. With a trial and error procedure the deposition parameters have been optimized to obtain the correct stoichiometry and crystalline phase. The evaporation conditions to achieve the correct cobalt oxide salt rock structure, when evaporating over amorphous silicon nitride, are: 525 K of substrate temperature, 2.5·10-4 mbar of oxygen partial pressure and 1 Å/s of evaporation rate. Once the parameters were optimized a set of ultra thin film ranging from samples of 1 nm of nominal thickness to 20nm thick and bulk samples were grown. With the aim to characterize the samples and study their microstructure and morphology, X-ray diffraction, transmission electron microscopy, electron diffraction, energy dispersive X-ray spectroscopy and quasi-adiabatic nanocalorimetry techniques are utilised. The final results show a size dependent effect of the antiferromagnetic transition. Its Néel temperature becomes depressed as the size of the grains forming the layer decreases.

Sobre els mapes de semblança quàntica molecular

Partint de les definicions usuals de Mesures de Semblança Quàntica (MSQ), es considera la dependència d'aquestes mesures respecte de la superposició molecular. Pel cas particular en qnè els sistemes comparats siguin una molècula i un Àtom i que les mesures es calculin amb l'aproximació EASA, les MSQ esdevenen funcions de les tres coordenades de l'espai. Mantenint fixa una de les tres coordenades, es pot representar fàcilment la variació del valor de semblança en un pla determinat, i obtenir els anomenats mapes de semblança. En aquest article, es comparen els mapes de semblança obtinguts amb diferents MSQ per a sistemes senzills

Modelling of Mercury’s surface composition and remote detection from the orbit with the BepiColombo Mercury Planetary Orbiter

It can be assumed that the composition of Mercury’s thin gas envelope (exosphere) is related to thecomposition of the planets crustal materials. If this relationship is true, then inferences regarding the bulkchemistry of the planet might be made from a thorough exospheric study. The most vexing of allunsolved problems is the uncertainty in the source of each component. Historically, it has been believedthat H and He come primarily from the solar wind, while Na and K originate from volatilized materialspartitioned between Mercury’s crust and meteoritic impactors. The processes that eject atoms andmolecules into the exosphere of Mercury are generally considered to be thermal vaporization, photonstimulateddesorption (PSD), impact vaporization, and ion sputtering. Each of these processes has its owntemporal and spatial dependence. The exosphere is strongly influenced by Mercury’s highly ellipticalorbit and rapid orbital speed. As a consequence the surface undergoes large fluctuations in temperatureand experiences differences of insolation with longitude. We will discuss these processes but focus moreon the expected surface composition and solar wind particle sputtering which releases material like Caand other elements from the surface minerals and discuss the relevance of composition modelling

Calorimetry of dehydrogenation and dangling-bond recombination in several hydrogenated amorphous silicon materials

Differential scanning calorimetry (DSC) was used to study the dehydrogenation processes that take place in three hydrogenated amorphous silicon materials: nanoparticles, polymorphous silicon, and conventional device-quality amorphous silicon. Comparison of DSC thermograms with evolved gas analysis (EGA) has led to the identification of four dehydrogenation processes arising from polymeric chains (A), SiH groups at the surfaces of internal voids (A'), SiH groups at interfaces (B), and in the bulk (C). All of them are slightly exothermic with enthalpies below 50 meV/H atoms , indicating that, after dissociation of any SiH group, most dangling bonds recombine. The kinetics of the three low-temperature processes [with DSC peak temperatures at around 320 (A),360 (A'), and 430°C (B)] exhibit a kinetic-compensation effect characterized by a linea relationship between the activation entropy and enthalpy, which constitutes their signature. Their Si-H bond-dissociation energies have been determined to be E (Si-H)0=3.14 (A), 3.19 (A'), and 3.28 eV (B). In these cases it was possible to extract the formation energy E(DB) of the dangling bonds that recombine after Si-H bond breaking [0.97 (A), 1.05 (A'), and 1.12 (B)]. It is concluded that E(DB) increases with the degree of confinement and that E(DB)>1.10 eV for the isolated dangling bond in the bulk. After Si-H dissociation and for the low-temperature processes, hydrogen is transported in molecular form and a low relaxation of the silicon network is promoted. This is in contrast to the high-temperature process for which the diffusion of H in atomic form induces a substantial lattice relaxation that, for the conventional amorphous sample, releases energy of around 600 meV per H atom. It is argued that the density of sites in the Si network for H trapping diminishes during atomic diffusion

Anharmonicity contributions to the vibrational second hyperpolarizability of conjugated oligomers

Restricted Hartree-Fock 6-31G calculations of electrical and mechanical anharmonicity contributions to the longitudinal vibrational second hyperpolarizability have been carried out for eight homologous series of conjugated oligomers - polyacetylene, polyyne, polydiacetylene, polybutatriene, polycumulene, polysilane, polymethineimine, and polypyrrole. To draw conclusions about the limiting infinite polymer behavior, chains containing up to 12 heavy atoms along the conjugated backbone were considered. In general, the vibrational hyperpolarizabilities are substantial in comparison with their static electronic counterparts for the dc-Kerr and degenerate four-wave mixing processes (as well as for static fields) but not for electric field-induced second harmonic generation or third harmonic generation. Anharmonicity terms due to nuclear relaxation are important for the dc-Kerr effect (and for the static hyperpolarizability) in the σ-conjugated polymer, polysilane, as well as the nonplanar π systems polymethineimine and polypyrrole. Restricting polypyrrole to be planar, as it is in the crystal phase, causes these anharmonic terms to become negligible. When the same restriction is applied to polymethineimine the effect is reduced but remains quantitatively significant due to the first-order contribution. We conclude that anharmonicity associated with nuclear relaxation can be ignored, for semiquantitative purposes, in planar π-conjugated polymers. The role of zero-point vibrational averaging remains to be evaluated

Basis set superposition error-counterpoise corrected potential energy surfaces : application to hydrogen peroxide ... X (X=F-, Cl-, Br-, Li+, Na+) complexes

Møller-Plesset (MP2) and Becke-3-Lee-Yang-Parr (B3LYP) calculations have been used to compare the geometrical parameters, hydrogen-bonding properties, vibrational frequencies and relative energies for several X- and X+ hydrogen peroxide complexes. The geometries and interaction energies were corrected for the basis set superposition error (BSSE) in all the complexes (1-5), using the full counterpoise method, yielding small BSSE values for the 6-311 + G(3df,2p) basis set used. The interaction energies calculated ranged from medium to strong hydrogen-bonding systems (1-3) and strong electrostatic interactions (4 and 5). The molecular interactions have been characterized using the atoms in molecules theory (AIM), and by the analysis of the vibrational frequencies. The minima on the BSSE-counterpoise corrected potential-energy surface (PES) have been determined as described by S. Simón, M. Duran, and J. J. Dannenberg, and the results were compared with the uncorrected PES

The mapping of the local contributions of Fermi and Coulomb correlation into intracule and extracule density distributions

The contributions of the correlated and uncorrelated components of the electron-pair density to atomic and molecular intracule I(r) and extracule E(R) densities and its Laplacian functions ∇2I(r) and ∇2E(R) are analyzed at the Hartree-Fock (HF) and configuration interaction (CI) levels of theory. The topologies of the uncorrelated components of these functions can be rationalized in terms of the corresponding one-electron densities. In contrast, by analyzing the correlated components of I(r) and E(R), namely, IC(r) and EC(R), the effect of electron Fermi and Coulomb correlation can be assessed at the HF and CI levels of theory. Moreover, the contribution of Coulomb correlation can be isolated by means of difference maps between IC(r) and EC(R) distributions calculated at the two levels of theory. As application examples, the He, Ne, and Ar atomic series, the C2-2, N2, O2+2 molecular series, and the C2H4 molecule have been investigated. For these atoms and molecules, it is found that Fermi correlation accounts for the main characteristics of IC(r) and EC(R), with Coulomb correlation increasing slightly the locality of these functions at the CI level of theory. Furthermore, IC(r), EC(R), and the associated Laplacian functions, reveal the short-ranged nature and high isotropy of Fermi and Coulomb correlation in atoms and molecules

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

Nuclear magnetic resonance chemical shifts with the statistical average of orbital-dependent model potentials in Kohn-Sham density functional theory

The performance of the SAOP potential for the calculation of NMR chemical shifts was evaluated. SAOP results show considerable improvement with respect to previous potentials, like VWN or BP86, at least for the carbon, nitrogen, oxygen, and fluorine chemical shifts. Furthermore, a few NMR calculations carried out on third period atoms (S, P, and Cl) improved when using the SAOP potential

Two complementary molecular energy decomposition schemes: the Mayer and Ziegler-Rauk methods in comparison

In the present paper we discuss and compare two different energy decomposition schemes: Mayer's Hartree-Fock energy decomposition into diatomic and monoatomic contributions [Chem. Phys. Lett. 382, 265 (2003)], and the Ziegler-Rauk dissociation energy decomposition [Inorg. Chem. 18, 1558 (1979)]. The Ziegler-Rauk scheme is based on a separation of a molecule into fragments, while Mayer's scheme can be used in the cases where a fragmentation of the system in clearly separable parts is not possible. In the Mayer scheme, the density of a free atom is deformed to give the one-atom Mulliken density that subsequently interacts to give rise to the diatomic interaction energy. We give a detailed analysis of the diatomic energy contributions in the Mayer scheme and a close look onto the one-atom Mulliken densities. The Mulliken density ρA has a single large maximum around the nuclear position of the atom A, but exhibits slightly negative values in the vicinity of neighboring atoms. The main connecting point between both analysis schemes is the electrostatic energy. Both decomposition schemes utilize the same electrostatic energy expression, but differ in how fragment densities are defined. In the Mayer scheme, the electrostatic component originates from the interaction of the Mulliken densities, while in the Ziegler-Rauk scheme, the undisturbed fragment densities interact. The values of the electrostatic energy resulting from the two schemes differ significantly but typically have the same order of magnitude. Both methods are useful and complementary since Mayer's decomposition focuses on the energy of the finally formed molecule, whereas the Ziegler-Rauk scheme describes the bond formation starting from undeformed fragment densities