# Suggestions for simulating transport through membranes with Molecular Dynamics (MD) and Monte Carlo (MC)?

I'm intending to perform a simulation of the transport of simple ions or water through perm-selective membranes, using Molecular Dynamics (MD) and/or Monte Carlo (MC). I've read around that, in general, with Monte Carlo the polymerization of the membrane can be modeled starting from linear chains and monomers, and get the pseudo-membrane well cross-linked. Once the polymer structure is optimized, MD is performed with explicit solvent, the membrane putted in between two different solutions, so transport properties through membranes can be explored.

I never used MD nor MC. This is where I stand in the knowledge:

Monte Carlo (stochastic approach): Useful to calculate thermodynamic macroscopic properties. The simulation follows the Markov Process, where the steps are independent one from each other (Jan Labanowski).

The stochastic approach, called Monte Carlo, is based on exploring the energy surface by randomly probing the geometry of the molecular system. The iterations are independent of one another (i.e., the system does not contain any "memory''). Is accepted more by physicists than by chemists, probably because MC is not a deterministic method and does not offer time evolution of the system in a form suitable for viewing. (But) for Markov chains, there are efficient methods for deriving time related quantities such as relaxation times (…) Many chemical problems in statistical mechanics are approached more efficiently with MC, and some (e.g., simulations of polymers chains on a lattice) can only be done efficiently with MC.

Molecular dynamics (deterministic approach): To find trajectories and get an insightful movie of the molecular movement, with parameters of interest. Normally the system composes of around 5000 atoms with explicit solvent (water molecules generally). Periodic conditions are employed. The Verlet algorithm and the SHAKE routine are commonly used; here the bond lengths are fixed and the whole molecules are allow to move between certain potential range given (in position and energy), this repeated successively through a minimum path over the potential surface. It's a leap-frog method.

Considering this, I basically don't know where to begin. I would like to know the basic howto's and softwares for MD (GROMACS, CHARMM, LAMMPS, AMBER, etc.) and ways to perform Monte Carlo for polymers. Also, besides the fundamentals, I would like to know if there is like a Handbook for the features of each MD software or modes of application (such us potentials/ functionals used, force fields), or the most important practical keys to be taken into account before simulating.

• I would recommend "Understanding Molecular Simulation [amazon.com/Understanding-Molecular-Simulation-Second-Edition/dp/… ]" by Frenkel and Smit for the theory part Commented Dec 3, 2015 at 18:52
• Another enlightening book is "Computer Simulation of Liquids", by Allen and Tildesley.
– user23061
Commented Dec 4, 2015 at 14:51
• Hi @Orr22, your link seems broken (in the comment above), could you fix that? Commented Jan 20, 2017 at 1:12
• Great thanks for the alert, finally I found this as good introduction amzn.to/2jK5RuS. But, and, especially the documentation of each molecular dynamics software, they are very self-explanatory. Commented Jan 21, 2017 at 14:03
• en.wikipedia.org/wiki/MARTINI Commented May 13, 2017 at 20:04