10KNpt/exp_resumeはだめです fix extra all print 1 "$i $s $d $g $l $f $t $v $e $p" file outp.txt fix nveall all nve minimize 1.0e-18 1.0e-18 100000 100000 unfix nveall fix npt1 all npt temp ${npt_temp} ${npt_temp} 1 iso ${atmos} ${atmos} 1 …

#fix dt/reset のXmax決める用 temp = 0 timing = 0 for i in range(M-1): if abs(x_detec1[i+1] - x_detec1[i])>temp: temp = abs(x_detec1[i+1] - x_detec[i]) timing = i print(temp) print(timing)

2) 10K以上の弾性波伝播で，detectorが波を受けるまでの間にゆっくりと座標を変化させている（右にずれている）。弾性波を流さずに長いタイムステップnveで放置したらどうなるか見て，この現象についても考える。10Knve_doNothing(nve_relax)など →前も書い…

前回の問題点 20K以上の昇温で，セルサイズは収束しているように見えるが，温度変化がおかしい。より長いタイムステップをかけて移ろいを見ないといけない。20Kexp_longなど →15000 temp/rescale, 15000 relaxで以下のグラフ。 とりあえず1K前後の温度変動は…

20K以上の昇温で，セルサイズは収束しているように見えるが，温度変化がおかしい。より長いタイムステップをかけて移ろいを見ないといけない。20Kexp_longなど 10K以上の弾性波伝播で，detectorが波を受けるまでの間にゆっくりと座標を変化させている（右に…

M1で初飲み会した。みんな研究モチベないよね～～という話。 tabelog.com 店員さん優しかった。 LAMMPSのif文でやられるのもうn回（nは十分大きい自然数）なので，いい加減メモ。 docs.lammps.org LAMMPSのif構文は論理演算子は使えても，（加減乗除などの）…