2025-11-10T02:33:53.699950

An explicit formula for the Laplace-Beltrami operator on the Stiefel manifold

Birtea, Casu, Comanescu

We derive an explicit formula for the Laplace-Beltrami operator on the orthogonal Stiefel manifold, viewed as a constraint submanifold of the Euclidean space of real matrices equipped with the Frobenius metric. Using the general framework of Laplace operators on constraint manifolds, we provide the formula for the Laplace-Beltrami operator in terms of the ambient Euclidean coordinates. The result extends previously known cases, recovering the formulas for the sphere and the special orthogonal group as particular instances.

academic

Stiefel多様体上のLaplace-Beltrami作用素の明示公式

基本情報

論文ID: 2509.18877
タイトル: An explicit formula for the Laplace-Beltrami operator on the Stiefel manifold
著者: Petre Birtea, Ioan Caşu, Dan Comănescu（ティミショアラ西大学）
分類: math.DG（微分幾何学）、math-ph（数学物理学）、math.MP（数学物理学）
発表日: 2025年10月14日（arXiv v2）
論文リンク: https://arxiv.org/abs/2509.18877

要約

本論文は、Frobenius計量を備えた実行列ユークリッド空間の制約部分多様体として見なされる直交Stiefel多様体上のLaplace-Beltrami作用素の明示公式を導出している。制約多様体上のLaplace作用素の一般的枠組みを利用して、著者は環境ユークリッド座標で表現されたLaplace-Beltrami作用素の公式を提供している。この結果は先行する既知の場合を一般化し、球面と特殊直交群の公式を特殊な実例として復元している。

研究背景と動機

問題定義

本論文が解決しようとする中心的問題は、直交Stiefel多様体 $St_p^n = \{U \in M_{n \times p}(\mathbb{R}) | U^T U = I_p\}$ 上（ただし $n \geq p \geq 1$ ）でLaplace-Beltrami作用素の明示公式を確立することである。

研究の重要性

理論的意義: Laplace-Beltrami作用素はリーマン幾何学の中心的対象であり、微分幾何学、数学物理学、確率論に広く応用されている
応用価値: Stiefel多様体は最適化理論、機械学習、信号処理などの分野で頻繁に現れ、特に直交制約を伴う最適化問題に関連している
統一的枠組み: 本研究は異なる種類の制約多様体上のLaplace作用素を扱うための統一的方法を提供している

既存方法の限界

先行研究は主に特殊な場合（球面 $p=1$ または特殊直交群 $p=n$ など）に集中していた
一般的なStiefel多様体上のLaplace-Beltrami作用素の明示表現が欠けていた
既存の射影方法は存在するが、直接的な計算公式が不足していた

核心的貢献

主要な理論的結果: 直交Stiefel多様体上のLaplace-Beltrami作用素の明示公式を確立（定理1.3）
既知結果の統一: 球面と特殊直交群の公式を特殊な場合として包含
計算的枠組み: 環境座標に基づく計算可能な表現を提供
幾何学的洞察: Stiefel多様体の接空間構造と計量性質の深い分析

方法の詳細

タスク定義

滑らかな関数 $\tilde{f}: St_p^n \to \mathbb{R}$ が与えられたとき、Stiefel多様体上のそのLaplace-Beltrami作用素 $\Delta_{St_p^n} \tilde{f}(U)$ （ただし $U \in St_p^n$ ）を計算することが目標である。

理論的枠組み

制約多様体上のLaplace作用素

著者は制約多様体上のLaplace作用素の一般理論を採用している。制約多様体 $S_c := F^{-1}(c)$ （ただし $F = (F_1, \ldots, F_k): M \to \mathbb{R}^k$ は制約関数）に対して、Laplace-Beltrami作用素は以下の公式で与えられる：

$\Delta_{S_c} \tilde{f} = \text{tr}\left((T(T^t T)^{-1} T^t)[\text{Hess } f]\right) - \sum_{\alpha=1}^k \sigma_\alpha \text{tr}\left((T(T^t T)^{-1} T^t)[\text{Hess } F_\alpha]\right)$

