Study of the oxygen vacancy influence on magnetic properties of Fe- and Co-doped SnO2 diluted alloys

Transition-metal (TM)-doped diluted magnetic oxides (DMOs) have attracted attention from both experimental and theoretical points of view due to their potential use in spintronics towards new nanostructured devices and new technologies. In the present work, we study the magnetic properties of Sn0.96TM0.04O2 and Sn0.96TM0.04O1.98(V (O))(0.02), where TM = Fe and Co, focusing in particular in the role played by the presence of O vacancies nearby the TM. The calculated total energy as a function of the total magnetic moment per cell shows a magnetic metastability, corresponding to a ground state, respectively, with 2 and 1 mu(B)/cell, for Fe and Co. Two metastable states, with 0 and 4 mu(B)/cell were found for Fe, and a single value, 3 mu(B)/cell, for Co. The spin-crossover energies (E (S)) were calculated. The values are E (S) (0/2) = 107 meV and E (S) (4/2) = 25 meV for Fe. For Co, E (S) (3/1) = 36 meV. By creating O vacancies close to the TM site, we show that the metastablity and E (S) change. For iron, a new state appears, and the state with zero magnetic moment disappears. The ground state is 4 mu(B)/cell instead of 2 mu(B)/cell, and the energy E (S) (2/4) is 30 meV. For cobalt, the ground state is then found with 3 mu(B)/cell and the metastable state with 1 mu(B)/cell. The spin-crossover energy E (S) (1/3) is 21 meV. Our results suggest that these materials may be used in devices for spintronic applications that require different magnetization states.

Changes in a nanoparticle's spectroscopic signal mediated by the local environment

Using a first-principles theoretical model the adsorption of a methyl radical on different sized silver nanoparticles is compared to the adsorption of the same radical on model surfaces. Calculations of our structural, dynamical and electronic properties indicated that small changes in the local environment will lead to small changes in infrared (IR) wavenumbers, but in dramatic changes in the IR signal. Our calculations indicate the lower the adsorption site coordination, the higher is the signal strength, suggesting that small changes in the electronic charge distribution will result in bigger changes in the polarizability and hence in the spectroscopic signal intensity. This effect explains, among others, the signal magnification observed for nanoparticles in surface enhanced Raman spectroscopic (SERS) experiments.

Digital magnetic heterostructures based on GaN using GGA-1/2 approach

We present ab-initio calculations of seven digital magnetic heterostructures, GaN delta-doped with V, Cr, Mn, Fe, Co, Ni, and Cu, forming two-dimensional systems. Only GaN delta-doped with V or Cr present a ferromagnetic ground state with high Curie temperatures. For both, to better describe the electronic properties, we used the GGA-1/2 approach. The ground state of GaN/Cr resulted in a two dimensional half-metal, with 100% spin polarization. For GaN/V, we obtained an insulating state: integer magnetic moment of 2.0 mu(B), a minority spin gap of 3.0 eV close to the gap of GaN, but a majority spin gap of 0.34 eV. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4751285]

The role of Ti as a catalyst for the dissociation of hydrogen on a Mg(0001) surface

In this paper, the dissociative chemisorption of hydrogen on both pure and Ti-incorporated Mg(0001) surfaces are studied by ab initio density functional theory (DFT) calculations. The calculated dissociation barrier of hydrogen molecule on a pure Mg(0001) surface (1.05 eV) is in good agreement with comparable theoretical studies. For the Ti-incorporated Mg(0001) surface, the activated barrier decreases to 0.103 eV due to the strong interaction between the molecular orbital of hydrogen and the d metal state of Ti. This could explain the experimentally observed improvement in absorption kinetics of hydrogen when transition metals have been introduced into the magnesium materials.

Ab initio studies of hydrogen desorption from low index magnesium hydride surface

The low index Magnesium hydride surfaces, MgH2(001) and MgH2(110), have been studied by ab intio Density Functional Theory (DFT) calculations. It was found that the MgH2(110) surface is more stable than MgH2(001) surface, which is in good agreement with the experimental observation. The H-2 desorption barriers vary depending on the crystalline surfaces that are exposed and also the specific H atom sites involved-they are found to be generally high, due to the thermodynamic stability of the MgH2, system, and are larger for the MgH2(001) surface. The pathway for recombinative desorption of one in-plane and one bridging H atom from the MgH2(110) surface was found to be the lowest energy barrier amongst those computed (172 KJ/mol) and is in good agreement with the experimental estimates. (c) 2006 Elsevier B.V. All rights reserved.

