Coarse grained modeling of a metal-organic framework/polymer composite and its gas adsorption at the nanoparticle level

TL;DR

This study develops a hybrid coarse-grained simulation approach to model ZIF-8/PVDF composites at the nanoparticle level, revealing how nanoparticle size and shape influence polymer penetration and CO2 adsorption, aligning with experimental observations.

Contribution

A novel hybrid MARTINI/Force Matching force field enables nanoparticle-level simulation of MOF/polymer composites, capturing detailed interfacial interactions and gas adsorption behavior.

Findings

Smaller nanoparticles show less polymer penetration.

All composites exhibit higher CO2 adsorption in the nanoparticle domain.

Polymer conformation varies with nanoparticle size and shape.

Abstract

Simulations have acted as a cornerstone to understand MOF/polymer interface structure, however, no molecular-level simulation has yet been performed at the nanoparticle scale. In this work, a hybrid MARTINI/Force Matching (FM) force field was developed and successfully implemented to model the ZIF-8/PVDF composite at a coarse grained resolution. Inter-phase interactions were modeled using FM potentials, which strive to reasonably reproduce the forces from an atomistic benchmark model, while intraphase interactions are modeled using the general-purpose MARTINI potentials. Systems made of a ZIF-8 nanoparticle embedded into a PVDF matrix were considered to evaluate the effect of nanoparticle size and morphology in the polymer structure and in the CO2 adsorption. Results show that simulations at the nanoparticle level are crucial for depicting the polymer penetration. Notably, the smallest nanoparticle exhibited the least extent of polymer penetration, while the cubic nanoparticle exhibited the highest amount. Polymer conformation and local density values change similarly in all ZIF-8/PVDF systems depending on whether the polymer lies inside or outside of the nanoparticle domain. All composite models present more significant CO2 adsorption in the nanoparticle domain than in the PVDF phase, in agreement with experiments. More remarkably, the small rhombic dodecahedron ZIF-8/PVDF system presents a larger equilibrium amount of gas adsorbed at ambient condition compared to the other two systems, in alignment with the observed polymer penetration trend. On the other hand, the amount of CO2 adsorbed at equilibrium is lower for the rhombic dodecahedron morphology than for the cubic one, contrary to the intuitive expectation founded in the polymer penetration trend. This result could be a reflection of a difference in the number of surface adsorption sites.