[논문 리뷰] Self-organized criticality in a discrete model for Smoluchowski's equation with limited aggregations

Mathieu Merle, Raoul Normand|arXiv (Cornell University)|2015. 09. 03.

Stochastic processes and statistical mechanics참고 문헌 25인용 수 4

한 줄 요약

이 논문은 제한된 응집을 고려한 스몰وخอ프스키의 응집 방정식을 위한 이산 확률 모델을 제안한다. 여기서 크기가 임계값을 초과하는 군집은 겔화되어 비활성화된다. 이 모델은 자가조직 임계성(self-organized criticality)을 나타낸다: 겔화 이전에는 군집 크기 분포가 초임계적(Galton-Watson) 분포를 따르고, 겔화 이후에는 임계적(Galton-Watson) 분포를 따르며, 시스템의 점근적 행동은 스몰وخ오프스키 방정식의 유체역학적 근사와 일치하며, 최종 농도에 대한 해석적 공식에 대한 명시적 설명을 제공한다.

ABSTRACT

We introduce and study a discrete random model for Smoluchowski's equation with limited aggregations. The latter is a model of coagulation introduced by Bertoin which may exhibit gelation. In our model, a large number of particles are initially given a prescribed number of arms. These arms are activated independently after exponential times and successively linked together. However, when the size of a cluster goes above a fixed threshold, it falls instantaneously into the gel, meaning that it no longer interacts with the other particles in solution. The concentrations in this model asymptotically obey Smoluchowski's equation with limited aggregations. In this article, we study the discrete features of this model. We are able to argue that it remains closely related with a configuration model. We specifically obtain explicit expressions for the parameters of this configuration model, which show that it is subcritical before gelation, but remains critical afterwards. As a consequence, the asymptotic distribution of a typical cluster in solution is that of a subcritical Galton-Watson tree before gelation, while it is that of a critical Galton-Watson tree after gelation. Our model therefore exhibits self-organized criticality. Our study relies heavily on the study of the configuration model in the critical window. One of the main results of the paper is an extension of a result of Janson and Luczak on the size of the largest components in this critical window. Finally, as a consequence of our explicit expressions, we provide an explanation of analytic formulas established by Normand and Zambotti regarding the limiting concentrations in Smoluchowski's equation.

연구 동기 및 목표

• 유한한 군집 크기 제한을 갖는 응집 과정을 모델링하여 물리적 겔화를 모방한다.
• 한계가 있는 응집을 따르는 스몰호프스키 방정식의 점근적 행동을 따르는 이산 확률 시스템을 수립한다.
• 군집 크기 분포 전이를 통해 수학적으로 엄밀한 방식으로 자가조직 임계성을 입증한다.
• 모델의 역학을 구성 모델(configuration model)과 연결하고, 장손과 루츠카의 임계 창(window) 결과를 확장한다.
• 노르맹드와 산보티가 유도한 최종 농도에 대한 해석적 공식에 대한 명시적 설명을 제공한다.

제안 방법

• 입자는 고정된 수의 팔을 갖는다. 이 팔들은 지수 분포를 따르는 기다림 시간 후에 독립적으로 활성화된다.
• 팔은 쌍으로 연결되어 군집을 형성하지만, 고정된 임계값을 초과하는 군집은 즉시 겔화되어 향후 응집에 참여하지 않는다.
• 모델은 구성 모델로의 사상(mapping)을 통해 분석되며, 기저의 무작위 그래프 과정에 대한 명시적 매개변수 표현식이 유도된다.
• 구성 모델의 임계 창 분석을 확장하여 겔화 이전과 이후의 최대 성분 크기를 기술한다.
• 유체역학적 근사를 확립하여 한계가 있는 응집을 갖는 스몰호프스키 방정식 수렴을 보여준다.
• 모델의 점근적 행동은 갈튼-워슨 분포 과정을 통해 보렐 분포와 연결되며, 기존의 해석적 해를 설명한다.

실험 결과

연구 질문

• RQ1 한계가 있는 응집을 고려한 스몰호프스키 방정식에 대한 이산 확률 모델은 자가조직 임계성을 나타내는가?
• RQ2 겔화 임계점에서 군집 크기 분포는 초임계적에서 임계적 행동으로 어떻게 변화하는가?
• RQ3 구성 모델 프레임워크는 이 응집 모델에서 임계 창 행동을 포괄하도록 확장될 수 있는가?
• RQ4 모델의 매개변수와 다항 핵을 갖는 스몰호프스키 방정식의 알려진 해석적 해 사이의 정확한 연결은 무엇인가?
• RQ5 겔화 시점에서의 시스템 전이가 군집 크기의 점근적 분포에 어떤 영향을 미치는가?

주요 결과

• 모델은 구성 모델과 밀접하게 관련되어 있으며, 겔화 시점에서 초임계적에서 임계적 행동로의 전이를 보여주는 명시적 매개변수 표현식을 갖는다.
• 겔화 이전에는 일반 군집 크기의 분포가 초임계적 갈튼-워슨 트리의 분포와 일치한다.
• 겔화 이후에는 일반 군집 크기의 분포가 임계적 갈튼-워슨 트리의 분포와 일치하여 자가조직 임계성이 확인된다.
• 논문은 구성 모델의 임계 창에서 최대 성분의 크기를 기술하는 장손과 루츠카의 결과를 이 응집 맥락으로 확장한다.
• 모델의 점근적 농도 프로파일은 스몰호프스키 방정정의 알려진 해석적 해와 일치하며, 특히 포아송 갈튼-워슨 과정의 총 후손 수에 대한 보렐 분포를 포함한다.
• 질량 분포의 尾행동은 겔화 이전에는 지수 감소를 보이다가 겔화 이후에는 지수 $-1/2$ 의 거듭제곱 감소를 보이며, 임계 시스템과 일치한다.