1
|
Di Pino S, Perez Sirkin YA, Morzan UN, Sánchez VM, Hassanali A, Scherlis DA. Water Self-Dissociation is Insensitive to Nanoscale Environments. Angew Chem Int Ed Engl 2023; 62:e202306526. [PMID: 37379226 DOI: 10.1002/anie.202306526] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2023] [Revised: 06/26/2023] [Accepted: 06/27/2023] [Indexed: 06/30/2023]
Abstract
Nanoconfinement effects on water dissociation and reactivity remain controversial, despite their importance to understand the aqueous chemistry at interfaces, pores, or aerosols. The pKw in confined environments has been assessed from experiments and simulations in a few specific cases, leading to dissimilar conclusions. Here, with the use of carefully designed ab initio simulations, we demonstrate that the energetics of bulk water dissociation is conserved intact to unexpectedly small length-scales, down to aggregates of only a dozen molecules or pores of widths below 2 nm. The reason is that most of the free-energy involved in water autoionization comes from breaking the O-H covalent bond, which has a comparable barrier in the bulk liquid, in a small droplet of nanometer size, or in a nanopore in the absence of strong interfacial interactions. Thus, dissociation free-energy profiles in nanoscopic aggregates or in 2D slabs of 1 nm width reproduce the behavior corresponding to the bulk liquid, regardless of whether the corresponding nanophase is delimited by a solid or a gas interface. The present work provides a definite and fundamental description of the mechanism and thermodynamics of water dissociation at different scales with broader implications on reactivity and self-ionization at the air-liquid interface.
Collapse
Affiliation(s)
- Solana Di Pino
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Buenos Aires, C1428EHA, Argentina
- Condensed Matter and Statistical Physics, International Centre for Theoretical Physics, I-34151, Trieste, Italy
| | - Yamila A Perez Sirkin
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Buenos Aires, C1428EHA, Argentina
| | - Uriel N Morzan
- Condensed Matter and Statistical Physics, International Centre for Theoretical Physics, I-34151, Trieste, Italy
| | - Verónica M Sánchez
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Buenos Aires, C1428EHA, Argentina
| | - Ali Hassanali
- Condensed Matter and Statistical Physics, International Centre for Theoretical Physics, I-34151, Trieste, Italy
| | - Damian A Scherlis
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Buenos Aires, C1428EHA, Argentina
- Condensed Matter and Statistical Physics, International Centre for Theoretical Physics, I-34151, Trieste, Italy
| |
Collapse
|
2
|
Pederson JP, McDaniel J. DFT-based QM/MM with Particle-Mesh Ewald for Direct, Long-Range Electrostatic Embedding. J Chem Phys 2022; 156:174105. [DOI: 10.1063/5.0087386] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a DFT-based, QM/MM implementation with long-range electrostatic embedding achieved by direct real-space integration of the particle mesh Ewald (PME) computed electrostatic potential. The key transformation is the interpolation of the electrostatic potential from the PME grid to the DFT quadrature grid, from which integrals are easily evaluated utilizing standard DFT machinery. We provide benchmarks of the numerical accuracy with choice of grid size and real-space corrections, and demonstrate that good convergence is achieved while introducing nominal computational overhead. Furthermore, the approach requires only small modification to existing software packages, as is demonstrated with our implementation in the OpenMM and Psi4 software. After presenting convergence benchmarks, we evaluate the importance of long-range electrostatic embedding in three solute/solvent systems modeled with QM/MM. Water and BMIM/BF4 ionic liquid were considered as ``simple' and ``complex' solvents respectively, with water and p-phenylenediamine (PPD) solute molecules treated at QM level of theory. While electrostatic embedding with standard real-space truncation may introduce negligible error for simple systems such as water solute in water solvent, errors become more significant when QM/MM is applied to complex solvents such as ionic liquids. An extreme example is the electrostatic embedding energy for oxidized PPD in BMIM/BF4 for which real-space truncation produces severe error even at 2-3 nm cutoff distances. This latter example illustrates that utilization of QM/MM to compute redox potentials within concentrated electrolytes/ionic media requires carefully chosen long-range electrostatic embedding algorithms, with our presented algorithm providing a general and robust approach.
Collapse
Affiliation(s)
| | - Jesse McDaniel
- Chemistry, Georgia Institute of Technology, United States of America
| |
Collapse
|
3
|
Ougaard Dohn A, Selli D, Fazio G, Ferraro L, Mortensen JJ, Civalleri B, Di Valentin C. Interfacing CRYSTAL/AMBER to Optimize QM/MM Lennard⁻Jones Parameters for Water and to Study Solvation of TiO₂ Nanoparticles. Molecules 2018; 23:molecules23112958. [PMID: 30428551 PMCID: PMC6278561 DOI: 10.3390/molecules23112958] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2018] [Revised: 11/01/2018] [Accepted: 11/06/2018] [Indexed: 11/28/2022] Open
Abstract
Metal oxide nanoparticles (NPs) are regarded as good candidates for many technological applications, where their functional environment is often an aqueous solution. The correct description of metal oxide electronic structure is still a challenge for local and semilocal density functionals, whereas hybrid functional methods provide an improved description, and local atomic function-based codes such as CRYSTAL17 outperform plane wave codes when it comes to hybrid functional calculations. However, the computational cost of hybrids are still prohibitive for systems of real sizes, in a real environment. Therefore, we here present and critically assess the accuracy of our electrostatic embedding quantum mechanical/molecular mechanical (QM/MM) coupling between CRYSTAL17 and AMBER16, and demonstrate some of its capabilities via the case study of TiO2 NPs in water. First, we produced new Lennard–Jones (LJ) parameters that improve the accuracy of water–water interactions in the B3LYP/TIP3P coupling. We found that optimizing LJ parameters based on water tri- to deca-mer clusters provides a less overstructured QM/MM liquid water description than when fitting LJ parameters only based on the water dimer. Then, we applied our QM/MM coupling methodology to describe the interaction of a 1 nm wide multilayer of water surrounding a spherical TiO2 nanoparticle (NP). Optimizing the QM/MM water–water parameters was found to have little to no effect on the local NP properties, which provide insights into the range of influence that can be attributed to the LJ term in the QM/MM coupling. The effect of adding additional water in an MM fashion on the geometry optimized nanoparticle structure is small, but more evident effects are seen in its electronic properties. We also show that there is good transferability of existing QM/MM LJ parameters for organic molecules–water interactions to our QM/MM implementation, even though these parameters were obtained with a different QM code and QM/MM implementation, but with the same functional.
Collapse
Affiliation(s)
- Asmus Ougaard Dohn
- Faculty of Physical Sciences and Science Institute, University of Iceland, 107 Reykjavík, Iceland.
| | - Daniele Selli
- Dipartimento di Scienza dei Materiali, Università di Milano-Bicocca, via Cozzi 55, 20125 Milano, Italy.
| | - Gianluca Fazio
- Dipartimento di Scienza dei Materiali, Università di Milano-Bicocca, via Cozzi 55, 20125 Milano, Italy.
| | - Lorenzo Ferraro
- Dipartimento di Scienza dei Materiali, Università di Milano-Bicocca, via Cozzi 55, 20125 Milano, Italy.
| | - Jens Jørgen Mortensen
- CAMD, Department of Physics, Technical University of Denmark, 2800 Kongens Lyngby, Denmark.
| | - Bartolomeo Civalleri
- Dipartimento di Chimica, Università di Torino and NIS Centre of Excellence, Via P. Giuria 7, I-10129 Torino, Italy.
| | - Cristiana Di Valentin
- Dipartimento di Scienza dei Materiali, Università di Milano-Bicocca, via Cozzi 55, 20125 Milano, Italy.
| |
Collapse
|
4
|
Sirkin YAP, Hassanali A, Scherlis DA. One-Dimensional Confinement Inhibits Water Dissociation in Carbon Nanotubes. J Phys Chem Lett 2018; 9:5029-5033. [PMID: 30113846 DOI: 10.1021/acs.jpclett.8b02183] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
The effect of nanoconfinement on the self-dissociation of water constitutes an open problem whose elucidation poses a serious challenge to experiments and simulations alike. In slit pores of width ≈1 nm, recent first-principles calculations have predicted that the dissociation constant of H2O increases by almost 2 orders of magnitude [ Muñoz-Santiburcio and Marx, Phys. Rev. Lett. 2017 , 119 , 056002 ]. In the present study, quantum mechanics-molecular mechanics simulations are employed to compute the dissociation free-energy profile of water in a (6,6) carbon nanotube. According to our results, the equilibrium constant Kw drops by 3 orders of magnitude with respect to the bulk phase value, at variance with the trend predicted for confinement in two dimensions. The higher barrier to dissociation can be ascribed to the undercoordination of the hydroxide and hydronium ions in the nanotube and underscores that chemical reactivity does not exhibit a monotonic behavior with respect to pore size but may vary substantially with the characteristic length scale and dimensionality of the confining media.
Collapse
Affiliation(s)
- Yamila A Perez Sirkin
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE, Facultad de Ciencias Exactas y Naturales , Universidad de Buenos Aires, Ciudad Universitaria , Buenos Aires C1428EHA , Argentina
| | - Ali Hassanali
- Condensed Matter and Statistical Physics , International Centre for Theoretical Physics , I-34151 Trieste , Italy
| | - Damián A Scherlis
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE, Facultad de Ciencias Exactas y Naturales , Universidad de Buenos Aires, Ciudad Universitaria , Buenos Aires C1428EHA , Argentina
| |
Collapse
|
5
|
Morzan UN, Alonso de Armiño DJ, Foglia NO, Ramírez F, González Lebrero MC, Scherlis DA, Estrin DA. Spectroscopy in Complex Environments from QM–MM Simulations. Chem Rev 2018; 118:4071-4113. [DOI: 10.1021/acs.chemrev.8b00026] [Citation(s) in RCA: 97] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Uriel N. Morzan
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Diego J. Alonso de Armiño
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Nicolás O. Foglia
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Francisco Ramírez
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Mariano C. González Lebrero
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Damián A. Scherlis
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Darío A. Estrin
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| |
Collapse
|