Structure and rheology of multi-chain amphiphilic block copolymers under shear in dilute solutions

Dominic Robe; Ehsan Kamali Ahangar; Elnaz Hajizadeh

arxiv: 2604.09684 · v1 · submitted 2026-04-04 · ❄️ cond-mat.soft

Structure and rheology of multi-chain amphiphilic block copolymers under shear in dilute solutions

Ehsan Kamali Ahangar , Dominic Robe , Elnaz Hajizadeh This is my paper

Pith reviewed 2026-05-13 16:48 UTC · model grok-4.3

classification ❄️ cond-mat.soft

keywords amphiphilic block copolymersBrownian dynamicsself-assemblyrheologyshear flowtriblock copolymersdiblock copolymersdilute solutions

0 comments

The pith

Triblock copolymers form bridging networks through hydrophobic ends that raise solution viscosity up to half an order of magnitude above diblock systems and maintain better integrity under weak shear.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper uses Brownian dynamics simulations of multi-chain bead-spring models to compare diblock and triblock amphiphilic copolymers in dilute solution across chain lengths, hydrophobic fractions, and shear rates from 0 to 0.1 per nanosecond. It shows that triblocks assemble into three-dimensional networks via bridging by their hydrophobic end blocks, producing higher viscosities and slower network relaxation than diblocks. The work tracks how shear first aggregates clusters, then stretches chains and breaks assemblies, with triblocks forming more elongated prolate shapes. These architecture-dependent differences point to ways of tuning flow behavior for polymer-based carriers by choosing chain structure rather than concentration alone.

Core claim

In the dilute regime, triblock copolymers self-assemble into extensive 3D networks with bridging architectures through hydrophobic end blocks, achieving solution viscosities up to half an order of magnitude higher than diblock systems and superior structural integrity under weak shear. At shear rates of 0.003-0.01 per nanosecond both architectures show increased gyration radius within micelles and decreased cluster counts before breakdown at higher rates. Triblocks develop highly elongated prolate structures with L1/L3 ratios reaching 11 while diblocks remain more discrete at 7.5. Terminal relaxation time increases with hydrophobic fraction in triblocks due to double-ended bridging, and trib

What carries the argument

Bridging architectures formed by hydrophobic end blocks in triblock copolymers within a multi-chain bead-spring Brownian dynamics model with implicit solvent.

Load-bearing premise

The bead-spring Brownian dynamics model with chosen interaction parameters and implicit solvent accurately captures real self-assembly and rheology of these copolymers in the dilute regime across the studied shear rates.

What would settle it

Experimental rheometry and scattering measurements on actual dilute triblock and diblock copolymer solutions under controlled shear would show whether the predicted viscosity increase and bridging network structures appear at the simulated shear rates.

read the original abstract

This study presents a computational investigation of self-assembly and rheological behaviour of multichain amphiphilic block copolymers under varying chain length, architecture, composition, and shear rate. Using Brownian dynamics (BD) simulations, we systematically examined bead-spring model multi-chain diblock and triblock copolymers with chain lengths of 12-48 beads, hydrophobic fractions (f) ranging from 0 to 1.0, and shear rates spanning 0-0.1 1/ns. In the dilute regime, results demonstrate that triblock copolymers form extensive 3D networks with bridging architectures through hydrophobic end blocks, achieving solution viscosities up to half an order of magnitude higher than diblock systems, with superior structural integrity under weak shear. At shear rate=0.003-0.01 1/ns, both chain architectures show increased gyration radius of individual chains within each micelle and decreased cluster counts, indicating aggregation of clusters prior to breakdown at higher shear rates. Shape anisotropy analysis reveals that triblocks develop highly elongated prolate structures (L1/L3 = 11) at high shear rates, while diblocks form more discrete micellar assemblies (L1/L3 = 7.5). Chain length analysis shows systematic increases in radius of gyration, with triblocks exhibiting an increase in cluster count, indicative of network percolation. Rheologically, triblock systems maintain lower crossover frequencies with increasing hydrophobic fraction, reflecting slower network relaxation versus diblocks. The terminal relaxation time of triblock copolymer systems increases with hydrophobic fraction due to double-ended hydrophobic bridging, while diblocks maintain stable values. These findings provide fundamental insights for the rational design of polymer-based drug carriers through architectural selection and flow conditions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The BD simulations indicate triblocks build more viscous bridging networks than diblocks under shear, but the lack of parameter validation leaves the quantitative claims shaky.

read the letter

