Erratum: Gauging van der Waals interactions in aqueous solutions of 2D MOFs: when water likes organic linkers more than open-metal sites (Physical Chemistry Chemical Physics (2021) 23 (3135-3143) DOI: 10.1039/D0CP05923D)

In the originally published article, a few parameters were missing from the force fields affecting the reported results for the classical molecular dynamics simulation data in Section 3.1. (i.e., Fig. 2 and Table 1). Specifically, a number of dihedral angles were missing from the force field. Also, errors were found in the implementation of the q-TIP4P/F water potential which are fixed The simulations were repeated to correct this error, along with optimisation of the simulation details, which have been outlined here. The corrected figure and table are reported below. The corresponding tables of force field parameters along with data related to the MD results have been updated in the ESI. The input and output files of all simulations have also been provided as part of the Supplementary Materials to the original article. For MD simulations of bulk models of 2D MOFs, larger 2x2x3 periodic cells, composed of 1512 atoms and 72 metal centers, were used along with a larger cutoff of 10 Å for treating long-range electrostatic interactions. This is opposed to the 2 x 2x 2 periodic cells in the original paper with a 6 Å cutoff. The simulations were carried out using our in-house modified software package coined as DL_POLY Quantum v1.0 which is publicly available through our GitHub page.1 The corrected Fig. 2 is provided below. The trends observed and the related discussions in the main text on the importance of hydrogen bond formation for the adsorption of water are not changed. The corrected Table 1 is given below. To determine the orientational relaxation time (treor), we used a bi-exponential function2 as in: (Farmula Presented) This is opposed to the single exponential (i.e., eqn (3)) used in the original article. The related graphs are reported in Fig. S11 of the updated Supplementary Materials. For the new data from the bi-exponential function, the final relaxation times were calculated from the weighted average of the fitting parameters as: (Farmula Presented) The first 100 fs of the simulations which correspond to the fast liberational motion of the water molecules2 were excluded from these fits. The mean square displacement (MSD) plots for calculating diffusion coefficients (D) are updated in Fig. S12 and S13 of the new Supplementary Materials. The reported trends and main conclusions on the freer nature of water in Cu₃(HTTP)₂compared to Cu3(HHTP)2 remain unchanged. The noticeable difference of the new results compared to the original data is the increasing trend of both D_z and D_xy with respect to water concentration which is in line with the increasing trend of D_tot. The Royal Society of Chemistry apologises for these errors and any consequent inconvenience to authors and readers. (Figure Presented) (Table Presented).