经典力学-太阳对光的吸引.docVIP

太阳对光的吸引 只考虑太阳的质量，和光子的能量，看看太阳能够让光子偏转多少？ 先看经典力学下的太阳对光的偏转。根据万有引力，可以得到运动方程 其中，u=1/r, l=mvr=mωr2。所以，上述方程还可以简化为： d2udθ2+u=GML2 其中，L是vr或ωr2，是常数。于是，在经典力学中，受重力影响而发生偏转，好像和质量无关了。即使光子在经典力学中没有质量，也是可能发生偏转的。 太阳 太阳半径：6.960*108m，太阳质量：1.989*1030kg。 所以：GM/L2=GM/(CR)2=3.048e-15 参考前面（另外一篇文章）对万有引力公式的推导，可以得到这个方程的解是： 1r=GML2+? cos?(θ+θ0) 根据初始条件：光线从远处过来，r-∞，θ-0。所以，?= - GML2。但是这种初始条件是无法考察光线偏转的。（实际看看就知道了） 要考察光线偏转，需要的目标是：光线正好从太阳表面掠过。也就是θ=90度的时候，r正好是太阳半径R。要满足这个条件，可以简化θ0=-π/2，这样简化后，边界条件就成为： 1R=GML2+?= ?= 1R-GML2 根据这个值，再计算入射光线和出射光线。可以计算θ=0和π时候的r值，得到： |r0|=|rπ|=L2/GM θ L2/GM L2/GM R 根据上图，可以得到： θ=2*arctan(R/(L2/GM))=2*arctan(GM/(RC2))=2*2.1e-6rad=2*00.875。 （1度=3600) 如果要再精确一些，可以算出在0度和180度的速度方向，再看偏转角度。绘制上面方程的曲线如下： 1r=GML2+? sin(θ) // scilab stacksize(max); G=6.672e-11; //万有引力 C=28; // 光速 SR=6.960e8; // radiu of sun SM=1.989e30; // mass of sun K=G*SM/(C*SR*C*SR); PS=1/SR-K; X=1; Y=K+PS*sin(sita); r=X./Y; polarplot(sita,r); 在太阳附件（π/2）局部放大图： 耐住性子，分析一下0度和180度附件的曲线的斜率，也就是速度方向。其实就是计算切线。 y=ρsin?(θ)x=ρcos?(θ)= dy=sinθdρ+ρcosθdθdx=cosθdρ-ρsinθdθ = dydx=sinθdρ+ρcosθdθcosθdρ-ρsinθdθ= sinθdρ/dθ+ρcosθcosθdρ/dθ-ρsinθ 令GML2=K，那么曲线1r=K+? sin(θ)的dr/dθ为： drdθ=-?cos?(θ)(K+?sinθ)2 所以在θ=0处的斜率为： dydx|θ=0=ρ|θ=0dρ/dθ|θ=0=-1/K?*K2=-K?=-GML21R-GML2=-3.048*10-151.437*10-9=-2.1*10-6 斜率角度为：0.438。 所以在θ=π处的斜率为： dydx|θ=π=ρ|θ=πdρ/dθ|θ=π=K? 所以不管怎么算，偏转角度还是2*0.438=0.876。 这是经典力学的分析结果。