스칼라장의 작용과 경계 변분

1. 좌표 목록에서 장으로

진동하는 줄에는 위치마다 변위가 있습니다. 이를 $\phi(x,t)$로 쓰면 $x$는 장소를 가리키는 표지이고 $\phi$는 그곳에서 변하는 양입니다. 온도 지도 역시 장이지만 진동하는 줄과 같은 운동 법칙을 따르지는 않습니다. 장은 기술 방식이며, 파동이나 양자 현상을 자동으로 뜻하지 않습니다.

먼저 $0\le x\le L$ 위의 실수장 하나를 다룹니다. 움직이는 표면을 그리기 전에 단위, 경계 조건, 초기 변위와 초기 속도를 정해야 합니다. 아래 계산은 기준 척도를 선택해 무차원화한 좌표와 계수를 사용합니다.

2. 국소 작용과 모형의 가정

학습용 모형을 다음과 같이 정합니다.

$$S[\phi]=\int_{t_0}^{t_1}\!dt\int_0^L\!dx\,\mathcal L, \qquad \mathcal L=\frac12\phi_t^2-\frac{c^2}{2}\phi_x^2-V(\phi), \qquad V(\phi)=\frac a2\phi^2+\frac b4\phi^4.$$

$c>0$, $b\ge0$이며 $a$는 부호를 가질 수 있는 계수입니다. 이를 물리적인 질량의 제곱이라고 자동으로 해석하지 않습니다. 외력과 감쇠가 없는 국소적·보존적 실수 스칼라 모형입니다. 학습한 신경장이 반드시 이 작용을 따른다는 뜻도 아닙니다.

작용은 한 시점의 그림이 아니라 시간에 걸친 이력 전체의 양입니다. 무차원화 전에는 에너지와 단위도 다릅니다.

3. 장을 바꾼다는 것은 무엇일까?

$\phi$와 $\phi+\epsilon\eta$를 비교합니다. 시험 함수 $\eta(x,t)$는 허용되는 경계 변화를 따라야 합니다. 시간의 양 끝에서 장을 고정하여 $\eta(x,t_0)=\eta(x,t_1)=0$으로 둡니다. $\epsilon=0$에서 미분하면

