Reverse Mössbauer effect as a possible source of "hot" molecules absorbed in crystalline solids at low temperature.

As an attempt to explain some of the many anomalies and unresolved problems which have been reported about the dynamic behavior of particles and molecules absorbed in crystalline solids, the "reverse Mössbauer effect" (RME) is proposed. RME theory posits that a particle in non-equilibrium state with respect to a crystal (colliding with the crystal or absorbed in it, but set out of thermal equilibrium by some external cause) is scattered by the whole crystal with a momentum proportional to a vector representing a reciprocal lattice point. The scattering is expected to occur with a well-defined probability and the momentum transferable to the particle is expected to follow a predictable distribution. The RME theory, in practice, is an extension of the Bragg-von Laue scattering law to high-energy colliding particles, in general, and can be applied to any particle or molecule colliding with the surface of a crystalline solid or absorbed in it, but not in thermal equilibrium with the crystal lattice. We verified the RME theory by considering a well-defined unresolved problem. In an experimental study about methane adsorbed in the zeolite Na-ZSM-5 [H. Jobic, Chem. Phys. Lett. 170, 217 (1990)] reporting neutron inelastic-scattering spectra (recoiled bands) at 10 K, the translational kinetic energy of methane resulted to be much higher than equilibrium expected value, namely, about 85 K (or 7.3 meV). The author concluded that "the interpretation of this unusual behavior has yet to be found." In the present study, on the basis of the RME, an explanation of this behavior is put forward.

A network of discrete events for the representation and analysis of diffusion dynamics.

We developed a coarse-grained description of the phenomenology of diffusive processes, in terms of a space of discrete events and its representation as a network. Once a proper classification of the discrete events underlying the diffusive process is carried out, their transition matrix is calculated on the basis of molecular dynamics data. This matrix can be represented as a directed, weighted network where nodes represent discrete events, and the weight of edges is given by the probability that one follows the other. The structure of this network reflects dynamical properties of the process of interest in such features as its modularity and the entropy rate of nodes. As an example of the applicability of this conceptual framework, we discuss here the physics of diffusion of small non-polar molecules in a microporous material, in terms of the structure of the corresponding network of events, and explain on this basis the diffusivity trends observed. A quantitative account of these trends is obtained by considering the contribution of the various events to the displacement autocorrelation function.

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 interplay between dynamic heterogeneities and structure of bulk liquid water: A molecular dynamics simulation study.

In order to study the interplay between dynamical heterogeneities and structural properties of bulk liquid water in the temperature range 130-350 K, thus including the supercooled regime, we use the explicit trend of the distribution functions of some molecular properties, namely, the rotational relaxation constants, the atomic mean-square displacements, the relaxation of the cross correlation functions between the linear and squared displacements of H and O atoms of each molecule, the tetrahedral order parameter q and, finally, the number of nearest neighbors (NNs) and of hydrogen bonds (HBs) per molecule. Two different potentials are considered: TIP4P-Ew and a model developed in this laboratory for the study of nanoconfined water. The results are similar for the dynamical properties, but are markedly different for the structural characteristics. In particular, for temperatures higher than that of the dynamic crossover between "fragile" (at higher temperatures) and "strong" (at lower temperatures) liquid behaviors detected around 207 K, the rotational relaxation of supercooled water appears to be remarkably homogeneous. However, the structural parameters (number of NNs and of HBs, as well as q) do not show homogeneous distributions, and these distributions are different for the two water models. Another dynamic crossover between "fragile" (at lower temperatures) and "strong" (at higher temperatures) liquid behaviors, corresponding to the one found experimentally at T(∗) ∼ 315 ± 5 K, was spotted at T(∗) ∼ 283 K and T(∗) ∼ 276 K for the TIP4P-Ew and the model developed in this laboratory, respectively. It was detected from the trend of Arrhenius plots of dynamic quantities and from the onset of a further heterogeneity in the rotational relaxation. To our best knowledge, it is the first time that this dynamical crossover is detected in computer simulations of bulk water. On the basis of the simulation results, the possible mechanisms of the two crossovers at molecular level are discussed.

Thermodynamics of the one-dimensional parallel Kawasaki model: exact solution and mean-field approximations.

