We derive two distinct asymptotic expansions for the zeros $j_{ν,k}^{(n)}$ of the $n$-th derivative of Bessel function $J_ν^{(n)}(x)$. The first is a McMahon-type expansion for the case when $k \to \infty$ with fixed $ν$, for which we also establish an explicit error bound. The second addresses the case when $ν\to \infty$ with fixed $k$ and it involves the zeros of Airy functions and their derivatives. These results extend and refine the classical work of Wong, Lang, and Olver on the zeros of Bessel functions. In the course of obtaining our main results, we also generalize several auxiliary results, which in turn provide a broader framework for the study of zeros of special functions.
논문 ID : 2510.12353제목 : Bessel 함수 도함수의 영점의 점근적 거동저자 : Árpád Baricz, Pranav Kumar, Saminathan Ponnusamy분류 : math.CA (고전 해석 및 상미분방정식)발표 시간 : 2025년 10월 14일논문 링크 : https://arxiv.org/abs/2510.12353v1 본 논문은 Bessel 함수의 n차 도함수 J ν ( n ) ( x ) J_\nu^{(n)}(x) J ν ( n ) ( x ) 의 영점 j ν , k ( n ) j_{\nu,k}^{(n)} j ν , k ( n ) 에 대한 두 가지 서로 다른 점근 전개식을 유도한다. 첫 번째는 k → ∞ k \to \infty k → ∞ 이고 ν \nu ν 가 고정일 때의 McMahon 형 전개이며, 명시적 오차 한계를 설정한다. 두 번째는 ν → ∞ \nu \to \infty ν → ∞ 이고 k k k 가 고정인 경우를 다루며, Airy 함수 및 그 도함수의 영점을 포함한다. 이러한 결과는 Bessel 함수 영점에 관한 Wong, Lang, Olver의 고전적 업적을 확장하고 개선한다. 주요 결과를 얻는 과정에서 여러 보조 결과를 일반화하여 특수 함수 영점 연구에 대한 더 광범위한 틀을 제공한다.
Bessel 함수 영점은 응용 수학 및 수학 물리학의 다양한 문제에서 중요한 역할을 한다. 특히 양자역학, 산란 이론, 파동 전파 등 관련 분야에서 핵심적이다. 광범위한 적용 가능성으로 인해 Bessel 함수 및 그 영점의 점근적 거동을 이해하는 것이 상당한 관심을 받고 있다.
고전적 McMahon 전개의 제한 : 초기 McMahon의 Bessel 함수 J ν ( x ) J_\nu(x) J ν ( x ) 영점에 관한 연구는 주로 함수 자체에 초점을 맞추었으며, 도함수는 다루지 않았다.Wong과 Lang의 업적 : 이들은 이러한 아이디어를 확장하여 J ν ′ ′ ( x ) J''_\nu(x) J ν ′′ ( x ) 영점의 점근 전개를 얻었으나, 2차 도함수로만 제한되었다.연구 공백 : n차 Bessel 함수 도함수 영점의 성질을 이해하는 데 있어 상당한 진전이 있었음에도 불구하고, 이러한 영점의 점근 전개는 문헌에서 아직 탐구되지 않았다.본 논문은 이러한 공백을 메우기 위해 McMahon 형 점근 전개를 제시하고 ν → ∞ \nu \to \infty ν → ∞ 이고 n이 고정인 경우를 분석하는 것을 목표로 한다. 연구는 McMahon이 도입한 표준 방법을 사용하여 J ν ( n ) ( x ) J_\nu^{(n)}(x) J ν ( n ) ( x ) 의 점근 전개를 유도하고 k → ∞ k \to \infty k → ∞ 일 때 j ν , k ( n ) j_{\nu,k}^{(n)} j ν , k ( n ) 의 점근적 거동을 설정한다.
McMahon 형 점근 전개 수립 : Bessel 함수 n차 도함수의 영점에 대해 k → ∞ k \to \infty k → ∞ 이고 ν \nu ν 가 고정일 때의 점근 전개식을 제시한다.명시적 오차 한계 제공 : McMahon 형 전개에 대한 명확한 오차 추정을 설정하여 근사의 정확도를 정량화한다.대 매개변수 점근 전개 유도 : ν → ∞ \nu \to \infty ν → ∞ 이고 k k k 가 고정인 경우를 다루며, Airy 함수 및 그 도함수의 영점을 포함한다.보조 결과 일반화 : 주요 결과를 얻는 과정에서 여러 기초 이론 결과를 일반화하여 특수 함수 영점 연구에 더 광범위한 틀을 제공한다.고전 이론 확장 : Wong, Lang, Olver의 고전적 업적을 Bessel 함수 자체에서 그 임의 차수 도함수로 확장한다.제1종 Bessel 함수 J ν ( x ) J_\nu(x) J ν ( x ) 의 n차 도함수 J ν ( n ) ( x ) J_\nu^{(n)}(x) J ν ( n ) ( x ) 의 양의 영점 j ν , k ( n ) j_{\nu,k}^{(n)} j ν , k ( n ) 의 점근적 거동을 연구한다. 여기서:
ν > 0 \nu > 0 ν > 0 은 Bessel 함수의 차수n ≥ 1 n \geq 1 n ≥ 1 은 도함수의 차수k k k 는 영점의 인덱스정리 2.1 : n ∈ N 0 n \in \mathbb{N}_0 n ∈ N 0 이고 x x x 가 클 때, Bessel 함수의 2n차 도함수는 다음과 같이 표현된다:
1 2 π x J ν ( 2 n ) ( x ) = cos ( x − ν π 2 − π 4 ) τ ν ( 2 n ) ( x ) + sin ( x − ν π 2 − π 4 ) θ ν ( 2 n ) ( x ) \sqrt{\frac{1}{2}\pi x} J_\nu^{(2n)}(x) = \cos\left(x - \frac{\nu\pi}{2} - \frac{\pi}{4}\right)\tau_\nu^{(2n)}(x) + \sin\left(x - \frac{\nu\pi}{2} - \frac{\pi}{4}\right)\theta_\nu^{(2n)}(x) 2 1 π x J ν ( 2 n ) ( x ) = cos ( x − 2 ν π − 4 π ) τ ν ( 2 n ) ( x ) + sin ( x − 2 ν π − 4 π ) θ ν ( 2 n ) ( x )
여기서:
τ ν ( 2 n ) ( x ) = ∑ m = 0 ∞ α 2 m , 2 n x 2 m , θ ν ( 2 n ) ( x ) = ∑ m = 0 ∞ α 2 m + 1 , 2 n x 2 m + 1 \tau_\nu^{(2n)}(x) = \sum_{m=0}^{\infty} \frac{\alpha_{2m,2n}}{x^{2m}}, \quad \theta_\nu^{(2n)}(x) = \sum_{m=0}^{\infty} \frac{\alpha_{2m+1,2n}}{x^{2m+1}} τ ν ( 2 n ) ( x ) = ∑ m = 0 ∞ x 2 m α 2 m , 2 n , θ ν ( 2 n ) ( x ) = ∑ m = 0 ∞ x 2 m + 1 α 2 m + 1 , 2 n
계수는 다음 점화식을 만족한다:
α 2 m , 2 n = α 2 m , 2 n − 1 − 4 m − 1 2 α 2 m − 1 , 2 n − 1 \alpha_{2m,2n} = \alpha_{2m,2n-1} - \frac{4m-1}{2}\alpha_{2m-1,2n-1} α 2 m , 2 n = α 2 m , 2 n − 1 − 2 4 m − 1 α 2 m − 1 , 2 n − 1 α 2 m + 1 , 2 n = − α 2 m − 1 , 2 n − 1 − 4 m + 1 2 α 2 m , 2 n − 1 \alpha_{2m+1,2n} = -\alpha_{2m-1,2n-1} - \frac{4m+1}{2}\alpha_{2m,2n-1} α 2 m + 1 , 2 n = − α 2 m − 1 , 2 n − 1 − 2 4 m + 1 α 2 m , 2 n − 1
역변환 기법을 통해 영점의 McMahon 형 전개를 얻는다:
짝수 차수 도함수 (J ν ( 2 m ) J_\nu^{(2m)} J ν ( 2 m ) 의 영점):
j ν , k + 1 ( 2 m ) = α + α 1 , 2 m α 0 , 2 m ⋅ 1 α + O ( 1 α 3 ) j_{\nu,k+1}^{(2m)} = \alpha + \frac{\alpha_{1,2m}}{\alpha_{0,2m}} \cdot \frac{1}{\alpha} + O\left(\frac{1}{\alpha^3}\right) j ν , k + 1 ( 2 m ) = α + α 0 , 2 m α 1 , 2 m ⋅ α 1 + O ( α 3 1 )
여기서 α = k π + ν π 2 − π 4 \alpha = k\pi + \frac{\nu\pi}{2} - \frac{\pi}{4} α = kπ + 2 ν π − 4 π
홀수 차수 도함수 (J ν ( 2 m + 1 ) J_\nu^{(2m+1)} J ν ( 2 m + 1 ) 의 영점):
j ν , k ( 2 m + 1 ) = β − α 1 , 2 m + 1 α 0 , 2 m + 1 ⋅ 1 β + O ( 1 β 3 ) j_{\nu,k}^{(2m+1)} = \beta - \frac{\alpha_{1,2m+1}}{\alpha_{0,2m+1}} \cdot \frac{1}{\beta} + O\left(\frac{1}{\beta^3}\right) j ν , k ( 2 m + 1 ) = β − α 0 , 2 m + 1 α 1 , 2 m + 1 ⋅ β 1 + O ( β 3 1 )
여기서 β = k π + ν π 2 − 3 π 4 \beta = k\pi + \frac{\nu\pi}{2} - \frac{3\pi}{4} β = kπ + 2 ν π − 4 3 π
정리 2.2 : 점근 전개의 오차항에 대해 다음과 같은 한계 추정이 성립한다:
2n차 도함수의 경우, ν ≥ − 2 n + 1 2 \nu \geq -2n + \frac{1}{2} ν ≥ − 2 n + 2 1 , n ∈ N n \in \mathbb{N} n ∈ N 이고 x x x 가 충분히 클 때:
∣ δ 2 n ( ν , x ) ∣ ≤ 4 ( ν + 2 n ) 2 − 1 4 x exp { 4 ( ν + 2 n ) 2 − 1 4 x } |\delta_{2n}(\nu,x)| \leq \frac{4(\nu+2n)^2-1}{4x} \exp\left\{\frac{4(\nu+2n)^2-1}{4x}\right\} ∣ δ 2 n ( ν , x ) ∣ ≤ 4 x 4 ( ν + 2 n ) 2 − 1 exp { 4 x 4 ( ν + 2 n ) 2 − 1 }
ν → ∞ \nu \to \infty ν → ∞ 이고 k k k 가 고정일 때, Olver의 균일 점근 전개 이론과 Airy 함수를 활용한다:
짝수 차수 도함수의 경우:
j ν , k ( 2 n ) = ν − a k 2 1 / 3 ν 1 / 3 + 2 1 / 3 a k 2 10 ν − 1 / 3 + O ( ν − 1 ) j_{\nu,k}^{(2n)} = \nu - \frac{a_k}{2^{1/3}}\nu^{1/3} + \frac{2^{1/3}a_k^2}{10}\nu^{-1/3} + O(\nu^{-1}) j ν , k ( 2 n ) = ν − 2 1/3 a k ν 1/3 + 10 2 1/3 a k 2 ν − 1/3 + O ( ν − 1 )
홀수 차수 도함수의 경우:
j ν , k ( 2 n + 1 ) = ν − a k ′ 2 1 / 3 ν 1 / 3 + 2 1 / 3 ( a k ′ ) 2 10 ν − 1 / 3 + O ( ν − 1 ) j_{\nu,k}^{(2n+1)} = \nu - \frac{a'_k}{2^{1/3}}\nu^{1/3} + \frac{2^{1/3}(a'_k)^2}{10}\nu^{-1/3} + O(\nu^{-1}) j ν , k ( 2 n + 1 ) = ν − 2 1/3 a k ′ ν 1/3 + 10 2 1/3 ( a k ′ ) 2 ν − 1/3 + O ( ν − 1 )
여기서 a k a_k a k 와 a k ′ a'_k a k ′ 는 각각 Airy 함수 Ai ( x ) \text{Ai}(x) Ai ( x ) 와 Ai ′ ( x ) \text{Ai}'(x) Ai ′ ( x ) 의 k번째 음의 영점이다.
통일된 점화 틀 : 모든 차수 도함수 계수의 통일된 점화식을 수립하여 각 차수별 계산을 피한다.개선된 오차 분석 : Hethcote 방법을 수정하여 고전적 결과보다 더 정확한 오차 한계를 얻는다.확장된 적용 범위 : 고전 이론을 특정 저차 도함수에서 임의 차수 도함수로 확장한다.균일 점근 이론의 응용 : Olver의 균일 점근 전개 이론을 고차 도함수 경우에 교묘하게 적용한다.순수 수학 이론 연구이므로 "실험"은 주로 이론 검증 및 수치 검증을 의미한다:
수학적 귀납법 검증 : 점화식의 정확성을 수학적 귀납법으로 증명점근 분석 : 극한 분석을 통해 점근 전개의 유효성 검증오차 한계 검증 : 구체적 계산을 통해 오차 한계의 타당성 검증McMahon의 원래 결과 (J ν ( x ) J_\nu(x) J ν ( x ) 의 영점) Wong과 Lang의 결과 (J ν ′ ′ ( x ) J''_\nu(x) J ν ′′ ( x ) 의 영점) Olver의 균일 점근 전개 이론 McMahon 형 전개의 수립 : 임의 차수 Bessel 함수 도함수 영점의 McMahon 형 점근 전개를 성공적으로 유도하며, 명확한 계수 표현식을 제공한다.오차 한계의 개선 : 기존 결과에 비해 새로운 오차 한계가 더 정확하며, 특히 고차 도함수의 경우에 그러하다.대 매개변수 전개의 완전성 : ν → ∞ \nu \to \infty ν → ∞ 인 경우 주요항과 수정항을 포함한 완전한 점근 급수를 제공한다.3차 및 4차 도함수의 영점에 대해 다음을 얻는다:
j ν , k ( 3 ) = ν − a k 2 1 / 3 ν 1 / 3 + 2 1 / 3 a k 2 10 ν − 1 / 3 + a k 3 700 ν − 1 + ⋯ j_{\nu,k}^{(3)} = \nu - \frac{a_k}{2^{1/3}}\nu^{1/3} + \frac{2^{1/3}a_k^2}{10}\nu^{-1/3} + \frac{a_k^3}{700}\nu^{-1} + \cdots j ν , k ( 3 ) = ν − 2 1/3 a k ν 1/3 + 10 2 1/3 a k 2 ν − 1/3 + 700 a k 3 ν − 1 + ⋯
j ν , k ( 4 ) = ν − a k ′ 2 1 / 3 ν 1 / 3 + 2 1 / 3 ( a k ′ ) 2 10 ν − 1 / 3 + ( a k ′ ) 3 700 ν − 1 + ⋯ j_{\nu,k}^{(4)} = \nu - \frac{a'_k}{2^{1/3}}\nu^{1/3} + \frac{2^{1/3}(a'_k)^2}{10}\nu^{-1/3} + \frac{(a'_k)^3}{700}\nu^{-1} + \cdots j ν , k ( 4 ) = ν − 2 1/3 a k ′ ν 1/3 + 10 2 1/3 ( a k ′ ) 2 ν − 1/3 + 700 ( a k ′ ) 3 ν − 1 + ⋯
McMahon 형 전개에 대해 조건을 만족할 때:
짝수 차수 도함수: ∣ j ν , k + 1 ( 2 n ) − ( k + ν 2 − 1 4 ) π ∣ ≤ 0.26 ( 4 ( ν + 2 n ) 2 − 1 ) π ( k + ν 2 − 0.314 ) |j_{\nu,k+1}^{(2n)} - (k + \frac{\nu}{2} - \frac{1}{4})\pi| \leq \frac{0.26(4(\nu+2n)^2-1)}{\pi(k + \frac{\nu}{2} - 0.314)} ∣ j ν , k + 1 ( 2 n ) − ( k + 2 ν − 4 1 ) π ∣ ≤ π ( k + 2 ν − 0.314 ) 0.26 ( 4 ( ν + 2 n ) 2 − 1 ) 홀수 차수 도함수: ∣ j ν , k ( 2 n + 1 ) − ( k + ν 2 − 1 4 ) π ∣ ≤ 0.26 ( 4 ( ν + 2 n + 1 ) 2 − 1 ) π ( k + ν 2 − 0.314 ) |j_{\nu,k}^{(2n+1)} - (k + \frac{\nu}{2} - \frac{1}{4})\pi| \leq \frac{0.26(4(\nu+2n+1)^2-1)}{\pi(k + \frac{\nu}{2} - 0.314)} ∣ j ν , k ( 2 n + 1 ) − ( k + 2 ν − 4 1 ) π ∣ ≤ π ( k + 2 ν − 0.314 ) 0.26 ( 4 ( ν + 2 n + 1 ) 2 − 1 ) McMahon (1895) : Bessel 함수 영점의 고전적 점근 전개 수립Watson (1944) : 이론 기초를 완성하고 더 엄격한 증명 제공Olver (1954) : 균일 점근 전개 이론 개발Wong & Lang (1990-1991) : 이론을 Bessel 함수의 2차 도함수로 확장Baricz 등 (2018) : n차 도함수 영점의 실성과 양성 증명본 논문은 이러한 발전 맥락에서 이론을 특정 저차 도함수 경우에서 임의 차수 도함수로 완전히 일반화하고, 더 정확한 오차 분석을 제공한다. 이는 해당 분야의 중요한 이론적 완성이다.
완전성 : Bessel 함수 임의 차수 도함수 영점의 완전한 점근 이론 수립정확성 : 기존 결과보다 더 정확한 오차 한계 제공통일성 : 모든 차수 도함수에 적용 가능한 통일된 이론 틀 제시공백 해소 : Bessel 함수 도함수 영점 점근 이론의 마지막 퍼즐 완성방법론 혁신 : 고차 도함수 점근 문제를 다루는 일반적 방법 제공응용 잠재력 : 관련 물리 및 공학 문제에 대한 이론적 도구 제공적용 조건 : 일부 결과는 특정 매개변수 범위 조건을 만족해야 함계산 복잡성 : 고차항의 계산이 상당히 복잡해질 수 있음수치 검증 : 이론 결과는 실용성 확인을 위해 더 많은 수치 검증 필요수정 Bessel 함수로의 확장 : 유사한 방법을 수정 Bessel 함수 및 그 도함수에 적용다른 특수 함수 : Wright 함수, 초 Bessel 함수 등으로 확장수치 알고리즘 : 이론 결과를 기반으로 효율적인 영점 계산 알고리즘 개발이론적 완전성 : 해당 문제에 대한 완전한 이론적 해결책 제공방법론의 엄밀성 : 성숙한 수학 분석 방법 사용, 증명 과정 엄격결과의 일반성 : 특정 문제 해결뿐 아니라 일반적 이론 틀 제공오차 분석 : 명확한 오차 한계 제공으로 결과의 실용성 강화수치 검증 부족 : 순수 이론 연구로서 충분한 수치 실험 검증 부족응용 사례 부족 : 구체적 응용 장면과 사례 제시 부족계산 복잡도 : 매우 고차 도함수의 경우 계산이 비현실적일 수 있음학술적 영향 : 특수 함수 이론의 중요 분야 완성실용적 가치 : Bessel 함수 도함수 영점이 필요한 응용에 이론적 기초 제공방법론적 가치 : 제시된 기술 방법이 다른 유사 문제에 적용 가능수학 물리학 : 파동 방정식, 열전도 등 문제의 해석 해공학 응용 : 신호 처리, 진동 분석 등 분야수치 계산 : 고정밀 영점 계산 알고리즘의 이론적 기초논문은 해당 분야의 핵심 문헌을 인용하며, 다음을 포함한다:
McMahon의 개척적 업적 (1895) Watson의 고전 교과서 (1944) Olver의 점근 이론 (1954, 1962, 1974) Wong & Lang의 확장 업적 (1990, 1991) 최근 관련 연구 (Baricz 등, 2018; Dimitrov & Lun, 2025) 이러한 인용은 연구의 역사적 연속성과 이론적 완전성을 반영한다.