[论文解读] Improved percolation thresholds for rods in three-dimensional boxes

Răzvan V. Florian, Zoltán Néda|arXiv (Cornell University)|Oct 3, 2001

Theoretical and Computational Physics被引用 4

一句话总结

本研究通过扩展的蒙特卡洛模拟，改进了三维各向同性取向棒体的渗透阈值估计。在细长棒极限（R/L → 0）下，验证了排除体积规则，但发现在较大长径比或棒长接近盒子尺寸时存在显著的有限尺寸效应及与该规则的偏离。

ABSTRACT

We improve our previous results for the percolation thresholds of isotropically oriented rods in three dimensional boxes. We prove again the applicability of the excluded volume rule in the slender-rod limit (radius/length -> 0). Other limits for the rod sizes are discussed and important finite-size effects are revealed.

研究动机与目标

• 通过高精度模拟，改进三维立方体中各向同性取向棒体的渗透阈值估计。
• 检验排除体积规则在细长棒极限（R/L → 0）下的有效性，并评估其在该区域之外的适用性。
• 量化棒长相对于系统盒子尺寸所引起的有限尺寸效应。
• 研究非细长棒体在L→0或L→1极限下对排除体积规则的偏离行为。
• 提供不同棒长径比和尺寸下渗透阈值的精确数值结果。

提出的方法

• 在单位立方体中模拟可渗透的带帽圆柱体（带半球形端帽的圆柱体），采用周期性边界条件以保持边缘附近的均匀性。
• 在扩展区域内采用均匀随机采样棒中心位置，以避免边缘效应。
• 对棒体采用各向同性取向采样，以确保取向的统计均匀性。
• 通过比较棒体轴线间的最小距离与2R来计算棒体交叠，对部分重叠的棒体采用分数贡献。
• 通过动态簇合并追踪簇的形成与渗透；渗透定义为从一个面贯穿到对面。
• 对每个(L,R)组合的5000次实现取Nc（渗透时的棒数）的平均值，以估计临界浓度qp = Nc/V。

实验结果

研究问题

• RQ1 对于带帽圆柱棒体，在细长棒极限（R/L → 0）下，排除体积规则（qp ∝ 1/Vex）是否成立？
• RQ2 当棒长L与盒子尺寸相当时，有限尺寸效应如何影响排除体积规则的适用性？
• RQ3 在R/L → 0和L → 0极限下，偏差s = qp·Vex − 1如何随R/L和L变化？
• RQ4 收敛到理论s = 0极限的过程如何依赖于棒长和长径比？
• RQ5 对于固定的R/L比值，能否通过L的线性外推来估计L → 0极限下的s值？

主要发现

• 在细长棒极限（R/L → 0）下，模拟证实了排除体积规则（s = 0），当L = 0.2且R/L < 0.05时，s趋近于零。
• 当R/L < 0.05且L = 0.15时，s随R/L按幂律变化，指数为0.58，表明收敛至s = 0的过程缓慢。
• 当L → 0时，s线性收敛至非零值：R/L = 0.5时为1.58，R/L = 0.25时为1.46，R/L = 0.1时为1.14，表明对排除体积规则存在持续偏离。
• 对于固定的R（R = 0.01和0.02），随着L增加趋近于1，s偏离零，表明当棒长接近盒子尺寸时规则失效。
• 当R/L比值较大且L增加时，收敛至s = 0的过程显著变慢，尤其在L ≈ 1附近更为严重。
• 即使R/L较小，有限尺寸效应也十分显著且不可忽略，尤其当L不远远小于盒子尺寸时，这削弱了排除体积规则的一般适用性。