1
|
Scheurer M, Dreuw A, Epifanovsky E, Head-Gordon M, Stauch T. Modeling Molecules under Pressure with Gaussian Potentials. J Chem Theory Comput 2021; 17:583-597. [PMID: 33350311 DOI: 10.1021/acs.jctc.0c01212] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
The computational modeling of molecules under high pressure is a growing research area that augments experimental high-pressure chemistry. Here, a new electronic structure method for modeling atoms and molecules under pressure, Gaussians On Surface Tesserae Simulate HYdrostatic Pressure (GOSTSHYP) approach, is introduced. In this method, a set of Gaussian potentials is distributed evenly on the van der Waals surface of the investigated chemical system, leading to a compression of the electron density and the atomic scaffold. Since no parameters other than pressure need to be specified, GOSTSHYP allows straightforward geometry optimizations and ab initio molecular dynamics simulations of chemical systems under pressure for nonexpert users. Calculated energies, bond lengths, and dipole moments under pressure fall within the range of established computational methods for high-pressure chemistry. A Diels-Alder reaction and the cyclotrimerization of acetylene showcase the ability of GOSTSHYP to model pressure-induced chemical reactions. The connection to mechanochemistry is pointed out.
Collapse
Affiliation(s)
- Maximilian Scheurer
- Interdisciplinary Center for Scientific Computing, Heidelberg University, D-69120 Heidelberg, Germany
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Heidelberg University, D-69120 Heidelberg, Germany
| | - Evgeny Epifanovsky
- Q-Chem Inc., 6601 Owens Dr, Suite 105, Pleasanton, California 94588, United States
| | - Martin Head-Gordon
- Pitzer Center for Theoretical Chemistry, University of California, Berkeley, South Dr, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Rd, Berkeley, California 94720, United States
| | - Tim Stauch
- Institute for Physical and Theoretical Chemistry, University of Bremen, Leobener Str. NW2, D-28359 Bremen, Germany.,Bremen Center for Computational Materials Science, University of Bremen, Am Fallturm 1, D-28359 Bremen, Germany.,MAPEX Center for Materials and Processes, University of Bremen, Bibliothekstr. 1, D-28359 Bremen, Germany
| |
Collapse
|
2
|
Hung PK, Ha NTT, Hong NV. Computer simulation of diffusion in silica liquid under temperature and pressure. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2013; 36:60. [PMID: 23793830 DOI: 10.1140/epje/i2013-13060-9] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/02/2012] [Accepted: 05/21/2013] [Indexed: 06/02/2023]
Abstract
We have studied the diffusion mechanism in silica liquid following a new approach where the diffusion rate is estimated via the rate of SiO(x) → SiO(x±1) and the mean square displacement of Si particles per SiO(x) → SiO(x±1). Molecular dynamics simulation has been conducted for a model consisting of 1998 particles over a wide range of temperatures (3000-4500 K) and pressure (from 0 to 25.75 GPa). Our results show that the rate of SiO(x) → SiO(x±1) increases either with increasing the temperature or pressure. Further, we find that SiO(x) → SiO(x±1) is heterogeneously distributed through the network structure of the liquid. In particular, it is concentrated on a small section of Si particles in a low-temperature regime and at ambient pressure. The spatial localisation of SiO(x) → SiO(x±1) originates from the fact that the stable unit in low- and high-pressure regime is SiO4 and SiO6, respectively. The major change in the diffusion mechanism under pressure or temperature concerns the change in the distribution of SiO(x) → SiO(x±1) through the network structure. It is finally shown that the spatial localisation of SiO(x) → SiO(x±1) is responsible for the dynamics heterogeneity and the diffusion anomaly for silica liquid. This finding supports the concept that as the temperature approaches the glass transition point, SiO(x) → SiO(x±1) spatially localises such that the diffusivity drops and the dynamics are anomalously slow.
Collapse
Affiliation(s)
- P K Hung
- Department of Computational Physics, Hanoi University of Science and Technology, 1 Dai Co Viet, Hanoi, Viet Nam.
| | | | | |
Collapse
|
3
|
Williams Q, Hemley RJ, Kruger MB, Jeanloz R. High-pressure infrared sepctra of ∝ -quartz, coesite, stishovite and silica glass. ACTA ACUST UNITED AC 2012. [DOI: 10.1029/93jb02171] [Citation(s) in RCA: 77] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
4
|
Rustad JR, Yuen DA, Spera FJ. The statistical geometry of amorphous silica at lower mantle pressures: Implications for melting slopes of silicates and anharmonicity. ACTA ACUST UNITED AC 2012. [DOI: 10.1029/91jb01925] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
5
|
Karki BB, Stixrude L. First-principles study of enhancement of transport properties of silica melt by water. PHYSICAL REVIEW LETTERS 2010; 104:215901. [PMID: 20867116 DOI: 10.1103/physrevlett.104.215901] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/04/2010] [Indexed: 05/29/2023]
Abstract
First-principles molecular dynamics simulations show that water (8.25 wt%) dramatically affects the transport properties of SiO2 liquid increasing the diffusivity and decreasing the viscosity by an order of magnitude. At 3000 K, the diffusivity of Si, O, and H, and the viscosity vary anomalously with pressure. Highly mobile protons make the hydrous liquid a potential superionic conductor. The predicted dynamical changes are associated with structural depolymerization and water speciation, which changes from being dominated by hydroxyls at low pressure to extended structures at high pressure.
Collapse
Affiliation(s)
- Bijaya B Karki
- Department of Computer Science, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | | |
Collapse
|
6
|
Du J, Corrales LR. Structure, dynamics, and electronic properties of lithium disilicate melt and glass. J Chem Phys 2006; 125:114702. [PMID: 16999496 DOI: 10.1063/1.2345060] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Ab initio molecular dynamics simulations within the framework of density functional theory have been performed to study the structural, dynamic, and electronic properties of lithium disilicate melt and the glass derived from quenching the melt. It is found that lithium ions have a much higher diffusion coefficient and show different diffusion mechanisms than the network forming silicon and oxygen ions in the melt. The simulated lithium disilicate glass structure has 100% four coordinated silicon, close to theoretical nonbridging oxygen to bridging oxygen ratio (2:3), and Q(n) distributions of 20.8%, 58.4%, and 20.8% for n=2,3,4, respectively. In the melt there are considerable amounts (10%-15%) of silicon coordination defects; however, the average silicon coordination number remains about 4, similar to that in the glass. The lithium ion coordination number increases from 3.7 in the glass to 4.4 in the melt mainly due to the increase of bridging oxygen in the first coordination shell. The bond length and bond angle distributions, vibrational density of states, and static structure factors of the simulated glass were determined where the latter was found to be in good agreement with experimental measurement. Atomic charges were obtained based on Bader and Hirshfeld population analyses [Atoms in Molecule: A Quantum Theory (Oxford University Press, Oxford, 1990); Theor. Chim. Acta 44, 129 (1977)]. The average Bader charges found in lithium disilicate glass were -1.729, 3.419, and 0.915 for oxygen, silicon, and lithium, respectively. The corresponding Hirshfeld charges were -0.307, 0.550, and 0.229. The electronic densities of states of the melt and glass were calculated and compared with those of crystalline lithium disilicate.
Collapse
Affiliation(s)
- Jincheng Du
- Chemical Science Division, Pacific Northwest National Laboratory, Richland, Washington 99352, USA.
| | | |
Collapse
|
7
|
Bouhifd MA, Jephcoat AP. Aluminium control of argon solubility in silicate melts under pressure. Nature 2006; 439:961-4. [PMID: 16495996 DOI: 10.1038/nature04583] [Citation(s) in RCA: 35] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2005] [Accepted: 01/12/2006] [Indexed: 11/08/2022]
Abstract
Understanding of the crystal chemistry of the Earth's deep mantle has evolved rapidly recently with the gradual acceptance of the importance of the effect of minor elements such as aluminium on the properties of major phases such as perovskite. In the early Earth, during its formation and segregation into rocky mantle and iron-rich core, it is likely that silicate liquids played a large part in the transport of volatiles to or from the deep interior. The importance of aluminium on solubility mechanisms at high pressure has so far received little attention, even though aluminium has long been recognized as exerting strong control on liquid structures at ambient conditions. Here we present constraints on the solubility of argon in aluminosilicate melt compositions up to 25 GPa and 3,000 K, using a laser-heated diamond-anvil cell. The argon contents reach a maximum that persists to pressures as high as 17 GPa (up to 500 km deep in an early magma ocean), well above that expected on the basis of Al-free melt experiments. A distinct drop in argon solubility observed over a narrow pressure range correlates well with the expected void loss in the melt structure predicted by recent molecular dynamics simulations. These results provide a process for noble gas sequestration in the mantle at various depths in a cooling magma ocean. The concept of shallow partial melting as a unique process for extracting noble gases from the early Earth, thereby defining the initial atmospheric abundance, may therefore be oversimplified.
Collapse
Affiliation(s)
- M Ali Bouhifd
- University of Oxford, Department of Earth Sciences, Parks Road, Oxford OX1 3PR, UK.
| | | |
Collapse
|
8
|
Micoulaut M, Guissani Y, Guillot B. Simulated structural and thermal properties of glassy and liquid germania. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2006; 73:031504. [PMID: 16605529 DOI: 10.1103/physreve.73.031504] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/18/2005] [Revised: 01/10/2006] [Indexed: 05/08/2023]
Abstract
Structural, dynamical, and thermal properties of germanium dioxide are investigated with classical molecular dynamics simulations from the amorphous to the liquid state. Pair correlation functions and coordination numbers are computed under pressure change and show the progressive conversion of the tetrahedral network into an octahedral network, in agreement with experiments. The thermodynamical behavior of the liquid is investigated by means of an equation of state that allows a precise estimation of the compressibility. At low temperature, the diffusion constant D shows an Arrhenius law that progressively deviates when the temperature is increased. The overall comparison with simulated silica permits finally to outline not only the differences in the physical behavior of these two similar systems but also to stress the limitation of the employed germania potential.
Collapse
Affiliation(s)
- M Micoulaut
- Laboratoire de Physique Théorique de la Matière Condensée, CNRS UMR 7600, Université Pierre et Marie Curie, Boite 121, 4, Place Jussieu, 75252 Paris Cedex 05, France
| | | | | |
Collapse
|
9
|
Huang L, Kieffer J. Molecular dynamics study of cristobalite silica using a charge transfer three-body potential: Phase transformation and structural disorder. J Chem Phys 2003. [DOI: 10.1063/1.1529684] [Citation(s) in RCA: 91] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
10
|
Shell MS, Debenedetti PG, Panagiotopoulos AZ. Molecular structural order and anomalies in liquid silica. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2002; 66:011202. [PMID: 12241346 DOI: 10.1103/physreve.66.011202] [Citation(s) in RCA: 89] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/11/2002] [Indexed: 05/23/2023]
Abstract
The present investigation examines the relationship between structural order, diffusivity anomalies, and density anomalies in liquid silica by means of molecular dynamics simulations. We use previously defined orientational and translational order parameters to quantify local structural order in atomic configurations. Extensive simulations are performed at different state points to measure structural order, diffusivity, and thermodynamic properties. It is found that silica shares many trends recently reported for water [J. R. Errington and P. G. Debenedetti, Nature 409, 318 (2001)]. At intermediate densities, the distribution of local orientational order is bimodal. At fixed temperature, order parameter extrema occur upon compression: a maximum in orientational order followed by a minimum in translational order. Unlike water, however, silica's translational order parameter minimum is broad, and there is no range of thermodynamic conditions where both parameters are strictly coupled. Furthermore, the temperature-density regime where both structural order parameters decrease upon isothermal compression (the structurally anomalous regime) does not encompass the region of diffusivity anomalies, as was the case for water.
Collapse
Affiliation(s)
- M Scott Shell
- Department of Chemical Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | | | | |
Collapse
|
11
|
Ehrman SH. Effect of Particle Size on Rate of Coalescence of Silica Nanoparticles. J Colloid Interface Sci 1999; 213:258-261. [PMID: 10191028 DOI: 10.1006/jcis.1999.6105] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
In high temperature processes, at high concentrations, aerosol particles grow by collisions and coalescence. The rate of coalescence is an important parameter for predicting the final primary particle size. For some materials, particularly silica, predictions of final primary particle size using coalescence rates based upon bulk material properties are not in agreement with experimental results. One explanation may be that the high internal pressure in very small particles (less than 10 nm in diameter) affects the mobility of material within the particles and hence the coalescence rate. In this Note, a new approach for estimating rates of coalescence of particles in the initial stages of growth is presented. Coalescence of liquid particles is assumed to be rate-limited by atomic mobility, and the effect of internal pressure on diffusivity is considered. Copyright 1999 Academic Press.
Collapse
Affiliation(s)
- SH Ehrman
- Laboratory for Radio and Environmental Chemistry, Paul Scherrer Institute, Villigen, PSI, CH-5232, Switzerland
| |
Collapse
|
12
|
Diefenbacher J, McMillan PF, Wolf GH. Molecular Dynamics Simulations of Na2Si4O9 Liquid at High Pressure. J Phys Chem B 1998. [DOI: 10.1021/jp973276a] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jason Diefenbacher
- Science & Engineering of Materials Program, Arizona State University, Tempe, Arizona 85287-1704
| | - Paul F. McMillan
- Center for Solid State Science, Arizona State University, Tempe, Arizona 85287-1704, and Department of Chemistry & Biochemistry, Arizona State University, Tempe, Arizona 85287-1604
| | - George H. Wolf
- Department of Chemistry & Biochemistry, Arizona State University, Tempe, Arizona 85287-1604
| |
Collapse
|
13
|
Barrat JL, Badro J, Gillet P. A Strong to Fragile Transition in a Model of Liquid Silica. MOLECULAR SIMULATION 1997. [DOI: 10.1080/08927029708024165] [Citation(s) in RCA: 37] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
14
|
Guillot B, Guissani Y. A Molecular Dynamics Study of the Vibrational Spectra of Silica Polyamorphs. MOLECULAR SIMULATION 1997. [DOI: 10.1080/08927029708024167] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
15
|
|
16
|
Venuti E. High-pressure densification of silica glass: A molecular-dynamics simulation. PHYSICAL REVIEW. B, CONDENSED MATTER 1996; 54:3809-3816. [PMID: 9986279 DOI: 10.1103/physrevb.54.3809] [Citation(s) in RCA: 36] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/12/2023]
|
17
|
Sarnthein J, Pasquarello A, Car R. Model of vitreous SiO2 generated by an ab initio molecular-dynamics quench from the melt. PHYSICAL REVIEW. B, CONDENSED MATTER 1995; 52:12690-12695. [PMID: 9980432 DOI: 10.1103/physrevb.52.12690] [Citation(s) in RCA: 79] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/12/2023]
|
18
|
Tsuneyuki S, Matsui Y. Molecular dynamics study of pressure enhancement of ion mobilities in liquid silica. PHYSICAL REVIEW LETTERS 1995; 74:3197-3200. [PMID: 10058136 DOI: 10.1103/physrevlett.74.3197] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/23/2023]
|
19
|
Tse JS, Klug DD. High-pressure densification of amorphous silica. PHYSICAL REVIEW. B, CONDENSED MATTER 1992; 46:5933-5938. [PMID: 10002276 DOI: 10.1103/physrevb.46.5933] [Citation(s) in RCA: 67] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/12/2023]
|
20
|
Meade C, Hemley RJ, Mao HK. High-pressure x-ray diffraction of SiO2 glass. PHYSICAL REVIEW LETTERS 1992; 69:1387-1390. [PMID: 10047204 DOI: 10.1103/physrevlett.69.1387] [Citation(s) in RCA: 145] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/23/2023]
|
21
|
Della Valle RG, Andersen HC. Molecular dynamics simulation of silica liquid and glass. J Chem Phys 1992. [DOI: 10.1063/1.463056] [Citation(s) in RCA: 95] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
22
|
Farber DL, Williams Q. Pressure-Induced Coordination Changes in Alkali-Germanate Melts: An in Situ Spectroscopic Investigation. Science 1992; 256:1427-30. [PMID: 17791610 DOI: 10.1126/science.256.5062.1427] [Citation(s) in RCA: 21] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2022]
Abstract
The structure of liquid Na(2)Ge(2)O(5).H(2)O, a silicate melt analog, has been studied with Raman spectroscopy to pressures of 2.2 gigapascals. Upon compression, a peak near approximately 240 wavenumbers associated with octahedral GeO(6) groups grows relative to a peak near approximately 500 wavenumbers associated with tetrahedral GeO(4) groups. This change corresponds to an increase in octahedral germanium in the liquid from near 0% at ambient pressures to >50% at a pressure of 2.2 gigapascals. Silicate liquids plausibly undergo similar coordination changes at depth in the Earth. Such structural changes may generate decreases in the fusion slopes of silicates at high pressures as well as neutrally buoyant magmas within the transition zone of the Earth's mantle.
Collapse
|
23
|
Tse JS, Klug DD. The structure and dynamics of silica polymorphs using a two‐body effective potential model. J Chem Phys 1991. [DOI: 10.1063/1.461198] [Citation(s) in RCA: 74] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
24
|
Rustad JR, Yuen DA, Spera FJ. Molecular dynamics of amorphous silica at very high pressures (135 GPa): Thermodynamics and extraction of structures through analysis of Voronoi polyhedra. PHYSICAL REVIEW. B, CONDENSED MATTER 1991; 44:2108-2121. [PMID: 9999759 DOI: 10.1103/physrevb.44.2108] [Citation(s) in RCA: 20] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/12/2023]
|
25
|
Stebbins JF. NMR evidence for five-coordinated silicon in a silicate glass at atmospheric pressure. Nature 1991. [DOI: 10.1038/351638a0] [Citation(s) in RCA: 155] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|