Ion-Gated Gas Separation through Porous Graphene

Similar documents
Universal Repulsive Contribution to the. Solvent-Induced Interaction Between Sizable, Curved Hydrophobes: Supporting Information

Supporting Information for Solid-liquid Thermal Transport and its Relationship with Wettability and the Interfacial Liquid Structure

Computational Predictions of 1-Octanol/Water Partition Coefficient for Imidazolium based Ionic Liquids.

WATER PERMEATION THROUGH GRAPHENE NANOSLIT BY MOLECULAR DYNAMICS SIMULATION

Supporting Information

Structure-Property Relationships of Porous Materials for Carbon Dioxide Separation and Capture

Supporting Information

Supporting Information for: Physics Behind the Water Transport through. Nanoporous Graphene and Boron Nitride

k θ (θ θ 0 ) 2 angles r i j r i j

Diffusion of Water and Diatomic Oxygen in Poly(3-hexylthiophene) Melt: A Molecular Dynamics Simulation Study

Supporting Information

Unusual Entropy of Adsorbed Methane on Zeolite Templated Carbon. Supporting Information. Part 2: Statistical Mechanical Model

Opening the Gate: Framework Flexibility in ZIF-8 Explored by Experiments and Simulations

STRUCTURE OF IONS AND WATER AROUND A POLYELECTROLYTE IN A POLARIZABLE NANOPORE

Introduction to molecular dynamics

Ionic Liquids simulations : obtention of structural and transport properties from molecular dynamics. C. J. F. Solano, D. Beljonne, R.

Supporting Information

Non-equilibrium molecular dynamics simulation study of the behavior of hydrocarbon-isomers in silicalite

SUPPLEMENTARY INFORMATION

Determination of Kamlet-Taft parameters for selected solvate ionic liquids.

Effect of polarizability of halide anions on the ionic salvation in water clusters

INVESTIGATION OF THE ABSORPTION OF CO 2 IN IONIC LIQUID. Kalyan Dhar 1 * and Syed Fahim 1

Supporting Information

Schwarzites for Natural Gas Storage: A Grand- Canonical Monte Carlo Study

Theoretical comparative study on hydrogen storage of BC 3 and carbon nanotubes

Multiscale Coarse-Graining of Ionic Liquids

All-atom Molecular Mechanics. Trent E. Balius AMS 535 / CHE /27/2010

PORE SIZE DISTRIBUTION OF CARBON WITH DIFFERENT PROBE MOLECULES

Unit Cell-Level Thickness Control of Single-Crystalline Zinc Oxide Nanosheets Enabled by Electrical Double Layer Confinement

Molecular Dynamics Simulation of a Nanoconfined Water Film

The change in free energy on transferring an ion from a medium of low dielectric constantε1 to one of high dielectric constant ε2:

Dynamics. capacities of. (A.I.1) (E b ), bond. E total stretching (A.I.2) (A.I.3) A.I.1. Ebond force (A.I.4) (A.I.5)

Molecular dynamics study on the microstructure of CH3COOLi solutions with different concentrations

Phase Equilibria of binary mixtures by Molecular Simulation and PR-EOS: Methane + Xenon and Xenon + Ethane

MOLECULAR DYNAMICS SIMULATIONS OF CH4 CLATHRATE HYDRATE DISSOCIATION ADJACENT TO HYDRATED SILICA SURFACES

Supporting Information

Dioxide Is Facilitated In Narrow Carbon. Nanopores

First principles molecular dynamics simulation of a taskspecific. ionic liquid based on silver-olefin complex:

Shear Properties and Wrinkling Behaviors of Finite Sized Graphene

Competitive I 2 Sorption by Cu-BTC from Humid Gas Streams

Supplementary Information for: Controlling Cellular Uptake of Nanoparticles with ph-sensitive Polymers

Glass-Transition and Side-Chain Dynamics in Thin Films: Explaining. Dissimilar Free Surface Effects for Polystyrene and Poly(methyl methacrylate)