Catalytic effects of subsurface carbon in the chemisorption of hydrogen on a Mg(0001) surface: an ab-initio study

Ab initio density functional theory (DFT) calculations are performed to explore possible catalytic effects on the dissociative chemisorption of hydrogen on a Mg(0001) surface when carbon is incorporated into Mg materials. The computational results imply that a C atom located initially on a Mg(0001) surface can migrate into the subsurface and occupy an fcc interstitial site, with charge transfer to the C atom from neighboring Mg atoms. The effect of subsurface C on the dissociation of H-2 on the Mg(0001) surface is found to be relatively marginal: a perfect sublayer of interstitial C is calculated to lower the barrier by 0.16 eV compared with that on a pure Mg(0001) surface. Further calculations reveal, however, that sublayer C may have a significant effect in enhancing the diffusion of atomic hydrogen into the sublayers through fcc channels. This contributes new physical understanding toward rationalizing the experimentally observed improvement in absorption kinetics of H2 when graphite or single walled carbon nanotubes (SWCNT) are introduced into the Mg powder during ball milling.

First-principle study of adsorption of hydrogen on Ti-doped Mg(0001) surface

Ab initio density functional theory (DFT) calculations are performed to study the adsorption of H-2 molecules on a Ti-doped Mg(0001) surface. We find that two hydrogen molecules are able to dissociate on top of the Ti atom with very small activation barriers (0.103 and 0.145 eV for the first and second H-2 molecules, respectively). Additionally, a molecular adsorption state of H-2 above the Ti atom is observed for the first time and is attributed to the polarization of the H-2 molecule by the Ti cation. Our results parallel recent findings for H-2 adsorption on Ti-doped carbon nanotubes or fullerenes. They provide new insight into the preliminary stages of hydrogen adsorption onto Ti-incorporated Mg surfaces.

Low valence cation doping of bulk Cr2O3: Charge compensation and oxygen vacancy formation

The different oxidation states of chromium allow its bulk oxide form to be reducible, facilitating the oxygen vacancy formation process, which is a key property in applications such as catalysis. Similar to other useful oxides such as TiO2, and CeO2, the effect of substitutional metal dopants in bulk Cr2O3 and its effect on the electronic structure and oxygen vacancy formation are of interest, particularly in enhancing the latter. In this paper, density functional theory (DFT) calculations with a Hubbard + U correction (DFT+U) applied to the Cr 3d and O 2p states, are carried out on pure and metal-doped bulk Cr2O3 to examine the effect of doping on the electronic and geometric structure. The role of dopants in enhancing the reducibility of Cr2O3 is examined to promote oxygen vacancy formation. The dopants are Mg, Cu, Ni, and Zn, which have a formal +2 oxidation state in their bulk oxides. Given this difference in host and, dopant oxidation states, we show that to predict the correct ground state two metal dopants charge compensated with an oxygen vacancy are required. The second oxygen atom removed is termed "the active" oxygen vacancy and it is the energy required to remove this atom that is related to the reduction process. In all cases, we find that substitutional doping improves the oxygen vacancy formation of bulk Cr2O3 by lowering the energy cost.

An understanding of chemoselective hydrogenation on crotonaldehyde over Pt(111) in the free energy landscape: The microkinetics study based on first-principles calculations

Microkinetic model is developed in the free energy landscape based on density functional theory (DFT) to quantitatively investigate the reaction mechanism of chemoselective partial hydrogenation of crotonaldehyde to crotyl alcohol over Pt(1 1 1) at the temperature of 353 K. Three different methods (mobile, immobile and collision theory models) were carried out to obtain free energy barrier of adsorption/desorption processes. The results from mobile and collision theory models are similar. The calculated TOFs from both models are close to the experiment value. However, for the immobile model, in which the free energy barrier of desorption approaches the energy barrier, the calculated TOF is 2 orders of magnitude lower than the other models. The difficulty of adsorption/ desorption may be overestimated in the immobile model. In addition, detailed analyses show that for the surface hydrogenation elementary steps, the entropy and internal energy effects are small under the reaction condition, while the zero-point-energy (ZPE) correction is significant, especially for the multi-step hydrogenation reaction. The total energy with the ZPE correction approaches to the full free energy calculation for the surface reaction under the reaction condition. (c) 2011 Elsevier B.V. All rights reserved.

