skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Effects of Defects and Presence of Open-Metal Sites on the Structure and Dynamics of Water in Hydrophobic Zeolitic Imidazolate Frameworks
Most of the chemistry in nanoporous materials with small pore sizes and windows takes place on the outer surface, which is in direct contact with the substrate/solvent, rather than within the pores and channels. Here, we report the results of our comprehensive atomistic molecular dynamics (MD) simulations to decipher the interaction of water with a realistic finite ∼5.1 nm nanoparticle (NP) model of ZIF-8, with edges containing undercoordinated Zn metal sites, vs a conventionally employed pristine crystalline bulk (CB) model. The hydrophobic interior surface of the CB model imparts significant dynamical behavior on water molecules with (i) increasing diffusivity from the surface toward the center of the pores and (ii) confined water, at low concentration, showing similar diffusivity to that of the bulk water. On the other hand, water molecules adsorbed on the surface of the NP model exhibit a range of characteristics, including “coordinated”, “confined”, and “bulk-like” behavior. Some of the water molecules form coordinative bonds with the undercoordinated Zn metal centers and act as nucleation sites for the water droplets to form, facilitating diffusion into the pores. However, diffusion of water molecules is limited to the areas near the surface and not all the way to the core of the NP model. Our atomistic MD simulations provide insights into the stability of ZIFs in aqueous solutions despite hydrolysis of their outer surface. Such insights are helpful in designing more robust nanoporous materials for applications in humid environments.  more » « less
Award ID(s):
2302617 2302618
PAR ID:
10533704
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
American Chemical Society
Date Published:
Journal Name:
Journal of Chemical Information and Modeling
Volume:
63
Issue:
22
ISSN:
1549-9596
Page Range / eLocation ID:
7097 to 7106
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Mismatched complex oxide thin films and heterostructures have gained significant traction for use as electrolytes in intermediate temperature solid oxide fuel cells, wherein interfaces exhibit variation in ionic conductivity as compared to the bulk. Although misfit dislocations present at interfaces in these structures impact ionic conductivity, the fundamental mechanisms responsible for this effect are not well understood. To this end, a kinetic lattice Monte Carlo (KLMC) model was developed to trace oxygen vacancy diffusion at misfit dislocations in SrTiO3/BaZrO3 heterostructures and elucidate the atomistic mechanisms governing ionic diffusion at oxide interfaces. The KLMC model utilized oxygen vacancy migration energy barriers computed using molecular statics. While some interfaces promote oxygen vacancy diffusion, others impede their transport. Fundamental factors such as interface layer chemistry, misfit dislocation structure, and starting and ending sites of migrating ions play a crucial role in oxygen diffusivity. Molecular dynamics (MD) simulations were further performed to support qualitative trends for oxygen vacancy diffusion. Overall, the agreement between KLMC and MD is quite good, though MD tends to predict slightly higher conductivities, perhaps a reflection of nuanced structural relaxations that are not captured by KLMC. The current framework comprising KLMC modeling integrated with molecular statics offers a powerful tool to perform mechanistic studies focused on ionic transport in thin film oxide electrolytes and facilitate their rational design. 
    more » « less
  2. Metallic zinc (Zn) has been considered one of the most promising anode materials for next-generation aqueous Zn batteries due to its low redox potential and high storage capacity. However, excessive dendrite formation in Zn metal, corrosion, the evolution of hydrogen gas during the cycling process, and the poor Zn-ions (Zn2+) transport from the electrolyte to the electrode limits its practical application. One of the most effective strategies to suppress Zn dendrite growth and promote Zn2+ transport is to introduce suitable protective layers between the Zn metal electrode and the electrolyte. Herein, we mathematically simulated the dynamic interactions between the Zn deposition on the anode and the resulting displacement of a protective layer that covers the anode, the latter of which can simultaneously inhibit Zn dendrite growth and enhance the Zn2+ transport through the interface between Zn anode and the protective layer. Our simulation results indicate that a protective layer of high Zn2+ diffusivity not only improves the deposition rate of the Zn metal but also prevents the dendrite growth by homogenizing the Zn2+ concentration at the anode surface. In addition, it is revealed that the anisotropic Zn2+ diffusivity in the protective layer influences the 2D diffusion of Zn2+. Higher Zn2+ diffusivity perpendicular to the Zn metal surface inhibits the dendrite growth, while higher diffusivity parallel to the Zn metal surface promotes dendrite growth. Our work thus provides a fundamental understanding and a design principle of controlling anisotropic Zn2+ diffusion in the protective layer for better suppression of dendrite growth in Zn metal batteries. 
    more » « less
  3. Abstract Micro- and nanoporous materials have gathered attention from the scientific community due to their size dependent properties, including but not limited to high specific surface area, surface diffusivity, bulk diffusivity and permeability, catalytic activity, and distinct optical properties. In this work, spherical nanoporous copper (np-Cu) powders, due to their nanosized porosity and low Cu2O content, show hemispherical total reflectance of 20% which is significantly lower than its bulk counterpart value for solid or molten copper of approximately 97% at wavelengths of most commercial Laser Powder Bed Fusion (L-PBF) commercial machines. The low-reflectance of np-Cu powders has the potential to be used in L-PBF to improve laser absorption, volumetric energy efficiency, and throughput of this additive manufacturing process. In fact, a prepared mixture of solid Cu powders containing only 5 wt.% of np-Cu powders reflects 34.8 % less than pure copper powders as shown in this paper. Np-Cu powders are fabricated via chemical dealloying of gas atomized CuAl alloy in a robust and scalable approach, and then mixed with pure copper powders to prepare hybrid feedstocks. Under this framework, the crucial role of deglomeration strategies to achieve homogeneity and flowability of np-Cu/Cu hybrid mixtures are evaluated via particle imaging to determine agglomerate size and composition with an eye at obtaining a high-quality print in L-PBF. In np-Cu powders fabrication, washing them in low-surface tension fluids upholds the highest degree of deglomeration in their fabrication process, and for hybrid feedstocks preparation, pre-mixing Cu and CuAl prior to dealloying yields the best homogeneity results with smallest size of agglomerates and good flowability. 
    more » « less
  4. Nanoporous carbons play an important role in different electrochemical applications such as being utilized as electrodes in supercapacitors. Application of electric potential to a porous electrode in electrolyte solution stimulates adsorption or desorption of ions on the electrode surface. Electrosorption causes appearance of solvation pressure in the pores and results in electrode deformation. In this work, using molecular dynamics simulations and the continuum theory based on the modified Poisson-Boltzmann equation, we studied the structure of the electrical double layer in slit graphitic micropores filled with a NaCl aqueous solution, and solvation pressure in these pores. We focused on the behavior of the solvation pressure as a function of pore width and surface charge density. Within molecular dynamics simulations, two different water models were used -- an explicit model based on SPC/E water molecules and an implicit model, i.e., structureless background with fixed dielectric permittivity. The latter allows us to relate molecular dynamics simulations to the continuum theory. Simulations with explicit water show a qualitatively different behavior of the solvation pressure in the 1 and 2 nm pores as a function of the surface charge density. We demonstrated that the value of the solvation pressure is defined by a delicate balance between Van der Waals and electrostatic contributions. We demonstrated that the theory predicts the dependence of the solvation pressure on the pore width, which matches the results of simulations using the implicit water model. Finally, we adapted the continuum theory, developed for adsorption-induced deformation to estimate the deformation of a carbon electrode due to electrosorption. Our results can be used in the further development of nanoporous actuators working based on electrosorption-induced deformation. 
    more » « less
  5. Governing partial differential equations (PDEs) play a critical role in materials research and applications, as they describe essential physics underlying materials behaviour. Traditionally, these equations are developed through phenomenological modelling of experimental results or first principle analysis based on conservation laws. In addition, molecular dynamics (MD) simulations capture atomistic-scale behaviour with detailed physics. However, translating atomistic insights into continuum-scale governing equations remains a significant challenge. Empowered by recent advances in data-driven modelling, we develop a computational framework to learn governing PDEs directly from atomistic simulation data. The framework integrates numerical differentiation of MD data with the identification of constitutive relationships. It proves effective and efficient in learning governing PDEs from noisy and limited MD datasets, without requiring prior knowledge of the final PDEs. Using this framework, we identify a nonlinear PDE governing solid-state diffusion in nickel–hydrogen alloys. This PDE reveals a highly concentration-dependent diffusivity that varies over an order of magnitude. Our data-driven computational framework paves the way for cross-scale constitutive modelling. 
    more » « less