Analysis of a carbon dimer bound to a vacancy in iron using density functional theory and a new tight binding model

Recent density functional theory (DFT) calculations by Foerst et al. have predicted that vacancies in both low and high carbon steels have a carbon dimer bound to them. This is likely to change the thinking of metallurgists in the kinetics of the development of microstructures. While the notion of a C2 molecule bound to a vacancy in Fe will potentially assume a central importance in the atomistic modeling of steels, neither a recent tight binding (TB) model nor existing classical interatomic potentials can account for it. Here we present a new TB model for C in Fe, based on our earlier work for H in Fe, which correctly predicts the structure and energetics of the carbon dimer at a vacancy in Fe. Moreover the model is capable of dealing with both concentrated and dilute limits of carbon in both bcc-Fe and fcc-Fe as comparisons with DFT show. We use both DFT and TB to make a detailed analysis of the dimer and to come to an understanding as to what governs the choice of its curious orientation within the vacancy

A fully quantum mechanical calculation of the diffusivity of hydrogen in iron using the tight binding approximation and path integral theory

We present calculations of free energy barriers and diffusivities as functions of temperature for the diffusion of hydrogen in bcc-Fe. This is a fully quantum mechanical approach since the total energy landscape is computed using a new self consistent, transferable tight binding model for interstitial impurities in magnetic iron. Also the hydrogen nucleus is treated quantum mechanically and we compare here two approaches in the literature both based in the Feynman path integral formulation of statistical mechanics. We find that the quantum transition state theory which admits greater freedom for the proton to explore phase space gives result in better agreement with experiment than the alternative which is based on fixed centroid calculations of the free energy barrier. We also find results in better agreement compared to recent centroid molecular dynamics (CMD) calculations of the diffusivity which employed a classical interatomic potential rather than our quantum mechanical tight binding theory. In particular we find first that quantum effects persist to higher temperatures than previously thought, and conversely that the low temperature diffusivity is smaller than predicted in CMD calculations and larger than predicted by classical transition state theory. This will have impact on future modeling and simulation of hydrogen trapping and diffusion

Vibrational Instability of Metal-Poor Low-Mass Main-Sequence Stars

We find that low-degree low-order g-modes become unstable in metal-poor low-mass stars due to the $\varepsilon$-mechanism of the pp-chain. Since the outer convection zone of these stars is limited only to the very outer layers, the uncertainty in the treatment of convection does not affect the result significantly. The decrease in metallicity leads to decrease in opacity and hence increase in luminosity of a star. This makes the star compact and results in decrease in the density contrast, which is favorable to the $\varepsilon$-mechanism instability. We find also instability for high order g-modes of metal-poor low-mass stars by the convective blocking mechanism. Since the effective temperature and the luminosity of metal-poor stars are significantly higher than those of Pop I stars, the stars showing $\gamma$ Dor-type pulsation are substantially less massive than in the case of Pop I stars. We demonstrate that those modes are unstable for about $1\,M_\odot$ stars in the metal-poor case.Comment: 4 pages, 4 figures, To be published in Astrophysics and Space Science Proceedings series (ASSP). Proceedings of the "20th Stellar Pulsation Conference Series: Impact of new instrumentation and new insights in stellar pulsations", 5-9 September 2011, Granada, Spai

Electronic structure and total energy of interstitial hydrogen in iron: Tight binding models

An application of the tight binding approximation is presented for the description of electronic structure and interatomic force in magnetic iron, both pure and containing hydrogen impurities. We assess the simple canonical d-band description in comparison to a non orthogonal model including s and d bands. The transferability of our models is tested against known properties including the segregation energies of hydrogen to vacancies and to surfaces of iron. In many cases agreement is remarkably good, opening up the way to quantum mechanical atomistic simulation of the effects of hydrogen on mechanical properties

Local volume effects in the generalized pseudopotential theory