Theoretical and experimental characterisation of energy in an electrostatic discharge

Electrostatic discharges have been identified as the most likely cause in a number of incidents of fire and explosion with unexplained ignitions. The lack of data and suitable models for this ignition mechanism creates a void in the analysis to quantify the importance of static electricity as a credible ignition mechanism. Quantifiable hazard analysis of the risk of ignition by static discharge cannot, therefore, be entirely carried out with our current understanding of this phenomenon. The study of electrostatics has been ongoing for a long time. However, it was not until the wide spread use of electronics that research was developed for the protection of electronics from electrostatic discharges. Current experimental models for electrostatic discharge developed for intrinsic safety with electronics are inadequate for ignition analysis and typically are not supported by theoretical analysis. A preliminary simulation and experiment with low voltage was designed to investigate the characteristics of energy dissipation and provided a basis for a high voltage investigation. It was seen that for a low voltage the discharge energy represents about 10% of the initial capacitive energy available and that the energy dissipation was within 10 ns of the initial discharge. The potential difference is greatest at the initial break down when the largest amount of the energy is dissipated. The discharge pathway is then established and minimal energy is dissipated as energy dissipation becomes greatly influenced by other components and stray resistance in the discharge circuit. From the initial low voltage simulation work, the importance of the energy dissipation and the characteristic of the discharge were determined. After the preliminary low voltage work was completed, a high voltage discharge experiment was designed and fabricated. Voltage and current measurement were recorded on the discharge circuit allowing the discharge characteristic to be recorded and energy dissipation in the discharge circuit calculated. Discharge energy calculations show consistency with the low voltage work relating to discharge energy with about 30-40% of the total initial capacitive energy being discharged in the resulting high voltage arc. After the system was characterised and operation validated, high voltage ignition energy measurements were conducted on a solution of n-Pentane evaporating in a 250 cm3 chamber. A series of ignition experiments were conducted to determine the minimum ignition energy of n-Pentane. The data from the ignition work was analysed with standard statistical regression methods for tests that return binary (yes/no) data and found to be in agreement with recent publications. The research demonstrates that energy dissipation is heavily dependent on the circuit configuration and most especially by the discharge circuit's capacitance and resistance. The analysis established a discharge profile for the discharges studied and validates the application of this methodology for further research into different materials and atmospheres; by systematically looking at discharge profiles of test materials with various parameters (e.g., capacitance, inductance, and resistance). Systematic experiments looking at the discharge characteristics of the spark will also help understand the way energy is dissipated in an electrostatic discharge enabling a better understanding of the ignition characteristics of materials in terms of energy and the dissipation of that energy in an electrostatic discharge.

Energy gap in the aetiology of body weight gain and obesity: A challenging concept with a complex evaluation and pitfalls

The concept of energy gap(s) is useful for understanding the consequence of a small daily, weekly, or monthly positive energy balance and the inconspicuous shift in weight gain ultimately leading to overweight and obesity. Energy gap is a dynamic concept: an initial positive energy gap incurred via an increase in energy intake (or a decrease in physical activity) is not constant, may fade out with time if the initial conditions are maintained, and depends on the 'efficiency' with which the readjustment of the energy imbalance gap occurs with time. The metabolic response to an energy imbalance gap and the magnitude of the energy gap(s) can be estimated by at least two methods, i.e. i) assessment by longitudinal overfeeding studies, imposing (by design) an initial positive energy imbalance gap; ii) retrospective assessment based on epidemiological surveys, whereby the accumulated endogenous energy storage per unit of time is calculated from the change in body weight and body composition. In order to illustrate the difficulty of accurately assessing an energy gap we have used, as an illustrative example, a recent epidemiological study which tracked changes in total energy intake (estimated by gross food availability) and body weight over 3 decades in the US, combined with total energy expenditure prediction from body weight using doubly labelled water data. At the population level, the study attempted to assess the cause of the energy gap purported to be entirely due to increased food intake. Based on an estimate of change in energy intake judged to be more reliable (i.e. in the same study population) and together with calculations of simple energetic indices, our analysis suggests that conclusions about the fundamental causes of obesity development in a population (excess intake vs. low physical activity or both) is clouded by a high level of uncertainty.

Free Energy Barriers for Escape of Water Molecules from Protein Hydration Layer

