Search

Formic acid (HCOOH) adsorption on rutile TiO2 (110) has been studied by s- and p-polarized infrared reflection-absorption spectroscopy (IRRAS) and spin-polarized density functional theory together with Hubbard U contributions (DFT+U) calculations. To compare with IRRAS spectra, the results from the DFT+U calculations were used to simulate IR spectra by employing a threelayer model, where the adsorbate layer was modelled using Lorentz oscillators with calculated dielectric constants. To account for the experimental observations, four possible formate adsorption geometries were calculated, describing both the perfect (110) surface, and surfaces with defects; either O vacancies or hydroxyls. The majority species seen in IRRAS was confirmed to be the bridging bidentate formate species with associated symmetric and asymmetric frequencies of the ν(OCO) modes measured to be at 1359 cm-1 and 1534 cm-1, respectively. The in-plane δ(C-H) wagging mode of this species couples to both the tangential and the normal component of the incident p-polarized light, which results in absorption and emission bands at 1374 cm-1 and 1388 cm-1. IRRAS spectra measured on surfaces prepared to be either reduced, stoichiometric, or to contain surplus O adatoms, were found to be very similar. By comparisons with computed spectra, it is proposed that in our experiments, formate binds as a minority species to an in-plane Ti5c atom and a hydroxyl, rather than to O vacancy sites, the latter to a large extent being healed even at our UHV conditions. Excellent agreement between calculated and experimental IRRAS spectra is obtained. The results emphasize the importance of protonation and reactive surface hydroxyls - even under UHV conditions - as reactive sites in e.g., catalytic applications. [ABSTRACT FROM AUTHOR]

We examine the real space structure and the electronic structure (particularly Ce4f electron localization) of oxygen vacancies in CeO2 (ceria) as a function of U in density functional theory studies with the rotationally invariant forms of the LDA+U and GGA+U functionals. The four nearest neighbor Ce ions always relax outwards, with those not carrying localized Ce4f charge moving furthest. Several quantification schemes show that the charge starts to become localized at U≈3 eV and that the degree of localization reaches a maximum at ∼6 eV for LDA+U or at ∼5.5 eV for GGA+U. For higher U it decreases rapidly as charge is transferred onto second neighbor O ions and beyond. The localization is never into atomic corelike states; at maximum localization about 80–90% of the Ce4f charge is located on the two nearest neighboring Ce ions. However, if we look at the total atomic charge we find that the two ions only make a net gain of (0.2–0.4)e each, so localization is actually very incomplete, with localization of Ce4f electrons coming at the expense of moving other electrons off the Ce ions. We have also revisited some properties of defect-free ceria and find that with LDA+U the crystal structure is actually best described with U=3–4 eV, while the experimental band structure is obtained with U=7–8 eV. (For GGA+U the lattice parameters worsen for U>0 eV, but the band structure is similar to LDA+U.) The best overall choice is U≈6 eV with LDA+U and ≈5.5 eV for GGA+U, since the localization is most important, but a consistent choice for both CeO2 and Ce2O3, with and without vacancies, is hard to find. [ABSTRACT FROM AUTHOR]

