[Paper Review] A physicist's guide to the solution of Kummer's equation and confluent hypergeometric functions

W. N. Mathews, M. A. Esrick+2 more|arXiv (Cornell University)|Nov 8, 2021

Quantum Mechanics and Non-Hermitian Physics49 references27 citations

TL;DR

This paper provides a comprehensive, case-by-case guide to solving Kummer's equation (the confluent hypergeometric equation) by systematically identifying the correct pair of linearly independent solutions for all parameter combinations of a and b, especially when a, b, or a−b are integers—cases where standard treatments fail. It resolves long-standing ambiguities in textbook physics by showing when standard functions like M(a,b,z), e^z M(1+a−b,2−b,z), and U(a,b,z) are insufficient, and derives explicit second solutions using logarithmic terms when necessary, validated through a detailed analysis of the hydrogen atom problem.

ABSTRACT

The confluent hypergeometric equation, also known as Kummer's equation, is one of the most important differential equations in physics, chemistry, and engineering. Its two power series solutions are the Kummer function, M(a,b,z), often referred to as the confluent hypergeometric function of the first kind, and z^{1-b}M(1+a-b,2-b,z), where a and b are parameters that appear in the differential equation. A third function, the Tricomi function, U(a,b,z), sometimes referred to as the confluent hypergeometric function of the second kind, is also a solution of the confluent hypergeometric equation that is routinely used. All three of these functions must be considered in a search for two linearly independent solutions of the confluent hypergeometric equation. There are situations, when a, b, and a - b are integers, where one of these functions is not defined, or two of the functions are not linearly independent, or one of the linearly independent solutions of the differential equation is different from these three functions. Many of these special cases correspond precisely to cases needed to solve physics problems. This leads to significant confusion about how to work with confluent hypergeometric equations, in spite of authoritative references such as the NIST Digital Library of Mathematical Functions. Here, we carefully describe all of the different cases one has to consider and what the explicit formulas are for the two linearly independent solutions of the confluent hypergeometric equation. Our results are summarized in Table I in Section 3. As an example, we use these solutions to study the bound states of the hydrogenic atom, going beyond the standard treatment in textbooks. We also briefly consider the cutoff Coulomb potential. We hope that this guide will aid physics instruction that involves the confluent hypergeometric differential equation.

Motivation & Objective

• To resolve widespread confusion in physics textbooks about the correct linearly independent solutions of the confluent hypergeometric equation when parameters a, b, or a−b are integers.
• To identify and systematically classify all cases where the standard solutions M(a,b,z), e^z M(1+a−b,2−b,z), and U(a,b,z) are not linearly independent or not defined.
• To provide a complete, user-friendly reference table (Table 1) that maps all parameter combinations to the correct pair of linearly independent solutions.
• To correct the standard textbook treatment of the hydrogen atom bound states by showing that the conventional solution is incomplete or incorrect in certain parameter regimes.
• To demonstrate the necessity of logarithmic terms in solutions when standard functions fail, using the hydrogen atom and cutoff Coulomb potential as concrete examples.

Proposed method

• The authors perform a systematic case analysis based on the integer or non-integer nature of a, b, and a−b, using the NIST Digital Library of Mathematical Functions (DLMF) as the primary reference for functional identities and analytic continuations.
• They derive explicit formulas for the second linearly independent solution in cases where the standard three functions (M, e^z M, U) fail to be independent or are undefined, particularly when a, b, or a−b are integers.
• The method involves identifying critical parameter thresholds (e.g., b ∈ Z≤0, a ∈ Z≤0, a−b = −(1+q)) and applying known DLMF identities (e.g., 13.2.7, 13.2.8, 13.2.9) to determine valid solution pairs.
• For cases requiring logarithmic terms (e.g., when U(a,b,z) contains ln z), they derive explicit series expressions with psi functions, such as in equations (A.1) and (A.2), which are not present in standard treatments.
• The solution strategy is summarized in a comprehensive decision table (Table 1) that maps parameter conditions to the correct solution pair.
• The method is validated by re-analyzing the hydrogen atom problem, showing that the standard textbook solution is incomplete when a or b are integers, and deriving the correct solution using the proposed framework.

Experimental results

Research questions

• RQ1 When do the standard confluent hypergeometric functions M(a,b,z), e^z M(1+a−b,2−b,z), and U(a,b,z) fail to provide two linearly independent solutions for Kummer’s equation?
• RQ2 What explicit form must the second solution take when the standard functions are not linearly independent or undefined, particularly in cases where a, b, or a−b are integers?
• RQ3 How does the correct solution to the confluent hypergeometric equation differ from the textbook treatment of the hydrogen atom, and what is the corrected bound state wavefunction?
• RQ4 What role do logarithmic terms (ln z) play in the solutions, and under what parameter conditions do they appear?
• RQ5 Can a unified, case-by-case procedure be constructed to always yield two linearly independent solutions for any real or complex a and b?

Key findings

• The paper identifies six main cases based on the integer or non-integer nature of a, b, and a−b, and provides a complete classification of which solution pairs are valid in each.
• In cases where b ∈ Z≤0 and a ∈ Z≤0 with m > n (i.e., a = −m, b = −n, m > n), the standard functions M and U fail to be independent, and a new solution involving logarithmic terms is required, explicitly derived in equation (A.1).
• For the hydrogen atom problem with integer quantum numbers, the standard textbook solution is incomplete; the correct solution requires the use of logarithmic terms when a = −n, b = 1+n, as shown in equation (A.2).
• The Tricomi function U(a,b,z) contains logarithmic terms only when b ∈ Z≤0 and a ∈ Z≤0, but not when b ∈ Z>0, which resolves a long-standing confusion in the literature.
• The paper demonstrates that in cases where a ∈ Z>0 and b ∈ Z>0 with m < n (i.e., a = 1+m, b = 1+n, m < n), the function U(a,b,z) is well-defined and free of logarithmic terms, contrary to some expectations.
• The authors derive explicit second solutions with psi functions and logarithmic terms in cases where the standard three functions are not independent, such as in equations (A.1) and (A.2), which are not present in standard references or textbooks.