[논문 리뷰] On variational solution of the four-body Santilli-Shillady model of $H_2$ molecule

A. K. Aringazin|ArXiv.org|2000. 01. 24.

Statistical Mechanics and Entropy참고 문헌 3인용 수 4

한 줄 요약

이 논문은 H₂ 분자의 4체 산틸리-실라디 이소화학적 모델에 릿츠 변분법을 적용하여 전자 상관관계와 이소전자우밍 형성을 기술하기 위해 단거리형 헐텐 포텐셜을 도입한다. 결과적으로 지수형 스크리닝된 쿠론트 포텐셜의 변형은 실험적 결합 에너지와 결합 길이를 1–6%의 이소전자우밍 상기 기여도로 잘 맞추며, 헐텐 포텐셜은 이소전자우밍 반지름에 대해 이산적인 값을 강제하여 양자화된 안정성을 나타낸다.

ABSTRACT

In this paper, we apply Ritz variational approach to a new isochemical model of $H_2$ molecule suggested by Santilli and Shillady. We studied Gaussian, $V_g$, and exponential, $V_e$, screened Coloumb potential {\it approximations}, as well as the original Hulten potential, $V_h$, case. Both the Coloumb and exchange integrals have been calculated only for $V_e$ owing to Gegenbauer expansion while for $V_g$ and $V_h$ cases we achieved analytical results only for the Coloumb integrals. We conclude that the $V_e$-based model is capable to fit experimental data on $H_2$ molecule in confirmation of the results of numerical HFR approach by Santilli and Shillady. Also, we achieved the energy-based estimation of the weight of the isoelectronium phase which is appeared to be of the order of 1%...6%, for the case of $V_e$. However, we note that this is {\it not} the result corresponding to the original Santilli-Shillady model, which is based on the Hulten potential $V_h$. An interesting result is that in order to prevent divergency of the Coloumb integral for $V_h$ the correlation length parameter $r_c$ should run discrete set of values. We used this condition in our $V_e$-based model.

연구 동기 및 목표

• 전삼체 산틸리-실라디 모델을 전자 상관관계를 헐텐 포텐셜로 포함한 완전한 4체 시스템으로 확장하기 위해.
• 4체 이소화학적 H₂ 모델에서 기초 상태 에너지와 결합 길이를 계산하기 위해 릿츠 변분법을 적용하기 위해.
• 이소전자우밍 상기의 역할을 에너지 기여도와 안정성 가중치를 추정하여 분석하기 위해.
• 헐텐 포텐셜이 이소전자우밍 효과적 반지름에 대해 이산적, 양자화된 값을 유도하는지 조사하기 위해.
• 실험적 H₂ 데이터에 대한 적합도를 비교하기 위해 헐텐 포텐셜과 지수형/가우시안 스크리닝된 쿠론트 포텐셜 간의 결과를 비교하기 위해.

제안 방법

• 총 에너지를 최소화하기 위해 γ와 ρ = γR로 매개변수화된 수소 유사 파동함수를 사용한 릿츠 변분 접근법을 적용하였다.
• 필요한 경우 게겐바우어 전개를 사용하여 헐텐, 지수형 스크리닝, 가우시안 스크리닝된 쿠론트 포텐셜에 대해 쿠론트 및 교환 적분을 해석적으로 계산하였다.
• 지수형 스크리닝된 포텐셜의 경우 변분 매개변수(γ, ρ, V₀, r_c)를 최적화하기 위해 수치적 피팅을 사용하였다.
• 헐텐 포텐셜을 단거리 상관 상호작용으로 간주하여 V₀와 r_c 매개변수를 설정하였으며, 이소전자우밍을 상관된 전자 쌍으로 모델링하였다.
• 스크리닝된 포텐셜 대비 쿠론트 포텐셜의 상대적 에너지 기여도를 통해 이소전자우밍 상기의 기여도를 추정하였다.
• 최적 에너지와 결합 길이를 결정하기 위해 r_c에 대해 60개의 값에 대해 수치 분석을 수행하였다.

실험 결과

연구 질문

• RQ1 헐텐 포텐셜을 포함한 4체 산틸리-실라디 모델이 H₂ 분자의 실험적 결합 에너지와 결합 길이를 재현할 수 있는가?
• RQ2 헐텐 포텐셜이 이소전자우밍 효과적 반지름 r_c에 대해 허용 가능한 이산적 값 집합을 유도하는가?
• RQ3 지수형 스크리닝된 포텐셜 모델 하에서 H₂ 분자에서 이소전자우밍 상기의 상대적 에너지 기여도(가중치)는 얼마인가?
• RQ4 지수형 스크리닝된 포텐셜의 결과와 헐텐 포텐셜의 결과는 실험 데이터 적합도 측면에서 어떻게 비교되는가?
• RQ5 이소전자우밍 상기는 분자의 기초 상태 에너지에 유의미하게 기여할 정도로 안정적인가?

주요 결과

• 지수형 스크리닝된 쿠론트 포텐셜 모델은 실험적 H₂ 결합 에너지(−1.174474 a.u.)와 결합 길이(1.4011 a.u.)를 잘 맞추며, R_opt ≈ 1.4013 a.u.에서 최소 에너지 −1.1432 a.u.를 달성한다.
• 지수형 스크리닝된 포텐셜의 경우 이소전자우밍 상기의 기여도는 스크리닝된 포텐셜의 상대적 에너지 기여도로 정의되며, 1–6%로 추정된다.
• 헐텐 포텐셜의 경우 변분 에너지는 이소전자우밍 반지름 r_c를 이산적인 값 집합으로 제한하여, 이는 양자화된 비연속적 안정성 조건을 나타낸다.
• 모든 경우에서 이소전자우밍의 효과적 반지름은 0.25 a.u. 이하로 제약을 받으며, 이는 밀도가 높고 국소화된 전자 쌍임을 시사한다.
• 헐텐 포텐셜의 경우 독특한 특징이 존재한다: 변분 과정에서 r_c는 반드시 특정한 이산값들만을 취할 수 있으며, 이는 지수형 또는 가우시안 스크리닝 모델에서는 관찰되지 않는 성질이다.
• 결과는 초보적이다. 완전한 4체 모델은 헐텐 포텐셜에 기반하고 있으며, 헐텐 및 가우시안 경우의 교환 적분은 여전히 해석적으로 다룰 수 없다.