보존되지 않는 자기조직화 임계성 모형
에너지가 보존되지 않는 시스템에서 임계성이 나타나는가라는 문제에 관해 재작년에 두 개의 글을 쓴 적이 있습니다.
모래쌓기 모형과 지진 모형의 차이(2007. 4. 5.)
그해 6월에는 명쾌한 결론이 나지 않은 이 문제에 대해 여러 논문을 찾아 공부한 결과를 정리하여 세미나에서 발표하기도 했습니다. 그런데 제 블로그에서도 여러 번 언급했던 미겔 무노즈(Miguel A. Munoz) 그룹에서 이 문제에 관한 준 리뷰급 논문을 올 봄에 냈습니다. 아직 저널에 게재된 건 아니고 아카이브에 올라온 걸 보실 수 있습니다: arXiv:0905.1799v2. 제가 재작년 공부하면서 봤던 논문들이 정리되어 있어서 비교적 쉽게 읽어내려갈 수 있었습니다.
이 논문의 논점을 이해하려면 꼭 봐야할 논문이 하나 더 생겼는데요, 지금 이 글에서 소개하려고 합니다. 프루스너와 젠센(줄여서 PJ로 부를게요.)이 2002년 <유로피직스 레터스(EPL)>에 낸 건데요, "A solvable non-conservative model of Self-Organised Criticality"라는 제목입니다. 이전 글들에서 옌센은 비보존 시스템에서도 SOC가 보인다는 주장을 한 그룹에 속했다는 걸 염두에 두어야겠습니다.
무노즈가 논문에서 말한대로 PJ 모형은 지진 모형(earthquake model)과 숲불 모형(forest-fire model; FFM)의 중간쯤에 있습니다. 모형을 보죠. N개의 자리가 있고 각 자리는 n개의 이웃과 연결되어 있습니다. 보존 변수는 α라고 합니다. 각 자리 i에는 에너지 zi라는 값을 갖습니다. 이 값에 따라 각 자리는 '안정(stable)', '불안정(susceptible)', '활성(active)' 상태로 나뉩니다.
아래첨자 act는 '활성', st는 '안정', c는 '불안정'을 뜻합니다. '윗줄'은 평균을 뜻하고요, pst와 pc는 각각 안정한 자리와 불안정한 자리의 밀도입니다. s는 사태의 크기이며 꺽쇠는 사태 크기의 평균을 뜻합니다. 위 좌변은 '풀기'에 의해 에너지가 줄어드는 양이며 우변의 첫째항은 '몰기'에 의해 에너지가 늘어나는 양, 둘째항은 '방아쇠 당기기'에 의해 에너지가 늘어나는 양입니다.
여기서 1/θ라는 변수가 발산한다면 우변이 발산하므로 좌변도 발산해야 하고 그건 사태 크기의 평균이 발산한다는 걸 뜻하는데요, 이는 곧 '임계성'을 나타낸다고 볼 수 있다고 합니다. 결론적으로 보존되지 않는 모형에서도 임계성이 나타날 수 있다는 거죠. 하지만 1/θ가 발산해야 한다는 조건은 특정 변수를 미세조정해야 한다는 말이므로 엄밀히 말해서 '자기조직화'는 아닙니다.
그럼 이 모형이 왜 숲불 모형과 연관되는지만 살펴보고 마치겠습니다. 숲불 모형에서 격자 위의 각 자리는 비어 있음(E), 나무 있음(T), 불이 남(F) 중 하나의 상태를 가질 수 있습니다. 빈 자리에 p의 확률로 나무가 자라고, 나무가 있는 자리에 f의 확률로 불이 납니다. 한 자리에서 불이 나면 '즉각' 이웃한 자리의 나무로 불이 옮겨붙습니다. 불이 난 자리는 빈 자리가 됩니다. 더이상 이웃한 나무가 없을 때까지 불이 번지다가 멈춥니다. 그리고 또 빈 자리에 나무가 자라고 또 나무에는 불이 붙고 퍼지기를 되풀이합니다.
p나 f나 둘 중 하나가 0이라면 나무가 모두 타버리거나 불이 전혀 나지 않거나겠죠. 임계상태가 되려면 둘 다 0보다는 커야 합니다. 그런데 p와 f가 모두 유한하다면, 나무가 매우 빨리 자랄테고 또 화재도 빈번히 일어나므로 나무가 덩어리로 생겼다 덩어리로 없어지는 활동적인 상태가 될 겁니다. 그래서 임계상태가 되려면 두 변수 모두 0에 매우 가까워야 합니다. 마지막으로, 빈 자리에 나무가 생기는 속도보다 화재가 나타나는 속도가 매우 작아야만 다양한 규모의 화재가 일어날 수 있고 그게 화재 크기의 거듭제곱 분포의 원인이 되어 임계점이라 할 수 있습니다.
PJ 모형에서 α가 0으로 가는 극한은