Ab Initio Calculations of the Main Crystal Surfaces of Baryte (BaSO4

Synopsis. We revisited the theoretical morphology of baryte for determining the character and different terminations of the crystal faces, and present...
0 downloads 0 Views 3MB Size
Subscriber access provided by Kaohsiung Medical University 4

Ab Initio Calculations of the Main Crystal Surfaces of Baryte (BaSO) Erica Bittarello, Marco Bruno, and Dino Aquilano Cryst. Growth Des., Just Accepted Manuscript • DOI: 10.1021/acs.cgd.8b00460 • Publication Date (Web): 24 Apr 2018 Downloaded from http://pubs.acs.org on May 8, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

Ab Initio Calculations of the Main Crystal Surfaces of Baryte (BaSO4) BITTARELLO Erica,*1 BRUNO Marco,2 AQUILANO Dino2

1

Dipartimento di Chimica, Università degli Studi di Torino, Via Pietro Giuria 7, 10125 Torino, Italy

2

Dipartimento di Scienze della Terra, Università degli Studi di Torino, Via Valperga Caluso 35, 10125 Torino, Italy *Corresponding author [email protected]

ABSTRACT

In this work, we revisited the theoretical morphology of baryte for determining the character (Flat, Stepped and Kinked) and the different profiles (or terminations) of the crystal faces, and presented an accurate ab initio study of the structures and surface energies at 0 K of the low-index (001), (011), (210), (101), (100), (010), (201), (012) and (211) baryte faces, by using the hybrid Hartree−Fock/density functional B3LYP Hamiltonian. According to the surface energy values, the stability order of the baryte faces was found to be (210) ≅ (001) < (211) < (101) < (010) < (011) < (201) ≅ (100) < (012). Then, the equilibrium shape of baryte was drawn and compared with the previous ones obtained at an empirical level. This is a preliminary work, which is necessary to successively evaluate the influence of water and/or impurities adsorption on the most important baryte surfaces, both at equilibrium and during growth. Indeed, to evaluate these adsorption effects on the crystal morphology, it is fundamental an in-depth analysis of the structure and energetics of the crystal faces acting as substrates.

1 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 35

1. INTRODUCTION Baryte, BaSO4, is an orthorhombic sulfate, named after the Greek βαρ́ ς (heavy), due to its unusual heaviness for a non-metallic mineral.1 Baryte is a common gangue mineral in low-temperature hydrothermal veins and is typical in deposits derived from weathered baryte-bearing limestones; it is also an accessory mineral in igneous rocks and a primary component of submarine volcanogenic sulfide deposits. Baryte is commonly in association with fluorite, calcite, dolomite, rhodochrosite, gypsum, sphalerite, galena and stibnite. Usually found as thin to thick tabular {001} crystals, bounded by {210} alone or in combination with {101}, {011} or other forms. It also can display flattened {001}, and elongated to prismatic [010] or [100], and more rarely prismatic [001], or equant. Often baryte appears as massive aggregates or clusters of tabular crystals with edges projecting into crest-like forms, or as rosettes-like aggregates. Baryte is isostructural with celestine (SrSO4) and anglesite (PbSO4), its crystal structure being first determined by James and Wood.2 Sahl,3 Colville and Staudhammer4 and Hill5 refined the structural model in space group Pnma. The unit-cell parameters5 for baryte are: a0 = 8.8842(12), b0 = 5.4559(8), c0 = 7. 1569(9) Å, V = 346.904Å3, Z = 4. The structure of the baryte-celestine series was successively refined by Miyake et al.6, Jacobsen et al.7 and Antao8 to determine the structural trends that are expected across the series and the structure refinements were carried out in space group Pbnm. In this work we choose to use the crystallographic parameters by Hill 5. Commonly, baryte is used in the industries for the production of barium-hydroxide for sugar refining and as a white pigment for paper, textiles and paint.9 It is also suitable for other purposes due to its high specific gravity, opaqueness to X-rays, inertness, and whiteness10. Recently, controlled morphology and crystal growth of massive production of BaSO4 micro- and nanoparticles have been closely studied for multiple applications in medical area, in oil industry,11 electronics, TV screen, glass, car filters, paint industry and ceramics.12 Moreover, baryte is also used as space filler in automotive brake friction materials.13 The amount, the morphology and size of the baryte crystals affect the tribological properties of these materials. Therefore, the ability to control the size and morphology of the baryte crystals is of fundamental importance to produce highly efficient materials for the industry. The rich literature on baryte can be roughly divided in two large domains: experimental and theoretical. The experimental works can be, in turn, sub-divided in: I. precipitation and growth, including an averaged mass morphology14-23 2 ACS Paragon Plus Environment

Page 3 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

II. detailed nucleation and growth processes, including surface morphology (during both growth and dissolution)24-32 III. influence by organic and inorganic additives,33-37 or surfactants38 on the growth features of baryte. The theoretical works are mainly devoted to predict the equilibrium and growth morphology of baryte, both at 0 K and room temperature, either with or without surface relaxation.39-48 Special attention should be deserved as well to a kind of papers in which an accurate surface microtopography has been integrated by theoretical aspects concerning both the overall growth morphology and the kinetics of the steps building up the surface patterns (spirals and/or 2Dnuclei).49-55 In this work we aim at revising the theoretical morphology of baryte, (i) by applying the Hartman and Perdok39 analysis to determine the character (Flat, Stepped and Kinked) and the different profiles (or terminations) of the crystal faces, and (ii) by performing quantum-mechanical calculations to obtain their optimized structure and surface energy at 0 K and in the vacuum. We employed a hybrid Hartree-Fock (HF)-density functional theory (DFT) approach that, notably, has never before been applied to the study of baryte surfaces; in particular, the B3LYP functional was used. The following crystal faces were taken into account: (001), (011), (210), (101), (100), (010), (201), (012) and (211). This preliminary work is necessary to successively evaluate the influence of water and/or impurities adsorption on the most important baryte surfaces, both at equilibrium and during growth. Indeed, to evaluate these adsorption effects on the crystal morphology, it is fundamental an in-depth analysis of the structure and energetics of the crystal faces acting as substrates. The optimized surfaces structures to be determined in this work will be used as well, in a forthcoming paper, to calculate the adsorption energies released when water molecules or specific impurities are placed in contact with the crystal faces. This article is structured as follows: section 2 provides an outline of the computational methodology. In sections 3: (i) the Hartman and Perdok39 analysis on the main crystal faces of barite is performed, and (ii) their optimized structures and surface energies obtained at ab initio level are described and compared with the results of previous computational studies.

3 ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 35

2. COMPUTATIONAL DETAILS The calculations were performed with the ab initio CRYSTAL14 code56-59 and at the DFT (Density Functional Theory) level with the B3LYP Hamiltonian.60-62 In CRYSTAL14, the multi-electronic wave-function is constructed as an anti-symmetrized product (Slater determinant) of mono-electronic crystalline orbitals (COs) which are linear combinations of local functions (i.e.: atomic orbitals, AOs) centred on each atom of the crystal. In turn, AOs are linear combinations of Gaussian-type functions (GTF, the product of a Gaussian times a real solid spherical harmonic to give s-, p- and d-type AOs). In this study, oxygen was described by (8s)(411sp)-(1d) contraction63: the exponents (in bohr−2 units) of the most diffuse sp-shells are 0.59 and 0.25; the exponent of the single Gaussian d-shell is 0.5. For barium atoms, we used the pseudopotential by Hay and Wadt 64 in combination with the tightest sp-contraction from Piskunov et al.,65 and three sp-shells with exponent 0.95, 0.45, 0.17, and two d-shells with exponents 0.4 and 0.15.66 In order to test the effect of the basis set of sulphur on the geometry and energy of the crystal surfaces, we performed the calculations by using: I.

BS1: a pseudopotential67 to describe the first ten electrons and (31sp)-(1d) contraction to describe the remaining six electrons; the exponent of the diffuse sp-shell is 0.13 and that of the d-shell is 0.21.

II.

BS2: an all Gaussian-type basis sets recently developed68; it consists of a TZVP contraction that has been systematically derived for solid state systems starting from def2-TZVP69-70 basis sets devised for molecular cases.

It is important to highlight that others basis set for Ba and O were tested before to choose those above reported. Unfortunately, not one was able to reproduce in a satisfactory way the bulk structure of barite: the distances in the SO4 tetrahedra obtained at quantum-mechanical level are on average ∼7% higher than those experimentally measured. Therefore, we choose a combination of basis set that better reproduces the bulk structure and, at the same time, allows to perform the calculations in a reasonable time, compatibly with the resources of calculus actually in our hand. Further computational details (e.g., thresholds controlling the accuracy of the calculations) are given as Supporting Information. The CRYSTAL14 output files, listing the optimized fractional coordinates and optimized 2D cell parameters of the (001), (011), (210), (101), (100), (010), (201), (012) and (211) slabs, are freely available at http://mabruno.weebly.com/download. All of these slabs are charge neutral and retain either the centre of inversion or the mirror plane parallel to the (hkl) face, to ensure that the dipole 4 ACS Paragon Plus Environment

Page 5 of 35 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

moment perpendicular to the slab is equal to zero. The calculations were performed by considering 2D slabs56-57 with a thickness sufficient to obtain an accurate description of the surfaces. The slab thickness is considered appropriate when the bulk-like properties are reproduced at the centre of the slab. Further details are given as Supporting Information. The specific surface energy γ (J/m2) at T = 0 K was calculated by means of the relation:56-57

γ = lim Es(n) = lim →∞

→∞

()     

(1)

where E(n) and Ebulk are the energy of a n-layer slab and of the bulk, respectively; A is the area of the primitive unit cell of the surface; the factor 2 in the denominator accounts for both upper and lower surfaces of the slab. Es(n) is thus the energy (per unit area) required for the formation of the surface from the bulk. When n → ∞, Es(n) will converge to the surface energy per unit area (γ). The number of dhkl layers to be considered in each slab, n, was set to values such that the relative difference between Es(n) and Es(n − 1) was