The generalized pseudopotential theory (GPT) is a powerful method for deriving real-space transferable interatomic potentials. Using a coarse-grained electronic structure, one can explicitly calculate the pair ion-ion and multi-ion interactions in simple and transition metals. While successful in determining bulk properties, in central force metals the GPT fails to describe crystal defects for which there is a significant local volume change. A previous paper [J. A. Moriarty and R. Phillips, Phys. Rev. Lett. 66, 3036 (1991)PRLTAO0031-900710.1103/PhysRevLett.66.3036] found that by allowing the GPT total energy to depend upon some spatially averaged local electron density, the energetics of vacancies and surfaces could be calculated within experimental ranges. In this paper, we develop the formalism further by explicitly calculating the forces and stress tensor associated with this total energy. We call this scheme the adaptive GPT (aGPT) and it is capable of both molecular dynamics (MD) and molecular statics. We apply the aGPT to vacancy formation, divacancy binding, and stacking faults in hcp Mg. We also calculate the local electron density corrections to the bulk elastic constants and phonon dispersion for which there is refinement over the baseline GPT treatment. In addition, we demonstrate aGPT-MD simulation through the calculation of thermal expansion in magnesium to 700 K.</p

Phoretic Motion of Spheroidal Particles Due To Self-Generated Solute Gradients

We study theoretically the phoretic motion of a spheroidal particle, which generates solute gradients in the surrounding unbounded solvent via chemical reactions active on its surface in a cap-like region centered at one of the poles of the particle. We derive, within the constraints of the mapping to classical diffusio-phoresis, an analytical expression for the phoretic velocity of such an object. This allows us to analyze in detail the dependence of the velocity on the aspect ratio of the polar and the equatorial diameters of the particle and on the fraction of the particle surface contributing to the chemical reaction. The particular cases of a sphere and of an approximation for a needle-like particle, which are the most common shapes employed in experimental realizations of such self-propelled objects, are obtained from the general solution in the limits that the aspect ratio approaches one or becomes very large, respectively.Comment: 18 pages, 5 figures, to appear in European Physical Journal

A tight binding model for water

We demonstrate for the first time a tight binding model for water incorporating polarizable anions. A novel aspect is that we adopt a "ground up" approach in that properties of the monomer and dimer only are fitted. Subsequently we make predictions of the structure and properties of hexamer clusters, ice-XI and liquid water. A particular feature, missing in current tight binding and semiempirical hamiltonians, is that we reproduce the almost two-fold increase in molecular dipole moment as clusters are built up towards the limit of bulk liquid. We concentrate on properties of liquid water which are very well rendered in comparison with experiment and published density functional calculations. Finally we comment on the question of the contrasting densities of water and ice which is central to an understanding of the subtleties of the hydrogen bond

Generalized stacking fault energetics and dislocation properties: compact vs. spread unit dislocation structures in TiAl and CuAu

We present a general scheme for analyzing the structure and mobility of dislocations based on solutions of the Peierls-Nabarro model with a two component displacement field and restoring forces determined from the ab-initio generalized stacking fault energetics (ie., the so-called $\gamma$-surface). The approach is used to investigate dislocations in L1$_{0}$ TiAl and CuAu; predicted differences in the unit dislocation properties are explicitly related with features of the $\gamma$-surface geometry. A unified description of compact, spread and split dislocation cores is provided with an important characteristic "dissociation path" revealed by this highly tractable scheme.Comment: 7 two columns pages, 2 eps figures. Phys. Rev. B. accepted November 199

The stabilizing role of itinerant ferromagnetism in inter-granular cohesion in iron

We present a simple, general energy functional for ferromagnetic materials based upon a local spin density extension to the Stoner theory of itinerant ferromagnetism. The functional reproduces well available ab initio results and experimental interfacial energies for grain boundaries in iron. The model shows that inter-granular cohesion along symmetric tilt boundaries in iron is dependent upon strong magnetic structure at the interface, illuminates the mechanisms underlying this structure, and provides a simple explanation for relaxation of the atomic structure at these boundaries.Comment: In review at Phys. Rev. Lett. Submitted 23 September 1997; revised 16 March 199