Explore the Glass Transition with Coarse-Grained Molecular Simulations

← Back to AIMS

By Charles Xie ✉

A plastic cup is made of polystyrene, long chains of styrene units tangled together. A single chain has thousands of atoms, and a piece of plastic holds many chains, far more than a web browser can simulate one atom at a time. AIMS models each styrene unit as a single bead instead, and an AI agent sets up the experiments. A student asks a plain-language question, and the agent cools the melt step by step, measures how far its beads move, and posts the measured data in the chat as a live chart. In the window below, twenty chains of twenty beads fill a box at 450 K, well above the glass transition. Press Play and watch the chains slide past one another.

Why Trade Atoms for Beads?

An atomistic simulation computes the force on every atom from every nearby atom, a femtosecond or two at a time. For a few hundred atoms a browser can do that fast enough to watch, and the phase-change and heat-conduction articles on this site work that way. What makes a glass is chains too tangled to slide past one another, which involves hundreds of chain segments and takes nanoseconds. A model that keeps every hydrogen atom spends its whole budget on detail the glass transition does not depend on.

Coarse-grained modeling keeps the degrees of freedom that matter and drops the rest. A landmark of the approach is the 1990 simulation by Kremer and Grest, who modeled polymer melts as chains of soft beads on springs and reproduced how the chains move and tangle. A bead has a mass, a size, and an attraction to other beads, and chains of beads have springs, bending stiffness, and a torsion term. One bead in this lab has the mass of a styrene unit, 104 g/mol, and an attraction fitted so that the melt stops flowing near the glass transition of polystyrene. Because the units are real, the thermometer in the lab reads kelvin. Companion articles apply the same approach to liquid crystals, protein folding, and soap micelles.

Running Experiments by Asking

AIMS is an AI-native environment. Building a melt, stepping its temperature, measuring how far its beads move, and plotting the result are the agent's work. A student asks "when does the melt become a glass?", and the agent recognizes the intent, configures the lab, runs the experiment, and posts a live chart of the results in the conversation.

A Melt and Its Glass Transition

A polymer melt is a tangle of long chains. Warm, the chains slide past one another and the material flows; cold, they lock in place and the material is a glass, hard and brittle and without any crystal order. The temperature at which that happens is the glass transition, and for polystyrene it is 373 K (100 °C), which is why a foam cup softens in boiling water. This lab packs twenty chains of twenty beads into a box, 400 beads in all, in the manner of Kremer and Grest, and the melt opens at 450 K. Press Play and open the analysis panel: the mean squared displacement of the beads grows steadily with time, which is a liquid diffusing, and the radius of gyration of each chain fluctuates as the chain changes shape.

Ask the agent to cool the melt from 550 K to 200 K in steps of 50 K and it measures the diffusion coefficient at each step: 0.291 Å²/ps at 550 K, 0.062 at 400 K, 0.045 at 350 K, and 0.018 at 200 K. A falling coefficient by itself leaves the glass transition undefined, because a bead boxed in by its neighbors still rattles, and every 10 ps window reads a value above zero. What tells a liquid from a glass is whether the displacement keeps growing with time: in a liquid the mean squared displacement rises in proportion to the time elapsed, and in a glass it rises for a picosecond and then levels off. The agent fits that growth exponent at every step and charts it too. It falls from 0.80 at 550 K to 0.27 at 400 K, 0.15 at 350 K, and 0.01 at 200 K, where each bead moves only 0.7 Å, 0.15 of its own diameter, inside its cage. The exponent drops below 0.3 at about 410 K, give or take 40 K, and that is the glass transition of the model on a 10 ps time scale.

Polystyrene's own glass transition is 373 K, and the attraction of the bead was chosen to land near it, so the agreement is a fit; repeat runs land tens of kelvin apart. A glass transition seen in picoseconds after cooling at a trillion kelvin per second is a kinetic transition: it moves with the cooling rate, the length of the window, the depth of the bead well, and the fact that nearly half of a 400-bead blob is surface, where a bead has fewer neighbors to cage it. Hold each step longer and the measured transition moves up. Warm the melt again and the motion returns.

Conclusion

This article uses a polymer melt to illustrate a new way to learn science: students converse with an AI agent that turns a plain-language question into a molecular dynamics experiment and hands back measured data. The glass transition of a plastic becomes something students measure for themselves in a browser, along with the reason the measured value depends on how fast the melt is cooled. By pairing a physically grounded coarse-grained simulation with an AI that can set it up and run it, AIMS turns a student's curiosity into a genuine scientific investigation of the plastics of everyday life.

← Back to AIMS