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.

Phonon spectra and thermodynamic properties of rare gas solids based on empirical and semi-empirical (ab initio) two-body potentials : a comparative study /

We study the phonon dispersion, cohesive and thermal properties of raxe gas solids Ne, Ar, Kr, and Xe, using a variety of potentials obtained from different approaches; such as, fitting to crystal properties, purely ab initio calculations for molecules and dimers or ab initio calculations for solid crystalline phase, a combination of ab initio calculations and fitting to either gas phase data or sohd state properties. We explore whether potentials derived with a certain approaxih have any obvious benefit over the others in reproducing the solid state properties. In particular, we study phonon dispersion, isothermal ajid adiabatic bulk moduli, thermal expansion, and elastic (shear) constants as a function of temperatiue. Anharmonic effects on thermal expansion, specific heat, and bulk moduli have been studied using A^ perturbation theory in the high temperature limit using the neaxest-neighbor central force (nncf) model as developed by Shukla and MacDonald [4]. In our study, we find that potentials based on fitting to the crystal properties have some advantage, particularly for Kr and Xe, in terms of reproducing the thermodynamic properties over an extended range of temperatiures, but agreement with the phonon frequencies with the measured values is not guaranteed. For the lighter element Ne, the LJ potential which is based on fitting to the gas phase data produces best results for the thermodynamic properties; however, the Eggenberger potential for Ne, where the potential is based on combining ab initio quantum chemical calculations and molecular dynamics simulations, produces results that have better agreement with the measured dispersion, and elastic (shear) values. For At, the Morse-type potential, which is based on M0ller-Plesset perturbation theory to fourth order (MP4) ab initio calculations, yields the best results for the thermodynamic properties, elastic (shear) constants, and the phonon dispersion curves.

Phonon spectra and thermal properties of some fcc metals using embedded-atom potentials /

By employing the embedded-atom potentials of Mei et ai.[l], we have calculated the dynamical matrices and phonon dispersion curves for six fee metals (Cu,Ag,Au,Ni,Pd and Pt). We have also investigated, within the quasiharmonic approximation, some other thermal properties of these metals which depend on the phonon density of states, such as the temperature dependence of lattice constant, coefficient of linear thermal expansion, isothermal and adiabatic bulk moduli, heat capacities at constant volume and constant pressure, Griineisen parameter and Debye temperature. The computed results are compared with the experimental findings wherever possible. The comparison shows a generally good agreement between the theoretical values and experimental data for all properties except the discrepancies of phonon frequencies and Debye temperature for Pd, Pt and Au. Further, we modify the parameters of this model for Pd and Pt and obtain the phonon dispersion curves which is in good agreement with experimental data.

Application of [Lambda] to the fourth perturbation theory in calculating the equation of state of rare gas solids and fcc metals

We have calculated the thermodynamic properties of monatomic fcc crystals from the high temperature limit of the Helmholtz free energy. This equation of state included the static and vibrational energy components. The latter contribution was calculated to order A4 of perturbation theory, for a range of crystal volumes, in which a nearest neighbour central force model was used. We have calculated the lattice constant, the coefficient of volume expansion, the specific heat at constant volume and at constant pressure, the adiabatic and the isothermal bulk modulus, and the Gruneisen parameter, for two of the rare gas solids, Xe and Kr, and for the fcc metals Cu, Ag, Au, Al, and Pb. The LennardJones and the Morse potential were each used to represent the atomic interactions for the rare gas solids, and only the Morse potential was used for the fcc metals. The thermodynamic properties obtained from the A4 equation of state with the Lennard-Jones potential, seem to be in reasonable agreement with experiment for temperatures up to about threequarters of the melting temperature. However, for the higher temperatures, the results are less than satisfactory. For Xe and Kr, the thermodynamic properties calculated from the A2 equation of state with the Morse potential, are qualitatively similar to the A 2 results obtained with the Lennard-Jones potential, however, the properties obtained from the A4 equation of state are in good agreement with experiment, since the contribution from the A4 terms seem to be small. The lattice contribution to the thermal properties of the fcc metals was calculated from the A4 equation of state, and these results produced a slight improvement over the properties calculated from the A2 equation of state. In order to compare the calculated specific heats and bulk moduli results with experiment~ the electronic contribution to thermal properties was taken into account~ by using the free electron model. We found that the results varied significantly with the value chosen for the number of free electrons per atom.

Thermodynamic properties and Debye-Waller factor of fee materials

