[论文解读] Eigenfunctions of the Multidimensional Linear Noise Fokker-Planck Operator via Ladder Operators

Todd K. Leen, Robert Friel+1 more|arXiv (Cornell University)|Aug 31, 2016

Statistical Mechanics and Entropy参考文献 7被引用 5

一句话总结

本文为具有线性漂移和恒定扩散的多维福克-普朗克算子推导出升降算符，实现了本征函数与本征值的显式构造。关键贡献在于证明：在特定坐标系下，这些本征函数可简化为埃尔米特函数的乘积之和，构成求解时变且非齐次福克-普朗克方程的双正交基。

ABSTRACT

The eigenfunctions and eigenvalues of the Fokker-Planck operator with linear drift and constant diffusion are required for expanding time-dependent solutions and for evaluating our recent perturbation expansion for probability densities governed by a nonlinear master equation. Although well-known in one dimension, for multiple dimensions the eigenfunctions are not explicitly given in the literature. We develop raising and lowering operators for the Fokker-Planck (FP) operator and its adjoint, and use them to obtain expressions for the corresponding eigenvalues and eigenfunctions. We show that the eigenfunctions for the forward and adjoint FP operators form a bi-orthogonal set, and that the eigenfunctions reduce to sums of products of Hermite functions in a particular coordinate system.

研究动机与目标

• 为多维前向与伴随福克-普朗克算子提供本征函数与本征值的显式、系统化构造。
• 解决多维情况下本征函数表达式缺失的问题，特别是当漂移矩阵与扩散矩阵不能同时对角化时。
• 为前向与伴随算子建立一个双正交本征函数基，以实现时变解与非线性主方程摄动解的展开。
• 证明在静止态协方差为球对称且方差为1/2的坐标系下，本征函数可简化为埃尔米特函数的乘积之和。

提出的方法

• 引入升降算符 $ V_I = -e_I abla $，其中 $ e_I $ 为漂移矩阵 $ A $ 的右特征向量，从静止态 $ f_0 $ 出发生成前向福克-普朗克算子 $ L $ 的本征函数。
• 利用 $ A $ 的左特征向量 $ w_I $ 定义伴随升降算符 $ ar{V}_I $，从常数函数 $ g_0 = 1 $ 出发生成伴随算子 $ L^lat $ 的本征函数。
• 推导对易关系 $[L, V_I] = \bar{\nu}_I V_I$ 与 $[L^\flat, \bar{V}_I] = \nu_I \bar{V}_I $，确认 $ V_I $ 与 $ \bar{V}_I $ 分别为 $ L $ 与 $ L^\flat $ 的升降算符。
• 将本征函数构造为 $ f_{n_1,\nu_2,\nu_3,\nu_4} = \big(\bigotimes_I (V_I)^{n_I}\big) f_0 $，对应本征值 $ \nu = \nu_I n_I $，类似地构造 $ g_K $，形成完备的双正交集。
• 变换至静止态协方差矩阵 $ \tilde{\nu} $ 与单位阵成比例的坐标系，表明在此框架下，本征函数可表示为埃尔米特函数乘积乘以 $ f_0 $。
• 利用双正交性与归一化关系 $ \braket{g_K | f_K} = \nu_K $，通过本征函数展开求解非齐次福克-普朗克方程。

实验结果

研究问题

• RQ1 当漂移矩阵与扩散矩阵不能同时对角化时，如何系统构造具有线性漂移与恒定扩散的多维福克-普朗克算子的本征函数与本征值？
• RQ2 前向与伴随福克-普朗克算子在多维空间中的本征函数结构如何？其与埃尔米特函数有何关联？
• RQ3 能否为前向与伴随福克-普朗克算子分别定义升降算符，以生成完备的双正交本征基？
• RQ4 在何种坐标变换下，福克-普朗克算子的本征函数可简化为一维埃尔米特函数的乘积？
• RQ5 如何利用本征函数基求解时变且非齐次的福克-普朗克方程？

主要发现

• 前向福克-普朗克算子 $ L $ 的本征函数构造为 $ f_{n_1,\nu_2,\nu_3,\nu_4} = \big(\bigotimes_I (V_I)^{n_I}\big) f_0 $，对应本征值 $ \nu = \nu_I n_I $，其中 $ V_I $ 为由漂移矩阵 $ A $ 的右特征向量导出的升降算符。
• 伴随算子 $ L^\flat $ 的本征函数通过 $ g_{n_1,\nu_2,\nu_3,\nu_4} = \big(\bigotimes_I (\bar{V}_I)^{n_I}\big) g_0 $ 生成，对应本征值 $ \nu = \nu_I n_I $，使用 $ A $ 的左特征向量。
• 本征函数 $ f_K $ 与 $ g_K $ 构成双正交集，满足 $ \braket{g_K | f_K} = \nu_K $，且归一化满足 $ \braket{g_0 | f_0} = 1 $。
• 在静止态协方差矩阵 $ \tilde{\nu} $ 与单位阵成比例（即 $ \tilde{\nu} = \frac{1}{2} I $）的坐标系下，本征函数简化为埃尔米特函数乘积之和乘以 $ f_0 $，推广了一维埃尔米特函数基。
• 福克-普朗克算子 $ L $ 可表示为 $ L = \frac{1}{2} \nu_I V_I \bar{V}_I^\flat $，且 $ L^\flat = \frac{1}{2} \nu_I^* \bar{V}_I V_I^\flat $，类似于量子谐振子哈密顿量，但使用非厄米算符。
• 非齐次福克-普朗克方程（如 $ L P^{(i)} = q^{(i)} $）的解为 $ P^{(i)}(x) = \frac{1}{\nu_K \braket{g_K | f_K}} \braket{g_K | q^{(i)}} f_K(x) $，对 $ K \neq 0 $ 求和，从而实现非线性主方程的微扰求解。