de Pablo Group

Our group investigates the physics and thermodynamics of complex materials using statistical mechanics, molecular simulations, and machine learning.  Using the results, we design new systems for technological applications.


Principal Investigator

Juan de Pablo

Sequence blockiness controls the structure of polyampholyte necklaces

Rumyantsev, A.M., Johner, A. and de Pablo, J.J., 2021. Sequence Blockiness Controls the Structure of Polyampholyte Necklaces. ACS macro letters, 10, pp.1048-1054.

Polyelectrolyte Complex Coacervation across a Broad Range of Charge Densities

Neitzel, A.E., Fang, Y.N., Yu, B., Rumyantsev, A.M., de Pablo, J.J. and Tirrell, M.V., 2021. Polyelectrolyte complex coacervation across a broad range of charge densities. Macromolecules, 54(14), pp.6878-6890.

Scaling theory of neutral sequence-specific polyampholytes

Rumyantsev, A.M., Jackson, N.E., Johner, A. and De Pablo, J.J., 2021. Scaling Theory of Neutral Sequence-Specific Polyampholytes. Macromolecules, 54(7), pp.3232-3246.

Harnessing Peptide Binding to Capture and Reclaim Phosphate

Fowler, Whitney C., et al. "Harnessing Peptide Binding to Capture and Reclaim Phosphate." Journal of the American Chemical Society 143.11 (2021): 4440-4450. Whitney C. Fowler, Chuting Deng, Gabriella M. Griffen, Tess Teodoro, Ashley Z. Guo, Michal Zaiden, Moshe Gottlieb*, Juan J. de Pablo, Matthew V. Tirrell

Role of Molecular Architecture on Ion Transport in Ethylene oxide-Based Polymer Electrolytes

Deng, Chuting, et al. "Role of Molecular Architecture on Ion Transport in Ethylene oxide-Based Polymer Electrolytes." Macromolecules 54.5 (2021): 2266-2276. Chuting Deng, Michael A. Webb, Peter Bennington, Daniel Sharon, Paul F. Nealey, Shrayesh N. Patel, Juan J. de Pablo

Role of solvation site segmental dynamics on ion transport in ethylene-oxide based side-chain polymer electrolytes

Bennington, Peter, et al. "Role of solvation site segmental dynamics on ion transport in ethylene-oxide based side-chain polymer electrolytes." Journal of Materials Chemistry A 9.15 (2021): 9937-9951. Peter Bennington, Chuting Deng, Daniel Sharon, Michael A. Webb, Juan J. de Pablo, Paul F. Nealey, Shrayesh N. Patel

Modeling the Binding Mechanism of Remdesivir, Favilavir, and Ribavirin to SARS-CoV-2 RNA-Dependent RNA Polymerase

Byléhn, F., Menéndez, C.A., Perez-Lemus, G.R., Alvarado, W. and De Pablo, J.J., 2021. Modeling the binding mechanism of remdesivir, favilavir, and ribavirin to SARS-CoV-2 RNA-dependent RNA polymerase. ACS central science, 7(1), pp.164-174.

Polyelectrolyte complex coacervates: Recent developments and new frontiers

Rumyantsev, A.M., Jackson, N.E. and De Pablo, J.J., 2021. Polyelectrolyte Complex Coacervates: Recent Developments and New Frontiers. Annual Review of Condensed Matter Physics, 12, pp.155-176.

Molecular mass dependence of interfacial tension in complex coacervation

Audus, D.J., Ali, S., Rumyantsev, A.M., Ma, Y., de Pablo, J.J. and Prabhu, V.M., 2021. Molecular Mass Dependence of Interfacial Tension in Complex Coacervation. Physical Review Letters, 126(23), p.237801.

Complex coacervation of statistical polyelectrolytes: Role of monomer sequences and formation of inhomogeneous coacervates

Yu, B., Rumyantsev, A.M., Jackson, N.E., Liang, H., Ting, J.M., Meng, S., Tirrell, M.V. and de Pablo, J.J., 2021. Complex coacervation of statistical polyelectrolytes: role of monomer sequences and formation of inhomogeneous coacervates. Molecular Systems Design & Engineering.

View All Publications

Liquid crystals (LCs) are a phase of matter that flows like a liquid, but the orientations of the molecules are highly ordered over a very long range. This presence of long-range orientation results in interesting behavior of systems that employ LCs. In our group, we model LCs on multiple scales in an effort to engineer new applications for the laboratory and industry. At the atomistic level, we investigate the behavior of LCs near surfaces to determine the types and strength of anchoring present at different surfaces. At a mesoscale, we study systems mixtures and determine the accessibility different phases of LCs. On the largest scales, we investigate the behavior of particles, from the nanometer to micron scale, and observe their behavior in an LC solvent; the presence of defects in the LC has a marked effect on particle behavior, so by controlling the defect with fields (flow, electric, magnetic,etc.), we can dictate particle behavior in a well controlled manner.  We also study how introducing active biological agents, such as bacteria or myosin motors, influences the dynamics of LC systems.

Our research group uses coarse grain models to study the biophysics of DNA and chromatin. Recent efforts have included:

  1. Developing hierarchical coarse-grained models to study DNA and chromatin at multiple length scales.
  2. Relating nucleosome interactions to chromatin fiber structure
  3. Exploring collective motions in chromatin using nonlinear manifold learning
  4. Examining the effect of DNA-protein interactions on the structure and dynamics of chromatin

With the help of modern polymer chemistry, macromolecules can now be synthesized with exquisite control over architecture, composition, and monomer sequence. These parameters work cooperatively to control the dynamics, self-assembly, and morphology of the resulting materials and determine how they may be applied in new technologies.  To understand and predict these properties, we develop and utilize multi-scale Monte Carlo and molecular dynamics simulations as well as novel theoretical treatments.  Our current ares of research include:

  • Controlling complex coacervation of polyelectrolytes with random charge sequences
  • Structure and rheology of "miktoarm" polymers and their blends
  • Using multi-block co-polymers to design complex nanostructures in solution
  • Dynamics of interlocking polymers


Molecular simulations are typically limited by the time scale of sampling. With reasonable amount of computational resources one can only simulate on the order of hundreds of nanoseconds. On the other hand in real complex systems most of the phenomena of interest occur at orders-of-magnitude longer time scales. To solve such problems, we develop new advanced sampling algorithms that can accelerate discovery in these systems in both Monte Carlo (MC) and Molecular Dynamics (MD) simulations, which have been incorporated into the Software Suite for Advanced General Ensemble Simulations (SSAGES) code.

Many soft matter systems exhibit a wide range of length scales, making them challenging to simulate and study.  For example, colloidal suspensions are naturally represented as discrete particles embedded in a continuum since the colloids are much larger than the individual molecules of the solvent.  To couple the dynamics of particles and continua, we have developed the Continuum-Particle Simulation Software (COPSS) code, which uses innovative numerical methods to efficiently simulate mesoscopic soft matter systems.

The SSAGES and COPSS codes are available at:


Machine learning provides a powerful set of methods for extracting useful information from large and complex data sets, such as those routinely generated by molecular simulations.  Our group uses these two complementary tools to gain insight into the properties of molecular systems and to make predictions for materials design.  Recent efforts include:

  • Using artificial neural networks (ANNs) for force-biasing during molecular dynamics simulations
  • Predicting electronic structure properties of conjugated polymers at coarse-grained resolution using ANNs
  • Inverse design of materials based on desired physical properties