10 окт. 2026 г.

Schrödinger: 100 Years Later

Schrodinger: 100 Years Later

How wave mechanics meets an extremal problem in approximation theory

In 1926, Erwin Schrodinger proposed understanding motion through the possible shapes of a wave. A hundred years later, the differential operator of the spherical rotor appears in a study of the sharp constant relating the norms of functions with a restricted spectrum. The physicist asks which states a system permits. The mathematician asks which shape produces the best possible inequality. Sometimes both questions lead to the same spectral family.

This note brings together Schrodinger's second paper, Quantisation as an Eigenvalue Problem, and the work of Andriy Bondarenko, Joaquim Ortega-Cerda, Danylo Radchenko and Kristian Seip on the Hormander-Bernhardsson extremal function. The modern preprint first appeared on 7 April 2025; the version considered here is dated 23 January 2026 [2]. Schrodinger's paper ends with a different date: the manuscript was received by the journal on 23 February 1926 [1].

The comparison requires precision. In Schrodinger's paper we find the free spherical rotor. In the modern problem we encounter a complex perturbation of that rotor. The full perturbed operator has already been studied in physics; what interests us here is its emergence from extremality in analysis.

1. A paper that explains a way of seeing

The second paper occupies pages 489-527 of volume 79 of Annalen der Physik. It both sets out a programme for a new mechanics and tests that programme on concrete examples. It has three main parts: the relationship between mechanics and optics; the passage from a geometrical description to a wave description; and applications to the oscillator and rotors [1].

The opening takes the reader back to Hamilton. The trajectories of a mechanical system and the rays of geometrical optics obey related variational principles. Yet a ray is an approximate description of light. When diffraction matters, we need a wave, with its amplitude, phase and spatial form. Schrodinger carries this question into mechanics: might a trajectory likewise be an approximation to a fuller wave description?

He does not conceal the fact that the familiar analogy does not yet determine a unique wave equation. In particular, choosing a second-order equation is motivated by simplicity and then tested on examples. This is an important feature of the paper: a powerful idea comes with a clear account of the choice being made.

The formulation of the problem itself begins to change. Instead of searching for an individual path, one must study the collection of admissible wave forms. The space on which those forms are defined becomes part of the physical content.

2. The space of possible configurations

Schrodinger works in configuration space: a point in that space specifies the complete configuration of the system. For a single body, this may be ordinary three-dimensional space. For several bodies, we need the coordinates of all of them. When constraints are present, the configuration space may become a surface or another geometrical object.

Consider two different atoms joined by a bond of fixed length, and for the moment disregard the motion of the molecule as a whole. Only the direction of the bond changes. Its possible directions form the sphere \(S^2\). This sphere describes the molecule's orientations; no material spherical surface is needed.

Kinetic energy determines the geometry of configuration space. The wave equation therefore contains the Laplacian associated with that geometry. This is one of the paper's main lessons: the differential operator expresses the structure of possible motion.

In modern notation, the free rotor has Hamiltonian

\[ \mathsf H_0=\frac{\hbar^2}{2I}(-\Delta_{S^2}), \]

where \(I\) is the moment of inertia and \(\hbar\) is Planck's constant divided by \(2\pi\). The prefactor sets the energy scale; the Laplacian describes how the wave form changes between nearby orientations.

3. From motion to form: three central examples

The wave packet. In the second part of the paper, Schrodinger discusses how a superposition of waves forms a localised group. At suitable scales, its motion approximates that of a classical system. He distinguishes phase velocity from group velocity and acknowledges the role of de Broglie's work. This account of the passage to a classical description should not be read as a claim that an arbitrary wave packet always retains its shape.

The oscillator. In the third part, Hermite functions appear, along with energy levels now written as \(E_k=\hbar\omega(k+\tfrac12)\). The ground-state energy is nonzero. Schrodinger discusses the shapes of the eigenfunctions, their nodes, and the interlacing of zeros of successive polynomials. This part is especially familiar to a reader of approximation theory: well-known orthogonal polynomials acquire a physical meaning.

The rotor. This example makes particularly clear why the number of available directions cannot be reduced arbitrarily. With a fixed axis, the configuration space is a circle and the energy levels are proportional to \(k^2\). With a free axis, it is a sphere and the levels are proportional to \(k(k+1)\). Even when an individual classical motion looks simple, the wave problem must take account of the entire space of permitted configurations.

Later, Schrodinger combines rotation and vibration in a diatomic molecule. At the end of the paper, he turns to perturbation theory and notes agreement between the first-order Stark effect and Epstein's formula. Together, these examples demonstrate that the new programme works, while preserving the distinction between exact solutions and approximations [1].

4. The crucial passage: the free rotor and the form of the equation

