The Flat Phase of Crystalline Membranes

TL;DR

This study uses high-statistics Monte Carlo simulations to confirm the existence of a stable flat phase in crystalline membranes, revealing how simple models can generate necessary elastic constants and determining critical exponents.

Contribution

It demonstrates that a simple spring and bending energy model can dynamically produce elastic constants, confirming the flat phase in crystalline membranes through extensive simulations.

Findings

Confirmed the flat phase in crystalline membranes.

Measured critical exponents: ν ≈ 0.95, ζ ≈ 0.64, η_u ≈ 0.50.

Results support theoretical scaling relations.

Abstract

We present the results of a high-statistics Monte Carlo simulation of a phantom crystalline (fixed-connectivity) membrane with free boundary. We verify the existence of a flat phase by examining lattices of size up to $128^2$. The Hamiltonian of the model is the sum of a simple spring pair potential, with no hard-core repulsion, and bending energy. The only free parameter is the the bending rigidity $\kappa$. In-plane elastic constants are not explicitly introduced. We obtain the remarkable result that this simple model dynamically generates the elastic constants required to stabilise the flat phase. We present measurements of the size (Flory) exponent $\nu$ and the roughness exponent $\zeta$. We also determine the critical exponents $\eta$ and $\eta_u$ describing the scale dependence of the bending rigidity ($\kappa(q) \sim q^{-\eta}$) and the induced elastic constants ($\lambda(q) \sim \mu(q) \sim q^{\eta_u}$). At bending rigidity $\kappa = 1.1$, we find $\nu = 0.95(5)$ (Hausdorff dimension $d_H = 2/\nu = 2.1(1)$), $\zeta = 0.64(2)$ and $\eta_u = 0.50(1)$. These results are consistent with the scaling relation $\zeta = (2+\eta_u)/4$. The additional scaling relation $\eta = 2(1-\zeta)$ implies $\eta = 0.72(4)$. A direct measurement of $\eta$ from the power-law decay of the normal-normal correlation function yields $\eta \approx 0.6$ on the $128^2$ lattice.