1
|
Gavazzoni C, Skaf MS. Adsorption of CO 2 and CH 4 in MIL-47 investigated by the 3D-RISM molecular theory of solvation. Phys Chem Chem Phys 2020; 22:13240-13247. [PMID: 32500908 DOI: 10.1039/d0cp01025a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Metal-organic frameworks (MOFs) comprise a class of highly porous nanomaterials formed by the assembly of organic molecular templates connected by metal ions. These materials exhibit a large diversity of pore size and geometry, topology, surface area, and chemical functionality. MOFs are particularly promising materials for developing new technologies for capture and storage of greenhouse gases such as methane and carbon dioxide. Here we apply the three dimensional reference interaction site model (3D-RISM) molecular theory of solvation to study the interactions of CO2 and CH4 with the metal-organic material MIL-47. The 3D-RISM integral equations were solved to determine the three dimensional density correlation functions of the gas (solvent) relative to the atomic positions of the MIL-47 framework, treated as static solute sites. The computed solvent spatial distributions inside MIL-47 pores were used to identify whether or not there exist preferable binding sites and the binding free energy landscape for the gas of interest at low computational costs compared with other molecular modeling techniques, such as grand canonical Monte Carlo and molecular dynamics simulations. The 3D-RISM formalism was applied to pure CO2, pure CH4, and binary mixtures of these gases of various compositions under different pressure conditions. The results indicate that both gases bind very weakly to MIL-47 and that this material exhibits nearly vanishing CO2/CH4 selectivity. The 3D-RISM computations presented here can be extended to investigate the physical adsorption of gases on other MOFs and nanoporous materials, providing an alternative low-cost computational approach to study gas capture and storage in nanoporous materials in general and, in particular, to determine the binding free-energy landscape in these systems.
Collapse
Affiliation(s)
- Cristina Gavazzoni
- Institute of Chemistry and Center for Computing in Engineering and Sciences, University of Campinas - Unicamp, Campinas, SP 13082-864, Brazil.
| | | |
Collapse
|
2
|
Munaò G, Saija F. Integral equation study of soft-repulsive dimeric fluids. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:115101. [PMID: 28155850 DOI: 10.1088/1361-648x/aa5306] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
We study fluid structure and water-like anomalies of a system constituted by dimeric particles interacting via a purely repulsive core-softened potential by means of integral equation theories. In our model, dimers interact through a repulsive pair potential of inverse-power form with a softened repulsion strength. By employing the Ornstein-Zernike approach and the reference interaction site model (RISM) theory, we study the behavior of water-like anomalies upon progressively increasing the elongation λ of the dimers from the monomeric case ([Formula: see text]) to the tangent configuration ([Formula: see text]). For each value of the elongation we consider two different values of the interaction potential, corresponding to one and two length scales, with the aim to provide a comprehensive description of the possible fluid scenarios of this model. Our theoretical results are systematically compared with already existing or newly generated Monte Carlo data: we find that theories and simulations agree in providing the picture of a fluid exhibiting density and structural anomalies for low values of λ and for both the two values of the interaction potential. Integral equation theories give accurate predictions for pressure and radial distribution functions, whereas the temperatures where anomalies occur are underestimated. Upon increasing the elongation, the RISM theory still predicts the existence of anomalies; the latter are no longer observed in simulations, since their development is likely precluded by the onset of crystallization. We discuss our results in terms of the reliability of integral equation theories in predicting the existence of water-like anomalies in core-softened fluids.
Collapse
Affiliation(s)
- Gianmarco Munaò
- Dipartimento di Scienze Matematiche e Informatiche, Scienze Fisiche e Scienze della Terra, Università degli Studi di Messina, Viale F. Stagno d'Alcontres 31, 98166 Messina, Italy
| | | |
Collapse
|
3
|
Borówko M, Rżysko W, Sokołowski S, Staszewski T. Integral equations theory for two-dimensional systems involving nanoparticles. Mol Phys 2016. [DOI: 10.1080/00268976.2016.1234079] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Małgorzata Borówko
- Department for the Modeling of Physico-Chemical Processes, Maria Curie-Skłodowska University, Lublin, Poland
| | - Wojciech Rżysko
- Department for the Modeling of Physico-Chemical Processes, Maria Curie-Skłodowska University, Lublin, Poland
| | - Stefan Sokołowski
- Department for the Modeling of Physico-Chemical Processes, Maria Curie-Skłodowska University, Lublin, Poland
| | - Tomasz Staszewski
- Department for the Modeling of Physico-Chemical Processes, Maria Curie-Skłodowska University, Lublin, Poland
| |
Collapse
|
4
|
Omelyan I, Kovalenko A. MTS-MD of Biomolecules Steered with 3D-RISM-KH Mean Solvation Forces Accelerated with Generalized Solvation Force Extrapolation. J Chem Theory Comput 2016; 11:1875-95. [PMID: 26574393 DOI: 10.1021/ct5010438] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
We developed a generalized solvation force extrapolation (GSFE) approach to speed up multiple time step molecular dynamics (MTS-MD) of biomolecules steered with mean solvation forces obtained from the 3D-RISM-KH molecular theory of solvation (three-dimensional reference interaction site model with the Kovalenko-Hirata closure). GSFE is based on a set of techniques including the non-Eckart-like transformation of coordinate space separately for each solute atom, extension of the force-coordinate pair basis set followed by selection of the best subset, balancing the normal equations by modified least-squares minimization of deviations, and incremental increase of outer time step in motion integration. Mean solvation forces acting on the biomolecule atoms in conformations at successive inner time steps are extrapolated using a relatively small number of best (closest) solute atomic coordinates and corresponding mean solvation forces obtained at previous outer time steps by converging the 3D-RISM-KH integral equations. The MTS-MD evolution steered with GSFE of 3D-RISM-KH mean solvation forces is efficiently stabilized with our optimized isokinetic Nosé-Hoover chain (OIN) thermostat. We validated the hybrid MTS-MD/OIN/GSFE/3D-RISM-KH integrator on solvated organic and biomolecules of different stiffness and complexity: asphaltene dimer in toluene solvent, hydrated alanine dipeptide, miniprotein 1L2Y, and protein G. The GSFE accuracy and the OIN efficiency allowed us to enlarge outer time steps up to huge values of 1-4 ps while accurately reproducing conformational properties. Quasidynamics steered with 3D-RISM-KH mean solvation forces achieves time scale compression of conformational changes coupled with solvent exchange, resulting in further significant acceleration of protein conformational sampling with respect to real time dynamics. Overall, this provided a 50- to 1000-fold effective speedup of conformational sampling for these systems, compared to conventional MD with explicit solvent. We have been able to fold the miniprotein from a fully denatured, extended state in about 60 ns of quasidynamics steered with 3D-RISM-KH mean solvation forces, compared to the average physical folding time of 4-9 μs observed in experiment.
Collapse
Affiliation(s)
- Igor Omelyan
- Department of Mechanical Engineering, University of Alberta , Mechanical Engineering Building 4-9, Edmonton, Alberta T6G 2G8, Canada.,National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada.,Institute for Condensed Matter Physics, National Academy of Sciences of Ukraine , 1 Svientsitskii Street, Lviv 79011, Ukraine
| | - Andriy Kovalenko
- Department of Mechanical Engineering, University of Alberta , Mechanical Engineering Building 4-9, Edmonton, Alberta T6G 2G8, Canada.,National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
| |
Collapse
|
5
|
Casanova D, Gusarov S, Kovalenko A, Ziegler T. Evaluation of the SCF Combination of KS-DFT and 3D-RISM-KH; Solvation Effect on Conformational Equilibria, Tautomerization Energies, and Activation Barriers. J Chem Theory Comput 2015; 3:458-76. [PMID: 26637028 DOI: 10.1021/ct6001785] [Citation(s) in RCA: 54] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The effect of solvation on conformational equilibria, tautomerization energies, and activation barriers in simple SN2 reactions is reproduced by using the self-consistent field coupling of the Kohn-Sham density functional theory (KS-DFT) for electronic structure and the three-dimensional reference interaction site model with the closure approximation of Kovalenko and Hirata (3D-RISM-KH) for molecular solvation structure. These examples are used in order to validate the implementation of the 3D-RISM-KH method in the Amsterdam Density Functional package. The computations of the free energy difference in the trans/gauche conformational equilibrium for 1,2-dichloroethane in different solvents; the relative tautomerization free energy for cytosine, isocytosine, and guanine; and the free energy activation barrier for a CH3X-type (X = F, Cl, Br) SN2 reaction exhibit agreement with the experimental data. The method is also applied to the electronic and hydration structure of carbon single-wall nanotubes of different diameters, including the effect of water located in the inner space of the nanotubes. A comparison with continuum models of solvation (including COSMO) as well as with other more precise and computationally more expensive calculations is made to demonstrate the accuracy and predictive capability of the new KS-DFT/3D-RISM-KH method.
Collapse
Affiliation(s)
- David Casanova
- Departament de Química Inorgànica and Centre de Recerca en Química Teòrica, Universitat de Barcelona, Martí i Franquès 1-11, 08028 Barcelona, Spain, National Institute for Nanotechnology, National Research Council of Canada, 421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada, Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta, T6G 2G8, Canada, and Department of Chemistry, University of Calgary, Calgary, Alberta, T2N 1N4, Canada
| | - Sergey Gusarov
- Departament de Química Inorgànica and Centre de Recerca en Química Teòrica, Universitat de Barcelona, Martí i Franquès 1-11, 08028 Barcelona, Spain, National Institute for Nanotechnology, National Research Council of Canada, 421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada, Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta, T6G 2G8, Canada, and Department of Chemistry, University of Calgary, Calgary, Alberta, T2N 1N4, Canada
| | - Andriy Kovalenko
- Departament de Química Inorgànica and Centre de Recerca en Química Teòrica, Universitat de Barcelona, Martí i Franquès 1-11, 08028 Barcelona, Spain, National Institute for Nanotechnology, National Research Council of Canada, 421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada, Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta, T6G 2G8, Canada, and Department of Chemistry, University of Calgary, Calgary, Alberta, T2N 1N4, Canada
| | - Tom Ziegler
- Departament de Química Inorgànica and Centre de Recerca en Química Teòrica, Universitat de Barcelona, Martí i Franquès 1-11, 08028 Barcelona, Spain, National Institute for Nanotechnology, National Research Council of Canada, 421 Saskatchewan Drive, Edmonton, Alberta, T6G 2M9, Canada, Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta, T6G 2G8, Canada, and Department of Chemistry, University of Calgary, Calgary, Alberta, T2N 1N4, Canada
| |
Collapse
|
6
|
Lyubimova O, Stoyanov SR, Gusarov S, Kovalenko A. Electric Interfacial Layer of Modified Cellulose Nanocrystals in Aqueous Electrolyte Solution: Predictions by the Molecular Theory of Solvation. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2015; 31:7106-7116. [PMID: 26053228 DOI: 10.1021/acs.langmuir.5b00680] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
The X-ray crystal structure-based models of Iα cellulose nanocrystals (CNC), both pristine and containing surface sulfate groups with negative charge 0-0.34 e/nm(2) produced by sulfuric acid hydrolysis of softwood pulp, feature a highly polarized "crystal-like" charge distribution. We perform sampling using molecular dynamics (MD) of the structural relaxation of neutral pristine and negatively charged sulfated CNC of various lengths in explicit water solvent and then employ the statistical mechanical 3D-RISM-KH molecular theory of solvation to evaluate the solvation structure and thermodynamics of the relaxed CNC in ambient aqueous NaCl solution at a concentration of 0.0-0.25 mol/kg. The MD sampling induces a right-hand twist in CNC and rearranges its initially ordered structure with a macrodipole of high-density charges at the opposite faces into small local spots of alternating charge at each face. This surface charge rearrangement observed for both neutral and charged CNC significantly affects the distribution of ions around CNC in aqueous electrolyte solution. The solvation free energy (SFE) of charged sulfated CNC has a minimum at a particular electrolyte concentration depending on the surface charge density, whereas the SFE of neutral CNC increases linearly with NaCl concentration. The SFE contribution from Na(+) counterions exhibits behavior similar to the NaCl concentration dependence of the whole SFE. An analysis of the 3D maps of Na(+) density distributions shows that these model CNC particles exhibit the behavior of charged nanocolloids in aqueous electrolyte solution: an increase in electrolyte concentration shrinks the electric interfacial layer and weakens the effective repulsion between charged CNC particles. The 3D-RISM-KH method readily treats solvent and electrolyte of a given nature and concentration to predict effective interactions between CNC particles in electrolyte solution. We provide CNC structural models and a modeling procedure for studies of effective interactions and the formation of ordered phases of CNC suspensions in electrolyte solution.
Collapse
Affiliation(s)
- Olga Lyubimova
- †National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
- ‡Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta T6G 2G8, Canada
| | - Stanislav R Stoyanov
- †National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
- ‡Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta T6G 2G8, Canada
- §Department of Chemical and Materials Engineering, University of Alberta, Edmonton, Alberta T6G 2V4, Canada
| | - Sergey Gusarov
- †National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
| | - Andriy Kovalenko
- †National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
- ‡Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta T6G 2G8, Canada
| |
Collapse
|
7
|
Munaò G, Urbic T. Structure and thermodynamics of core-softened models for alcohols. J Chem Phys 2015; 142:214508. [DOI: 10.1063/1.4922164] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Gianmarco Munaò
- Dipartimento di Fisica e di Scienze della Terra, Università degli Studi di Messina, Viale F. Stagno d’Alcontres 31, 98166 Messina, Italy
| | - Tomaz Urbic
- Department of Chemistry and Chemical Technology, Chair of Physical Chemistry, University of Ljubljana, Večna pot 113, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
8
|
Huš M, Munaò G, Urbic T. Properties of a soft-core model of methanol: an integral equation theory and computer simulation study. J Chem Phys 2014; 141:164505. [PMID: 25362323 DOI: 10.1063/1.4899316] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Thermodynamic and structural properties of a coarse-grained model of methanol are examined by Monte Carlo simulations and reference interaction site model (RISM) integral equation theory. Methanol particles are described as dimers formed from an apolar Lennard-Jones sphere, mimicking the methyl group, and a sphere with a core-softened potential as the hydroxyl group. Different closure approximations of the RISM theory are compared and discussed. The liquid structure of methanol is investigated by calculating site-site radial distribution functions and static structure factors for a wide range of temperatures and densities. Results obtained show a good agreement between RISM and Monte Carlo simulations. The phase behavior of methanol is investigated by employing different thermodynamic routes for the calculation of the RISM free energy, drawing gas-liquid coexistence curves that match the simulation data. Preliminary indications for a putative second critical point between two different liquid phases of methanol are also discussed.
Collapse
Affiliation(s)
- Matej Huš
- Department of Chemistry and Chemical Technology, University of Ljubljana, Chair of Physical Chemistry, Večna pot 113, SI-1000 Ljubljana, Slovenia
| | - Gianmarco Munaò
- Dipartimento di Fisica e di Scienze della Terra, Università degli Studi di Messina, Viale F. Stagno d'Alcontres 31, 98166 Messina, Italy
| | - Tomaz Urbic
- Department of Chemistry and Chemical Technology, University of Ljubljana, Chair of Physical Chemistry, Večna pot 113, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
9
|
Kobryn AE, Nikolić D, Lyubimova O, Gusarov S, Kovalenko A. Dissipative Particle Dynamics with an Effective Pair Potential from Integral Equation Theory of Molecular Liquids. J Phys Chem B 2014; 118:12034-49. [DOI: 10.1021/jp503981p] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Alexander E. Kobryn
- National
Institute for Nanotechnology, National Research Council of Canada, 11421
Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
| | - Dragan Nikolić
- National
Institute for Nanotechnology, National Research Council of Canada, 11421
Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
- Department
of Mechanical Engineering, University of Alberta, Edmonton, Alberta T6G 2G8, Canada
| | - Olga Lyubimova
- National
Institute for Nanotechnology, National Research Council of Canada, 11421
Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
- Department
of Mechanical Engineering, University of Alberta, Edmonton, Alberta T6G 2G8, Canada
| | - Sergey Gusarov
- National
Institute for Nanotechnology, National Research Council of Canada, 11421
Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
| | - Andriy Kovalenko
- National
Institute for Nanotechnology, National Research Council of Canada, 11421
Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
- Department
of Mechanical Engineering, University of Alberta, Edmonton, Alberta T6G 2G8, Canada
| |
Collapse
|
10
|
Kovalenko A. Predictive Multiscale Modeling of Nanocellulose Based Materials and Systems. ACTA ACUST UNITED AC 2014. [DOI: 10.1088/1757-899x/64/1/012040] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
11
|
Omelyan I, Kovalenko A. Multiple time step molecular dynamics in the optimized isokinetic ensemble steered with the molecular theory of solvation: accelerating with advanced extrapolation of effective solvation forces. J Chem Phys 2014; 139:244106. [PMID: 24387356 DOI: 10.1063/1.4848716] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
We develop efficient handling of solvation forces in the multiscale method of multiple time step molecular dynamics (MTS-MD) of a biomolecule steered by the solvation free energy (effective solvation forces) obtained from the 3D-RISM-KH molecular theory of solvation (three-dimensional reference interaction site model complemented with the Kovalenko-Hirata closure approximation). To reduce the computational expenses, we calculate the effective solvation forces acting on the biomolecule by using advanced solvation force extrapolation (ASFE) at inner time steps while converging the 3D-RISM-KH integral equations only at large outer time steps. The idea of ASFE consists in developing a discrete non-Eckart rotational transformation of atomic coordinates that minimizes the distances between the atomic positions of the biomolecule at different time moments. The effective solvation forces for the biomolecule in a current conformation at an inner time step are then extrapolated in the transformed subspace of those at outer time steps by using a modified least square fit approach applied to a relatively small number of the best force-coordinate pairs. The latter are selected from an extended set collecting the effective solvation forces obtained from 3D-RISM-KH at outer time steps over a broad time interval. The MTS-MD integration with effective solvation forces obtained by converging 3D-RISM-KH at outer time steps and applying ASFE at inner time steps is stabilized by employing the optimized isokinetic Nosé-Hoover chain (OIN) ensemble. Compared to the previous extrapolation schemes used in combination with the Langevin thermostat, the ASFE approach substantially improves the accuracy of evaluation of effective solvation forces and in combination with the OIN thermostat enables a dramatic increase of outer time steps. We demonstrate on a fully flexible model of alanine dipeptide in aqueous solution that the MTS-MD/OIN/ASFE/3D-RISM-KH multiscale method of molecular dynamics steered by effective solvation forces allows huge outer time steps up to tens of picoseconds without affecting the equilibrium and conformational properties, and thus provides a 100- to 500-fold effective speedup in comparison to conventional MD with explicit solvent. With the statistical-mechanical 3D-RISM-KH account for effective solvation forces, the method provides efficient sampling of biomolecular processes with slow and/or rare solvation events such as conformational transitions of hydrated alanine dipeptide with the mean life times ranging from 30 ps up to 10 ns for "flip-flop" conformations, and is particularly beneficial for biomolecular systems with exchange and localization of solvent and ions, ligand binding, and molecular recognition.
Collapse
Affiliation(s)
- Igor Omelyan
- National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
| | - Andriy Kovalenko
- National Institute for Nanotechnology, 11421 Saskatchewan Drive, Edmonton, Alberta T6G 2M9, Canada
| |
Collapse
|
12
|
Multiscale modeling of solvation in chemical and biological nanosystems and in nanoporous materials. PURE APPL CHEM 2013. [DOI: 10.1351/pac-con-12-06-03] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Statistical–mechanical, 3D-RISM-KH molecular theory of solvation (3D reference interaction site model with the Kovalenko–Hirata closure) is promising as an essential part of multiscale methodology for chemical and biomolecular nanosystems in solution. 3D-RISM-KH explains the molecular mechanisms of self-assembly and conformational stability of synthetic organic rosette nanotubes (RNTs), aggregation of prion proteins and β-sheet amyloid oligomers, protein-ligand binding, and function-related solvation properties of complexes as large as the Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC) and GroEL/ES chaperone. Molecular mechanics/Poisson–Boltzmann (generalized Born) surface area [MM/PB(GB)SA] post-processing of molecular dynamics (MD) trajectories involving SA empirical nonpolar terms is replaced with MM/3D-RISM-KH statistical–mechanical evaluation of the solvation thermodynamics. 3D-RISM-KH has been coupled with multiple time-step (MTS) MD of the solute biomolecule driven by effective solvation forces, which are obtained analytically by converging the 3D-RISM-KH integral equations at outer time-steps and are calculated in between by using solvation force coordinate extrapolation (SFCE) in the subspace of previous solutions to 3D-RISM-KH. The procedure is stabilized by the optimized isokinetic Nosé–Hoover (OIN) chain thermostatting, which enables gigantic outer time-steps up to picoseconds to accurately calculate equilibrium properties. The multiscale OIN/SFCE/3D-RISM-KH algorithm is implemented in the Amber package and illustrated on a fully flexible model of alanine dipeptide in aqueous solution, exhibiting the computational rate of solvent sampling 20 times faster than standard MD with explicit solvent. Further substantial acceleration can be achieved with 3D-RISM-KH efficiently sampling essential events with rare statistics such as exchange and localization of solvent, ions, and ligands at binding sites and pockets of the biomolecule. 3D-RISM-KH was coupled with ab initio complete active space self-consistent field (CASSCF) and orbital-free embedding (OFE) Kohn–Sham (KS) density functional theory (DFT) quantum chemistry methods in an SCF description of electronic structure, optimized geometry, and chemical reactions in solution. The (OFE)KS-DFT/3D-RISM-KH multi-scale method is implemented in the Amsterdam Density Functional (ADF) package and extensively validated against experiment for solvation thermochemistry, photochemistry, conformational equilibria, and activation barriers of various nanosystems in solvents and ionic liquids (ILs). Finally, the replica RISM-KH-VM molecular theory for the solvation structure, thermodynamics, and electrochemistry of electrolyte solutions sorbed in nanoporous materials reveals the molecular mechanisms of sorption and supercapacitance in nanoporous carbon electrodes, which is drastically different from a planar electrical double layer.
Collapse
|
13
|
Munaò G, Costa D, Giacometti A, Caccamo C, Sciortino F. Structure and phase behavior of colloidal dumbbells with tunable attractive interactions. Phys Chem Chem Phys 2013; 15:20590-9. [DOI: 10.1039/c3cp52425f] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- G Munaò
- Dipartimento di Fisica, Università di Roma "La Sapienza", Piazzale Aldo Moro 2, 00185 Roma, Italy.
| | | | | | | | | |
Collapse
|
14
|
da Costa LM, Hayaki S, Stoyanov SR, Gusarov S, Tan X, Gray MR, Stryker JM, Tykwinski R, Carneiro JWDM, Sato H, Seidl PR, Kovalenko A. 3D-RISM-KH molecular theory of solvation and density functional theory investigation of the role of water in the aggregation of model asphaltenes. Phys Chem Chem Phys 2012; 14:3922-34. [DOI: 10.1039/c2cp23131j] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
15
|
Munaó G, Costa D, Sciortino F, Caccamo C. Simulation and theory of a model for tetrahedral colloidal particles. J Chem Phys 2011; 134:194502. [DOI: 10.1063/1.3582904] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
|
16
|
Stumpe MC, Blinov N, Wishart D, Kovalenko A, Pande VS. Calculation of local water densities in biological systems: a comparison of molecular dynamics simulations and the 3D-RISM-KH molecular theory of solvation. J Phys Chem B 2010; 115:319-28. [PMID: 21174421 DOI: 10.1021/jp102587q] [Citation(s) in RCA: 66] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Water plays a unique role in all living organisms. Not only is it nature's ubiquitous solvent, but it also actively takes part in many cellular processes. In particular, the structure and properties of interfacial water near biomolecules such as proteins are often related to the function of the respective molecule. It can therefore be highly instructive to study the local water density around solutes in cellular systems, particularly when solvent-mediated forces such as the hydrophobic effect are relevant. Computational methods such as molecular dynamics (MD) simulations seem well suited to study these systems at the atomic level. However, due to sampling requirements, it is not clear that MD simulations are, indeed, the method of choice to obtain converged densities at a given level of precision. We here compare the calculation of local water densities with two different methods: MD simulations and the three-dimensional reference interaction site model with the Kovalenko-Hirata closure (3D-RISM-KH). In particular, we investigate the convergence of the local water density to assess the required simulation times for different levels of resolution. Moreover, we provide a quantitative comparison of the densities calculated with MD and with 3D-RISM-KH and investigate the effect of the choice of the water model for both methods. Our results show that 3D-RISM-KH yields density distributions that are very similar to those from MD up to a 0.5 Å resolution, but for significantly reduced computational cost. The combined use of MD and 3D-RISM-KH emerges as an auspicious perspective for efficient solvent sampling in dynamical systems.
Collapse
Affiliation(s)
- Martin C Stumpe
- Department of Bioengineering, Stanford University, 318 Campus Drive West, Stanford, California 94305, USA
| | | | | | | | | |
Collapse
|
17
|
Kaminski JW, Gusarov S, Wesolowski TA, Kovalenko A. Modeling Solvatochromic Shifts Using the Orbital-Free Embedding Potential at Statistically Mechanically Averaged Solvent Density. J Phys Chem A 2010; 114:6082-96. [DOI: 10.1021/jp100158h] [Citation(s) in RCA: 77] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Jakub W. Kaminski
- Département de Chimie Physique, Université de Genève, 30 quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland, National Research Council of Canada. National Institute for Nanotechnology, 421 Saskatchewan Drive T6G 2M9 Edmonton, Canada, and Department of Mechanical Engineering, University of Alberta, T6G 2G8 Edmonton, Canada
| | - Sergey Gusarov
- Département de Chimie Physique, Université de Genève, 30 quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland, National Research Council of Canada. National Institute for Nanotechnology, 421 Saskatchewan Drive T6G 2M9 Edmonton, Canada, and Department of Mechanical Engineering, University of Alberta, T6G 2G8 Edmonton, Canada
| | - Tomasz A. Wesolowski
- Département de Chimie Physique, Université de Genève, 30 quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland, National Research Council of Canada. National Institute for Nanotechnology, 421 Saskatchewan Drive T6G 2M9 Edmonton, Canada, and Department of Mechanical Engineering, University of Alberta, T6G 2G8 Edmonton, Canada
| | - Andriy Kovalenko
- Département de Chimie Physique, Université de Genève, 30 quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland, National Research Council of Canada. National Institute for Nanotechnology, 421 Saskatchewan Drive T6G 2M9 Edmonton, Canada, and Department of Mechanical Engineering, University of Alberta, T6G 2G8 Edmonton, Canada
| |
Collapse
|
18
|
Munaò G, Costa D, Saija F, Caccamo C. Simulation and reference interaction site model theory of methanol and carbon tetrachloride mixtures. J Chem Phys 2010; 132:084506. [DOI: 10.1063/1.3314296] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
19
|
Abstract
We propose a site site generalization of the Lovett-Mow-Buff-Wertheim integro-differential equation for the one-particle density distributions to polyatomic fluids. The method provides microscopic description of liquid interfaces of molecular fluids and solutions. It uses the inhomogeneous site-site direct correlation function of molecular fluid consistently constructed by nonlinear interpolation between the homogeneous ones. The site site correlations of the coexisting bulk phases are obtained from the reference interaction site model (RISM) integral equation with our closure approximation. For illustration, we calculated the structure of the planar liquid-vapor as well as liquid-liquid interfaces of n-hexane and methanol at ambient conditions.
Collapse
Affiliation(s)
- Andriy Kovalenko
- National Institute for Nanotechnology, National Research Council of Canada and Department of Mechanical Engineering, University of Alberta, W6-025, ECERF Bldg., 9107-116 Str., Edmonton, AB, T6G 2V4, Canada.
| | | |
Collapse
|
20
|
Stoyanov SR, Gusarov S, Kovalenko A. Multiscale modelling of asphaltene disaggregation. MOLECULAR SIMULATION 2008. [DOI: 10.1080/08927020802411711] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Stanislav R. Stoyanov
- a National Institute for Nanotechnology , National Research Council of Canada , Edmonton, AB, Canada , T6G 2M9
- b Department of Mechanical Engineering , University of Alberta , Edmonton, AB, Canada , T6G 2G8
| | - Sergey Gusarov
- a National Institute for Nanotechnology , National Research Council of Canada , Edmonton, AB, Canada , T6G 2M9
| | - Andriy Kovalenko
- a National Institute for Nanotechnology , National Research Council of Canada , Edmonton, AB, Canada , T6G 2M9
| |
Collapse
|
21
|
Woelki S, Kohler HH, Krienke H, Schmeer G. Improvements of DRISM calculations: symmetry reduction and hybrid algorithms. Phys Chem Chem Phys 2008; 10:898-910. [DOI: 10.1039/b712306j] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
22
|
Costa D, Munaó G, Saija F, Caccarno C. Reference interaction site model and molecular dynamics study of structure and thermodynamics of methanol. J Chem Phys 2007; 127:224501. [DOI: 10.1063/1.2803059] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
|
23
|
Woelki S, Kohler HH, Krienke H. A Singlet-RISM Theory for Solid/Liquid Interfaces Part I: Uncharged Walls. J Phys Chem B 2007; 111:13386-97. [DOI: 10.1021/jp068998t] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Stefan Woelki
- Institute of Analytical Chemistry, Chemo- and Biosensors, University of Regensburg, D-93040 Regensburg, Germany, and Institute of Physical and Theoretical Chemistry, University of Regensburg, D-93040 Regensburg, Germany
| | - Hans-Helmut Kohler
- Institute of Analytical Chemistry, Chemo- and Biosensors, University of Regensburg, D-93040 Regensburg, Germany, and Institute of Physical and Theoretical Chemistry, University of Regensburg, D-93040 Regensburg, Germany
| | - Hartmut Krienke
- Institute of Analytical Chemistry, Chemo- and Biosensors, University of Regensburg, D-93040 Regensburg, Germany, and Institute of Physical and Theoretical Chemistry, University of Regensburg, D-93040 Regensburg, Germany
| |
Collapse
|
24
|
Tanimura A, Kovalenko A, Hirata F. Structure of electrolyte solutions sorbed in carbon nanospaces, studied by the replica RISM theory. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2007; 23:1507-17. [PMID: 17241081 DOI: 10.1021/la061617i] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
The replica RISM theory is used to investigate the structure of electrolyte solutions confined in carbonized polyvinylidene chloride (PVDC) nanoporous material, compared to bulk electrolyte solution. Comparisons are made between the models of electrolyte solution sorbed in the carbonized PVDC material and a single carbon nanosphere in bulk electrolyte solution. Particular attention is paid to the chemical potential balance between the species of the sorbed electrolyte solution and the bulk solution in contact with the nanoporous material. As a result of the strong hydrophobicity of the carbonized PVDC material in the absence of activating chemical groups, the densities of water and ions sorbed in the material are remarkably low compared to those in the ambient bulk solution. The interaction between water molecules and cations becomes strong in nanospaces. It turns out that, in carbon nanopores, a cation adsorbed at the carbon surface is fully surrounded by the hydration shell of water molecules which separates the cation and the surface. Distinctively, an anion is adsorbed in direct contact with the carbon surface, which squeezes a part of its hydration shell out. The tendency increases toward smaller cations, which are characterized as "positive hydration" ions. In the bulk, cations are not hydrated so strongly and behave similarly to anions. The results suggest that the specific capacitance of an electric double-layer supercapacitor with nanoporous electrodes is intimately related to the solvation structure of electrolyte solution sorbed in nanopores, which is affected by the microscopic structure of the nanoporous electrode.
Collapse
Affiliation(s)
- A Tanimura
- Institute for Molecular Science, Okazaki 444-8585, Japan
| | | | | |
Collapse
|
25
|
Vörtler HL, Kettler M. Efficient simulation of chemical potentials and phase equilibria in associating fluids: monomer/dimer insertion versus gradual particle insertion in primitive water models. Mol Phys 2006. [DOI: 10.1080/00268970500404455] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
26
|
Livadaru L, Kovalenko A. Self-consistent molecular theory of polymers in melts and solutions. J Phys Chem B 2005; 109:10631-9. [PMID: 16852290 DOI: 10.1021/jp044654w] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We propose a self-consistent molecular theory of conformational properties of flexible polymers in melts and solutions. The method employs the polymer reference interaction site model for the intermolecular correlations and the Green function technique for the intramolecular correlations. We demonstrate this method on n-alkane molecules in different environments: water, hexane, and in melt, corresponding to poor, good, and theta condition, respectively. The numerical results of the intramolecular correlation function, the radius of gyration, and the characteristic ratio of a polymer chain are indicative of conformational changes from one environment to another and are in agreement with other findings in the literature. Scaling laws for the chain size with respect to the number of monomers are discussed. We show results for the intra- and intermolecular correlation functions and the medium-induced potential. We also extract the Kuhn length and the characteristic ratio for the infinite chain limit for melts. The latter is compared to the experimental results and computer simulation. The conformational free energy per monomer in different solvents is calculated. Our treatment can be generalized readily to other polymer-solvent systems, for example, those containing branched copolymers and polar solvents.
Collapse
Affiliation(s)
- Lucian Livadaru
- National Institute for Nanotechnology, National Research Council of Canada, W6-010, ECERF Building, 9107-116 Street, Edmonton, AB, Canada, T6G 2V4
| | | |
Collapse
|
27
|
Microscopic description of a liquid–vapor interface by an inhomogeneous integral equation theory. Chem Phys Lett 2004. [DOI: 10.1016/j.cplett.2004.08.114] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
28
|
Livadaru L, Kovalenko A. Molecular description of the collapse of hydrophobic polymer chains in water. J Chem Phys 2004; 121:4449-52. [PMID: 15332873 DOI: 10.1063/1.1789916] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We propose a self-consistent molecular theory of conformational properties of flexible polymers in solution. It is applied to the collapse of a hydrophobic polymer chain in water, and can be readily generalized to any polymer-solvent system (e.g., copolymers with high complexity). We stress the potential of this method for a variety of problems, such as protein folding.
Collapse
Affiliation(s)
- Lucian Livadaru
- National Institute for Nanotechnology, National Research Council of Canada, 9107-116 St., W6-010, ECERF Bldg., Edmonton, AB, T6G 2V4, Canada
| | | |
Collapse
|
29
|
Tanimura A, Kovalenko A, Hirata F. Molecular theory of an electrochemical double layer in a nanoporous carbon supercapacitor. Chem Phys Lett 2003. [DOI: 10.1016/s0009-2614(03)01336-8] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|