Five-site Model Research Articles

Comparison of Eight Classical Lennard-Jones-Based H2 Molecular Models in the Gas Phase at Temperatures and Pressures Relevant to Hydrogen On-Board Storage Tanks.

This work compares eight classical H2 molecular models in the gas phase taken from the existing literature. All models are based on Lennard-Jones (LJ) 12-6 terms for the van der Waals interactions and hence easier to transfer to multiphase molecular simulations than more sophisticated potentials. The H2 potentials tested include one-site, two-site, three-site, and five-site models, with the sites being either the H atoms, the center-of-mass of the H2 molecule, or massless sites. For the multisite models, high-frequency H-H stretching modes can lead to poor equipartition of the kinetic energy, and the timestep for molecular dynamics (MD) simulations should be reduced to maintain a stable numerical integration of the equations of motion. As such, only those models with rigid bonds are considered. In the present case, 600 MD simulations of H2 gas were carried out over a large range of temperatures (-50 to +90 °C) and at densities corresponding to a pressure range of 50 to 2000 bar, which include the operating conditions of on-board storage tanks in hydrogen-fueled vehicles. Most of the models under study were found to reproduce reasonably well the experimental pVT phase diagram as well as the solubility. Discrepancies only became significant at the highest densities tested, and these could be used to rank the different models. All model diffusion coefficients were essentially indistinguishable from experimental results, and as such, kinetically dominated dynamic properties could not be used as a criterion for the choice of model. Among the eight models tested, two of them, i.e., the two-site model of Yang and Zhong and the one-site model derived from Buch performed very well over the range of conditions tested. They represent a good compromise between realism, simplicity, and computational efficiency.

Five-Site Model for Brownian Dynamics Simulations of a Molecular Walker in Three Dimensions.

Motor proteins play an important role in many biological processes and have inspired the development of synthetic analogues. Molecular walkers, such as kinesin, dynein, and myosin V, fulfill a diverse set of functions including transporting cargo along tracks, pulling molecules through membranes, and deforming fibers. The complexity of molecular motors and their environment makes it difficult to model the detailed dynamics of molecular walkers over long time scales. In this work, we present a simple, three-dimensional model for a molecular walker on a bead-spring substrate. The walker is represented by five spherically symmetric particles that interact through common intermolecular potentials and can be simulated efficiently in Brownian dynamics simulations. The movement of motor protein walkers entails energy conversion through ATP hydrolysis while artificial motors typically rely on a local conversion of energy supplied through external fields. We model energy conversion through rate equations for mechanochemical states that couple positional and chemical degrees of freedom and determine the walker conformation through interaction potential parameters. We perform Brownian dynamics simulations for two scenarios: In the first, the model walker transports cargo by walking on a substrate whose ends are fixed. In the second, a tethered motor pulls a mobile substrate chain against a variable force. We measure relative displacements and determine the effects of cargo size and retarding force on the efficiency of the walker. We find that, while the efficiency of our model walker is less than for the biological system, our simulations reproduce trends observed in single-molecule experiments on kinesin. In addition, the model and simulation method presented here can be readily adapted to biological and synthetic systems with multiple walkers.

Hybrid In Silico and TR-FRET-Guided Discovery of Novel BCL-2 Inhibitors.

B-Cell lymphoma 2 (BCL-2) regulates cell death in humans. In this study, combined multiscale in silico approaches and in vitro studies were employed. A small-molecule library that includes more than 210 000 compounds was used. The predicted therapeutic activity value (TAV) of the compounds in this library was computed with the binary cancer quantitative structure-activity relationships (QSAR) model. The molecules with a high calculated TAV were used in 26 individual toxicity QSAR models. As a result of this screening protocol, 288 nontoxic molecules with high predicted TAV were identified. These selected hits were then screened against the BCL-2 target protein using hybrid docking and molecular dynamics (MD) simulations. The interaction energies of identified compounds were compared with two known BCL-2 inhibitors. Then, the short MD simulations were carried out by initiating the best docking poses of 288 molecules. Average MM/GBSA energies were computed, and long MD simulations were employed to selected hits. The same calculations were also applied for two known BCL-2 inhibitors. Moreover, a five-site (AHRRR) structure-based pharmacophore model was constructed, and this model was used in the screening of the same database. On the basis of hybrid data-driven ligand identification study, final hits were selected and used in in vitro studies. Based on results of the time-resolved fluorescence resonance energy transfer (TR-FRET) analysis, further filtration was carried out for the U87-MG cell line tests. MTT cell proliferation assay analysis results showed that selected three potent compounds were significantly effective on glioma cells.

