Quark gluon plasma

Recent trends in the research of quark gluon plasma (QGP) are surveyed and the current experimental and theoretical status regarding the properties and signals of QGP is reported. We hope that the experiments commencing at relativistic heavy-ion collider (RHIC) in 2000 will provide a glimpse of the QGP formation.     

Status of the quark gluon plasma search

A selection of results are discussed that support the conclusion that strongly interacting quark gluon plasma is produced in heavy-ion collisions at the Relativistic Heavy Ion Collider at BNL.     

Hard loop effective theory of the (anisotropic) quark gluon plasma

The generalization of the hard thermal loop effective theory to anisotropic plasmas is described with a detailed discussion of anisotropic dispersion laws and plasma instabilities. The numerical results obtained in real-time lattice simulations of the hard loop effective theory are reviewed, both for the stationary anisotropic case and for a quark-gluon plasma undergoing boost-invariant expansion.     

Instability of the plasma oscillations in finite temperature perturbative QCD

We use linear-response dielectric theory to show that the baryon-poor QCD plasma based on the perturbative vacuum is unstable, even at a high temperature. If deconfinement occurs in nuclear collisions or the early universe, it is not accompanied by the restoration of the perturbative vacuum.     

Viscous quark‐gluon plasma in the early universe

In the present work a study is given for the evolution of a flat, isotropic and homogeneous Universe, which is filled with a causal bulk viscous cosmological fluid. We describe the viscous properties by an ultra‐relativistic equation of state, and bulk viscosity coefficient obtained from recent lattice QCD calculations. The basic equation for the Hubble parameter is derived by using the energy equation obtained from the assumption of the covariant conservation of the energy‐momentum tensor of the matter in the Universe. By assuming a power law dependence of the bulk viscosity coefficient, temperature and relaxation time on the energy density, we derive the evolution equation for the Hubble function. By using the equations of state from recent lattice QCD simulations and heavy‐ion collisions we obtain an approximate solution of the field equations. In this treatment for the viscous cosmology, no evidence for singularity is observed. For example, both the Hubble parameter and the scale factor are finite at t = 0, where t is the comoving time. Furthermore, their time evolution essentially differs from the one associated with non‐viscous and ideal gas. Also it is noticed that the thermodynamic quantities, like temperature, energy density and bulk pressure remain finite. Particular solutions are also considered in order to prove that the free parameter in this model does qualitatively influence the final results.     

The extent of strangeness equilibration in quark gluon plasma

The evolution and production of strangeness from chemically equilibrating and transversely expanding quark gluon plasma which may be formed in the wake of relativistic heavy-ion collisions is studied with initial conditions obtained from the self screened parton cascade (SSPC) model. The extent of partonic equilibration increases almost linearly with the square of the initial energy density, which can then be scaled with the number of participants.     

瞬发中子密度衰减法计算中子代时间

采用蒙特卡罗程序MCNP计算了西安脉冲堆中子代时间。使用MCNP程序模拟了反应堆瞬发中子通量密度衰减,基于忽略缓发中子项的点堆动力学方程计算出中子代时间。在微次临界下,研究了次临界度、源的分布、计数区域等对西安脉冲堆中子代时间计算结果的影响。计算分析表明:采用瞬发中子密度衰减法计算中子代时间时,微次临界度、源分布、计数区域等对计算结果影响都很小;误差产生的主要原因是忽略缓发中子项的点堆动力学方程并不能较好地反应瞬发中子通量密度的衰减规律。     

瞬发中子密度衰减法计算中子代时间

采用蒙特卡罗程序MCNP计算了西安脉冲堆中子代时间。使用MCNP程序模拟了反应堆瞬发中子通量密度衰减,基于忽略缓发中子项的点堆动力学方程计算出中子代时间。在微次临界下,研究了次临界度、源的分布、计数区域等对西安脉冲堆中子代时间计算结果的影响。计算分析表明：采用瞬发中子密度衰减法计算中子代时间时,微次临界度、源分布、计数区域等对计算结果影响都很小;误差产生的主要原因是忽略缓发中子项的点堆动力学方程并不能较好地反应瞬发中子通量密度的衰减规律。     