Electron impact ionization cross sections measured close to threshold are reported for both the monomers and dimers of ozone and OC1O using a new high resolution electron impact apparatus. The present appearance energies AE(O[sup +, sub 3]-/O[sub 3])=12.70±0.02 eV, AE (OC1O[sup +]/OC1O)= 10.55±0.02 and AE(C1O[sup +]/OC10 = 13.37±0.03 eV derived from the measured ionization cross sections are in excellent agreement with the vertical threshold values determined for these ions by high resolution PES and PIMS photoionization studies. The corresponding appearance energies determined for the dimer ions, 10.10±0.3 eV for (O[sub 3])[sup +, sub 2] and 9.87±0.2 eV for (OC1O)[sup +, sub 2], are both red shifted with respect to the monomer case. The bond energy (0.70-0.3+0.5) eV of (OC1O)[sup +, sub 2] estimated from these data is similar to that of other dimer ions, whereas the bond energy of (O[sub 3]-O[sup +, sub 3]) with (2.55-0.4+0.6) eV is rather large suggesting an unusual structure for the cationic ozone dimer ion. Based on quantum chemical calculations on various levels we are led to the conclusion that the ion produced by ionization of the ozone dimer is no longer a conventional dimer ion where the two monomer units are still present (as is the case for the OC1O system), but rather an ion of form O[sub 2]...O[sup +, sub 4] or a twisted ring structure of (O[sub 6])[sup +]. [ABSTRACT FROM AUTHOR]

Mg-doped VO2 thin films with thermochromic properties were made by reactive DC magnetron co-sputtering onto heated substrates, and spectral absorption was recorded at room temperature in the 0.5 < ħω < 3.5 eV energy range. Clear evidence was found for a widening of the main band gap from 1.67 to 2.32 eV as the Mg/(V + Mg) atomic ratio went from zero to 0.19, thereby significantly lowering the luminous absorption. This technologically important effect could be reconciled with spin-polarized density functional theory calculations using the Heyd-Scuseria-Ernzerhof [Heyd et al., J. Chem. Phys. 118, 8207 (2003); ibid. 124, 219906 (2006)] hybrid functional. Specifically, the calculated luminous absorptance decreased when the Mg/(V + Mg) ratio was increased from 0.125 to 0.250. [ABSTRACT FROM AUTHOR]

It is known that, in the proximity of a neutral wall, liquids experience diffusion enhancement relative to their bulk diffusion, but the origin of this phenomenon is still unknown. We report a molecular dynamics simulation investigating the dynamics of a simple liquid in the proximity to a non-interacting smooth confining wall, which exhibits a strong diffusion enhancement within the liquid layers adjacent to the wall. We present an analysis of these results, demonstrating that the observed diffusion enhancement can be accounted for, with numerical accuracy, using the universal scaling law that relates the liquid diffusion rate to the excess entropy. These results show that the scaling law, which has so far only been used for the description of the bulk liquid diffusion, can be successfully used to describe the diffusion in liquids under nano-scale confinement., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International (CC BY-NC-ND) license (https://creativecommons.org/licenses/by-nc-nd/4.0/).)

High-accuracy molecular force field models suited for hot gases and plasmas are not as abundant as those geared toward ambient pressure and temperature conditions. Here, we present an improved version of our previous electron-density based force field model that can now account for polarization effects by adjusting the atomic valence electron contributions to match ab initio calculated Mulliken partial charges. Using a slightly modified version of the Hohenberg-Kohn theorem, we also include an improved theoretical formulation of our model when applied to systems with degenerate ground states. We present two variants of our polarizable model, fitted from ab initio reference data calculated at CCSD(T)/cc-pVTZ and CCSD(T)/CEP-31G levels of theory, that both accurately model water dimer interaction energies. Further improvements include the additional interaction components with fictitious non-spherically symmetric, yet atom-centered, electron densities and fitting the exchange and correlation coefficients against analytical expressions. The latter removes all unphysical oscillations that are observed in the previous non-polarizable variant of our force field., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution-NonCommercial 4.0 International (CC BY-NC) license (https://creativecommons.org/licenses/by-nc/4.0/).)

The H2O and H2O2 molecules resemble each other in a multitude of ways as has been noted in the literature. Here, we present density functional theory (DFT) calculations for the H2O2(s) and H2O2·2H2O(s) crystals and make selected comparisons with ice polymorphs. The performance of a number of dispersion-corrected density functionals-both self-consistent and a posteriori ones-are assessed, and we give special attention to the D3 correction and its effects. The D3 correction to the lattice energies is large: for H2O2(s) the D3 correction constitutes about 25% of the lattice energy using PBE, much more for RPBE, much less for SCAN, and it primarily arises from non-H-bonded interactions out to about 5 Å.The large D3 corrections to the lattice energies are likely a consequence of several effects: correction for missing dispersion interaction, the ability of D3 to capture and correct various other kinds of limitations built into the underlying DFT functionals, and finally some degree of cell-contraction-induced polarization enhancement. We find that the overall best-performing functionals of the twelve examined are optPBEvdW and RPBE-D3. Comparisons with DFT assessments for ices in the literature show that where the same methods have been used, the assessments largely agree., (© 2023 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

Theoretical and experimental vibrational signatures of H2O and OH- (dissociated water) adsorbed on stoichiometric ceria{111} surfaces are compared. The experimental ones were collected from low-coverage experiments in the literature, and the theoretical anharmonic frequencies were generated using density functional theory calculations employing the optPBE-vdW functional for coverages from 0.5 to a few monolayers. It is found that (i) the experiments and our calculations overall agree well, lending credibility to both; (ii) the calculations manage to resolve the large class of H-bonded motifs into frequency classes that can guide experimental assignments; (iii) it is possible to find a geometrical H-bond definition that also captures the OH vibrational frequency downshifts well: R(H⋯O) ≤2.5 Å and the O-H⋯O angle θ ≥ 100°; and (iv) the frequency vs electric field relations for water and hydroxides (i.e., dissociated water) follow different and well-separated curves., (© 2023 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

A computational protocol is developed for efficient studies of partially reduced redox-active oxides using the self-consistent charge density functional tight-binding method. The protocol is demonstrated for ceria, which is a prototypical reducible oxide material. The underlying idea is to achieve a consistent (and harmonized) set of Slater-Koster (SK) tables with connected repulsive potentials that enable switching on and off the in-valence description of the Ce 4f states without serious loss of accuracy in structure and energetics. The implicit treatment of the Ce 4f states, with the use of f-in-core SK-tables, is found to lead to a significant decrease in computational time. More importantly, it allows for explicit control of the oxidation states of individual Ce atoms. This makes it possible to "freeze" the electronic configuration, thereby allowing the exploration of the energetics for various meta-stable configurations. We anticipate that the outlined strategy can help to shed light on the interplay between the size, shape, and redox activity for nanoceria and other related materials.

We report an ab initio molecular dynamics (MD) simulation investigating the effect of a fully hydrated surface of TiO 2 on the water dynamics. It is found that the universal relation between the rotational and translational diffusion characteristics of bulk water is broken in the water layers near the surface with the rotational diffusion demonstrating progressive retardation relative to the translational diffusion when approaching the surface. This kind of rotation-translation decoupling has so far only been observed in the supercooled liquids approaching glass transition, and its observation in water at a normal liquid temperature is of conceptual interest. This finding is also of interest for the application-significant studies of the water interaction with fully hydrated nanoparticles. We note that this is the first observation of rotation-translation decoupling in an ab initio MD simulation of water.

Thin film structures of water on the CeO 2 (111) surface for coverages between 0.5 and 2.0 water monolayers have been optimized and analyzed using density functional theory (optPBE-vdW functional). We present a new 1.0 ML structure that is both the lowest in energy published and features a hydrogen-bond network extending the surface in one-dimension, contrary to what has been found in the literature, and contrary to what has been expected due to the large bulk ceria cell dimension. The adsorption energies for the monolayer and multilayered water structures agree well with experimental temperature programmed desorption results from the literature, and we discuss the stability window of CeO 2 (111) surfaces covered with 0.5-2.0 ML of water.

Combining experimental spectroscopy and hybrid density functional theory calculations, we show that the incorporation of fluoride ions into a prototypical reducible oxide surface, namely, ceria(111), can induce a variety of nontrivial changes to the local electronic structure, beyond the expected increase in the number of Ce 3+ ions. Our resonant photoemission spectroscopy results reveal new states above, within, and below the valence band, which are unique to the presence of fluoride ions at the surface. With the help of hybrid density functional calculations, we show that the different states arise from fluoride ions in different atomic layers in the near surface region. In particular, we identify a structure in which a fluoride ion substitutes for an oxygen ion at the surface, with a second fluoride ion on top of a surface Ce 4+ ion giving rise to F 2p states which overlap the top of the O 2p band. The nature of this adsorbate F - -Ce 4+ resonant enhancement feature suggests that this bond is at least partially covalent. Our results demonstrate the versatility of anion doping as a potential means of tuning the valence band electronic structure of ceria.

Unraveling the atomistic details of solid/liquid interfaces, e.g., by means of vibrational spectroscopy, is of vital importance in numerous applications, from electrochemistry to heterogeneous catalysis. Water-oxide interfaces represent a formidable challenge because a large variety of molecular and dissociated water species are present at the surface. Here, we present a comprehensive theoretical analysis of the anharmonic OH stretching vibrations at the water/ZnO(101¯0) interface as a prototypical case. Molecular dynamics simulations employing a reactive high-dimensional neural network potential based on density functional theory calculations have been used to sample the interfacial structures. In the second step, one-dimensional potential energy curves have been generated for a large number of configurations to solve the nuclear Schrödinger equation. We find that (i) the ZnO surface gives rise to OH frequency shifts up to a distance of about 4 Å from the surface; (ii) the spectrum contains a number of overlapping signals arising from different chemical species, with the frequencies decreasing in the order ν(adsorbed hydroxide) > ν(non-adsorbed water) > ν(surface hydroxide) > ν(adsorbed water); (iii) stretching frequencies are strongly influenced by the hydrogen bond pattern of these interfacial species. Finally, we have been able to identify substantial correlations between the stretching frequencies and hydrogen bond lengths for all species.

In this work, a range of van der Waals type density functionals are applied to the H 2 O/NaCl(001) and H 2 O/MgO(001) interface systems to explore the effect of an explicit dispersion treatment. The functionals we use are the self-consistent vdW functionals vdW-DF, vdW-DF2, optPBE-vdW, optB88-vdW, optB86b-vdW, and vdW-DF-cx, as well as the dispersion-corrected PBE-TS and PBE-D2 methods; they are all compared with the standard PBE functional. For both NaCl(001) and MgO(001), we find that the dispersion-flavoured functionals stabilize the water-surface interface by approximately 20%-40% compared to the PBE results. For NaCl(001), where the water molecules remain intact for all overlayers, the dominant contribution to the adsorption energy from "density functional theory dispersion" stems from the water-surface interactions rather than the water-water interactions. The optPBE-vdW and vdW-DF-cx functionals yield adsorption energies in good agreement with available experimental values for both NaCl and MgO. To probe the strengths of the perturbations of the adsorbed water molecules, we also calculated water dipole moments and found an increase up to 85% for water at the MgO(001) surface and 70% at the NaCl(001) surface, compared to the gas-phase dipole moment.

We investigate the effects of anion doping with fluorine impurities on the chemistry of the CeO 2 (111) facet, using the results of DFT + U calculations. We consider three prototypical processes: the formation of oxygen vacancies, the adsorption of O 2 and H 2 O molecules, and the re-oxidation of the surface with fragments of the two molecules. We find that the first two of these processes are not strongly affected, but that the presence of F lowers the energy gained in the re-oxidation of the surface in comparison to the healing of an oxygen vacancy, by 1.47 eV in the case of O 2 (provided that the F is part of a cluster) and by 0.92 eV in the case of H 2 O. Based on these results, we suggest that F could enhance the redox chemistry of ceria by toggling between being in the surface and on the surface, effectively facilitating the release of lattice O by acting as a "place holder" for it. Finally, we find that the desorption of F as either 12F 2 or HF is energetically unfavourable, suggesting that F doped ceria should be stable in the presence of O 2 and H 2 O.

We have calculated the anharmonic OH(-)(aq) vibrational spectrum in aqueous solution with a "classical Monte Carlo simulation + QM/MM + vibrational" sequential approach. A new interaction model was used in the Monte Carlo simulations: a modified version of the charged-ring hydroxide-water model from the literature. This spectrum is compared with experiment and with a spectrum based on CPMD-generated structures, and the hydration structures and H-bonding for the two models are compared. We find that: (i) the solvent-induced frequency shift as well as the absolute OH(-) frequency are in good agreement with experiment using the two models; (ii) the Raman and IR bands are very similar, in agreement with experiment; (iii) the hydration structure and H-bonding around the ion are very different with the two ion-water interaction models (charged-ring and CPMD); (iv) a cancellation effect between different regions of the hydration shell makes the total spectra similar for the two interaction models, although their hydration structures are different; (v) the net OH(-) frequency shift is a blueshift of about +80 cm(-1) with respect to frequency of the gas-phase ion.

We present periodic "DFT+U" studies of single oxygen vacancies on the CeO(2)(110) surface using a number of different supercells, finding a range of different local minimum structures for the vacancy and its two accompanying Ce(III) ions. We find three different geometrical structures in combination with a variety of different Ce(III) localization patterns, several of which have not been studied before. The desired trapping of electrons was achieved in a two-stage optimization procedure. We find that the surface oxygen nearest to the vacancy either moves within the plane towards the vacancy, or rises out of the surface into either a symmetric or an unsymmetric bridge structure. Results are shown in seven slab geometry supercells, p(2 × 1), p(2 × 2), p(2 × 3), p(3 × 2), p(2 × 4), p(4 × 2), and p(3 × 3), and indicate that the choice of supercell can affect the results qualitatively and quantitatively. An unsymmetric bridge structure with one nearest and one next-nearest neighbour Ce(III) ion (a combination of localizations not previously found) is the ground state in all (but one) of the supercells studied here, and the relative stability of other structures depends strongly on supercell size. Within any one supercell the formation energies of the different vacancy structures differ by up to 0.5 eV, but the same structure can vary by up to ∼1 eV between supercells. Furthermore, finite size scaling suggests that the remaining errors (compared to still larger supercells) can also be ∼1 eV for some vacancy structures.

Results from first-principles calculations present a rather clear picture of the interaction of SO(2) with unreduced and partially reduced (111) and (110) surfaces of ceria. The Ce(3+)∕Ce(4+) redox couple, together with many oxidation states of S, give rise to a multitude of SO(x) species, with oxidation states from +III to +VI. SO(2) adsorbs either as a molecule or attaches via its S-atom to one or two surface oxygens to form sulfite (SO(3)(2-)) and sulfate (SO(4)(2-)) species, forming new S-O bonds but never any S-Ce bonds. Molecular adsorption is found on the (111) surface. SO(3)(2-) structures are found on both the (111) and (110) surfaces of both stoichiometric and partially reduced ceria. SO(4)(2-) structures are observed on the (110) surface together with the formation of two reduced Ce(3+) surface cations. SO(2) can also partially heal the ceria oxygen vacancies by weakening a S-O bond, when significant electron transfer from the surface (Ce4f) into the lowest unoccupied molecular orbital of the SO(2) adsorbate takes place and oxidizes the surface Ce(3+) cations. Furthermore, we propose a mechanism that could lead to monodentate sulfate formation at the (111) surface.

The anharmonic OH stretching vibrational frequencies, ν(OH), for the first-shell water molecules around the Li(+), Ca(2+), Mg(2+), and Al(3+) ions in dilute aqueous solutions have been calculated based on classical molecular dynamics (MD) simulations and quantum-mechanical (QM) calculations. For Li(+)(aq), Ca(2+)(aq), Mg(2+)(aq), and Al(3+)(aq), our calculated IR frequency shifts, Δν(OH), with respect to the gas-phase water frequency, are about -300, -350, -450, and -750 cm(-1), compared to -290, -290, -420, and -830 cm(-1) from experimental infrared (IR) studies. The agreement is thus quite good, except for the order between Li(+) and Ca(2+). Given that the polarizing field from the Ca(2+) ion ought to be larger than that from Li(+)(aq), our calculated result seems reasonable. Also the absolute OH frequencies agree well with experiment. The method we used is a sequential four-step procedure: QM(electronic) to make a force field+MD simulation+QM(electronic) for point-charge-embedded M(n+) (H(2)O)(y) (second shell) (H(2)O)(z) (third shell) clusters+QM(vibrational) to yield the OH spectrum. The many-body Ca(2+)-water force-field presented in this paper is new. IR intensity-weighting of the density-of-states frequency distributions was carried out by means of the squared dipole moment derivatives.

Anharmonic OH phonon dispersion curves have been calculated for the Mg(OH)(2) crystal. A crystal Hamiltonian was set up for the vibrational problem, where the coordinates consists of the bond lengths of two hydroxide ions in the central unit cell. Its two-dimensional potential energy surface was constructed from first principle calculations within the density functional theory approximation. Dispersion curves were calculated by diagonalizing the Hamiltonian in a basis of singly excited crystal functions. The single particle functions used to construct the crystal states were taken from a Morse oscillator basis set. These well chosen functions made it possible to restrict calculations to include only very few functions, which greatly contributed to a transparent presentation of the underlying theory. All calculations could be done analytically except for the calculation of a few integrals. We have compared our results with those of a series of harmonic lattice dynamics calculations and have found that the anharmonicity shifts the IR and Raman dispersion curves downward appreciably and slightly changes the energy differences between both curves. From an analysis of the harmonic results we conclude that incorporating the coupling between OH stretching motion and the motion of their centers of mass will appreciably change the overall features of the dispersion curves. Extension of the anharmonic model along these lines will cause no problem to the theoretical approach presented in this paper.

In this paper we evaluate the performance of density functional theory with the B3LYP functional for calculations on ceria (CeO(2)) and cerium sesquioxide (Ce(2)O(3)). We demonstrate that B3LYP is able to describe CeO(2) and Ce(2)O(3) reasonably well. When compared to other functionals, B3LYP performs slightly better than the hybrid functional PBE0 for the electronic properties but slightly worse for the structural properties, although neither performs as well as LDA+U(U=6 eV) or PBE+U(U=5 eV). We also make an extensive comparison of atomic basis sets suitable for periodic calculations of these cerium oxides. Here we conclude that there is currently only one type of cerium basis set available in the literature that is able to give a reasonable description of the electronic structure of both CeO(2) and Ce(2)O(3). These basis sets are based on a 28 electron effective core potential (ECP) and 30 electrons are attributed to the valence space of cerium. Basis sets based on 46 electron ECPs fail for these materials.

A two-dimensional quantum-mechanical vibrational model has been used to calculate the anharmonic OH vibrational frequencies in the layered Mg(OH)(2) (brucite) crystal. The underlying potential energy surface was generated by density functional theory (DFT) calculations. The resulting OH frequencies are upshifted (blueshifted) by about +75 cm(-1) with respect to the gas-phase OH frequency (+120 cm(-1) in experiments; the discrepancy is mainly due to inadequacies in the DFT and pseudopotential models). The Raman-IR split is about 50 cm(-1), both in the calculations and in experiments. We find that the blueshift phenomenon in brucite can qualitatively be explained by a parabolalike "OH frequency versus electric field" correlation curve pertaining to an OH(-) ion exposed to an electric field. We also find that it is primarily the neighbors within the Mg(OH)(2) layer that induce the blueshift while the interlayer interaction gives a smaller (and redshifting) contribution.

Anharmonic vibrational frequencies for the Raman-active (A(1g)) and the IR-active (A(2u)) modes have been calculated for the LiOH crystal within a plane-wave density functional theory (DFT) framework. We find that a two-dimensional quantum-mechanical vibrational approach, allowing for anharmonic coupling between symmetric and antisymmetric OH stretching modes, produces OH frequencies--both absolute frequencies and gas-to-solid frequency shifts--in good agreement with experiment. Remaining errors in the absolute frequencies are largely a consequence of the DFT model chosen. A one-dimensional normal-mode following vibrational treatment, on the other hand, fails to reproduce both absolute anharmonic frequencies and gas-to-solid frequency shifts.

We have combined an embedded-cluster model with an extension of the method of increments to treat the adsorption of molecules on a surface. In this way we are able to investigate the physisorption of CO on CeO(2)(110) at the MP2, MP4(SDTQ), and CCSD(T) levels with only moderate computational costs. We find that, at the CCSD(T) level, 25% of the adsorption energy originates from electron correlation. The interactions of the CO molecule with its five nearest cerium and oxygen neighbors in the surface layer make the largest contributions to the electron correlation. Approximately 97% of the adsorption-induced electron correlation energy part of the adsorption energy is recovered by the method of increments (in our chosen expansion), at the MP2 level.

Using the DFT+U method, i.e., first principles density functional theory calculations with the inclusion of on-site Coulomb interaction, the effects of Pd doping on the O vacancy formation energy (E(vac)) in CeO(2) has been studied. We find that E(vac) is lowered from 3.0 eV in undoped ceria to 0.6 eV in the Pd-doped compound. Much of this decrease can be attributed to emerging Pd-induced gap states above the valence band and below the empty Ce 4f states. These localized defect states involve the Pd ion and its nearest neighbors, which are also the main acceptors of the extra electrons left on reduction. The effect of the Pd dopant on the geometric structure is very modest for CeO(2) but considerable for CeO(2-x).

Ab initio embedded-cluster calculations have been performed for the CeO2(110) surface using temperature induced structures from molecular dynamics (MD) snapshots. As a first step towards understanding how temperature induced distortions of the surface structure influence the surface oxygen reactivity, the energy cost of removing an O atom from the surface was calculated for 41 snapshots from the MD simulation at 300 K. The quantum mechanical embedded-cluster calculations show that already at 300 K the dynamics causes significant fluctuations (root mean square of 0.37 eV) in the O vacancy formation energy (Evac) while the distribution of the two excess electrons associated with the vacancy is virtually unaffected by the surface dynamics and remains localized on the two Ce ions close to the vacancy. It is also found that the quantum mechanical Evac fluctuations can be reproduced by oxygen vacancy calculations using only the relaxed shell-model force field (FF) itself and the MD geometries. Using the FF as the interaction model, the effect of raising the temperature to 750 K and the effect of doping with Ca were investigated for the oxygen vacancy formation.

The Zr doping in CeO(2) may change the reduction properties and therefore the redox properties of CeO(2). Using first-principles density functional theory with the inclusion of on-site Coulomb interaction for a 96-atom supercell, these effects are studied by comparing the differences in atomic structures, electronic structures, and reduction energies of the doped CeO(2) and those of the nondoped CeO(2). It is found that (1) Zr doping of the ceria structure results in important modifications involving nonequivalent O atoms; (2) the oxygen anions (still four-coordinated) next to the doping center show considerably lower reduction energies (by 0.6 eV) and larger displacements ("higher mobilities"); (3) an O vacancy is most easily created close to the Zr centers, therefore the Zr-doping centers might serve as nucleation centers for vacancy clustering; and (4) the electrons left by the released oxygen localize on two Ce cations neighboring the vacancy, which results in the reduction of two Ce(4+) ions.

The optimized geometry, energetics, and vibrational properties of Al(D2O)n3+ clusters, with n = 1,2,4, and 6, have been studied using plane waves, different local basis sets, different methodologies [density-functional theory, MP2, CCSD(T)], and different functionals (BLYP, PBE). Moreover, Car-Parrinello molecular-dynamics (MD) simulations using the BLYP functional, plane waves, and the Vanderbilt ultrasoft pseudopotentials have been performed for an aqueous Al3+ solution with 1 ion and 32 D2O molecules in a periodic box at room temperature, studied for 10 ps. The cluster calculations were performed to pinpoint possible shortcomings of the electronic structure description used in the Car-Parinello MD (CPMD) simulation. For the clusters, the hydration structure and interaction energies calculated with the 'BLYP/plane-wave' approach agree well with high-level ab initio methods but the exchange-correlation functional introduces errors in the OD stretching frequencies (both in the absolute values and in the ion-induced shifts). For the aqueous solution, the CPMD simulation yields structural properties in good agreement with experimental data. The CPMD-simulated OD stretching vibrational band for the first-shell water molecules around Al3+ is strongly downshifted by the influence of the ion and is compared with experimental data from the literature. To make such a comparison meaningful, the influences of a number of systematic effects have been addressed, such as the exchange-correlation functional, the fictitious electron mass, anharmonicity effects, and the small box size in the simulation. Each of these factors (except the last one) is found to affect the OD frequency by 100 cm(-1) or more. The final "corrected" frequencies agree with experiment within approximately 30 cm(-1) for bulk water but are too little downshifted for the first-shell Al3+ (aq) water molecules (by approximately 200 cm(-1)).

The atomic and electronic structure of (111), (110), and (100) surfaces of ceria (CeO2) were studied using density-functional theory within the generalized gradient approximation. Both stoichiometric surfaces and surfaces with oxygen vacancies (unreduced and reduced surfaces, respectively) have been examined. It is found that the (111) surface is the most stable among the considered surfaces, followed by (110) and (100) surfaces, in agreement with experimental observations and previous theoretical results. Different features of relaxation are found for the three surfaces. While the (111) surface undergoes very small relaxation, considerably larger relaxations are found for the (110) and (100) surfaces. The formation of an oxygen vacancy is closely related to the surface structure and occurs more easily for the (110) surface than for (111). The preferred vacancy location is in the surface layer for CeO2(110) and in the subsurface layer (the second O-atomic layer) for CeO2(111). For both surfaces, the O vacancy forms more readily than in the bulk. An interesting oscillatory behavior is found for the vacancy formation energy in the upper three layers of CeO2(111). Analysis of the reduced surfaces suggests that the additional charge resulting from the formation of the oxygen vacancies is localized in the first three layers of the surface. Furthermore, they are not only trapped in the 4f states of cerium., ((c) 2004 American Institute of Physics)

A combined molecular dynamics (MD)+quantum mechanics (QM) method for studying processes on ionic surfaces is presented. Through the combination of classical MD and ab initio embedded-cluster calculations, this method allows the modeling of surface processes involving both the structural and dynamic features of the substrate, even for large-scale systems. The embedding approach used to link the information from the MD simulation to the cluster calculation is presented, and rigorous tests have been carried out to ensure the feasibility of the method. The electrostatic potential and electron density resulting from our embedded-cluster model have been compared with periodic slab results, and confirm the satisfying quality of our embedding scheme as well as the importance of applying embedding in our combined MD+QM approach. We show that a highly accurate representation of the Madelung potential becomes a prerequisite when the embedded-cluster approach is applied to temperature-distorted surface snapshots from the MD simulation.

Many-body potentials for the aqueous Li(+), Na(+), Mg(2+), and Al(3+) ions have been constructed from ab initio cluster calculations. Pure pair, effective pair, effective three-body, and effective polarizable models were created and used in subsequent molecular dynamics simulations. The structures of the first and second solvation shells were studied using radial distribution functions and angular-radial distribution functions. The effective three-body and polarizable potentials yield similar first-shell structures, while the contraction of the O-O distances between the first and second solvation shells is more pronounced with the polarizable potentials. The definition of the tilt angle of the water molecules around the ions is discussed. When a proper definition is used, it is found that for Li(+), Mg(2+), and Al(3+) the water molecules prefer a trigonal orientation, but for Na(+) a tetrahedral orientation (ion in lone-pair direction) is preferred. The self-diffusion coefficients for the water molecules and the ions were calculated; the ionic values follow the order obtained from experiment, although the simulated absolute values are smaller than experiment for Mg(2+) and Al(3+).