Automated Multiscale Approach To Predict Self-Diffusion from a Potential Energy Field

Senja Dominique Barthel, Amber Kashan Mace, Berend Smit
2019
Article

Résumé

For large-scale screening studies there is a need to estimate the diffusion of gas molecules in nanoporous materials more efficiently than (brute force) molecular dynamics. In particular for systems with low diffusion coefficients molecular dynamics can be prohibitively expensive. An alternative is to compute the hopping rates between adsorption sites using transition state theory. For large-scale screening this requires the automatic detection of the transition states between the adsorption sites along the different diffusion paths. Here an algorithm is presented that analyzes energy grids for the moving particles. It detects the energies at which diffusion paths are formed, together with their directions. This allows for easy identification of nondiffusive systems. For diffusive systems, it partitions the grid coordinates assigned to energy basins and transitions states, permitting a transition state theory based analysis of the diffusion. We test our method on CH4 diffusion in zeolites, using a standard kinetic Monte Carlo simulation based on the output of our grid analysis. We find that it is accurate, fast, and rigorous without limitations to the geometries of the diffusion tunnels or transition states.

Official source

https://infoscience.epfl.ch/record/287152?ln=fr

À propos de ce résultat

Cette page est générée automatiquement et peut contenir des informations qui ne sont pas correctes, complètes, à jour ou pertinentes par rapport à votre recherche. Il en va de même pour toutes les autres pages de ce site. Veillez à vérifier les informations auprès des sources officielles de l'EPFL.

Concepts associés

Chargement

Publications associées

Chargement

Granular materials are large sets of macroscopic particles that interact solely via contact forces. The static behavior depends on the contact network and on the surface friction forces between grains; when they are set in motion (typically by vibrations) their dynamics is dominated by inelastic collisions. For these reasons granular media show an extremely rich phenomenology, ranging from fluid-like properties (if strongly vibrated), to "jamming", glassy, behavior (if weakly vibrated), to aging and hysteretical phenomena observed when they become trapped in frozen, amorphous states. The objective of this work is to study these states and transitions, and to characterize the analogies found between the dynamic behavior of vibrated granular media and the glass transition observed in thermal glass-formers. These analogies justify the interest in granular materials, because granular media can be seen as simplified model systems useful in the study of out of equilibrium thermodynamics, and, in general, to the larger framework known as "complexity". The granular materials considered here are composed of spheric, polished glass spheres. Since the surface state plays an important role in the grain-grain interaction, some measurements were also performed with acid etched beads, having different surface roughness. The samples are vertically vibrated to achieve vibrofluidization. Different kinds of vibration are used, to highlight different properties of the system. We first consider the transition between the fluid and the subcooled glassy phase, using different experimental techniques. The most important one is a torsion oscillator, that interacts with the granular media via immersed probes. The torsion oscillator can be used in forced mode. A torque is applied on the probe, and we measure the mechanical response function (complex susceptibility). In general, a relaxation is found and it is interpreted as the signature of the irreversible energy loss (damping) in granular collisions. This relaxation has an intrinsic time scale, and systematic analysis of it shows that a clear parallel can be traced to the behavior of "strong" glasses. In particular, it is found that (i) the relaxation time is a function of a unified control parameter, proportional to the square root of the average vibration, and phenomenologically equivalent to an effective temperature; (ii) the functional form with which the relaxation times approach the final "frozen" state has an Arrhenius, or Vögel-Fulcher-Tamman (VFT) behavior. The same torsion oscillator is employed in free mode. In this case, no external torque is applied, and the probe moves adapting its position under the effect of the continuous rearrangements in the sample. The system is studied by computing the power spectral density of the (angular position) time series. The resulting spectra represent a "configurational noise" as the system randomly hops from one configuration to the following. This allows to define, using a completely different approach, the same intrinsic time scale observed in forced mode measurements. The comparison of the two techniques allows to obtain a more complete and detailed picture of the dynamics in the jamming region. From this comparison, it was inferred that the system is also influenced by an effective vibration frequency, and that the relaxation time has indeed a non-Arrhenius behavior as a function of a control parameter defined as as = √ Γ/ωs. A model was developed combining rheological observations to a statistic approach describing extremal phenomena. This model justifies the appearance of both the control parameter and the VFT evolution of the relaxation. Furthermore, the model is predictive and exposes the effect of a few other rheologic properties of granular system. The effect of surface roughness are considered, showing that the static and dynamic surface friction coefficients are well described by the model. A second relevant part of this work is devoted to an explicit verification that macroscopic probes act as Brownian objects. This fact is often used to interpret experimental data (also in the present work) and to propose theoretical model. However, no explicit evidence has ever been discussed. This is hard to do, using a constrained system such as the torsion oscillator, because the restoring coefficient influences the dynamics of diffusion. To overcome the problem we built a different apparatus, called "Brownian motor", where the probes are mounted on ball bearings, so that they are free to turn without constraint. The properties of the time series of the position of the free turning probe and of the torsionally constrained oscillator can finally be analyzed and compared with simple simulations. The data show an overall diffusion-like behavior, that is influenced by the presence of constraints. Using fractal analysis we estimate the diffusion, or Hurst exponent. This allows to verify that a "macroscopic" object (the probe) immersed in the "microscopic" granular medium indeed behaves as a Brownian object, and that its dynamics can be studied in detail, showing that it undergoes anomalous diffusion. This work is concluded with a discussion on a few possible developments. The most promising idea is a novel approach to the study of the geometrical properties of the contact network of granular assemblies, that is responsible for many of the properties of the granular sample. By using Magnetic Resonance Imaging, the static 3-D structure of granular media can be reconstructed with unprecedented accuracy, resolution and ease of reproducibility. From the spatial information we can extract all the properties of static granular media: the compaction factor, the grain-grain correlation function, the free volume and other observables. Systematic studies could allow experimental confirmations of the many theoretical models that have been proposed in the last years and that still lack a thorough comparison with experiments. This idea does not conclude the perspectives of this work, that are vast and intriguing. A few promising subjects are reviewed more into detail in the corresponding Perspective section. To name a few we cite: measurements of induced aging in non-vibrated samples, the Brownian motor, stick and slip phenomena and their comparison with earthquakes.

