제한된 자원과 랭뮐러 동역학을 포함한 비보존 비대칭 배제 과정의 정상 밀도와 상전이

초록

본 연구는 양쪽 끝이 유한한 입자 저장소와 연결된 랭뮐러 동역학(Lk)을 갖는 TASEP 모델을 분석한다. 평균장 이론과 대규모 몬테카를로 시뮬레이션을 통해 정상 밀도 프로파일과 상도표를 도출했으며, 기존 개방형 TASEP‑Lk 모델과는 전혀 다른 상구조가 나타난다. 특히 일부 기존 상이 사라지고, 반대로 링‑결함 모델보다 더 많은 상이 존재한다는 점이 주요 결과이다.

상세 분석

본 논문은 전통적인 개방형 TASEP에 랭뮐러 동역학(Lk)을 도입하고, 입출구 양쪽을 유한한 입자 저장소와 연결함으로써 ‘자원 제한’이라는 새로운 비보존 요소를 추가한 모델을 제시한다. 모델의 핵심 변수는 입구 효율 α, 출구 효율 β, 그리고 전체 시스템에 걸친 랭뮐러 첨가·제거 비율을 나타내는 Ω이다. 입구·출구 속도는 저장소에 남아 있는 입자 수 N_R에 비례·반비례하는 피드백 함수 f(N_R)=N_R/L 로 정의되어, 저장소가 고갈될 경우 입구가 급격히 감소하고 출구가 증가하는 비선형 효과를 만든다. 이는 기존의 무한 저장소 가정과는 근본적으로 다르며, 시스템 전체의 질량 보존을 약하게 깨뜨린다.

평균장 방정식은 연속 방정식 형태로 도출되며, 대규모(L→∞) 한계에서 랭뮐러 속도 ω=Ω/L 로 스케일링한다. 이 스케일링은 입자 이동(속도 1)과 랭뮐러 첨착·분리 과정이 동등한 시간 척도를 갖게 하여, 두 메커니즘 간의 경쟁이 정상 상태 밀도에 미치는 영향을 명확히 드러낸다. 평균장 해석을 통해 정상 밀도 ρ(x)와 전류 J는 구간별(저밀도 LD, 고밀도 HD, 최대 전류 MC) 상에 따라 서로 다른 형태의 해를 갖는다. 특히, Ω가 0에 가까울 때는 기존 개방형 TASEP와 유사한 LD‑HD‑MC 3상 구조가 나타나지만, Ω가 증가하면 새로운 혼합 상(LD‑MC‑HD, Pure MC 등)이 등장한다. 이는 랭뮐러 동역학이 시스템 전체에 걸쳐 균일하게 입자를 공급·제거함으로써, 전통적인 경계‑주도 상전이 외에 ‘부피‑주도’ 상전이가 발생하기 때문이다.

상도표는 (α,β) 평면과 (α,Ω) 평면 두 축으로 제시된다. (α,β) 평면에서는 Ω값에 따라 LD‑HD 경계가 비선형적으로 변형되고, MC 영역이 확대되며, 일부 기존의 DW(도메인 월) 상이 사라진다. 특히, Ω가 충분히 크면 β가 매우 큰 경우에도 LD‑HD 경계가 존재하지 않고 Pure MC 상만 남는다. (α,Ω) 평면에서는 β가 작을 때 LD‑HD 상이 넓게 퍼지고, β가 중간값이면 LD‑MC‑HD 혼합 상이 넓어지며, β가 매우 큰 경우에는 Ω가 증가함에 따라 순수 MC 상이 지배한다. 이러한 상구조의 차이는 저장소 피드백이 입구·출구 속도를 동시에 조절함으로써, 전류와 밀도가 서로 얽힌 비선형 관계를 만들기 때문이다.

Monte Carlo 시뮬레이션 결과는 평균장 이론이 예측한 상경계와 거의 일치한다. 시뮬레이션은 특히 경계 근처에서 발생하는 미세한 밀도 진동과, 랭뮐러 첨착·제거에 의해 발생하는 작은 충격파를 포착했으며, 이는 평균장 근사에서 무시된 고차 상관 효과를 보여준다. 또한, 저장소 크기 L에 대한 스케일링 검증을 통해 N_R/L 비율이 일정할 때 정상 상태가 L에 독립적임을 확인하였다.

결론적으로, 본 연구는 ‘자원 제한 + 랭뮐러 동역학’이라는 두 가지 비보존 메커니즘이 결합될 때, 전통적인 TASEP의 상구조가 크게 변형되고 새로운 상이 등장함을 입증한다. 이는 세포 내 모터 단백질의 제한된 ATP·리보솜 풀, 혹은 교통 시스템에서 제한된 차량 수와 도로 진입·출구 제어와 같은 실제 현상을 모델링하는 데 중요한 이론적 토대를 제공한다.