Contrasting spectroscopic response of human hemoglobin in presence of graphene oxides and its reduced form: Comparative approach with carbon quantum dots.

Recently, a considerable amount of research is being directed towards study of graphene oxide (GO) and its reduced form (RGO) since their exposed functional groups make them better candidates in nanobiotechnolgy. In order to assess their biocompatibility, the nature of interactions between Human Hemoglobin (HHb) and GO/RGO are monitored since a comparative spectroscopic approach towards understanding their nature of interactions has not been investigated previously. UV-vis spectroscopy reveals hyperchromicity for HHb-GO system and hypochromicity for HHb-RGO system in the region of absorption of tryptophan/tyrosine residues. Notably, although steady-state fluorescence static quenching of HHb for GO and enhancement of fluorescence for RGO is noticed, but average fluorescence-lifetime is remaining unchanged in presence of GO/RGO. Calorimetric data illustrates three-site and five-site binding model to be the best-fit model for GO and RGO respectively. Also, synchronous fluorescence quenching corresponding to alterations in microenvironment of tryptophan/ tyrosine residues is observed only in presence of GO. Likewise FTIR spectroscopy elucidates involvement of both amide I and amide II bond of HHb backbone through H-bonding interaction only for GO. Furthermore RLS spectra demonstrate an increase and a decrease in signal for GO and RGO respectively. Surprisingly, secondary structure of HHb is maintained upon interaction with both GO/RGO, as revealed by CD spectroscopy, thus supporting their potential application in biological microenvironment. Thus it appears that the spectroscopic properties of HHb upon interaction with GO is altered upon its reduction to RGO. Furthermore the role of HHb as good candidate for bimolecular interaction has been highlighted.

Structural features of high-local-density water molecules: Insights from structure indicators based on the translational order between the first two molecular shells.

The two-liquids scenario for liquid water assumes the existence of two competing preferential local molecular structural states characterized by either low or high local density. While the former is expected to present good local order thus involving privileged structures, the latter is usually regarded as conforming a high-entropy unstructured state. A main difference in the local arrangement of such "classes" of water molecules can be inferred from the degree of translational order between the first and second molecular shells. This is so, since the low-local-density molecules present a clear gap between the first two shells while in the case of the high-local-density ones, one or more molecules from the second shell have collapsed toward the first one, thus populating the intershell region. Some structural indicators, like the widely employed local structure index and the recently introduced ζ index, have been devised precisely on the basis of this observation, being successful in detecting well-structured low-local-density molecules. However, the nature of the high-local-density state has been mainly disregarded over the years. In this work we employ molecular dynamics simulations for two water models (the extended simple point charge model and the five-site model) at the liquid and supercooled regimes combined with the inherent dynamics approach (energy minimizations of the instantaneous configurations) in order to both rationalize the detailed structural and topological information that these indicators provide and to advance in our understanding of the high-density state.

Investigating the Effects of Linker Extension on H2 Sorption in the rht-Metal–Organic Framework NU-111 by Molecular Simulations