氮等离子体枪中电子输运行为的蒙特卡罗模拟

采用蒙特卡罗模拟,对氮等离子体枪的两个环形电极间的氮气辉光放电过程中电子的输运过程进行研究.计算在不同平均电场与粒子数密度的比值(E/N)下,电子与氮分子发生不同碰撞的概率、出射电子的平均能量、方向角分布和电子的能量分布.结果表明,电子能量近似服从玻尔兹曼分布.随着E/N的升高,电子平均能量升高,发生激发、离化、电离和离化电离碰撞的概率增大;非均匀分布的电场使分子获得更高的离化率,同时显著增强出射电子的能量.模拟结果为等离子体应用设计提供了参考依据.     

多群输运-燃耗耦合蒙特卡罗计算

介绍了三维多群中子输运-燃耗耦合P3近似蒙特卡罗程序MCMG-BURN，该程序是在连续截面蒙特卡罗程序MCNP和反应堆栅元均匀化程序WIMS基础上发展而来的，使用多群截面为模拟临界实验堆和高通量工程实验堆，取得与MCNP和实验一致的结果，MCMG-BURN具有与MCNP相同的精度，但计算时间较MCNP要少得多。     

非均匀磁场约束条件下的电子束空气等离子体特性

为了揭示大气环境电子束等离子体的性质,基于蒙特卡罗程序包Geant4建立了一个包含电离、激发以及轫致辐射等物理过程的计算模型,用于模拟非均匀磁场约束条件下高能强流稳态电子束的输运特性、以及大气环境等离子体的性质。结果表明:非均匀磁场可以有效控制电子束在空气中的输运轨迹,显著降低电子束的发散;随着电子束在空气中行程的增加,电子束能谱开始展宽并向低能区移动;输运装置出口能量损失比电子束射程末端高2个量级,且随着电子束输运距离的增加,等离子体密度降低;等离子体密度的高低与电子束能量直接相关。     

蒙特卡罗模拟在辉光放电鞘层离子输运研究中的应用

利用蒙特卡罗方法模拟氩气直流辉光放电鞘层内离子的运动过程.模拟基于离子与中性原子的电荷转移和弹性散射两种主要的散射过程，考虑了碰撞截面与能量相关和不相关两种情况，在弹性散射中采用了势场相互作用模型和刚性球碰撞两种模型.通过模拟得到不同气压和不同放电电压下离子入射阴极的能量分布和角度分布，并对几种模型的模拟结果进行了比较和讨论. 关键词： 辉光放电 等离子体鞘层 蒙特卡罗模拟     

流注放电低温等离子体中电子输运系数计算的蒙特卡罗模型

流体或者粒子-流体混合数值仿真是研究流注放电基本物理机制的常用手段,而精确的电子输运系数是保证其仿真正确性的必要前提.鉴于现有电子输运系数求解工具存在一定缺陷,本文开发了采用蒙特卡罗方法求解低温等离子体中电子输运系数的仿真工具,测试表明其准确性和精确度均较高.研究了氮氧气体混合比及大气压下三体碰撞吸附对电子输运系数的影响.氮气中流注放电仿真表明,流体仿真中采用本模型改进后的电子输运系数可显著改善流注通道内部的等离子体参数分布.     

Measurement of temperature of a dusty plasma from its configuration

A new concept called “configurational temperature” is introduced in the context of dusty plasma, where the temperature of the dust particles submerged in the plasma can be measured directly from the positional information of the individual dust particles and the interaction potential between the dust grains. This method does not require the velocity information of individual particles, which is a key parameter to measure the dust temperature in the conventional method. The technique is initially tested using two-dimensional (2D) OpenMP parallel molecular dynamics and Monte Carlo simulation and then compared with the temperature evaluated from experimental data. The experiment have been carried out in the Dusty Plasma Experimental (DPEx) device, where a 2D stationary plasma crystal of melamine formaldehyde particles is formed in the cathode sheath of a DC glow discharge argon plasma. The kinetic temperature of the dust is calculated using the standard particle image velocimetry technique at different pressures. An extended simulation result for the three-dimensional case is also presented, which can be employed for the temperature measurement of a three-dimensional dust crystal in laboratory devices.     

