First principles LDA+U and GGA+U study of cerium oxides: Dependence on the effective U-parameter

The electronic structure and properties of cerium oxides (CeO2 and Ce2O3) have been studied in the framework of the LDA+U and GGA(PW91)+U implementations of density functional theory. The dependence of selected observables of these materials on the effective U parameter has been investigated in detail. The examined properties include lattice constants, bulk moduli, density of states, and formation energies of CeO2 and Ce2O3. For CeO2, the LDA+U results are in better agreement with experiment than the GGA+U results whereas for the computationally more demanding Ce2O3 both approaches give comparable accuracy. Furthermore, as expected, Ce2O3 is much more sensitive to the choice of the U value. Generally, the PW91 functional provides an optimal agreement with experiment at lower U energies than LDA does. In order to achieve a balanced description of both kinds of materials, and also of nonstoichiometric CeO2¿x phases, an appropriate choice of U is suggested for LDA+U and GGA+U schemes. Nevertheless, an optimum value appears to be property dependent, especially for Ce2O3. Optimum U values are found to be, in general, larger than values determined previously in a self-consistent way.

Theoretical approach to the magnetostructural correlations in spin-Peierls compound CuGeO3

A theoretical density-functional study has been carried out to analyze the exchange coupling in the chains of CuGeO3 using discrete models. The results show a good agreement with the experimental exchange coupling constant (J) together with a strong dependence of J with the Cu-O-Cu angle. The calculation of the J values for a distorted model indicates a larger degree of dimerization than those reported previously.

Computational Approach to the Study of Thermal Spin Crossover Phenomena

The key parameters associated to the thermally induced spin crossover process have been calculated for a series of Fe(II) complexes with mono-, bi-, and tridentate ligands. Combination of density functional theory calculations for the geometries and for normal vibrational modes, and highly correlated wave function methods for the energies, allows us to accurately compute the entropy variation associated to the spin transition and the zero-point corrected energy difference between the low- and high-spin states. From these values, the transition temperature, T 1/2, is estimated for different compounds.

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.

Correlation of the EPR properties of perchlorotriphenylmethyl radical and their efficiency as DNP polarizers

Water soluble perchlorinated trityl (PTM) radicals were found to be effective 95 GHz DNP (dynamic nuclear polarization) polarizers in ex situ (dissolution) 13C DNP (Gabellieri et al., Angew Chem., Int. Ed. 2010, 49, 3360). The degree of the nuclear polarization obtained was reported to be dependent on the position of the chlorine substituents on the trityl skeleton. In addition, on the basis of the DNP frequency sweeps it was suggested that the 13C NMR signal enhancement is mediated by the Cl nuclei. To understand the DNP mechanism of the PTM radicals we have explored the 95 GHz EPR characteristics of these radicals that are relevant to their performance as DNP polarizers. The EPR spectra of the radicals revealed axially symmetric g-tensors. A comparison of the spectra with the 13C DNP frequency sweeps showed that although the solid effect mechanism is operational the DNP frequency sweeps reveal some extra width suggesting that contributions from EPR forbidden transitions involving 35,37Cl nuclear flips are likely. This was substantiated experimentally by ELDOR (electron-electron double resonance) detected NMR measurements, which map the EPR forbidden transitions, and ELDOR experiments that follow the depolarization of the electron spin upon irradiation of the forbidden EPR transitions. DFT (density functional theory) calculations helped to assign the observed transitions and provided the relevant spin Hamiltonian parameters. These results show that the 35,37Cl hyperfine and nuclear quadrupolar interactions cause a considerable nuclear state mixing at 95 GHz thus facilitating the polarization of the Cl nuclei upon microwave irradiation. Overlap of Cl nuclear frequencies and the 13C Larmor frequency further facilitates the polarization of the 13C nuclei by spin diffusion. Calculation of the 13C DNP frequency sweep based on the Cl nuclear polarization showed that it does lead to an increase in the width of the spectra, improving the agreement with the experimental sweeps, thus supporting the existence of a new heteronuclear assisted DNP mechanism.

Effect of pressure on La2(WO4)3 with a modulated scheelite-type structure

We have studied the effect of pressure on the structural and vibrational properties of lanthanum tritungstate La2(WO4)3. This compound crystallizes under ambient conditions in the modulated scheelite-type structure known as the α phase. We have performed x-ray diffraction and Raman scattering measurements up to a pressure of 20 GPa, as well as ab initio calculations within the framework of the density functional theory. Up to 5 GPa, the three methods provide a similar picture of the evolution under pressure of α-La2(WO4)3. At 5 GPa, we begin to observe some structural changes, and above 6 GPa we find that the x-ray patterns cannot be indexed as a single phase. However, we find that a mixture of two phases with C2/c symmetry accounts for all diffraction peaks. Our ab initio study confirms the existence of several C2/c structures, which are very close in energy in this compression range. According to our measurements, a state with medium-range order appears at pressures above 9 and 11 GPa, from x-ray diffraction and Raman experiments, respectively. Based upon our theoretical calculations we propose several high-pressure candidates with high cationic coordinations at these pressures. The compound evolves into a partially amorphous phase at pressures above 20 GPa.