Grand canonical Monte Carlo (GCMC) simulations were performed to investigate hydrogen sorption in the rht-metal–organic framework (MOF) NU-111, which is composed of Cu2+ ions coordinated to 1,3,5-tris[(1,3-carboxylate-5-(4-(ethynyl)phenyl))butadiynyl]-benzene linkers. This MOF was shown to have a high Brunauer–Emmett–Teller surface area of 5000 ± 80 m2 g–1 and a hydrogen uptake of 20.8 mg g–1 at 77 K and 1 bar. It has a similar structure to that of PCN-61, which was studied previously by our group [Forrest, K. A.; et al. J. Phys. Chem. C 2012, 116, 15538−15549]; however, NU-111 has an extra triple bond on each arm in the linker body. The purpose of this study was to investigate the effects of linker extension on the MOF–sorbate interactions using PCN-61 as a basis for comparison. Three hydrogen models of varying complexities were used in the simulations: a single-site model, which only includes van der Waals interactions, a five-site model, which includes van der Waals and permanent electrostatic interact...

Force Field Benchmark of Amino Acids: I. Hydration and Diffusion in Different Water Models.

Thermodynamic and kinetic properties are of critical importance for the applicability of computational models to biomolecules such as proteins. Here we present an extensive evaluation of the Amber ff99SB-ILDN force field for modeling of hydration and diffusion of amino acids with three-site (SPC, SPC/E, SPC/Eb, and TIP3P), four-site (TIP4P, TIP4P-Ew, and TIP4P/2005), and five-site (TIP5P and TIP5P-Ew) water models. Hydration free energies (HFEs) of neutral amino acid side chain analogues have little dependence on the water model, with a root-mean-square error (RMSE) of ∼1 kcal/mol from experimental observations. On the basis of the number of interacting sites in the water model, HFEs of charged side chains can be putatively classified into three groups, of which the group of three-site models lies between those of four- and five-site water models; for each group, the water model dependence is greatly eliminated when the solvent Galvani potential is considered. Some discrepancies in the location of the first hydration peak ( RRDF) in the ion-water radial distribution function between experimental and calculated observations were detected, such as a systematic underestimation of the acetate (Asp side chain) ion. The RMSE of calculated diffusion coefficients of amino acids from experiment increases linearly with the increasing diffusion coefficients of the solvent water models at infinite dilution. TIP3P has the fastest diffusivity, in line with literature findings, while the "FB" and "OPC" water model families as well as TIP4P/2005 perform well, within a relative error of 5%, and TIP4P/2005 yields the most accurate estimate for the water diffusion coefficient. All of the tested water models overestimate amino acid diffusion coefficients by approximately 40% (TIP4P/2005) to 200% (TIP3P). Scaling of protein-water interactions with TIP4P/2005 in the Amber ff99SBws and ff03ws force fields leads to more negative HFEs but has little influence on the diffusion of amino acids. The most recent FF/water combinations of ff14SB/OPC3, ff15ipq/SPC/Eb, and fb15/TIP3P-FB do not show obvious improvements in accuracy for the tested quantities. These findings here establish a benchmark that may aid in the development and improvement of classical force fields to accurately model protein dynamics and thermodynamics.

Building better water models using the shape of the charge distribution of a water molecule.

The unique properties of liquid water apparently arise from more than just the tetrahedral bond angle between the nuclei of a water molecule since simple three-site models of water are poor at mimicking these properties in computer simulations. Four- and five-site models add partial charges on dummy sites and are better at modeling these properties, which suggests that the shape of charge distribution is important. Since a multipole expansion of the electrostatic potential describes a charge distribution in an orthogonal basis set that is exact in the limit of infinite order, multipoles may be an even better way to model the charge distribution. In particular, molecular multipoles up to the octupole centered on the oxygen appear to describe the electrostatic potential from electronic structure calculations better than four- and five-site models, and molecular multipole models give better agreement with the temperature and pressure dependence of many liquid state properties of water while retaining the computational efficiency of three-site models. Here, the influence of the shape of the molecular charge distribution on liquid state properties is examined by correlating multipoles of non-polarizable water models with their liquid state properties in computer simulations. This will aid in the development of accurate water models for classical simulations as well as in determining the accuracy needed in quantum mechanical/molecular mechanical studies and ab initio molecular dynamics simulations of water. More fundamentally, this will lead to a greater understanding of how the charge distribution of a water molecule leads to the unique properties of liquid water. In particular, these studies indicate that p-orbital charge out of the molecular plane is important.