铪离子等离子体源离子注入铜基体的数值模拟

通过SRIM软件对铪离子等离子体源离子注入铜进行了模拟。模拟了铪离子注入铜的核阻止本领、电子阻止本领、入射深度随能量的变化,以及在不同注入条件下铪离子的摩尔浓度分布,并对模拟结果进行了分析。结果显示:能量低于6 MeV时核阻止本领占主导地位,高于6 MeV时电子阻止本领成为主要的能量损失,并且离子注入过程中会出现能量沉积的Bragg峰和质量沉积区域较集中的现象,入射深度随能量的增加而增加。     

增强辉光放电等离子体离子注入的三维PIC／MC模拟

采用三维粒子模拟／蒙特卡洛模型自洽地模拟了增强辉光放电等离子体离子注入过程中离子产生和注入,获得了放电空间的离子总数、电势分布、等离子体密度分布和离子入射剂量等信息．模拟结果表明,5μs时鞘层达到稳定扩展,15μs时离子的产生与注入达到平衡,证实了增强辉光放电等离子体离子注入能在一定条件下实现白持的辉光放电．注入过程中,在点状阳极正下方存在一个高密度的等离子体区域,证实了电子聚焦效应．除靶台边缘外,离子的注入速率稳定且入射剂量均匀．脉冲负偏压提高时注入速率增加但入射剂量的均匀性变差．     

A Monte Carlo study of the classical two-dimensional one-component plasma

We present results from extensive Monte Carlo simulations of the fluid phase of the two-dimensional classical one-component plasma (OCP). The difficulties associated with the infinite range of the logarithmic Coulomb interaction are eliminated by confining the particles to the surface of a sphere. The results are compared to those obtained for a planar system with screened Coulomb interactions and periodic boundary conditions; in this case the infinite tail of the Coulomb interaction is treated as a perturbation. The exact simulation results are used to test various approximate theories, including a semiempirical modification of the hypernetted-chain (HNC) integral equation. The OCP freezing transition is located at a coupling= e²/k_BT–140.     

On surface properties of a weakly-coupled classical one-component plasma: The size effect in finite systems

The surface properties of a weakly coupled classical one-component plasma of finite size are calculated exactly within the Poisson-Boltzmann (PB) approximation scheme. It is found that the ion density profile and the surface energy for a spherical system show strong size dependence. The surface energy also strongly depends on the position of the hard wall introduced for achieving an appropriate equilibrium ion configuration. These results indicate that the recent Monte Carlo simulation data for a spherical system must be interpreted, at least in the weak-coupling regime, as including substantial size effects and cannot be directly compared with the theoretical calculations for the planar surface. For a slab, on the other hand, such size effects are found to be very small if the hard wall is placed at sufficiently distant position from the surface. The dominant contribution to the surface energy which is omitted in the PB approximation is also estimated by the perturbation calculations.     

Approximating the number of monomer-dimer coverings of a lattice

We study the problem of counting the number of coverings of ad-dimensional rectangular lattice by a specified number of monomers and dimers. This problem arises in several models in statistical physics, and has been widely studied. A classical technique due to Fisher, Kasteleyn, and Temperley solves the problem exactly in two dimensions when the number of monomers is zero (the dimer covering problem), but is not applicable in higher dimensions or in the presence of monomers. This paper presents the first provably polynomial-time approximation algorithms for computing the number of coverings with any specified number of monomers ind-dimensional rectangular lattices with periodic boundaries, for any fixed dimensiond, and in two-dimensional lattices with fixed boundaries. The algorithms are based on Monte Carlo simulation of a suitable Markov chain, and, in constrast to most Monte Carlo algorithms in statistical physics, have rigorously derived performance guarantees that do not rely on any assumptions. The method generalizes to counting coverings of any finite vertex-transitive graph, a class which includes most natural finite lattices with periodic boundary conditions.