[Paper Review] Comment on "Four-component relativistic theory for NMR parameters: Unified formulation and numerical assessment of different approaches" [J. Chem. Phys. 130, 144102 (2009)]

Radosław Szmytkowski, Patrycja Stefańska|arXiv (Cornell University)|Feb 9, 2011

Advanced Chemical Physics Studies2 references5 citations

TL;DR

This comment provides a closed-form analytical solution for the magnetic dipole shielding constant σ in relativistic hydrogen-like atoms by simplifying an infinite series of gamma functions in a prior four-component relativistic NMR theory. Using hypergeometric function identities and gamma function properties, the authors derive a compact expression that matches previously known results, significantly simplifying the original complex series representation in Eq. (3).

ABSTRACT

In the paper commented on [J. Chem. Phys. 130 (2009) 144102], Cheng et al. derived a formula for the magnetic dipole shielding constant $σ$ for the Dirac one-electron atom in its ground state. That formula involves an infinite series of ratios of the Euler's gamma functions. We show that with some algebra the series may be expressed in terms of elementary functions. This leads to a simple closed-form expression for the shielding constant.

Motivation & Objective

• To simplify the complex infinite series representation of the magnetic dipole shielding constant σ in the four-component relativistic NMR theory.
• To demonstrate that the series in Eq. (3) of the referenced paper can be summed in closed form using hypergeometric function identities.
• To provide a more compact and analytically tractable expression for σ that matches known results from earlier studies.
• To validate the equivalence of the simplified expression with previously derived results by Moore, Pyper, Zhang, and Ivanov et al.

Proposed method

• Utilizes the identity ζΓ(ζ) = Γ(ζ+1) to re-express the gamma function terms in the series of Eq. (3).
• Applies the generalized hypergeometric series identity ∑[Γ(n+a₁)Γ(n+a₂)/Γ(n+b)] zⁿ/n! = [Γ(a₁)Γ(a₂)/Γ(b)] × ₂F₁(a₁,a₂;b;z) for |z| ≤ 1.
• Employs Gauss’s summation formula for hypergeometric functions at z=1: ₂F₁(a,b;c;1) = Γ(c)Γ(c−a−b)/[Γ(c−a)Γ(c−b)] under Re(c−a−b) > 0.
• Applies the identity γ₂² = γ₁² + 3 to simplify the resulting gamma function ratios.
• Performs algebraic manipulation of the resulting expressions to derive a closed-form expression for σ₊₂.
• Combines the simplified σ₊₂ with the original σ₋₁ to obtain the final closed-form expression for σ.

Experimental results

Research questions

• RQ1 Can the infinite series representation of the shielding constant σ₊₂ in the four-component relativistic NMR theory be summed in closed form?
• RQ2 Does the simplified expression for σ₊₂ match previously derived analytical results for the magnetic dipole shielding constant?
• RQ3 Can the use of hypergeometric function identities and gamma function properties lead to a more compact and analytically useful expression for σ?
• RQ4 Is the final closed-form expression for σ equivalent to those reported by Moore, Pyper and Zhang, and Ivanov et al.?
• RQ5 What is the analytical structure of the shielding constant σ when expressed in terms of elementary functions rather than special functions?

Key findings

• The infinite series in Eq. (3) for σ₊₂ is analytically summable using hypergeometric function identities and Gauss’s formula.
• The simplified expression for σ₊₂ is σ₊₂ = (2Zα²/27) × (γ₁ + 2)/(γ₁ + 1), which is significantly simpler than the original series.
• The final closed-form expression for the total shielding constant is σ = −(2Zα²/27) × (4γ₁³ + 6γ₁² − 7γ₁ − 12)/[γ₁(γ₁ + 1)(2γ₁ − 1)], derived by combining σ₋₁ and σ₊₂.
• This closed-form expression for σ is identical to results previously derived by Moore, Pyper and Zhang, and Ivanov et al., confirming consistency across different approaches.
• The derivation confirms that the original series in Eq. (3) reduces to an elementary function expression, eliminating the need for numerical summation.
• The result provides a more efficient and analytically transparent representation of the relativistic magnetic dipole shielding constant in hydrogen-like atoms.