공초점 슬라이스에서 본 탄성 변동

초록

이 논문은 3차원 콜로이드 결정의 탄성 변동을 단일 공초점 평면으로 투사했을 때 나타나는 비표준 분산 관계를 이론적으로 분석한다. 투사된 2차원 슬라이스는 ω∝√q 형태의 이상적인 스케일을 보이며, 이를 통해 2차원 실험 데이터만으로 3차원 탄성 상수를 역산할 수 있는 구체적인 식을 제시한다.

상세 분석

본 연구는 3차원 입방체 결정의 작은 변형을 기술하는 변위 텐서 u_i와 변위 구배 텐서 u_{ij}를 시작점으로 한다. 라메 상수 λ, 전단계수 μ, 그리고 비등방성 파라미터 ν를 포함하는 탄성 매트릭스 D̂_{ik}(k) (식 1)를 정의하고, 그 역행렬인 정적 탄성 그린 함수 Ĝ_{jk}(k) (식 2)를 구한다. 열역학적 평균 ⟨u_i(k)u_j(−k)⟩=β^{−1}Ĝ_{ij}(k) (식 4)를 이용해 파동벡터 k에 대한 고유값 d_i(k)와 보조 변수 ω_i^2(k)=d_i(k)를 도입한다.

그 다음, 실험적으로 관찰 가능한 2차원 공초점 슬라이스를 고려한다. 3차원 그린 함수를 실공간으로 역변환한 뒤, 평면 (1,1,1) 에 수직인 단위벡터 N을 이용해 2차원 좌표 x를 3차원 벡터 \bar{x}=∑{α=1}^2 x_α r^{(α)} 로 확장한다. 투사된 2차원 그린 함수 Q{αβ}(x)는 거의 등방성 형태 Q_{αβ}(x)= (C_0/4πx)δ_{αβ}+ (C_{00}/4πx^3) x_α x_β 로 표현된다 (식 8). C_0, C_{00}는 3차원 그린 함수의 텐서 성분을 δ_{ij}−N_iN_j 및 \bar{x}_i\bar{x}_j/x^2 로 각각 적분해 얻는다 (식 9).

Fourier 변환을 수행하면 2차원 탄성 매트릭스 D̂_{αβ}(q)의 역인 Q̂_{αβ}(q)는 E_0(q)δ_{αβ}+E_{00}(q)q_αq_β/q^2 형태가 된다 (식 10). 여기서 중요한 점은 q의 크기에 대한 의존성이 1/q가 되며, 이는 전통적인 3차원 탄성 이론에서 나타나는 q^2와는 근본적으로 다르다는 것이다.

실제 계산을 단순화하기 위해 FedoroV의 등방성 근사를 도입한다. 비등방성 파라미터 ν를 반영한 유효 라메 상수 λ̃=λ+ν/5, μ̃=μ+ν/5 (식 11)를 정의하고, 이를 이용해 등방성 그린 함수 Ĝ_{ij}(k)= (1/μ̃k^2)δ_{ij}−(λ̃+μ̃)/(μ̃(λ̃+2μ̃)k^4)k_i k_j 를 얻는다 (식 12). 이 식으로부터 2차원 투사된 그린 함수의 스칼라 계수들을 직접 읽어낼 수 있으며, 최종적으로 D̂_{αβ}(q)의 고유값, 즉 효과적인 분산 관계를

ω_⊥^2 = 2 μ̃ q (전단 모드)
ω_∥^2 = 4 μ̃(λ̃+2μ̃)/(λ̃+3μ̃) q (압축 모드)

(식 14) 로 얻는다. 따라서 ω ∝ √q 스케일이 나타나며, 이는 “비데비”(non‑Debye) 형태라 부른다.

이론 검증을 위해 160×160×162 입자를 갖는 면심입방체(FCC) 구조의 분자동역학 시뮬레이션을 수행하였다. 3차원 전체에서의 분산 곡선은 ω∝k 스케일을 보였고, 이를 통해 λ, μ, ν = 42.8, 51.8, −53.8 (단위: 입자 직경, 질량, k_BT) 를 추정했다. 동일한 시스템을 (1,1,1) 평면으로 투사한 뒤 2차원 변위 상관함수를 계산하고, 식 (14)와 (11)의 파라미터를 사용해 예측 곡선을 그렸다. 시뮬레이션 결과와 이론 곡선은 3 % 이내의 오차로 일치했으며, 이는 비등방성 효과와 FedoroV 근사의 제한을 제외하고도 충분히 정확함을 보여준다.

또한, 2차원 원판(디스크) 격자를 선상으로 투사한 별도 시뮬레이션을 수행해 등방성 2차원 탄성 이론을 그대로 적용했을 때도 동일한 ω^2∝q 관계가 관찰되었다. 이는 투사된 시스템이 일반적으로 |q| 의 선형 스케일을 갖는 비표준 분산을 보인다는 점을 뒷받침한다.

결론적으로, 본 연구는 3차원 탄성 매질을 2차원 공초점 슬라이스로 관찰할 때 나타나는 비표준 분산 관계를 정확히 예측하는 이론적 프레임워크를 제공한다. 실험적으로 측정된 2차원 변위 상관함수만으로도 3차원 라메 상수를 역산할 수 있는 실용적인 식을 제시함으로써, 콜로이드 결정 및 유리 전이와 같은 연성 물질의 미세역학적 특성을 보다 효율적으로 탐구할 수 있게 된다.