Concentration-dependent shear response of multi-chain amphiphilic block copolymer self-assemblies

Dominic Robe; Ehsan Kamali Ahangar; Elnaz Hajizadeh

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

Concentration-dependent shear response of multi-chain amphiphilic block copolymer self-assemblies

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

Pith reviewed 2026-05-10 16:34 UTC · model grok-4.3

classification ❄️ cond-mat.soft

keywords amphiphilic block copolymersself-assemblyBrownian dynamicsshear responseviscosity inversionbridging networksdrug deliverymicellar morphology

0 comments

The pith

Diblock copolymers exhibit higher viscosity at equilibrium while triblocks maintain superior viscosity under flow through bridging networks.

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

The paper uses Brownian dynamics simulations to examine how amphiphilic diblock and triblock copolymers self-assemble into micelles and other structures across dilute and semi-dilute concentrations under varying shear. Morphologies shift from spheres to cigars, cylinders, worms, sheets, or fragments depending on hydrophobic fraction and shear rate, with collective reorganization in denser solutions. Rheological results identify a consistent inversion: diblocks are more viscous when at rest, yet triblocks sustain higher viscosity during flow because their architecture allows chains to bridge between aggregates. Aggregation numbers scale differently in each regime, confirming a concentration-driven change, while viscoelastic response follows the same non-terminal power law in all cases due to micellar relaxation times. These patterns offer guidance on controlling flow properties for applications such as drug delivery.

Core claim

Brownian dynamics simulations of multi-chain diblock and triblock copolymers across dilute and semi-dilute regimes establish that equilibrium viscosity is higher for diblocks while triblocks retain higher viscosity under shear via bridging networks. In the dilute regime, quiescent micelles evolve from spherical to cigar-like at moderate shear before fragmenting, with shape depending on hydrophobic fraction f. In semi-dilute conditions, shear drives reorganization toward sheets before fragmentation, again f-dependent. Aggregation number exponents of 0.833 (dilute) and 1.07 (semi-dilute) mark the regime transition, and viscoelasticity shows universal non-terminal power-law scaling governed by,

What carries the argument

Bridging networks formed by triblock copolymer chains that connect separate micelles and sustain viscosity under applied shear.

If this is right

Diblock copolymers would be preferred for formulations that must remain stable and viscous when static.
Triblock copolymers would be preferred when materials must resist flow or maintain structure under shear.
Aggregation number scaling confirms a clear shift from dilute to semi-dilute behavior at the studied concentrations.
Universal non-terminal viscoelastic scaling implies micellar relaxation dynamics control flow response regardless of topology or concentration.
These patterns directly inform how to adjust block copolymer architecture for desired injectability in drug delivery.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Formulations could combine both architectures to achieve rest stability plus shear resistance in a single mixture.
The same inversion may appear in other self-assembling amphiphiles beyond the simulated hydrophobic fractions.
Higher-concentration or entangled regimes might show additional network effects that further enhance triblock performance.
Experimental validation in solvent conditions matching the simulation potentials would strengthen the link to practical use.

Load-bearing premise

The chosen Brownian dynamics parameters, interaction potentials, and shear rates produce self-assembly morphologies and rheological behavior that match those of actual physical block copolymer solutions.

What would settle it

Laboratory rheological measurements on real diblock and triblock copolymer solutions that directly compare zero-shear viscosity against viscosity under comparable shear rates, testing whether the predicted inversion appears.

read the original abstract

Amphiphilic block copolymers self-assemble into diverse nanoscale morphologies with significant implications for drug delivery. This work presents systematic Brownian dynamics simulations of multi-chain diblock and triblock copolymers across dilute and semi-dilute unentangled regimes, hydrophobic fractions, f of 0-1, and shear rates of 0-0.1 1/ns. In the dilute regime, quiescent conditions yield spherical micelles evolving to cigar-like structures at shear rate ~0.01 1/ns and fragmenting at higher shear; varying f produces dispersed chains (f=0), cigar-like (f=0.25), short cylindrical (f=0.5), and gnarled or worm-like (f=0.75) micelles, culminating in sheet-like phase-separated structures (f=1). While, in the semi-dilute regime, shear drives collective reorganisation toward sheet-like morphologies at moderate rates before fragmentation; the f-dependent progression yields cigar-like (f=0.25), sheet-like (f=0.5), and necklace micelles (f=0.75), with larger phase-separated domains at f=1. Rheological characterisation reveals a universal architectural inversion between equilibrium and flow conditions: diblocks show higher equilibrium viscosity while triblocks maintain superior viscosity under flow via bridging networks. Aggregation number scaling exponents of alpha=0.833 in dilute, consistent with star-to-crew-cut bounds of 0.8 to 1.0, and alpha=1.07 in semi-dilute confirm the concentration-driven transition between regimes. Viscoelastic analysis establishes universal non-terminal power-law scaling across all conditions, governed by micellar relaxation dynamics independent of concentration or topology. These findings provide valuable insights into tailoring the injectability and flow behaviour of block copolymers in drug delivery formulations.

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 simulations map a viscosity inversion where diblocks win at rest but triblocks do better under shear via bridging, backed by morphology sweeps and scaling exponents, though the physical mapping remains unanchored.