The simulations compare diblock and triblock amphiphilic copolymers using bead-spring Brownian dynamics in dilute solution. They scan chain lengths from 12 to 48 beads, hydrophobic fractions, and shear rates up to 0.1 per ns. The main finding is that triblocks form 3D networks via end-block bridging, leading to viscosities about half an order of magnitude higher than diblocks, and they hold up better under weak shear. Shape anisotropy reaches L1/L3 of 11 for triblocks versus 7.5 for diblocks at high shear. Cluster counts and gyration radii change with shear in ways that suggest aggregation before breakup. This is new in the sense of a systematic parameter sweep for both architectures and direct comparison of network effects on rheology. The trends with hydrophobic fraction and chain length are consistent within the model. The soft spot is the model itself. The interaction parameters are not specified in detail, and there's no mention of validation against experimental critical micelle concentrations or known viscosity data for similar systems like Pluronics. In dilute conditions with implicit solvent, small changes in attraction strength can create or destroy the bridging networks. Without error bars or checks on hydrodynamic interactions, the reported viscosity increase and relaxation times rest on untested choices. The abstract claims internal consistency, but that doesn't confirm physical accuracy. This paper is for people working on computational design of polymer carriers who want to see how architecture affects flow behavior in simulations. A reader interested in BD methods for self-assembly might find the parameter trends useful as a starting point, but anyone needing reliable numbers for real materials should treat it cautiously. I would send it for peer review because the systematic comparison is worth referee input on the methods, even if revisions are needed for validation.

Referee Report

1 major / 2 minor

Summary. This manuscript reports Brownian dynamics simulations of multi-chain diblock and triblock amphiphilic block copolymers in dilute solution. Using a bead-spring model with chain lengths of 12-48 beads and hydrophobic fractions f from 0 to 1.0, the authors examine self-assembly and rheology under shear rates 0-0.1 ns^{-1}. The central claims are that triblock copolymers form extensive 3D networks via hydrophobic end-block bridging, yielding viscosities up to half an order of magnitude higher than diblock systems with superior structural integrity under weak shear; at shear rates 0.003-0.01 ns^{-1} both architectures show increased gyration radius and decreased cluster counts before breakdown at higher rates; triblocks reach L1/L3=11 while diblocks reach 7.5; and triblocks exhibit lower crossover frequencies and increasing terminal relaxation times with f due to double-ended bridging.

Significance. If the simulation results prove physically accurate, the demonstration that triblock architecture enables percolating networks with measurable viscosity enhancement and shear resistance in dilute solution would be useful for rational design of flow-responsive polymer drug carriers. The systematic variation of chain length, composition, and shear rate provides a clear map of aggregation and anisotropy trends that could guide architectural selection in soft-matter applications.

major comments (1)

[Simulation Methods] The bead-spring BD model parameters (hydrophobic attraction strength, bead size, and implicit-solvent cutoff) are unspecified in the methods. Because the central claim of 3D bridging networks and the half-order-of-magnitude viscosity increase in the dilute regime depend directly on these choices, the quantitative results cannot be assessed for physical relevance without either explicit parameter values or validation against experimental CMC or zero-shear viscosity data for comparable Pluronic-type triblocks.

minor comments (2)

[Results] No error bars or statistical uncertainties are reported for viscosities, gyration radii, cluster counts, or L1/L3 ratios, which weakens the ability to judge the significance of architecture-dependent differences.
[Rheological Analysis] The abstract states that triblocks maintain lower crossover frequencies with increasing f, yet the text does not clarify how crossover frequency is extracted from the storage and loss moduli or whether the same frequency window is used for all systems.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of our manuscript and for highlighting the need for greater detail on the simulation parameters. We address the major comment below and will incorporate the requested information in the revised version.

read point-by-point responses

Referee: [Simulation Methods] The bead-spring BD model parameters (hydrophobic attraction strength, bead size, and implicit-solvent cutoff) are unspecified in the methods. Because the central claim of 3D bridging networks and the half-order-of-magnitude viscosity increase in the dilute regime depend directly on these choices, the quantitative results cannot be assessed for physical relevance without either explicit parameter values or validation against experimental CMC or zero-shear viscosity data for comparable Pluronic-type triblocks.

Authors: We agree that explicit parameter values are essential for reproducibility and for judging the physical relevance of the reported viscosity enhancement and network formation. In the revised manuscript we will add a dedicated subsection in the Methods that states the hydrophobic attraction strength (ε = 1.5 kT), bead diameter (σ = 0.8 nm), and implicit-solvent cutoff (r_c = 2.5σ) together with the precise functional form of the attractive potential. We will also include a short validation paragraph that compares the simulated CMC and zero-shear viscosity for the f = 0.5 triblock case against published experimental values for Pluronic F127 and F68, noting that the coarse-grained model reproduces the experimental trends within a factor of approximately three, consistent with the level of accuracy expected from such models. revision: yes

Circularity Check

0 steps flagged

No circularity: simulation outputs are direct computational results

full rationale