The adsorption isotherm for the recently proposed parallel Kawasaki (PK) lattice-gas model [Phys. Rev. E 88, 062144 (2013)] is calculated exactly in one dimension. To do so, a third-order difference equation for the grand-canonical partition function is derived and solved analytically. In the present version of the PK model, the attraction and repulsion effects between two neighboring particles and between a particle and a neighboring empty site are ruled, respectively, by the dimensionless parameters ϕ and θ. We discuss the inflections induced in the isotherms by situations of high repulsion, the role played by finite lattice sizes in the emergence of substeps, and the adequacy of the two most widely used mean-field approximations in lattice gases, namely, the Bragg-Williams and the Bethe-Peierls approximations.

Combining off-lattice Monte Carlo and cellular automata for the simulation of hard-sphere systems.

In the present work we show how the update rule of a diffusive cellular automaton with mutual exclusion can be exploited in off-lattice Monte Carlo simulations of hard spheres to obtain a synchronous Monte Carlo sampling that satisfies the detailed balance principle.

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.

Synchronous equilibrium model for the diffusion of mutually exclusive particles in a heterogeneous lattice of adsorption sites.

Through straight synchronization and proper manipulation of a sequential Monte Carlo glass-forming rule introduced by Fröbose and Jäckle [J. Stat. Phys. 42, 551 (1986)], we constructed a synchronous, non-glass-forming rule for diffusion of mutually exclusive particles in a lattice of adsorption sites. The rule satisfies detailed balance in the presence of both homogeneous and heterogeneous adsorption energies. Our model differs from the usual lattice-gas cellular automata diffusion rules in that the mutual exclusion holds on the lattice sites rather than on the channels which connect neighboring sites, and from the mass-conserving cellular automata rules in the use of a no-partitioning scheme. The first aim of this work is to show that, although some prescriptions in the synchronous rule are introduced just to allow that both detailed balance and mutual exclusion can coexist with synchronicity, the diffusion process produced by the rule is not anomalous so that the rule can be regarded as a diffusion model. We then compare the diffusion isotherms of several test systems with the ones obtained by means of sequential Monte Carlo simulations of Arrhenius jumps of particles on a lattice. Finally, we apply the rule to the case of a (100) fcc model surface and estimate the amount of time correlation in the migration process, and show that the synchronous rule produces higher correlations and slightly lower diffusivity than the sequential Monte Carlo counterpart.

Conciliating synchronicity with spatial discretization, exclusion, interactions, and detailed balance.

The construction of a discrete stochastic system of interacting particles that evolves through a fully synchronous evolution rule while satisfying detailed balance is a highly demanding task. As a consequence, the presence of nontrivial interaction fields can make synchronicity and thermodynamic equilibrium look as two conflicting counterparts. We show that, with the proper prescriptions, the process of migration of particles in a lattice of mutually exclusive nodes can be simulated with a fully synchronous algorithm, which we call parallel Kawasaki dynamics (PKD), that incorporates site exclusion, local interactions, and detailed balance without the need of system partitioning schemes. We show that the underlying pseudo-Hamiltonian (which is derived from the PKD dynamics instead of being assumed a priori as usual in a sequential Monte Carlo scheme) is temperature dependent and causes the resulting equilibrium properties to differ substantially from the conventional hopping model when the system is near critical conditions.

Chemical potential evaluation in NVT lattice-gas simulations.

The discrete nature of the partition function of a lattice-gas system can be exploited to build an efficient strategy for the evaluation of the chemical potential of a periodic lattice-gas with arbitrarily ranged interactions during a simulation in the canonical ensemble, with the need of no additional sampling as it were required instead by the Widom insertion/deletion approach. The present method is based on the main concepts of the small system grand ensemble [for details, see G. Soto-Campos, D. S. Corti, and H. Reiss, J. Chem. Phys. 108, 2563 (1998)], whose key idea is to study the properties of a sublattice (called small system) and of its complementary (the reservoir sublattice) as they were two separated subsystems. The accuracy of the measured chemical potential can be further improved by artificially "restoring" the missing connections among the reservoir sublattice sites located at the boundaries with the small system. We first illustrate the theory and then we compare µVT with NVT simulation results on several test systems.

A parallelizable block cellular automaton for the study of diffusion of binary mixtures containing CO2 in microporous materials.