Promoting effect of ethanol on dewetting transition in the confined region of melittin tetramer

MD simulation of methane in nanochannels

Agency, Honcho, Kawaguchi, Saitama (Japan), University, Tsushima, Kita-ku, Okayama (Japan),

Supporting Information for. Hydrogen Bonding Structure at Zwitterionic. Lipid/Water Interface

Molecular Dynamics Investigation of the ω-current in the Kv1.2 Voltage Sensor Domains

Molecular Dynamics Study of Carbon Dioxide Storage in Carbon-Based Organic Nanopores

Supplementary Material. Physisorption of Hydrophobic and Hydrophilic 1-alkyl-3- methylimidazolium Ionic Liquids on the Graphite Plate Surface

Molecular dynamics simulation of limiting conductances for LiCl, NaBr, and CsBr in supercritical water

New Six-site Acetonitrile Model for Simulations of Liquid Acetonitrile and its Aqueous Mixtures

Supporting Information for. Ab Initio Metadynamics Study of VO + 2 /VO2+ Redox Reaction Mechanism at the Graphite. Edge Water Interface

Electronegativity is a very useful concept for the explanation or understanding of chemical reactivity throughout the periodic table.

Why Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar.

CARBON 2004 Providence, Rhode Island. Adsorption of Flexible n-butane and n-hexane on Graphitized Thermal Carbon Black and in Slit Pores

SUPPLEMENTARY INFORMATION

Supporting Information

Lecture 11: Potential Energy Functions

arxiv: v1 [cond-mat.soft] 26 Nov 2018

Electronic Supporting Information Topological design of porous organic molecules

Thomas Roussel, Roland J.-M. Pellenq, Christophe Bichara. CRMC-N CNRS, Campus de Luminy, Marseille, cedex 09, France. Abstract.

Explanation of Dramatic ph-dependence of Hydrogen Binding on Noble Metal Electrode: Greatly Weakened Water Adsorption at High ph.

Hyeyoung Shin a, Tod A. Pascal ab, William A. Goddard III abc*, and Hyungjun Kim a* Korea

ATOMISTIC MODELLING OF CROSSLINKED EPOXY POLYMER

Supplementary Figure 1 Irregular arrangement of E,E-8-mer on TMA. STM height images formed when

Developing Monovalent Ion Parameters for the Optimal Point Charge (OPC) Water Model. John Dood Hope College

Tailoring the Properties of Quadruplex Nucleobases for Biological and Nanomaterial Applications

Comparative Study of the Water Response to External Force at Nanoscale and Mesoscale

Supplemntary Infomation: The nanostructure of. a lithium glyme solvate ionic liquid at electrified. interfaces

Atomic and molecular interaction forces in biology

ONETEP PB/SA: Application to G-Quadruplex DNA Stability. Danny Cole

Supporting Information: Improved Parametrization. of Lithium, Sodium, Potassium, and Magnesium ions. for All-Atom Molecular Dynamics Simulations of

Mechanical Properties of Tetra-Polyethylene and Tetra-Polyethylene Oxide Diamond Networks via Molecular Dynamics Simulations

Supplementary Information for Atomistic Simulation of Spinodal Phase Separation Preceding Polymer Crystallization

Supporting information for: Anomalous Stability of Two-Dimensional Ice. Confined in Hydrophobic Nanopore

Molecular Dynamics Simulations. Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia

Peptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia

MARTINI simulation details

The micro-properties of [hmpy+] [Tf 2 N-] Ionic liquid: a simulation. study. 1. Introduction

arxiv: v2 [cond-mat.stat-mech] 23 Sep 2009

Structural transition and solid-like behavior of alkane films confined in nano-spacing

Supporting Information Soft Nanoparticles: Nano Ionic Networks of Associated Ionic Polymers

Proceedings of the ASME 2009 International Mechanical Engineering Congress & Exposition

Supplementary Information. Surface Microstructure Engenders Unusual Hydrophobicity in. Phyllosilicates