EPFL2012

Chargement

What Are the Best Materials To Separate a Xenon/Krypton Mixture?

Rocio Mercado, Sondre Kvalvag Schnell, Cory Matthew Simon, Berend Smit

Accelerating progress in the discovery and deployment of advanced nanoporous materials relies on chemical insight and structure-property relationships for rational design. Because of the complexity of this problem, trial-and-error is heavily involved in the laboratory today. A cost-effective route to aid experimental materials discovery is to construct structure models of nanoporous materials in silico and use molecular simulations to rapidly test them and elucidate data-driven guidelines for rational design. For example, highly tunable nanoporous materials have shown promise as adsorbents for separating an industrially relevant gaseous mixture of xenon and krypton. In this work, we characterize, screen, and analyze the Nanoporous Materials Genorne, a database of about 670 000 porous material structures, for candidate adsorbents for xenon/krypton separations. For over half a million structures, the computational resources required for a brute-force screening using grand-canonical Monte Carlo simulations of Xe/Kr adsorption are prohibitive. To overcome the computational cost, we used a hybrid approach combining machine learning algorithms (random forests) with molecular simulations. We compared the results from our large-scale screening with simple pore models to rationalize the strong link between pore size and selectivity. With this insight, we then analyzed the anatomy of the binding sites of the most selective materials. These binding sites can be constructed from tubes, pockets, rings, or cages and are often composed of nondiscrete chemical fragments. The complexity of these binding sites emphasizes the importance of high-throughput computational screenings to identify optimal materials for a given application. Interestingly, our screening study predicts that the two most selective materials in the database are an aluminophosphate zeolite analogue and a calcium based coordination network, both of which have already been synthesized but not yet tested for Xe/Kr separations.

American Chemical Society2015

Chargement

, , ,

For applications of metal-organic frameworks (MOFs) such as gas storage and separation, flexibility is often seen as a parameter that can tune material performance. In this work we aim to determine the optimal flexibility for the shape selective separation of similarly sized molecules (e.g., Xe/Kr mixtures). To obtain systematic insight into how the flexibility impacts this type of separation, we develop a simple analytical model that predicts a material's Henry regime adsorption and selectivity as a function of flexibility. We elucidate the complex dependence of selectivity on a framework's intrinsic flexibility whereby performance is either improved or reduced with increasing flexibility, depending on the material's pore size characteristics. However, the selectivity of a material with the pore size and chemistry that already maximizes selectivity in the rigid approximation is continuously diminished with increasing flexibility, demonstrating that the globally optimal separation exists within an entirely rigid pore. Molecular simulations show that our simple model predicts performance trends that are observed when screening the adsorption behavior of flexible MOFs. These flexible simulations provide better agreement with experimental adsorption data in a high-performance material that is not captured when modeling this framework as rigid, an approximation typically made in high-throughput screening studies. We conclude that, for shape selective adsorption applications, the globally optimal material will have the optimal pore size/chemistry and minimal intrinsic flexibility even though other nonoptimal materials' selectivity can actually be improved by flexibility. Equally important, we find that flexible simulations can be critical for correctly modeling adsorption in these types of systems.

American Chemical Society2017

Chargement

What Are the Best Materials To Separate a Xenon/Krypton Mixture?

Rocio Mercado, Sondre Kvalvag Schnell, Cory Matthew Simon, Berend Smit

American Chemical Society2015

EPFL2012