배제 큐잉 프로세스 동역학 분석

초록

본 논문은 배제 효과가 포함된 큐잉 프로세스(EQP)를 평행 업데이트 TASEP 형태로 모델링하고, 입자 수 ⟨Nₜ⟩와 시스템 길이 ⟨Lₜ⟩의 시간 의존성을 생성함수와 도메인 월 이론으로 분석한다. 두 가지 수렴/발산 상이 존재하며, 각각 고밀도와 최대 전류 서브상으로 세분된다. 수렴 상에서는 도메인 월 예측이 Monte Carlo 시뮬레이션과 정량적으로 일치하지만, 발산 상에서는 ⟨Nₜ⟩에만 일치한다.

상세 분석

본 연구는 기존에 정적 해만 구해진 배제 큐잉 프로세스(EQP)를 동적 관점에서 확장한다. EQP는 입자들이 한 줄로 이어지는 일종의 비정상적인 TASEP이며, 시스템 길이가 입자 수에 따라 가변한다는 점이 특징이다. 저자들은 평행 업데이트 규칙을 채택해, 매 시간 단계마다 모든 입자가 동시에 이동하거나 정지하도록 설정하였다. 이때 입자들은 서로 겹칠 수 없으며, 새로운 입자는 왼쪽 경계에서 확률 α로, 오른쪽 경계에서는 확률 β로 탈출한다. 이러한 규칙은 전통적인 개방형 TASEP와 유사하지만, 시스템 길이가 동적으로 변한다는 점에서 차별화된다.

동적 특성을 분석하기 위해 두 가지 접근법을 사용한다. 첫 번째는 생성함수 기법이다. 입자 수 Nₜ와 시스템 길이 Lₜ에 대한 확률분포를 이중 생성함수 G(z, w; t)=∑_{n,l}P(Nₜ=n, Lₜ=l)zⁿwˡ 로 정의하고, 마스터 방정식에 적용해 재귀 관계를 도출한다. 이후 복소 적분과 특이점 분석을 통해 장기 한계에서 ⟨Nₜ⟩와 ⟨Lₜ⟩가 선형적으로 증가하거나 감소하는 두 가지 거시적 거동을 확인한다. 특히, α와 β의 조합에 따라 고밀도(HD)와 최대 전류(MC) 서브상이 나타나며, 각각 ⟨Nₜ⟩∝t, ⟨Lₜ⟩∝t 혹은 ⟨Nₜ⟩∝t^{1/2}, ⟨Lₜ⟩∝t^{1/2} 형태의 스케일링을 보인다.

두 번째는 도메인 월(Domain Wall) 이론이다. 시스템을 고밀도 영역과 저밀도 영역으로 구분하고, 두 영역 사이의 경계(도메인 월)가 시간에 따라 확산하거나 이동한다는 가정 하에, 월의 평균 속도와 확산계수를 α, β, 그리고 내부 흐름 J에 대한 함수로 표현한다. 이 접근법은 특히 수렴 상(⟨Nₜ⟩,⟨Lₜ⟩가 유한값으로 수렴)에서 정확한 예측을 제공한다. 도메인 월이 왼쪽으로 이동하면 시스템 길이가 감소하고, 오른쪽으로 이동하면 증가한다. 저자들은 이 이론을 바탕으로 각 서브상에서 월의 속도와 전류를 계산하고, 이를 Monte Carlo 시뮬레이션 결과와 비교하였다.

시뮬레이션은 시스템 크기를 동적으로 조정하면서 10⁶ 시간 단계 이상을 수행했으며, 평균값과 분산을 충분히 수집하였다. 결과는 다음과 같이 요약된다. (1) 수렴 상에서는 도메인 월 이론이 ⟨Nₜ⟩와 ⟨Lₜ⟩ 모두에 대해 오차가 1 % 이내로 매우 정확했다. (2) 발산 상에서는 ⟨Nₜ⟩에 대한 예측은 여전히 정확했지만, ⟨Lₜ⟩는 시뮬레이션보다 과소평가되는 경향을 보였다. 이는 월이 무한히 멀리 이동하면서 시스템 길이의 비선형 성장에 대한 추가적인 상호작용을 무시했기 때문으로 해석된다. 마지막으로, 고밀도와 최대 전류 서브상 사이의 경계는 α와 β의 임계값 α_c=β/(1+β)와 β_c=α/(1+α) 로 명시적으로 표현된다.

이러한 분석을 통해 EQP의 동적 위상도는 전통적인 TASEP와 유사하면서도, 가변 길이라는 새로운 자유도가 복잡한 스케일링과 상 전이를 야기한다는 중요한 통찰을 얻을 수 있었다.