粒子模拟程序的发展及其在激光等离子体相互作用研究中的-计算物理.PDFVIP

计 算 物 理 第25卷第l期 V01．25．No．I CHINESE OFCOMPUTA，rIONALPHYSICS 2008年1月 JOURNAL Jan．．2008 【文章编号】1001．246X(2008)Ol-0043．08 粒子模拟程序的发展及其在激光等离子体 相互作用研究中的应用 陈 民， 盛政明， 郑 君， 马燕云， 张 杰 (中国科学院物理研究所北京凝聚态物理国家实验室，北京 100080) [摘 要】 介绍针对超短超强激光脉冲与等离子体相互作用研究的多维粒子模拟程序KLAP．在其一维程序 KIAPID中。考虑场电离、碰撞离化及两体碰撞效应后，程序可以用于研究短脉冲激光与中性物质的相互作用． 在其三维程序KLAP3D中，为了研究加速能量达GeV的长距离激光尾波场加速问题，程序采用移动窗口技术。使 得模拟尺度可以达到厘米量级．同时介绍了利用KLAP程序得到的有关THz辐射、激光与中性气体相互作用中的 脉冲及离化波前演化、激光固体靶作用中表面电子加速及激光尾波场加速的研究实例． [关键词】 粒子模拟；等离子体；超短超强激光；并行计算；移动窗口 [中图分类号】 0411．3；0532+．13[文献标识码】 A O 引言 近年来，随着超短超强激光脉冲技术的发展，相对论强度的激光脉冲与等离子体的相互作用引起了人们 的极大兴趣．它具有广泛的应用前景，如快点火激光核聚变…、新型粒子加速器乜]、高能离子源及其医学诊断 和治疗¨。51、新型辐射源№1等．与以往的长脉冲非相对论强度激光脉冲与等离子体作用不同，超短超强激光 脉冲与等离子体的相互作用具有作用时间短、强动力学效应、高度非线性、多过程耦合等特征．粒子之间在有 限短的时间内无法通过碰撞达到热力学平衡，不能用一个简单的流体形式描述该类粒子．此外在有些领域， 少数粒子的动力学特征更是人们研究的对象(如激光尾波场加速中的电子注入和捕获)．在这种情形下，把等 离子体视为单流或多流的流体力学行为已不再合适，通常的流体力学模拟也就难以胜任，人们便将模拟方法 转向描述更为细致的动力学模拟． 在目前应用最为广泛的激光等离子体动力学模拟中，一种是基于求解各种粒子在相空间的分布函数的 化为用其低阶矩方程所组成的流体力学模拟’9】．如果研究体系中粒子在相空间的分布较为广泛，且人们对 激光等离子体相互作用领域中另一动力学模拟方法HJ训． 由于采用有限大小的粒子云模型，常规的粒子模拟方法也忽略粒子之间的近程碰撞作用，这使得它有点 解法的复杂性，使得处理二维或三维空间问题时，计算时间远远超过PIC模拟；其次，当粒子在相空间的分布 具有锐的边界时(如激光与固体靶相互作用，存在锐的等离子体真空边界)，会造成分布函数的这种锐特征的 扩散效应，其频谱中的短波长成分会以错误的速度扩散．为了避免这一效应，通常是增大相空间中的分辨率 (Ax，△t，)，但会使模拟时间进一步增大】．由于粒子模拟方法采用的是求解较为简单的粒子运动方程，且数 值方法非常成熟，并行效率高，因而在目前的激光等离子体相互作用领域大有作为． ’ 目前国际上发展了很多粒子模拟程序．如美国加州大学洛杉矶分校等单位联合开发的OSIRIS程序¨¨， 【收稿日期】2006—10—19；[修回日期】2007—03—27 和科学院知识创新工程资助项目 【作者简介】陈民(1979一)，安徽马鞍山，博士生，主要从事相对论激光等离子体相互作用中的粒子加速及高能粒子输运等方面的研究． 计 算 物 理 第25卷 它是一个全相对论，全并行，全隐格式，且包含场及碰撞离化的大型程序．应用该程序在尾波场加速方面取得 LaserPlasma 了丰富的结果；Pukh