방사형 디랙 방정식 수치 해법

초록

본 논문은 방사형 디랙 방정식을 선형 스칼라 퍼텐셜, 쿨롱 벡터 퍼텐셜, 그리고 두 퍼텐셜의 혼합 형태에 대해 단계별 수치적 방법으로 풀이하고, 구체적인 에너지 스펙트럼과 파동함수를 제시한다.

상세 요약

본 연구는 방사형 디랙 방정식의 수치해법을 체계적으로 정리함으로써, 양자색역학(QCD)에서의 쿼크 구속 모델과 원자·뮤온 원자와 같은 전자기계 시스템을 동시에 다룰 수 있는 통합적 프레임워크를 제공한다. 첫 번째 단계에서는 질량이 0인 u, d 쿼크에 대해 선형 스칼라 퍼텐셜 Vₛ(r)=kr를 적용한다. 이 경우 질량항이 사라져 전자가 완전 상대론적 행동을 보이며, 상·하 컴포넌트가 강하게 결합한다. 저자들은 초기 추정값을 이용해 경계조건을 맞춘 뒤, 레이놀즈-쿠타(RK4)와 같은 4차 Runge‑Kutta 적분을 수행하고, 에너지 E를 변분적으로 조정하는 이터레이션을 도입한다. 수렴 조건은 상위 컴포넌트 g(r)와 하위 컴포넌트 f(r)의 연속성 및 정규화 조건을 동시에 만족시키는 지점에서 정의된다. 이 과정에서 근접한 고유값을 찾기 위해 이분법적(bracketing) 탐색과 뉴턴‑랩슨 보정이 결합되어 효율성을 높인다.

다음으로 질량이 있는 s, c, b 쿼크에 대해 동일한 Vₛ(r) 모델을 적용한다. 질량항 m가 추가되면서 하위 컴포넌트 f(r)의 진폭이 억제되고, 에너지 스펙트럼이 비선형적으로 변한다. 저자들은 질량에 따른 파라미터 스케일링을 제시하고, 실제 하드론 스펙트럼과 비교해 모델의 정량적 타당성을 검증한다.

세 번째 섹션에서는 시간‑같은 성분을 갖는 벡터 퍼텐셜 Vᵥ(r)=−α/r(쿨롱) 를 도입한다. 여기서는 해밀턴ian이 전통적인 수소 원자와 동일해, 분석적으로 알려진 에너지 공식 Eₙ=mc²