Molecular simulation adsorption models

Molecular simulation methods can be a complement to surface complexation modeling on metal-bacteria adsorption reactions, which provides a more detailed and atomistic information of how metal cations interact with specific functional groups within bacterial cell wall. Johnson et al., (2006) applied molecular dynamics (MD) simulations to analyze equilibrium structures, coordination bond distances of metal-ligand complexes. 

It is important to propose molecular and theoretical models to describe the forces, energy, structure and dynamics of water near mineral surfaces. Our understanding of experimental results concerning hydration forces, the hydrophobic effect, swelling, reaction kinetics and adsorption mechanisms in aqueous colloidal systems is rapidly advancing as a result of recent Monte Carlo (MC) and molecular dynamics (MO) models for water properties near model surfaces. This paper reviews the basic MC and MD simulation techniques, compares and contrasts the merits and limitations of various models for water-water interactions and surface-water interactions, and proposes an interaction potential model which would be useful in simulating water near hydrophilic surfaces. In addition, results from selected MC and MD simulations of water near hydrophobic surfaces are discussed in relation to experimental results, to theories of the double layer, and to structural forces in interfacial systems. 

In this review, we introduce another approach to study the multiscale structures of polymer materials based on a lattice model. We first show the development of a Helmholtz energy model of mixing for polymers based on close-packed lattice model by combining molecular simulation with statistical mechanics. Then, holes are introduced to account for the effect of pressure. Combined with WDA, this model of Helmholtz energy is further applied to develop a new lattice DFT to calculate the adsorption of polymers at solid-liquid interface. Finally, we develop a framework based on the strong segregation limit (SSL) theory to predict the morphologies of micro-phase separation of diblock copolymers confined in curved surfaces. 

To establish the molecular thermodynamic model for uniform systems based on concepts from statistical mechanics, an effective method by combining statistical mechanics and molecular simulation has been recommended (Hu and Liu, 2006). Here, the role of molecular simulation is not limited to be a standard to test the reliability of models. More directly, a few simulation results are used to determine the analytical form and the corresponding coefficients of the models. It retains the rigor of statistical mechanics, while mathematical difficulties are avoided by using simulation results. The method is characterized by two steps (1) based on a statistical-mechanical derivation, an analytical expression is obtained first. The expression may contain unknown functions or coefficients because of mathematical difficulty or sometimes because of the introduced simplifications. (2) The form of the unknown functions or unknown coefficients is then determined by simulation results. For the adsorption of polymers at interfaces, simulation was used to test the validity of the weighting function of the WDA in DFT. For the meso-structure of a diblock copolymer melt confined in curved surfaces, we found from MC simulation that some more complex structures exist. From the information provided by simulation, these complex structures were approximated as a combination of simple structures. Then, the Helmholtz energy of these complex structures can be calculated by summing those of the different simple structures. 

As expected, the total interaction energies depend strongly on the van der Waals radii (of both sorbate and sorbent atoms) and the surface densities. This is true for both HK type models (Saito and Foley, 1991 Cheng and Yang, 1994) and more detailed statistical thermodynamics (or molecular simulation) approaches (such as Monte Carlo and density functional theory). Knowing the interaction potential, molecular simulation techniques enable the calculation of adsorption isotherms (see, for example, Razmus and Hall, (1991) and Cracknell etal. (1995)). 

Molecular simulation of the adsorption of gases by the ALPOs was pioneered by Cracknell and Gubbins (1993), who pointed out that the aluminophosphates should be easier to model than the aluminosilicates. There are two important advantages first, the charge neutrality of the framework means that there are no exchangeable cations to be taken into account (this is, of course, also true for pure Silicalite) and second, the modelling is simpler because the pores are unidirectional with no interconnections. 

The Horvath-Kawazoe (HK) method is capable of generating model isotherms more efficiently than either molecular simulation (MS) or density functional theory (DFT) to characterize the pore size distribution (PSD) of microporous solids. A two-stage HK method is introduced that accounts for monolayer adsorption in mesopores prior to capillary condensation. PSD analysis results from the original and two-stage HK models are evaluated. 

In a combined simulation and experimental project we set out to assess the adsorption properties of a series of zeolites. In the present work the adsorption properties of n-butane and iso-butane on MFI are being studied. The experimental part consists in the validation of the molecular simulation model, by confirming its results. The experiments were performed in a constructed in-house manometric apparatus coupled with a NIR spectrometer (Perkin Elmer, FT-IR system, GX Spectrum). Figure 1 is a scheme of the experimental set-up. 

Fig.7 Adsorption equilibrium for IPA-TCE acid site model system (molecular simulation)...

In this study, equilibria and isosteric heat of adsorption for the system of chlorinated hydrocarbons and Y-type zeolite were obtained with gravimetric method and chromatographic method. By comparing an experiment result with a molecular simulation result, the validity of forcefield parameters and zeohte model was exartuned... 

Figure 3 Comparison of PSDs obtained using the Dubinin-Stoeckli (DS), Horvalh-Kawazoe (HK), and density fitnctional theory (DFT) methods to interpret an isotherm generated from molecular simulation of nitrogen adsorption in a model carbon that has an Gaussian distribution of slit pore widths (18]. Results are shown for mean pore widths of 8.9 A (left) and 16.9 A (right).

Figure 6 PSDs obtained for methane adsorption in square model carbon pores using molecular simulation to interpret an activated carbon isotherm. PSD results are shown for regularization smoothing parameter values of 1 (solid line), 10 (open circles), 100 (open diamonds), 600 (filled circles), and 800 (filled diamonds) [25].

In this paper, a modified HK method is presented which accounts for spatial variations in the density profile of a fluid (argon) adsorbed within a carbon slit pore. We compare the pore width/filling pressure correlations predicted by the original HK method, the modified HK method, and methods based upon statistical thermodynamics (density functional theory and Monte Carlo molecular simulation). The inclusion of the density profile weighting in the HK adsorption energy calculation improves the agreement between the HK model and the predictions of the statistical thermodynamics methods. Although the modified Horvath-Kawazoe adsorption model lacks the quantitative accuracy of the statistical thermodynamics approaches, it is numerically convenient for ease of application, and it has a sounder molecular basis than analytic adsorption models derived from the Kelvin equation. 

The principal drawback of the DFT method is that it is computationally intensive relative to the classical adsorption models, although it is still much less compute-intensive than full Monte Carlo molecular simulation. A semianalytic adsorption model that retains computational efficiency while accounting for gas-solid potential interactions in micropores was originally proposed by Horvath and Kawazoe [12], In the Horvath-Kawazoe or HK method, a pore filling correlation is obtained by calculating the mean heat of adsorption (/> required to transfer an adsorbate molecule from the gas phase to the condensed phase in a slit pore of width // ... 

The basis for thermodynamic calculations is the adsorption isotherm, which gives the amount of gas adsorbed in the nanopores as a function of the external pressure. Adsorption isotherms are measured experimentally or calculated from theory using molecular simulations. Potential functions are used to constmct a detailed molecular model for atom-atom interactions and a distribution of point charges is used to reproduce the polarity of the solid material and the adsorbing molecules. Recently, ab initio quantum chemistry has been applied to the theoretical determination of these potentials, as discussed in another chapter of this book. 

---