Preparation of Model

c 2010 ANJAN V. RAGHUNATHAN

Atomistic Modeling of Cross-linked Epoxy Polymer

CE 530 Molecular Simulation

arxiv: v1 [physics.chem-ph] 8 Mar 2010

What Does an Ionic Liquid Surface Really Look Like? Unprecedented Details from Molecular Simulations

An Informal AMBER Small Molecule Force Field :

Supporting Information

Morphology of Supported Polymer Electrolyte. Ultra-thin Films: a Numerical Study. Supporting Information

Site dependent hydrogenation in Graphynes: A Fully Atomistic Molecular Dynamics Investigation

Supplementary Figure 1. Schematic of rapid thermal annealing process: (a) indicates schematics and SEM cross-section of the initial layer-by-layer

Interface Resistance and Thermal Transport in Nano-Confined Liquids


Amino Acids and Proteins at ZnO-water Interfaces in Molecular Dynamics Simulations: Electronic Supplementary Information

SUPPLEMENTARY INFORMATION

Supporting information

Transcription:

Online Supporting Information for: Ion-Gated Gas Separation through Porous Graphene Ziqi Tian, Shannon M. Mahurin, Sheng Dai,*,, and De-en Jiang *, Department of Chemistry, University of California, Riverside, California 92521, United States Chemical Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States Department of Chemistry, The University of Tennessee, Knoxville, Tennessee 37996, United States *Corresponding authors: dais@ornl.gov; djiang@ucr.edu Table of Contents 1. Computational details 2. Pore-size modulation by the anion as reflected by the pore-anion distance 3. Gas permeation through the IL/graphene composite membranes of different pore sizes 4. Gas permeation through the [emim]bf 4 monolayer without the graphene layer 5. Effect of ionic liquid thickness on gas permeation 6. Force field parameters 7. References 1. Computational details Classical molecular dynamics (CMD) simulations were performed with the LAMMPS package 1 in the canonical (NVT) ensemble with two-dimensional periodic boundary condition (PBC) in xy directions. The porous graphene was fixed in the box, with dimensions of 49.2 Å by 42.6 Å in x and y directions, respectively. Forty ionic liquid pairs were coated above the surface in the case of a monolayer. Fifty gas molecules were placed in the gas phase above the ionic liquid layer on the feed side of the bichamber setup. The initial height of the feed-side chamber was adjusted to make the initial gas pressure at 10 atm, determined by the Peng Robinson equation. A vacuum layer was added on the other side of the porous graphene in the bi-chamber setup as the permeate side, with a thickness of 90 Å. During simulation, the size of the simulation box was fixed. The non-polarizable OPLS-AA type force field 2 was employed for the ionic liquids. Lennard- Jones parameters for the [emim] cation and all parameters for the BF 4 - and PF 6 - anions were from Liu et al. 3 The parameters for the bonded terms of the [emim] cation were from Borodin. 4 The partial atomic charges for the [emim] cation were obtained by fitting to the electrostatic potential at the B3LYP/6-31g(d) level from the RESP method of the Merz Kollman scheme. This set of parameters for [emim][bf 4 ] and [emim][pf 6 ] could well reproduce their bulk dynamic properties, as shown in our previous work. 5 The atomic charges of porous graphene were obtained based on DFT-derived electrostatic potential, also known as the Repeating Electrostatic Potential Extracted Atomic charges (REPEAT) method. 6 Regarding the Lennard-Jones terms for the carbon atoms in the graphene layer, we chose 0.086 kcal/mol for ε and 3.4 Å for σ as recommended, 7-8 while ε of 0.015 kcal/mol and σ of 2.45 Å were used for the terminating hydrogen atoms around the graphene pore rim. For CO 2 and N 2, three-site models were adopted, 9, 10 as our previous simulation. 11 All-atom model was used for methane. 2 Lennard-Jones potential terms were evaluated via the Lorentz-Berthelot mixing rule with a cutoff of 12 Å. S1

