Understanding the phase behavior of polyelectrolyte coacervation
is crucial for many applications, including consumer formulations,
wet adhesives, processed food, and drug delivery. However, in most
cases, modeling coacervation is not easily accessed by molecular simulation
methods due to the long-range nature of electrostatic forces and the
typically high molecular weights of the species involved. We present
a modeling strategy to study complex coacervation leveraging the strengths
of both particle simulations and polymer field theory. Field theory
is uniquely suited to capture larger-length scales that are inaccessible
to particle simulations, but its predictive capability is limited
by the need to specify emergent parameters. Using model coacervate-forming
systems consisting of poly(acrylic acid) and poly(allylamine hydrochloride),
we show an original way to use small-scale, all-atom simulations to
parameterize field-theoretic models via the relative entropy coarse-graining
approach. The dependence of coacervation on the salt concentration,
molecular weight, and charge stoichiometry is predicted without fitting
to experimental data and is consistent with experimental trends including
asymmetric phase behavior from non-stoichiometric mixtures of polyelectrolytes.
This demonstrates a unique simulation approach to study phase behavior
in coacervate-forming systems, which is particularly useful when chemical
specificity is of interest.