skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Exact analytical solution of the Flory–Huggins model and extensions to multicomponent systems
The Flory–Huggins theory describes the phase separation of solutions containing polymers. Although it finds widespread application from polymer physics to materials science to biology, the concentrations that coexist in separate phases at equilibrium have not been determined analytically, and numerical techniques are required that restrict the theory’s ease of application. In this work, we derive an implicit analytical solution to the Flory–Huggins theory of one polymer in a solvent by applying a procedure that we call the implicit substitution method. While the solutions are implicit and in the form of composite variables, they can be mapped explicitly to a phase diagram in composition space. We apply the same formalism to multicomponent polymeric systems, where we find analytical solutions for polydisperse mixtures of polymers of one type. Finally, while complete analytical solutions are not possible for arbitrary mixtures, we propose computationally efficient strategies to map out coexistence curves for systems with many components of different polymer types.  more » « less
Award ID(s):
2011750 2245850
PAR ID:
10585836
Author(s) / Creator(s):
;
Publisher / Repository:
J. Chemical Physics
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
161
Issue:
4
ISSN:
0021-9606
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We introduce a lattice framework that incorporates elements of Flory–Huggins solution theory and the q-state Potts model to study the phase behavior of polymer solutions and single-chain conformational characteristics. Without empirically introducing temperature-dependent interaction parameters, standard Flory–Huggins theory describes systems that are either homogeneous across temperatures or exhibit upper critical solution temperatures. The proposed Flory–Huggins–Potts framework extends these capabilities by predicting lower critical solution temperatures, miscibility loops, and hourglass-shaped spinodal curves. We particularly show that including orientation-dependent interactions, specifically between monomer segments and solvent particles, is alone sufficient to observe such phase behavior. Signatures of emergent phase behavior are found in single-chain Monte Carlo simulations, which display heating- and cooling-induced coil–globule transitions linked to energy fluctuations. The framework also capably describes a range of experimental systems. Importantly, and by contrast to many prior theoretical approaches, the framework does not employ any temperature- or composition-dependent parameters. This work provides new insights regarding the microscopic physics that underpin complex thermoresponsive behavior in polymers. 
    more » « less
  2. Polymer/ionic liquid systems are being increasingly explored, yet those exhibiting lower critical solution temperature (LCST) phase behavior remain poorly understood. Poly(benzyl methacrylate) in certain ionic liquids constitute unusual LCST systems, in that the second virial coefficient (A2) in dilute solutions has recently been shown to be positive, indicative of good solvent behavior, even above phase separation temperatures, where A2 < 0 is expected. In this work, we describe the LCST phase behavior of poly(benzyl methacrylate) in 1-butyl-1-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide for three different molecular weights (32, 63, and 76 kg/mol) in concentrated solutions (5–40% by weight). Turbidimetry measurements reveal a strong concentration dependence to the phase boundaries, yet the molecular weight is shown to have no influence. The critical compositions of these systems are not accessed, and must therefore lie above 40 wt% polymer, far from the values (ca. 10%) anticipated by Flory-Huggins theory. The proximity of the experimental cloud point to the coexistence curve (binodal) and the thermo-reversibility of the phase transitions, are also confirmed at various heating and cooling rates. 
    more » « less
  3. Phase separation in multicomponent mixtures is of significant interest in both fundamental research and technology. Although the thermodynamic principles governing phase equilibria are straightforward, practical determination of equilibrium phases and constituent compositions for multicomponent systems is often laborious and computationally intensive. Here, we present a machine-learning workflow that simplifies and accelerates phase-coexistence calculations. We specifically analyze capabilities of neural networks to predict the number, composition, and relative abundance of equilibrium phases of systems described by Flory-Huggins theory. We find that incorporating physics-informed material constraints into the neural network architecture enhances the prediction of equilibrium compositions compared to standard neural networks with minor errors along the boundaries of the stable region. However, introducing additional physics-informed losses does not lead to significant further improvement. These errors can be virtually eliminated by using machine-learning predictions as a warm-start for a subsequent optimization routine. This work provides a promising pathway to efficiently characterize multicomponent phase coexistence. 
    more » « less
  4. Phase separation in multicomponent mixtures is of significant interest in both fundamental research and technology. Although the thermodynamic principles governing phase equilibria are straightforward, practical determination of equilibrium phases and constituent compositions for multicomponent systems is often laborious and computationally intensive. Here, we present a machine-learning workflow that simplifies and accelerates phase-coexistence calculations. We specifically analyze capabilities of neural networks to predict the number, composition, and relative abundance of equilibrium phases of systems described by Flory-Huggins theory. We find that incorporating physics-informed material constraints into the neural network architecture enhances the prediction of equilibrium compositions compared to standard neural networks with minor errors along the boundaries of the stable region. However, introducing additional physics-informed losses does not lead to significant further improvement. These errors can be virtually eliminated by using machine-learning predictions as a warm-start for a subsequent optimization routine. This work provides a promising pathway to efficiently characterize multicomponent phase coexistence. 
    more » « less
  5. The excess free energy of mixing Δ G ex governs the phase behavior of mixtures and controls material properties. It is challenging, however, to measure Δ G ex in simulations. Previously, we developed a method that combines molecular dynamics (MD) simulations with thermodynamic integration along the path of transformation of chains to predict the Flory Huggins interaction parameter χ for polymer mixtures and block copolymers. However, this method is best applied when the constituent molecules of the blends are structurally related. To overcome this limitation, we have developed a new method to predict Δ G ex for mixtures. We perform simulations to induce phase separation within a mixture by gradually weakening the interaction between different species. To compute Δ G ex we measure the thermodynamic work required to modify the interactions and the interfacial energy between the separated phases. We validate our method by applying it first to equimolar mixtures of labeled and unlabeled Lennard-Jones (LJ) beads, and labeled and unlabeled benzene, which results in good agreement with ideal solution theory. Then we compute the excess free energy of mixing for equimolar mixtures of benzene and pyridine, using both united-atom (UA) and all-atom (AA) potentials. Our results using UA potentials predict a value for Δ G ex about four times the experimental value, whereas using AA potentials gives results consistent with experiment, highlighting the need for good potentials to faithfully represent mixture behavior. 
    more » « less