LIAM MCGRATH / PHYSICSOPENMM · MOLECULAR EXPLORER

MATTER, IN MOTION

DipoleAtlas.

From locked-in directions to a cloud of possibilities. Explore the changing geometry of molecular order.

μ

REAL SIMULATION
ILLUSTRATIVE MATERIALS

A model, not a material prediction. Short runs can retain their starting structure. These visuals do not certify a phase transition.

Controls select saved simulations—not a live force-field calculation. Changing the diagnostic or display does not change the dynamics. Link to this selection ↗ · Theory & design ↗

MULTIPOLAR ORDER · O₆—
FCC REFERENCE ORDER—
DENSITY · g/cm³—
POTENTIAL · kJ/mol molecule—
EMPIRICAL POTENTIAL · selected T/PUnavailable

Model energy ↔ experimental context

Literature vaporization enthalpy

Reference source ↗

Same selected T and P?

Density anchor

Density source ↗

ΔHᵥₐₚ is not absolute potential energy. No missing value is filled with zero or extrapolated. Reference states, assumptions and comparison limits ↗

The molecular box

Drag to rotate · charged-model arrows point −q → +q; uncharged arrows show the geometric axis · site separation/arrows enlarged for visibility, not to scale · color encodes azimuth φ · box normalized to view (check density for expansion) · sampled frames, not physical playback speed

TEMPERATURE × PRESSURE

Response surface

Beyond orientation

Fluctuation responses and MSD slope are exploratory estimates from short, correlated trajectories—not converged heat capacities, diffusion coefficients or transition temperatures. How to interpret them ↗

Download figure book ↗All records · ZIPSummary · CSVRun metadata · JSONRead theory ↗Theory · PDF

READ THE MODEL

Order is a clue, not a verdict.

Why more than a single angle?

Six narrow directional populations can have almost no net polarization. We retain spherical-harmonic degrees 1–6 using pairwise Legendre polynomials. O₆ combines their positive, self-pair-corrected powers. It detects several kinds of orientational order, not an exact number of Gaussian modes.

What would suggest melting?

Loss of crystal order alongside sustained molecular displacement and a change in energy or density. A dilute state additionally needs strong density loss. Dipole disorder alone can also be a rotating solid. Short trajectories and hysteresis can be misleading.

What this run cannot claim

Qualitative two-site polar Lennard-Jones surrogates, not atomistic models of the named chemicals. Interaction scales are illustrative; hydrogen bonds, molecular flexibility, chemical reactions and decomposition are absent. Short finite periodic-box heating/cooling runs do not establish melting points, boiling points, vapour pressures or phase boundaries. High-temperature results are model behavior, not predictions for stable real molecules.

MSD removes affine box scaling and resets at each plateau. Pooled angular KL has positive finite-sample bias. The multipolar display score also has a small positive noise floor; signed spectra are retained in CSV. Browser preview is capped at 120 frames and 256 molecules per condition; full saved samples remain in the ZIP.