Andrea Gabrieli

Partial Charges in Periodic Systems: Improving Electrostatic Potential (ESP) Fitting via Total Dipole Fluctuations and Multiframe Approaches

Two major improvements to the state-of-the-art Repeating Electrostatic Potential Extracted Atomic (REPEAT) method, for generating accurate partial charges for molecular simulations of periodic structures, are here developed. The first, D-REPEAT, consists in the simultaneous fit of the electrostatic potential (ESP), together with the total dipole fluctuations (TDF) of the framework. The second, M-REPEAT, allows the fit of multiple ESP configurations at once. When both techniques are fused into one, DM-REPEAT method, the resulting charges become remarkably stable over a large set of fitting regions, giving a robust and physically sound solution to the buried atoms problem. The method capabilities are extensively studied in ZIF-8 framework, and subsequently applied to IRMOF-1 and ITQ-29 crystal structures. To our knowledge, this is the first time that this approach is proposed in the context of periodic systems.

The central cell model: a mesoscopic hopping model for the study of the displacement autocorrelation function.

On the mesoscale, the molecular motion in a microporous material can be represented as a sequence of hops between different pore locations and from one pore to the other. On the same scale, the memory effects in the motion of a tagged particle are embedded in the displacement autocorrelation function (DACF), the discrete counterpart of the velocity autocorrelation function (VACF). In this paper, a mesoscopic hopping model, based on a lattice-gas automata dynamics, is presented for the coarse-grained modeling of the DACF in a microporous material under conditions of thermodynamic equilibrium. In our model, that we will refer to as central cell model, the motion of one tagged particle is mimicked through probabilistic hops from one location to the other in a small lattice of cells where all the other particles are indistinguishable; the cells closest to the one containing the tagged particle are simulated explicitly in the canonical ensemble, whereas the border cells are treated as mean-field cells in the grand-canonical ensemble. In the present paper, numerical simulation of the central cell model are shown to provide the same results as a traditional lattice-gas simulation. Along with this a mean-field theory of self-diffusion which incorporates time correlations is discussed.

A coarse-grained method based on the analysis of short molecular dynamics trajectories for the simulation of non-Markovian dynamics of molecules adsorbed in microporous materials

We developed a coarse-grained model suitable for the study of adsorbed molecules in microporous materials. A partition of the space available to the motion of adsorbed molecules was carried out, which allows to formulate the dynamics in terms of jumps between discrete regions. The probabilities of observing given pairs of successive jumps were calculated from Molecular Dynamics (MD) simulations, performed on small systems, and used to drive the motion of molecules in a lattice-gas model. Dynamics is thus reformulated in terms of event-space dynamics and this allows to treat the system despite its inherent non markovity. Despite the assumptions enforced in the algorithm, results show that it can be applied to various spherical molecules adsorbed in the all-silica zeolite ITQ-29, establishing a suitable direct bridge between MD simulation results and coarse-grained models.

Local free energies for the coarse-graining of adsorption phenomena: The interacting pair approximation

We investigate the coarse-graining of host-guest systems under the perspective of the local distribution of pore occupancies, along with the physical meaning and actual computability of the coarse-interaction terms. We show that the widely accepted approach, in which the contributions to the free energy given by the molecules located in two neighboring pores are estimated through Monte Carlo simulations where the two pores are kept separated from the rest of the system, leads to inaccurate results at high sorbate densities. In the coarse-graining strategy that we propose, which is based on the Bethe-Peierls approximation, density-independent interaction terms are instead computed according to local effective potentials that take into account the correlations between the pore pair and its surroundings by means of mean-field correction terms without the need for simulating the pore pair separately. Use of the interaction parameters obtained this way allows the coarse-grained system to reproduce more closely the equilibrium properties of the original one. Results are shown for lattice-gases where the local free energy can be computed exactly and for a system of Lennard-Jones particles under the effect of a static confining field.