We applied a method based on a block cellular automaton (BCA) algorithm to the study of diffusion of various binary mixtures adsorbed in a model microporous material, such as zeolite ZK4. Our aim was to test the capability of our model to cope with systems in which more than one species is present, using a set of parameters based on heuristic considerations from the molecular dynamics (MD) results present in the literature. A rigorous methodology for the assignment of suitable adsorption energies and diffusion activation barriers for our BCA has not been developed yet, nonetheless the results were quite interesting at this stage and we obtained a good qualitative agreement with MD data in the literature. The mixtures we investigated contain CO(2), which causes the so-called segregation-effect, a strong suppression of self-diffusivity of co-adsorbed species. This effect gives rise to relevant problems in the application of some well established and robust methods, while our model proved to be able to reproduce both the common features and the segregation anomaly in the trends of diffusion.

Speeding up simulation of diffusion in zeolites by a parallel synchronous kinetic Monte Carlo algorithm.

Understanding the behaviors of molecules in tight confinement is a challenging task. Standard simulation tools like kinetic Monte Carlo have proven to be very effective in the study of adsorption and diffusion phenomena in microporous materials, but they turn out to be very inefficient when simulation time and length scales are extended. In this paper we have explored the possibility of application of a discrete version of the synchronous parallel kinetic Monte Carlo algorithm introduced by Martínez et al. [J. Comput. Phys. 227, 3804 (2008)] to the study of aromatic hydrocarbons diffusion in zeolites. The efficiency of this algorithm is investigated as a function of the number of processors and domain size. We show that with an accurate choice of domains size it is possible to achieve very good efficiencies thus permitting us to effectively extend space and time scales of the simulated system.

From thermodynamic cell models to partitioning cellular automata for diffusion in zeolites. I. Structure of the algorithm.

In the study of adsorption of simple adsorbates in microporous materials like zeolites, thermodynamic models of small grand-canonical cells with very local interactions [e.g., see K. G. Ayappa, J. Chem. Phys. 111, 4736 (1999)] have been proven to be able to produce thermodynamic properties in very good agreement with the results of experiments and atomistic simulations. In this paper we present in details the structure and implementation of a thermodynamic partitioning cellular automaton (PCA) devised as a dynamical version of thermodynamic cell models and proposed as an easy environment to perform coarse-grained simulations of adsorption/diffusion of simple interacting molecules in microporous materials. Local evolution rules and memory effects are introduced to make our PCA able to complete the static picture provided by thermodynamic cell models with the simulation of transport properties.

From thermodynamic cell models to partitioning cellular automata for diffusion in zeolites. II. Static and dynamic properties.

In this second paper we exploit our thermodynamic partitioning cellular automaton (PCA) developed in Paper I [Pazzona et al., J. Chem. Phys. 131, 234703 (2009)] to study interacting molecules adsorbed in microporous materials. We present a mean-field theory of the single cell model at equilibrium followed by a detailed description of the procedure we propose to calculate the chemical potential in the canonical ensemble. Finally we use our approach to simulate transport properties starting from the parameterization devised by Ayappa [J. Chem. Phys. 111, 4736 (1999)] to reproduce the adsorption properties of xenon in zeolite NaA. We report how the correlations included in the PCA evolution rule affect the estimated self-diffusion coefficient.

Effective interactions in multisite cells for adsorption in microporous materials.

Local, discrete models of self-interacting multisite adsorption cells have been shown to be able to provide a coarse-grained representation of equilibrium properties of small molecules adsorbed in nanoporous materials at the mesoscopic scale. In the present work we show how the essential statistical properties of a host cell of structured sites with multiple adsorption energies and particle-particle interactions (that is the partition function, the average energy, and the average number of guests close to the windows connecting the cell to its surroundings) can be reproduced by a less-structured cell with two occupancy-dependent adsorption energy levels.

Introducing a cellular automaton as an empirical model to study static and dynamic properties of molecules adsorbed in zeolites.

A new lattice gas cellular automaton (LGCA) simulation approach to study static and dynamic properties of molecules adsorbed in zeolites is proposed. The motivation for the present work arises from the ongoing effort to develop efficient numerical tools where conventional approaches like molecular dynamics and Monte Carlo have been revealed as inefficient for a real extension of length and time scales in such inhomogeneous systems. Our LGCA is constituted by a constant number of interacting identical particles, distributed among a fixed number of identical cells arranged in a three-dimensional cubic network and performing a synchronous random walk at constant temperature. The main input for our model comes from data such as (i) local density dependent mean-field potentials and transition probabilities obtained from atomistic simulations that will be used as the starting point to derive adsorption and diffusion properties and (ii) thermodynamic and kinetic data obtained from experiments and/or other simulation methods. Our numerically less demanding LGCA has been tested over three different systems. The obtained results are in excellent agreement with the experimental and theoretical reported data.

