Modelling Diffusion in Zeolites by Molecular Dynamics Simulations

Abstract

The diffusion of molecules sorbed in zeolites is of growing interest for understanding the mechanisms of chemical processes with regard to selectivity and reactivity [1].

MD simulations give insight into physical systems on the molecular level allowing to study and visualize the motion of molecules even beyond the possibilities of experiments [2,3]. Single system parameters can easily be varied to study their influence, also those parameters that are fixed in reality (e.g., the size of particles). We present a cross section of our recent work to illustrate the capabilities of MD: The self diffusion coefficients (D) of a mixture of methane and xenon in silicalite show remarkable deviations from those of the pure species. This is shown and confirmed by PFG NMR experiments [4].

Simulating ethane in zeolite A the mechanism of diffusion has been studied. The effects of rotation on the diffusion lead to cases where D decreases with growing temperature [5].

The independence of self diffusion on lattice vibrations is proven even for zeolites with windows of guest particle size comparing simulations with rigid and vibrating zeolite lattice [6].     

Dissipative particle dynamics simulation of the soft micro actuator using polymer chain displacement in electro-osmotic flow

ABSTRACT

In this research, the numerical simulation of a soft polymer micro actuator performance has been investigated using the dissipative particle dynamics method in electro-osmotic flow. Effective factors including electro-osmotic flow and polymer chain parameters have been studied. First of all, considering a wide range of electro-osmotic parameters, the validation of analytical results is carried out in a simple micro channel. The electric field and zeta potential changes are linearly related to the flow rate, and the kh parameter behaves nonlinearly to around the kh?=?10. In the following, a convergent–divergent channel is used for the soft micro actuator simulation in which a polymer chain as a heart of actuation is embedded in the middle. As the main control parameter, the direction of the electric field is changed every 4?s, and it leads to a reciprocating motion. The numerical results indicate that the displacement of the soft polymer chain will be increased by enhancing the electric field, the number of beads, decreasing the harmonic bond coefficient and also exposing more length of a polymer chain in front of fluid flow. The results of this study may be useful for some future applications such as artificial fibres and muscles.     

Hydrogen bond analysis of confined water in mesoporous silica using the reactive force field

ABSTRACT

The structural and dynamical properties of water confined in nanoporous silica with a pore diameter of 2.7?nm were investigated by performing large-scale molecular dynamics simulations using the reactive force field. The radial distribution function and diffusion coefficient of water were calculated, and the values at the centre of the pore agreed well with experimental values for real water. In addition, the pore was divided into thin coaxial layers, and the average number of hydrogen bonds, hydrogen bond lifetime and hydrogen bond strength were calculated as a function of the radial distance from the pore central axis. The analysis showed that hydrogen bonds involving silanol (Si–OH) have a longer lifetime, although the average number of hydrogen bonds per atom does not change from that at the pore centre. The longer lifetime, as well as smaller diffusion coefficient, of these hydrogen bonds is attributed to their greater strength.     

Configurational Temperature for Brownian Dynamics

Abstract

Expressions for the configurational temperature are evaluated in Brownian dynamics simulations of a Lennard-Jones fluid and compared with the input temperature which is used to generate the random displacements. It is found that the two temperatures agree in the limit of large numbers of particles, and even for moderate system sizes the configurational temperature is a useful check on the correctness of the simulation algorithm. Investigation of the autocorrelation functions shows that for Lennard-Jones and power-law fluids, the correlation time of the configurational temperature is shorter than other typical thermodynamic quantities, and it generally increases with the range of the potential.     

Molecular Dynamics Simulations of Proteins in Water Without the Truncation of Long-range Coulomb Interactions

