Theory
Full stack materials design requires understanding the relationship between fine-grained chemical detail and resulting macroscopic properties. The computational workflow at the BioPacific MIP provides a suite of tools to do exactly that.
There are three levels of simulation capabilities in toolbox: all-atom resolution molecular dynamics, coarse-grained molecular dynamics, and field theoretic simulations. The key to making this all work is the careful construction of coarse-grained molecular models that are simultaneously simulatable in field theoretic simulations while as faithful as possible to the underlying all atom molecular dynamics simulations.
Many of the details, nuances, and theoretical considerations of model selection are outlined in the best practices page.
On this page we give a high level overview of the mathematical transformations underlying the development of coarse grained models from all atom simulations, and the field theoretic, continuum simulation of particle based models.
From all atom to coarse-grained models: the relative entropy
Systematic bottom-up coarse-graining refers to the systematic elimination of chemical degrees of freedom. Most commonly, this is effected by replacing collections of atoms with pseudo-atoms, or beads, typically representing the center of mass of the atoms being replaced.
Given a particular mapping from all atom coordinates
The optimal potential of mean force
However, this object is multi-body in nature and intractable to determine. In practice, pairwise potentials are more practical, and several procedures exist to find pairwise approximations of
Schematically, in 1D, the relative entropy can be visualized as thus (image taken from Shell, Adv. Chem. Phys. 2016):
In practice, the minimization of the relative entropy can become a high dimensional optimization problem achieved by solving for:
where
This is handled by the in-house coarse graining simulation software at BioPacific MIP.
From coarse-grained particle models to continuum: field theoretic models
There are several exact mathematical transformations from particle to continuum descriptions. The one employed in the BioPACIFIC simulation software suite is the auxiliary field approach.
A good series of lectures on field theory can be found courtesy of Professor Rob Riggleman.
In short, the first step is to define particle density
where the sum is over all particles in the system of species
This then allows for the critical mathematical trick of using the Hubbard Stratonovich transformation, an exact transformation that allows for the transformation of any pairwise, positive definite interaction
Definint the shorthand:
we express the Hubbard Stratonovich transformation:
where we have introduced functional integrals over a newly introduced "auxiliary" field
Note that the density fields only appear linearly in the argument of the exponential, or the new "effective Hamiltonian". This decouples all configuration integrals from one another and is what allows field theories to efficiently treate the configurational space of large molecules. Analytic continuation of the functional integral over
The field theory enables full field theoretic simulation and useful approximations like self-consistent field theory. It also provides helpful analytical tools for understanding polymeric behavior in disordered solutions, using mean-field analysis and the random phase approximation.