read the letter

The simulations find that diblock copolymers have higher viscosity at rest, but triblocks take over under shear because of bridging networks that form. This inversion comes out of the Brownian dynamics runs across different concentrations and hydrophobic fractions. They cover a good range: dilute and semi-dilute, f from 0 to 1, shear from 0 to 0.1 per ns. The morphology sequences are laid out clearly, with things like spherical micelles turning into cigars then fragments, or in semi-dilute going to sheets then necklaces. The aggregation number scales as 0.833 in dilute, which fits the expected range, and 1.07 in semi-dilute, marking the regime change. The viscoelastic response shows the same power-law behavior no matter the setup. What stands out is the systematic comparison between the two topologies under both equilibrium and flow. That gives practical pointers for tuning injectability in drug delivery systems. The main concern is whether these viscosities and transitions hold up in real materials. The abstract gives no comparison to experiments, no details on how the interaction potentials were chosen or coarse-grained, and no conversion of simulation units to physical ones. At shear rates of 0.1 per nanosecond, you have to wonder if you're still in the linear regime or seeing something the model exaggerates. Without that, the universal inversion might not translate. This work is aimed at computational soft matter folks and formulation scientists who use simulations to guide polymer design. It has enough new data points on the architecture effect to merit a serious look from referees, who could push for validation or more analysis on the shear regime. I would recommend sending it for peer review rather than desk rejecting it.

Referee Report

2 major / 2 minor

Summary. The manuscript uses Brownian dynamics simulations to study multi-chain diblock and triblock amphiphilic block copolymers in dilute and semi-dilute unentangled regimes. It reports f-dependent morphology transitions (spherical to cigar-like to cylindrical to sheet-like) under quiescent and shear conditions (0-0.1 1/ns), aggregation-number scaling exponents of 0.833 (dilute) and 1.07 (semi-dilute), a rheological inversion (diblocks higher equilibrium viscosity; triblocks higher under flow via bridging), and universal non-terminal power-law viscoelastic scaling independent of concentration or topology.

Significance. If the simulated viscosities and morphologies map to physical block-copolymer solutions, the architectural inversion and scaling results would provide concrete guidance for tuning injectability and flow behavior in drug-delivery formulations.

major comments (2)

[Abstract] Abstract (rheological characterisation paragraph): the central claim of a 'universal architectural inversion' (diblocks higher equilibrium viscosity, triblocks superior under flow via bridging) is load-bearing but rests on unspecified interaction potentials, Brownian dynamics parameters, and shear rates up to 0.1 1/ns with no reported unit conversion, Lees-Edwards stress-tensor extraction details, or comparison to experimental viscosities of comparable systems; this leaves open whether the inversion is physical or an artifact of the model regime.
[Abstract] Abstract (aggregation number scaling sentence): the exponents α=0.833 (dilute, consistent with 0.8-1.0 bounds) and α=1.07 (semi-dilute) are stated without error bars, the concentration range fitted, the precise definition of aggregation number, or the fitting procedure, which is required to substantiate the claimed concentration-driven regime transition.

minor comments (2)

[Abstract] The abstract does not specify the number of chains or system size used in the multi-chain simulations, which is needed to evaluate finite-size effects on self-assembly and bridging statistics.
[Abstract] Shear rates are given only in absolute units (1/ns) without Weissenberg numbers or discussion of the linear-response regime for the micellar relaxation times.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful and constructive review of our manuscript. We address each major comment point by point below, providing clarifications on methodology and results. Revisions have been made to strengthen the presentation of key claims.

read point-by-point responses

Referee: [Abstract] Abstract (rheological characterisation paragraph): the central claim of a 'universal architectural inversion' (diblocks higher equilibrium viscosity, triblocks superior under flow via bridging) is load-bearing but rests on unspecified interaction potentials, Brownian dynamics parameters, and shear rates up to 0.1 1/ns with no reported unit conversion, Lees-Edwards stress-tensor extraction details, or comparison to experimental viscosities of comparable systems; this leaves open whether the inversion is physical or an artifact of the model regime.

Authors: We acknowledge that the abstract's brevity omits full methodological specifics, which are detailed in the Methods section: interaction potentials follow standard Lennard-Jones for hydrophobic/hydrophilic beads and FENE for chain connectivity; Brownian dynamics uses a friction coefficient of 1.0 and time step of 0.001 in reduced units. Shear rates are reported in simulation units (1/ns) as conventional; we have added a clarifying sentence on unit mapping to the revised abstract and discussion. Lees-Edwards boundary conditions are applied for simple shear, with the stress tensor obtained from the virial expression (full protocol in SI). Direct experimental viscosity comparisons are not performed, as the study emphasizes mechanistic trends, but we have expanded the discussion to note consistency with bridging effects observed in experimental triblock systems. The inversion arises robustly from triblock bridging networks under flow versus diblock behavior at equilibrium, supported by morphology and scaling data across regimes, indicating it is physical within the model rather than an artifact. revision: partial
Referee: [Abstract] Abstract (aggregation number scaling sentence): the exponents α=0.833 (dilute, consistent with 0.8-1.0 bounds) and α=1.07 (semi-dilute) are stated without error bars, the concentration range fitted, the precise definition of aggregation number, or the fitting procedure, which is required to substantiate the claimed concentration-driven regime transition.