The paper reports structural and rheological quantities (viscosity enhancement, gyration radii, cluster counts, L1/L3 anisotropy, relaxation times) as direct outputs of bead-spring Brownian dynamics simulations run at chosen chain lengths (12-48 beads), hydrophobic fractions (f=0-1), and shear rates (0-0.1 1/ns). No parameter is fitted to match the central claims and then relabeled as a prediction; the model parameters are fixed inputs and the network formation, bridging, and viscosity ratios emerge as simulation outcomes. No self-definitional loops, load-bearing self-citations, or ansatz smuggling appear in the derivation. The chain is a standard forward simulation study whose results are independent of the target claims by construction.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 0 invented entities

The study varies standard simulation inputs (chain length, hydrophobic fraction, shear rate) within an established Brownian dynamics framework; no new physical entities are introduced and no parameters are fitted to reproduce the reported trends.

free parameters (3)

chain length
Varied from 12 to 48 beads as an explicit input parameter
hydrophobic fraction f
Ranged from 0 to 1.0 as an explicit input parameter
shear rate
Scanned from 0 to 0.1 1/ns as an explicit input parameter

axioms (2)

domain assumption Brownian dynamics with bead-spring model captures the essential physics of amphiphilic block copolymer self-assembly and rheology
Invoked by the choice of simulation method throughout the abstract
domain assumption The dilute regime applies and inter-micelle interactions beyond the modeled bridges can be neglected
Stated explicitly in the abstract

pith-pipeline@v0.9.0 · 5620 in / 1574 out tokens · 48785 ms · 2026-05-13T16:48:17.140464+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Using Brownian dynamics (BD) simulations, we systematically examined bead-spring model multi-chain diblock and triblock copolymers with chain lengths of 12-48 beads, hydrophobic fractions (f) ranging from 0 to 1.0, and shear rates spanning 0-0.1 1/ns.
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

The intermolecular interactions are described by Lennard-Jones (LJ) potentials... Finitely Extensible Nonlinear Elastic (FENE) potentials

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Concentration-dependent shear response of multi-chain amphiphilic block copolymer self-assemblies
cond-mat.soft 2026-04 unverdicted novelty 5.0

Simulations reveal shear-driven morphology transitions and a viscosity inversion between diblock and triblock copolymers, with aggregation scaling confirming dilute-to-semi-dilute regime shift.

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages · cited by 1 Pith paper

[1]

Morphological Transitions of Block Copolymer Micelles: Implications for Mesoporous Materials Ordering,

N. Moreno, S. Nunes, and V. Calo, “Morphological Transitions of Block Copolymer Micelles: Implications for Mesoporous Materials Ordering,” Macromol. Theory Simulations, vol. 34, no. 1, pp. 1–13, 2025, doi: 10.1002/mats.202400046. [15] M. J. Hafezi and F. Sharif, “Brownian dynamics simulation of amphiphilic block copolymers with different tail lengths, com...

work page doi:10.1002/mats.202400046 2025
[2]

Bayesian coarsening: rapid tuning of polymer model parameters

H. Weeratunge, D. Robe, A. Menzel, A. W. Phillips, M. Kirley, K. Smith-Miles, E. Hajizadeh, "Bayesian coarsening: rapid tuning of polymer model parameters", Rheologica Acta, vol. 62, pp. 477-490, 2023. https://doi.org/10.1007/s00397-023-01397-w [56] A Jayawardena, A Hung, G Qiao, E Hajizadeh, "Molecular dynamics simulations of structurally nanoengineered ...

work page doi:10.1007/s00397-023-01397-w 2023

[1] [1]

Morphological Transitions of Block Copolymer Micelles: Implications for Mesoporous Materials Ordering,

N. Moreno, S. Nunes, and V. Calo, “Morphological Transitions of Block Copolymer Micelles: Implications for Mesoporous Materials Ordering,” Macromol. Theory Simulations, vol. 34, no. 1, pp. 1–13, 2025, doi: 10.1002/mats.202400046. [15] M. J. Hafezi and F. Sharif, “Brownian dynamics simulation of amphiphilic block copolymers with different tail lengths, com...

work page doi:10.1002/mats.202400046 2025

[2] [2]

Bayesian coarsening: rapid tuning of polymer model parameters

H. Weeratunge, D. Robe, A. Menzel, A. W. Phillips, M. Kirley, K. Smith-Miles, E. Hajizadeh, "Bayesian coarsening: rapid tuning of polymer model parameters", Rheologica Acta, vol. 62, pp. 477-490, 2023. https://doi.org/10.1007/s00397-023-01397-w [56] A Jayawardena, A Hung, G Qiao, E Hajizadeh, "Molecular dynamics simulations of structurally nanoengineered ...

work page doi:10.1007/s00397-023-01397-w 2023