STRUCTURE AND PROPERTIES OF MD SIMULATED Na 2 0.Si0 2 MELT COMPARISON OF THE BORN-MAYER-HUGGINS AND PAULING INTERIONIC POTENTIALS
|
|
- Amberlynn Phelps
- 5 years ago
- Views:
Transcription
1 Ceramics - Silikaty 37, s (1993) 83 STRUCTURE AND PROPERTIES OF MD SIMULATED Na 2 0.Si0 2 MELT COMPARISON OF THE BORN-MAYER-HUGGINS AND PAULING INTERIONIC POTENTIALS BEATA HATALOVA, MAREK LISKA Institute of Inorganic Chemistry of the Slo.vak Academy of Sciences, Dubravska cesta 9, Bratislava Received MD.simulations of Na 20.Si0 2 melt were performed using Born-Mayer-Huggins and Pauling pair potential functions in the temperature range of K. Thermodynamic properties as total energy, pressure and heat capacity were compared for the two simulated systems as well as their structural characteristics. Thermodynamic behaviour of both systems is similar and the calculated values of the heat capacity are close to the experimental one. There are some differences between the simulated structures: Born-Mayer-Huggins pair potential functions support occurence of some silicon atoms coordinated by five oxygens and more branched structure of silicate polyanions with many rather small rings. On the other hand, the structure simulated using Pauling pair potential functions has more chain-like character with few large rings and there are no 5-coordinated silicon atoms present. INTRODUCTION In recent years the method of molecular dynamics (MD) has been widely used to elucidate the random structure of silica and silicate melts and glasses as well as some of their properties [14,15]. All these simulations are based on ionic approximation and two types of interionic potentials have been mostly applied, the Born-Mayer-Huggins potential (BMH) parametrized by Woodcock [1] and the Pau\ing potential (PA) parametrized by Mitra[2]. All parameters were adjusted to fit the properties of Si0 2 and both potentials do work well enough to reproduce the local tetrahedral network structure of amorphous silica. Some problems occure when alkalis are added to Si0 2 and become more obvious with increasing concentration of alkali oxide. They concern the presence of 5-coordinated silicon atoms and incre ing of pressure in simulated systems. Different authors reported the results for alkali silicates [3,4,16.::..19] but always using only one of the mentioned potentials. Comparison of these results is rather difficult because not only the content of alkali oxides was different, but the cooling procedures varied seriously, as well. The aim of this work is to show the differences in simulated Na 2 0.Si0 2 structure and properties using BMH and PA pair potential functions (PPFs) with parameters that were employed in simulations of sodium silicate systems by Soules [3] and Mitra [4], respectively. METHOD The system of 378 particles (63 Si, and 126 Na) was simulated applying periodic boundary conditions. The density of the system was kept constant. It was obtained by linear extrapolation from experimental data [5] for temperature 1873 K, that is 2.11 g/ cm 3 Simulations were performed for two models of interionic interactions; the first represented by BMH PPFs (-r ) z z e2 ½ =Ai exp ----=! J J {! 411" otij ( u + u ) A;j=(l+ z;/v;+zi/vi)bexp {! 3 (1) ( ) and the second by PA PPFs defined as interparticle force: -- - ZiZje 2 Tij [1. (. ) n <r;+o'j ( ).L' J + sign z,z - 411'Eor 3 3 ;j Tij ] ( 3 ) where r;i is the distance between the i-th and j-th particles, z the electronic charge, u the ionic radius, v th number of valence shell electrons, b and {! are constants and n is the measure of the repulsion force. Parameters of the individual ion pairs [3,4] are summarized in Table I. The simulations started from a random configuration at 5000 K. After 2000 time steps (2ps) the system was cooled down to 3000, then 2500 and 2000 K by scaling the velocities of all particles. Each dicrease of temperature was followed by 2ps long equilibration. From 2000 to 1000 K the cooliug procedure comprised sudden decrease of the temperature by 100 K followed by lps equilibration (with numeric control of the set temperature) and another 1 )>S of measuring phase. Individual simulation runs were carried out to obtain results for comparison with X-ray diffraction experiment [6] at 1873, 1673 and 1523 K. They started from the configurations reached in the previous simulations for the closest temperature and employed the same cooling procedure except the temperature was set to the required values. All calculations were performed using the P3M3DP program [7] that comprises the calculation of total Coulomb force by particle-particle particle-mesh
2 84 B. Hatalova, M. Liska Table I Parameters of the used pair potential functions Potential type Parameter Si 0 Na Zi Ui [A] Born-Mayer Vi Huggins b [J] 3.38x10-20 a [A] 0.29 Zi Pauling u; [A] n 10 method [8]. The numerical integration time step was s. RESULTS AND DISCUSSION Thermodynamic properties As it is obvious from the values of ionic charges the total energy of the system simulated with BMH PPFs is about 2.8 times greater than with PA PPFs. Using the linear least squares method the values of isochoric -4.15E+007 heat capacity cv and thermal pressure coefficient 'YV were calculated from the temperature dependances of total energy and pressure of the simulated systems, respectively. These are shown in Figs. 1 and 2. The isobaric heat capacity C p of the simulated systems was calculated according to the following equation c p = cv + <Xp"fvT (4) where cv and 'YV were taken from the simulation results (as mentioned above), the thermal expansivity <Xp was calculated from the density measurements [5], and temperature 1500 K was taken as the mean value of the temperature interval K for which cv was determined. The results and comparison with experiment are shown in Table II. The calculation of pressure from the virial theorem is not accurate enough, so in MD simulations the pressure up to 0.5 GPa is considered the normal pressure [2]. In both simulated systems the pressure was higher. It ranged from 0.8 to 3 GPa. As it has been shown in other works, pressure would lower if the density of simulated system was lowered. For example, MD study of the molar volume of Na 2 0.3Si0 2 using BMH PPFs for description of interparticle interactions [9] gave the molar volume about 22 cm 3 /mol for the melt at 2000 K. That is the density of 1.57 g/ cm 3 only, comparing with the value g/cm 3 extrapolated from the density measurements [10]. Similar problems with the density of sodium trisilicate were reported by Murray et al. who used combined Keat E+007 :;,, -4.35E+007 i -1.21E "i::.c t ;,s.oo -1.22E BMH PPFB + PA PPFs + BMH PPF8 PA PPFa -1.23E+00 8s 0 f500 temper«ture 2 00 (K) t-.....,......,-,-.,...,..,...,...,,-,-..., temperature 2600 (K) Fig. 1. Temperature dependance of the total energy in Na2 O.Si0 2 melt simulated using Born-Mayer-Huggins (BMH) and Pauling (PA} pair potential functions (PPFs). Fig. 2. Temperature dependance of the pressure in Na 20.Si0 2 melt simulated using Born-Mayer-Huggins {BMH) and Pauling {PA) pair potential functions (PPFs).
3 Structure and Properties of MD Simulated Na2 O.Si02 Melt Table II Calculated and experimental values of the thermal pressure coefficient 'YV, thermal expansivity coefficient <Xp, isochoric anq. isobaric heat capacity cv and cp, respectively Potential 'YV <X p cv Cp type [MPa K- 1 ] [m a kg-1k-1] (J kg- 1 I( -1 ] (J kg- 1 K-1] BMH PA Experiment x10-8 * ** *Ref. (5] **Ref. (6] ing and Lenard-Jones potential [11] for modelling its structure. Such high pressure was not observed in cases when Si0 2 was simulated in accordance with the small thermal expansivity of silica. Application of the potential functions parametrized for Si0 2 to alkali silicate systems may lead to mentioned increase of pressure, as well. Despite the high pressure in the simulated systems the values of the heat capacity and thermal pressure coefficient are reasonable comparing with experiment. Struct ure The main features of the structure of simulated Na20.Si02 melt for both types of PPFs are similar. In agreement with general opinion the structure is built of Si04 tetrahedra that are connected through their corners. The tetrahedra are only slightly distorted, as the mean OSiO bond angle is close to the value for regular tetrahedron The mean values and distributions of the SiOSi bond angle for both simulated systems refl.exing the arrangement of tetrahedra are in Table IV and Fig. 6, as well as the results for intertetrahedral bond angle OSiO. Sodium ions are placed in intermediate sites of the silicate framework. Their coor.dination sphere created by oxygens is not regular and the number of neighbouring oxygens vary. Fig. 9. The structure of Na2 O.Si02 melt simulated using BMH PPFs at 1679 K. Fig. 4. The structure of Na2 O.Si02 melt simulated using PA PPFs at 1679 K. Ceramics - Silika.ty c. 2, 1993
4 86 B. Hatalova, M. Liska Table III Temperature dependance of the parameters of short range order, average distance Tij [10-10 m] and coordination number N;j, in molten Na20.Si02 a:n,d its MD models simulated using Born-Mayer-Huggins (BMH) and Pauling (PA) potential functions T [K] SiO NaO 0 0 Si Si Tij N;i Tij N;i Tij N;i Tij N;i BMH PA Exp [6] Comparison of the short range structure characteristics, i.e. bond lengths and coordination numbers, with X-ray diffraction experiment [6] is presented in Table III. The bond lengths represented by positions of the first maxima in the pair radial distribution functions (RDFs) agree quite well for all pairs. Significant differences between the simulations and experiment occured in coordination numbers. Though the experimental error for coordination numbers is 0.3, the differences for 00 and SiSi pairs are much bigger. The values for the MD modeled melts are by lower than those determined from the X-ray measurement. However, from the theoretical point of view Na20.Si02 is built of tetrahedra with two bridging and two non-bridging oxygens that are connected into chains. That is the structure of crystalline sodium silicate, as well,_and this chain-like structure supports the coordination number SiSi equal 2 rather than from the difraction results (Table III). The 2 9 Si NMR measurements of sodium silicate glass support this value, too, because only tetrahedra with two bridging oxygens were observed [12]. The most distinct differences in the coordination number for NaO pairs can be explained by the character ofnao pair RDF (Fig. 5). The coordinatio.n sphere of Na ion is not well defined, so it is difficult to set the cutoff distance for its first coordination sphere: If the cut-off distance is set to the first flat minimum after first maximum in the NaO RDF (these results), the coordination number comprises perhaps only the non-bridging oxygen neighbours. The short range order parameters are constant. in the studied temperature interval (Table III). Comparison of the structures simulated using BMH and PA PPFs yields the following results: 1) BMH PPFs support occurence of some silicon atoms that are coordinated by 5 oxygens. This can be identified by rising of the mean coordination number over the value of 4 (Table III), as well as in the 0Si0 bond angle distribution function (Fig. 6). A small shoulder.appears at the lower angle side of the main peak together with a very small and flat peak at They represent the angles in a deformed trigonal bipyramid which arises from a coordination tetrahedron when the fifth oxygen enters coordination sphere of a silicon atom [13]. The presence of Table IV The mean bond angles fo.r-the simulated structure at 1673 K. The numbers in brackets are the standard deviations Potential Bond angle [ 0 ] type OSiO SiOSi BMH {10.2) (11.2) PA (7.7) (13.6) '
5 Structure and Properties of MD Simulated Na2 O.Si02 Melt BMH potentia,l P.A potentia.i Si0 0Si0 I :,'\ {\ i \ 4.0 f \ Si0Si 1 I \ I I j \ I \ ; \ ) \. \./"11 1 -' ,...:,..., angle (degree).s 2.0 i \ ; \ / i i i i :\ ( \./ Fig. 6. Bond angle distributions of the simulated Na2 O.Si02 melts. Table V Percentage of species in the simulated melts at 1673 K Potential type Coordination species no n1 n2 n3 n4 ns 20.0 BMH PA coordination might be the result of inadequate ratio of the anion/cation radii, that is 1.23 for BMH parametrization, though it should be at least 2.44 [20]. k. 2) The two structures differ in the bridging to non- bridging oxygen ratio and in the distribution of O n species that represent tetrahedra with n bridging oxygens [12]. The species rls, present in the structure simulated using BMH PPFs, is the case of a 5- coordinated silicon with all bridging oxygens. While the theoretical value of the bridging to non-bridging oxygen ratio for Na 2 0.Si0 2 is 0.5, we obtained 0.59 and 0.55 for BMH and PA PPFs, respectively. The distribution of O n species is in Table V, and shows, that the structure obtained using BMH PPFs is more branched than the other one. It consists of many smaller rings including 4-6 tetrahedra while the structure obtained using PA PPFs has more chain-like character with small number of large rings with 8 to 15 tetrahedra. Comparing these facts with 29Si NMR i;esults [12] that confirm only presence of 0 2 species and admit less than 10% of the others, the structure obtained using PA PPFs seems more realistic NaO 0.0 -t.tt"rnri-t.,.-j"'r1""1'"t"1-frtt"rrrri..,...,-tt"t"trri-rt.,...,.,-.-, o.o 'o e.o Fig. 5. Pair radial dtfll1'1l'bht1(/rfj;mjtions (RDF) of Na20.Si02 melt simulated using BMH PPFs at 1673 K. 00
6 88 B. Hatalova, M. Liska CONCLUSIONS The MD models of sodium silicate melt obtained using Born-Mayer-Ituggins and Pauling PPFs [3,4] for description of interparticle interactions agree in the main features rather well with experimental findings, though there are some differences between both simulated structures. Thermodynamic properties of the systems, except the total energy, are similar for both potential types and their parametrizations. The heat capacities and thermal pressure coefficients have reasonable values, but the pressure is quite high. The first step to improvement of MD results should be reparametrization of the present PPFs because their parameters, as in general all fitted ones, may depend on chemical composition of a simulated system. References [1] Woodcock L. V., Angell C. A., Cheeseman P.: J. Chem. Phys. 65, 1565 (1976). [2] Mitra S. K: Phil. Mag. B45, 529 (1982). [3] Soules T. F.: J. Chem. Phys. 71, 4570 (1979). [4] Mitra S. K., Hockney R. W. in: The Structure of Non Crystalline Materials,. (ed. by Gaskell P. H., Parker, J. M.), p. 316, Taylor and Francis Ltd., New York [5] Coenen M. in: Svojstva Stekol i Stekloobrazujushtchich Rasplavov I., ( ed. by Mazurin 0. V., Strelcina M. V., Schvajko-Schvajkovskaja T. P.), p. 210, N auka, Leningrad [6] Waseda Y., Toguri J. M.: Transactions ISIJ 17, 601 (1977). [7] Eastwood J. W., Hockney R. W., Lawrence D. N.: Computer Phys. Commun. 19, 215 (1980). [8] Eastwood J. W. in: Computational Methods in Classical and Quantum Physics (ed.by Hooper M. B.), p. 206, Advance Publications, London [9] Soules T. F.: J. Non-Cryst. Splids 49, 29 (1982). [10] Heidtkamp G., Endell K. in: Svojstva Stekol i Stekloobrazujushtchich Rasplavov I., (ed. by Mazurin 0. V., Strelcina M. V., Schvajko-Schvajkovskaja T. P.), p. 204, N auka, Leningrad [11] Murray R. A., Song L. W., Ching W. Y.: J. Non Cryst.. Solids 94, 133 (1987). [12] Dupree R., Holland D., McMillan P. W., Pettifer R. F.: J. Non-Cryst. Solids 68, 399 (1984). [13] Hatalova B., Liska M: J. Non-Cryst. Solids 146, 218 (1992). [14] Soules T. F. in: Glass, Science and Technology, Vol. 4A, Structure, Microstructure and Properties, (ed. by Uhhnan D. R., Kreidl N. J.), p. 267, Academic Press, New York, [15] Hatalova B., Liska M: Ceramics-Silikaty 32, 359 (1988). [16] Angel C. A., Cheeseman P.A., Tamaddon S.: Science 218, 885 (1982). [17] Tesar A. A., Varshneya A. K.: J. Chem. Phys. 87, 2986 (1987). [18] Newell R. G., Feuston B. P., Garofalini S. H.: J. Mater. Res. 4, 434 (1989). [19] Inoue H., Yasui I.: Phys. Chem. Glasses 28, 63 (1987). [20] Dekker A. J.: Solid State Physics, p. 127, Prentice Hall, Englewood Cliffs Submitted in English by the authors MD MODEL TAVENINY Na 2 0.Si0 2 POROVNANIE BORN-MAYER-HUGGINSOVYCH A PAULINGOVYCH INTERAKCNYCH POTENCIJ\LOV BEATA HATALOVA, MAREK LISKA Ustav anorganickej chemie SA V, Dubravska cesta 9, Bratislava V teplotnom intervale K bola met6dou molekulovej dynamiky simulovana tavenina N a2 0.Si02. V jednom subore boli medzicasticove sily reprezentovane Born-Mayer-Hugginsovymi (BMH) a v druhom Paulingovymi (PA) parovymi potencialovymi funkciami. Odhliadnuc od celkovej energie, ktora je v dosledku vysokych nabojov v pripade BMH potencialov 2.8-krat vii.csia, SU termodynamicke vlastnosti oboch suborov velhti podobne. V oboch je prilis vysoky tlak, co vsak nema vplyv na hoduoty izochorickej teplotnej rozpinavosti a merneho tepla. Napriek dobrej zhode oboch struktur s vysledkami rontgeuovej difrakcnej analyzy, navzajom Sa V detailoch odlisuju. BMH potencialove funkcie vedu k rozvetvenej strukture s velkym poctom pomerne malych cyklov vytvorenych tetraedrami Si04, V ktorej sa vsak nachadzaju aj at6my kremika koordinovane piatimi kyslikmi. Struktura, ktora vznika pri simulacii s pouzitfm PA poteucialovych funkcif, je charakterizovana refazcami prepojenymi leu do niekolkych velkych cyklov a vsetky at6my kremika Sa nachadzaju V tetraedrickej koordinacii kyslfkom. Obr. 1. Teplotna zavislosf celkovej energie v simulovanych systemoch Na2 O.Si02 s pouiitfrn Born-Mayer Hugginsovych (BMH) a Paulingovych (PA} parovych potenicalovych funkci( (PPFs). Obr. 2. Teplotna zavislose tlaku v simulovanych systemoch Na20.Si02 s pouiium BMH a PA PPFs. Obr. 3. Struktura taveniny Na2 O.Si02 simulvoanej s pouiium BMH PPFs pri teplote 1673 K. Obr. 4. Struktura taveniny Na2 O.Si02 simulvoanej s pouiium PA PPFs pri teplote 1673 I(. Obr. 5. Parove radialne distribucne funkcie pre taveninu Na20.Si02 simulovanu s pouiitim BMH PPFs pri 1673 I(. Obr. 6. Distribucie viizbovych uhlov v simulovanej tavenine Na20.Si02.
The structure of sodium trisilicate glass via molecular dynamics employing three-body potentials
The structure of sodium trisilicate glass via molecular dynamics employing three-body potentials R. G. ewell, B.P. Feuston, and S.H. Garofalini Department of Ceramic Engineering, Rutgers University, Piscataway,
More informationA Molecular Dynamics Simulation of a Homogeneous Organic-Inorganic Hybrid Silica Membrane
A Molecular Dynamics Simulation of a Homogeneous Organic-Inorganic Hybrid Silica Membrane Supplementary Information: Simulation Procedure and Physical Property Analysis Simulation Procedure The molecular
More informationTHE STRUCTURE OF THE FIRST CO-ORDINATION SHELL OF Pb ATOMS IN LEAD-SILICATE GLASSES: A MOLECULAR DYNAMICS STUDY
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 5, 67-74 (1999) THE STRUCTURE OF THE FIRST CO-ORDINATION SHELL OF Pb ATOMS IN LEAD-SILICATE GLASSES: A MOLECULAR DYNAMICS STUDY ANNA RYBICKA, JAROSŁAW RYBICKI,
More informationTHE STRUCTURE OF POROUS AND SPONTANEOUSLY DENSIFIED AMORPHOUS PbSiO 3 : A MOLECULAR DYNAMICS STUDY
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 10(1), 21-38 (2004) THE STRUCTURE OF POROUS AND SPONTANEOUSLY DENSIFIED AMORPHOUS PbSiO 3 : A MOLECULAR DYNAMICS STUDY KATARZYNA CHOMENKO 1, GRZEGORZ BERGMAŃSKI
More informationShell-model molecular dynamics calculations of modified silicate glasses
Shell-model molecular dynamics calculations of modified silicate glasses Antonio Tilocca, 1, * Nora H. de Leeuw, 1,2 and Alastair N. Cormack 1,3 1 Department of Chemistry, University College London, London
More informationMolecular Dynamics Simulation of the Thermodynamic and Structural Properties for the CaO SiO 2 System
ISIJ International, Vol. 44 (4), No. 11, pp. 1817 185 Molecular Dynamics Simulation of the Thermodynamic and Structural Properties for the Ca Si System Won-Gap SE and Fumitaka TSUKIHASHI Department of
More informationAtomistic structure of calcium silicate intergranular films between prism and basal planes in silicon nitride: A molecular dynamics study
Atomistic structure of calcium silicate intergranular films between prism and basal planes in silicon nitride: A molecular dynamics study Xiaotao Su and Stephen H. Garofalini Department of Ceramics and
More informationEarth Solid Earth Rocks Minerals Atoms. How to make a mineral from the start of atoms?
Earth Solid Earth Rocks Minerals Atoms How to make a mineral from the start of atoms? Formation of ions Ions excess or deficit of electrons relative to protons Anions net negative charge Cations net
More informationChapter 3. The structure of crystalline solids 3.1. Crystal structures
Chapter 3. The structure of crystalline solids 3.1. Crystal structures 3.1.1. Fundamental concepts 3.1.2. Unit cells 3.1.3. Metallic crystal structures 3.1.4. Ceramic crystal structures 3.1.5. Silicate
More informationT=293 K is shown in Fig. 1. The figure shows that in. Crystalline fragments in glasses ' ' PHYSICAL REVIEW B VOLUME 45, NUMBER 13
), PHYSICAL REVIEW B VOLUME 5, NUMBER 13 1 APRIL 1992-I Crystalline fragments in glasses Giomal A. Antonio Departamento de Fisica, Universidade Fstadual Paulista Bauru, Caixa Postal 73, 1733 Bauru, So
More informationTHE HIGH TEMPERATURE RAMAN SPECTROSCOPY AND ION CLUSTER THEORY FOR METALLURGICAL SLAG AND MAGMA
THE HIGH TEMPERATURE RAMAN SPECTROSCOPY AND ION CLUSTER THEORY FOR METALLURGICAL SLAG AND MAGMA Guo-Chang Jiang, Yong-Quan Wu & Jing-Lin You Shanghai University, China ABSTRACT This paper presents the
More informationStructural properties of low-density liquid alkali metals
PRAMANA c Indian Academy of Sciences Vol. 65, No. 6 journal of December 2005 physics pp. 1085 1096 Structural properties of low-density liquid alkali metals A AKANDE 1, G A ADEBAYO 1,2 and O AKINLADE 2
More informationIonic Coordination and Silicate Structures
Ionic Coordination and Silicate Structures Pauling s Rules A coordination polyhedron of anions forms around a cation Ionic distance determined by radii Coordination number determined by radius ratio. May
More informationDensity, Volume, and Packing: Part 2 Thursday, September 4, 2008 Steve Feller Coe College Physics Department
Density, Volume, and Packing: Part 2 Thursday, September 4, 2008 Steve Feller Coe College Physics Department Lecture 3 of Glass Properties Course Some Definitions x = molar fraction of alkali or alkaline-earth
More informationStructural properties of low density liquid alkali metals
PRAMANA c Indian Academy of Sciences Vol. xx, No. x journal of xxxxxx 2005 physics pp. 1 12 Structural properties of low density liquid alkali metals A AKANDE 1, G A ADEBAYO 1,2 and O AKINLADE 2 1 The
More informationLecture 8 Chemical/Electronic Structure of Glass
Lecture 8 Chemical/Electronic Structure of Glass Syllabus Topic 6. Electronic spectroscopy studies of glass structure Fundamentals and Applications of X-ray Photoelectron Spectroscopy (XPS) a.k.a. Electron
More informationConduction Modeling in Mixed Alkali Borate Glasses
International Journal of Pure & Applied Physics ISSN 0973-1776 Vol.1 No.2 (2005), pp. 191-197 Research India Publications http://www.ripub lication.com/ijpap.htm Conduction Modeling in Mixed Alkali Borate
More information5. STRUCTURES AND DEFECTS IN AMORPHOUS SOLIDS
62 5. STRUCTURES AND DEFECTS IN AMORPHOUS SOLIDS 5.1 Review/Background: In Chapter 4, we discussed the origin of crystal structures and Bravais lattices based on Euler relationship. In this chapter, we
More informationThe dative covalent bond acts like an ordinary covalent bond when thinking about shape so in NH 4. the shape is tetrahedral
1.3 Bonding Definition Ionic bonding is the electrostatic force of attraction between oppositely charged ions formed by electron transfer. Metal atoms lose electrons to form ve ions. Non-metal atoms gain
More informationSolids. properties & structure
Solids properties & structure Determining Crystal Structure crystalline solids have a very regular geometric arrangement of their particles the arrangement of the particles and distances between them is
More informationCOMPUTER SIMULATIONS OF LIQUID SILICA: WATER-LIKE THERMODYNAMIC AND DYNAMIC ANOMALIES, AND THE EVIDENCE FOR POLYAMORPHISM
COMPUTER SIMULATIONS OF LIQUID SILICA: WATER-LIKE THERMODYNAMIC AND DYNAMIC ANOMALIES, AND THE EVIDENCE FOR POLYAMORPHISM IVAN SAIKA-VOIVOD 1 and PETER H. POOLE 2 1 Department of Physics and Physical Oceanography,
More informationCHAPTER 4. Crystal Structure
CHAPTER 4 Crystal Structure We can assume minerals to be made of orderly packing of atoms or rather ions or molecules. Many mineral properties like symmetry, density etc are dependent on how the atoms
More informationLITHIUM COORDINATION MEASUREMENT IN THE LITHIUM DISILICATE GLASSES
Romanian Reports in Physics, Volume 56, No. 3, P. 413-418, 2004 LITHIUM COORDINATION MEASUREMENT IN THE LITHIUM DISILICATE GLASSES D. RADU 1, C. MAZILU 2, B. GRIGORAS 2 1 University Politehnica Bucharest,
More informationMolecular dynamics simulation of lithium diffusion in Li 2 O Al 2 O 3 SiO 2 glasses
Solid State Ionics 166 (2004) 365 373 www.elsevier.com/locate/ssi Molecular dynamics simulation of lithium diffusion in Li 2 O Al 2 O 3 SiO 2 glasses Weiqun Li, Stephen H. Garofalini* Department of Ceramic
More information2.2.2 Bonding and Structure
2.2.2 Bonding and Structure Ionic Bonding Definition: Ionic bonding is the electrostatic force of attraction between oppositely charged ions formed by electron transfer. Metal atoms lose electrons to form
More informationAccurate computation of shear viscosity from equilibrium molecular dynamics simulations. D. NEVINS* and F. J. SPERA
Accurate computation of shear viscosity from equilibrium molecular dynamics simulations D. NEVINS* and F. J. SPERA Department of Earth Science, University of California, Santa Barbara, CA 93106, USA *dean@nevins.org
More informationA NEW CONCEPT ABOUT THE STATIC STRUCTURE AND ALKALI NON-BRIDGING OXYGEN DEPENDENCE OF THE MIXED ALKALI EFFECT IN BISMUTHATE GLASSES
Bulg. J. Phys. 30 (2003) 141 151 A NEW CONCEPT ABOUT THE STATIC STRUCTURE AND ALKALI NON-BRIDGING OXYGEN DEPENDENCE OF THE MIXED ALKALI EFFECT IN BISMUTHATE GLASSES SH. EL-DESOUKI, E.E. SHAISHA, A.A. BAHGAT
More informationLecture 4! ü Review on atom/ion size! ü Crystal structure (Chap 4 of Nesseʼs book)!
Lecture 4! ü Review on atom/ion size! ü Crystal structure (Chap 4 of Nesseʼs book)! 15 C 4+ 42 Si 4+ Size of atoms! Hefferan and O Brien, 2010; Earth Materials Force balance! Crystal structure (Chap. 4)!
More informationHigh Temperature Materials. By Docent. N. Menad. Luleå University of Technology ( Sweden )
Course KGP003 Ch. 12 High Temperature Materials By Docent. N. Menad Dept. of Chemical Engineering and Geosciences Div. Of process metallurgy Luleå University of Technology ( Sweden ) Ceramic materials
More informationProperties of calcium fluoride up to 95 kbar: A theoretical study
Bull. Mater. Sci., Vol. 33, No. 4, August 2010, pp. 413 418. Indian Academy of Sciences. Properties of calcium fluoride up to 95 kbar: A theoretical study CHUN-SHENG WANG School of Traffic and Transportation,
More informationComputer modelling of molten halides using diffraction data
D. K. BELASHCHENKO and O. I. OSTROVSKI. Computer modelling of molten halides using diffraction data. VII International Conference on Molten Slags Fluxes and Salts, The South African Institute of Mining
More information2011 CHEM 120: CHEMICAL REACTIVITY
2011 CHEM 120: CHEMICAL REACTIVITY INORGANIC CHEMISTRY SECTION Lecturer: Dr. M.D. Bala Textbook by Petrucci, Harwood, Herring and Madura 15 Lectures (4/10-29/10) 3 Tutorials 1 Quiz 1 Take-home test https://chemintra.ukzn.ac.za/
More informationCHEM Principles of Chemistry II Chapter 10 - Liquids and Solids
CHEM 1212 - Principles of Chemistry II Chapter 10 - Liquids and Solids 10.1 Intermolecular Forces recall intramolecular (within the molecule) bonding whereby atoms can form stable units called molecules
More informationBonding in Solids. What is the chemical bonding? Bond types: Ionic (NaCl vs. TiC?) Covalent Van der Waals Metallic
Bonding in Solids What is the chemical bonding? Bond types: Ionic (NaCl vs. TiC?) Covalent Van der Waals Metallic 1 Ions and Ionic Radii LiCl 2 Ions (a) Ions are essentially spherical. (b) Ions may be
More informationModeling & Simulation of Glass Structure
Modeling & Simulation of Glass Structure VCG Lecture 19 John Kieffer Department of Materials Science and Engineering University of Michigan 1 Overview Historical perspective Simulation methodologies Theoretical
More informationChapter 10 Liquids and Solids. Problems: 14, 15, 18, 21-23, 29, 31-35, 37, 39, 41, 43, 46, 81-83, 87, 88, 90-93, 99, , 113
Chapter 10 Liquids and Solids Problems: 14, 15, 18, 21-23, 29, 31-35, 37, 39, 41, 43, 46, 81-83, 87, 88, 90-93, 99, 104-106, 113 Recall: Intermolecular vs. Intramolecular Forces Intramolecular: bonds between
More informationEdexcel Chemistry A-level
Edexcel Chemistry A-level Topic 2 - Bonding and Structure Flashcards What are ions? What are ions? Charged particles that is formed when an atom loses or gains electrons What is the charge of the ion when
More informationChapter 12: Structures & Properties of Ceramics
Chapter 12: Structures & Properties of Ceramics ISSUES TO ADDRESS... Bonding and structure of ceramic materials as compared with metals Chapter 12-1 Atomic Bonding in Ceramics Bonding: -- Can be ionic
More informationIonic Bonding. Chem
Whereas the term covalent implies sharing of electrons between atoms, the term ionic indicates that electrons are taken from one atom by another. The nature of ionic bonding is very different than that
More informationDECRIRE LA TRANSITION VITREUSE: Simulations moléculaires et approches topologiques. Matthieu Micoulaut, LPTMC, Paris Sorbonne Universités UPMC
DECRIRE LA TRANSITION VITREUSE: Simulations moléculaires et approches topologiques Matthieu Micoulaut, LPTMC, Paris Sorbonne Universités UPMC C. Yildrim, O. Laurent, B. Mantisi, M. Bauchy «Chasseurs d
More information4. Interpenetrating simple cubic
2 1. The correct structure t of CsClCl crystal is 1. Simple cubic 2. Body centered cubic 3. Face centered cubic 4. Interpenetrating simple cubic If corner as well as the particle at the center are same
More informationCHAPTER 2. Atomic Structure And Bonding 2-1
CHAPTER 2 Atomic Structure And Bonding 2-1 Structure of Atoms ATOM Basic Unit of an Element Diameter : 10 10 m. Neutrally Charged Nucleus Diameter : 10 14 m Accounts for almost all mass Positive Charge
More informationChapter 2 Atoms, chemical bonding, material structure, and physical properties Homework Solutions
Chapter 2 Atoms, chemical bonding, material structure, and physical properties Homework Solutions Concept questions 1. The Pauli exclusion principle says that no two electrons that occupy the same space
More informationIndentation of Silicon: Phase Transition? Indentation of Silicon: Phase Transition?
Indentation of Silicon: Phase Transition? Kallman et al., Phys. Rev. B. 47, 7705 (1993). Smith et al., Acta. Mater. 49, 4089 (2001). MD: NVT 350,000 atoms Cut away of crystalline Si indented with a tetrahedral
More informationRadiation Damage Modeling of Fused Silica in Fusion Systems
1 Radiation Damage Modeling of Fused Silica in Fusion Systems F. Mota 1), M.J. Caturla 2), J.M. Perlado 1), A. Ibarra 3), M. León 3), J.Mollá 3) 1) Instituto de Fusion Nuclear (DENIM) / ETSII / Universidad
More informationCooling rate effects in amorphous Silica: A Computer Simulation Study. Abstract
August 14, 2018 Cooling rate effects in amorphous Silica: A Computer Simulation Study Katharina Vollmayr, Walter Kob and Kurt Binder Institut für Physik, Johannes Gutenberg-Universität, Staudinger Weg
More informationMetallic & Ionic Solids. Crystal Lattices. Properties of Solids. Network Solids. Types of Solids. Chapter 13 Solids. Chapter 13
1 Metallic & Ionic Solids Chapter 13 The Chemistry of Solids Jeffrey Mack California State University, Sacramento Crystal Lattices Properties of Solids Regular 3-D arrangements of equivalent LATTICE POINTS
More informationChapter 2. Atomic Structure
Chapter 2 Atomic Structure 2 6 (a) Aluminum foil used for storing food weighs about 0. g per square cm. How many atoms of aluminum are contained in one 6.25 cm 2 size of foil? (b) Using the densities and
More informationEarth and Planetary Materials
Earth and Planetary Materials Spring 2013 Lecture 3 2013.01.14 14 1 Close Packed Anion Arrays Closest Packing Coordination number (C.N.) : number of anions bonded to a cation larger cation, higher C.N.
More informationThe Periodic Table and Chemical Reactivity
The and Chemical Reactivity Noble gases Less electronegative elements More electronegative elements Then what is electronegativity? The tendency of an atom to attract an electron (or electron density)
More informationPeriodicity SL (answers) IB CHEMISTRY SL
(answers) IB CHEMISTRY SL Syllabus objectives 3.1 Periodic table Understandings: The periodic table is arranged into four blocks associated with the four sublevels s, p, d, and f. The periodic table consists
More informationMolecular dynamics study on surface structure of a-sio 2 by charge equilibration method
Journal of Non-Crystalline Solids 318 (2003) 131 141 www.elsevier.com/locate/jnoncrysol Molecular dynamics study on surface structure of a-sio 2 by charge equilibration method C. Wang, N. Kuzuu *, Y. Tamai
More informationCrystal Structure and Chemistry
Crystal Structure and Chemistry Controls on Crystal Structure Metallic bonding closest packing Covalent bonding depends on orbital overlap and geometry Ionic bonding Pauling s Rules Coordination Principle
More informationSilicate Structures. Silicate Minerals: Pauling s s Rules and. Elemental Abundance in Crust. Elemental Abundance in Crust: Pauling s s Rules
Silicate Minerals: Pauling s s Rules and Silicate Structures February 6, 2007 Elemental Abundance in Crust Fe Ion O 2- Si 4+ Al 3+, 3+ Ca Na + K + Mg mol % 2.6 1.4 mol% x charge 4.8 3.8 2.6 1.4 3.8 Sum
More informationAdvanced Vitreous State - Physical Properties of Glass
Advanced Vitreous State - Physical Properties of Glass Lecture 25: Charge Conduction Properties of Glass: Steve W. Martin Ionic Conduction in Glass - Part 1 Relationship to Glass Structure and Composition
More informationPY2N20 Material Properties and Phase Diagrams
PY2N20 Material Properties and Phase Diagrams Lecture 10 P. Stamenov, PhD School of Physics, TCD PY2N20-10 Modern CMOS pair structure Photolithographic Process CMOS Processing Steps Cu Damascene Process
More informationMolecular Dynamics Simulation of Ionic Melts
Molecular Dynamics Simulation of Ionic Melts David K. Belashchenko* and Oleg I. Ostrovski** *Department of Physical Chemistry, Moscow Steel and Alloys Institute, Moscow, Russian Federation. Telephone:
More informationStructural order in glassy water
Structural order in glassy water Nicolas Giovambattista, 1 Pablo G. Debenedetti, 1 Francesco Sciortino, 2 and H. Eugene Stanley 3 1 Department of Chemical Engineering, Princeton University, Princeton,
More informationStudy of Mechanisms of Ion Transport in Ion Conducting Glasses
Study of Mechanisms of Ion Transport in Ion Conducting Glasses P. Bury a), P. Hockicko a), M. Jamnický b) and I. Jamnický a) a) Department of Physics, Žilina University, 010 26 Žilina, Slovakia (bury@fel.utc.sk)
More informationDefinition: An Ionic bond is the electrostatic force of attraction between oppositely charged ions formed by electron transfer.
3 Bonding Definition An Ionic bond is the electrostatic force of attraction between oppositely charged ions formed by electron transfer. Metal atoms lose electrons to form +ve ions. on-metal atoms gain
More information2. Calculation.Methods The MD simulations have been performed by using a modified central force potential [14] and
T he E ffect o f P re ssu re on th e H y d ro g e n B o n d S tru c tu re o f Liquid W a te r G. Pálinkás Central Research Institute for Chemistry of the Hungarian Academy of Sciences, Budapest. Hungary
More informationEarth Materials I Crystal Structures
Earth Materials I Crystal Structures Isotopes same atomic number, different numbers of neutrons, different atomic mass. Ta ble 1-1. Su mmar y of quantu m num bers Name Symbol Values Principal n 1, 2,
More informationMODELLING THE STRUCTURE OF. NANOPARTICLE EMBEDDING MATRICES : MOLECULAR DYNAMICS IN Li 2 B 4 O 7
MODELLING THE STRUCTURE OF NANOPARTICLE EMBEDDING MATRICES : MOLECULAR DYNAMICS IN Li 2 B 4 O 7 Alain MARBEUF and Janis KLIAVA Centre de Physique Moléculaire Optique et Hertzienne UMR5798, Université Bordeaux
More informationLiquid-Liquid Phase Transitions and Water-Like Anomalies in Liquids
1 Liquid-Liquid Phase Transitions and Water-Like Anomalies in Liquids Erik Lascaris Final oral examination 9 July 2014 2 Outline Anomalies in water and simple models Liquid-liquid phase transition in water
More informationGlass V.M. Sglavo GlassEng - UNITN 2017
Glass Definition material with no long range order matarial characterized by the glass transition phenomenon amorphous solid without long range order and ordered atomic structure, which shows the glass
More information3.012 PS Issued: Fall 2004 Due: pm
3.012 PS 2 3.012 Issued: 09.15.04 Fall 2004 Due: 09.22.04 5pm Graded problems: 1. In discussing coordination numbers and deriving the permitted range of radius ratio, R A / R B, allowed for each ( where
More informationCHAPTER-3 CLASSIFICATION OF ELEMENTS AND PERIODICITY IN PROPERTIES OF ELEMENTS
CHAPTER-3 CLASSIFICATION OF ELEMENTS AND PERIODICITY IN PROPERTIES OF ELEMENTS Mandeleev s Periodic Law:- The properties of the elements are the periodic function of their atomic masses. Moseley, the English
More informationAtoms, Molecules and Minerals
Atoms, Molecules and Minerals Atoms Matter The smallest unit of an element that retain its properties Molecules - a small orderly group of atoms that possess specific properties - H 2 O Small nucleus surrounded
More informationAtomic Arrangement. Primer in Materials Spring
Atomic Arrangement Primer in Materials Spring 2017 30.4.2017 1 Levels of atomic arrangements No order In gases, for example the atoms have no order, they are randomly distributed filling the volume to
More informationBonding. Bringing the atoms together
Bonding Bringing the atoms together More than one atom Until now, we have been consumed with describing individual atoms of elements. However, isolating individual atoms in most elements is an arduous
More informationBonding and Packing: building crystalline solids
Bonding and Packing: building crystalline solids The major forces of BONDING Gravitational forces: F = G m m 1 2 F = attractive forces between 2 bodies G = universal graviational constant (6.6767 * 10
More informationCHEM1902/ N-2 November 2014
CHEM1902/4 2014-N-2 November 2014 The cubic form of boron nitride (borazon) is the second-hardest material after diamond and it crystallizes with the structure shown below. The large spheres represent
More informationIntroduction to Condensed Matter Physics
Introduction to Condensed Matter Physics Crystalline Solids - Introduction M.P. Vaughan Overview Overview of course Crystal solids Crystal structure Crystal symmetry The reciprocal lattice Band theory
More informationSOLID STATE CHEMISTRY
SOLID STATE CHEMISTRY Crystal Structure Solids are divided into 2 categories: I. Crystalline possesses rigid and long-range order; its atoms, molecules or ions occupy specific positions, e.g. ice II. Amorphous
More informationUNIVERSITY COLLEGE LONDON. University of London EXAMINATION FOR INTERNAL STUDENTS. For The Following Qualifications:-
UNIVERSITY COLLEGE LONDON University of London EXAMINATION FOR INTERNAL STUDENTS For The Following Qualifications:- B.Sc. M.Sci. Physics 1B28: Thermal Physics COURSE CODE : PHYSIB28 UNIT VALUE : 0.50 DATE
More informationChapter 10: Liquids and Solids
Chapter 10: Liquids and Solids Chapter 10: Liquids and Solids *Liquids and solids show many similarities and are strikingly different from their gaseous state. 10.1 Intermolecular Forces Intermolecular
More informationAtomic Structure of the Glassy State (Introduction to Glass Science) CERAMIC ENGINEERING 103
Atomic Structure of the Glassy State (Introduction to Glass Science) CERAMIC ENGINEERING 103 Instructor: Office Hours : Textbook: Richard K. Brow, 222 McNutt Hall, 341-6812, brow@umr.edu Tuesdays and Thursdays,
More informationElectron configurations follow the order of sublevels on the periodic table.
Electron configurations follow the order of sublevels on the periodic table. 1 The periodic table consists of sublevel blocks arranged in order of increasing energy. Groups 1A(1)-2A(2) = s level Groups
More informationDeveloping Monovalent Ion Parameters for the Optimal Point Charge (OPC) Water Model. John Dood Hope College
Developing Monovalent Ion Parameters for the Optimal Point Charge (OPC) Water Model John Dood Hope College What are MD simulations? Model and predict the structure and dynamics of large macromolecules.
More information3. Bonding Ionic Bonding
3. Bonding Ionic Bonding Metal atoms lose electrons to form +ve ions. on-metal atoms gain electrons to form -ve ions. Mg goes from 1s 2 2s 2 2p 6 3s 2 to Mg 2+ 1s 2 2s 2 2p 6 goes from 1s 2 2s 2 2p 4 to
More informationCovalent Bonding. In nature, only the noble gas elements exist as uncombined atoms. All other elements need to lose or gain electrons
In nature, only the noble gas elements exist as uncombined atoms. They are monatomic - consist of single atoms. All other elements need to lose or gain electrons To form ionic compounds Some elements share
More informationChapter 11. Intermolecular Forces and Liquids & Solids
Chapter 11 Intermolecular Forces and Liquids & Solids The Kinetic Molecular Theory of Liquids & Solids Gases vs. Liquids & Solids difference is distance between molecules Liquids Molecules close together;
More informationTheoretical approach to the Poisson s ratio behaviour during structural changes in metallic glasses. Abstract
Theoretical approach to the Poisson s ratio behaviour during structural changes in metallic glasses Eloi Pineda Dept. de Física i Enginyeria Nuclear, ESAB, Universitat Politècnica de Catalunya, Campus
More informationStructure of Cement Phases from ab initio Modeling Crystalline C-S-HC
Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC Sergey V. Churakov sergey.churakov@psi.ch Paul Scherrer Institute Switzerland Cement Phase Composition C-S-H H Solid Solution Model
More informationarxiv:cond-mat/ v1 [cond-mat.dis-nn] 27 Mar 1997
arxiv:cond-mat/9703237v1 [cond-mat.dis-nn] 27 Mar 1997 Molecular Dynamics Computer Simulation of the Dynamics of Supercooled Silica J. Horbach, W. Kob 1 and K. Binder Institute of Physics, Johannes Gutenberg-University,
More informationChemical bonds. In some minerals, other (less important) bond types include:
Chemical bonds Chemical bond: force of attraction between two or more atoms/ions Types of bonds in crystals: Ionic bond: electrostatic attraction between two oppositely charged ions. This type of bond
More information2 Bonding and structure Answers to Exam practice questions
Pages 77 80 Exam practice questions 1 a) Bonding between atoms in a metal is the result of electrostatic attractions between positive metal ions in a lattice and delocalised electrons in the outer shell
More informationIonic Bonding. Example: Atomic Radius: Na (r = 0.192nm) Cl (r = 0.099nm) Ionic Radius : Na (r = 0.095nm) Cl (r = 0.181nm)
Ionic Bonding Ion: an atom or molecule that gains or loses electrons (acquires an electrical charge). Atoms form cations (+charge), when they lose electrons, or anions (- charge), when they gain electrons.
More informationEverything starts with atomic structure and bonding
Everything starts with atomic structure and bonding not all energy values can be possessed by electrons; e- have discrete energy values we call energy levels or states. The energy values are quantized
More informationP. W. Anderson [Science 1995, 267, 1615]
The deepest and most interesting unsolved problem in solid state theory is probably the theory of the nature of glass and the glass transition. This could be the next breakthrough in the coming decade.
More informationUniversity of California, Los Angeles, CA , U.S.A. California, Los Angeles, CA , U.S.A.
Cooling-Rate Effects in Sodium Silicate Glasses: Bridging the Gap between Molecular Dynamics Simulations and Experiments Xin Li 1, Weiying Song 1, Kai Yang 1, N M Anoop Krishnan 1, Bu Wang 1, Gaurav Sant
More informationChapter 12. Insert picture from First page of chapter. Intermolecular Forces and the Physical Properties of Liquids and Solids
Chapter 12 Insert picture from First page of chapter Intermolecular Forces and the Physical Properties of Liquids and Solids Copyright McGraw-Hill 2009 1 12.1 Intermolecular Forces Intermolecular forces
More informationUNIT-1 SOLID STATE. Ans. Gallium (Ga) is a silvery white metal, liquid at room temp. It expands by 3.1% on solidifica-tion.
UNIT-1 SOLID STATE 1 MARK QUESTIONS Q. 1. Name a liquefied metal which expands on solidification. Ans. Gallium (Ga) is a silvery white metal, liquid at room temp. It expands by 3.1% on solidifica-tion.
More informationMetallic and Ionic Structures and Bonding
Metallic and Ionic Structures and Bonding Ionic compounds are formed between elements having an electronegativity difference of about 2.0 or greater. Simple ionic compounds are characterized by high melting
More informationINTERMOLECULAR FORCES: LIQUIDS, SOLIDS & PHASE CHANGES (Silberberg, Chapter 12)
INTERMOLECULAR FORCES: LIQUIDS, SOLIDS & PASE CANGES (Silberberg, Chapter 12) Intermolecular interactions Ideal gas molecules act independently PV=nRT Real gas molecules attract/repulse one another 2 n
More information2. Bonding Ionic Bonding
2. Bonding Ionic Bonding Metal atoms lose electrons to form +ve ions. on-metal atoms gain electrons to form -ve ions. Mg goes from 1s 2 2s 2 2p 6 3s 2 to Mg 2+ 1s 2 2s 2 2p 6 goes from 1s 2 2s 2 2p 4 to
More informationImplementation of consistent embedding for a larger system Amorphous silica
Journal of Computer-Aided Materials Design (2006) 13:61 73 Springer 2006 DOI 10.1007/s10820-006-9015-z Implementation of consistent embedding for a larger system Amorphous silica KRISHNA MURALIDHARAN a,,
More informationCeramic Bonding. CaF 2 : large SiC: small
Recall ceramic bonding: - Mixed ionic and covalent. - % ionic character ( f ) increases with difference in electronegativity Large vs small ionic bond character: Ceramic Bonding CaF : large SiC: small
More informationAdvanced Vitreous State - Physical Properties of Glass
Advanced Vitreous State - Physical Properties of Glass Lecture 26: Charge Conduction Properties of Glass: Ionic Conduction in Glass - Part 2 Activation Energies in Glass Steve W. Martin Department of Materials
More informationIntermolecular Forces, Liquids, Solids. IM Forces and Physical Properties
Intermolecular Forces, Liquids, Solids Interactions Between Molecules: What does it take to separate two (or more) molecules from one another? or What holds molecules close to one another? Structure/Property
More information