fork download
  1. PROGRAM LennardJones
  2. implicit none
  3.  
  4. ! 変数の定義
  5. INTEGER, PARAMETER :: Nplot = 101 ! プロット数
  6. REAL(8), PARAMETER :: R_ini = 0.95d0 ! 粒子間距離の最小値
  7. REAL(8), PARAMETER :: R_end = 2.0d0 ! 粒子間距離の最大値
  8. REAL(8), PARAMETER :: Eps = 1.0d0 ! Potential parameter 1
  9. REAL(8), PARAMETER :: Sigma = 1.0d0 ! Potential parameter 2
  10.  
  11. INTEGER i
  12. REAL(8) dR, R, pot
  13.  
  14. ! 距離の増分
  15. dR = (R_end - R_ini) / real(Nplot-1,8)
  16.  
  17.  
  18. ! print header
  19. write(*,'(a)') '#R, Potential'
  20.  
  21.  
  22. ! 粒子間距離をスキャンしながらLJポテンシャルを計算&出力
  23. do i= 1, Nplot
  24.  
  25. ! 距離の計算
  26. R = R_ini + dR * real(i-1,8)
  27.  
  28. ! ポテンシャルエネルギーの計算
  29. pot = 4.d0 * Eps * ( (Sigma/R)**12 - (Sigma/R)**6 )
  30.  
  31. ! 出力
  32. write(*,'(e12.6, e14.6)') R, pot
  33.  
  34. enddo ! i
  35.  
  36. END PROGRAM LennardJones
  37.  
Success #stdin #stdout 0.01s 5280KB
stdin
Standard input is empty
stdout
#R,  Potential
0.950000E+00  0.196097E+01
0.960500E+00  0.139351E+01
0.971000E+00  0.921669E+00
0.981500E+00  0.530452E+00
0.992000E+00  0.207244E+00
0.100250E+01 -0.585941E-01
0.101300E+01 -0.276039E+00
0.102350E+01 -0.452675E+00
0.103400E+01 -0.594912E+00
0.104450E+01 -0.708175E+00
0.105500E+01 -0.797057E+00
0.106550E+01 -0.865448E+00
0.107600E+01 -0.916645E+00
0.108650E+01 -0.953443E+00
0.109700E+01 -0.978214E+00
0.110750E+01 -0.992970E+00
0.111800E+01 -0.999415E+00
0.112850E+01 -0.998997E+00
0.113900E+01 -0.992941E+00
0.114950E+01 -0.982288E+00
0.116000E+01 -0.967918E+00
0.117050E+01 -0.950574E+00
0.118100E+01 -0.930887E+00
0.119150E+01 -0.909388E+00
0.120200E+01 -0.886523E+00
0.121250E+01 -0.862669E+00
0.122300E+01 -0.838141E+00
0.123350E+01 -0.813203E+00
0.124400E+01 -0.788073E+00
0.125450E+01 -0.762933E+00
0.126500E+01 -0.737933E+00
0.127550E+01 -0.713196E+00
0.128600E+01 -0.688820E+00
0.129650E+01 -0.664886E+00
0.130700E+01 -0.641456E+00
0.131750E+01 -0.618579E+00
0.132800E+01 -0.596293E+00
0.133850E+01 -0.574625E+00
0.134900E+01 -0.553592E+00
0.135950E+01 -0.533207E+00
0.137000E+01 -0.513475E+00
0.138050E+01 -0.494397E+00
0.139100E+01 -0.475970E+00
0.140150E+01 -0.458185E+00
0.141200E+01 -0.441036E+00
0.142250E+01 -0.424508E+00
0.143300E+01 -0.408591E+00
0.144350E+01 -0.393268E+00
0.145400E+01 -0.378524E+00
0.146450E+01 -0.364342E+00
0.147500E+01 -0.350707E+00
0.148550E+01 -0.337600E+00
0.149600E+01 -0.325004E+00
0.150650E+01 -0.312902E+00
0.151700E+01 -0.301276E+00
0.152750E+01 -0.290110E+00
0.153800E+01 -0.279385E+00
0.154850E+01 -0.269086E+00
0.155900E+01 -0.259197E+00
0.156950E+01 -0.249701E+00
0.158000E+01 -0.240583E+00
0.159050E+01 -0.231828E+00
0.160100E+01 -0.223422E+00
0.161150E+01 -0.215350E+00
0.162200E+01 -0.207599E+00
0.163250E+01 -0.200156E+00
0.164300E+01 -0.193008E+00
0.165350E+01 -0.186143E+00
0.166400E+01 -0.179550E+00
0.167450E+01 -0.173216E+00
0.168500E+01 -0.167132E+00
0.169550E+01 -0.161286E+00
0.170600E+01 -0.155669E+00
0.171650E+01 -0.150272E+00
0.172700E+01 -0.145084E+00
0.173750E+01 -0.140098E+00
0.174800E+01 -0.135305E+00
0.175850E+01 -0.130697E+00
0.176900E+01 -0.126266E+00
0.177950E+01 -0.122004E+00
0.179000E+01 -0.117906E+00
0.180050E+01 -0.113963E+00
0.181100E+01 -0.110169E+00
0.182150E+01 -0.106519E+00
0.183200E+01 -0.103007E+00
0.184250E+01 -0.996254E-01
0.185300E+01 -0.963706E-01
0.186350E+01 -0.932367E-01
0.187400E+01 -0.902190E-01
0.188450E+01 -0.873126E-01
0.189500E+01 -0.845130E-01
0.190550E+01 -0.818159E-01
0.191600E+01 -0.792171E-01
0.192650E+01 -0.767126E-01
0.193700E+01 -0.742986E-01
0.194750E+01 -0.719716E-01
0.195800E+01 -0.697280E-01
0.196850E+01 -0.675645E-01
0.197900E+01 -0.654779E-01
0.198950E+01 -0.634652E-01
0.200000E+01 -0.615234E-01