蒙特卡罗软件简介

蒙特卡罗简介蒙特卡罗是一款模拟实验室中粒子碰撞的软件，可以代替实验操作节约成本。我们要做的是正确写好输入卡（也就是构造一个和实验室相近的实验条件）。首先建立一个新建文本文档，然后重命名，把文本文档的扩展名去掉。这样文件就建立好了，点开写输入卡就可以了。下面介绍输入卡如何书写和如何运行。一、mcnp输入文件的编写输入文件包括几何区域（cell）定义块、面（surface）定义块、物质定义块、粒子源定义块。1.cell定义块区域号物质号物质密度区域定义（1）区域号从1开始到最后（2）物质序号从1开始到最后，但是真空用0代替，跳过真空部分继续。（3）物质密度：单位是g/cm3（4）区域定义：向左（—），向右（+），向下（—），向上（+）例子：如下图如何编写cell卡（图和卡对应）11-5.323(14-131-2):(15-14-2)$区域一GE晶体22-5.323(2-3-1315):(-316-15)$区域二GE死层33-8.9-114-13$区域三铜棒40(-413-12):(-417-16):(3-416-13)$区域四真空54-2.70(-512-11):(4-517-12):(3-5-1719)$区域五铝层65-6.5-17-318$区域六铍层76-0.0012-319-18$区域七空气87-2.2(-2021-7):(-2223-7):(6-7-2122)$区域八[CF2CF2]n109-11.3(-927-1011):(8-9-1125):(28-9-2526)$区域十铅109-11.31110-0.0012-24#1#2#3#4#5#6#7#8#9#10$区域十空气12024 $1110-0.0012-24#1#2#3#4#5#6#7#8#9#10$区域十空气12024 $区域十一真空2.面（surface）定义块几何区域定义完还要把所需的面进行定义CZ：代表以z为对称轴的柱面PZ：代表以Z为对称轴的平面我们需要用到的只有柱面和平面，其他的要在仪器手册上查找例子（还是上面的图形）1cz0.52cz2.73cz2.84cz3.935cz46cz3.57cz3.758cz209cz29.710Pz19.5311pz13.5312pz13.463.物质定义块（1）MODE定义模式（定义模拟实验中用的是什么粒子）格式：MODEP（质子俘获）。在输入卡上就用MODEP。其他的要在手册上查询。（2）定义区域例子重要性（解决粒子在区域中是否传输问题）格式：IMP:P1110（模拟光子在区域1,2,3中传输，进入区域4便结束对其模拟）1代表重要也就是可以传输，0代表不重要也就是截止传输模拟注意：咱们前面对应的cell有N部分，那么1的个数就是（N-1）个，最后一个填写0。源的说明（对测量的放射源的性质进行说明）1)SDEF卡SDEF源的位置，能量，能谱，方向，类型，分布等等。这部分内容比较复杂，只做简单介绍。CEL=出射的几何快序号ERG=出射的能量POS=XYZ源的几何位置RAD=抽样位置离开POS或者离开AXS的径向距离，相当于圆柱形的半径。已乂丁二沿AXS离开AXS的距离。相当于圆柱的高。（2）源信息卡Sin：源信息SPn：源的概率L：光子能谱D：相应间隔概率（离散型概率）Sin和SPn可以在中另一个decay程序中查出。我已经整理出镭、钍、钾的衰变图纲（文件夹里有）。如需要其他的自己可以查找。但是注意SPn的总和为1•也就是概率归一化。例子：K粒子源定义（还是上面的彩图对应）SDEFCEL=9ERG=D1POS=00-3・5RAD=D2EXT=D3L1.460800SP1D103.5SP2-210-3.253.25SP3-210物质说明格式：MmZZZAAA核素分之比解释：M是固定格式不变；m是几何快中的物质序号，要和cell中的物质序号对应上。ZZZ是原子数；AAA是质量数。（1）如是单质M9820001铅单质（单质核素分之比为1）（2）如是混合物例子：有机玻璃成分C5H8O2（甲基丙烯酸甲酯）C分之比= 乞12 二0.612x5+8x1+16x2H分之比= 8x1 二0.0812x5+8x1+16x2二0.32或者用1-0.6-0.08=0.32O分之比= 二0.32或者用1-0.6-0.08=0.3212x5+8x1+16x2其他（无需改动）F8:P1e801e-517.932E-317.9331E-38190I1.8691ft8geb0.9748e-30.8415e-3NPS1e7P