力梯度辛方法在圆型限制性三体问题中的应用

旋转坐标系下的圆型限制性三体问 题因含非惯性系所附加的影响部分使得动能不是动量的严格二次型, 可能导致力梯度辛积分算法的应用遇到困难. 从Lie算子运算出发, 严格论证了力梯度算子在这种情形下的物理意义 仍然像质心惯性坐标系下的圆型限制性三体问题那样是引力的梯度, 而不是引力与非惯性力所得合力的梯度, 表明了力梯度辛方法适合求解旋转坐标系下的圆型限制性三体问题. 通过应用四阶力梯度辛方法、最优化四阶力梯度辛方法和Forest-Ruth 辛方法分别求解该问题, 进行了数值对比研究, 结果显示最优化型力梯度算法能够取得最好精度. 还应用最优化型算法计算两邻近轨道的Lyapunov指数和快速Lyapunov指标, 确保高精度辛方法能够贯穿于这些混沌指标计算的全过程, 以便准确刻画此系统的动力学定性性质. 关键词： 辛积分器 圆型限制性三体问题 混沌 Lyapunov指数     

Dynamics of particles around a pseudo-Newtonian Kerr black hole with halos

The regular and chaotic dynamics of test particles in a superposed field between a pseudo-Newtonian Kerr black hole and quadrupolar halos is detailed.In particular,the dependence of dynamics on the quadrupolar parameter of the halos and the spin angular momentum of the rotating black hole is studied.It is found that the small quadrupolar moment,in contrast with the spin angular momentum,does not have a great effect on the stability and radii of the innermost stable circular orbits of these test particles.In addition,chaos mainly occurs for small absolute values of the rotating parameters,and does not exist for the maximum counter-rotating case under some certain initial conditions and parameters.This means that the rotating parameters of the black hole weaken the chaotic properties.It is also found that the counter-rotating system is more unstable than the co-rotating one.Furthermore,chaos is absent for small absolute values of the quadrupoles,and the onset of chaos is easier for the prolate halos than for the oblate ones.     

旋转致密双星的引力波特征

研究了轨道和旋转效果到2.5阶后牛顿旋转致密双星拉格朗日动力学与引力波的关系, 分析了有序和混沌轨道的引力波特征.发现当加速度不考虑辐射项时, 有序双星系统辐射的引力波具有周期或拟周期的特征, 而混沌双星系统辐射的引力波却具有明显的混沌特征.当加速度含有辐射项贡献时, 双星必会出现并合现象.此时, 原保守有序双星系统需较长时间才能完成并合过程, 引力波形在双星并合前仍保留拟周期的基本特点;然而, 原保守混沌双星系统仅在较短时间内就会并合, 但因并合时间太短, 无法获取足够的动力学信息导致引力波形的特征不易分辨.     

两个三阶最优化力梯度辛积分器的对称组合

利用已存在的三阶最优化力梯度辛格式以对称组合方法获得两个新的四阶力梯度辛积分器.它们在求解摄动Kepler混沌问题的能量精度和一维定态Schrö,dinger方程的能量本征值精度方面比Forest-Ruth四阶非力梯度辛积分器要好得多,甚至还要明显优越于已有的四阶最优化力梯度辛积分器.     

Performance of an irreversible quantum refrigeration cycle

A new model of a quantum refrigeration cycle composed of two adiabatic and two isomagnetic field processes is established. The working substance in the cycle consists of many non-interacting spin-1/2 systems. The performance of the cycle is investigated, based on the quantum master equation and semi-group approach. The general expressions of several important performance parameters, such as the coefficient of performance, cooling rate, and power input, are given. It is found that the coefficient of performance of this cycle is in the closest analogy to that of the classical Carnot cycle. Furthermore, at high temperatures the optimal relations of the cooling rate and the maximum cooling rate are analysed in detail. Some performance characteristic curves of the cycle are plotted, such as the cooling rate versus the maximum ratio between high and low ``temperatures' of the working substances, the maximum cooling rate versus the ratio between high and low ``magnetic fields' and the ``temperature' ratio between high and low reservoirs. The obtained results are further generalized and discussed, so that they may be directly applied to describing the performance of the quantum refrigerator using spin-$J$ systems as the working substance. Finally, the optimum characteristics of the quantum Carnot and Ericsson refrigeration cycles are derived by analogy.     

Dynamics of Charged Particles Around a Magnetized Black Hole with Pseudo–Newtonian Potential

The linear stability of equilibria of charged particles moving near a compact object with a dipole magnetic field and a pseudo-Newtonian potential is analyzed detailedly.An optimal fourth-order force gradient symplectic method,as a global symplectic integrator that can simultaneously solve both the equations of motion and the variational equations,is used to calculate fast Lyapunov indicators.In this way,dynamical structures are described,and parameter domains for causing chaos are found.     

半隐Euler法和隐中点法嵌入混合辛积分器的比较

当Hamilton函数分解为可积和不可积两部分时,前者能用分析方法给出解析解,而后者可借助一阶半隐Euler法或二阶隐中点法等数值求解,将这种解析和数值解法组合能构造二阶混合辛积分器.理论分析表明Euler嵌入法的稳定区要小于中点嵌入法的.再分别以圆形限制性三体问题和相对论自旋致密双星后牛顿Hamilton构型为例,详细比较了两嵌入法的性能特点.二者的数值精度、稳定性及计算效率与Hamilton的分解方式和轨道类型有关.就圆形限制性三体问题而言,当Hamilton采用势能和含坐标与动量混合项在内的动能分解 关键词： 辛积分器 后牛顿近似 自旋致密双星 混沌