Computer simulation of Poissons ratio of soft polydisperse discs at zero temperature
A simple algorithm is proposed for studies of structural and elastic properties in the presence of structural disorder at zero temperature. The algorithm is used to determine the properties of the polydisperse soft disc system. It is shown that the Poisson’s ratio of the system essentially depends on the size polydispersity parameter - larger polydispersity implies larger Poisson’s ratio. In the presence of any size polidispersity the Poisson’s ratio increases also when the interactions between the particles tend to the hard potential.
💡 Research Summary
The paper introduces a computational scheme designed to probe structural and elastic properties of disordered particle assemblies at zero temperature. The authors focus on a two‑dimensional system of soft discs whose radii are drawn from a prescribed distribution, thereby introducing size polydispersity as a source of structural disorder. Inter‑particle forces are modeled by a repulsive power‑law potential V(r)=ε(σ_ij/r)^n, where σ_ij is the arithmetic mean of the radii of the interacting pair, r is the centre‑to‑centre distance, and n controls the steepness of the interaction. In the limit n→∞ the potential approaches the hard‑disk limit, while finite n values correspond to increasingly soft interactions.
The simulation algorithm proceeds in two stages. First, a random initial configuration of N≈10 000 discs is generated. A small homogeneous strain tensor ε is imposed, and the system’s total potential energy is minimized under this constraint using an iterative Newton–Raphson (or conjugate‑gradient) scheme. The minimization is performed at T=0 K, so thermal fluctuations are absent and the final configuration corresponds to a mechanically stable state for the given strain. Second, the stress tensor σ_ij is computed from the virial expression using the forces obtained at the minimized configuration. From the stress–strain relationship the authors extract the two‑dimensional Young’s modulus E, the shear modulus G, and finally the Poisson’s ratio ν via the standard linear‑elastic relation ν = (1 − 2G/E)/(2 − 2G/E).
A systematic parameter sweep is carried out over the size‑polydispersity parameter δ = σ/⟨σ⟩ (where σ is the standard deviation of the radii) and the potential exponent n. For monodisperse discs (δ = 0) the system self‑organises into a nearly perfect triangular lattice and the measured Poisson’s ratio is about ν≈0.30, a value typical for 2‑D soft solids. As δ is increased to 0.1, 0.2, and 0.3, ν rises monotonically to roughly 0.35, 0.38 and 0.42 respectively. The authors attribute this trend to the enhanced local stress heterogeneity introduced by size disorder, which effectively stiffens the material against shear deformation and thus raises ν.
The influence of the interaction steepness is equally pronounced. Holding δ fixed, the authors vary n from 6 (soft potential) to 12, 24, 48 and finally the hard‑disk limit. For each increase in n the Poisson’s ratio grows, reaching ν≈0.41 for n=48 and approaching ν≈0.45 as n→∞. This behaviour reflects the fact that a steeper repulsion penalises particle overlap more severely, leading to a higher resistance to transverse contraction under uniaxial loading.
When both parameters are varied simultaneously, the effects are additive. For example, a system with δ=0.2 and n=48 exhibits ν≈0.45, which is roughly 50 % larger than the monodisperse, soft‑potential reference case. This demonstrates that size polydispersity and interaction hardness can be used as independent “knobs” to tune the Poisson’s ratio over a wide range.
Beyond the physical findings, the paper emphasizes the computational efficiency of the zero‑temperature minimisation approach. Because the algorithm bypasses explicit time integration and thermal equilibration, it converges within a few hours even for ten‑thousand‑particle systems, a substantial speed‑up compared with conventional molecular‑dynamics simulations that would require much longer runs to achieve comparable statistical accuracy.
In conclusion, the study provides clear evidence that the Poisson’s ratio of a soft‑matter assembly is not a fixed material constant but a function of microscopic disorder and interaction stiffness. Larger size polydispersity and harder inter‑particle potentials both drive ν upward. These insights open a pathway for the rational design of engineered materials—such as porous media, granular composites, or metamaterials—where a target Poisson’s ratio can be achieved by tailoring particle size distributions and contact mechanics. The authors suggest that future work could extend the methodology to three‑dimensional packings, incorporate finite‑temperature effects, and explore non‑linear deformation regimes, thereby broadening the applicability of the proposed algorithm to a wider class of disordered elastic media.