主要讨论了Cox-Ross-Rubinstein（CRR）模型和广义的CRR模型，并研究了如何基于CRR模型和广义的CRR模型利用Monte Carlo模拟计算资产价格以及期权价值。

【doc】 PbTiO3薄膜生长的Monte Carlo模拟Ⅰ：模型与算法.doc

运用Monte Carlo方法对线型高分子链格点模型进行模拟，研究了高分子链的动态行为随体系内组分间相互作用能的变化情况。讨论了随着溶剂分子与高分子链段间的相互作用能εPS与高分子链段间的相互作用能εPP的增大，链动力学物理量gcm（t），gm （t），ge（t）和gr（t）的变化情况。结果显示，εPP对高分子链的动态行为较εPS影响较小，但体系缠结的特征不完全取决于εPS和εPP的变化，主要受体系链长度的影响。

这组 MATLAB 例程模拟具有周期性边界条件的三维硬球堆积。 它初始化一个均匀晶格，然​​后使用 Metropolis 算法对粒子位置进行多次迭代退火。 然后计算 $g_2(r)$，或成对径向分布函数，在集合中的许多中心上取平均值。 对于大的初始粒子密度，结果 $g_2$ 与 Percus-Yevick 近似给出的解析近似匹配，给定足够高的粒子密度的任意精度。

该 Ising 模型用于通过应用 Metropolis 算法-蒙特卡洛方法模拟磁系统（正、负或随机自旋）。 运行主文件，输入晶格大小（最好是 100），然后选择一个输入自旋作为初始配置。 设置了两个不同的温度（T=2.0 和 T=2.5）。 例如，在 T=2（低温）下使用选定的正自旋进行初始化，大多数自旋是黑色的，这是因为翻转自旋的机会很小，并且材料具有铁磁性。 对于 T=2.5 的情况，由于温度较高，它会产生自旋翻转的趋势。 因此，材料失去磁化。 自旋似乎没有对齐，由此产生的配置似乎是随机无序的。 那是因为顺磁行为。 模拟的下一部分是可观测值计算：平均磁化、平均能量、平均磁化率和比热。 为了计算平均能量和磁化强度，我们必须找到能量和磁化值变化很小的时间依赖性（能量和磁化强度随时间增加变化很小的时间）因此，我们选择精度 p 并检查间隔（应满足精度的时间步数）。 这些间隔应取决于初始配

本文件是用Mathematica代码编写的Monte carlo模拟。 包含均匀分布，指数分布，以及均匀分布在圆和球面上的点。

IsingModel-蒙特卡洛 日语中的Ising模型的蒙特卡洛（MC）仿真简介。 Ising2d_MC 用日语记录Ising模型和MC模拟（Metropolis方法和热浴方法）。 使用Python的各向同性Ising mocdel的MC模拟代码（计算磁化强度，比热，Binder累积量和相关函数）。 Python代码中的注释为英文。 Ising2d_各向异性 使用Python（计算磁化强度）的各向异性Ising模型的MC模拟代码。

用MonteCairo方法通过matlab对薄膜外延生长进行了计算机模拟。模型中引入Morse势描述粒子间的相互作用，考虑粒子的沉积、吸附粒子的扩散和蒸发3个过程。研究了粒子间相互作用范围α和允许粒子行走的最大步数对薄膜生长形貌的影响。结果表明：在不同的6c值下，随粒子行走步数的增加，薄膜的生长都经历了从分散、分形、混合到团聚的过程；在同一最大步数下，α值越大，薄膜越易趋向于分散生长。

本资源是2019年全国大学数学竞赛C题省一等奖论文内含R语言代码。有用的伙伴下载学习交流。 声明：此论文只供自己学习使用，内容切勿用于商用。

使用 Geant4 工具包对闪烁粒子检测器进行蒙特卡罗建模 正在对闪烁中子探测器进行大量研究，作为全球短缺的氦 3 的替代品。 这个特殊的应用是设计一个超冷中子探测器，用于在超流氦中以 0.5 K 运行的中子 EDM 实验。 这项工作是我一部分。 Monte Carlo 模拟在原型探测器的设计和优化中起着至关重要的作用。 闪烁体是一种吸收高能粒子并发射光子的材料，然后可以使用标准光检测器（例如 PMT 或 CCD）对其进行计数。 为了最大限度地提高检测效率，尤其是对低能粒子的检测效率，必须最大限度地提高收集到的光量。 光学模拟包括产生闪烁光，光子的不同体积内的跟踪，并用所述光收集系统和转换成电子信号的相互作用。 要回答的一个重要问题涉及估计的光产额值（测量的光子与产生的光子的百分比），这是优化探测器效率和性能的关键参数。 这将取决于闪烁体的选择、几何形状、包裹探测器的反射材料的选择、机