Inversion of Diffraction Data for Amorphous Materials
Total Page:16
File Type:pdf, Size:1020Kb
The University of Southern Mississippi The Aquila Digital Community Faculty Publications 9-22-2016 Inversion of Diffraction Data for Amorphous Materials Anup Pandey Ohio University Parthapratim Biswas University of Southern Mississippi, [email protected] D.A. Drabold Ohio University Follow this and additional works at: https://aquila.usm.edu/fac_pubs Part of the Physics Commons Recommended Citation Pandey, A., Biswas, P., Drabold, D. (2016). Inversion of Diffraction Data for Amorphous Materials. Scientific Reports, 6, 1-8. Available at: https://aquila.usm.edu/fac_pubs/16741 This Article is brought to you for free and open access by The Aquila Digital Community. It has been accepted for inclusion in Faculty Publications by an authorized administrator of The Aquila Digital Community. For more information, please contact [email protected]. www.nature.com/scientificreports OPEN Inversion of diffraction data for amorphous materials Anup Pandey1, Parthapratim Biswas2 & D. A. Drabold3 The general and practical inversion of diffraction data–producing a computer model correctly Received: 06 July 2016 representing the material explored–is an important unsolved problem for disordered materials. Such Accepted: 31 August 2016 modeling should proceed by using our full knowledge base, both from experiment and theory. In this Published: 22 September 2016 paper, we describe a robust method to jointly exploit the power of ab initio atomistic simulation along with the information carried by diffraction data. The method is applied to two very different systems: amorphous silicon and two compositions of a solid electrolyte memory material silver-doped GeSe3. The technique is easy to implement, is faster and yields results much improved over conventional simulation methods for the materials explored. By direct calculation, we show that the method works for both poor and excellent glass forming materials. It offers a means to adda priori information in first- principles modeling of materials, and represents a significant step toward the computational design of non-crystalline materials using accurate interatomic interactions and experimental information. On the eve of the First World War, William Lawrence Bragg and his father, William Henry Bragg, exposed crys- talline solids to X-rays and discovered what we now call “Bragg diffraction”, strong reflection at particular inci- dent angles and wavelengths. These “Bragg peaks” were sharply defined and, when analyzed with a wave theory of the X-rays, led to clear evidence of order in the crystalline state1. By analyzing the diffraction angles at which the peaks appeared and the wavelength of the X-rays, the full structure of the crystal could be ascertained. In the language of modern solid state physics, the X-ray structure factor of a single crystal consists of a sequence of sharp spikes, which are broadened in a minor way by thermal effects. The information obtained from this palisade of delta functions, arising from a crystal, is sufficient for the determination of lattice structure. The rapid devel- opment of X-ray Crystallography in the past several decades had made it possible to successfully determine the structure of complex protein molecules, with more than 105 atoms, leading to the formation of a new branch of protein crystallography in structural biology2. In contrast with crystals, amorphous materials and liquids have structure factors that are smooth, and thus contain far less specific information about structure. The lack of sharp peaks principally originates from the presence of local atomic ordering in varying length scales, and no long-range order in the amorphous state. The resulting structure factor is one-dimensional and is effectively a sum rule that must be satisfied by the three-dimensional amorphous solids. This presents a far more difficult problem of structural determination of amorphous solids that requires the development of new tools and reasoning. A natural approach to address the problem is to carry out computer simulations, either employing molecular dynamics or Monte Carlo, with suit- able interatomic potentials. We have called this approach the “simulation paradigm”3 elsewhere. By contrast, the other limit is to attempt to invert the diffraction data by “Reverse Monte Carlo” (RMC) or otherwise without using any interatomic potential but information only4,5. This we have called the “information paradigm”3. The information paradigm in its purest form produces models reproducing the data using a random process. These models tend to be maximally disordered and chemically unrealistic. The information paradigm is closely related to the challenge of Materials by Design6,7, for which one imposes external constraints to incorporate additional information on a model to enable a set of preferred physical properties that are of technological utility. Neither paradigm is ideal, or even adequate. The simulation paradigm is plagued by severe size and time-scale limitations that misrepresent the real process of forming a glass, not to mention imperfect interatomic interac- tions. For amorphous materials with no or weak glass-forming ability, either approach is rather desperate, and leads to the formation of unrealistic models with too many structural defects in the networks. In this paper we 1Department of Physics and Astronomy, Condensed Matter Surface Science Program, Ohio University, Athens OH 45701, USA. 2Department of Physics and Astronomy, The University of Southern Mississippi, Hattiesburg MS 39406, USA. 3Department of Physics and Astronomy, Ohio University, Athens OH 45701, USA. Correspondence and requests for materials should be addressed to D.A.D. (email: [email protected]) SCIENTIFIC REPORTS | 6:33731 | DOI: 10.1038/srep33731 1 www.nature.com/scientificreports/ introduce ab initio Force Enhanced Atomic Refinement (AIFEAR). A preliminary trial of the algorithm using only empirical potentials recently appeared8. Others have undertaken related approaches9–17. By including ‘uniformity’ as a constraint for the refinement of models, Goodwin and coworkers showed their Invariant Environment Refinement Technique11 to produce improved models of a-Si and other systems. A liquid-quench procedure, combined with a hybrid Reverse Monte Carlo approach, which incorporates both experimental and energy-based constraints has been employed by Opletal and coworkers12. A similar approach via hybrid RMC with empirical bonded and non-bonded forces was used by Gereben and Pusztai to study liquid dimethyl trisulfide18. Likewise, by refining the initial interatomic empirical potential-energy function and fitting the input experimental structure-factor data, empirical potential structure refinement has been quite successful in predicting liquid structures14. An alternative approach, experi- mentally constrained molecular relaxation, which incorporates experimental information in first-principles mod- eling of materials in a ‘self-consistent’ manner16 was discussed in refs 15 and 16. Recently, a means for including electronic a priori information has also appeared17. These methods have all contributed significantly to the field, yet they have limitations such as employing empirical potentials of limited reliability8,12, or unacceptable con- vergence properties15,16. A general and successful framework for inverting solid state diffraction data does not exist. AIFEAR is a major step toward this important goal. 2 We begin with some definitions. If V(X1 … Xn) is the energy functional for atomic coordinates {Xi} and χ measures the discrepancy between diffraction experiment and theory, we seek to find a set of atomic coordi- 2 nates {Xi} with the property that V = minimum and χ is within experimental error. AIFEAR is a simple iter- ative process consisting of (i) producing a structural model at random (at a sensible but not necessarily exact density, which may not be available), (ii) invoking N accepted moves within conventional RMC19 followed by M conjugate-gradient (CG) steps using ab initio interactions. We then iterate (ii) until convergence. The final results do not depend heavily on the numerical values of N and M, which were chosen to be 1000 and 10, respec- tively, for the present work. For the examples discussed here, we find that significantly fewerab initio force calls are needed for AIFEAR than for an ab initio ‘melt-quench’ simulation. In addition, AIFEAR avoids the problem of relative weighting of V and χ2 in a penalty or target energy functional as in hybrid approaches developed elsewhere12,20. If the density of the material is unknown, it is straightforward to carry out the simulation at zero pressure (with variable cell geometries) in the CG loop, and simply pass the modified supercell vectors back to the RMC loop. To illustrate the efficacy of this new approach, we begin with a persistently vexing problem: the structure of amorphous Si which is particularly difficult because the network is over-constrained21,22 and it is not a glass former. The only methods that yield really satisfactory results are the Wooten-Weaire-Winer (WWW)23 and Activation Relaxation Technique24 methods. Structural and electronic experiments reveal that coordination defects in good quality material have a concentration less than a part in 1000. As such, a satisfactory model should have at most a few percent (or less) defects. Inversion methods like RMC and ab initio melt-quench both produce unsatisfactory models with far too many coordination defects compared to experiments. In this illustra- tion, we employ