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

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

基本信息

论文ID: 2509.18877
标题: An explicit formula for the Laplace-Beltrami operator on the Stiefel manifold
作者: Petre Birtea, Ioan Caşu, Dan Comănescu (West University of Timişoara)
分类: math.DG (微分几何), math-ph (数学物理), math.MP (数学物理)
发表时间: 2025年10月14日 (arXiv v2)
论文链接: https://arxiv.org/abs/2509.18877

摘要

本文推导了正交Stiefel流形上Laplace-Beltrami算子的显式公式，将其视为配备Frobenius度量的实矩阵欧几里得空间的约束子流形。利用约束流形上Laplace算子的一般框架，作者提供了用环境欧几里得坐标表示的Laplace-Beltrami算子公式。该结果推广了先前已知的情况，恢复了球面和特殊正交群作为特殊实例的公式。

研究背景与动机

问题定义

本文要解决的核心问题是在正交Stiefel流形 $St_p^n = \{U \in M_{n \times p}(\mathbb{R}) | U^T U = I_p\}$ 上建立Laplace-Beltrami算子的显式公式，其中 $n \geq p \geq 1$ 。

研究重要性

理论意义: 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$ 矩阵，定义为： $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}$$$

关键计算步骤

变换矩阵的计算： $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))$