matlab二维势阱简谐振动程序,二维简谐势阱中的理想气体玻色凝聚的数值分析及模拟...

1

2

[/(2)]c T N ζθ= (10) ()x ζ为黎曼函数，当1000N =时，24.66c T θ=，通过对若干点的计算，确定处于基态的粒子的个数，如表1所示，(0N 去除小数部分)

/c T T

0.6 0.83 0.86

0.890.920.950.98 1.0 1.01 1.04 1.07 1.1 1.20N 595 258 206 155 107 66 37 26 22 14 10 8 4

表1.各温度下处于基态粒子数(1000N =)

做出相关图像，如图5

图5.基态占有率随温度变化的关系

当100N =时7.8c T θ=，同样的做出表2。 表2.各温度下处于基态粒子数(100N =)

为便于后面的讨论，将1000N =和100N =时0//N N T θ?的图像做于同一坐标中，如图6所示

/T θ 4 5 6

7.8 9 10 15 20 0N 68 52 33 8 3 2 0.58940.2873