Redfield Treatment of Multipathway Electron Transfer in Artificial Photosynthetic Systems.

Coherence effects on electron transfer in a series of symmetric and asymmetric two-, three-, four-, and five-site molecular model systems for photosystem I in cyanobacteria and green plants were studied. The total site energies of the electronic Hamiltonian were calculated using the density functional theory (DFT) formalism and included the zero point vibrational energies of the electron donors and acceptors. Site energies and couplings were calculated using a polarizable continuum model to represent various solvent environments, and the site-to-site couplings were calculated using fragment charge difference methods at the DFT level of theory. The Redfield formalism was used to propagate the electron density from the donors to the acceptors, incorporating relaxation and dephasing effects to describe the electron transfer processes. Changing the relative energies of the donor, intermediate acceptor, and final acceptor molecules in these assemblies has profound effects on the electron transfer rates as well as on the amplitude of the quantum oscillations observed. Increasing the ratio of a particular energy gap to the electronic coupling for a given pair of states leads to weaker quantum oscillations between sites. Biasing the intermediate acceptor energies to slightly favor one pathway leads to a general decrease in electron transfer yield.

A single-site multipole model for liquid water.

Accurate and efficient empirical potential energy models that describe the atomistic interactions between water molecules in the liquid phase are essential for computer simulations of many problems in physics, chemistry, and biology, especially when long length or time scales are important. However, while models with non-polarizable partial charges at four or five sites in a water molecule give remarkably good values for certain properties, deficiencies have been noted in other properties and increasing the number of sites decreases computational efficiency. An alternate approach is to utilize a multipole expansion of the electrostatic potential due to the molecular charge distribution, which is exact outside the charge distribution in the limits of infinite distances or infinite orders of multipoles while partial charges are a qualitative representation of electron density as point charges. Here, a single-site multipole model of water is presented, which is as fast computationally as three-site models but is also more accurate than four- and five-site models. The dipole, quadrupole, and octupole moments are from quantum mechanical-molecular mechanical calculations so that they account for the average polarization in the liquid phase, and represent both the in-plane and out-of-plane electrostatic potentials of a water molecule in the liquid phase. This model gives accurate thermodynamic, dynamic, and dielectric properties at 298 K and 1 atm, as well as good temperature and pressure dependence of these properties.

Molecular simulations and density functional theory calculations of bromine in clathrate hydrate phases.

Bromine forms a tetragonal clathrate hydrate structure (TS-I) very rarely observed in clathrate hydrates of other guest substances. The detailed structure, energetics, and dynamics of Br2 and Cl2 in TS-I and cubic structure I (CS-I) clathrate hydrates are studied in this work using molecular dynamics and quantum chemical calculations. X-ray diffraction studies show that the halogen-water-oxygen distances in the cages of these structures are shorter than the sum of the van der Waals radii of halogen and oxygen atoms. This suggests that the stabilizing effects of halogen bonding or other non-covalent interactions (NCIs) may contribute to the formation of the unique tetragonal bromine hydrate structure. We performed molecular dynamics simulations of Br2 and Cl2 clathrate hydrates using our previously developed five-site charge models for the dihalogen molecules [Dureckova et al. Can. J. Chem. 93, 864 (2015)] which reproduce the computed electrostatic potentials of the dihalogens and account for the electropositive σ-hole of the halogen bond donor (the dihalogen). Analysis of the radial distribution functions, enthalpies of encapsulation, velocity and orientation autocorrelation functions, and polar angle distributions are carried out for Br2 and Cl2 guests in various cages to contrast the properties of these guests in the TS-I and CS-I phases. Quantum chemical partial geometry optimizations of Br2 and Cl2 guests in the hydrate cages using the M06-2X functional give short halogen-water distances compatible with values observed in X-ray diffraction experiments. NCI plots of guest-cage structures are generated to qualitatively show the relative strength of the non-bonding interactions between dihalogens and water molecules. The differences between behaviors of Br2 and Cl2 guests in the hydrate cages may explain why bromine forms the unique TS-I phase.