For our comparison, the relevant passage is Section 3, item 3, "The rigid rotor with a free axis," on pages 520-522 of the original. The main formulas are on page 521, equations (33)-(34'). Here Schrodinger obtains the spherical wave equation and the spectrum \(k(k+1)\) [1].

In modern notation, after dividing by the energy scale,

\[ -\Delta_{S^2}\psi=\lambda\psi, \qquad \lambda=\frac{2IE}{\hbar^2}. \tag{1} \]

Restrict attention to axisymmetric states: they depend on the polar angle \(\theta\), but not on the azimuthal angle. Equation (1) then becomes

\[ -\frac1{\sin\theta}\frac{d}{d\theta} \left(\sin\theta\frac{d\psi}{d\theta}\right) =\lambda\psi. \tag{2} \]

With \(\xi=\cos\theta\) and \(g(\xi)=\psi(\theta)\), we obtain

\[ \boxed{ -\frac{d}{d\xi}\left((1-\xi^2)g'(\xi)\right) =\lambda g(\xi),\qquad -1<\xi<1. } \tag{3} \]

It matters that we retain this divergence form: first a change in \(g\), then the geometrical factor \(1-\xi^2\), then another differentiation. This form explains both integration by parts and the appearance of Legendre polynomials.

The factor vanishes at \(\xi=\pm1\), but these endpoints are the poles of a sphere, not the walls of a container. The state must be regular there; imposing zero values at both endpoints would define a different problem. The regular solutions of (3) are Legendre polynomials, with eigenvalues \(k(k+1)\).

At the end of this discussion, Schrodinger explicitly observes that the level \(k=0\) corresponds to a nonzero constant function on the sphere. Today we would describe the corresponding normalised state as giving a uniform probability distribution over orientations. That is a modern explanation; the probabilistic interpretation of the wave function should not be attributed to this paper.

5. An equation arising from extremality

The modern problem concerns entire functions of exponential type at most \(\pi\) that are integrable on the real axis. Subject to \(f(0)=1\), one seeks to minimise \(\int_{\mathbb R}|f(t)|\,dt\). Equivalently, one seeks the sharp constant \(\mathscr C\) in the inequality \(|f(0)|\leq\mathscr C\|f\|_1\).

This comparison of norms belongs to the Akhiezer-Krein line of extremal problems; the particular extremal function was studied by Hormander and Bernhardsson. In [2], it is factorised as \(\varphi(z)=\Phi(z)\Phi(-z)\), with alternating real zeros assigned to the two factors. Theorem 1.1 gives a second-order equation for \(\Phi\).

We write it in the form

\[ z^2\Phi''+(2z-\alpha)\Phi' +(b^2z^2-\lambda)\Phi=0, \qquad \alpha=\frac1{2\mathscr C},\quad b=\frac\pi2. \tag{4} \]

Here \(\Phi(0)=1\) and \(\lambda=-\alpha\Phi'(0)\); the latter relation follows immediately by setting \(z=0\) in (4). Thus the parameters of the equation are tied to the extremal function.

At first sight, (4) seems far removed from (3): its variable is complex, and the coefficient of the highest derivative is \(z^2\). However, in Section 4.2 the authors pass to a Fourier representation and obtain a perturbed Legendre operator [2].

6. The exact comparison in the Fourier representation

To see how the transformation works, use the representation

\[ \Phi(z)=\int_{-1}^{1}g(\xi)e^{ibz\xi}\,d\xi. \tag{5} \]

Any normalising factor can be absorbed into \(g\). Differentiation with respect to \(z\) becomes multiplication by \(ib\xi\). Multiplication by \(z\) is transferred to differentiation with respect to \(\xi\).

The pair of terms \(z^2\Phi''+2z\Phi'\) is especially significant: together, they transform into a divergence-form operator. For regular \(g\), the boundary term contains \(1-\xi^2\) and vanishes. The resulting equation is

\[ \boxed{ -\frac{d}{d\xi}\left((1-\xi^2)g'(\xi)\right) -ia\xi g(\xi)=\lambda g(\xi), \qquad a=\alpha b. } \tag{6} \]

This is an exact transformation of the equation, with no limiting procedure or asymptotic replacement. The representation and the operator domain are justified in the spectral analysis of [2].

Compare the two forms. Equation (3) contains the free spherical Laplacian. Equation (6) adds a linear imaginary potential. Since \(\xi=\cos\theta\), that potential is proportional to the cosine of the angle between an orientation and a chosen axis. For a real dipole in a real field, the same angular dependence occurs with a real coefficient; an imaginary coefficient takes the model into non-Hermitian quantum mechanics.

The full operator (6) already appears in a 2013 paper by Katherine Jones-Smith and Rudolf Kalveks. Their Hamiltonian \(-\Delta_{S^2}+i\gamma\cos\theta\), restricted to the axisymmetric sector, coincides with (6) when \(\gamma=-a\). The regularity conditions at the poles also agree [3].

The zeros of the extremal function \(\varphi\) and its factor \(\Phi\) must not, however, be identified with the zeros of the wave function \(g\): these functions are related by Fourier transformation and factorisation.

Reflection from north to south reverses the sign of \(\xi\), while complex conjugation reverses the sign of \(i\). Together, the two operations leave the operator invariant: this is the relevant \(PT\) symmetry. It provides a framework for studying regions in which the spectrum of a complex operator is real. Physicists vary the coupling strength; the extremal problem ties the selected value to a sharp constant.

7. The variational meaning: what connects the two problems

At the very beginning of his second paper, Schrodinger emphasises the connection between the wave equation and a variational problem. Here, the "principle of least action" should be understood as the principle of stationary action: the first variation vanishes, although a minimum need not exist.

The free spectral problem arises from stationarity of a quadratic energy subject to a normalisation constraint. For (6), one can write the complex quadratic functional

\[ \mathcal Q[g]=\frac12\int_{-1}^{1} \left((1-\xi^2)(g')^2-ia\xi g^2-\lambda g^2\right)\,d\xi. \tag{7} \]

Integration by parts in the first variation of (7) gives (6). There are no complex conjugates in this expression: the functional is complex, and we are discussing stationarity rather than the minimum of a real physical energy.

The original approximation problem uses a different functional: the integral of the absolute value. For admissible real variations, its first variation is determined by the sign of the extremal function. A signature, orthogonality and nodes at which the sign changes arise naturally here.

The common foundation is the compatibility of a function with all permitted variations. A full explanation of the connection must nevertheless trace how the condition of \(L^1\) extremality becomes a spectral condition for the factor of half the exponential type. Merely writing down the quadratic functional (7) does not accomplish this: it does not yet explain why this particular extremal solution and this particular parameter \(a\) are selected.

This distinction is especially important in the study of Favard's problem on \(S^3\). The original operator acts on a three-dimensional sphere, whereas the auxiliary rotor equation acts on \(S^2\). Coincidence of the auxiliary operators does not replace a justification of the spherical reduction or the sharpness of the constant. It gives us a language in which to seek that justification.

For an accessible introduction to this sphere problem, its geometry, harmonics, Green's function approach and signatures, see Yuriy Kryakin's Favard on the Sphere \(S^3\). Part I: To Hear the Music of the Spheres. That note introduces the question of estimating a function through its Laplacian after the low harmonics have been removed, and outlines the subsequent parts of the discussion.

8. What is worth remembering a hundred years later

Schrodinger's second paper leaves us with several guiding ideas that are particularly close to our concerns. The geometry of configuration space determines the operator. A wave form requires us to consider every available direction. Regularity and global compatibility select the eigenvalues. Orthogonal polynomials describe concrete states. A variational principle connects form with equation. A simple model allows a powerful idea to be tested and its limits to become visible.

For readers wishing to locate these passages in the original: geometry and the optical analogy are in Section 1; wave groups and the passage to the equation are in Section 2; the oscillator is discussed on pages 514-519; the rotors on pages 519-522; and the diatomic molecule and concluding remarks on perturbations on pages 522-527. These are the journal's printed page numbers, not the page numbers of the PDF file.

A hundred years later, the paper retains its appeal through its way of asking a question: first understand the space of possible states, then find the equation and conditions that select compatible forms. Approximation theory reaches a similar question through the demand for the best possible inequality.

These general observations also apply to contemporary problems in approximation theory. The geometry of the domain and the underlying measure determine the operator, the moments and the admissible variations. In sharp inequalities, best \(L^1\) approximation and extremal problems for functions with a restricted spectrum, a signature can certify optimality, while a spectral equation can help describe the extremal form. The analogy becomes useful when we identify the actual transformation, the boundary conditions and the criterion selecting the solution. Those details allow ideas to move between problems and also show where a proposed transfer still requires proof.

The meeting between mathematics and physics is particularly convincing where an exact transformation can be displayed. The extremal problem gives equation (4); its Fourier representation reveals operator (6); and the free part of that operator is already present in Schrodinger's treatment of the rotor. In this way, the history of an idea acquires a verifiable mathematical content.

Sources

[1] E. Schrodinger. Quantisierung als Eigenwertproblem (Zweite Mitteilung) [Quantisation as an Eigenvalue Problem (Second Communication)]. Annalen der Physik, 79 (1926), 489-527; volume 384 in the journal's continuous numbering. DOI: 10.1002/andp.19263840602. Original paper. The equations in this note are written in modern notation.

[2] A. Bondarenko, J. Ortega-Cerda, D. Radchenko, K. Seip. The Hormander-Bernhardsson extremal function. arXiv:2504.05205v2, 23 January 2026; first version, 7 April 2025. Theorem 1.1 and Section 4.2 are central to the comparison. Full text.

[3] K. Jones-Smith, R. Kalveks. Vector Models in PT Quantum Mechanics. arXiv:1304.5692, 2013. See Section 2, equations (5), (15), (17): a particle on a sphere and its axisymmetric sector.

[4] Y. Kryakin. Favard on the Sphere \(S^3\). Part I: To Hear the Music of the Spheres. An introduction to the spherical Favard problem and the accompanying series of notes.

Historical statements, the exact coincidence of operators, and the open question of transferring variational conditions have different statuses in this note. Its purpose is to make the connection understandable and verifiable, while leaving room for further investigation.

=========================================================

Historical note. An important precursor of the spherical model is the work of Carl M. Bender and R. J. Kalveks, Extending PT Symmetry from Heisenberg Algebra to E2 Algebra (2011; preprint 2010). They studied a quantum rotor on the circle with an imaginary cosine potential and explicitly proposed extensions to higher Euclidean algebras. Jones-Smith and Kalveks subsequently developed the spherical model in 2013. In the preface to his dissertation, Kalveks credits Bender with introducing him to PT-symmetric quantum mechanics and supervising his research, while tracing the proposal to investigate Euclidean algebras to discussions with Chris Isham.

======================================================

# Шрёдингер: 100 лет спустя

*Как волновая механика встречается с экстремальной задачей теории приближений*

В 1926 году Эрвин Шрёдингер предложил смотреть на движение через возможные формы волны. Сто лет спустя дифференциальный оператор сферического ротатора появляется в исследовании точной константы между нормами функций с ограниченным спектром. Физик спрашивает, какие состояния допускает система. Математик спрашивает, какая форма обеспечивает наилучшую оценку. Иногда оба вопроса приводят к одному спектральному семейству.

Повод для этой заметки - вторая статья Шрёдингера «Квантование как задача на собственные значения» и работа Андрия Бондаренко, Жоакима Ортеги-Серды, Даниила Радченко и Кристиана Сейпа об экстремальной функции Хёрмандера-Бернхардссона. Современный препринт впервые появился 7 апреля 2025 года; рассматриваемая версия датирована **23 января 2026 года** [2]. В конце статьи Шрёдингера стоит другая дата: рукопись получена редакцией **23 февраля 1926 года** [1].

Сопоставление требует точности. У Шрёдингера мы найдём свободный сферический ротатор. В современной задаче появляется его комплексное возмущение. Полный возмущённый оператор уже изучался в физике; интерес представляет его возникновение из экстремальности в анализе.

## 1. Статья, в которой объясняется способ видеть

Вторая статья занимает страницы 489-527 в 79-м томе *Annalen der Physik*. Это одновременно программа новой механики и проверка программы на конкретных примерах. Она состоит из трёх больших частей: связи механики с оптикой; перехода от геометрического описания к волновому; применения к осциллятору и ротаторам [1].

Начало возвращает читателя к Гамильтону. Траектории механической системы и лучи геометрической оптики подчиняются родственным вариационным принципам. Но луч - приближённое описание света. Когда существенна дифракция, нужна волна со своей амплитудой, фазой и пространственной формой. Шрёдингер переносит этот вопрос в механику: не является ли траектория таким же приближением к более полному волновому описанию?

Он не скрывает, что знакомая аналогия ещё не определяет единственного волнового уравнения. В частности, выбор уравнения второго порядка обосновывается стремлением к простоте и затем проверяется на примерах. Это важная черта работы: сильная идея сопровождается ясным указанием на сделанный выбор.

Здесь начинает меняться сама постановка задачи. Вместо поиска отдельного пути приходится изучать множество допустимых волновых форм. Пространство, на котором они определены, становится частью физического содержания.

## 2. Пространство возможных положений

Шрёдингер работает в конфигурационном пространстве: его точка задаёт полное положение системы. Для одного тела это может быть привычное трёхмерное пространство. Для нескольких тел нужны координаты всех тел. Для системы с ограничениями пространство может стать поверхностью или другой геометрической фигурой.

Возьмём два различных атома, соединённых связью постоянной длины, и пока забудем о перемещении молекулы как целого. Меняется только направление связи. Возможные направления образуют сферу $S^2$. Она описывает ориентации молекулы; поверхность материального шара здесь не требуется.

Кинетическая энергия задаёт геометрию конфигурационного пространства. Поэтому в волновом уравнении возникает лапласиан, соответствующий этой геометрии. Это один из главных уроков статьи: дифференциальный оператор выражает устройство возможного движения.

В современной записи свободный ротатор имеет гамильтониан

\[
\mathsf H_0=\frac{\hbar^2}{2I}(-\Delta_{S^2}),
\]

где $I$ - момент инерции, а $\hbar$ - постоянная Планка, делённая на $2\pi$. Множитель устанавливает масштаб энергии; лапласиан описывает изменение волновой формы между близкими ориентациями.

## 3. От движения к форме: три главных примера

**Волновой пакет.** Во второй части статьи Шрёдингер обсуждает, как суперпозиция волн образует локализованную группу. При подходящих масштабах её движение приближает движение классической системы. Он различает фазовую и групповую скорости и отмечает роль исследований де Бройля. Это объяснение перехода к классическому описанию не следует воспринимать как утверждение, что произвольный пакет всегда сохраняет свою форму.

**Осциллятор.** В третьей части появляются функции Эрмита и уровни, которые сегодня записывают как $E_k=\hbar\omega(k+\tfrac12)$. Нулевая энергия оказывается ненулевой. Шрёдингер обсуждает форму собственных функций, их узлы и разделение нулей соседних полиномов. Читателю теории приближений эта часть особенно близка: известные ортогональные полиномы приобретают физическое содержание.

**Ротатор.** Здесь особенно ясно видно, почему нельзя произвольно сокращать число доступных направлений. При закреплённой оси остаётся окружность и уровни, пропорциональные $k^2$. При свободной оси возникает сфера и уровни, пропорциональные $k(k+1)$. Даже если отдельное классическое движение выглядит простым, волновая задача должна учитывать всё пространство разрешённых положений.

Позднее Шрёдингер соединяет вращение и колебания двухатомной молекулы. В конце статьи он обращается к теории возмущений и отмечает согласие результата для эффекта Штарка первого порядка с формулой Эпштейна. Вместе эти примеры показывают работоспособность новой программы, сохраняя различие между точным решением и приближением [1].

## 4. Главное место: свободный ротатор и форма уравнения

Для нашего сопоставления нужно открыть **§ 3, пункт 3 «Жёсткий ротатор со свободной осью», страницы 520-522 оригинала**. Основные формулы находятся на странице 521, уравнения (33)-(34′). Здесь Шрёдингер получает сферическое волновое уравнение и спектр $k(k+1)$ [1].

В современных обозначениях после деления на масштаб энергии:

\[
-\Delta_{S^2}\psi=\lambda\psi,
\qquad
\lambda=\frac{2IE}{\hbar^2}.
\tag{1}
\]

Выберем осесимметричные состояния: они зависят от полярного угла $\theta$, но не от азимутального угла. Тогда (1) превращается в

\[
-\frac1{\sin\theta}\frac{d}{d\theta}
\left(\sin\theta\frac{d\psi}{d\theta}\right)
=\lambda\psi.
\tag{2}
\]

После замены $\xi=\cos\theta$, $g(\xi)=\psi(\theta)$ получаем

\[
\boxed{
-\frac{d}{d\xi}\left((1-\xi^2)g'(\xi)\right)
=\lambda g(\xi),\qquad -1<\xi<1.
}
\tag{3}
\]

Важно сохранить именно дивергентную форму: сначала изменение $g$, затем геометрический множитель $1-\xi^2$, затем ещё одно дифференцирование. Она объясняет и интегрирование по частям, и появление полиномов Лежандра.

Множитель исчезает при $\xi=\pm1$, но это полюса сферы, а не стенки сосуда. Требуется регулярность состояния; задавать нулевые значения на обоих концах было бы другой задачей. Регулярные решения (3) - полиномы Лежандра, а собственные значения равны $k(k+1)$.

В конце этого пункта Шрёдингер отдельно отмечает: уровню $k=0$ соответствует постоянная ненулевая функция на сфере. Сегодня мы описали бы такое нормированное состояние как равномерное распределение вероятностей по ориентациям. Это современное объяснение; закон вероятностной интерпретации волновой функции не следует приписывать данной статье.

## 5. Уравнение, пришедшее из экстремальности

Современная задача ставится для целых функций экспоненциального типа не выше $\pi$, интегрируемых на вещественной оси. При условии $f(0)=1$ нужно минимизировать $\int_{\mathbb R}|f(t)|\,dt$. Эквивалентно ищется точная константа $\mathscr C$ в оценке $|f(0)|\leq\mathscr C\|f\|_1$.

Это задача о соотношении норм, принадлежащая линии Ахиезера-Крейна; рассматриваемую экстремальную функцию изучали Хёрмандер и Бернхардссон. В работе [2] она факторизуется как $\varphi(z)=\Phi(z)\Phi(-z)$ с распределением чередующихся вещественных нулей между множителями. Теорема 1.1 даёт для $\Phi$ уравнение второго порядка.

Запишем его в форме

\[
z^2\Phi''+(2z-\alpha)\Phi'
+(b^2z^2-\lambda)\Phi=0,
\qquad
\alpha=\frac1{2\mathscr C},\quad b=\frac\pi2.
\tag{4}
\]

Здесь $\Phi(0)=1$ и $\lambda=-\alpha\Phi'(0)$; последняя связь сразу следует из (4) при $z=0$. Таким образом, параметры уравнения связаны с экстремальной функцией.

На первый взгляд (4) далеко от (3): переменная комплексная, а коэффициент при старшей производной равен $z^2$. Однако в § 4.2 авторы переходят к представлению Фурье и получают возмущённый оператор Лежандра [2].

## 6. Точное сопоставление в представлении Фурье

Чтобы увидеть механику перехода, используем запись

\[
\Phi(z)=\int_{-1}^{1}g(\xi)e^{ibz\xi}\,d\xi.
\tag{5}
\]

Нормировочный множитель можно включить в $g$. Дифференцирование по $z$ превращается в умножение на $ib\xi$. Умножение на $z$ переносится в дифференцирование по $\xi$.

Особенно существенна пара членов $z^2\Phi''+2z\Phi'$: вместе они переходят в дивергентный оператор. При регулярных $g$ граничный член содержит $1-\xi^2$ и исчезает. Получается

\[
\boxed{
-\frac{d}{d\xi}\left((1-\xi^2)g'(\xi)\right)
-ia\xi g(\xi)=\lambda g(\xi),
\qquad a=\alpha b.
}
\tag{6}
\]

Это точное преобразование уравнения, без предельного перехода или асимптотической замены. Представление и область определения оператора обосновываются в спектральном анализе [2].

Сравним формы. В (3) стоит свободный сферический лапласиан. В (6) к нему добавлен линейный мнимый потенциал. При $\xi=\cos\theta$ он пропорционален косинусу угла между ориентацией и выбранной осью. У реального диполя в реальном поле такая угловая зависимость возникает с вещественным коэффициентом; мнимый коэффициент переводит модель в неэрмитову квантовую механику.

Полный оператор (6) уже присутствует в работе Кэтрин Джонс-Смит и Рудольфа Кальвекса 2013 года. Их гамильтониан $-\Delta_{S^2}+i\gamma\cos\theta$ в осесимметричном секторе совпадает с (6) при $\gamma=-a$. Совпадают и условия регулярности в полюсах [3].

При этом нули экстремальной функции $\varphi$ и её множителя $\Phi$ нельзя отождествлять с нулями волновой функции $g$: эти функции связаны преобразованием Фурье и факторизацией.

Отражение север-юг меняет знак $\xi$, а комплексное сопряжение меняет знак $i$. Вместе эти операции сохраняют оператор: это соответствующая $PT$-симметрия. Она позволяет исследовать области вещественного спектра комплексного оператора. Физики варьируют силу воздействия; экстремальная задача связывает выбранное значение с точной константой.

## 7. Вариационный смысл: что объединяет две задачи

В самом начале второй статьи Шрёдингер подчёркивает связь волнового уравнения с вариационной задачей. Принцип «наименьшего действия» при этом следует понимать как принцип стационарного действия: первая вариация исчезает, хотя минимум может отсутствовать.

Свободная спектральная задача получается из стационарности квадратичной энергии при фиксированной нормировке. Для (6) можно записать комплексный квадратичный функционал

\[
\mathcal Q[g]=\frac12\int_{-1}^{1}
\left((1-\xi^2)(g')^2-ia\xi g^2-\lambda g^2\right)\,d\xi.
\tag{7}
\]

Интегрирование по частям в первой вариации (7) даёт (6). В этой записи нет комплексных сопряжений: функционал комплексный, и речь идёт о стационарности, а не о минимуме вещественной физической энергии.

У исходной задачи приближения функционал другой: интеграл абсолютной величины. Для вещественных допустимых вариаций его первая вариация определяется знаком экстремальной функции. Здесь естественно появляются сигнатура, ортогональность и узлы изменения знака.

Общая основа - согласованность функции со всеми разрешёнными изменениями. Но полноценное объяснение связи должно проследить, как условие $L^1$-экстремальности преобразуется в спектральное условие для половинного множителя. Одного наличия квадратичного функционала (7) для этого недостаточно: такой функционал ещё не объясняет выбор именно экстремального решения и параметра $a$.

В исследовании Фавара на $S^3$ это различие особенно важно. Исходный оператор действует на трёхмерной сфере, а вспомогательное уравнение ротатора - на $S^2$. Совпадение вспомогательного оператора не заменяет обоснования сферической редукции и точности константы. Оно даёт язык, в котором можно искать это обоснование.

## 8. Что стоит запомнить через сто лет

Вторая статья Шрёдингера оставляет несколько особенно близких нам ориентиров. Геометрия конфигурационного пространства определяет оператор. Волновая форма требует учитывать все доступные направления. Регулярность и глобальная согласованность выделяют собственные значения. Ортогональные полиномы описывают конкретные состояния. Вариационный принцип связывает форму с уравнением. Простая модель позволяет проверить сильную идею и увидеть границы её применимости.

Для читателя, который хочет найти эти места в оригинале: геометрия и оптическая аналогия находятся в § 1; волновые группы и переход к уравнению - в § 2; осциллятор - на страницах 514-519; ротаторы - на страницах 519-522; двухатомная молекула и заключительные замечания о возмущениях - на страницах 522-527. Здесь указана печатная пагинация журнала, а не номер страницы файла.

Через сто лет привлекательность работы сохраняется в её способе ставить вопрос: сначала понять пространство возможных состояний, затем найти уравнение и условия, которые выделяют согласованные формы. Теория приближений приходит к похожему вопросу через требование наилучшей оценки.

Их встреча становится особенно убедительной там, где можно показать точное преобразование. Из экстремальной задачи приходит уравнение (4); представление Фурье открывает в нём оператор (6); свободная часть этого оператора стоит в разделе о ротаторе у Шрёдингера. Так история идеи получает проверяемое математическое содержание.

### Источники

**[1]** E. Schrödinger. *Quantisierung als Eigenwertproblem (Zweite Mitteilung).* Annalen der Physik, 79 (1926), 489-527; в сквозной нумерации томов - 384. DOI: [10.1002/andp.19263840602](https://doi.org/10.1002/andp.19263840602). [Оригинал статьи](https://jeti.uni-freiburg.de/vorles_Theo_IV/schroedinger_II.pdf). Формулы в заметке приведены в современных обозначениях.

**[2]** A. Bondarenko, J. Ortega-Cerdà, D. Radchenko, K. Seip. *The Hörmander-Bernhardsson extremal function.* [arXiv:2504.05205v2](https://arxiv.org/abs/2504.05205v2), 23 января 2026; первая версия - 7 апреля 2025. Для сопоставления существенны теорема 1.1 и § 4.2. [Полный текст](https://arxiv.org/html/2504.05205v2).

**[3]** K. Jones-Smith, R. Kalveks. *Vector Models in PT Quantum Mechanics.* [arXiv:1304.5692](https://arxiv.org/abs/1304.5692), 2013. См. § 2, уравнения (5), (15), (17): частица на сфере и её осесимметричный сектор.

*Исторические формулировки, точное совпадение операторов и открытый вопрос о переносе вариационных условий имеют в этой заметке разный статус. Задача текста - сделать связь понятной и проверяемой, сохранив пространство для дальнейшего исследования.*
 

9 окт. 2026 г.

S^3 Part I To hear the music of the spheres

 

To Hear the Sphere

Favard on the Sphere \(S^3\). I. Geometry, Oscillations and the Idea

The circle has a familiar language: sines and cosines. Any sufficiently regular oscillation can be decomposed into simple components, and we can ask which of them are slow, which are fast, and which matter for approximation. The sphere has such a language too. It is more complex, but preserves the essential point: geometry determines the admissible oscillations, while a differential operator measures their frequencies.

We want to understand how large a function on the three-dimensional sphere can be when its Laplacian is known and its low harmonics have been removed. The search for a sharp coefficient will lead us to Green's function, approximation in the integral norm, and a special pattern of sign changes. We begin with geometry: the reader needs only an intuitive understanding of distance, averages and oscillation.

How to Imagine \(S^3\)

The familiar surface of a ball is the two-dimensional sphere \(S^2\): one can move along it in two independent directions. It lies in three-dimensional space. Similarly, \(S^3\) is the surface of the unit ball in four-dimensional space. It has three independent tangent directions:

\[ S^3=\{(x_0,x_1,x_2,x_3):x_0^2+x_1^2+x_2^2+x_3^2=1\}. \]

The fourth coordinate helps us define the sphere, while the intrinsic dimension of its surface is three. It is difficult to depict the whole sphere in an ordinary drawing, but the geometry we need can be described through distances and layers.

Choose a north pole. Let \(\theta\) denote the distance from it to another point; on the unit sphere, this distance equals the central angle. It ranges from zero to \(\pi\). Points at the same distance form an ordinary two-dimensional sphere of radius \(\sin\theta\). Near the pole this layer is small, at the equator it is largest, and then it contracts again towards the opposite pole.

A function that has the same value at every point of a layer is called zonal. It can be specified as a function of the single variable \(\theta\). To integrate it, however, we must take the size of the layer into account. For the usual surface measure,

\[ \int_{S^3}A(\theta)\,d\sigma =4\pi\int_0^\pi A(\theta)\sin^2\theta\,d\theta. \]

The factor \(4\pi\) is the area of the unit two-dimensional sphere. The same formula gives the area of \(S^3\), namely \(2\pi^2\). The weight \(\sin^2\theta\) expresses the geometry, which must not be lost when we pass to a single variable.

This will be the decisive difference between dimensions. On \(S^2\), the radial weight is \(\sin\theta\). There, the substitution \(x=\cos\theta\) turns the measure into one with constant density. On \(S^3\), the same substitution leaves the weight \(\sqrt{1-x^2}\). Almost identical geometrical pictures lead to different conditions on the extremal nodes.

Let us fix the normalisations at once: the sphere has unit radius, and we use the usual surface measure. The uniform norm measures the maximum absolute value of a function; the integral norm measures the integral of its absolute value. Replacing surface measure by probability measure does not change the ratio of the norms of two functions, provided the same norm is used for both, but it changes the coefficients of integral kernels. Normalisation therefore matters particularly when we invert the operator.

Harmonics: the Individual Voices of the Sphere

On the circle, a sine or cosine of a single frequency remains the same oscillation after two differentiations: only its multiplier changes. Spherical harmonics play the corresponding role on the sphere. They are grouped by degree \(k\).

A harmonic of degree \(k\) can be defined by taking a homogeneous polynomial in four-dimensional space whose Euclidean Laplacian vanishes, and restricting it to the sphere. For our account, however, another property is more convenient: it is an eigenfunction of the spherical Laplacian.

The Laplace-Beltrami operator measures the combined variation of a function in tangent directions. It takes curvature into account and is independent of coordinates. Write the positive operator as \(L=-\Delta_{S^3}\). For a harmonic of degree \(k\),

\[ LY_k=k(k+2)Y_k. \]

Constants form degree zero. The four coordinate functions have degree one and eigenvalue three. In degree two, the eigenvalue is eight. There are \((k+1)^2\) independent harmonics of degree \(k\): a single frequency on the sphere admits an entire family of shapes.

Different degrees are orthogonal: the integral of the product of two harmonics of different degrees is zero. Together, they provide a complete language for decomposing square-integrable functions. Removing low harmonics is therefore as natural a condition as removing low frequencies from a Fourier series.

Let \(\Pi_{n-1}\) denote the space spanned by all harmonics of degrees less than \(n\). The notation \(f\perp\Pi_{n-1}\) means that the integral of \(fP\) vanishes for every polynomial \(P\) in this space. When \(n=1\), only a zero mean is required. When \(n=2\), the four coordinate moments must also vanish. This condition makes sense for integrable functions as well, since the test polynomials are bounded.

What coefficient is needed to estimate \(\|f\|\) in terms of \(\|Lf\|\) for such functions? The first admissible spectral level is \(n(n+2)\). In the quadratic norm, orthogonality immediately gives the exact answer \(1/[n(n+2)]\). In the uniform and integral norms, different harmonics can combine to create high peaks. For these norms, we must study the shape of the function, as well as the individual eigenvalues.

This is where Favard's problem arises: to find a coefficient of order \(n^{-2}\), determine its sharpness, and describe the extremal behaviour. A single harmonic gives an initial lower bound; a coordinated combination of harmonics may determine the true answer.

Averaging over rotations that fix the pole extracts the zonal part. It preserves the degree of a polynomial and does not increase the integral error. Many of the extremal questions we need therefore become one-dimensional. This is an exact reduction, although the weight itself continues to retain the dimension of the sphere.

From a Sharp Constant to the Right Sign

The history of the method begins with the inversion of differentiation. On the circle, this leads to periodic Bernoulli kernels. Favard's work in the 1930s showed how best approximation of the kernel yields sharp direct estimates for classes of smooth functions. On the sphere, Green's function of the Laplacian plays the same role.

Akhiezer and Krein developed the duality of extremal problems with linear constraints and prescribed moments. In this theory, minimising the integral norm is linked to a bounded function that certifies sharpness. When such a function takes the values \(+1\) and \(-1\), it becomes a signature: the sign of the extremal residual.

The signature has two duties. It must be orthogonal to every admissible approximating polynomial, and it must agree with the sign of the error. Integration against it then gives a lower bound equal to the upper bound already found. The approximation thus receives a certificate of its optimality.

Geronimus representations allow alternating nodes to be described through analytic polynomials. In Babenko's method, the search for such a sign becomes the starting point for difficult extremal problems. This language was also developed in joint studies by Babenko, Deikalova, Gorbachev and Kryakin. Its strength is that it reveals a common construction behind different formulations.

For our problem, an unexpected connection emerged with the comparison of norms: how narrow can the peak of a polynomial of a given degree be? After transforming the even spherical problem, this question selects the same signature as the approximation of Green's function. Different functionals use the same nodes. This coincidence requires proof and rests on the positivity of the interpolation remainder.

The modern description of the Hormander-Bernhardsson extremal function, obtained by Bondarenko, Ortega-Cerda, Radchenko and Seip, supplies an equation and a spectral mechanism for locating the nodes. In our working derivation, this leads to the limiting constant \(F_{1,3}\approx2.85099942\). It concerns best approximation of the kernel and the uniform-norm estimate; the sharpness of the corresponding integral operator norm on the constrained space remains a separate problem.

The full route through the series is as follows:

  1. Chapter II. "Favard: Recovering a Function from Its Oscillations." The circle and \(S^2\), exact signatures and two explicit series.
  2. Chapter III. "The Sphere Puts a Question to Its Nodes." The formulation on \(S^3\), Green's function and a close Gegenbauer candidate.
  3. Chapter IV. "How We Learned to Find the Extremal Form." The history of the norm problem and the ideas leading from moments to an equation.
  4. Chapter V. "A Bridge to Entire Functions." The limiting profile, a common signature and the formula for the spherical constant.
  5. Chapter VI. "When Two Constants Meet on the Sphere." The direct route through the circle, a simple estimate and the remaining questions.

The chapters are self-contained: the necessary notation is repeated. Detailed checks of the ongoing investigation are collected in the proof notes. Three questions are useful while reading: which norm are we estimating, with respect to which measure do we check orthogonality, and what certifies sharpness?

We shall seek a number and reconstruct the form that determines it. On the circle, that form follows a uniform rhythm. On the sphere, the rhythm changes with the measure. To hear it is to understand how geometry chooses the oscillations and their sign.

 

 # Услышать сферу

*Фавар на сфере S³. I. Геометрия, колебания и замысел*

У окружности есть знакомый язык: синусы и косинусы. Любое достаточно регулярное колебание можно разложить на простые составляющие и спросить, какие из них медленные, какие быстрые, какие важны для приближения. У сферы тоже есть такой язык. Он сложнее, но сохраняет главное: геометрия определяет допустимые колебания, а дифференциальный оператор измеряет их частоты.

Мы хотим понять, насколько велика функция на трёхмерной сфере, если известен её лапласиан и удалены низкие гармоники. Поиск точного коэффициента приведёт нас к функции Грина, интегральному приближению и особому рисунку перемен знака. Начнём с геометрии: для чтения достаточно представлений о расстоянии, среднем значении и колебании.

## Как представить S³

Обычная поверхность шара — двумерная сфера \(S^2\): по ней можно двигаться в двух независимых направлениях. Она лежит в трёхмерном пространстве. Аналогично \(S^3\) — поверхность единичного шара в четырёхмерном пространстве. На ней есть три независимых касательных направления:

\[
S^3=\{(x_0,x_1,x_2,x_3):x_0^2+x_1^2+x_2^2+x_3^2=1\}.
\]

Четвёртая координата помогает задать сферу, а внутренняя размерность её поверхности равна трём. Целиком изобразить её обычным рисунком трудно, но нужную нам геометрию можно описать через расстояния и слои.

Выберем северный полюс. Расстояние от него до другой точки обозначим \(\theta\); на единичной сфере оно совпадает с центральным углом. Оно изменяется от нуля до \(\pi\). Точки с одинаковым расстоянием образуют обычную двумерную сферу радиуса \(\sin\theta\). Возле полюса этот слой мал, у экватора — наибольший, затем снова сжимается к противоположному полюсу.

Функция, одинаковая во всех точках слоя, называется зональной. Её достаточно задать как функцию одной переменной \(\theta\). Однако интегрировать её следует с учётом размера слоя. Для обычной поверхностной меры

\[
\int_{S^3}A(\theta)\,d\sigma
=4\pi\int_0^\pi A(\theta)\sin^2\theta\,d\theta.
\]

Множитель \(4\pi\) — площадь единичной двумерной сферы. Из той же формулы получается площадь \(S^3\), равная \(2\pi^2\). Вес \(\sin^2\theta\) выражает геометрию, которую нельзя потерять при переходе к одной переменной.

Это будет решающим различием размерностей. На \(S^2\) радиальный вес равен \(\sin\theta\). Замена \(x=\cos\theta\) там превращает меру в постоянную. На \(S^3\) после такой замены остаётся \(\sqrt{1-x^2}\). Почти одинаковые геометрические картины приводят к разным условиям на экстремальные узлы.

Нормировки зафиксируем сразу: радиус сферы равен единице, мера — обычная поверхностная. Равномерная норма измеряет максимум абсолютной величины функции; интегральная — интеграл её абсолютной величины. Замена площади на вероятностную меру не изменяет отношения этих норм у двух функций, но меняет коэффициенты интегральных ядер. Поэтому нормировка важна именно при обращении оператора.

## Гармоники: отдельные голоса сферы

На окружности синус или косинус одной частоты остаётся тем же колебанием после двукратного дифференцирования: меняется только множитель. На сфере аналогичную роль играют сферические гармоники. Их объединяют в группы по степени \(k\).

Можно определить гармонику степени \(k\), взяв однородный полином в четырёхмерном пространстве, чей евклидов лапласиан равен нулю, и ограничив его на сферу. Но для нашего рассказа удобнее другое свойство: это собственная функция сферического лапласиана.

Лапласиан Лапласа—Бельтрами измеряет суммарное изменение функции по касательным направлениям. Он учитывает кривизну и не зависит от координат. Положительный оператор обозначим \(L=-\Delta_{S^3}\). Для гармоники степени \(k\)

\[
LY_k=k(k+2)Y_k.
\]

Постоянные образуют степень нуль. Четыре координатные функции имеют степень один и собственное значение три. В степени два собственное значение равно восьми. Всего независимых гармоник степени \(k\) имеется \((k+1)^2\): одна частота на сфере допускает целое семейство форм.

Разные степени ортогональны: их попарные интегралы произведений равны нулю. Вместе они дают полный язык разложения квадратично интегрируемых функций. Поэтому удалить низкие гармоники — столь же естественное условие, как удалить низкие частоты из ряда Фурье.

Обозначим через \(\Pi_{n-1}\) все гармоники степеней меньше \(n\). Запись \(f\perp\Pi_{n-1}\) означает, что интеграл \(fP\) равен нулю для каждого полинома \(P\) из этого пространства. При \(n=1\) требуется только нулевое среднее. При \(n=2\) дополнительно исчезают четыре координатных момента. Это условие осмысленно и в интегральной норме, поскольку тестовые полиномы ограничены.

Какой коэффициент нужен для оценки \(\|f\|\) через \(\|Lf\|\) на таких функциях? Первый допустимый уровень равен \(n(n+2)\). В квадратичной норме ортогональность сразу даёт точный ответ \(1/[n(n+2)]\). В равномерной и интегральной нормах разные гармоники могут совместно создавать высокие пики. Для них нужно изучить форму функции, а не только отдельные собственные значения.

Отсюда возникает задача Фавара: найти коэффициент порядка \(n^{-2}\), понять его точность и описать экстремальное поведение. Единичная гармоника даёт начальную нижнюю границу; согласованное сочетание гармоник может определять настоящий ответ.

Усреднение по вращениям вокруг полюса выделяет зональную часть. Оно сохраняет степень полинома и не увеличивает интегральную ошибку. Поэтому многие необходимые нам экстремальные вопросы становятся одномерными. Это точное сокращение задачи, хотя сам вес продолжает хранить размерность сферы.

## От точной константы к точному знаку

История метода начинается с обращения дифференцирования. На окружности оно приводит к периодическим ядрам Бернулли. Работы Фавара 1930-х годов показали, как наилучшее приближение ядра даёт точные прямые оценки для классов гладких функций. На сфере ту же роль играет функция Грина лапласиана.

Ахиезер и Крейн развивали двойственность экстремальных задач с линейными условиями и заданными моментами. В этой теории минимизация интегральной нормы связана с ограниченной функцией, удостоверяющей точность. Когда такая функция принимает значения \(+1\) и \(-1\), она становится сигнатурой — знаком экстремального остатка.

У сигнатуры две обязанности. Она должна быть ортогональна всем разрешённым приближающим полиномам и совпадать со знаком ошибки. Тогда интегрирование против неё даёт нижнюю границу, равную найденной верхней. Так приближение получает свидетельство своей оптимальности.

Представления Геронимуса позволяют описывать чередующиеся узлы через аналитические полиномы. В методе Бабенко поиск такого знака становится отправной точкой для трудных экстремальных задач. Этот язык развивался и в совместных исследованиях Бабенко, Дейкаловой, Горбачева и Крякина. Его достоинство — возможность увидеть общую конструкцию за различными формулировками.

Для нашей задачи неожиданным оказался выход к сравнению норм: насколько узким может быть пик полинома заданной степени? После преобразования чётной сферической задачи этот вопрос выбирает ту же сигнатуру, что и приближение функции Грина. Разные функционалы используют одни узлы. Совпадение требует доказательства и опирается на положительность интерполяционного остатка.

Современное описание экстремальной функции Хёрмандера—Бернхардссона, полученное Бондаренко, Ортегой-Сердой, Радченко и Сейпом, даёт уравнение и спектральный механизм поиска узлов. В нашем рабочем выводе это приводит к предельной константе \(F_{1,3}\approx2{,}85099942\). Она относится к точному приближению ядра и равномерной оценке; точность соответствующей ограниченной интегральной операторной нормы остаётся отдельной задачей.

Полный маршрут цикла таков:

1. **Глава II. «Фавар: восстановить функцию по её колебаниям».** Окружность и \(S^2\), точные сигнатуры и два явных ряда.
2. **Глава III. «Сфера задаёт вопрос своим узлам».** Постановка на \(S^3\), функция Грина и близкий кандидат Гегенбауэра.
3. **Глава IV. «Как научились находить экстремальную форму».** История нормовой задачи и идеи, ведущие от моментов к уравнению.
4. **Глава V. «Мост к целым функциям».** Предельный профиль, общая сигнатура и формула сферической константы.
5. **Глава VI. «Когда две константы встречаются на сфере».** Прямой круговой путь, простая оценка и оставшиеся вопросы.

Главы автономны: необходимые обозначения повторяются. Подробные проверки рабочего исследования собраны в доказательных заметках. Для чтения полезны три вопроса: какую норму мы оцениваем, с какой мерой проверяем ортогональность и чем удостоверена точность.

Мы будем искать число и восстанавливать определяющую его форму. На окружности её задаёт равномерный ритм. На сфере он меняется вместе с мерой. Услышать его — значит понять, как геометрия выбирает колебания и их знак.

S^3 Part II: In the beginning were $S^1$ and $S^2$

# Favard: Recovering a Function from Its Oscillations

*Favard on the Sphere \(S^3\). II. The Circle, the Sphere, and Exact Signatures*

How well can a polynomial of bounded degree approximate a function? Smoothness helps by reconciling neighboring values. But how can this help be measured precisely?

At the core of our story lies Favard’s theorem. Its mechanism inverts differentiation, reducing a problem about an entire class of functions to a single kernel. Once the best approximation of this kernel is understood, direct error bounds become available for every function in the class. On the sphere, differentiation is replaced by the Laplacian, and the Bernoulli kernel is replaced by Green’s function.

We begin with the circle, where the exact nodes are known. Next, we examine \(S^2\), where the geometry unexpectedly preserves the simplicity of the signs. These two cases will prepare us for \(S^3\), where the familiar method encounters a new weight.

## The Circle: A Kernel in Place of a Function Class

On the unit circle, a function is periodic with period \(2\pi\). A trigonometric polynomial of degree less than \(n\) consists of a constant term alongside sines and cosines of frequencies \(1, \ldots, n-1\). The absence of these frequencies means that the function contains no low-order harmonics. We denote this condition by \(f \perp \mathcal{T}_{n-1}\).

For the basic Laplacian \(L = -d^2/d\theta^2\), the Bohr–Favard inequality takes the form

\[
\boxed{\displaystyle
\|f\|_\infty \le \frac{\pi^2}{8n^2} \|Lf\|_\infty,
\qquad f \perp \mathcal{T}_{n-1}.}
\]

The constant is sharp. A pure harmonic of frequency \(n\) would yield only \(1/n^2\). The extremal function utilizes multiple frequencies in concert: its second derivative approaches an alternating step function (a square wave). Smooth approximations to this function show that the constant cannot be improved.

For an arbitrary smooth periodic function, the same constant bounds the approximation error: \(e_{n-1}(f)_\infty \le [\pi^2/(8n^2)] \|f''\|_\infty\). Here, \(e_{n-1}\) denotes the best uniform approximation error by trigonometric polynomials of degree less than \(n\). The lower frequencies are absorbed by the approximating polynomial.

Why does an exact constant emerge? On functions with zero mean, \(L\) is inverted by convolution with a periodic kernel. Under the normalization of convolution by the measure \(d\theta/(2\pi)\), this kernel is given by the Fourier series \(2\sum_{k\ge1} \cos(k\theta)/k^2\). Subtracting a low-degree polynomial from the kernel yields a remainder whose \(L^1\)-norm bounds the reconstruction error. The best remainder has the sign signature \(\operatorname{sign}\cos(n\theta)\), which is orthogonal to lower frequencies.

Its nodes are equally spaced. The sign function itself serves as a dual certificate: integrating against it yields a lower bound that matches the norm of the remainder. This is precisely how sharpness is established, rather than merely an upper bound. The computation yields

\[
F_{1,1} = \frac{4}{\pi} \sum_{j=0}^{\infty} \frac{(-1)^j}{(2j+1)^3} = \frac{\pi^2}{8} \approx 1.2337005501.
\]

The kernel contains the second power of the frequency; the third power appears after integration against the step-like sign function. In classical notation, this is the Favard constant \(K_2\); the index \(1\) in \(F_{1,1}\) indicates the first power of the Laplacian.

## The Sphere: Dimension Enters via the Measure

On the unit sphere \(S^2 \subset \mathbb{R}^3\), a spherical harmonic of degree \(k\) has eigenvalue \(k(k+1)\) with respect to \(L = -\Delta_{S^2}\). Spherical polynomials of degree less than \(n\) form the space \(\Pi_{n-1}\). The condition \(f \perp \Pi_{n-1}\) again signifies the absence of low-frequency components.

The Laplacian combines second derivatives along tangential directions and is coordinate-free. Its inversion kernel depends solely on the geodesic distance \(\theta\) between points. Under the standard surface measure, where the area of the sphere is \(4\pi\), the normalized Green's function is

\[
G_2(\theta) = \frac{-1 - 2\log\sin(\theta/2)}{4\pi}.
\]

It has zero mean. The logarithmic singularity at the pole is integrable—a local footprint of the two-dimensional Laplacian. Just as on the circle, any zero-mean function \(f\) can be recovered via convolution: \(f = G_2 * Lf\).

Let \(\mathcal{E}_n(S^2)\) denote the best \(L^1\)-error in approximating \(G_2\) by spherical polynomials of degree less than \(n\). This yields the bound

\[
\|f\|_p \le \mathcal{E}_n(S^2) \|Lf\|_p,
\qquad f \perp \Pi_{n-1}, \quad 1 \le p \le \infty.
\]

For the uniform norm (\(p = \infty\)), the constant \(\mathcal{E}_n\) is sharp. For the \(L^1\)-norm, this inequality serves as an upper bound; its sharpness on the subspace orthogonal to low-degree harmonics is a separate question. This distinction persists on \(S^3\) as well.

Averaging around the pole allows one to seek the approximating polynomial among zonal polynomials—those depending only on \(\theta\). The radial measure is \(2\pi\sin\theta\,d\theta\). Under the change of variables \(x = \cos\theta\), it becomes \(2\pi\,dx\): the weight drops out. Therefore,

\[
\mathcal{E}_n(S^2) = \frac{1}{2} \inf_{\deg q < n} \int_{-1}^{1} |-\log(1-x) - q(x)|\,dx.
\]

This is an exact algebraic problem on a standard interval. It connects directly with the classical Korkine–Zolotarev theorem: among all monic polynomials of degree \(n\), the minimal \(L^1\)-norm is attained by \(2^{-n}U_n\), where \(U_n\) is the Chebyshev polynomial of the second kind.

Its sign is orthogonal to all polynomials of lower degree. In terms of the angular variable, it takes a particularly simple form:

\[
\operatorname{sign}U_n(\cos\theta) = \operatorname{sign}\sin((n+1)\theta), \qquad \theta_j = \frac{j\pi}{n+1}.
\]

All derivatives of the logarithm are positive. Consequently, the interpolation error at these nodes possesses precisely the required sign, and duality certifies that this yields the best approximation. Note a subtle point: the zonal harmonics on \(S^2\) are the Legendre polynomials, whereas the extremal nodes are governed by a different family, \(U_n\). The spectrum and the signature play distinct roles.

## Two Exact Constants and the Path to the Third

On \(S^2\), not only the asymptotic order of the error is known, but also its exact series expansion:

\[
\boxed{\displaystyle
\mathcal{E}_n(S^2) = 2\sum_{j=0}^{\infty} \frac{1}{(2j+1)\bigl[(2j+1)^2(n+1)^2 - 1\bigr]}.}
\]

From this expression, one readily reads off both the finite problem and its limit. For instance, \(\mathcal{E}_1(S^2) = \log 2\). As the degree grows,

\[
\lim_{n\to\infty} n^2 \mathcal{E}_n(S^2) = F_{1,2} = 2\sum_{j=0}^{\infty} \frac{1}{(2j+1)^3} = \frac{7}{4}\zeta(3) \approx 2.1035995805.
\]

The value \(\zeta(3) = \sum_{k\ge1} k^{-3}\) appears because the sum of the reciprocal cubes of the odd integers equals \(\tfrac{7}{8}\zeta(3)\). On the circle, alternating signs produced a power of \(\pi\); on the two-dimensional sphere, a strictly positive series led to a zeta value. These are two manifestations of the same underlying mechanism: integrating the kernel against the exact signature.

One can also dispense with limits: the denominator of each term is bounded below by \((2j+1)^3 n(n+2)\). Hence, for all \(n \ge 1\),

\[
\mathcal{E}_n(S^2) < \frac{F_{1,2}}{n(n+2)} \le \frac{F_{1,2}}{n^2}.
\]

The denominator \(n(n+2)\) arises here from the exact series and should not be confused with the first admissible eigenvalue \(n(n+1)\). The asymptotic formula demonstrates that in the general \(n^{-2}\) bound for the uniform norm, the constant \(F_{1,2}\) cannot be decreased.

Let us return to the role of Favard's approach in approximation theory. Approximating the kernel and convolving the resulting polynomial with \(Lf\) yields \(e_{n-1}(f)_p \le \mathcal{E}_n(S^2)\|Lf\|_p\), where \(e_{n-1}\) is the error of approximating \(f\) by spherical polynomials.

In this way, a single calculation for Green's function governs an entire class of smooth functions. When smoothness is prescribed in terms of finite differences, one first constructs an averaging operator connecting the differences to the Laplacian, and then applies Favard’s inequality. In the operator-theoretic approach to direct theorems, this is the central link; constructing the averaging operator and analyzing its geometry requires dedicated effort.

On the circle, the same nodes serve all even-order Bernoulli kernels. This allows kernels to be combined while preserving sign compatibility. On the sphere, however, such compatibility requires verifying the remainder for each family of operators individually.

Thus, as we transition to \(S^3\), we have two guiding principles. We know what to look for: the inversion kernel, its exact sign, and its \(L^1\)-approximation error. And we know what cannot be guessed from the spectrum alone: the distribution of the extremal nodes. On \(S^3\), after the substitution \(x = \cos\theta\), a residual weight \(\sqrt{1-x^2}\) remains. This breaks the off-the-shelf Korkine–Zolotarev framework and poses a new problem, which is the subject of the next chapter.

 

# Фавар: восстановить функцию по её колебаниям

*Фавар на сфере S³. II. Окружность, сфера и точные сигнатуры*

Насколько хорошо полином ограниченной степени может описать функцию? Гладкость помогает, согласуя соседние значения. Но как измерить эту помощь точно?

В центре нашего рассказа стоит теорема Фавара. Её механизм обращает дифференцирование и переносит задачу о целом классе функций на одно ядро. Если понять наилучшее приближение этого ядра, становятся доступны прямые оценки для всех функций класса. На сфере производную заменяет лапласиан, а ядро Бернулли — функция Грина.

Начнём с окружности, где точные узлы известны. Затем посмотрим на \(S^2\): там геометрия неожиданно сохраняет простоту знаков. Эти два случая подготовят нас к \(S^3\), где знакомый метод встретится с новым весом.

## Окружность: ядро вместо класса функций

На единичной окружности функция периодична с периодом \(2\pi\). Полином степени меньше \(n\) содержит постоянную составляющую и синусы и косинусы частот \(1,\ldots,n-1\). Отсутствие этих частот означает, что у функции нет низких гармоник. Это условие обозначим \(f\perp\mathcal T_{n-1}\).

Для первого лапласиана \(L=-d^2/d\theta^2\) неравенство Бора—Фавара принимает вид

\[
\boxed{\displaystyle
\|f\|_\infty\le\frac{\pi^2}{8n^2}\|Lf\|_\infty,
\qquad f\perp\mathcal T_{n-1}.}
\]

Коэффициент точен. Разрешённая гармоника частоты \(n\) дала бы только \(1/n^2\). Экстремальная функция использует несколько частот согласованно: её вторая производная приближается к чередующейся ступеньке. Гладкие приближения этой функции показывают, что коэффициент уменьшить нельзя.

Для произвольной гладкой периодической функции тот же коэффициент оценивает ошибку приближения: \(e_{n-1}(f)_\infty\le[\pi^2/(8n^2)]\|f''\|_\infty\). Здесь \(e_{n-1}\) — наименьшая равномерная ошибка среди полиномов степени меньше \(n\). Низкие частоты принимает на себя приближающий полином.

Почему возникает точная константа? На функциях с нулевым средним \(L\) обращается свёрткой с периодическим ядром. При нормировке свёртки мерой \(d\theta/(2\pi)\) оно имеет ряд \(2\sum_{k\ge1}\cos(k\theta)/k^2\). Вычитая из ядра полином низкой степени, получаем остаток, чья интегральная норма оценивает ошибку восстановления. Наилучший остаток имеет знак \(\operatorname{sign}\cos(n\theta)\), ортогональный низким частотам.

Его узлы равномерны. Сам знак служит двойственным свидетельством: интегрирование против него даёт нижнюю границу, совпадающую с нормой остатка. Именно так доказывается точность, а не только верхняя оценка. Вычисление даёт

\[
F_{1,1}=\frac4\pi\sum_{j=0}^{\infty}
\frac{(-1)^j}{(2j+1)^3}=\frac{\pi^2}{8}
\approx1{,}2337005501.
\]

В ядре стоит вторая степень частоты; третья появляется после интегрирования против ступенчатого знака. В классической нумерации это константа \(K_2\); индекс \(1\) в \(F_{1,1}\) означает первую степень лапласиана.

## Сфера: размерность входит через меру

На единичной сфере \(S^2\subset\mathbb R^3\) гармоника степени \(k\) имеет собственное значение \(k(k+1)\) для \(L=-\Delta_{S^2}\). Сферические полиномы степени меньше \(n\) образуют пространство \(\Pi_{n-1}\). Условие \(f\perp\Pi_{n-1}\) снова означает отсутствие низких уровней.

Лапласиан объединяет вторые производные по касательным направлениям и не зависит от выбора координат. Его ядро обращения зависит только от расстояния \(\theta\) между точками. При обычной поверхностной мере, когда площадь сферы равна \(4\pi\), нормированная функция Грина равна

\[
G_2(\theta)=\frac{-1-2\log\sin(\theta/2)}{4\pi}.
\]

Она имеет нулевое среднее. Логарифмическая особенность у полюса интегрируема; это местный след двумерного лапласиана. Как и на окружности, \(f\) с нулевым средним восстанавливается свёрткой \(G_2*Lf\).

Обозначим через \(\mathcal E_n(S^2)\) наименьшую интегральную ошибку приближения \(G_2\) сферическими полиномами степени меньше \(n\). Она даёт оценку

\[
\|f\|_p\le\mathcal E_n(S^2)\|Lf\|_p,
\qquad f\perp\Pi_{n-1},\quad1\le p\le\infty.
\]

Для равномерной нормы коэффициент \(\mathcal E_n\) точен. Для интегральной нормы здесь утверждается верхняя оценка; её точность на подпространстве без низких гармоник — отдельный вопрос. Это различие сохранится и на \(S^3\).

Усреднение вокруг полюса позволяет искать приближающий полином среди зональных — зависящих только от \(\theta\). Радиальная мера равна \(2\pi\sin\theta\,d\theta\). При замене \(x=\cos\theta\) она становится \(2\pi\,dx\): вес исчезает. Поэтому

\[
\mathcal E_n(S^2)=\frac12\inf_{\deg q<n}
\int_{-1}^{1}|-\log(1-x)-q(x)|\,dx.
\]

Это точная алгебраическая задача на обычном отрезке. Она встречается с классическим результатом Коркина—Золотарёва: минимальную интегральную норму среди монических полиномов степени \(n\) имеет \(2^{-n}U_n\), где \(U_n\) — полином Чебышёва второго рода.

Его знак ортогонален всем меньшим степеням. В угловой переменной он особенно прост:

\[
\operatorname{sign}U_n(\cos\theta)
=\operatorname{sign}\sin((n+1)\theta),
\qquad\theta_j=\frac{j\pi}{n+1}.
\]

У логарифма все производные положительны. Поэтому интерполяционный остаток в этих узлах имеет именно нужный знак, и двойственность удостоверяет наилучшее приближение. Заметим тонкость: зональные гармоники \(S^2\) — полиномы Лежандра; экстремальные узлы задаёт другое семейство, \(U_n\). Спектр и сигнатура выполняют разные роли.

## Две точные константы и путь к третьей

На \(S^2\) известен не только порядок ошибки, но и её точный ряд:

\[
\boxed{\displaystyle
\mathcal E_n(S^2)=2\sum_{j=0}^{\infty}
\frac{1}{(2j+1)\bigl[(2j+1)^2(n+1)^2-1\bigr]}.}
\]

Здесь легко увидеть одновременно конечную задачу и её предел. Например, \(\mathcal E_1(S^2)=\log2\). При росте степени

\[
\lim_{n\to\infty}n^2\mathcal E_n(S^2)
=F_{1,2}=2\sum_{j=0}^{\infty}\frac1{(2j+1)^3}
=\frac74\zeta(3)\approx2{,}1035995805.
\]

Функция \(\zeta(3)=\sum_{k\ge1}k^{-3}\) появляется потому, что сумма кубических обратных степеней нечётных чисел равна \(\tfrac78\zeta(3)\). У окружности чередование знаков привело к степени \(\pi\); у двумерной сферы положительный ряд привёл к дзета-значению. Это разные результаты одного механизма: интегрирования ядра против точной сигнатуры.

Можно обойтись и без предельного символа: знаменатель каждого слагаемого не меньше \((2j+1)^3n(n+2)\). Поэтому для всех \(n\ge1\)

\[
\mathcal E_n(S^2)<\frac{F_{1,2}}{n(n+2)}
\le\frac{F_{1,2}}{n^2}.
\]

Знаменатель \(n(n+2)\) здесь возникает из точного ряда. Его не следует путать с первым допустимым собственным значением \(n(n+1)\). Предельная формула показывает, что в общей оценке с \(n^{-2}\), для равномерной нормы, коэффициент \(F_{1,2}\) нельзя уменьшить.

Вернёмся к месту Фавара в теории приближений. Приближая ядро и свёртывая найденный полином с \(Lf\), получаем \(e_{n-1}(f)_p\le\mathcal E_n(S^2)\|Lf\|_p\). Здесь \(e_{n-1}\) — ошибка приближения функции сферическим полиномом.

Так одно вычисление для функции Грина управляет целым классом гладких функций. Если гладкость задана разностями, сначала строят усреднение, связывающее разности с лапласианом, а затем применяют оценку Фавара. В операторном подходе к прямым теоремам это центральное звено; выбор усреднения и его геометрия требуют собственной работы.

На окружности одни узлы служат всем чётным порядкам ядер Бернулли. Это позволяет сочетать ядра, сохраняя согласованность знаков. На сфере такое согласование требует проверки остатка для каждого семейства операторов.

Итак, к переходу на \(S^3\) у нас есть два ориентира. Мы знаем, что искать: ядро обращения, его точный знак и интегральную ошибку. И знаем, чего нельзя угадывать по одному спектру: расположение экстремальных узлов. На \(S^3\) после замены \(x=\cos\theta\) останется вес \(\sqrt{1-x^2}\). Он нарушит готовую схему Коркина—Золотарёва и поставит новую задачу, которой посвящена следующая глава.

---

S^3 III Sphere and nodes

 # The Sphere Poses a Question to Its Nodes

*Favard on the Sphere \(S^3\). III. The Laplacian, Green's Function, and the Search for the Signature*

Favard's theorem begins with a simple observation: if a function possesses no slow oscillations, its magnitude can be controlled in terms of its derivative. On the circle, oscillations are measured by frequencies. On the sphere, their place is taken by spherical harmonics, and the natural differential operator is the Laplacian.

Our goal now is to formulate this problem on the three-dimensional sphere and identify where its exact constant resides. Geometry provides a weight in the integral, the Laplacian supplies the kernel, and the approximation problem demands finding the correct sign-change nodes. These very nodes will turn out to be the central unknown.

## Which Oscillations Do We Exclude?

The unit sphere \(S^3\) lies in four-dimensional space. Its intrinsic dimension is three, and its total volume (surface area) is \(2\pi^2\). We work with the standard Riemannian volume measure \(d\sigma\). The norms \(L^1\) and \(L^\infty\) measure the integral of the absolute value and its supremum, respectively.

We denote the positive Laplacian by \(L = -\Delta_{S^3}\). Spherical harmonics of degree \(k\) are its eigenfunctions, satisfying

\[
L Y_k = k(k+2)Y_k.
\]

Constants correspond to degree zero. Let \(\Pi_{n-1}\) denote the space of all spherical harmonics of degrees strictly less than \(n\). The condition \(f \perp \Pi_{n-1}\) means that \(\int_{S^3} f P\,d\sigma = 0\) for every \(P\) in this space. We exclude not only the constant component, but all low-frequency levels up to degree \(n-1\).

To begin, we consider smooth functions. We seek the smallest constant \(B_{n,p}\) such that

\[
\|f\|_p \le B_{n,p} \|Lf\|_p,
\qquad f \perp \Pi_{n-1}, \quad n \ge 1.
\]

The natural scale is \(n^{-2}\), as the first permitted eigenvalue is \(n(n+2)\). In the \(L^2\)-norm, the answer is immediately furnished by the orthogonality of harmonics: \(B_{n,2} = 1/[n(n+2)]\). However, the integral and uniform norms detect function concentration in different ways. Here, the first eigenvalue alone is no longer sufficient.

We invert the operator. If \(h = Lf\), then \(h\) is likewise orthogonal to low-degree harmonics. The function \(f\) is recovered via the Green's kernel:

\[
f(x) = \int_{S^3} G(d(x,y)) h(y)\,d\sigma(y),
\qquad
G(\theta) = \frac{(\pi - \theta)\cot\theta - 1/2}{4\pi^2}.
\]

Here, \(\theta\) denotes the geodesic distance to the pole. The kernel has zero mean and satisfies \(LG = \delta_{\mathrm{pole}} - 1/(2\pi^2)\). Near the pole, \(G(\theta) \sim 1/(4\pi\theta)\); at the antipode, the singularity is removable. This serves simultaneously as a local inverse of the Laplacian and a global correction for the compactness of the sphere.

Orthogonality provides freedom: one may subtract any low-degree polynomial from the kernel without altering the recovered function. In this way, an inequality for a differential operator transforms into a problem of best approximation.

## When the Sign Certifies Optimality

Let us fix a pole. Averaging over rotations around this pole does not increase the approximation error. It is therefore sufficient to consider polynomials that depend solely on \(x = \cos\theta\). The zonal harmonics are expressed in terms of Chebyshev polynomials of the second kind:

\[
U_k(\cos\theta) = \frac{\sin((k+1)\theta)}{\sin\theta}.
\]

These span all algebraic polynomials of the relevant degrees. The radial measure is \(4\pi\sin^2\theta\,d\theta\), which, after the change of variables, gives rise to the weight \(w(x) = \sqrt{1-x^2}\). Setting \(g_0(x) = G(\arccos x)\), the target approximation error of the kernel takes the form

\[
\boxed{\displaystyle
E_n = 4\pi \inf_{\deg q < n} \int_{-1}^{1} |g_0(x) - q(x)| w(x)\,dx.}
\]

This is an exact one-dimensional problem. It yields \(B_{n,p} \le E_n\) for all \(1 \le p \le \infty\). In the uniform norm (\(p = \infty\)), equality is attained. For the original integral norm (\(p = 1\)), equality requires a separate proof: bounding the norm of the kernel does not by itself establish it.

Why does the sign of the error become decisive? Suppose \(s(x)\) takes values in \(\pm 1\) and is orthogonal to all polynomials of degree less than \(n\) with respect to the weight \(w\). Then, for any admissible \(q\),

\[
\int_{-1}^1 |g_0 - q| w\,dx \ge \int_{-1}^1 (g_0 - q) s w\,dx = \int_{-1}^1 g_0 s w\,dx.
\]

If the sign of the remainder coincides with \(s\), the inequality turns into an equality. The signature is precisely such a certifying sign: it converts a candidate approximation into a provably optimal one.

To obtain it, we minimize \(\int_{-1}^1 |Q| w\,dx\) among monic polynomials of degree \(n\), where the leading coefficient is fixed to one. Variational arguments with respect to the remaining coefficients demonstrate that the sign of the minimizer \(Q_n\) is orthogonal to all lower degrees. It has exactly \(n\) simple roots in the interior \((-1, 1)\); otherwise, a polynomial vanishing at the sign-change points would contradict this orthogonality.

It remains to verify that this sign function suits the Green's function specifically. If \(q_{n-1}\) interpolates \(g_0\) at the roots of \(Q_n\), then

\[
g_0(x) - q_{n-1}(x) = \frac{Q_n(x)}{4\pi^2} \int_0^\infty \frac{\cosh\eta\,d\eta}{(\cosh\eta - x)Q_n(\cosh\eta)}.
\]

For \(-1 < x < 1\), the integral is strictly positive. Consequently, the sign of the remainder is indeed \(\operatorname{sign}Q_n\). This is a crucial stroke of good fortune: the minimal monic polynomial selects the nodes, and the positivity of the remainder transfers them directly into the Favard problem.

Setting \(h(y) = \operatorname{sign}Q_n(\cos d(x,y))\), we obtain an admissible right-hand side of unit norm and the value \(L^{-1}h(x) = E_n\). Smoothing via the heat flow extends this sharpness to the class of smooth functions. Therefore, \(B_{n,\infty} = E_n\).

## Why Gegenbauer Is Almost the Answer

The lowest degrees allow us to verify the construction without heavy machinery. For \(n = 1\), symmetry gives \(Q_1(x) = x\). The signature is positive on one hemisphere and negative on the other. The best constant polynomial is \(G(\pi/2) = -1/(8\pi^2)\); direct integration of the remainder gives \(E_1 = 1/2\).

The next case already reveals the difficulty. For \(n = 2\), we have \(Q_2(x) = x^2 - a^2\). The odd moment vanishes by symmetry; the requirement that the zeroth moment vanish demands that the central interval contain half of the total weight. Thus,

\[
a\sqrt{1-a^2} + \arcsin a = \frac{\pi}{4},
\qquad a \approx 0.4039727533.
\]

A very close candidate is provided by the Gegenbauer polynomials \(C_n^{(2)}\). However, for \(C_2^{(2)}(x) = 12x^2 - 2\), the roots are \(\pm 1/\sqrt{6} \approx \pm 0.4082482905\). Substitution reveals that the zeroth moment of its sign is not zero. The discrepancy is small, but for exact extremality, it is essential.

This distinction has a straightforward explanation: ordinary orthogonality of a polynomial is a linear condition on its values, whereas in our problem, it is its **sign** that must be orthogonal. Replacing values with signs alters the equations governing the nodes. A nearby family of classical orthogonal polynomials is under no obligation to solve the new problem.

Nonetheless, the Gegenbauer candidate remains valuable. Interpolation at its roots yields a similarly positive remainder, providing a constructive upper bound. Even a relatively crude estimate leads to

\[
n(n+2)E_n < 12, \qquad n \ge 1.
\]

Thus, Favard’s inequality is established for all \(n\) with the correct order of magnitude. The lower bound \(B_{n,p} \ge 1/[n(n+2)]\) follows by testing a single harmonic of degree \(n\). Between these bounds lies the question of the exact constant and the limit of \(n^2 E_n\).

For even degrees, another avenue opens. If \(n = 2m\), symmetry allows one to fold the two halves of the radial integral and set \(t = 2\theta\). The weight becomes \(1 - \cos t\). The identical sign conditions arise in minimizing

\[
\int_0^\pi |P(t)|\,dt,
\qquad \deg P \le m, \quad P(0) = 1,
\]

where \(P\) is a cosine polynomial. Indeed, any admissible variation vanishing at zero takes the form \((1 - \cos t)q(\cos t)\). In this way, spherical nodes meet the classical problem of comparing uniform and integral norms.

The question is now sharply posed: determine the signature, characterize its nodes as the degree grows, and evaluate the Green's functional on this sign function.

# Сфера задаёт вопрос своим узлам

*Фавар на сфере S³. III. Лапласиан, функция Грина и поиск сигнатуры*

Теорема Фавара начинается с простого наблюдения: если у функции нет медленных колебаний, её величину можно контролировать через производную. На окружности колебания измеряются частотами. На сфере их место занимают сферические гармоники, а естественным оператором становится лапласиан.

Теперь наша цель — поставить эту задачу на трёхмерной сфере и увидеть, где скрывается её точная константа. Геометрия даст вес в интеграле, лапласиан даст ядро, а задача приближения потребует найти правильные точки перемены знака. Именно эти точки окажутся главным неизвестным.

## Какие колебания мы исключаем

Единичная сфера \(S^3\) лежит в четырёхмерном пространстве. Её собственная размерность равна трём, площадь — \(2\pi^2\). Работаем с обычной поверхностной мерой \(d\sigma\). Нормы \(L^1\) и \(L^\infty\) измеряют соответственно интеграл абсолютной величины и её максимум.

Положительный лапласиан обозначим \(L=-\Delta_{S^3}\). Гармоники степени \(k\) — его собственные функции, причём

\[
LY_k=k(k+2)Y_k.
\]

Постоянные соответствуют степени нуль. Пространство всех гармоник степеней меньше \(n\) обозначим \(\Pi_{n-1}\). Условие \(f\perp\Pi_{n-1}\) означает, что \(\int fP\,d\sigma=0\) для каждого \(P\) из этого пространства. Мы исключаем не только постоянную составляющую, но и следующие низкие уровни — вплоть до степени \(n-1\).

Для начала рассматриваем гладкие функции. Ищем наименьший коэффициент \(B_{n,p}\), при котором

\[
\|f\|_p\le B_{n,p}\|Lf\|_p,
\qquad f\perp\Pi_{n-1},\quad n\ge1.
\]

Естественный масштаб — \(n^{-2}\): первый разрешённый собственный уровень равен \(n(n+2)\). В норме \(L^2\) ответ сразу даёт ортогональность гармоник: \(B_{n,2}=1/[n(n+2)]\). Но интегральная и равномерная нормы иначе чувствуют сосредоточение функции. Здесь одного собственного значения уже недостаточно.

Обратим оператор. Если \(h=Lf\), то \(h\) также ортогональна низким гармоникам. Функция восстанавливается по ядру Грина:

\[
f(x)=\int_{S^3}G(d(x,y))h(y)\,d\sigma(y),
\qquad
G(\theta)=\frac{(\pi-\theta)\cot\theta-1/2}{4\pi^2}.
\]

Здесь \(\theta\) — расстояние до полюса. Ядро имеет нулевое среднее и удовлетворяет \(LG=\delta_{\mathrm{pole}}-1/(2\pi^2)\). Возле полюса \(G(\theta)\sim1/(4\pi\theta)\); в антиподе особенность устранима. Это одновременно локальное обращение лапласиана и глобальная поправка на замкнутость сферы.

Ортогональность даёт свободу: из ядра можно вычитать полином низкой степени, не меняя восстановленную функцию. Так неравенство для дифференциального оператора превращается в задачу наилучшего приближения.

## Когда знак удостоверяет точность

Зафиксируем полюс. Усреднение по вращениям вокруг него не увеличивает ошибку приближения. Поэтому достаточно полиномов, зависящих только от \(x=\cos\theta\). Зональные гармоники выражаются через полиномы Чебышёва второго рода:

\[
U_k(\cos\theta)=\frac{\sin((k+1)\theta)}{\sin\theta}.
\]

Они порождают все алгебраические полиномы нужной степени. Радиальная мера равна \(4\pi\sin^2\theta\,d\theta\), а после замены переменной возникает вес \(w(x)=\sqrt{1-x^2}\). Обозначим \(g_0(x)=G(\arccos x)\). Искомое расстояние ядра имеет вид

\[
\boxed{\displaystyle
E_n=4\pi\inf_{\deg q<n}\int_{-1}^{1}|g_0(x)-q(x)|w(x)\,dx.}
\]

Это точная одномерная задача. Она даёт \(B_{n,p}\le E_n\) для всех \(1\le p\le\infty\). В равномерной норме будет равенство. Для исходной интегральной нормы равенство требует отдельного доказательства: оценка нормы ядра сама по себе его не устанавливает.

Почему решающим становится знак ошибки? Пусть \(s(x)\) принимает значения \(\pm1\) и ортогональна всем полиномам степени меньше \(n\) с весом \(w\). Тогда для любого допустимого \(q\)

\[
\int|g_0-q|w\,dx\ \ge\ \int(g_0-q)sw\,dx
=\int g_0sw\,dx.
\]

Если знак остатка совпал с \(s\), неравенство превращается в равенство. Сигнатура — это такой удостоверяющий знак. Она превращает найденное приближение в доказанно наилучшее.

Чтобы получить её, минимизируем \(\int|Q|w\,dx\) среди монических полиномов степени \(n\): старший коэффициент закреплён и равен единице. Вариации остальных коэффициентов показывают, что знак минимизатора \(Q_n\) ортогонален всем меньшим степеням. У него ровно \(n\) простых внутренних корней. Иначе полином по точкам смены знака нарушил бы эту ортогональность.

Остаётся проверить, что найденный знак подходит именно функции Грина. Если \(q_{n-1}\) интерполирует её в корнях \(Q_n\), то

\[
g_0(x)-q_{n-1}(x)=\frac{Q_n(x)}{4\pi^2}
\int_0^\infty\frac{\cosh\eta\,d\eta}
{(\cosh\eta-x)Q_n(\cosh\eta)}.
\]

При \(-1<x<1\) интеграл положителен. Следовательно, знак остатка действительно равен \(\operatorname{sign}Q_n\). Это важная удача: минимальный монический полином выбирает узлы, а положительный остаток переносит их в задачу Фавара.

Взяв \(h(y)=\operatorname{sign}Q_n(\cos d(x,y))\), получаем допустимую правую часть нормы единица и значение \(L^{-1}h(x)=E_n\). Тепловое сглаживание даёт тот же результат для гладкого класса. Поэтому \(B_{n,\infty}=E_n\).

## Почему Гегенбауэр — почти ответ

Первые степени позволяют проверить конструкцию без большой техники. При \(n=1\) симметрия даёт \(Q_1(x)=x\). Сигнатура положительна на одной полусфере и отрицательна на другой. Наилучший постоянный полином равен \(G(\pi/2)=-1/(8\pi^2)\); прямое интегрирование остатка даёт \(E_1=1/2\).

Следующий случай уже выявляет трудность. При \(n=2\) имеем \(Q_2(x)=x^2-a^2\). Момент нечётной степени исчезает по симметрии; постоянный момент требует, чтобы центральный интервал содержал половину веса. Поэтому

\[
a\sqrt{1-a^2}+\arcsin a=\frac\pi4,
\qquad a\approx0{,}4039727533.
\]

Очень близкий кандидат дают многочлены Гегенбауэра \(C_n^{(2)}\). Но у \(C_2^{(2)}(x)=12x^2-2\) корни равны \(\pm1/\sqrt6\approx\pm0{,}4082482905\). Подстановка показывает: постоянный момент его знака не равен нулю. Расхождение невелико, однако для точной экстремальности оно существенно.

Это различие имеет простую причину. Обычная ортогональность полинома — линейное условие на его значения. В нашей задаче ортогонален его **знак**. Замена значений знаками меняет уравнения для узлов. Близкое семейство ортогональных полиномов не обязано решать новую задачу.

При этом кандидат Гегенбауэра остаётся полезным. Интерполяция в его корнях имеет такой же положительный остаток и даёт конструктивную верхнюю оценку. Уже сравнительно грубый расчёт приводит к

\[
n(n+2)E_n<12,\qquad n\ge1.
\]

Тем самым неравенство Фавара получено для всех \(n\), с правильным порядком. Нижняя оценка \(B_{n,p}\ge1/[n(n+2)]\) следует из одной гармоники степени \(n\). Между этими границами остаётся вопрос о точном коэффициенте и о пределе \(n^2E_n\).

Для чётных степеней появляется ещё одна дверь. Если \(n=2m\), симметрия позволяет сложить две половины радиального интеграла и положить \(t=2\theta\). Вес превращается в \(1-\cos t\). Те же условия на знак возникают при минимизации

\[
\int_0^\pi|P(t)|\,dt,
\qquad \deg P\le m,\quad P(0)=1,
\]

где \(P\) — косинус-полином. Действительно, допустимая вариация, исчезающая в нуле, имеет вид \((1-\cos t)q(\cos t)\). Так сферические узлы встречаются с задачей о сравнении равномерной и интегральной норм.

Теперь вопрос поставлен точно: найти сигнатуру, описать её узлы при росте степени и вычислить функционал Грина на этом знаке.

S^3 Part IV Exstreme form

 # How the Extremal Form Came to Be Found

*Favard on the Sphere \(S^3\). IV. Akhiezer–Krein, Geronimus, Babenko, and the Modern Solution*

An exact constant often conceals an underlying function. Determining a few of its decimal digits is frequently easier than uncovering the structure of the function itself: where its zeros lie, why the sign changes precisely there, and what equation governs the entire landscape.

For the spherical Favard problem, this distinction proved decisive. Once the Laplacian is inverted, the problem reduces to approximating Green's function in the \(L^1\)-norm. Such an approximation demands the correct error sign—the signature. The quest for this signature weaves together several distinct traditions. Below, we trace the essential ideas, maintaining a clear distinction between the general framework and the specific solution.

## The Norm Problem and Bohr’s Question

The uniform norm measures the maximum absolute value of a function; the integral norm sums its absolute values. For a polynomial of bounded degree, a tall peak cannot be arbitrarily narrow. Consequently, an inequality connects these two norms. We are interested in its sharp constant and the polynomial that achieves it.

One can fix the value of a polynomial at a point—for instance, setting \(P(0) = 1\)—and minimize its integral norm. Alternatively, rather than fixing a point value, one can prescribe a coefficient, a derivative, or another linear functional. This yields the general question: how small can the norm be under a prescribed linear constraint?

In their 1936 paper and 1938 monograph, Akhiezer and Krein investigated minimum problems with prescribed initial moments. Their framework relates such a problem to its dual: instead of minimizing the integral norm of a polynomial, one seeks a bounded function with prescribed initial Fourier coefficients. In the extremal case, this dual function is proportional to the sign of the extremal polynomial.

For trigonometric polynomials, this theory, further developed by Geronimus, provides a finite-dimensional computational apparatus: the constant is determined by a determinantal equation. This yields an exact solution in implicit form. However, it does not yet make transparent the asymptotic behavior of the nodes as the degree grows, nor the shape of the limiting function. These questions require deeper analysis.

Another line of inquiry originates with Bohr’s inequality: a function devoid of low frequencies can be bounded in terms of its derivative. This is the very principle at the heart of Favard’s theorem. The absence of slow oscillations allows one to invert differentiation while maintaining control over the norm.

In 1993, Hörmander and Bernhardsson studied a two-dimensional extension of this inequality for the Cauchy–Riemann operator. Their work was also motivated by estimates related to small denominators. Rotational averaging and duality led them to a norm-comparison problem for entire functions of exponential type with bounded spectrum.

They obtained extremely sharp numerical bounds. Under the normalization of exponential type \(\pi\), the constant comparing the point value to the integral norm is approximately \(0.540928822\). Yet identifying the analytic form of the extremal function from these numerical calculations proved elusive. It is therefore important to distinguish the achievements: the numerical value was known to high precision, while the analytic structure of the extremal remained an open problem.

## Babenko: Finding the Signature First

In Babenko’s method, the starting point is a sign function orthogonal to a given subspace. Its sign changes select the nodes. Once found, this signature can be employed across different problems to evaluate their respective functionals.

This strategy emerged in the joint work of Babenko and Kryakin on the \(L^1\)-approximation of step functions. It likewise proved fruitful for comparing the uniform and integral norms of polynomials. The underlying unifying principle is duality: the correct sign certifies extremality.

Babenko’s research note from July 2, 2019, provides a detailed record of this search. Among cosine polynomials of degree at most \(m\), he considered the problem

\[
P(0) = 1, \qquad \int_0^\pi |P(t)|\,dt \longrightarrow \min.
\]

The admissible variations take the form \((1 - \cos t)q(\cos t)\). Consequently, the sign of the minimizer must be orthogonal to lower-degree polynomials with respect to the weight \(1 - \cos t\). It is precisely this weight that later emerges when folding the even-degree problem on \(S^3\).

Geronimus’s representation splits the alternating nodes into two groups. Their products form half-degree factors; an associated analytic polynomial has its roots inside the unit disk. The geometry of real nodes is thus translated into the language of complex phase. This construction was already present in Babenko’s 2019 work and should not be attributed to exposure to the modern solution.

The note sets up the moment equations, works out low-degree cases, and performs extensive computations. For a single sign change, the problem reduces to the familiar equation \(\omega = \cos\omega\); the node is \(\pi/2 + \omega\). Numerical comparisons between approximate models and the target nodes were pushed up to degree 5000.

Yet no simple recurrence relation characteristic of classical orthogonal polynomials emerged; numerical tests indicated a divergence from such standard recurrences. A candidate model for the limiting Fourier transform was proposed as a circular arc, but its extremality remained conjectural. Close agreement in plots and roots was not yet an exact proof.

The role of Gorbachev was substantial. His 2005 paper linked Konyagin’s integral problem (as he referred to the Akhiezer–Krein problem) to Nikolskii-type norm-comparison constants. Babenko’s 2019 note outlines a plan suggested by Gorbachev: construct an interpolation formula using the extremal zeros and deduce a differential equation from the values of the derivatives.

The joint program of Babenko, Deykalova, Gorbachev, and Kryakin is documented in their 2021 lecture on the sharp Jackson–Nikolskii inequality. This represents ongoing work on the method; the 2019 note reflects the stage reached at that time rather than a definitive description of the limiting extremal.

## From Nodes to an Equation

The modern work of Bondarenko, Ortega-Cerdà, Radchenko, and Seip, released as a preprint in 2025 and updated in January 2026, provides an analytic description of the Hörmander–Bernhardsson function. The two lines of inquiry converge on a closely related construction.

The extremal \(\varphi\) has only real zeros, \(\pm\tau_k\). The authors distribute these zeros between the factors \(\Phi(z)\) and \(\Phi(-z)\), so that \(\varphi(z) = \Phi(z)\Phi(-z)\). The moment conditions yield a second-order differential equation for \(\Phi\). An auxiliary functional relation selects the correct solution and governs the zeros.

The computational problem thus becomes a spectral one. In an appropriate orthogonal basis, the resulting matrix is tridiagonal: each coefficient couples only to its immediate neighbors. An eigenvalue alone, however, is not enough; one also requires an eigenvector condition to pick out the extremal branch. This determines both the norm constant and the parameter \(\lambda\), an eigenvalue of an auxiliary operator.

In the scaling used here, the zeros are denoted by \(u_k = \pi\tau_k\). For comparison, one introduces the uniform grid \(v_k = \pi(k + 1/2)\). The parameter \(\lambda\) controls the leading deviation of the true nodes from this grid. The full asymptotic expansion allows expressions involving the nodes to be converted into convergent formulas.

We now introduce the notation used in the spherical chapters. Let \(H\) denote the extremal entire function of exponential type 1, normalized by \(H(0) = 1\). Integrability and this exponential type mean that its spectrum is supported on \([-1, 1]\). We define

\[
\boxed{\displaystyle
B = \frac{1}{2} \int_{\mathbb{R}} |H(u)|\,du \approx 2.903887282752.}
\]

This is half of the minimal \(L^1\)-norm. In the normalization of the modern paper, \(H(u) = \varphi(u/\pi)\) and \(B = \pi/(2\mathscr{C})\), where \(\mathscr{C}\) is the point-evaluation constant. The number \(B\) should not be confused with the spherical Favard constant: both problems share the same underlying shape, but evaluate different quantities.

For our purposes, the primary contribution of the modern solution is the transition from sign constraints to an equation that characterizes the extremal and renders it computable. We deliberately make use of this mechanism. The language of signatures and node splitting, foundational to Babenko’s approach, clarifies its connection to our problem.

This explains why this narrative belongs in the Favard series: the norm problem determines the sign. One must then verify that Green's function conforms to this sign and evaluate the corresponding functional. The next chapter carries out this transition to \(S^3\). The same extremal form will find a new mathematical purpose.

---

*Sources:* N. I. Akhiezer, M. G. Krein, **On Two Minimum Problems Connected with the Problem of Moments** (1936); **Some Questions in the Theory of Moments** (1938). Ya. L. Geronimus, papers on extremal problems (1937–1938). L. Hörmander, B. Bernhardsson, **An extension of Bohr’s inequality** (1993), pp. 179–194. A. G. Babenko, draft **Signum Function**, July 2, 2019, §§3.1–3.6. D. V. Gorbachev, [**Konyagin’s integral problem and Nikolskii’s \((C, L)\)-constants**](https://www.mathnet.ru/rus/timm191) (2005). [Lecture by Babenko and coauthors](https://kma.kmath.ru/Science/Conferences/SBS2021/SBS2021_report.pdf), S. B. Stechkin School, 2021, p. 15 of the report. A. Bondarenko, J. Ortega-Cerdà, D. Radchenko, K. Seip, [**The Hörmander–Bernhardsson extremal function**](https://arxiv.org/abs/2504.05205v2).

 # Как научились находить экстремальную форму

*Фавар на сфере S³. IV. Ахиезер—Крейн, Геронимус, Бабенко и современное решение*

Точная константа часто скрывает за собой функцию. Найти несколько её десятичных знаков бывает легче, чем понять форму этой функции: где находятся нули, почему знак меняется именно там и какое уравнение управляет всей картиной.

Для сферического Фавара это различие оказалось решающим. После обращения лапласиана задача сводится к приближению функции Грина в интегральной норме. А такое приближение требует правильного знака ошибки — сигнатуры. История поиска этого знака соединяет несколько традиций. Проследим идеи, которые нам нужны, сохраняя различие между общей схемой и конкретным решением.

## Нормовая задача и вопрос Бора

Равномерная норма измеряет наибольшую абсолютную величину функции; интегральная суммирует абсолютные величины. У полинома ограниченной степени высокий пик не может быть произвольно узким. Поэтому между этими нормами существует неравенство. Нас интересуют его наилучший коэффициент и полином, на котором он достигается.

Можно зафиксировать значение полинома в точке, например \(P(0)=1\), и минимизировать интегральную норму. Можно вместо значения задать коэффициент, производную или другой линейный функционал. Так возникает общий вопрос: насколько мала норма при заданном линейном условии?

В работе 1936 года и книге 1938 года Ахиезер и Крейн исследовали минимальные задачи с заданными начальными моментами. Их схема связывает такую задачу с двойственной: вместо минимизации интегральной нормы полинома ищется ограниченная функция с предписанными начальными коэффициентами Фурье. В экстремальном случае она пропорциональна знаку полинома.

Для тригонометрических полиномов эта теория, развитая Геронимусом, даёт конечномерный механизм вычисления: константа определяется уравнением для определителя. Это точное решение в неявной форме. Однако оно ещё не делает прозрачными поведение узлов при росте степени и форму предельной функции. Эти вопросы требуют дальнейшего анализа.

Другая линия начинается с неравенства Бора: функцию без низких частот можно оценить через её производную. Это тот же принцип, который лежит в центре теоремы Фавара. Отсутствие медленных колебаний позволяет обратить дифференцирование с контролем нормы.

В 1993 году Хёрмандер и Бернхардссон изучали двумерное расширение этого неравенства для оператора Коши—Римана. Мотивом были также оценки, связанные с малыми знаменателями. Усреднение по вращениям и двойственность привели их к нормовой задаче для целых функций с ограниченным спектром.

Они получили чрезвычайно точные численные границы. В нормировке экспоненциального типа \(\pi\) коэффициент сравнения значения в точке с интегральной нормой равен примерно \(0{,}540928822\). Но распознать экстремальную функцию по вычислениям не удалось. Поэтому важно различать достижения: численное значение уже было известно с большой точностью; аналитическое устройство экстремали оставалось предметом поиска.

## Бабенко: сначала найти сигнатуру

В методе Бабенко отправной объект — знак, ортогональный заданному пространству. Его перемены знака выбирают узлы. После этого разные задачи могут использовать найденную сигнатуру для вычисления своих функционалов.

Такой способ работы проявился в совместных исследованиях Бабенко и Крякина об интегральном приближении ступенчатых функций. Он оказался полезен и для сравнения равномерной и интегральной норм полиномов. Общность состоит в двойственности: правильный знак удостоверяет экстремальность.

Записка Бабенко от 2 июля 2019 года даёт подробный снимок поиска. Среди косинус-полиномов степени не выше \(m\) рассматривалась задача

\[
P(0)=1,\qquad
\int_0^\pi|P(t)|\,dt\longrightarrow\min.
\]

Вариации имеют вид \((1-\cos t)q(\cos t)\). Поэтому знак минимизатора должен быть ортогонален полиномам меньшей степени с весом \(1-\cos t\). Именно этот вес позднее возникает при свёртке чётной задачи на \(S^3\).

Представление Геронимуса разделяет чередующиеся узлы на две группы. Их произведения образуют половинные множители; связанный аналитический полином имеет корни внутри единичного круга. Геометрия действительных узлов переводится в язык комплексной фазы. Эта конструкция уже присутствует у Бабенко в 2019 году и не должна приписываться знакомству с современным решением.

В записке выписаны моментные уравнения, разобраны малые степени и проведены большие вычисления. Для одной точки смены знака задача сводится к знакомому уравнению \(\omega=\cos\omega\); узел равен \(\pi/2+\omega\). Сравнение приближённых моделей с искомыми узлами доведено до степени 5000.

При этом простая рекурсия привычных ортогональных полиномов не обнаружилась. Вычисления показали расхождение с такими рекурсиями. Предлагалась красивая модель предельного преобразования Фурье в виде дуги окружности, но её экстремальность оставалась гипотезой. Близость графиков и корней ещё не давала точного ответа.

Существенна роль Горбачева. Его работа 2005 года связала интегральную задачу Конягина (так он называл задачу Ахиезера-Крейна)  с константами сравнения норм. А в записке 2019 года приведён предложенный им план: построить интерполяционную формулу по экстремальным нулям и через значения производных прийти к дифференциальному уравнению.

Совместное направление Бабенко, Дейкаловой, Горбачева и Крякина документировано докладом 2021 года о точном неравенстве Джексона—Никольского. Перед нами продолжавшаяся работа над методом; записка 2019 года фиксирует достигнутый тогда этап, а не окончательное описание предельной экстремали.

## От узлов к уравнению

Современная работа Бондаренко, Ортеги-Серды, Радченко и Сейпа, опубликованная как препринт в 2025 году и уточнённая в версии января 2026 года, даёт аналитическое описание функции Хёрмандера—Бернхардссона. Две линии исследования встречаются в близкой конструкции.

Экстремаль \(\varphi\) имеет только вещественные нули \(\pm\tau_k\). Авторы распределяют их между множителями \(\Phi(z)\) и \(\Phi(-z)\), так что \(\varphi(z)=\Phi(z)\Phi(-z)\). Моментные условия дают дифференциальное уравнение второго порядка для \(\Phi\). Дополнительное функциональное соотношение выделяет нужное решение и помогает описать нули.

Вычислительная задача становится спектральной. В подходящем ортогональном базисе матрица трёхдиагональна: каждый коэффициент связан только с соседними. Но собственного значения недостаточно. Нужно ещё условие на собственный вектор, выбирающее экстремальную ветвь. Так определяются нормовая константа и параметр \(\lambda\), собственное значение вспомогательного оператора.

В нашем масштабе нули обозначаются \(u_k=\pi\tau_k\). Для сравнения вводятся точки равномерной сетки \(v_k=\pi(k+1/2)\). Параметр \(\lambda\) управляет первым отклонением настоящих узлов от этой сетки. Полная асимптотика позволяет превращать выражения по узлам в сходящиеся формулы.

Теперь введём обозначение, используемое в сферических главах. Пусть \(H\) — экстремальная целая функция типа единица, нормированная условием \(H(0)=1\). Интегрируемость и этот тип означают, что её частоты ограничены отрезком \([-1,1]\). Положим

\[
\boxed{\displaystyle
B=\frac12\int_{\mathbb R}|H(u)|\,du
\approx2{,}903887282752.}
\]

Это половина минимальной интегральной нормы. В нормировке современной статьи \(H(u)=\varphi(u/\pi)\) и \(B=\pi/(2\mathscr C)\), где \(\mathscr C\) — коэффициент оценки значения в точке. Число \(B\) не следует смешивать с константой сферического Фавара: задачи используют общую форму, но вычисляют разные величины.

Для нас главный вклад современного решения — переход от условий на знак к уравнению, которое определяет экстремаль и делает её вычислимой. Мы сознательно используем этот механизм. Язык сигнатур и разделения узлов, существенный в методе Бабенко, позволяет увидеть его связь с нашей задачей.

Так объясняется место этой истории в цикле о Фаваре. Нормовая задача выбирает знак. Затем требуется проверить, что функция Грина подчиняется ему, и вычислить соответствующий функционал. Следующая глава проведёт этот переход к \(S^3\). Одна и та же экстремальная форма получит другое математическое назначение.

---

*Источники.* Н. И. Ахиезер, М. Г. Крейн, **О двух minimum-проблемах, связанных с проблемами моментов** (1936); **О некоторых вопросах теории моментов** (1938). Я. Л. Геронимус, работы об экстремальных задачах (1937–1938). L. Hörmander, B. Bernhardsson, **An extension of Bohr’s inequality** (1993), с. 179–194. А. Г. Бабенко, черновик **Signum Function**, 2 июля 2019 года, §§3.1–3.6. Д. В. Горбачев, [**Интегральная задача Конягина и (C,L)-константы Никольского**](https://www.mathnet.ru/rus/timm191) (2005). [Доклад Бабенко с соавторами](https://kma.kmath.ru/Science/Conferences/SBS2021/SBS2021_report.pdf), школа С. Б. Стечкина, 2021, с. 15 отчёта. A. Bondarenko, J. Ortega-Cerdà, D. Radchenko, K. Seip, [**The Hörmander–Bernhardsson extremal function**](https://arxiv.org/abs/2504.05205v2).