Stiefel多様体の制約表現

Stiefel多様体は以下の制約関数により定義される：

$F_{aa}(U) = \frac{1}{2}\|u_a\|^2$ 、 $1 \leq a \leq p$ （列ベクトルの長さ制約）
$F_{bc}(U) = \langle u_b, u_c \rangle$ 、 $1 \leq b < c \leq p$ （列ベクトルの直交性制約）

ここで $u_1, \ldots, u_p$ は行列 $U$ の列ベクトルである。

接空間分析

接空間の明示的記述

$U \in St_p^n$ に対して、接空間は以下の形式を持つ： $T_U St_p^n = \{UA + (I_n - UU^t)C | A \in M_{p \times p}(\mathbb{R}), A = -A^t, C \in M_{n \times p}(\mathbb{R})\}$

基の構成

著者は接空間の正規直交基 $B_U = B'_U \cup B''_U$ を構成した：

反対称部分 $B'_U$ ： $\Delta'_{ab}(U) = UA_{ab}, \quad 1 \leq a < b \leq p$ ここで $A_{ab}$ は標準的な反対称行列基である
射影部分 $B''_U$ ： $\Delta''_{ic}(U) = (I_n - UU^t)C_{ic}, \quad i \in \{p+1, \ldots, n\}, c \in \{1, \ldots, p\}$

主要な結果

定理1.3（主要結果）

滑らかな関数 $\tilde{f}: St_p^n \to \mathbb{R}$ に対して、Laplace-Beltrami作用素は以下の明示公式を持つ：

$\Delta_{St_p^n} \tilde{f}(U) = \Delta f(U) - \left(\frac{n-p+1}{2}\right)\text{tr}(U^t \nabla f(U)) - \frac{1}{2}\text{tr}\left((I_p \otimes (UU^t) + \Lambda(U))\text{Hess } f(U)\right)$

ここで：

$f: M_{n \times p}(\mathbb{R}) \to \mathbb{R}$ は $\tilde{f}$ の滑らかな拡張
$\Lambda(U)$ は $np \times np$ 行列で、以下のように定義される： $\Lambda(U) = \begin{bmatrix} u_1 u_1^t & \cdots & u_p u_1^t \\ \vdots & \ddots & \vdots \\ u_1 u_p^t & \cdots & u_p u_p^t \end{bmatrix}$

主要な計算ステップ

変換行列の計算： $T^t T = \begin{bmatrix} 2I_{\frac{p(p-1)}{2}} & O \\ O & I_p \otimes (I_{n-p} - U_2 U_2^t) \end{bmatrix}$
射影作用素の導出： $T(T^t T)^{-1} T^t = I_{np} - \frac{1}{2}I_p \otimes (UU^t) - \frac{1}{2}\Lambda(U)$
Lagrange乗数関数： $\Sigma(U) = \frac{1}{2}(\nabla f(U)^t U + U^t \nabla f(U))$

実験設定

理論的検証

本論文は主に理論的研究であり、検証は以下の方法で行われている：

特殊な場合の検証：
- $p=1$ のとき、球面上のLaplace-Beltrami作用素公式を復元
- $p=n$ のとき、特殊直交群上の公式を復元
既知結果との比較：
- 文献11の二パラメータ族リーマン計量の場合との等価性を検証
- 射影方法から得られた結果との一致を確認

計算複雑度分析

主要な計算は $np \times np$ 行列演算を含む
Hessian行列の計算複雑度は $O(n^2 p^2)$
全体的なアルゴリズム複雑度は $O(n^2 p^2)$

実験結果

理論的検証結果

球面の場合（ $p=1$ ）

$p=1$ のとき、Stiefel多様体は単位球面 $S^{n-1}$ に退化し、公式は既知の球面Laplace作用素に簡約される： $\Delta_{S^{n-1}} \tilde{f}(u) = \Delta f(u) - \frac{n-1}{2}\langle u, \nabla f(u) \rangle - \frac{1}{2}\text{tr}((uu^t)\text{Hess } f(u))$