[发明专利]模拟方法及程序

权利要求书

1.一种模拟方法，其特征在于，具有：

(a)工序，在考虑包含N个粒子、各粒子质量为m、且能够以表示对粒子间距离的依赖性的无因次函数f与相互作用系数ε之积εf来表示粒子间的相互作用势能的粒子系统S时，使用配置有该粒子系统S的空间的维数d确定大于1的第1重正化因子α、0以上d以下的第2重正化因子γ及0以上的第3重正化因子δ，并通过转换式N′＝N/α^d求出重正化的粒子数N′，通过转换式m′＝mα^δ/α^γ求出重正化的粒子质量m′，通过转换式ε′＝εα^γ求出重正化的相互作用系数ε′；及

(b)工序，对包含重正化的粒子数为N′个的粒子、各粒子具有重正化的粒子质量m′，且由所述无因次函数f和重正化的相互作用系数ε′之积ε′f表示粒子间的相互作用势能的粒子系统S′执行分子动力学计算。

2.如权利要求1所述的模拟方法，其特征在于，还具有(c)工序，所述(c)工序如下：在将所述工序(b)中用分子动力学计算求出的所述粒子系统S′中所含的某一粒子的位置矢量表示为q′，将所述工序(b)中用分子动力学计算求出的所述粒子系统S′中所含的某一粒子的运动量矢量表示为p′，将所述工序(b)中用分子动力学计算求出的所述粒子系统S′中所含的某一粒子的速度矢量表示为v′，将所述工序(b)中的分子动力学计算中的所述粒子系统S′中所含的某一粒子的某一时间间隔表示为t′时，使用在所述工序(a)中确定的所述第1重正化因子α、第2重正化因子γ及第3重正化因子δ中的至少一个，执行通过转换式q＝q′α求出所述粒子系统S中所含的某一粒子的位置矢量q的计算、通过转换式p＝p′/αδ/2求出所述粒子系统S中所含的某一粒子的运动量矢量p的计算、通过转换式v＝v′α^(-γ+δ/2)求出所述粒子系统S中所含的某一粒子的速度矢量v的计算、及通过转换式t＝t′α^(1+γ-δ/2)求出假设对所述粒子系统S执行的分子动力学计算中的某一时间间隔t的计算中的至少1个计算。

3.如权利要求1或2所述的模拟方法，其特征在于，

所述第2重正化因子γ为0。

4.如权利要求1～3中任一项所述的模拟方法，其特征在于，

在将n设为正整数时，所述第1重正化因子为2ⁿ。

5.如权利要求1～4中任一项所述的模拟方法，其特征在于，

所述维数d为3。

6.一种模拟方法，其特征在于，具有：

(a)工序，在考虑包含N个粒子、各粒子质量为m、且能够以表示对粒子间距离的依赖性的无因次函数f与相互作用系数ε之积εf来表示粒子间的相互作用势能的粒子系统S时，使用配置有该粒子系统S的空间的维数d确定大于1的第1重正化因子α和0以上的第3重正化因子δ，并通过转换式N′＝N/α^d求出重正化的粒子数N′，通过转换式m′＝mα^δ求出重正化的粒子质量m′，及

(b)工序，对包含重正化的粒子数为N′个的粒子、各粒子具有重正化的粒子质量m′、且由所述无因次函数f与所述相互作用系数ε之积εf表示粒子间的相互作用势能的粒子系统S′执行分子动力学计算。

7.如权利要求6所述的模拟方法，其特征在于，还具有(c)工序，所述(c)工序如下：在将所述工序(b)中用分子动力学计算求出的所述粒子系统S′中所含的某一粒子的位置矢量表示为q′，将所述工序(b)中用分子动力学计算求出的所述粒子系统S′中所含的某一粒子的运动量矢量表示为p′，将所述工序(b)中用分子动力学计算求出的所述粒子系统S′中所含的某一粒子的速度矢量表示为v′时，使用在所述工序(a)中确定的所述第1重正化因子α和所述第3重正化因子δ中的至少一方，执行通过转换式q＝q′α求出所述粒子系统S中所含的某一粒子的位置矢量q的计算、通过转换式p＝p′/α^δ/2求出所述粒子系统S中所含的某一粒子的运动量矢量p的计算、及通过转换式v＝v′α^δ/2求出所述粒子系统S中所含的某一粒子的速度矢量v的计算中的至少1个计算。