Adsorption energy and spin state of first-row transition metals adsorbed on MgO(100)

Slab and cluster model spin-polarized calculations have been carried out to study various properties of isolated first-row transition metal atoms adsorbed on the anionic sites of the regular MgO(100) surface. The calculated adsorption energies follow the trend of the metal cohesive energies, indicating that the changes in the metal-support and metal-metal interactions along the series are dominated by atomic properties. In all cases, except for Ni at the generalized gradient approximation level, the number of unpaired electron is maintained as in the isolated metal atom. The energy required to change the atomic state from high to low spin has been computed using the PW91 and B3LYP density-functional-theory-based methods. PW91 fails to predict the proper ground state of V and Ni, but the results for the isolated and adsorbed atom are consistent within the method. B3LYP properly predicts the ground state of all first-row transition atom the high- to low-spin transition considered is comparable to experiment. In all cases, the interaction with the surface results in a reduced high- to low-spin transition energy.

First principles analysis of the stability and diffusion of oxygen vacancies in meetal oxides

Oxygen vacancies in metal oxides are known to determine their chemistry and physics. The properties of neutral oxygen vacancies in metal oxides of increasing complexity (MgO, CaO, alpha-Al2O3, and ZnO) have been studied using density functional theory. Vacancy formation energies, vacancy-vacancy interaction, and the barriers for vacancy migration are determined and rationalized in terms of the ionicity, the Madelung potential, and lattice relaxation. It is found that the Madelung potential controls the oxygen vacancy properties of highly ionic oxides whereas a more complex picture arises for covalent ZnO.

Adsoption-induced restructuring of gold nanochains

The chemical properties of single-atomic chains of gold atoms are investigated using density functional calculations. The nanochains are shown to be unusually chemically active with strong chemisorption of oxygen atoms and carbon monoxide. The chemisorption energies vary significantly with the strain/stress conditions for the chain. Oxygen atoms are found to energetically prefer to get incorporated into a chain forming a new type of gold-oxygen nanochain with a conductance of one quantum unit. We suggest that the long bond lengths observed in electron microscopy investigations of gold chains can be due to oxygen incorporation.

Theoretical modeling of photon and electron-stimulated Na and K desorption from SiO2

The mechanism of generation of atomic Na and K from SiO2 samples has been studied using explicitly correlated wave function and density functional theory cluster calculations. Possible pathways for the photon and electron stimulated desorption of Na and K atoms from silicates are proposed, thus providing new insight on the generation of the tenuous Na and K atmosphere of the Moon.

MgO/Ag(001) interface structure and STM images from first principles

The interface of MgO/Ag(001) has been studied with density functional theory applied to slabs. We have found that regular MgO films show a small adhesion to the silver substrate, the binding can be increased in off-stoichiometric regimes, either by the presence of O vacancies at the oxide film or by a small excess of O atoms at the interface between the ceramic to the metal. By means of theoretical methods, the scanning tunneling microscopy signatures of these films is also analyzed in some detail. For defect free deposits containing 1 or 2 ML and at low voltages, tunnelling takes place from the surface Ag substrate, and at large positive voltages Mg atoms are imaged. If defects, oxygen vacancies, are present on the surface of the oxide they introduce much easier channels for tunnelling resulting in big protrusions and controlling the shape of the image, the extra O stored at the interface can also be detected for very thin films.

Oxygen vacancies as active sites for water dissociation on rutile TiO2(110)

Through an interplay between scanning tunneling microscopy experiments and density functional theory calculations, we determine unambiguously the active surface site responsible for the dissociation of water molecules adsorbed on rutile TiO2(110). Oxygen vacancies in the surface layer are shown to dissociate H2O through the transfer of one proton to a nearby oxygen atom, forming two hydroxyl groups for every vacancy. The amount of water dissociation is limited by the density of oxygen vacancies present on the clean surface exclusively. The dissociation process sets in as soon as molecular water is able to diffuse to the active site.

Bonding of gold nanoclusters to oxygen vacancies on rutile TiO2(110)

Through an interplay between scanning tunneling microscopy (STM) and density functional theory (DFT) calculations, we show that bridging oxygen vacancies are the active nucleation sites for Au clusters on the rutile TiO2(110) surface. We find that a direct correlation exists between a decrease in density of vacancies and the amount of Au deposited. From the DFT calculations we find that the oxygen vacancy is indeed the strongest Au binding site. We show both experimentally and theoretically that a single oxygen vacancy can bind 3 Au atoms on average. In view of the presented results, a new growth model for the TiO2(110) system involving vacancy-cluster complex diffusion is presented.

When Langmuir is too simple: H2 dissociation on Pd(111) at high coverage

Recent experiments of H2 adsorption on Pd(111) [T. Mitsui et al., Nature (London) 422, 705 (2003)] have questioned the classical Langmuir picture of second order adsorption kinetics at high surface coverage requiring pairs of empty sites for the dissociative chemisorption. Experiments find that at least three empty sites are needed. Through density functional theory, we find that H2 dissociation is favored on ensembles of sites that involve a Pd atom with no direct interaction with adsorbed hydrogen. Such active sites are formed by aggregation of at least 3 H-free sites revealing the complex structure of the "active sites."