To calculate the long-range electrostatic interaction in the 2D periodic simulation cell, we used the 3D Particle-Particle-Particle-Mesh (PPPM) method with a slab correction for our main simulation runs. In the z-direction, a fictitious empty volume was inserted between the 2D slabs to correct the interslab dipolar interactions. Since the PPPM method is time-consuming for the periodic 2D system, 12 we also used a much faster truncation method for our parallel runs to obtain statistics on gas permeance, whereby a cutoff of 15 Å was used for the electrostatic interactions. Our tests showed that these two methods give very similar gas permeation data for the same system. All the data reported in this work were averaged over 20 parallel simulations from various initial velocity distributions. For each simulation, at the beginning, the ionic liquid layer was heated up to 1000 K for 1 ns and quenched to 300 K in 1 ns, to disperse the ionic liquid on porous graphene evenly. In most cases, ionic liquid formed a uniform layer, covering all the nanopores. After ionic liquid coating, 25 ns simulation was carried out in NVT ensemble. Temperature of ionic liquid and gas were kept at 300 K 13, 14 with the Nose-Hoover algorithm. 2. Pore-size modulation by the anion as reflected by the pore-anion distance Figure S1. Probability distribution of the distance from the pore center to the center of the closest BF 4 - anion for the [emim][bf 4 ]/6.0-Å-porous-graphene system (Figure 1c in the main text): with and without the CO 2 gas (w/o gas). S2

3. Gas permeation through the IL/graphene composite membranes of different pore sizes Figure S2. (a) A larger graphene pore of 9.6 Å in size; (b) a smaller graphene pore of 4.2 Å in size; (c) permeation of CO 2 and CH 4 through the [emim][bf 4 ]/9.6-Å-porous-graphene system. 4. Gas permeation through the [emim][bf 4 ] monolayer without the graphene layer Figure S3. Comparison of CO 2 and CH 4 permeation through a hypothetical unsupported ionic liquid layer at the same simulation conditions used in Figure 4 in the text. The centers of mass for both cations and anions are fixed to prevent the collapse of the membrane. S3

5. Effect of ionic liquid thickness on gas permeation Figure S4. Permeation of CO 2 and CH 4 through the [emim][bf 4 ]/6.0-Å-porous-graphene system for different thickness of the [emim][bf 4 ] ionic liquid (IL) layer. 6. Force field parameters 6.1 Gas molecules Table S1. Partial atomic charges and Lennard-Jones parameters for gas molecules ε / kcal/mol σ / Å q / e CO 2 C 0.0559 2.757 0.6512 O 0.1600 2.565-0.3256 N 2 N 0.0728 3.318-0.4048 Center of Mass 0 0 0.8096 CH 4 C 0.0664 3.500-0.2400 H 0.0300 2.500 0.0600 6.2 Ionic liquids The OPLS force field was used for the ionic liquids The non-bonded term includes the Lennard-Jones and the Columbic interaction; the bonded term is expressed below: E bonded = K b (R R 0 ) 2 + K θ (θ θ 0 ) 2 + K φ [1 + d cos(nφ)] 2 bonds angles dihedrals impropers S4

Scheme 1. Partial atomic charges on the [emim], BF 4, and PF 6 ions. Atoms are labelled for the bonded-term parameters. Table S2. Lennard-Jones parameters for the ionic liquids ε / kcal/mol σ / Å Emim C (sp2; on the ring) 0.0860 3.400 C (sp3; off the ring) 0.0110 3.400 N 0.1700 3.250 H (-CH 3 ) 0.0157 2.500 H (-CH 2 -) 0.0157 2.650 H (on the ring) 0.0150 2.450 Anions B 0.0950 3.581 P 0.2000 3.742 F 0.0610 3.118 Bonded terms for the emim, BF 4, and PF 6 ion (K R, K θ, Kφ all in kcal/mol) Bond K R R 0 / Å C-Cm 309.0 1.520 H-Cm 327.5 1.100 C-N 369.5 1.472 Cm-N 369.5 1.472 H-C 327.5 1.100 Cc-N 650.0 1.372 Cc-Cc 309.0 1.353 H-Cc 327.5 1.083 B-F 290.0 1.389 P-F 190.0 1.600 Angle K θ θ 0 / H-Cm-C 43.0 110.5 H-Cm-H 38.5 107.7 Cm-C-N 70.0 109.0 H-C-N 35.0 104.6 C-N-Cc 70.0 124.3 S5