Structural and configurational properties of nanoconfined monolayer ice from first principles.

Understanding the structural tendencies of nanoconfined water is of great interest for nanoscience and biology, where nano/micro-sized objects may be separated by very few layers of water. Here we investigate the properties of ice confined to a quasi-2D monolayer by a featureless, chemically neutral potential, in order to characterize its intrinsic behaviour. We use density-functional theory simulations with a non-local van der Waals density functional. An ab initio random structure search reveals all the energetically competitive monolayer configurations to belong to only two of the previously-identified families, characterized by a square or honeycomb hydrogen-bonding network, respectively. We discuss the modified ice rules needed for each network, and propose a simple point dipole 2D lattice model that successfully explains the energetics of the square configurations. All identified stable phases for both networks are found to be non-polar (but with a topologically non-trivial texture for the square) and, hence, non-ferroelectric, in contrast to previous predictions from a five-site empirical force-field model. Our results are in good agreement with very recently reported experimental observations.

Transport properties of carbon dioxide and methane from molecular dynamics simulations.

Transport properties of carbon dioxide and methane are predicted for temperatures between (273.15 and 573.15) K and pressures up to 800 MPa by molecular dynamics simulations. Viscosities and thermal conductivities were obtained through the Green-Kubo formalism, whereas the Einstein relation was used to provide self-diffusion coefficient estimates. The differences in property predictions due to the force field nature and parametrization were investigated by the comparison of seven different CO2 models (two single-site models, three rigid three-site models, and two fully flexible three-site models) and three different CH4 models (two single-site models and one fully flexible five-site model). The simulation results show good agreement with experimental data, except for thermal conductivities at low densities. The molecular structure and force field parameters play an important role in the accuracy of the simulations, which is within the experimental deviations reported for viscosities and self-diffusion coefficients considering the most accurate CO2 and CH4 models studied. On the other hand, the molecular flexibility does not seem to improve accuracy, since the explicit account of vibrational and bending degrees of freedom in the CO2 flexible models leads to slightly less accurate results. Nonetheless, the use of a correctional term to account for vibrational modes in rigid models generally improves estimations of thermal conductivity values. At extreme densities, the caging effect observed with single-site representations of the molecules restrains mobility and leads to an unphysical overestimation of viscosities and, conversely, to the underestimation of self-diffusion coefficients. This result may help to better understand the limits of applicability of such force fields concerning structural and transport properties of dense systems.

Molecular Dynamics Study on Configuration Energy and Radial Distribution Functions of Ammonium Dihydrogen Phosphates Solution

Molecular dynamics simulations were carried out to study the configuration energy and radial distribution functions of mmonium dihydrogen phosphate solution at different temperatures. The dihydrogen phosphate ion was treated as a seven-site model and the ammonium ion was regarded as a five-site model, while a simple-point-charge model for water molecule. An unusually local particle number density fluctuation was observed in the system at saturation temperature. It can be found that the potential energy increases slowly with the temperature from 373 K to 404 K, which indicates that the ammonium dihydrogen phosphate has partly decomposed. The radial distribution function between the hydrogen atom of ammonium cation and the oxygen atom of dihydrogen phosphate ion at three different temperatures shows obvious difference, which indicates that the average H-bond number changes obviously with the temperature. The temperature has an influence on the combination between hydrogen atoms and phosphorus atoms of dihydrogen phosphate ion and there are much more growth units at saturated solutions.

A Polarizable and Transferable PHAST N2 Potential for Use in Materials Simulation.