Free energy barriers separating interfacial water molecules from the hydration layer at the surface of a protein to the bulk are obtained by using the umbrella sampling method of free energy calculation. We consider hydration layer of chicken villin head piece (HP-36) which has been studied extensively by molecular dynamics simulations. The free energy calculations reveal a strong sensitivity to the secondary structure. In particular, we find a region near the junction of first and second helix that contains a cluster of water molecules which are slow in motion, characterized by long residence times (of the order of 100 ps or more) and separated by a large free energy barrier from the bulk water. However, these ``slow'' water molecules constitute only about 5-10% of the total number of hydration layer water molecules. Nevertheless, they play an important role in stabilizing the protein conformation. Water molecules near the third helix (which is the important helix for biological function) are enthalpically least stable and exhibit the fastest dynamics. Interestingly, barrier height distributions of interfacial water are quite broad for water surrounding all the three helices (and the three coils), with the smallest barriers found for those near the helix-3. For the quasi-bound water molecules near the first and second helices, we use well-known Kramers' theory to estimate the residence time from the free energy surface, by estimating the friction along the reaction coordinate from the diffusion coefficient by using Einstein relation. The agreement found is satisfactory. We discuss the possible biological function of these slow, quasi-bound (but transient) water molecules on the surface.

Trends in elasticity and electronic structure of 5d transition metal diborides: first-principles calculations

We investigate the cohesive energy, heat of formation, elastic constant and electronic band structure of transition metal diborides TMB2 (TM = Hf, Ta, W, Re, Os and Ir, Pt) in the Pmmn space group using the ab initio pseudopotential total energy method. Our calculations indicate that there is a relationship between elastic constant and valence electron concentration (VEC): the bulk modulus and shear modulus achieve their maximum when the VEC is in the range of 6.8-7.2. In addition, trends in the elastic constant are well explained in terms of electronic band structure analysis, e.g., occupation of valence electrons in states near the Fermi level, which determines the cohesive energy and elastic properties. The maximum in bulk modulus and shear modulus is attributed to the nearly complete filling of TM d-B p bonding states without filling the antibonding states. On the basis of the observed relationship, we predict that alloying W and Re in the orthorhombic structure OsB2 might be harder than alloying the Ir element. Indeed, the further calculations confirmed this expectation.

Investigation of numerical atomic orbitals for first-principles calculations of the electronic and transport properties of silicon nanowire structures

This thesis is focused on the application of numerical atomic basis sets in studies of the structural, electronic and transport properties of silicon nanowire structures from first-principles within the framework of Density Functional Theory. First we critically examine the applied methodology and then offer predictions regarding the transport properties and realisation of silicon nanowire devices. The performance of numerical atomic orbitals is benchmarked against calculations performed with plane waves basis sets. After establishing the convergence of total energy and electronic structure calculations with increasing basis size we have shown that their quality greatly improves with the optimisation of the contraction for a fixed basis size. The double zeta polarised basis offers a reasonable approximation to study structural and electronic properties and transferability exists between various nanowire structures. This is most important to reduce the computational cost. The impact of basis sets on transport properties in silicon nanowires with oxygen and dopant impurities have also been studied. It is found that whilst transmission features quantitatively converge with increasing contraction there is a weaker dependence on basis set for the mean free path; the double zeta polarised basis offers a good compromise whereas the single zeta basis set yields qualitatively reasonable results. Studying the transport properties of nanowire-based transistor setups with p+-n-p+ and p+-i-p+ doping profiles it is shown that charge self-consistency affects the I-V characteristics more significantly than the basis set choice. It is predicted that such ultrascaled (3 nm length) transistors would show degraded performance due to relatively high source-drain tunnelling currents. Finally, it is shown the hole mobility of Si nanowires nominally doped with boron decreases monotonically with decreasing width at fixed doping density and increasing dopant concentration. Significant mobility variations are identified which can explain experimental observations.

Average ground-state energy of finite Fermi systems

Semiclassical theories such as the Thomas-Fermi and Wigner-Kirkwood methods give a good description of the smooth average part of the total energy of a Fermi gas in some external potential when the chemical potential is varied. However, in systems with a fixed number of particles N, these methods overbind the actual average of the quantum energy as N is varied. We describe a theory that accounts for this effect. Numerical illustrations are discussed for fermions trapped in a harmonic oscillator potential and in a hard-wall cavity, and for self-consistent calculations of atomic nuclei. In the latter case, the influence of deformations on the average behavior of the energy is also considered.