Cc-N-Cc 70.0 109.6 Cc-Cc-N 72.0 106.7 H-Cc-N 35.0 122.7 H-Cc-Cc 36.0 130.6 N-Cc-N 70.0 111.4 H-Cm-N 35.0 108.3 F-B-F 49.9 109.5 F-P-F 80.0 90.0/180.0 Dihedral Kφ d n N-Cc-N-Cc 10.000-1 2 N-Cc-N-Cm 2.000-1 2 H-Cc-N-Cc 1.500-1 2 H-Cc-N-Cm 1.500-1 2 Cc-Cc-N-Cc 10.000-1 2 Cc-Cc-N-Cm 2.000-1 2 N-Cc-Cc-H 1.500-1 2 N-Cc-Cc-N 10.000-1 2 H-Cc-Cc-H 1.500-1 2 N-C-Cm-H 0.160 1 3 H-C-Cm-H 0.150 1 3 H-Cm-N-Cc 0.164 1 3 Improper Kφ d n N-N-Cc-H 1.1-1 2 Cc-N-Cc-H 1.1-1 2 Cc-Cc-N-C 2.0-1 2 6.3 Porous graphene Table S3. Lennard-Jones parameters for C and H atoms on the porous graphene ε / kcal/mol σ / Å C 0.0860 3.400 H 0.0150 2.450 Cartesian coordinates (Å) and partial atomic charges on the 6.0-Å porous graphene Rectangular unit cell: a=24.60 Å, b=21.30 Å q / e x y H 0.1480 9.490680 8.388491 H 0.1490 13.475880 7.446969 H 0.1490 16.285200 10.420645 H 0.1490 8.314800 12.299428 H 0.1490 11.124120 15.275234 H 0.1490 15.109320 14.333712 H 0.1490 11.124120 7.446969 H 0.1480 15.109320 8.388491 H 0.1490 8.314800 10.420645 H 0.1490 9.490680 14.333712 H 0.1490 13.475880 15.275234 H 0.1490 16.285200 12.299428 S6

C 0.0000 0.000000-0.002130 C -0.0010 1.230000 2.128010 C 0.0000 0.000000 1.420803 C 0.0010 1.230000 3.548813 C 0.0000 2.462460 0.000000 C -0.0020 3.692460 2.125880 C 0.0000 2.462460 1.420803 C 0.0020 3.690000 3.548813 C -0.0010 4.922460 0.000000 C -0.0020 6.154920 2.130140 C 0.0000 4.922460 1.420803 C 0.0080 6.154920 3.553074 C -0.0030 7.384920-0.002130 C 0.0050 8.612460 2.134400 C -0.0020 7.384920 1.422934 C -0.0060 8.614920 3.557334 C 0.0000 9.842460 0.000000 C -0.0060 11.070000 2.125880 C 0.0010 9.842460 1.420803 C 0.0150 11.070000 3.548813 C 0.0030 12.300000 0.000000 C -0.0060 13.530000 2.125880 C 0.0030 12.300000 1.420803 C 0.0150 13.530000 3.548813 C 0.0000 14.757540 0.000000 C 0.0050 15.987540 2.134400 C 0.0010 14.757540 1.420803 C -0.0060 15.985080 3.557334 C -0.0030 17.215080-0.002130 C -0.0030 18.445080 2.130140 C -0.0020 17.215080 1.422934 C 0.0080 18.445080 3.553074 C -0.0010 19.677540 0.000000 C -0.0020 20.907540 2.125880 C 0.0000 19.677540 1.420803 C 0.0020 20.910000 3.548813 C 0.0000 22.137540 0.000000 C -0.0010 23.370000 2.128010 C 0.0000 22.137540 1.420803 C 0.0010 23.370000 3.548813 C -0.0010 0.000000 4.256020 C -0.0010 1.227540 6.386160 C 0.0000 0.000000 5.678953 C 0.0010 1.227540 7.809093 C -0.0010 2.460000 4.256020 C 0.0020 3.685080 6.384030 C 0.0000 2.457540 5.676823 C -0.0020 3.680160 7.809093 C -0.0010 4.917540 4.258150 C 0.0120 6.145080 6.386160 C -0.0050 4.912620 5.674693 S7

