양자 확산 제한 응집의 프랙탈 특성

초록

본 연구는 고전적 확산‑제한 응집(DLA)과 양자역학적 파동함수 전개를 이용한 양자 DLA(QDLA)를 비교한다. 동일한 격자와 초기 조건에서 시뮬레이션을 수행한 결과, 양자 입자 역시 프랙탈 구조를 형성하지만, 프랙탈 차원은 초기 파동팩킷의 폭에 따라 1.43에서 2 사이로 변한다. 평균적으로는 고전적 DLA와 거의 동일한 차원(≈1.69)을 보였다.

상세 분석

이 논문은 확산‑제한 응집(DLA) 현상을 고전적 랜덤 워크와 양자역학적 파동함수 전개 두 가지 관점에서 정량적으로 비교한다. 고전적 DLA는 입자가 격자상의 네 방향으로 확률 ¼씩 이동하는 이산적 랜덤 워크로 모델링되며, 연속적인 형태로는 확산 방정식 ∂Φ/∂t = D∇²Φ 로 기술된다. 저자는 이 확산 방정식을 수치적으로 풀 때 시간 간격 Δt<1을 선택해 안정성을 확보하고, D=¼와 격자 간격을 1로 두어 고전적 랜덤 워크와 동등한 동역학을 구현한다.

양자 DLA(QDLA)는 입자를 파동함수 ψ(x,y,t) 로 기술하고, 자유 입자에 대한 시간 의존 슈뢰딩거 방정식 iħ∂ψ/∂t = -(ħ²/2m)∇²ψ 를 명시적으로 전진-후진 차분법으로 풀어 전파한다. 여기서 확산 계수는 허수값이므로, 실제 확산 방정식과는 달리 수치적 불안정성이 오히려 억제된다. 파동함수의 절댓값 제곱 |ψ|²가 입자 존재 확률이 되며, 격자상의 이미 형성된 응집체(시드)와 겹치는 위치는 확률을 0으로 강제하고 전체 파동함수를 정규화한다.

시뮬레이션은 256×256, 512×512 격자에 시드를 중앙에 두고, 입자를 원형 경계 위에 가우시안 분포(σx=σy=10)로 초기화한다. 각 입자는 최대 500 000 타임스텝까지 확산·전파되며, 매 Δt⁻¹ 단계마다 시드 인접 셀을 검사해 확률이 가장 큰 위치에 응집시킨다. 병렬 MPI 구현을 통해 16코어에서 효율적으로 계산하고, 동일 파라미터 하에 13번씩 반복해 평균값과 표준편차를 구한다.

프랙탈 차원은 질량-반경 관계 M(r)=k rᵈ의 지수 d를 로그‑로그 선형 구간에서 최소제곱법으로 추정한다. 저자는 자동으로 선형 구간을 탐지하는 Kroll 방법을 적용했으며, 인간 주관을 배제하기 위해 다중 시뮬레이션 평균을 사용한다.

결과적으로 고전적 DLA는 d=1.67±0.04, QDLA는 d=1.69±0.03을 보였으며, 초기 파동팩킷 폭을 크게 하면 d≈1.45, 작게 하면 d≈1.91로 변한다. 이는 Pietronero의 이론적 예측( n=1 → d≈2, n=2 → d≈1.43)과 일치하는 경향을 보이나, 평균값이 고전적 결과와 거의 동일한 점이 특이하다. 저자는 양자 파동함수가 고전적 확산과 유사하게 시드 주변에서 소멸·재정규화되기 때문에, 입자가 가지는 회절 효과가 충분히 큰 간격을 만들지 못하고, 결과적으로 프랙탈 차원이 크게 변하지 않는다고 해석한다.

이 연구는 양자역학적 확산‑제한 응집이 실험적 재료 과학이나 나노구조 형성에 적용될 가능성을 제시하면서, 초기 파동함수의 스케일이 프랙탈 밀도에 미치는 영향을 정량화한다는 점에서 의미가 크다. 또한, 고전적 확산 방정식과 양자 슈뢰딩거 방정식이 동일한 라플라스 형태를 갖는다는 수학적 연결고리를 이용해 두 현상을 직접 비교한 최초의 시도라 할 수 있다.