The liquid structure of tetrachloroethene: Molecular dynamics simulations and reverse Monte Carlo modeling with interatomic potentials

Orsolya Gereben, L. Pusztai

Research output: Contribution to journalArticle

2 Citations (Scopus)

Abstract

The liquid structure of tetrachloroethene has been investigated on the basis of measured neutron and X-ray scattering structure factors, applying molecular dynamics simulations and reverse Monte Carlo (RMC) modeling with flexible molecules and interatomic potentials. As no complete all-atom force field parameter set could be found for this planar molecule, the closest matching all-atom Optimized Potentials for Liquid Simulations (OPLS-AA) intra-molecular parameter set was improved by equilibrium bond length and angle parameters coming from electron diffraction experiments [I. L. Karle and J. Karle, J. Chem. Phys. 20, 63 (1952)]. In addition, four different intra-molecular charge distribution sets were tried, so in total, eight different molecular dynamics simulations were performed. The best parameter set was selected by calculating the mean square difference between the calculated total structure factors and the corresponding experimental data. The best parameter set proved to be the one that uses the electron diffraction based intra-molecular parameters and the charges qC = 0.1 and q Cl = -0.05. The structure was further successfully refined by applying RMC computer modeling with flexible molecules that were kept together by interatomic potentials. Correlation functions concerning the orientation of molecular axes and planes were also determined. They reveal that the molecules closest to each other exclusively prefer the parallel orientation of both the molecular axes and planes. Molecules forming the first maximum of the center-center distribution have a preference for 60°axis orientation and >60°molecular plane arrangement. A second coordination sphere at ∼11 Å and a very small third one at ∼16 Å can be found as well, without preference for any axis or plane orientation.

Original languageEnglish
Article number164509
JournalThe Journal of Chemical Physics
Volume139
Issue number16
DOIs
Publication statusPublished - Oct 28 2013

Fingerprint

Tetrachloroethylene
Molecular dynamics
molecular dynamics
Molecules
Computer simulation
Liquids
liquids
simulation
Electron diffraction
molecules
electron diffraction
Atoms
Charge distribution
Bond length
X ray scattering
Neutrons
charge distribution
field theory (physics)
atoms
neutrons

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry

Cite this

@article{30fdceaaabae450e98ece8b765af14e1,
title = "The liquid structure of tetrachloroethene: Molecular dynamics simulations and reverse Monte Carlo modeling with interatomic potentials",
abstract = "The liquid structure of tetrachloroethene has been investigated on the basis of measured neutron and X-ray scattering structure factors, applying molecular dynamics simulations and reverse Monte Carlo (RMC) modeling with flexible molecules and interatomic potentials. As no complete all-atom force field parameter set could be found for this planar molecule, the closest matching all-atom Optimized Potentials for Liquid Simulations (OPLS-AA) intra-molecular parameter set was improved by equilibrium bond length and angle parameters coming from electron diffraction experiments [I. L. Karle and J. Karle, J. Chem. Phys. 20, 63 (1952)]. In addition, four different intra-molecular charge distribution sets were tried, so in total, eight different molecular dynamics simulations were performed. The best parameter set was selected by calculating the mean square difference between the calculated total structure factors and the corresponding experimental data. The best parameter set proved to be the one that uses the electron diffraction based intra-molecular parameters and the charges qC = 0.1 and q Cl = -0.05. The structure was further successfully refined by applying RMC computer modeling with flexible molecules that were kept together by interatomic potentials. Correlation functions concerning the orientation of molecular axes and planes were also determined. They reveal that the molecules closest to each other exclusively prefer the parallel orientation of both the molecular axes and planes. Molecules forming the first maximum of the center-center distribution have a preference for 60°axis orientation and >60°molecular plane arrangement. A second coordination sphere at ∼11 {\AA} and a very small third one at ∼16 {\AA} can be found as well, without preference for any axis or plane orientation.",
author = "Orsolya Gereben and L. Pusztai",
year = "2013",
month = "10",
day = "28",
doi = "10.1063/1.4826320",
language = "English",
volume = "139",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "16",

}

TY - JOUR

T1 - The liquid structure of tetrachloroethene

T2 - Molecular dynamics simulations and reverse Monte Carlo modeling with interatomic potentials

AU - Gereben, Orsolya

AU - Pusztai, L.

PY - 2013/10/28

Y1 - 2013/10/28

N2 - The liquid structure of tetrachloroethene has been investigated on the basis of measured neutron and X-ray scattering structure factors, applying molecular dynamics simulations and reverse Monte Carlo (RMC) modeling with flexible molecules and interatomic potentials. As no complete all-atom force field parameter set could be found for this planar molecule, the closest matching all-atom Optimized Potentials for Liquid Simulations (OPLS-AA) intra-molecular parameter set was improved by equilibrium bond length and angle parameters coming from electron diffraction experiments [I. L. Karle and J. Karle, J. Chem. Phys. 20, 63 (1952)]. In addition, four different intra-molecular charge distribution sets were tried, so in total, eight different molecular dynamics simulations were performed. The best parameter set was selected by calculating the mean square difference between the calculated total structure factors and the corresponding experimental data. The best parameter set proved to be the one that uses the electron diffraction based intra-molecular parameters and the charges qC = 0.1 and q Cl = -0.05. The structure was further successfully refined by applying RMC computer modeling with flexible molecules that were kept together by interatomic potentials. Correlation functions concerning the orientation of molecular axes and planes were also determined. They reveal that the molecules closest to each other exclusively prefer the parallel orientation of both the molecular axes and planes. Molecules forming the first maximum of the center-center distribution have a preference for 60°axis orientation and >60°molecular plane arrangement. A second coordination sphere at ∼11 Å and a very small third one at ∼16 Å can be found as well, without preference for any axis or plane orientation.

AB - The liquid structure of tetrachloroethene has been investigated on the basis of measured neutron and X-ray scattering structure factors, applying molecular dynamics simulations and reverse Monte Carlo (RMC) modeling with flexible molecules and interatomic potentials. As no complete all-atom force field parameter set could be found for this planar molecule, the closest matching all-atom Optimized Potentials for Liquid Simulations (OPLS-AA) intra-molecular parameter set was improved by equilibrium bond length and angle parameters coming from electron diffraction experiments [I. L. Karle and J. Karle, J. Chem. Phys. 20, 63 (1952)]. In addition, four different intra-molecular charge distribution sets were tried, so in total, eight different molecular dynamics simulations were performed. The best parameter set was selected by calculating the mean square difference between the calculated total structure factors and the corresponding experimental data. The best parameter set proved to be the one that uses the electron diffraction based intra-molecular parameters and the charges qC = 0.1 and q Cl = -0.05. The structure was further successfully refined by applying RMC computer modeling with flexible molecules that were kept together by interatomic potentials. Correlation functions concerning the orientation of molecular axes and planes were also determined. They reveal that the molecules closest to each other exclusively prefer the parallel orientation of both the molecular axes and planes. Molecules forming the first maximum of the center-center distribution have a preference for 60°axis orientation and >60°molecular plane arrangement. A second coordination sphere at ∼11 Å and a very small third one at ∼16 Å can be found as well, without preference for any axis or plane orientation.

UR - http://www.scopus.com/inward/record.url?scp=84903365540&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84903365540&partnerID=8YFLogxK

U2 - 10.1063/1.4826320

DO - 10.1063/1.4826320

M3 - Article

AN - SCOPUS:84903365540

VL - 139

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 16

M1 - 164509

ER -