A new program package (COSMOS90) for molecular dynamics simulations was developed to simulate large molecular systems consisting of more than tens of thousands of atoms without the truncation of long-range coulomb interactions. This program package was based on a new approximation scheme (PPPC) for calculating efficiently the coulomb interactions without sacrificing accuracy. In this approximation scheme, the group of charges at a long distance from each atom was represented by a total charge and total dipole moment of the group. In order to assess the accuracy of PPPC and the ability of COSMOS90, molecular dynamics simulations were carried out for a large system consisting of 16108 atoms (human lysozyme in water) for 50 ps using this program package. The coulomb energy per solute atom was calculated with only five percent of the error found in the 10 Å cut-off approximation (about 0.9 kcal/mol versus 18 kcal/mol, respectively). The molecular dynamics simulations using COSMOS90 require no more CPU time than the simulations based on the 10 Å cut-off approximation of the conventional programs for macromolecular simulations.     

Molecular dynamics study on perfect and defective graphene/calcium-silicate-hydrate composites under tensile loading

ABSTRACT

Atomic models of graphene/calcium-silicate-hydrate (G/C-S-H) are constructed by embedding perfect or defective graphene in molecular structures of amorphous C-S-H. Molecular dynamics (MD) simulation is utilised to study mechanical properties of the G/C-S-H and the enhancing effect of perfect and defective graphene is compared. The effects of temperature and strain rate on perfect and defective G/C-S-H are also investigated and compared. The results from present simulations show that (i) the defective graphene has better enhancing effect in C-S-H than perfect one and it grows with the increase of defect sizes; (ii) the tensile strength of G/C-S-H decreases with the increase of temperature and the defective G/C-S-H is more susceptible to temperature than the perfect one; (iii) the ultimate strength and the failure strain increase significantly with the increase of strain rate and the effects of strain rate on perfect and defective G/C-S-H are similar. These findings provide important atomic insights for understanding the mechanical behaviours of G/C-S-H composite.     

Molecular dynamics simulation of human interleukin-4: comparison with NMR data and effect of pH,counterions and force field on tertiary structure stability

The human protein interleukin-4 (IL-4) has been simulated at two different pH values 2 and 6, with different amounts of counterions present in the aqueous solution, and with two different force-field parameter sets using molecular dynamics simulation with the aim of validation of force field and simulation set-up by comparison to experimental nuclear magnetic resonance data, such as proton–proton nuclear Overhauser effect (NOE) distance bounds, ³ J(HN,HC_α) coupling constants and backbone N–H order parameters. Thirteen simulations varying in the length from 3 to 7 ns are compared.

At pH 6 both force-field parameter sets used do largely reproduce the NOE's and order parameters, the GROMOS 45A3 set slightly better than the GROMOS 53A6 set. ³ J values predicted from the simulation agree less well with experimental values. At pH 2 the protein unfolds, unless counterions are explicitly present in the system, but even then the agreement with experiment is worse than at pH 6. When simulating a highly charged protein, such as IL-4 at pH 2, the inclusion of counterions in the simulation seems mandatory.     

A Leap-frog Algorithm for Stochastic Dynamics

Abstract

A third-order algorithm for stochastic dynamics (SD) simulations is proposed, identical to the powerful molecular dynamics leap-frog algorithm in the limit of infinitely small friction coefficient γ. It belongs to the class of SD algorithms, in which the integration time step Δt is not limited by the condition Δt ≤ γ^?1, but only by the properties of the systematic force. It is shown how constraints, such as bond length or bond angle constraints, can be incorporated in the computational scheme. It is argued that the third-order Verlet-type SD algorithm proposed earlier may be simplified without loosing its third-order accuracy. The leap-frog SD algorithm is proven to be equivalent to the verlet-type SD algorithm. Both these SD algorithms are slightly more economical on computer storage than the Beeman-type SD algorithm.     

Dissipative particle dynamics of non-spherical particles using a Gaussian density model

The Gaussian density molecular model has been adapted for dissipative particle dynamics. The model, when combined with a soft potential, is shown to be a very flexible mesoscale model exhibiting a wide range of phase behaviour. The soft potential allows relatively large time steps to be used and hence a more rapid equilibration. In addition, the model can be used to study both uniaxial and biaxial systems. We have undertaken a number of pilot studies and have demonstrated that the Gaussian model is able to identify nematic–isotropic phase transitions in liquid crystals and the formation of ordered discotic phases.