[Paper Review] Scaling Relations for Contour Lines of Rough Surfaces

M. A. Rajabpour, S. Mehdi Vaez Allaei|ArXiv.org|Jul 5, 2009

Advanced Numerical Analysis Techniques3 references3 citations

TL;DR

This paper proposes scaling relations for contour lines of self-affine rough surfaces using field theory and numerical simulations, revealing Zipf-like scaling in ranked loop perimeters and area/size scaling with fractal dimension D=(3−H)/2. The key contribution is a robust, system-size-dependent scaling law validated across multiple H values, offering a method to extract fractal dimensions from contour statistics.

ABSTRACT

Equilibrium and non-equilibrium growth phenomena, e.g., surface growth, generically yields self-affine distributions. Analysis of statistical properties of these distributions appears essential in understanding statistical mechanics of underlying phenomena. Here, we analyze scaling properties of the cumulative distribution of iso-height loops (i.e., contour lines) of rough self-affine surfaces in terms of loop area and system size. Inspired by the Coulomb gas methods, we find the generating function of the area of the loops. Interestingly, we find that, after sorting loops with respect to their perimeters, Zipf-like scaling relations hold for ranked loops. Numerical simulations are also provided in order to demonstrate the proposed scaling relations.

Motivation & Objective

• To establish scaling laws for the cumulative distribution of contour line areas and perimeters in self-affine rough surfaces.
• To investigate how system size influences scaling behavior in non-conformally invariant rough surfaces.
• To derive and validate Zipf-like scaling relations for ranked contour loops using field-theoretic methods and numerical simulations.
• To provide a quantitative method for estimating the fractal dimension D of contour lines from statistical properties of loop size and rank.

Proposed method

• Formulates a generating function for loop areas using Coulomb gas-inspired field theory for self-affine surfaces.
• Applies fractional calculus and inhomogeneous translational invariance to define currents related to Wilson loops, enabling derivation of scaling exponents.
• Uses a Gaussian ensemble of self-affine surfaces with power-law Fourier spectrum S(q) ∝ |q|⁻²⁽¹⁺ᴴ⁾ to model roughness.
• Performs numerical simulations on large lattices (up to 4000²) to cut surfaces at various heights and extract contour line statistics.
• Employs rank-based analysis of loop perimeters and areas, fitting to power laws lₙ ∝ n⁻ξ and Aₙ ∝ L²/n²/ᵈ to estimate scaling exponents.
• Verifies scaling relations across multiple realizations and system sizes, with error analysis for estimated exponents.

Experimental results

Research questions

• RQ1 How do the statistical properties of contour lines on self-affine rough surfaces scale with system size?
• RQ2 Can Zipf-like scaling be observed in the ranked distribution of contour loop perimeters, and what is its exponent?
• RQ3 What is the relationship between the fractal dimension D of contour lines and the roughness exponent H in non-conformally invariant surfaces?
• RQ4 How do the average area and radius of gyration of ranked contour loops scale with loop rank and system size?
• RQ5 To what extent do scaling relations hold when cutting surfaces at heights other than the mean?

Key findings

• A new scaling relation is established: ranked loop perimeters scale as lₙ ∝ n⁻ξ, with ξ = D/d, where D = (3−H)/2 and d is the spatial dimension.
• For H = 0.3, the scaling relation holds over more than two orders of magnitude in n, with numerical exponent 1.38 ± 0.03, close to theoretical 1.35.
• The average area of ranked loops scales as Aₙ ∝ L²/n²/ᵈ, and the radius of gyration as Rₙ ∝ L/n¹/ᵈ, both confirmed numerically with small errors.
• The fractal dimension D = (3−H)/2 is consistently recovered across different H values, with higher H values showing lower statistical accuracy due to fewer loops.
• Scaling relations remain robust when cutting surfaces at heights other than the mean (e.g., h = 0.1σ to 0.9σ), with no significant deviation observed.
• The study confirms that large momenta do not affect scaling, while small momenta and system size are critical, indicating non-trivial dependence on system extent.