zhenzhubay.com

珍珠湾全球网

 找回密码
 用户注册

tag 标签: 行列式

相关帖子

版块 作者 回复/查看 最后发表

没有相关内容

相关日志

分享 量子物理中的行列式--WICK 定理新证(2)
热度 2 岳东晓 2017-3-24 06:38
好久没有写科技博文了。因为写科技博文其实不是一件容易的事。写完 量子物理的初中代数基础--配平方(1) 》中后,后续博文我一直没写,因为感觉写出来好像是在重复书上的内容。一篇博文至少得有点新东西吧。如果仅仅是重复书上有的东西,那读者还不如去 买本书看,一般来说,书上写得还严谨些。昨天终于想到了点新东西。本博文虽然标题说量子物理,但是内容完全是数学。至于这点数学为什么对物理重要是后续需要讲的。所以,现在不要去管物理,只看数学。 在前一篇里,我们推导了如下公式: $\int_{-\infty}^{+\infty} e^{- x A x } d^{n}x = \sqrt{\frac{\pi^n}{\det A}}$ $\int e^{- x A x + bx } d^{n}x =\exp( \frac{1}{4} b A^{-1} b ) \sqrt{\frac{\pi^n}{\det A}}$ 其中 A 是一个 n x n 的对称矩阵,b 为一个向量,积分上下限为正负无穷。 但在物理中,我们经常需要计算下面这样的积分: $\int x_i x_j e^{- x A x } d^{n}x $ 怎么办呢? 有一个常见的技巧是,将上面第二个积分对 b 取导数: $\int x_i x_j e^{- x A x + bx } d^{n}x = \int \frac{\partial}{\partial b_i} \frac{\partial}{\partial b_j} e^{- x A x + bx } d^{n}x = \sqrt{\frac{\pi^n}{\det A}} \frac{\partial}{\partial b_i} \frac{\partial}{\partial b_j} \exp(\frac{1}{4}b A^{-1} b) $ 因为积分变量 x 与 b 无关,积分-微分顺序可以交换。最后令上面的b=0,就会得到我们要的结果。 写开来看得更清楚一点: \frac{\partial}{\partial b_i} \frac{\partial}{\partial b_j} \exp(\frac{1}{4}b A^{-1} b) |_{b=0} = \frac{\partial}{\partial b_i} \frac{\partial}{\partial b_j} \exp (\frac{1}{4} b_k \ {A^{-1}}_{kl} \ b_l ) |_{b=0} = \frac{1}{2} \ {A^{-1}}_{ij} 上面的步骤目视可见,对b 进行一次微分后,拉下另一个 b。再进行一次微分时,这个拉下的 b 变成 1,剩下的另一半因为含有 b ,在 b=0 时就消失了。1/4 成了 1/2 是因为有两项,一项 $x_i A_{ij} x_j$, 一项 $x_j A_{ji} x_i$。这样,我们得到: \int x_i x_j e^{- x A x } d^{n}x = \frac{1}{2} \sqrt{\frac{\pi^n}{\det A}} \ {A^{-1}}_{ij} 如果要计算下面这样的积分呢? $\int x_i x_j x_k x_l e^{- x A x } d^{n}x $ 同样的方法,只需要对 b 进行四次偏微分最后令 b=0 即可: $\int x_i x_j x_k x_l e^{- x A x + bx } d^{n}x = \int \frac{\partial}{\partial b_i} \frac{\partial}{\partial b_j} \frac{\partial}{\partial b_k} \frac{\partial}{\partial b_l} e^{- x A x + bx } d^{n}x = \sqrt{\frac{\pi^n}{\det A}} \frac{\partial}{\partial b_i} \frac{\partial}{\partial b_j} \frac{\partial}{\partial b_k} \frac{\partial}{\partial b_l}\exp(\frac{1}{4} b A^{-1} b )$ 心算一下就知道,上面的结果里会把i,j,k,l 两两组合,出现 (i,j)(k,l), (i,k)(j,l) 与 (i,l) (j,k) 三个组合,最后是 $ {A^{-1}}_{ij} \ {A^{-1}}_{kl} + {A^{-1}}_{ik} \ {A^{-1}}_{jl} + {A^{-1}}_{il} \ {A^{-1}}_{jk}$ 注意,上面的 ${A^{-1}}_{ij}$ 是 A 的逆矩阵的ij 元素,而不是 $\frac{1}{A_{ij}}$。另外,注意,如果有奇数个 x, 这个积分由于对称性为零。 更多个 x 的情况,类似的两两组合就算出结果了。这叫做 WICK 定理,与其说是定理,还不如说是一个快速计算规则。它看似简单,却大有用场。量子物理里的很多计算都要用到。 上面的推导过程基本上是书上的,也是我们学过的。但读者可能想,加入一个 b*x 项,求导数后再令 b 为零,好像兜了一个圈子,巧妙确实巧妙,但为什么不能更为直接地计算呢。直接将那个高斯积分对 $A_{ij} $ 求导,不也生成出 $x_i x_j$吗?这条路肯定应该是可以的。且让我们来试试吧。 $ \int x_i x_j \exp({- x \ A \ x })\ d^{n}x = - \frac{\partial}{\partial A_{ij}} \int \exp({- x \ A \ x })\ d^{n}x = - \frac{\partial}{\partial A_{ij}} \sqrt{\frac{\pi^n}{\det A}}$ 摆在面前的问题是,我们需要计算一个矩阵的行列式对其一个矩阵元的偏导数。如果以这个矩阵元所在的行展开,行列式中求和中的一项是这个矩阵元乘以其余子式,而余子式矩阵除以其行列式为逆矩阵。 $\frac{\partial}{\partial A_{ij}} \det A = adj(A)_{ji} = (A^{-1})_{ij} \det A$ 利用这个关系,我们完全不需要加入 bx 项,直接就算出来了。 $- \frac{\partial}{\partial A_{ij}} \sqrt{\frac{\pi^n}{\det A}} = - \sqrt{\pi^n} \ (-\frac{1}{2}) (\det A) ^{-\frac{3}{2}} \frac{\partial}{\partial A_{ij}} \det A = \frac{1}{2} \sqrt{\frac{\pi^n}{\det A}}\ (A^{-1})_{ij} $ 这正是上面我们得出的$x_i x_j$ 的结果。可谓殊途同归。但是我们的问题还没完,还得推广到更多变量乘积的情况。 我们有 \int x_i x_j x_k x_l \ \exp({- x \ A \ x })\ d^{n}x = \frac{\partial}{\partial A_{kl}} \left 后面大括号里的我们已经算出来了,为 $ -\sqrt{1/\det A}\ (A^{-1})_{ij}$,现在的任务仅仅 是把这个结果对 $A_{kl}$ 求导数。$\det A $ 求导数我们知道怎么算了,问题是 $ (A^{-1})_{ij}$ 对 $A_{kl}$ 求导数 --- 逆矩阵的一个矩阵元对原矩阵元求偏导,看起来很不简单。但是稍微思考一下,发现有办法。 $A^{-1} A = 1$ 所以, $\partial A^{-1} \ A + A^{-1} \ \partial A =0$ 因此, $\partial A^{-1} = - A^{-1} \ \partial A \ A^{-1} $ 终于看到了曙光! $- \frac{\partial}{\partial A_{kl}} (A^{-1})_{ij} = (A^{-1})_{im} \ \frac{\partial A_{mn} }{\partial A_{kl} }\ (A^{-1} )_{nj} \\ = (A^{-1})_{im} \ \delta_{m,k} \delta_{n,l}\ (A^{-1} )_{nj} + (A^{-1})_{im} \ \delta_{m,l} \delta_{n,k}\ (A^{-1} )_{nj}\\ = (A^{-1})_{ik} \ (A^{-1} )_{lj} + (A^{-1})_{il} \ (A^{-1} )_{kj} $ 上面的负号正好跟前面的负号抵消。有了这个结果以及上面行列式求导的结果,相关的计算就是套下面这两个公式了: $\frac{\partial \det A }{\partial A_{ij}} = (A^{-1})_{ij} \det A$ $- \frac{\partial\ (A^{-1})_{ij}}{\partial A_{kl}} = (A^{-1})_{ik} \ (A^{-1} )_{lj} + (A^{-1})_{il} \ (A^{-1} )_{kj}$ 条条道路通罗马。至此,我们没有引入 bx 项作为生成器就得出了WICK 定理的结果。我个人认为,我上面的这个推导更为直接。
个人分类: 科普|6334 次阅读|1 个评论
分享 量子物理中的行列式--WICK 定理新证(2)
岳东晓 2017-3-21 14:52
好久没有写科技博文了。因为写科技博文其实不是一件容易的事。写完 量子物理的初中代数基础--配平方(1) 》中后,后续博文我一直没写,因为感觉写出来好像是在重复书上的内容。一篇博文至少得有点新东西吧。如果仅仅是重复书上有的东西,那读者还不如去 买本书看,一般来说,书上写得还严谨些。昨天终于想到了点新东西。本博文虽然标题说量子物理,但是内容完全是数学。至于这点数学为什么对物理重要是后续需要讲的。所以,现在不要去管物理,只看数学。 在前一篇里,我们推导了如下公式: $\int_{-\infty}^{+\infty} e^{- x A x } d^{n}x = \sqrt{\frac{\pi^n}{\det A}}$ $\int e^{- x A x + bx } d^{n}x =\exp( \frac{1}{4} b A^{-1} b ) \sqrt{\frac{\pi^n}{\det A}}$ 其中 A 是一个 n x n 的对称矩阵,b 为一个向量,积分上下限为正负无穷。 但在物理中,我们经常需要计算下面这样的积分: $\int x_i x_j e^{- x A x } d^{n}x $ 怎么办呢? 有一个常见的技巧是,将上面第二个积分对 b 取导数: $\int x_i x_j e^{- x A x + bx } d^{n}x = \int \frac{\partial}{\partial b_i} \frac{\partial}{\partial b_j} e^{- x A x + bx } d^{n}x = \sqrt{\frac{\pi^n}{\det A}} \frac{\partial}{\partial b_i} \frac{\partial}{\partial b_j} \exp(\frac{1}{4}b A^{-1} b) $ 因为积分变量 x 与 b 无关,积分-微分顺序可以交换。最后令上面的b=0,就会得到我们要的结果。 写开来看得更清楚一点: \frac{\partial}{\partial b_i} \frac{\partial}{\partial b_j} \exp(\frac{1}{4}b A^{-1} b) |_{b=0} = \frac{\partial}{\partial b_i} \frac{\partial}{\partial b_j} \exp (\frac{1}{4} b_k \ {A^{-1}}_{kl} \ b_l ) |_{b=0} = \frac{1}{2} \ {A^{-1}}_{ij} 上面的步骤目视可见,对b 进行一次微分后,拉下另一个 b。再进行一次微分时,这个拉下的 b 变成 1,剩下的另一半因为含有 b ,在 b=0 时就消失了。1/4 成了 1/2 是因为有两项,一项 $x_i A_{ij} x_j$, 一项 $x_j A_{ji} x_i$。这样,我们得到: \int x_i x_j e^{- x A x } d^{n}x = \frac{1}{2} \sqrt{\frac{\pi^n}{\det A}} \ {A^{-1}}_{ij} 如果要计算下面这样的积分呢? $\int x_i x_j x_k x_l e^{- x A x } d^{n}x $ 同样的方法,只需要对 b 进行四次偏微分最后令 b=0 即可: $\int x_i x_j x_k x_l e^{- x A x + bx } d^{n}x = \int \frac{\partial}{\partial b_i} \frac{\partial}{\partial b_j} \frac{\partial}{\partial b_k} \frac{\partial}{\partial b_l} e^{- x A x + bx } d^{n}x = \sqrt{\frac{\pi^n}{\det A}} \frac{\partial}{\partial b_i} \frac{\partial}{\partial b_j} \frac{\partial}{\partial b_k} \frac{\partial}{\partial b_l}\exp(\frac{1}{4} b A^{-1} b )$ 心算一下就知道,上面的结果里会把i,j,k,l 两两组合,出现 (i,j)(k,l), (i,k)(j,l) 与 (i,l) (j,k) 三个组合,最后是 $ {A^{-1}}_{ij} \ {A^{-1}}_{kl} + {A^{-1}}_{ik} \ {A^{-1}}_{jl} + {A^{-1}}_{il} \ {A^{-1}}_{jk}$ 注意,上面的 ${A^{-1}}_{ij}$ 是 A 的逆矩阵的ij 元素,而不是 $\frac{1}{A_{ij}}$。另外,注意,如果有奇数个 x, 这个积分由于对称性为零。 更多个 x 的情况,类似的两两组合就算出结果了。这叫做 WICK 定理,与其说是定理,还不如说是一个快速计算规则。它看似简单,却大有用场。量子物理里的很多计算都要用到。 上面的推导过程基本上是书上的,也是我们学过的。但读者可能想,加入一个 b*x 项,求导数后再令 b 为零,好像兜了一个圈子,巧妙确实巧妙,但为什么不能更为直接地计算呢。直接将那个高斯积分对 $A_{ij} $ 求导,不也生成出 $x_i x_j$吗?这条路肯定应该是可以的。且让我们来试试吧。 $ \int x_i x_j \exp({- x \ A \ x })\ d^{n}x = - \frac{\partial}{\partial A_{ij}} \int \exp({- x \ A \ x })\ d^{n}x = - \frac{\partial}{\partial A_{ij}} \sqrt{\frac{\pi^n}{\det A}}$ 摆在面前的问题是,我们需要计算一个矩阵的行列式对其一个矩阵元的偏导数。如果以这个矩阵元所在的行展开,行列式中求和中的一项是这个矩阵元乘以其余子式,而余子式矩阵除以其行列式为逆矩阵。 $\frac{\partial}{\partial A_{ij}} \det A = adj(A)_{ji} = (A^{-1})_{ij} \det A$ 利用这个关系,我们完全不需要加入 bx 项,直接就算出来了。 $- \frac{\partial}{\partial A_{ij}} \sqrt{\frac{\pi^n}{\det A}} = - \sqrt{\pi^n} \ (-\frac{1}{2}) (\det A) ^{-\frac{3}{2}} \frac{\partial}{\partial A_{ij}} \det A = \frac{1}{2} \sqrt{\frac{\pi^n}{\det A}}\ (A^{-1})_{ij} $ 这正是上面我们得出的$x_i x_j$ 的结果。可谓殊途同归。但是我们的问题还没完,还得推广到更多变量乘积的情况。 我们有 \int x_i x_j x_k x_l \ \exp({- x \ A \ x })\ d^{n}x = \frac{\partial}{\partial A_{kl}} \left 后面大括号里的我们已经算出来了,为 $ -\sqrt{1/\det A}\ (A^{-1})_{ij}$,现在的任务仅仅 是把这个结果对 $A_{kl}$ 求导数。$\det A $ 求导数我们知道怎么算了,问题是 $ (A^{-1})_{ij}$ 对 $A_{kl}$ 求导数 --- 逆矩阵的一个矩阵元对原矩阵元求偏导,看起来很不简单。但是稍微思考一下,发现有办法。 $A^{-1} A = 1$ 所以, $\partial A^{-1} \ A + A^{-1} \ \partial A =0$ 因此, $\partial A^{-1} = - A^{-1} \ \partial A \ A^{-1} $ 终于看到了曙光! $- \frac{\partial}{\partial A_{kl}} (A^{-1})_{ij} = (A^{-1})_{im} \ \frac{\partial A_{mn} }{\partial A_{kl} }\ (A^{-1} )_{nj} \\ = (A^{-1})_{im} \ \delta_{m,k} \delta_{n,l}\ (A^{-1} )_{nj} + (A^{-1})_{im} \ \delta_{m,l} \delta_{n,k}\ (A^{-1} )_{nj}\\ = (A^{-1})_{ik} \ (A^{-1} )_{lj} + (A^{-1})_{il} \ (A^{-1} )_{kj} $ 上面的负号正好跟前面的负号抵消。有了这个结果以及上面行列式求导的结果,相关的计算就是套下面这两个公式了: $\frac{\partial \det A }{\partial A_{ij}} = (A^{-1})_{ij} \det A$ $- \frac{\partial\ (A^{-1})_{ij}}{\partial A_{kl}} = (A^{-1})_{ik} \ (A^{-1} )_{lj} + (A^{-1})_{il} \ (A^{-1} )_{kj}$ 条条道路通罗马。至此,我们没有引入 bx 项作为生成器就得出了WICK 定理的结果。我个人认为,我上面的这个推导更为直接。
61 次阅读|0 个评论
分享 从行列式到量子物理:中国古代数学的新生命
热度 2 岳东晓 2017-1-18 09:49
爱斯基摩人的主要食物是海豹,但据说有一支爱斯基摩人的语言里没有海豹一词;每当爱斯基摩猎人相遇交换海豹信息的时候,需要说一大堆描述性的话才能说明是海豹这个动物。可以想象,这个表达的笨拙大大地约束了爱斯基摩的文学发展。假如他们要歌颂两个海豹猎人的传说,将会非常难以表达、记忆、传承。类似的,数学符号相当于数学语言的单词。古代中国人数学水平确实不低,但后来发展缓慢,为什么呢?我认为,原因之一是中国没有意识到符号表达数学关系的威力,没有发明各种的符号扩充数学语言的词汇量。如果把数学看成一个语言艺术领域,中国数学作品的相对贫乏 就不可避免了 。 行列式(Det: DETERMINANT)中学就学了。一个方形矩阵计算行列式,算大的矩阵要先算小的,符号还得交替变换。知其然不知其所以然,一头雾水。后来读西方写的数学史,才知道中国人早在《九章算术》(公元前2世纪)里就用到了行列式,而西方数学在两千年之后才用到。我们得承认,我们的古代祖先确实聪明,但后人就不行了,到唐朝时《九章算术》很多题目已经没人会做了。而西方则把这些数学步步推进,行列式在数学、物理经常出现。这个 det 符号往往令人望而生畏。它到底是什么? 先看一个简单的二元一次方程组: a x + b y = e c x + d y = f 简单的消元 ca x + cb y = c e ac x + ad y = a f 得出 (ad - cb ) y = af - ce 类似的: (da - bc ) x = de - bf 对比原来的方程组,这样对角相乘、然后相减的数 ad -bc 出现了,同一行、同一列的数字不会出现。这个数字我们称之为行列式。如果扩展到三元一次方程组,也有类似的数出现,但公式就复杂多了,因为有9个数字进行三个相乘(但同一行、同一列数字不会出现在一个乘积里)。九章算术里面解多元方程就是这么列阵进行。继续计算,可以发现中学数学里学到的递归计算规律。这是一个(线性)代数的理解。 数学与物理中,我们往往发现仅仅是符号的变化就能大大的简化各种推导与表达,甚至使很多看似复杂的结果变得非常显然。麦克斯韦尔的电磁学 论文发表时还没有矢量符号与偏微分符号,他的论文列出20多个分量方程,看得人眼花缭乱,根本不像现在这么优美简洁。狄拉克发明了 bra - ket 之后,量子力学的各种计算都大大简化,几乎成了机械套用。杨振宁当年为了找到规范场的数学 表达费了很长时间摸索,用 differential form 看,就是几分钟的事情 。在 这篇博文里 ,我做了一个CHERN- SIMONS场方程的推导,仅仅是几行;其中一行 $(\partial^a\epsilon_{ba\nu} - \sigma g_{b\nu})(\partial_{\mu} \epsilon^{f\mu\nu}-\sigma g^{\nu f}) B_f=0$,如果不是使用所谓 Einstein 的重复求和等规则,而是把各个分量摊开,这么多上标、下标,估计需要写满好几张纸,会看得我头昏眼花。 从几何上理解行列式,将每一列视为一个向量,行列式就是这些向量构成的多维多面体的体积(有朝向的)。在数学上,与向量结合而产生数字并且满足线性关系的量叫着张量。这个把多个向量吃入,输出体积的张量就是我们经常遇到的Levi-Civita 符号 $\epsilon{...}$ 。这是一个全反对称张量,当下标是 12... 时,数值为1,交换两个下标数字符号相反。 二维的情况: $\epsilon_{12}=-\epsilon_{21} =1, \epsilon_{11}=\epsilon_{22}=0$ 。 三维情况: $\epsilon_{123} = - \epsilon_{132} = - \epsilon_{213}= \epsilon_{231}= - \epsilon_{321} = \epsilon_{312}=1$,其余均为零。 用这个 $\epsilon_{...}$ 符号, n x n的矩阵 $A$ 的行列式 为 $\det A = \epsilon_{i_1 .. i_n} A_{1\ i_1}\cdots A_{n\ i_n}$ 这个表达适用于任意的 n 。以 2 x 2 矩阵为例, $\det A = \epsilon_{12} A_{11} A_{22} + \epsilon_{21} A_{12} A_{21} = A_{11} A_{22} - A_{12} A_{21}$ ,对角相乘减去斜对角相乘,这是中学数学都学了的。有兴趣地读者可以试试 3 x 3 矩阵的情况。 中国人发明的行列式在物理中运用极为广泛,但也是一种大大扩展了的运用。古希腊人把行星的运动用圆来解释,不行的话大圆加上小圆。类似的,物理中最基本构件是简谐振子。最基本的物理规律可以说就是没有规律。据 DYSON 回忆,费曼曾对他说一个粒子的运动其实是任意的,选择任何路径,可以跑到月球然后再回来,你把这些路径的几率(复数振幅)加起来就得到了实际结果。DYSON当时回答说:你疯了!当然我们知道,费曼没有疯。恰恰相反,他发现了最深刻的自然规律:路径积分。用费曼的路径积分来分析弹簧振子,把弹簧振子所有可能的路径加起来,包括飞出银河系再回来,有的读者可能会说,弹簧超光速拉伸出银河系早崩 断 -- 疯了。不是这样,费曼怎么是天才呢?总之,我们用费曼方法应该能得到这个经典力学弹簧的结果:在时间为弹簧振子的经典周期的时候,它应该有很大的几率复位。弹簧振子的拉格朗日为$L = \frac{1}{2} m v^2 - \frac{1}{2} k \ x^2$。剩下的这个路径积分是一个数学问题,简言之就是将所有可能的路径的作用量作为几率相角,然后将几率相加。这听起来可能令数学家们头皮发麻,但是理论物理却是家常便饭。下面我略加演示,请大家注意这个 Det 的出现。 首先, $L = \frac{1}{2} m (\frac{dx}{dt})^2 - \frac{1}{2} k \ x ^2 = \frac{1}{2} m - \frac{1}{2} k \ x^2\\= - \frac{1}{2}\ x\ \ x +\frac{1}{2} m \frac{d}{dt} (x \frac{dx}{dt})$, 令 $A = - m\ \frac{d^2}{dt^2} - k $, 则从 t=0, x=0, 到 t=T, x=0 的路径积分为 (自然单位 $\hbar=1$), $G = \int \mathcal {D}x \ e^ {i \int L(x)\ dt} = \int \mathcal{ D} x\ \exp(\frac{i}{2}\int dt \ x \ A\ x ) = \frac{C}{\sqrt{\det A}}$ 因为边界条件,L中的全微分项没有了。上面的积分是个高斯积分,C是一个常数。如果A 是一个矩阵,应该不难理解。把 A 对角化,对角相乘就是该矩阵的行列式。但我们的 A 不是矩阵,而是 一个微分算符。也就是说,上面的结果是 $G \propto \left ^{-\frac{1}{2}}$ 最早出现于九章算术的行列式 det 在现代物理学中大放光彩。
个人分类: 科普|6171 次阅读|0 个评论

Archiver|手机版|珍珠湾全球网

GMT+8, 2024-6-1 07:20 , Processed in 0.042401 second(s), 10 queries , Apc On.

Powered by Discuz! X2.5

回顶部