非谐势阱中玻色-爱因斯坦凝聚的数值计算

从G-P平均势场理论出发,探讨了玻色-爱因斯坦凝聚(BEC)的G-P方程的一维形式,用数值计算方法研究了非谐势阱中非理想玻色凝聚气体的基态和第一激发态解.给出了能量随非线性系数的变化规律.     

三维非谐势阱中玻色－爱因斯坦凝聚的数值计算

本文从G－P平均势场理论出发,探讨了三维球对称非谐势阱中玻色－爱因斯坦凝聚(BEC)的G－P方程;用数值计算方法研究了三维球对称非谐势阱中原子间有相互作用的玻色－爱因斯坦凝聚气体的基态解;分析了非谐振势能项对玻色－爱因斯坦凝聚体的分布、能量和化学势的影响。     

非谐势阱中二维G-P方程的数值计算

通过能量泛函的方法得到描述囚禁在非谐势阱中玻色-爱因斯坦凝聚体的二维G-P方程数值解,讨论原子间相互作用和非谐振势能项对玻色-爱因斯坦凝聚体的分布、能量和化学势的影响.     

三维简谐势阱中玻色-爱因斯坦凝聚的边界效应

在定义特征长度的基础上,应用Euler–MacLaurin公式,研究了理想玻色气体在三维简谐势阱中玻色-爱因斯坦凝聚的边界效应.结果表明:粒子的凝聚分数由于有限尺度和有限粒子数效应而减小,修正的凝聚分数和凝聚温度由于边界效应存在一个极大值,选择优化的最佳势阱参数,可以有效提高凝聚分数和凝聚温度;热容量的跃变存在边界效应和粒子数效应,选择合理的势阱参数时,热容量的跃变存在一个极小值.导出了简谐势阱中有限理想玻色气体的状态方程,揭示了压强的各向异性(或各向同性)取决于简谐势频率的各向异性(或各向同性).     

谐振子势阱囚禁玻色气体的玻色-爱因斯坦凝聚

基于Thomas-Fermi半经典近似方法研究了谐振子势阱约束下任意维理想玻色气体的玻色-爱因斯坦凝聚(BEC)．导出了玻色气体的BEC转变温度、基态粒子占据比例、内能和热容量等物理量的解析表达式,讨论了空间维度和谐振子势阱的影响．以二维和三维玻色系统为例,数值计算了上述热力学量,并与解析结果进行了对比,二者获得了较好的吻合．     

谐振子势阱囚禁玻色气体的玻色-爱因斯坦凝聚

基于Thomas-Fermi半经典近似研究了谐振子势阱约束下任意维理想玻色气体的玻色-爱因斯坦凝聚(BEC).导出了玻色气体的BEC转变温度、基态粒子占据比例、内能和热容量等物理量的解析表达式,讨论了空间维度和谐振子势阱的影响.以二维和三维玻色系统为例,数值计算了上述热力学量,并与解析结果进行了对比,二者获得了较好的吻合.     

球对称非谐势阱中玻色-爱因斯坦凝聚Gross-Pitaevskii方程的精确解

考虑了描述玻色 爱因斯坦凝聚的Gross-Pitaevskii(GP)方程, 得到了在球对称非谐势阱中玻色-爱因斯坦凝聚GP方程的精确亮孤子解。In this paper, we analyze Gross Pitaevskii equation which describes the dynamics of a bright soliton in trapped atomic Bose Einstein condensates, and obtain the exact bright soliton solution of Gross Pitaevskii equation in spherically symmetric non harmonic trap.     

势场中玻色-爱因斯坦凝聚的临界温度

给出了不同于文献的势场中玻色-爱因斯坦凝聚临界温度表达式.结果揭示了势场中理想玻色子气体凝聚的临界温度与势场之间的关系，表明势场中临界温度正比于无势场情况下临界温度T０ｃ，还给出了势场的有效性判据.势场的有效性是势场与玻尔兹曼常数k和无势场情况下临界温度Ｔ０ｃ乘积ｋＴ０ｃ的比较.当势场接近或大于ｋＴ０c时，临界温度会有效增加;当势场远小于ｋＴ０c时, 势场是无效的. 关键词： 玻色-爱因斯坦凝聚 临界温度 势阱     

谐振势阱中有弱相互作用的玻色爱因斯坦凝聚

应用数值计算的方法计算了谐振势阱中有弱相互作用的玻色气体凝聚的临界温度和基态占据库，计算结果表明，二者都随着散射长度的增大而减小；但与理想玻色气体相比仅差约０．４％。     

三势阱中玻色-爱因斯坦凝聚的开关特性

应用平均场近似的方法,研究了弱耦合的三势阱中玻色-爱因斯坦凝聚的开关效应.当粒子置于左阱时,可以通过在中间势阱中加入少量粒子控制左阱粒子向右阱的隧穿,从而呈现出明显的导通与截止行为.对中间势阱的深度和相对相位的影响也进行了讨论,并指出了该理论模型的一些潜在应用前景. 关键词： 玻色-爱因斯坦凝聚 开关效应 三势阱 平均场近似     

玻色-爱因斯坦凝聚在中国科学院上海光机所实现

文章报道了在稀薄87Rb原子气体中观测到的玻色-爱因斯坦凝聚(BEC)现象.在四极矩和Ioffe组合磁阱(QUIC)中装载了1×10\\+8个原子,经过19s蒸发冷却达到了相变条件.在高原子密度的情况下,文章作者观察到了BEC对探测光的衍射光环.这时降低磁阱势垒和绝热的放开冷原子样品,我们拍摄到冷原子和BEC的吸收像.根据数据拟合满足双高斯分布,表明发生了BEC相变.相变温度约215nK,凝聚的原子数约为5×10\\+4.     