C -0.0170 6.132780 7.800573 C -0.0080 7.382460 4.264540 C 0.0780 8.607540 6.415982 C -0.0170 7.380000 5.687474 C -0.3010 8.570640 7.804833 C -0.0180 9.844920 4.264540 C -0.3030 11.079840 6.358468 C 0.0810 9.864600 5.695994 C -0.0230 12.300000 4.245369 C -0.3030 13.520160 6.358468 C 0.1910 12.300000 5.668303 C -0.0180 14.755080 4.264540 C 0.0780 15.992460 6.415982 C 0.0810 14.735400 5.695994 C -0.3010 16.029360 7.804833 C -0.0080 17.217540 4.264540 C 0.0120 18.454920 6.386160 C -0.0170 17.220000 5.687474 C -0.0180 18.467220 7.800573 C -0.0010 19.682460 4.258150 C 0.0020 20.914920 6.384030 C -0.0060 19.687380 5.674693 C -0.0020 20.919840 7.809093 C -0.0010 22.140000 4.256020 C 0.0000 23.372460 6.386160 C 0.0000 22.142460 5.676823 C 0.0010 23.372460 7.809093 C -0.0010 0.000000 8.516300 C -0.0020 1.230000 10.650700 C 0.0000 0.000000 9.943494 C -0.0010 1.230000 12.071503 C -0.0020 2.457540 8.516300 C -0.0080 3.692460 10.644310 C 0.0080 2.464920 9.941363 C -0.0070 3.692460 12.077894 C 0.0080 4.912620 8.516300 C 0.0800 6.172140 10.633659 C -0.0150 4.920000 9.937103 C 0.0800 6.172140 12.088545 C 0.1890 7.365240 8.509909 C -0.3020 7.348020 9.915802 C 0.1890 17.234760 8.509909 C 0.0800 18.427860 10.633659 C -0.3020 17.251980 9.915802 C 0.0800 18.427860 12.088545 C 0.0080 19.687380 8.516300 C -0.0080 20.907540 10.644310 C -0.0150 19.680000 9.937103 C -0.0070 20.907540 12.077894 C -0.0020 22.142460 8.516300 C -0.0020 23.370000 10.650700 S8

