[发明专利]一种利用模拟分子动力学方法确定后装源驻留时间的方法

权利要求书

1.一种利用模拟分子动力学方法确定后装源驻留时间的方法，其特征在于，包括以下步骤；

(I)根据驻留点和剂量节制点的位置关系，采用式(1)算出每个驻留点在剂量节制点处的剂量：

dj=Σi=1Naijti---(1)]]>

(II)根据剂量节制点处的计算剂量与处方剂量的偏差和相邻驻留点的时间差构成评价函数：

O(t)=Σj=1Mωj(Σi=1Ndjiti-dj0)2+Σiα(ti+1-ti)2---(2)]]>

(III)根据式(2)对时间tk求偏导数，得到式(3)：

∂O(t)∂tk=2Σj=1MωjΣi=1Ndjitidjk-2Σj=1Mωjdj0djk+2αT---(3)]]>

(IV)将式(3)类似于一个分子动力学方程，分子间的相互作用力f_k表示如下：

-fk=∂O(t)∂tk=2Σj=1MωjΣi=1Ndjitidjk-2Σj=1Mωjdj0djk+2αT---(4)]]>

(V)采用下式(5)和(6)进行迭代，迭代过程中，通过不断改变τ的值，得出不同的驻留时间：

vk(τ+δτ)=vk(τ)+12mkδτ[fk(τ)+fk(τ+δτ)]---(5)]]>

tk(τ+δτ)=tk(τ)+δτvk(τ)+12mkδτ2fk(τ)---(6)]]>

(VI)当O(t)小于设定值时或者迭代次数达到设定值时退出迭代过程，得出最终的驻留时间解；

其中，M是剂量节制点的数目，N是驻留点的数目，ω_j是第j个剂量节制点的权重，t_i是第i个驻留点的驻留时间，是第i驻留点对第j个剂量节制点的剂量率，是第j个剂量节制点的节制剂量，α表示驻留时间平滑因子，v_k(τ)是第k个原子在τ时刻的速度，δτ是时间增量。