Mobility Heterogeneity in a 2D Gaussian Lattice Polymer: A Dynamic Monte Carlo Study

This study demonstrates through dynamic Monte Carlo simulations that while introducing mobility heterogeneity via different update rates in a two-block 2D Gaussian lattice polymer alters internal relaxation dynamics and block-resolved mean-squared displacement, the center-of-mass diffusion coefficient retains the standard ideal chain scaling of $D_{\rm cm} \sim N^{-1}$.

The Gist

Imagine a long, flexible snake made of beads, slithering across a grid floor. This is a polymer chain, a molecule found in everything from plastics to DNA. In this study, the author, Arpan Dey, uses a computer simulation to watch how this snake moves.

Here is the story of what he found, explained simply:

1. The Rules of the Game (The "Dictionary")

First, the author needed a set of rules for how the snake could move. He created a "move dictionary."

• The Grid: The snake lives on a square grid (like graph paper).
• The Constraint: The beads are connected by strings of fixed length. A bead can only move if it stays connected to its neighbors.
• The Moves:
  • Ends: The head and tail beads can wiggle to any empty spot next to them.
  • Middle: A bead in the middle is stuck between two neighbors. It can only move if it's at a "corner" of the grid, allowing it to flip to the opposite corner without breaking the strings.
• The Benchmark: When every bead gets an equal chance to try moving, the snake behaves exactly as physics predicts for a "perfect" chain (called the Rouse model). It wiggles locally, but the whole snake drifts slowly, and longer snakes drift even slower.

2. The Experiment: The "Lazy" vs. "Energetic" Snake

Next, the author wanted to see what happens if the snake isn't uniform. He split the snake into two halves:

• Block A (The Energetic Half): These beads get to try moving more often.
• Block B (The Lazy Half): These beads get to try moving less often.

Think of it like a relay race where the first half of the team is told to run as fast as they can, while the second half is told to jog slowly. The rules for how they move (the dictionary) stay the same; only the frequency of their attempts changes.

3. What Happened?

The results were a mix of "obvious" and "surprising."

The Obvious Part (Internal Chaos):
As expected, the "Energetic Half" (Block A) wiggled much more than the "Lazy Half" (Block B). If you measured how far each half traveled, the energetic side was clearly more active. The snake became asymmetric; one side was doing all the work while the other dragged its feet.

The Surprising Part (The Whole Snake):
Here is the big twist. Even though one half was frantic and the other was lazy, the speed of the entire snake's center did not change its fundamental rule.

In physics, there is a rule that says: The longer the snake, the slower it moves as a whole. Specifically, if you double the length of the snake, it moves half as fast.

• The Finding: Even with the "Energetic" and "Lazy" halves, the entire snake still followed this exact rule. Whether the snake was short or long, and whether the halves were equally active or very different, the overall speed still dropped in perfect proportion to the length.

4. Why Did This Happen? (The Analogy)

The author explains this with a simple logic:
Imagine the snake is a team of people pulling a heavy cart.

• If everyone pulls at the same speed, the cart moves at a certain pace.
• If half the team pulls twice as hard and the other half pulls half as hard, the total effort of the team changes slightly, but the relationship between the team size and the speed stays the same.

The "friction" (resistance to moving) of the whole snake is just the sum of the friction of all its parts. Because the snake is still one connected object, the internal differences (one side fast, one side slow) cancel out in a way that preserves the overall scaling law. The "Energetic" half doesn't drag the "Lazy" half fast enough to break the rule that "longer chains move slower."

5. The Bottom Line

The paper concludes that mobility heterogeneity (having parts of a molecule that are more active than others) changes how the molecule wiggles internally, but it does not change the fundamental law of how fast the whole molecule drifts through space.

• Internal motion: Changes drastically (one side moves more).
• Overall drift: Stays on the same predictable path ($D \sim 1/N$).

