고밀도 단일파일 운송에서 나타나는 비평형 위상 전이와 솔리톤 전류

초록

본 연구는 고밀도에서 주기적 포텐셜을 가로지르는 입자들의 단일파일 확산을 모델링한 BASEP에서, 입자 밀도가 임계값을 초과하면 열 활성화가 아닌 솔리톤(단일 클러스터 파동)에 의해 전류가 급격히 증가하는 비평형 위상 전이를 발견하였다. 전류‑밀도 관계와 전류 변동의 보편성 클래스(KPZ↔EW) 변화가 동시에 나타난다.

상세 분석

이 논문은 Brownian Asymmetric Simple Exclusion Process (BASEP)를 이용해, 직경 σ인 강체 입자들이 외부 구동력 f와 주기적 포텐셜 U(x)=U₀/2 cos(2πx/λ)를 따라 이동하는 1차원 단일파일 시스템을 연구한다. 입자들은 과잉감쇠 브라운 운동을 하며, 하드-스피어 제약 |x_i−x_j|≥σ을 만족한다. 시스템은 주기적 경계조건을 갖고, 입자 밀도 ρ=N/L으로 정의한다. 주요 관측량은 순간 전류 밀도 J(t)= (1/L)∑_i ẋ_i(t)이며, 장시간 평균 J(ρ)=⟨J(t)⟩을 분석한다.

0노이즈(D=0) 한계에서, ρ가 특정 임계값 ρ_c 이하이면 입자들은 포텐셜 우물에 고정돼 전류가 0이다. ρ>1이면 여러 입자가 같은 우물에 겹쳐 클러스터를 형성하고, 이 클러스터가 포텐셜 장벽을 ‘전달’하면서 전체 시스템을 전진시킨다. 핵심은 ‘기본 클러스터’ n_b(σ,f)로, 이는 포텐셜 최소점에 안정적으로 존재할 수 있는 최대 입자 수이며, Euler의 토시엔트 함수와 최소 잔여 자유공간 원칙을 통해 결정된다. 임계 밀도는

 ρ_c = n_b ⌈n_b σ⌉⁻¹

이며, 여기서 ⌈·⌉는 올림함수이다. ρ>ρ_c이면 솔리톤이 주기적으로 생성·소멸하면서 한 주기 τ_sol당 전체 입자 변위 Δ=λ를 만든다. 솔리톤 수 N_sol는 (ρ−ρ_c)·⌈n_b σ⌉에 비례하므로 전류는

 J = (ρ−ρ_c) v_sol, v_sol = ⌈n_b σ⌉/τ_sol

와 같이 선형적으로 증가한다. 작은 σ(예: σ=3/25)에서는 ⌈n_b σ⌉가 2 이상이 되어 ρ가 ρ* = ρ_c + 1/⌈n_b σ⌉²를 초과하면 전류 기울기가 급격히 변하는 ‘kink’가 나타난다.

노이즈(D>0)가 존재하면 열 활성화 전이로 인해 ρ<ρ_c에서도 미세한 전류가 흐르지만, 고밀도 영역에서는 솔리톤 매개 전류가 지배적이다. D가 커질수록 전류‑밀도 곡선이 부드러워지고, 특히 작은 σ에서는 전이 신호가 사라진다.

전류 변동의 보편성 클래스는 전류‑밀도 곡선의 곡률 J’’(ρ)에 따라 달라진다. J’’≈0이면 KPZ(α=4/3) 스케일링, J’’≠0이면 EW(α=3/2) 스케일링이 나타난다. 시뮬레이션 결과, ρ≈3.3(σ=1/4, f=0.1, D=0.05)에서는 J’’≈0에 가까워 EW 클래스(‑t⁻³/₂)로 전환되고, ρ≈3.05에서는 J’’>0이므로 KPZ(‑t⁻⁴/₃)로 돌아간다. 이는 솔리톤 전파가 전류 변동을 ‘무편향 확산’처럼 만들면서 비평형 KPZ 특성을 소멸시키는 메커니즘을 보여준다.

전반적인 위상도는 σ‑ρ, f‑ρ 평면에서 ‘jammed(빨강)’와 ‘current-carrying(초록)’ 영역으로 구분된다. σ가 0.35~0.8 구간에서는 n_b가 불규칙하게 변해 ρ_c에 ‘딥’(dip) 구조가 나타나며, 큰 구동력 f에서는 n_b가 단조 감소해 ρ_c가 단계적으로 감소한다.

이러한 결과는 동일한 포텐셜과 경계조건을 갖는 1차원 시스템에서도, 외부 구동과 입자 크기에 따라 전혀 다른 집합적 동역학(정지, 열 활성화, 솔리톤 전파)이 나타날 수 있음을 입증한다. 특히, 솔리톤에 의한 고전류 상태는 전통적인 ‘오픈’ 시스템에서 요구되는 입구·출구와 같은 비균일성 없이도 순수히 내부 상호작용만으로 비평형 위상 전이를 일으킬 수 있음을 강조한다.