Numerical solution of three-dimensional heterogeneous solid propellants

We have described, for the first time, a fully coupled low Mach number numerical algorithm which can be used to investigate the combustion of composite propellants. The code uses a body fitted grid along the moving interface by means of a mapping technique. Specifically, the surface is assumed to be single valued, thus allowing the use of a mapping function that effectively maps the propagating corrugated surface into a stationary flat surface. The transformed connection conditions are then applied along the stationary flat surface, allowing second-order one-sided derivatives to be used. This has proven to be a very robust and efficient way to treat the surface and the connection conditions. Numerical tests are performed and the scheme is shown to be second-order accurate in the spatial directions as well as in time. Selected results are presented for heterogeneous propellants.     

Oscillations of a bose-einstein condensate rotating in a harmonic plus quartic trap

We study the normal modes of a two-dimensional rotating Bose-Einstein condensate confined in a quadratic plus quartic trap. Hydrodynamic theory and sum rules are used to derive analytical predictions for the collective frequencies in the limit of high angular velocities Omega where the vortex lattice produced by the rotation exhibits an annular structure. We predict a class of excitations with frequency sqrt[6]Omega in the rotating frame, irrespective of the mode multipolarity m, as well as a class of low energy modes with frequency proportional to |m|/Omega. The predictions are in good agreement with results of numerical simulations based on the 2D Gross-Pitaevskii equation. The same analysis is also carried out at even higher angular velocities, where the system enters the giant vortex regime.     

简谐势阱中非理想气体玻色-爱因斯坦凝聚转变温度的数值研究

用数值计算的方法研究了简谐势阱中非理想气体的玻色 爱因斯坦凝聚情况.通过计算能级的分布，用图形的方式给出了考虑相互作用时，简谐势阱中体系转变温度相对无相互作用体系转变温度的变化趋势 对于存在排斥作用的体系，转变温度升高；对于存在吸引作用的体系，转变温度降低. 关键词： 玻色－爱因斯坦凝聚 相互作用 能级 转变温度     

Effect of quantum fluctuations on bose-einstein condensation of bilayer atomic gases in rapid rotation

We investigate vortex states and quantum fluctuations of bilayer atomic gases in rapid rotation. Among mean-field solutions of the Gross-Pitaevskii equation, we consider two types of vortex configurations: the vortex core positions of the layers are coincident or staggered. It is found that the coincident type is energetically preferred in the practical parameter regime. We also calculate the dispersion relations of collective modes and the filling factor for bosons outside the condensates. In the double layer system, quantum depletion is found to be suppressed due to interlayer tunneling. This means that Bose-Einstein condensation is stabilized compared to the single layer case.     

数值求解三维含时Schrodinger方程及其应用

发展了一套高精度、高效率的伪谱方法,以非微扰的方式求解真实原子三维含时Schrodinger方程.该方法选用二阶劈裂算符作为时间演化算子,分别选择能谱表象和坐标表象作为含时波函数演化的两个表象.在坐标表象下波函数的径向部分使用库仑波函数离散变量表象来离散;角向波函数展开在两维的Gauss-Legendre-Fourier格点上.以H原子的光激发和光电离过程为例,进行了数值计算并和解析解进行了比对.结果表明二者符合很好.该方法很好地处理了库仑奇点问题.还计算了强激光辐照H原子的多光子电离过程,并和其他的数值方案进行了比较.结果表明,在计算收敛的前提下本方法计算效率更高.

Abstract：

We present an accurate and effective pseudospectral method for solving the three-dimensional time-dependent Schrodinger equation involving the Coulomb potential. In this method, the Hamiltonian is evaluated by exploiting the two representations of the wave function. One is a grid representation, in which the angular dependence of the wave function is expanded in a two-dimensional Gauss-Legendre-Fourier grid in the coordinate space of polar and azimuthal angles. The radial coordinate is discretized using a discrete variable representation constructed from the Coulomb wave function (CWDVR) . The other is a spectral representation, in which the wave function is expanded in a set of square integrable functions chosen as the eigenfunctions of a zero-order Hamiltonian. The time of propagation of the wave function is calculated using the well-known second-order split-operator method implemented through the transform between the grid and spectral representations. Calculations on the photo-absorption strength of hydrogen atom are presented to demonstrate the accuracy of present method in low energy limit by the time-dependent wave-packet propagation method. As another example, the present method is applied to multiphoton ionization of H atom. For a wide range of field parameters, ionization rates calculated using the present method are in excellent agreement with those from other accurate numerical calculations. The new algorithm will be found more efficient than the close coupled wave packet method using CWDVR and/or methods based on evenly spaced grids.     

Thermodynamic evidence for nanoscale bose-einstein condensation in (4)He confined in nanoporous media

We report the measurements of the heat capacity of (4)He confined in nanoporous Gelsil glass that has nanopores of 2.5-nm diameter at pressures up to 5.3 MPa. The heat capacity has a broad peak at a temperature much higher than the superfluid transition temperature obtained using the torsional oscillator technique. The peak provides definite thermodynamic evidence for the formation of localized Bose-Einstein condensates on nanometer length scales. The temperature dependence of the heat capacity is described well by the excitations of phonons and rotons, supporting the existence of localized Bose-Einstein condensates.     

Bose-Einstein condensation of an ideal gas in a parabolic trap

The temperature-dependent velocity distribution function is found for the case of Bose-Einstein condensation of a finite number of noninteracting atoms trapped in a three-dimensional anisotropic parabolic trap. It is shown that at a temperature T of the order of the condensation temperature T ₀ the velocity distribution consists of an anisotropic part, reflecting the population of the ground state, and an isotropic part above the condensate. Pis’ma Zh. éksp. Teor. Fiz. 66, No. 8, 559–563 (25 October 1997)