Atomistic modeling of thermodynamic equilibrium and polymorphism of iron

We develop two new modified embedded-atom method (MEAM) potentials for elemental iron, intended to reproduce the experimental phase stability with respect to both temperature and pressure. These simple interatomic potentials are fitted to a wide variety of material properties of bcc iron in close agreement with experiments. Numerous defect properties of bcc iron and bulk properties of the two close-packed structures calculated with these models are in reasonable agreement with the available first-principles calculations and experiments. Performance at finite temperatures of these models has also been examined using Monte Carlo simulations. We attempt to reproduce the experimental iron polymorphism at finite temperature by means of free energy computations, similar to the procedure previously pursued by Müller et al (2007 J. Phys.: Condens. Matter 19 326220), and re-examine the adequacy of the conclusion drawn in the study by addressing two critical aspects missing in their analysis: (i) the stability of the hcp structure relative to the bcc and fcc structures and (ii) the compatibility between the temperature and pressure dependences of the phase stability. Using two MEAM potentials, we are able to represent all of the observed structural phase transitions in iron. We discuss that the correct reproductions of the phase stability among three crystal structures of iron with respect to both temperature and pressure are incompatible with each other due to the lack of magnetic effects in this class of empirical interatomic potential models. The MEAM potentials developed in this study correctly predict, in the bcc structure, the self-interstitial in the (110) orientation to be the most stable configuration, and the screw dislocation to have a non-degenerate core structure, in contrast to many embedded-atom method potentials for bcc iron in the literature.     

Sterically Active Electron Pairs in Lead Sulfide? An Investigation of the Electronic and Vibrational Properties of PbS in the Transition Region Between the Rock Salt and the α-GeTe-Type Modifications

Recently, we have investigated the energy landscape of PbS for many different pressures on the ab initio level by using Hartree-Fock and density functional theory to globally search for possible thermodynamically stable and metastable structures. The perhaps most fascinating observation was that besides the experimentally known modification exhibiting the rock salt structure a second minimum exists close-by on the landscape showing the low-temperature α-GeTe-type structure. In the present study, we investigate the possible reasons for the existence of this metastable modification; in particular we address the question, whether the α-GeTe-type modification might be stabilized (and conversely the rock salt modification destabilized) by steric effects of the non-bonding electron pair.     

How deeply are core electrons perturbed when valence electrons are spin polarized? The case study of transition metal compounds

The ferromagnetic and antiferromagnetic wave functions of the KMnF₃ perovskite have been evaluated quantum-mechanically by using an all electron approach and, for comparison, pseudopotentials on the transition metal and the fluorine ions. It is shown that the different number of α and β electrons in the d shell of Mn perturbs the inner shells, with shifts between the α and β eigenvalues that can be as large as 6 eV for the 3s level, and is far from negligible also for the 2s and 2p states. The valence electrons of F are polarized by the majority spin electrons of Mn, and in turn, spin polarize their 1s electrons. When a pseudopotential is used, such a spin polarization of the core functions of Mn and F can obviously not take place. The importance of such a spin polarization can be appreciated by comparing (i) the spin density at the Mn and F nuclear position, and then the Fermi contact constant, a crucial quantity for the hyperfine coupling, and (ii) the ferromagnetic–antiferromagnetic energy difference, when obtained with an all electron or a pseudopotential scheme, and exploring how the latter varies with pressure. This difference is as large as 50% of the all electron datum, and is mainly due to the rigid treatment of the F ion core. The effect of five different functionals on the core spin polarization is documented.     

The polymorphism of poly(vinylidene fluoride). I. The effect of head-to-head structure

The phenomenon of polymorphism in poly(vinylidene fluoride) has been observed recently by several authors. It has also been reported that high-resolution NMR measurements demonstrate the presence in this polymer of head-to-head linkages, resulting from the “backward” addition of from 5-6% of the monomer units. Since the van der Waals radii of fluorine (1.35 Å) and hydrogen (1.1-1.2 Å) are similar, the cocrystallization in a polymer chain of units that differ only by the substitution of fluorine atoms for hydrogen atoms is not unexpected. The two polymorphic forms of poly(vinylidene fluoride), examined in this investigation, have different chain conformations. Chains in phase I have a planar zigzag conformation, while chains in phase II are assumed to exhibit a 2₁ helical conformation. The incorporation into the polymer chain of small amounts of tetrafluoroethylene or trifluoroethylene comonomer favored the crystallization of phase I. This is in accord with the relative abilities, deduced from consideration of atomic size, of these comonomers to cocrystallize with vinylidene fluoride units in the two indicated chain conformations of the polymer. Since tetrafluoroethylene units are present in the head-to-head structure in the homopolymer, it can be concluded that the elimination of the head-to-head structure will eliminate or restrict crystallization in phase I.