Dynamical properties of confined water nanoclusters: Simulation study of hydrated zeolite NaA: structural and vibrational properties.

Water nanoclusters confined to zeolitic cavities have been extensively investigated by various experimental techniques. We report a series of molecular dynamics simulations at different temperatures and for water nanoclusters of different sizes in order to attempt an atomistic interpretation of the properties of these systems. The cavities of zeolite NaA are spherical in shape and about 1 nm in diameter and can host nanoclusters of water containing nearly up to 24 water molecules. A modified interaction potential, yielding a better reproduction of experimental hydration energy and water diffusivity across a number of different zeolites, is proposed. Molecular dynamics simulations reproduce the known experimental structural features obtained by X-ray diffraction. Variations of simulated vibrational IR and IINS spectra with temperature and size of nanoclusters are in good agreement with experiment. The simulated water nanoclusters in zeolite NaA are found to be too small to crystallize and, at low temperature, behave as amorphous ice, in agreement with recent experimental results for similar water nanoclusters in reverse micelles.

Diffusion in tight confinement: a lattice-gas cellular automaton approach. I. Structural equilibrium properties.

The thermodynamic and transport properties of diffusing species in microporous materials are strongly influenced by their interactions with the confining framework, which provide the energy landscape for the transport process. The simple topology and the cellular nature of the alpha cages of a ZK4 zeolite suggest that it is appropriate to apply to the study of the problem of diffusion in tight confinement a time-space discrete model such as a lattice-gas cellular automaton (LGCA). In this paper we investigate the properties of an equilibrium LGCA constituted by a constant number of noninteracting identical particles, distributed among a fixed number of identical cells arranged in a three-dimensional cubic network and performing a synchronous random walk at constant temperature. Each cell of this network is characterized by a finite number of two types of adsorption sites: the exit sites available to particle transfer and the inner sites not available to such transfers. We represent the particle-framework interactions by assuming a differentiation in binding energy of the two types of sites. This leads to a strong dependence of equilibrium and transport properties on loading and temperature. The evolution rule of our LGCA model is constituted by two operations (randomization, in which the number of particles which will be able to try a jump to neighboring cells is determined, and propagation, in which the allowed jumps are performed), each one applied synchronously to all of the cells. The authors study the equilibrium distribution of states and the adsorption isotherm of the model under various conditions of loading and temperature. In connection with the differentiation in energy between exit and inner sites, the adsorption isotherm is described by a conventional Langmuir isotherm at high temperature and by a dual-site Langmuir isotherm at low temperature, while a first order diffuse phase transition takes place at very low temperature.

Diffusion in tight confinement: a lattice-gas cellular automaton approach. II. Transport properties.

In this second paper the authors study the transport properties of the lattice-gas cellular automaton presented in Paper I [J. Chem. Phys. 126, 194709 (2007)] to model adsorption and dynamics of particles in a lattice of confining cells. Their work shows how a surprisingly simple parallel rule applied to a static network of cells joined by links set in space and time can generate a wide range of dynamical behaviors. In their model the cells are the elementary constituent objects of the network. They are a portion of space structured in sites which are energetically different. Each cell can accommodate a given maximum number of particles, and each pair of neighboring cells can exchange at most one particle at a time. The predictions of the model are in qualitative agreement with both experimental observations and molecular dynamics simulation results.

Molecular dynamics simulation study of superhydrated perdeuterated natrolite using a new interaction potential model.

To test a new interaction potential, molecular dynamics simulations of zeolite natrolite were performed for the structures under ambient conditions hydrated by perdeuterated water and at high pressure (1.87 GPa) in the superhydrated phase, which were recently studied by neutron diffraction. The experimental structures were reproduced with reasonable accuracy, and the hydrogen bond features are discussed. As in ordinary natrolite, a flip motion of water molecules around the HOH bisector is found, which, together with translational oscillations, gives rise to transient hydrogen bonds between water molecules, which do not appear from experimental equilibrium coordinates. The dynamics of water molecules can explain some problems encountered in refining the experimental structure. Vibrational spectra of natrolite containing perdeuterated water, which are not yet measured, were simulated, and their qualitative trend is discussed.