A polarizable and transferable intermolecular potential energy function, potentials with high accuracy, speed, and transferability (PHAST), has been developed from first principles for molecular nitrogen to be used in the modeling of heterogeneous processes such as materials sorption and separations. A five-site (van der Waals and point charge) anisotropic model, that includes many-body polarization, is proposed. It is parametrized to reproduce high-level electronic structure calculations (CCSD(T) using Dunning-type basis sets extrapolated to the CBS limit) for a representative set of dimer potential energy curves. Thus it provides a relatively simple yet robust and broadly applicable representation of nitrogen. Two versions are developed, differing by the type of mixing rules applied to unlike Lennard-Jones potential sites. It is shown that the Waldman-Hagler mixing rules are more accurate than Lorentz-Berthelot. The resulting potentials are demonstrated to be effective in modeling neat nitrogen but are designed to also be useful in modeling N2 interactions in a large array of environments such as metal-organic frameworks and zeolites and at interfaces. In such settings, capturing anisotropic forces and interactions with (open and coordinated) metals and charged/polar environments is essential. In developing the potential, it was found that adding a seemingly redundant dimer orientation, slip-parallel (S), improved the transferability of the potential energy surface (PES). Notably, one of the solid phases of nitrogen was not as accurately represented energetically without including S in the representative set. Liquid simulations, however, were unaffected and worked equally well for both potentials. This suggests that accounting for a wide variety of configurations is critical in designing a potential that is intended for use in heterogeneous environments where many orientations, including those not commonly explored in the bulk, are possible. Testing and validation of the potential are achieved via simulations of a thermal distribution of trimer geometries compared to analogous high level electronic structure calculations and molecular simulations of bulk pressure-density isotherms across the vapor, supercritical, and liquid phases. Crystal lattice parameters and energetics of the α-N2 and γ-N2 solid phases are also evaluated and determined to be in good agreement with experiment. Thus the proposed potential is shown to be efficacious for gas, liquid, and solid use, representing both disordered and ordered configurations.

Modeling and Simulation of Ethylene Polymerization in Industrial Slurry Reactor Series

A five-site comprehensive mathematical model was developed to simulate the steady-state behavior of industrial slurry polymerization of ethylene in multistage continuous stirred tank reactors. More specifically, the effects of various operating conditions (i.e., inflow rates of catalyst, hydrogen and comonomer) on the molecular structure and properties of polyethylene (i.e., Mw, Mn, polydispersity index (IPD), melt index, density, etc.) are fully assessed. It is shown that the proposed comprehensive model is capable of simulating the steady-state operation of an industrial slurry stirred tank reactor series. It is demonstrated that changing the catalyst flow rate, changes simultaneously the mean residence-time in both reactors, which plays a significant role on the establishment of polyethylene architecture properties such as molecular mass and IPD. The melt index and density of polyethylene are mainly controlled by hydrogen and comonomer concentration, respectively.

Thermodynamic analysis of Bacillus subtilis endospore protonation using isothermal titration calorimetry

Bacterial proton and metal adsorption reactions have the capacity to affect metal speciation and transport in aqueous environments. We coupled potentiometric titration and isothermal titration calorimetry (ITC) analyses to study Bacillus subtilis spore-proton adsorption. We modeled the potentiometric data using a four and five-site non-electrostatic surface complexation model (NE-SCM). Heats of spore surface protonation from coupled ITC analyses were used to determine site specific enthalpies of protonation based on NE-SCMs. The five-site model resulted in a substantially better model fit for the heats of protonation but did not significantly improve the potentiometric titration model fit. The improvement observed in the five-site protonation heat model suggests the presence of a highly exothermic protonation reaction circa pH 7 that cannot be resolved in the less sensitive potentiometric data. From the log Ks and enthalpies we calculated corresponding site specific entropies. Log Ks and site concentrations describing spore surface protonation are statistically equivalent to B. subtilis cell surface protonation constants. Spore surface protonation enthalpies, however, are more exothermic relative to cell based adsorption suggesting a different bonding environment. The thermodynamic parameters defined in this study provide insight on molecular scale spore-surface protonation reactions. Coupled ITC and potentiometric titrations can reveal highly exothermic, and possibly endothermic, adsorption reactions that are overshadowed in potentiometric models alone. Spore-proton adsorption NE-SCMs derived in this study provide a framework for future metal adsorption studies.

