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.

Attention:

The DOI auto-population feature in the Public Access Repository (PAR) will be unavailable from 4:00 PM ET on Tuesday, July 8 until 4:00 PM ET on Wednesday, July 9 due to scheduled maintenance. We apologize for the inconvenience caused.


Title: Assessment of the partial saddle point approximation in field-theoretic polymer simulations
Field-theoretic simulations are numerical treatments of polymer field theory models that go beyond the mean-field self-consistent field theory level and have successfully captured a range of mesoscopic phenomena. Inherent in molecularly-based field theories is a “sign problem” associated with complex-valued Hamiltonian functionals. One route to field-theoretic simulations utilizes the complex Langevin (CL) method to importance sample complex-valued field configurations to bypass the sign problem. Although CL is exact in principle, it can be difficult to stabilize in strongly fluctuating systems. An alternate approach for blends or block copolymers with two segment species is to make a “partial saddle point approximation” (PSPA) in which the stiff pressure-like field is constrained to its mean-field value, eliminating the sign problem in the remaining field theory, allowing for traditional (real) sampling methods. The consequences of the PSPA are relatively unknown, and direct comparisons between the two methods are limited. Here, we quantitatively compare thermodynamic observables, order-disorder transitions, and periodic domain sizes predicted by the two approaches for a weakly compressible model of AB diblock copolymers. Using Gaussian fluctuation analysis, we validate our simulation observations, finding that the PSPA incorrectly captures trends in fluctuation corrections to certain thermodynamic observables, microdomain spacing, and location of order-disorder transitions. For incompressible models with contact interactions, we find similar discrepancies between the predictions of CL and PSPA, but these can be minimized by regularization procedures such as Morse calibration. These findings mandate caution in applying the PSPA to broader classes of soft-matter models and systems.  more » « less
Award ID(s):
2104255
PAR ID:
10516157
Author(s) / Creator(s):
; ;
Publisher / Repository:
American Institute of Physics
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
159
Issue:
16
ISSN:
0021-9606
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Chiral building blocks have the ability to self-assemble and transfer chirality to larger hierarchical length scales, which can be leveraged for the development of novel nanomaterials. Chiral block copolymers, where one block is made completely chiral, are prime candidates for studying this phenomenon, but fundamental questions regarding the self-assembly are still unanswered. For one, experimental studies using different chemistries have shown unexplained diverging shifts in the order–disorder transition temperature. In this study, particle-based molecular simulations of chiral block copolymers in the disordered melt were performed to uncover the thermodynamic behavior of these systems. A wide range of helical models were selected, and several free energy calculations were performed. Specifically, we aimed to understand (1) the thermodynamic impact of changing the conformation of one block in chemically identical block copolymers and (2) the effect of the conformation on the Flory–Huggins interaction parameter, χ, when chemical disparity was introduced. We found that the effective block repulsion exhibits diverging behavior, depending on the specific conformational details of the helical block. Commonly used conformational metrics for flexible or stiff block copolymers do not capture the effective block repulsion because helical blocks are semiflexible and aspherical. Instead, pitch can quantitatively capture the effective block repulsion. Quite remarkably, the shift in χ for chemically dissimilar block copolymers can switch sign with small changes in the pitch of the helix. 
    more » « less
  2. The self-assembly of block polymers into well-ordered nanostructures underpins their utility across fundamental and applied polymer science, yet only a handful of equilibrium morphologies are known with the simplest AB-type materials. Here, we report the discovery of the A15 sphere phase in single-component diblock copolymer melts comprising poly(dodecyl acrylate)− block −poly(lactide). A systematic exploration of phase space revealed that A15 forms across a substantial range of minority lactide block volume fractions ( f L = 0.25 − 0.33) situated between the σ-sphere phase and hexagonally close-packed cylinders. Self-consistent field theory rationalizes the thermodynamic stability of A15 as a consequence of extreme conformational asymmetry. The experimentally observed A15−disorder phase transition is not captured using mean-field approximations but instead arises due to composition fluctuations as evidenced by fully fluctuating field-theoretic simulations. This combination of experiments and field-theoretic simulations provides rational design rules that can be used to generate unique, polymer-based mesophases through self-assembly. 
    more » « less
  3. Flocking phase transitions arise in many aligning active soft matter systems, and an interesting question concerns the role of “topological” vs. “metric” interactions on these transitions. While recent theoretical work suggests that the order–disorder transition in these polar aligning models is universally first order, numerical studies have suggested that topological models may instead have a continuous transition. Some recent simulations have found that some variations of topologically interacting flocking agents have a discontinuous transition, but unambiguous observations of phase coexistence using common Voronoi-based alignment remains elusive. In this work, we use a custom GPU-accelerated simulation package to perform million-particle-scale simulations of a Voronoi–Vicsek model in which alignment interactions stem from an XY-like Hamiltonian. By accessing such large systems on appropriately long time scales and in the time-continuous limit, we are able to show a regime of stable phase coexistence between the ordered and disordered phases, confirming the discontinuous nature of this transition in the thermodynamic limit. 
    more » « less
  4. Abstract We consider simple mean field continuum models for first order liquid–liquid demixing and solid–liquid phase transitions and show how the Maxwell construction at phase coexistence emerges on going from finite-size closed systems to the thermodynamic limit. The theories considered are the Cahn–Hilliard model of phase separation, which is also a model for the liquid-gas transition, and the phase field crystal model of the solid–liquid transition. Our results show that states comprising the Maxwell line depend strongly on the mean density with spatially localized structures playing a key role in the approach to the thermodynamic limit. 
    more » « less
  5. Experimental quantum simulators have become large and complex enough that discovering new physics from the huge amount of measurement data can be quite challenging, especially when little theoretical understanding of the simulated model is available. Unsupervised machine learning methods are particularly promising in overcoming this challenge. For the specific task of learning quantum phase transitions, unsupervised machine learning methods have primarily been developed for phase transitions characterized by simple order parameters, typically linear in the measured observables. However, such methods often fail for more complicated phase transitions, such as those involving incommensurate phases, valence-bond solids, topological order, and many-body localization. We show that the diffusion map method, which performs nonlinear dimensionality reduction and spectral clustering of the measurement data, has significant potential for learning such complex phase transitions unsupervised. This method works for measurements of local observables in a single basis and is thus readily applicable to many experimental quantum simulators as a versatile tool for learning various quantum phases and phase transitions. 
    more » « less