一維的電漿模型常被用作研究真實電漿的熱力學性質或遍歷行為;一維的電漿薄板模型 (plasma sheet model, PSM) 即為一個多粒子系統的典型例子。此篇論文所論及的電漿薄板模型,將薄板視作電子的類比,而離子則被簡化成固定不動的均勻背景,所有的薄板只能在該背景之中活動。熱鬆弛是熱力學系統由非熱平衡狀態趨向熱平衡的過程。倘若僅考慮雙粒子之間的相關性,由若干相同的一維粒子所構成的系統,並不具有熱鬆弛的現象。然如先前所述的電漿薄板模型,由於薄板與薄板之間或薄板與背景之間的庫倫作用力的緣故,當薄板的初始速度的機率密度分布,呈現均勻狀或拋物線形,其速度分布均會演變成馬士威爾分布,即該系統會趨向熱平衡。此庫倫作用力可當作三體粒子之間的相關性。此篇論文涵蓋的議題如下:一、 演算法對系統總能量的影響。 程式採取的演算法是將時間先行離散化,再計算各個時間區段的薄板軌跡。 因為解出的軌跡是近似解,因此有必要分析區段的長度與總能量變化的關 係,結果是總能量變化正比於區段長度的三次方。二、 達成熱平衡時,電位能與動能的比值。三、 由薄板的速度計算而得的卡方統計量,隨時間的變化。 卡方統計量除了針對薄板的速度分布會演變成馬士威爾分布的論點給出 統計學的依據之外,其隨時間的變化尚有助於測量熱鬆弛的特徵時間。卡方 統計量的漸進行為是指數衰減。四、 熱鬆弛的特徵時間與nD的相依性。 利用卡方統計量所測得的熱鬆弛的特徵時間與nD 2 大致成正比,同Dawson 的結論相吻合。鑒於兩種初始速度分布,不同於Dawson僅使用均勻分布, 遂採用次方擬合以解析特徵時間與nD 的某次方的正比關係,與初始速度分 布之間的關聯。拋物線形對應的次方比均勻狀的更加接近2,原因或許為前 者較後者更像馬士威爾分布。One-dimensional plasma models exhibit some thermodynamic properties or ergodic behavior of real plasmas; of many-body systems, a one-dimensional plasma model, consisting of a number of identical charged sheets embedded in a uniform fixed neutralizing background, is such a representative one. Thermal relaxation is an evolution during which a thermodynamic system approaches to thermal equilibrium. If a system of identical one-dimensional particles takes into account only two-particle correlations, no thermal relaxation takes place. Rather, due to the long-range forces among the background and the sheets themselves, as the initial velocity of sheets is either square or quadratic, the velocity of sheets will relax to Maxwellian; the system will attain to thermal equilibrium.The topics in this thesis:1. Energy conservation in the numerical schemeThe program deals with the discretized time to solve the approximated traces of sheets; therefore, it is necessary to examine the influence of the period ?t of a time interval on the total energy per sheet. It turns out that the total energy per sheet is dependent on the cube of the time increment.2. The ratio of the electric energy per sheet to the kinetic energy per sheet in thermal equilibrium3. Time development of the ??2 statistic constructed by the velocities of sheetsIn addition to provide statistical evidence that the velocity distribution of sheets relaxes to Maxwellian, the chi-squared statistic measures the thermal relaxation time, decaying exponentially asymptotically.4. The dependence of the thermal relaxation time ?R on nDBy use of the quadratic fitting, ?R seems to depend on the square of nD, and it is consistent with Dawson’s conclusion. With the power fitting, it is found that the exponent for the quadratic profile is closer to 2 than that for the square profile.