3D-QSAR and pharmacophore model study on aryl diphenolic azoles as estrogen receptor-β ligands

Pharmacophore and three-dimensional quantitative structure–activity relationship (3D-QSAR) model of 106 aryl diphenolic azoles as estrogen receptor-β ligands were proposed. The best pharmacophore model, five-site model including two hydrogen bond donors and three aromatic rings, has been established. The biological activity values (pIC50) and classification (active or inactive) of compounds were predicted by 3D-QSAR model. The model correlation coefficient (R 2) is 0.9623. The classification accuracies in predicting for training, test, and total datasets are 97.5, 84.6, and 94.3 %, respectively. The 3D-QSAR model maps revealed that hydrogen bond donor, hydrophobic, and electron-withdrawing fields would be the essential factors for designing new ligands with higher biological activity. The results indicate that the combination of pharmacophore and 3D-QSAR analyses could be a powerful modeling tool for finding novel estrogen receptor-β ligands.

Adsorption on Edge-Functionalized Bilayer Graphene Nanoribbons: Assessing the Role of Functional Groups in Methane Uptake

A combination of ab initio and classical Monte Carlo simulations is used to investigate the effects of functional groups on methane binding. Using Møller–Plesset (MP2) calculations, we obtain the binding energies for benzene functionalized with NH2, OH, CH3, COOH, and H2PO3 and identify the methane binding sites. In all cases, the preferred binding sites are located above the benzene plane in the vicinity of the benzene carbon atom attached to the functional group. Functional groups enhance methane binding relative to benzene (−6.39 kJ/mol), with the largest enhancement observed for H2PO3 (−8.37 kJ/mol) followed by COOH and CH3 (−7.77 kJ/mol). Adsorption isotherms are obtained for edge-functionalized bilayer graphene nanoribbons using grand canonical Monte Carlo simulations with a five-site methane model. Adsorbed excess and heats of adsorption for pressures up to 40 bar and 298 K are obtained with functional group concentrations ranging from 3.125 to 6.25 mol % for graphene edges functionalized with OH, NH2, and COOH. The functional groups are found to act as preferred adsorption sites, and in the case of COOH the local methane density in the vicinity of the functional group is found to exceed that of bare graphene. The largest enhancement of 44.5% in the methane excess adsorbed is observed for COOH-functionalized nanoribbons when compared to H terminated ribbons. The corresponding enhancements for OH- and NH2-functionalized ribbons are 10.5% and 3.7%, respectively. The excess adsorption across functional groups reflects the trends observed in the binding energies from MP2 calculations. Our study reveals that specific site functionalization can have a significant effect on the local adsorption characteristics and can be used as a design strategy to tailor materials with enhanced methane storage capacity.

Prediction of the Temperature-Dependent Thermal Conductivity and Shear Viscosity for Rigid Water Models

The temperature-dependent thermal conductivity and shear viscosity of liquid water between 283 and 363 K are evaluated for eight rigid models with reverse nonequilibrium molecular dynamics (RNEMD). In comparison with experimental data, five-site models (TIP5P and TIP5P-Ew) have apparent advantages in estimating thermal conductivities than other rigid water models that overestimate the value by tens of percent. For shear viscosity, no single model can reproduce all experimental data; instead, five- and four-site models show their own strength in a certain temperature range. Meanwhile, all of the current rigid models obtain lower values than experimental data when temperature is lower than 298 K, while the TIP5P and TIP5P-Ew models can relatively accurately predict the values over others at a temperature range from 298 to 318 K. At a higher temperature range shear viscosity of liquid water can be reproduced with a four-site model (TIP4P-2005, TIP4P-Ew) fairly well.