The author notes that this was tested on a "Gaussian" (ideal, non-sticky) snake. He tried to test it on a "sticky" snake (where beads can't overlap), but the simulation got too stuck to give clear results. So, this specific finding applies to the ideal, non-sticky version of the model.

In short: You can make one half of a polymer chain frantic and the other half lazy, and while the snake will look very uneven on the inside, its overall journey across the floor will still follow the same old, predictable rules.

Technical Summary

Technical Summary: Mobility Heterogeneity in a 2D Gaussian Lattice Polymer

Problem Statement
Polymer dynamics are governed by the interplay between stochastic monomer motion and connectivity constraints. While the Rouse model provides a standard description for ideal chains, real-world systems often exhibit mobility heterogeneity due to varying local environments, temperatures, or active forcing. Existing studies have explored heterogeneous Rouse models and active chains, but there is a need to isolate the specific effects of rate-induced mobility heterogeneity—where different segments of a chain are updated at different frequencies—without altering the underlying local move geometry or introducing energetic interactions. This paper addresses how such rate heterogeneity affects internal relaxation versus global center-of-mass transport in a minimal Gaussian setting.

Methodology
The study employs Dynamic Monte Carlo (DMC) simulations on a two-dimensional square lattice.

• Model: The polymer is modeled as a Gaussian chain (monomer overlap allowed) of length $N$. Connectivity is enforced via discrete, bond-preserving moves rather than harmonic springs.
• Move Dictionary: A computationally efficient three-monomer move dictionary is utilized. This pre-tabulated dictionary enumerates all allowed local configurations (bent, straight, and overlapping generators) for internal monomers based on their neighbors' positions, while end monomers follow standard end-bond rotation rules.
• Homogeneous Benchmark: Initially, a homogeneous chain is simulated where all monomers have equal selection probabilities. This validates the lattice algorithm against continuum Rouse predictions, specifically the mean-squared displacement (MSD) crossover and the scaling of the center-of-mass diffusion coefficient ($D_{cm} \sim N^{-1}$).
• Heterogeneous Model: The chain is divided into two equal blocks (Block A and Block B). While the local move dictionary remains identical for both, the attempt rates differ. Block A is selected with rate $\omega_A$ and Block B with $\omega_B$. The rate ratio is defined as $\rho = \omega_A / \omega_B$.
• Metrics: The study analyzes block-resolved MSDs, a normalized MSD asymmetry ($A_{MSD}$), and the chain-length dependence of $D_{cm}$ for various $\rho$ values (1, 2, and 4) and chain lengths ($N=20$ to $100$).

Key Results

1. Benchmark Validation: The homogeneous lattice simulation successfully reproduces Rouse-like behavior. The monomer MSD exhibits the expected crossover from subdiffusive motion (slope $\sim 1/2$) to diffusive motion (slope $\sim 1$) at times scaled by $N^2$. The center-of-mass diffusion coefficient scales as $D_{cm} \sim N^{-1}$.
2. Internal Dynamics: In the heterogeneous model ($\rho > 1$), the more frequently updated block (Block A) exhibits a larger MSD at early and intermediate times compared to Block B. This results in a positive normalized MSD asymmetry ($A_{MSD} > 0$). However, the qualitative Rouse-like crossover structure is preserved for both blocks.
3. Global Transport: Crucially, despite the internal mobility imbalance, the center-of-mass diffusion coefficient retains the Rouse scaling $D_{cm} \sim N^{-1}$ for all studied rate ratios. While the rate contrast $\rho$ alters the prefactor of $D_{cm}$ (specifically reducing it as $\rho$ increases), it does not change the scaling exponent with respect to chain length $N$.

Analytical Explanation
The paper provides a coarse-grained Rouse argument to explain the preservation of the $N^{-1}$ scaling. In the continuum limit, the attempt rate $q_i$ is inversely proportional to the effective friction $\gamma_i$ of a monomer. Summing the equations of motion for the entire chain, the internal spring forces cancel out, leaving a total friction $\Gamma_{tot} = \sum \gamma_i$. Since $D_{cm} \propto 1/\Gamma_{tot}$, and the total friction is a linear sum of the frictions of the two blocks (each proportional to $N/2$), the resulting diffusion coefficient scales as $1/N$. The rate ratio $\rho$ modifies the numerical prefactor (increasing total effective friction relative to the homogeneous case) but leaves the $N^{-1}$ dependence intact. This logic extends to polymers with more than two blocks, provided the fraction of monomers in each block remains fixed.

Significance and Claims
The paper claims that in a minimal Gaussian setting, rate-induced mobility heterogeneity modifies internal relaxation timescales without altering the fundamental Rouse scaling of center-of-mass transport.

• The study isolates the effect of update frequency from other sources of heterogeneity (such as excluded volume or force fields), demonstrating that the $D_{cm} \sim N^{-1}$ scaling is robust against variations in local attempt rates.
• The authors note that extending this specific construction to self-avoiding chains is non-trivial; preliminary tests showed slow relaxation and potential arrest of the junction monomer, preventing a clean extraction of scaling laws in that regime.
• The work serves as a controlled test case for understanding how local dynamic heterogeneity propagates (or fails to propagate) to global transport properties, offering a baseline for more complex models involving active polymers or heterogeneous environments.