Authors: We agree this information strengthens the claim and have revised the manuscript accordingly. Aggregation number is defined as the average number of copolymer chains per micellar aggregate, computed from cluster analysis with a cutoff distance of 1.5σ. The dilute fit uses volume fractions 0.005–0.05 and the semi-dilute fit uses 0.1–0.3; linear regression was performed on log-log plots of aggregation number versus concentration. Error bars from the fits are now reported: α = 0.833 ± 0.015 (dilute) and α = 1.07 ± 0.04 (semi-dilute). These details have been added to the Results section and referenced in the abstract to substantiate the dilute-to-semi-dilute transition. revision: yes

Circularity Check

0 steps flagged

No circularity: results are direct simulation outputs with no analytical derivation or self-referential reduction

full rationale

The paper consists entirely of Brownian dynamics simulation results across specified regimes, shear rates, and hydrophobic fractions. No analytical derivation, equations, or predictive model is presented whose outputs reduce to fitted inputs by construction. Claims such as the 'universal architectural inversion' and aggregation number scaling are reported as direct observations from the simulations (e.g., viscosity extracted from stress tensor under Lees-Edwards shear). No self-citations, uniqueness theorems, or ansatzes are invoked in a load-bearing way. The absence of an analytical chain means no opportunity for self-definitional or fitted-input circularity exists. This is the expected non-finding for a purely computational study.

Axiom & Free-Parameter Ledger

3 free parameters · 1 axioms · 0 invented entities

Central claims rest on the fidelity of the coarse-grained Brownian dynamics model to real amphiphilic copolymers and on the assumption that the chosen interaction strengths and chain lengths produce representative self-assembly.

free parameters (3)

hydrophobic fraction f
Varied parametrically from 0 to 1 as an input to control morphology.
shear rate
Scanned from 0 to 0.1 1/ns as an external driving parameter.
interaction potentials
Model-specific bead-bead attractions and repulsions that determine micelle stability; values not specified in abstract.

axioms (1)

domain assumption Brownian dynamics with the selected time step and thermostat accurately captures the overdamped dynamics of polymer chains in solution.
Invoked for all reported trajectories and rheological measurements.

pith-pipeline@v0.9.0 · 5639 in / 1326 out tokens · 29364 ms · 2026-05-10T16:34:40.474766+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages · 1 internal anchor

[1]

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

S. Liu and R. Sureshkumar, “Deformation, Rupture, and Morphology Hysteresis of Copolymer Nanovesicles in Uniform Shear Flow,” Langmuir, vol. 41, no. 8, pp. 5083–5096, 2025, doi: 10.1021/acs.langmuir.4c04200. [24] R. H. Colby, D. C. Boris, W. E. Krause, and S. Dou, “Shear thinning of unentangled flexible polymer liquids,” Rheol. Acta, vol. 46, no. 5, pp. 5...

work page internal anchor Pith review Pith/arXiv arXiv doi:10.1021/acs.langmuir.4c04200 2025
[2]

Constrained Bayesian accelerated design of acoustic polyurethane coatings with metamaterial features under hydrostatic pressure

P. Ding, D. Robe, M. Kirley, E. Hajizadeh, " Constrained Bayesian accelerated design of acoustic polyurethane coatings with metamaterial features under hydrostatic pressure", Structural and Multidisciplinary Optimization, vol. 68, no. 8, pp. 159, 2025, https://doi.org/10.1007/s00158-025-04107-7 [41] H. Weeratunge, Z. Shireen, S. Iyer, A. Menzel, A. W. Phi...

work page doi:10.1007/s00158-025-04107-7 2025

[1] [1]

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

S. Liu and R. Sureshkumar, “Deformation, Rupture, and Morphology Hysteresis of Copolymer Nanovesicles in Uniform Shear Flow,” Langmuir, vol. 41, no. 8, pp. 5083–5096, 2025, doi: 10.1021/acs.langmuir.4c04200. [24] R. H. Colby, D. C. Boris, W. E. Krause, and S. Dou, “Shear thinning of unentangled flexible polymer liquids,” Rheol. Acta, vol. 46, no. 5, pp. 5...

work page internal anchor Pith review Pith/arXiv arXiv doi:10.1021/acs.langmuir.4c04200 2025

[2] [2]

Constrained Bayesian accelerated design of acoustic polyurethane coatings with metamaterial features under hydrostatic pressure

P. Ding, D. Robe, M. Kirley, E. Hajizadeh, " Constrained Bayesian accelerated design of acoustic polyurethane coatings with metamaterial features under hydrostatic pressure", Structural and Multidisciplinary Optimization, vol. 68, no. 8, pp. 159, 2025, https://doi.org/10.1007/s00158-025-04107-7 [41] H. Weeratunge, Z. Shireen, S. Iyer, A. Menzel, A. W. Phi...

work page doi:10.1007/s00158-025-04107-7 2025