$$\delta S=\int dt\,dx\,[\phi_t\eta_t-c^2\phi_x\eta_x-V'(\phi)\eta]$$

입니다. 이는 함수에 작용하는 범함수의 방향 미분입니다. 슬라이더 하나를 바꾸는 것은 가능한 방향 하나일 뿐입니다. 정지 조건은 모든 허용 방향에서 미분이 사라져야 한다는 뜻입니다.

4. 경계항을 남겨 놓고 계산하기

부분적분하면

$$\delta S=\int dt\,dx\,[-\phi_{tt}+c^2\phi_{xx}-V'(\phi)]\eta +\left[\int dx\,\phi_t\eta\right]_{t_0}^{t_1} -\int dt\,[c^2\phi_x\eta]_{0}^{L}$$

입니다. 시간 경계항은 끝점 조건으로 사라집니다. 공간 끝점을 고정하면 $\eta=0$입니다. 주기 경계에서는 장과 변화량이 모두 주기적일 때 양 끝의 공간 항이 상쇄됩니다. 양 끝을 독립적으로 자유롭게 바꾸면 자연 경계 조건 $\phi_x=0$이 나옵니다. 서로 다른 실험이므로 경계항을 버리는 이유도 구분해야 합니다.

5. 모든 시험 방향에서 운동 방정식으로

지정한 경계 조건과 매끄러운 해를 가정하면

$$\phi_{tt}=c^2\phi_{xx}-a\phi-b\phi^3$$

을 얻습니다. 첫 항은 이웃 변위의 영향을 전달하고, 둘째 항은 국소적인 선형 반응, 셋째 항은 비선형 반응입니다. 부호는 변분 계산에서 결정됩니다. 애니메이션이 안정적으로 보이도록 임의로 선택하는 값이 아닙니다.

6. 정지한다는 말은 최소라는 뜻이 아니다

정적인 퍼텐셜의 극값과 동적인 작용의 정지는 다른 명제입니다. 자유 이차 모형의 이차 변분은

$$\delta^2S=\int dt\,dx\,[\eta_t^2-c^2\eta_x^2-a\eta^2]$$

입니다. 시간 방향으로 빠르게 변하는 허용 변화량은 이를 양수로 만들 수 있고, 공간 방향으로 빠르게 변하는 변화량은 음수로 만들 수 있습니다. 따라서 정지 이력은 안장점일 수 있습니다. 방정식 잔차의 제곱을 최소화하는 문제는 또 다른 최적화 문제이며 경계 조건과 식별 가능성을 따로 검토해야 합니다.

7. 공간 격자로 유한 모형 정하기

$N\ge3$, 간격 $h=L/N$, $q_{i+N}=q_i$인 주기 격자에서

$$L_h=h\sum_i\left[\frac{v_i^2}{2}-V(q_i)\right] -\frac{c^2}{2h}\sum_i(q_{i+1}-q_i)^2$$

로 정의합니다. 시간은 아직 연속입니다. 운동 에너지의 질량 행렬은 $hI$이며, $h=1$이 아니면 $I$가 아닙니다. 이 유한 라그랑지안을 미분하면

$$\dot q_i=v_i,\qquad \dot v_i=\frac{c^2}{h^2}(q_{i-1}-2q_i+q_{i+1})-aq_i-bq_i^3$$

입니다. 공간 이산화만으로 시간 적분 알고리즘이나 수렴 증명이 주어진 것은 아닙니다.

경계 종류에 따라 격자의 기하도 바뀝니다. 균일 Dirichlet 경계에서는 내부 $N$개 점 바깥의 값을 0으로 고정하고 $L=(N+1)h$로 둡니다. 내부 점과 고정 끝점을 연결하는 두 변도 에너지에 포함합니다. 무유속 Neumann 이산화에서는 길이 $L=Nh$의 구간에 $(i+1/2)h$인 $N$개 셀 중심을 두고 내부 $N-1$개 변만 사용합니다. 이는 유한 변분 모형의 정의이며, 모든 연속 경계 근사의 정확도를 자동으로 보장하지 않습니다.

$a=b=0$, 모든 $q_i=1$이면 주기·무유속 모형의 가속도는 0입니다. 그러나 고정값 0인 모형의 양 끝 내부 점은 $-c^2/h^2$로 가속합니다. 경계값 0과 경계 흐름 0을 단순히 ‘0 경계’라는 같은 조작으로 처리하면 안 됩니다.

8. 세 격자점을 정확히 계산하기

$N=3$, $h=c=1$, $a=b=0$, $q=(1,0,-1)$, $v=(1,0,0)$으로 둡니다. 이웃 차이는 $(-1,-1,2)$이므로 기울기 에너지는 $3$, 운동 에너지는 $1/2$입니다. 전체 에너지는 $H_h=7/2$입니다.

가속도는 $(-3,0,3)$입니다. 각 점은 가속하지만 가속도의 합은 0입니다. 이웃 상호작용은 운동을 재분배하며 각 좌표가 개별적으로 보존된다는 뜻은 아닙니다.

9. 힘의 부호를 다른 계산으로 검산하기

공간 퍼텐셜 에너지를 $U$라 하고 앞 예제에서 $q_0$만 바꿉니다.

$$U(q+\epsilon e_0)=3+3\epsilon+\epsilon^2.$$

방향 미분은 $3$입니다. 그 음수를 질량 $h=1$로 나누면 가속도 $-3$을 얻습니다. $\epsilon=1/10$일 때 두 에너지는 $331/100$과 $271/100$이며 중심 차분은 정확히 $3$입니다. 사차 퍼텐셜에서는 일반적으로 유한 차분의 절단 오차가 남으므로 같은 검산에서도 간격 의존성을 보존해야 합니다.

이 사차 격자 에너지에서는 유한 간격 오차를 정확히 제거할 수 있습니다. 좌표 방향을 $e_i$, 0이 아닌 진폭 변화량을 $\delta$라고 두면

$$\frac{\partial U}{\partial q_i} =\frac{U(q+\delta e_i)-U(q-\delta e_i)}{2\delta}-hbq_i\delta^2,$$ $$\frac{\partial^2 U}{\partial q_i^2} =\frac{U(q+\delta e_i)+U(q-\delta e_i)-2U(q)}{\delta^2}-\frac{hb\delta^2}{2}$$

입니다. $(q_i\pm\delta)^4$를 전개하여 두 보정항을 직접 확인해보세요. 이웃 기울기 에너지는 이차식이므로 여기서는 더 높은 차수의 오차를 만들지 않습니다. 두 좌표의 혼합 차분은 이웃 사이의 쌍선형 결합을 분리하며, 이 모형에서 인접하지 않은 좌표의 혼합 미분은 0입니다. 진폭 변화량 $\delta$와 공간 간격 $h$의 역할도 다릅니다.

따라서 기호 미분이나 자동 미분으로 얻은 힘·헤시안을 에너지 변화로 독립 검산할 수 있습니다. 이 정확한 항등식은 명시한 다항 에너지에 의존합니다. 임의의 학습된 퍼텐셜에 같은 보정을 적용해서는 안 됩니다.

10. 보존 법칙과 실제 시간 계산 구분하기

유한 모형에서

$$H_h=\frac h2\sum_i v_i^2+U(q),\qquad \frac{dH_h}{dt}=\sum_i\frac{\partial U}{\partial q_i}v_i+h\sum_i v_i\dot v_i=0$$

입니다. 세 점 예제의 두 항은 각각 $3$과 $-3$입니다. 상쇄 여부를 독립적인 수치 진단으로 사용할 수 있습니다. 다만 시간 적분 알고리즘은 에너지 오차를 만들 수 있습니다. 연속 시간 모형의 보존 법칙만으로 구현의 보존을 보장하지 않습니다.

연속장에서는 $e=\phi_t^2/2+c^2\phi_x^2/2+V$, $j=-c^2\phi_t\phi_x$에 대해 $\partial_te+\partial_xj=0$입니다. 따라서 구간의 에너지 변화율은 $j(0)-j(L)$이며, 경계로 드나드는 흐름을 확인해야 합니다.

11. 대칭·영모드·안정성

$a=b=0$이고 주기 경계일 때 모든 $q_i$에 같은 상수를 더해도 퍼텐셜이 변하지 않습니다. 균일 모드의 복원력은 0이고 $h\sum_i v_i$가 보존됩니다. 고정 경계는 일반적으로 이 이동 대칭을 깨뜨립니다. 이 합은 장의 값을 균일하게 옮기는 변환에 대한 정준 보존량입니다. 공간 표지 x를 이동시키는 변환에 대응하는 운동량과 자동으로 같은 양은 아닙니다. 연속장에서는 두 변환을 구분해야 합니다.

균일한 $q_i=Q$의 가속도는 $-aQ-bQ^3$입니다. $a=-1,b=1$이면 평형은 $Q=0,\pm1$입니다. 국소 퍼텐셜 곡률은 0에서 $-1$, 나머지 두 평형에서 $2$입니다. 이는 해당 퍼텐셜의 성질이며, 세 점짜리 유한 시연으로 열역학적 상전이를 입증한 것은 아닙니다.

국소 퍼텐셜의 하한과 유한 문제 전체의 안정성도 다릅니다. $a<0,b=0$이면 국소 퍼텐셜에는 하한이 없습니다. 그러나 내부 점 3개인 Dirichlet 모형에서 $h=1,c=2,a=-1$이면 전체 퍼텐셜 헤시안의 대각 성분은 $7$, 인접한 비대각 성분은 $-4$입니다. 선도 주행렬식이 $7,33,119$로 모두 양수이므로 양의 정부호입니다. 경계의 강성이 이 유한 이차 문제를 안정화하지만, 임의 구간 길이나 연속 극한의 안정성을 증명한 것은 아닙니다.

12. 시각화에서 드러내야 하는 것

변위, 속도, 이웃 기울기, 국소 힘, 에너지 항, 경계 종류를 따로 표시합니다. 모든 점을 같은 양만큼 이동시킨 경우와 비균일 교란을 비교합니다. 같은 변위 그림에서도 속도가 다르면 미래가 달라지므로 정지 화면만으로 상태를 결정할 수 없습니다.

띠나 표면은 $\phi$의 표현입니다. 세로축은 장의 값이며 새로운 물리적 공간 차원으로 오해하지 않도록 표시합니다. 움직임 외에 읽을 수 있는 표와 키보드 조작을 제공합니다. 잔차와 에너지 오차에는 사용한 이산화 방식과 시간 간격을 명시해야 합니다.

13. 이해 확인과 해설

초기 변위만으로 미래가 정해질까요? 일반적으로 아닙니다. 시간에 대한 이계 모형이므로 초기 속도와 호환되는 경계 조건도 필요합니다.

균일 변위이면 가속도가 0일까요? $aQ+bQ^3=0$일 때만 그렇습니다. 공간 미분이 사라져도 국소 퍼텐셜의 힘은 남을 수 있습니다.

수치 에너지가 보존되면 연속장 근사가 정확할까요? 아닙니다. 공간·시간 세분화, 의도한 경계 문제, 연구 주장에 필요한 관측량을 함께 검사해야 합니다.

고전 스칼라장이 곧 양자장일까요? 아닙니다. 양자화, 상태, 관측량, 정칙화에 별도 구성과 가정이 필요합니다.

14. AI 연구 연결과 참고 자료

장 모형은 구조를 가진 대리 동역학과 변분 목적함수를 설계하는 출발점입니다. 표본 궤적을 맞추는 것, 선택한 이산 방정식을 만족하는 것, 의도한 연속 문제를 근사하는 것을 구분합니다. 역문제에서는 서로 다른 퍼텐셜이나 경계 입력이 같은 희소 관측을 설명할 수 있습니다.

일반적인 장의 작용·경계 변분·Euler–Lagrange 구성은 David Tong의 Classical Field Theory 1.1·1.3절을 참고합니다. 세 점 주기 격자 계산은 이 교재의 예제입니다.