【手続補正書】
【提出日】令和1年6月11日(2019.6.11)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
くりこみ回数に依存するくりこみ因子λに基づいて、各々が一次元状に繋がった複数の
モノマー粒子からなる複数のポリマーで構成されたシミュレーション対象の粒子系Ｓに対
してくりこみ変換処理を行う工程と、
くりこまれた粒子系Ｓ’について分子動力学計算を実行することにより、くりこまれた
粒子系Ｓ’のモノマー粒子の位置ベクトル及び運動量ベクトルを求める工程と
を有し、
粒子系Ｓ及び粒子系Ｓ’のモノマー間の距離をｒで表し、空間の次元数をｄで表し、
粒子系Ｓのモノマーの質量をｍで表し、１つのポリマーを構成しているモノマー粒子の
数をＮｍで表し、ポリマーの数をＮｐで表し、モノマー粒子の座標をｑで表し、モノマー
粒子間の相互作用ポテンシャルφ（ｒ）として、
任意のモノマー粒子間においては、
【数１】
2017182231.tif 000001
が適用され、ポテンシャルＵ０（ｒ）は、エネルギの次元を持つパラメータεと、長さの
次元を持つパラメータσと、無次元関数ｆを用いて
【数２】
2017182231.tif 000002
と表すことができ、
同一ポリマー内の隣接するモノマー粒子間においては、パラメータとしてポテンシャル
Ｕ０（ｒ）で用いられたパラメータε、σを持つ有限伸長非線形弾性ポテンシャルＵｃｈ
（ｒ：ε，σ）がポテンシャルＵ０（ｒ）に追加されて
【数３】
2017182231.tif 000003
が適用され、
粒子系Ｓ’のモノマー粒子の質量をｍ’で表し、１つのポリマーを構成しているモノマ
ー粒子の数をＮｍ’で表し、ポリマーの数をＮｐ’で表し、モノマー粒子の座標をｑ’で
表したとき、下記の変換則
【数４】
2017182231.tif 000004
を適用し、
粒子系ＳにおけるポテンシャルＵ０（ｒ）、Ｕｃｈ（ｒ：ε，σ）に対応する粒子系Ｓ
’におけるポテンシャルＵ０’（ｒ）、Ｕｃｈ’（ｒ：ε，σ）として、
【数５】
2017182231.tif 000005
を用いて、粒子系Ｓ’について前記分子動力学計算を実行するシミュレーション方法。
【請求項２】
粒子系ＳにおけるポテンシャルＵ０（ｒ）は、レナード・ジョーンズ型ポテンシャルの
形で表される請求項１に記載のシミュレーション方法。
【請求項３】
有限伸長非線形弾性ポテンシャルＵｃｈ（ε，σ，ｒ）は、パラメータε、σに依存す
るパラメータｋ及びＲ０を用いて、
【数６】
2017182231.tif 000006
と表される請求項１または２に記載のシミュレーション方法。
【請求項４】
請求項１乃至３のいずれか１項に記載のシミュレーション方法をコンピュータで実行す
るためのコンピュータプログラム。
【請求項５】
請求項４に記載のコンピュータプログラムが、コンピュータで読み取り可能に記録され
た記録媒体。
【請求項６】
各々が一次元状に繋がった複数のモノマー粒子からなる複数のポリマーで構成されたシ
ミュレーション対象の粒子系Ｓを特徴づけるパラメータである空間の次元数ｄ、モノマー
粒子の質量ｍ、１つのポリマーを構成しているモノマー粒子の数Ｎｍ、ポリマーの数Ｎｐ
の値、及び粒子系Ｓの初期条件が入力される入力装置と、
前記入力装置から入力されたデータに基づいてシミュレーション処理を実行する処理ユ
ニットと
を有し、
前記処理ユニットは、
粒子系Ｓに対して、くりこみ回数に依存するくりこみ因子λに基づいてくりこみ変換を
行う処理と、
くりこまれた粒子系Ｓ’について分子動力学計算を実行することにより、くりこまれた
粒子系Ｓ’のモノマー粒子の位置及び運動量を求める処理と
を実行し、
前記処理ユニットは、
粒子系Ｓ及び粒子系Ｓ’のモノマー間の距離をｒで表し、
粒子系Ｓのモノマー粒子の座標をｑで表し、モノマー粒子間の相互作用ポテンシャルφ
（ｒ）として、
任意のモノマー粒子間においては、
【数７】
2017182231.tif 000007
が適用され、ポテンシャルＵ０（ｒ）は、エネルギの次元を持つパラメータεと、長さの
次元を持つパラメータσと、無次元関数ｆを用いて
【数８】
2017182231.tif 000008
と表すことができ、
同一ポリマー内の隣接するモノマー粒子間においては、パラメータとしてポテンシャル
Ｕ０（ｒ）で用いられたパラメータε、σを持つ有限伸長非線形弾性ポテンシャルＵｃｈ
（ｒ：ε，σ）がポテンシャルＵ０（ｒ）に追加されて
【数９】
2017182231.tif 000009
が適用され、
粒子系Ｓ’のモノマー粒子の質量をｍ’で表し、１つのポリマーを構成しているモノマ
ー粒子の数をＮｍ’で表し、ポリマーの数をＮｐ’で表し、モノマー粒子の座標をｑ’で
表したとき、下記の変換則
【数１０】
2017182231.tif 000010
を適用し、
粒子系ＳにおけるポテンシャルＵ０（ｒ）、Ｕｃｈ（ε，σ，ｒ）に対応する粒子系Ｓ
’におけるポテンシャルＵ０’（ｒ）、Ｕｃｈ’（ε，σ，ｒ）として、
【数１１】
2017182231.tif 000011
を用いて、粒子系Ｓ’について前記分子動力学計算を実行するシミュレーション装置。
【請求項７】
くりこみ回数に依存するくりこみ因子λに基づいて、各々が一次元状に繋がった複数の
モノマー粒子からなる複数のポリマーで構成されたシミュレーション対象の粒子系Ｓに対
してくりこみ変換処理を行う工程と、
くりこまれた粒子系Ｓ’について分子動力学計算を実行することにより、くりこまれた
粒子系Ｓ’のモノマー粒子の位置ベクトル及び運動量ベクトルを求める工程と
を有し、
粒子系Ｓ及び粒子系Ｓ’のモノマー間の距離をｒで表し、空間の次元数をｄで表し、
粒子系Ｓのモノマーの質量をｍで表し、１つのポリマーを構成しているモノマー粒子の
数をＮｍで表し、ポリマーの数をＮｐで表し、モノマー粒子の座標をｑで表し、モノマー
粒子間の相互作用ポテンシャルφ（ｒ）として、
任意のモノマー粒子間においては、
【数１２】
2017182231.tif 000012
が適用され、
同一ポリマー内の隣接するモノマー粒子間においては、パラメータとしてポテンシャル
Ｕ０（ｒ）で用いられたパラメータε、σを持つ有限伸長非線形弾性ポテンシャルＵｃｈ
（ｒ：ε，σ）がポテンシャルＵ０（ｒ）に追加されて
【数１３】
2017182231.tif 000013
が適用され、
ポテンシャルφ（ｒ）は、エネルギの次元を持つパラメータεと、長さの次元を持つパ
ラメータσと、無次元関数ｆを用いて
【数１４】
2017182231.tif 000014
と表すことができ、
粒子系Ｓ’のモノマー粒子の質量をｍ’で表し、１つのポリマーを構成しているモノマ
ー粒子の数をＮｍ’で表し、ポリマーの数をＮｐ’で表し、モノマー粒子の座標をｑ’で
表したとき、下記の変換則
【数１５】
2017182231.tif 000015
を適用し、
粒子系Ｓにおけるポテンシャルφ（ｒ）に対応する粒子系Ｓ’におけるポテンシャルφ
’（ｒ）として、
【数１６】
2017182231.tif 000016
を用いて、粒子系Ｓ’について前記分子動力学計算を実行するシミュレーション方法。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００１１
【補正方法】変更
【補正の内容】
【００１１】
本発明の一観点によると、
くりこみ回数に依存するくりこみ因子λに基づいて、各々が一次元状に繋がった複数の
モノマー粒子からなる複数のポリマーで構成されたシミュレーション対象の粒子系Ｓに対
してくりこみ変換処理を行う工程と、
くりこまれた粒子系Ｓ’について分子動力学計算を実行することにより、くりこまれた
粒子系Ｓ’のモノマー粒子の位置ベクトル及び運動量ベクトルを求める工程と
を有し、
粒子系Ｓ及び粒子系Ｓ’のモノマー間の距離をｒで表し、空間の次元数をｄで表し、
粒子系Ｓのモノマーの質量をｍで表し、１つのポリマーを構成しているモノマー粒子の
数をＮｍで表し、ポリマーの数をＮｐで表し、モノマー粒子の座標をｑで表し、モノマー
粒子間の相互作用ポテンシャルφ（ｒ）として、
任意のモノマー粒子間においては、
【数１】
2017182231.tif 000017
が適用され、ポテンシャルＵ０（ｒ）は、エネルギの次元を持つパラメータεと、長さの
次元を持つパラメータσと、無次元関数ｆを用いて
【数２】
2017182231.tif 000018
と表すことができ、
同一ポリマー内の隣接するモノマー粒子間においては、パラメータとしてポテンシャル
Ｕ０（ｒ）で用いられたパラメータε、σを持つ有限伸長非線形弾性ポテンシャルＵｃｈ
（ｒ：ε，σ）がポテンシャルＵ０（ｒ）に追加されて
【数３】
2017182231.tif 000019
が適用され、
粒子系Ｓ’のモノマー粒子の質量をｍ’で表し、１つのポリマーを構成しているモノマ
ー粒子の数をＮｍ’で表し、ポリマーの数をＮｐ’で表し、モノマー粒子の座標をｑ’で
表したとき、下記の変換則
【数４】
2017182231.tif 000020
を適用し、
粒子系ＳにおけるポテンシャルＵ０（ｒ）、Ｕｃｈ（ｒ：ε，σ）に対応する粒子系Ｓ
’におけるポテンシャルＵ０’（ｒ）、Ｕｃｈ’（ｒ：ε，σ）として、
【数５】
2017182231.tif 000021
を用いて、粒子系Ｓ’について前記分子動力学計算を実行するシミュレーション方法が提
供される。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００３４
【補正方法】変更
【補正の内容】
【００３４】
くりこまれた粒子系Ｓ’の各粒子の位置ベクトルｑ’と、シミュレーション対象の粒子
系Ｓの各粒子の位置ベクトルｑとの対応関係は、式（１３）に示したとおりである。くりこまれた粒子系Ｓ’における運動量ベクトルｐ’と、シミュレーション対象の粒子系Ｓにおける運動量ベクトルｐとは、下記の関係を有する。
【数２０】
2017182231.tif 000022