An official website of the United States government
Abstract A general algorithm is introduced to compute single‐chain partition functions in field‐theoretic simulations of polymers with nested tree‐like topologies, including self‐consistent field theory simulations that invoke the mean‐field approximation. The algorithm is an extension of a method used in a number of recent studies on the phase behavior of bottlebrush block copolymers. In those studies, the computational cost of computing single‐chain partition functions is reduced by aggregating the statistical weight of degenerate side arms. By extending this method to chains with arbitrary degrees of branching, the computational cost is reduced to scale with the total length of unique segments in the chain instead of the total length/mass of the entire chain. The method is first validated on a model dendrimer system by comparing results to coarse‐grained molecular dynamics simulations and also demonstrate its advantage over more conventional approaches to compute single‐chain partition functions. The algorithm is subsequently used to analyze the phase behavior of a molecularly informed field‐theoretic model of poly(butyl acrylate)‐graft‐poly(dodecyl acrylate) (pBA‐graft‐pDDA) copolymers in a dodecane solvent. The methodology can help advance field‐theoretic investigations of branched polymers by leveraging degeneracy in the chain to reduce computational cost and avoid the need to develop architecture‐specific algorithms.
more »
« less
This content will become publicly available on June 21, 2026
Improved algorithm for identifying partial saddle-points in polymer field theoretic simulations
Field-theoretic simulations that rely on a partial saddle-point approximation have become powerful tools for studying complex polymer materials. The computational cost of such simulations depends critically upon the efficiency of the iterative algorithm used to identify a partial saddle-point field configuration during each step of a stochastic simulation. We introduce a new algorithm for this purpose that relies on a physically motivated approximation in which the linear response of the density to a small change in a pressure-like field is approximated by the response of a hypothetical homogeneous system. The computational cost of the resulting algorithm is significantly less than that of the commonly used Anderson mixing algorithm.
more »
« less
- Award ID(s):
- 2103627
- PAR ID:
- 10629284
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 162
- Issue:
- 23
- ISSN:
- 0021-9606
- Page Range / eLocation ID:
- 234901
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We present a symmetric formulation of polymer field theory for incompressible systems containing any number M of monomer types, in which all monomers are treated on an equal footing. This is proposed as an alternative to the multispecies exchange formulation, which imposes incompressibility by eliminating one monomer type. The symmetric formulation is shown to correspond to the incompressible limit of a corresponding compressible model, and to reduce in the case M = 2 to the usual formulation of field theory for incompressible AB systems. An analysis of ABC systems (M = 3) identifies ranges of interaction parameter values in which a fully fluctuating field theory requires one, two or three imaginary-valued fields. ABC systems with parameters that satisfy the Hildebrand solubility parameter approximation are shown to require only one imaginary pressure-like field, much like AB systems. Generalization of the partial saddle-point approximation to M > 2 is discussed.more » « less
-
The behavior of complex-Langevin field-theoretic simulations (CL-FTSs) of polymer liquids is sensitive to the nature of saddle-point field configurations, which are solutions of self-consistent field theory (SCFT). Recent work [Kang et al. Macromolecules 2024, 57, 3850] has shown that SCFT saddle-points with real fields are generally not isolated solutions but rather members of a low-dimensional family of continuously-connected complex-valued saddle-points sharing the same Hamiltonian value. We show that this behavior is a natural consequence of the analyticity and translational invariance of the Hamiltonian, which together demand its invariance under generalized translations by displacements with complex components. We also present a numerical algorithm that minimizes the deleterious effects of this generalized symmetry on the stability of CL-FTSs.more » « less
-
We consider optimal experimental design (OED) for Bayesian nonlinear inverse problems governed by partial differential equations (PDEs) under model uncertainty. Specifically, we consider inverse problems in which, in addition to the inversion parameters, the governing PDEs include secondary uncertain parameters. We focus on problems with infinite-dimensional inversion and secondary parameters and present a scalable computational framework for optimal design of such problems. The proposed approach enables Bayesian inversion and OED under uncertainty within a unified framework. We build on the Bayesian approximation error (BAE) approach, to incorporate modeling uncertainties in the Bayesian inverse problem, and methods for A-optimal design of infinite-dimensional Bayesian nonlinear inverse problems. Specifically, a Gaussian approximation to the posterior at the maximuma posterioriprobability point is used to define an uncertainty aware OED objective that is tractable to evaluate and optimize. In particular, the OED objective can be computed at a cost, in the number of PDE solves, that does not grow with the dimension of the discretized inversion and secondary parameters. The OED problem is formulated as a binary bilevel PDE constrained optimization problem and a greedy algorithm, which provides a pragmatic approach, is used to find optimal designs. We demonstrate the effectiveness of the proposed approach for a model inverse problem governed by an elliptic PDE on a three-dimensional domain. Our computational results also highlight the pitfalls of ignoring modeling uncertainties in the OED and/or inference stages.more » « less
-
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