C 0.0080 22.135080 9.941363 C -0.0010 23.370000 12.071503 C -0.0010 0.000000 12.778710 C 0.0010 1.227540 14.913110 C 0.0000 0.000000 14.205904 C -0.0030 1.227540 16.333914 C 0.0070 2.464920 12.780840 C -0.0040 3.680160 14.913110 C 0.0000 2.457540 14.203774 C 0.0020 3.685080 16.338174 C -0.0180 4.920000 12.785100 C -0.0230 6.132780 14.919501 C 0.0130 4.912620 14.203774 C 0.0150 6.145080 16.336044 C -0.3020 7.348020 12.804272 C -0.3040 8.570640 14.917370 C 0.1930 7.365240 14.212294 C 0.0800 8.607540 16.304092 C -0.3020 11.079840 16.363735 C -0.3020 13.520160 16.363735 C -0.3040 16.029360 14.917370 C 0.0800 15.992460 16.304092 C -0.3020 17.251980 12.804272 C -0.0230 18.467220 14.919501 C 0.1930 17.234760 14.212294 C 0.0150 18.454920 16.336044 C -0.0180 19.680000 12.785100 C -0.0040 20.919840 14.913110 C 0.0130 19.687380 14.203774 C 0.0020 20.914920 16.338174 C 0.0070 22.135080 12.780840 C 0.0010 23.372460 14.913110 C 0.0000 22.142460 14.203774 C -0.0030 23.372460 16.333914 C 0.0030 0.000000 17.041120 C -0.0010 1.230000 19.173390 C 0.0000 0.000000 18.466184 C 0.0000 1.230000 20.594194 C 0.0020 2.457540 17.045380 C -0.0010 3.690000 19.171260 C 0.0000 2.460000 18.464054 C 0.0000 3.692460 20.594194 C -0.0060 4.912620 17.045380 C 0.0060 6.154920 19.167000 C 0.0010 4.917540 18.464054 C -0.0010 6.154920 20.592063 C -0.0200 7.380000 17.032599 C -0.0100 8.614920 19.164870 C -0.0040 7.382460 18.457663 C 0.0080 8.612460 20.587803 C 0.0790 9.864600 17.026209 S9

C 0.0130 11.070000 19.173390 C -0.0160 9.844920 18.457663 C -0.0050 11.070000 20.596324 C 0.1900 12.300000 17.051771 C 0.0130 13.530000 19.173390 C -0.0220 12.300000 18.476834 C -0.0050 13.530000 20.596324 C 0.0790 14.735400 17.026209 C -0.0100 15.985080 19.164870 C -0.0160 14.755080 18.457663 C 0.0080 15.987540 20.587803 C -0.0200 17.220000 17.032599 C 0.0060 18.445080 19.167000 C -0.0040 17.217540 18.457663 C -0.0010 18.445080 20.592063 C -0.0060 19.687380 17.045380 C -0.0010 20.910000 19.171260 C 0.0010 19.682460 18.464054 C 0.0000 20.907540 20.594194 C 0.0020 22.142460 17.045380 C -0.0010 23.370000 19.173390 C 0.0000 22.140000 18.464054 C 0.0000 23.370000 20.594194 7. References (1) Plimpton, S. J. Comput. Phys. 1995, 117, 1-19. (2) Jorgensen, W. L.; Maxwell, D. S.; TiradoRives, J. J. Am. Chem. Soc. 1996, 118, 11225-11236. (3) Liu, Z. P.; Huang, S. P.; Wang, W. C. J. Phys. Chem. B 2004, 108, 12978-12989. (4) Borodin, O. J. Phys. Chem. B 2009, 113, 11463-11478. (5) Babarao, R.; Dai, S.; Jiang, D. E. J. Phys. Chem. B 2011, 115, 9789-9794. (6) Campana, C.; Mussard, B.; Woo, T. K. J. Chem. Theory Comput. 2009, 5, 2866-2878. (7) Sidorenkov, A. V.; Kolesnikova, S. V.; Saletsky, A. M. Eur Phys J B 2016, 89, 220. (8) Wang, J. M.; Cieplak, P.; Kollman, P. A. J. Comput. Chem. 2000, 21, 1049-1074. (9) Murthy, C. S.; Singer, K.; Klein, M. L.; Mcdonald, I. R. Mol. Phys. 1980, 41, 1387-1399. (10) Harris, J. G.; Yung, K. H. J. Phys. Chem. 1995, 99, 12021-12024. (11) Liu, H. J.; Dai, S.; Jiang, D. E. Nanoscale 2013, 5, 9984-9987. (12) Yeh, I. C.; Berkowitz, M. L. J. Chem. Phys. 1999, 111, 3155-3162. (13) Nose, S. J. Chem. Phys. 1984, 81, 511-519. (14) Hoover, W. G. Phys. Rev. A 1985, 31, 1695-1697. S10