We have calculated the equation of state and the various thermodynamic properties of monatomic fcc crystals by minimizing the Helmholtz free energy derived in the high temperature limit for the quasiharmonic theory, QH, and the lowest-order (cubic and quartic), 'A2, anharmonic terms of the perturbation theory, PT. The total energy in each case is obtained by adding the static energy. The calculation of the thermal properties was carried out for a nearest-neighbour central-force model of the fcc lattice by means of the appropriate thermodynamic relations. We have calculated the lattice constant, the thermal expansion, the coefficient of volume expansion, the specific heat at constant volume and at constant pressure, the isothermal and adiabatic bulk moduli, and the Griineisen parameter, for the rare-gas solids Kr and Xe, and gold. Morse potential and modified Morse potential were each used to represent the atomic interaction for the three fcc materials. For most of the calculated thermodynamic properties from the QH theory, the results for Kr and Xe with the modified Morse potential show an improvement over the results for the Morse potential when compared with the experimental data. However, the results of the 'A 2 equation of state with the modified Morse potential are in good agreement with experiment only in the case of the specific heat at constant volume and at constant pressure. For Au we have calculated the lattice contribution from the QH and 'A 2 PT and the electronic contribution to the thermal properties. The electronic contribution was taken into account by using the free electron model. The results of the thermodynamic properties calculated with the modified Morse potential were similar to those obtained with the Morse potential. U sing the minimized equation of state we also calculated the Mossbauer recoilless fraction for Kr and Xe and the Debye-Waller factor (DWF) for Pb, AI, eu, Ag, and Au. The Mossbauer recoilless fraction was obtained for the above two potentials and Lennard-Jones potential. The L-J potential gives the best agreement with experiment for Kr. No experimental data exists for Xe. At low temperature the calculated DWF results for Pb, AI, and eu show a good agreement with experimental values, but at high temperature the experimental DWF results increase very rapidly. For Ag the computed values were below the expected results at all temperatures. The DWF results of the modified Morse potential for Pb, AI, eu and Ag were slightly better than those of the Morse potential. In the case of Au the calculated values were in poor agreement with experimental results. We have calculated the quasiharmonic phonon dispersion curves for Kr, Xe, eu, Ag, and Au. The calculated and experimental results of the frequencies agree quite well for all the materials except for Au where the longitudinal modes show serious discrepancies with the experimental results. In addition, the two lowest-order anharmonic contributions to the phonon frequency were derived using the Green's function method. The A 2 phonon dispersion curves have been calculated only for eu, and the results were similar to those of the QH dispersion curves. Finally, an expression for the Griineisen parameter "( has been derived from the anharmonic frequencies, and calculated for these materials. The "( results are comparable with those obtained from the thermodynamic definition.

Phonon spectra and temperature variation of thermodynamic properties of fcc metals via Finnis-Sinclair type many body potentials: Sutton-Chen and improved Sutton-Chen models

Volume(density)-independent pair-potentials cannot describe metallic cohesion adequately as the presence of the free electron gas renders the total energy strongly dependent on the electron density. The embedded atom method (EAM) addresses this issue by replacing part of the total energy with an explicitly density-dependent term called the embedding function. Finnis and Sinclair proposed a model where the embedding function is taken to be proportional to the square root of the electron density. Models of this type are known as Finnis-Sinclair many body potentials. In this work we study a particular parametrization of the Finnis-Sinclair type potential, called the "Sutton-Chen" model, and a later version, called the "Quantum Sutton-Chen" model, to study the phonon spectra and the temperature variation thermodynamic properties of fcc metals. Both models give poor results for thermal expansion, which can be traced to rapid softening of transverse phonon frequencies with increasing lattice parameter. We identify the power law decay of the electron density with distance assumed by the model as the main cause of this behaviour and show that an exponentially decaying form of charge density improves the results significantly. Results for Sutton-Chen and our improved version of Sutton-Chen models are compared for four fcc metals: Cu, Ag, Au and Pt. The calculated properties are the phonon spectra, thermal expansion coefficient, isobaric heat capacity, adiabatic and isothermal bulk moduli, atomic root-mean-square displacement and Gr\"{u}neisen parameter. For the sake of comparison we have also considered two other models where the distance-dependence of the charge density is an exponential multiplied by polynomials. None of these models exhibits the instability against thermal expansion (premature melting) as shown by the Sutton-Chen model. We also present results obtained via pure pair potential models, in order to identify advantages and disadvantages of methods used to obtain the parameters of these potentials.

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.

LDA or GGA? A combined experimental inelastic neutron scattering and ab initio lattice dynamics study of alkali metal hydrides

In a previous work, we carried out inelastic neutron scattering (INS) spectroscopy experiments and preliminary first principles calculations on alkali metal hydrides. The complete series of alkali metal hydrides, LiH, NaH, KH, RbH and CsH was measured in the high-resolution TOSCA INS spectrometer at ISIS. Here, we present the results of ab initio electronic structure calculations of the properties of the alkali metal hydrides using both the local density approximation (LDA) and the generalized gradient approximation (GGA), using the Perdew–Burke–Ernzerhof (PBE) parameterization. Properties calculated were lattice parameters, bulk moduli, dielectric constants, effective charges, electronic densities and inelastic neutron scattering (INS) spectra. We took advantage of the currently available computer power to use full lattice dynamics theory to calculate thermodynamic properties for these materials. For the alkali metal hydrides (LiH, NaH, KH, RbH and CsH) using lattice dynamics, we found that the INS spectra calculated using LDA agreed better with the experimental data than the spectra calculated using GGA. Both zero-point effects and thermal contributions to free energies had an important effect on INS and several thermodynamic properties.

Crystal engineering using functionalized adamantane

We performed a first-principles investigation on the structural, electronic and optical properties of crystals made of chemically functionalized adamantane molecules. Several molecular building blocks, formed by boron and nitrogen substitutional functionalizations, were considered to build zinc blende and wurtzite crystals, and the resulting structures presented large bulk moduli and cohesive energies, wide and direct bandgaps, and low dielectric constants (low-kappa materials). Those properties provide stability for such structures up to room temperature, superior to those of typical molecular crystals. This indicates a possible road map for crystal engineering using functionalized diamondoids, with potential applications ranging from space filling between conducting wires in nanodevices to nano-electromechanical systems.

Boron and nitrogen functionalized diamondoids: A first principles investigation

Chemically functionalized adamantane molecules have been investigated by first principles total energy calculations. Boron and nitrogen functionalized molecules were found to be very stable, consistent with available experimental data. Two hypothetical molecular crystals, involving functionalized adamantane, were investigated. These molecular crystals presented direct electronic bandgaps and large bulk moduli, which suggested a possible road for molecular self-assembly using functionalized diamondoids. (C) 2010 Elsevier B.V. All rights reserved.

Spin-strain coupling in a 3-D transition metal oxides

ab-initio Hartree Fock (HF), density functional theory (DFT) and hybrid potentials were employed to compute the optimized lattice parameters and elastic properties of perovskite 3-d transition metal oxides. The optimized lattice parameters and elastic properties are interdependent in these materials. An interaction is observed between the electronic charge, spin and lattice degrees of freedom in 3-d transition metal oxides. The coupling between the electronic charge, spin and lattice structures originates due to localization of d-atomic orbitals. The coupling between the electronic charge, spin and crystalline lattice also contributes in the ferroelectric and ferromagnetic properties in perovskites. The cubic and tetragonal crystalline structures of perovskite transition metal oxides of ABO3 are studied. The electronic structure and the physics of 3-d perovskite materials is complex and less well considered. Moreover, the novelty of the electronic structure and properties of these perovskites transition metal oxides exceeds the challenge offered by their complex crystalline structures. To achieve the objective of understanding the structure and property relationship of these materials the first-principle computational method is employed. CRYSTAL09 code is employed for computing crystalline structure, elastic, ferromagnetic and other electronic properties. Second-order elastic constants (SOEC) and bulk moduli (B) are computed in an automated process by employing ELASTCON (elastic constants) and EOS (equation of state) programs in CRYSTAL09 code. ELASTCON, EOS and other computational algorithms are utilized to determine the elastic properties of tetragonal BaTiO3, rutile TiO2, cubic and tetragonal BaFeO3 and the ferromagentic properties of 3-d transition metal oxides. Multiple methods are employed to crosscheck the consistency of our computational results. Computational results have motivated us to explore the ferromagnetic properties of 3-d transition metal oxides. Billyscript and CRYSTAL09 code are employed to compute the optimized geometry of the cubic and tetragonal crystalline structure of transition metal oxides of Sc to Cu. Cubic crystalline structure is initially chosen to determine the effect of lattice strains on ferromagnetism due to the spin angular momentum of an electron. The 3-d transition metals and their oxides are challenging as the basis functions and potentials are not fully developed to address the complex physics of the transition metals. Moreover, perovskite crystalline structures are extremely challenging with respect to the quality of computations as the latter requires the well established methods. Ferroelectric and ferromagnetic properties of bulk, surfaces and interfaces are explored by employing CRYSTAL09 code. In our computations done on cubic TMOs of Sc-Fe it is observed that there is a coupling between the crystalline structure and FM/AFM spin polarization. Strained crystalline structures of 3-d transition metal oxides are subjected to changes in the electromagnetic and electronic properties. The electronic structure and properties of bulk, composites, surfaces of 3-d transition metal oxides are computed successfully.

High-Pressure Behavior and Phase Stability of Al5BO9, a Mullite-Type Ceramic Material

Phase stability, elastic behavior, and pressure-induced structural evolution of synthetic boron-mullite Al5BO9 (a = 5.6780(7), b = 15.035(6), and c =7.698(3) Å, space group Cmc21, Z = 4) were investigated up to 25.6(1) GPa by in situ single-crystal synchrotron X-ray diffraction with a diamond anvil cell (DAC) under hydrostatic conditions. No evidence of phase transition was observed up to 21.7(1) GPa. At 25.6(1) GPa, the refined unit-cell parameters deviated significantly from the compressional trend, and the diffraction peaks appeared broader than at lower pressure. At 26.7(1) GPa, the diffraction pattern was not indexable, suggesting amorphization of the material or a phase transition to a high-pressure polymorph. Fitting the P–V data up to 21.7(1) GPa with a second-order Birch–Murnaghan Equation-of-State, we obtained a bulk modulus KT0 = 164(1) GPa. The axial compressibilities, here described as linearized bulk moduli, are as follows: KT0(a) = 244(9), KT0(b) = 120(4), and KT0(c) = 166(11) GPa (KT0(a):KT0(b):KT0(c) = 2.03:1:1.38). The structure refinements allowed a description of the main deformation mechanisms in response to the applied pressure. The stiffer crystallographic direction appears to be controlled by the infinite chains of edge-sharing octahedra running along [100], making the structure less compressible along the a-axis than along the b- and c-axis.

Synthesis of M2AC (M = Ti, V, Cr, Nb, Zr, A = Al, Sn, S) (MAX) carbide phases and study of their physical behavior under extreme conditions

Materials known as Mn+1AXn phases, where n is 1, 2, or 3, and M represents an early transition metal, A an A-group element, and X is either Carbon and/or Nitrogen [1], are fast becoming technologically important materials due to the interesting combination of unique properties. However, a lot of important information about the high temperature and high pressure behavior of many of these compounds is still missing, which needs to be determined systematically. ^ In this dissertation the synthesis of M2AC (M = Ti, V, Cr, Nb, Zr) and A = (Al, Sn, S) compounds by arc melting, vacuum sintering and piston cylinder synthesis is presented along with the synthesis of Zr 2SC, which has been synthesized for first time in bulk form, by piston cylinder technique. The microstructural analysis by electron microscopy and phase analysis by x-ray diffraction is presented next. Finally, a critical analysis of the behavior of these compounds under the application of extreme pressure (as high as 50 GPa) and temperature (≈ 1000°C) is presented. ^ The high pressure studies, up to 50 GPa, showed that these compounds were structurally intact and their bulk moduli ranged from 140 to 190 GPa. The high temperature studies in the inert atmosphere showed that the M 2SnC compounds were unstable above 650°C and the expansion along the a-axis was higher than that along the c-axis, unlike the other phases. M2SC compounds on the other hand showed negligible difference in the thermal expansion along the two axes. The oxidation study revealed that Ti2AC (Al, S) compounds had highest resistance to oxidation while the M2SnC compounds had the least. Furthermore, from the oxidation study of these compounds, which were short time oxidation experiments, it was found that all of these compounds oxidized to their respective binary oxides. ^

Applications of synchrotron x-ray powder diffraction in hydrated cements: high-resolution and high-pressure studies

The main aim of this study is to apply synchrotron radiation techniques for the study of hydrated cement pastes. In particular, the tetracalcium aluminoferrite phase, C4AF in cement nomenclature, is the major iron-containing phase in Ordinary Portland Cement (OPC) and in iron rich belite calcium sulfoaluminate cements. In a first study, the hydration mechanism of pure tetracalcium aluminoferrite phase with water-to-solid ratio of 1.0 has been investigated by HR-SXRPD (high resolution synchrotron X-ray powder diffraction). C4AF in the presence of water hydrates to form mainly an iron-containing hydrogarnet-type (katoite) phase, C3A0.84F0.16H6, as single crystalline phase. Its crystal structure and stoichiometry were determined by the Rietveld method and the final disagreement factors were RWP=8.1% and RF=4.8% [1]. As the iron content in the product is lower than that in C4AF, it is assumed that part of the iron also goes to an amorphous iron rich gel, like the hydrated alumina-type gel, as hydration proceeds. Further results from the high-resolution study will be discussed. In a second study, the behavior of pure and iron-containing katoites (C3AH6 and C3A0.84F0.16H6) under pressure have been analyzed by SXRPD using a diamond anvil cell (DAC) and then their bulk moduli were determined. The role of the pressure transmitting medium (PTM) has also been studied. In this case, silicone oil as well as methanol/ethanol mixtures have been used as PTM. Some “new peaks” were detected in the pattern for C3A0.84F0.16H6 as pressure increases, when using ethanol/methanol as PTM. These new peaks were still present at ambient pressure after releasing the applied pressure. They may correspond to crystalline nordstrandite or doyleite from the crystallization of amorphous aluminium hydroxide. The results from the high-pressure study will also be discussed.