| 查看: 1080 | 回复: 13 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】求一模拟程序（fortran）
|
|||
|
本人初次编程解决实际问题，感觉力不从心！！请求各位高手帮忙！ 问题如下：用fortran语言编程即时计算和展现十个氢原子在模拟盒中的运动情况。初速度满足麦克斯韦分布，采取硬球完全弹性碰撞 [ Last edited by lei0736 on 2009-11-25 at 11:58 ] |
回复此楼» 猜你喜欢
求调剂
已经有23人回复
-
301求调剂
已经有15人回复
-
304求调剂（085602，过四级，一志愿985）
已经有17人回复
-
302分求调剂 一志愿安徽大学085601
已经有12人回复
-
288环境专硕,求调材料方向
已经有23人回复
-
环境专硕调剂
已经有6人回复
-
22408 调剂材料
已经有6人回复
-
285求调剂
已经有12人回复
-
求调剂
已经有6人回复
-
085600材料与化工301分求调剂院校
已经有19人回复
zyy816
至尊木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 9698.1
- 散金: 7
- 红花: 2
- 帖子: 1987
- 在线: 320.1小时
- 虫号: 162947
- 注册: 2006-01-11
- 专业: 理论和计算化学
2楼2009-04-12 15:18:40
3楼2009-04-12 15:26:12
4楼2009-04-12 15:32:44
bay__gulf
金虫 (著名写手)
刘苏州
- 模拟EPI: 8
- 应助: 9 (幼儿园)
- 贵宾: 4.85
- 金币: 2332.8
- 红花: 1
- 帖子: 1344
- 在线: 271小时
- 虫号: 592012
- 注册: 2008-09-03
- 专业: 理论和计算化学
- 管辖: 分子模拟
★ ★ ★ ★ ★
jianchaoyv2651(金币+3,VIP+0): 4-12 16:02
lei0736(金币+2,VIP+0):谢谢 4-20 17:33
jianchaoyv2651(金币+3,VIP+0): 4-12 16:02
lei0736(金币+2,VIP+0):谢谢 4-20 17:33
|
http://www.mdbbs.org/viewthread.php?tid=6423 这里有个模拟256个Ar原子，LJ势的供学习用的frotran代码以及一些解释。 H原子因为强的量子力学效应不适合MD |
5楼2009-04-12 15:38:05
jianchaoyv2651(金币+11,VIP+0): 4-13 11:39
lei0736(金币+0,VIP+0):马甲行为 自己给自己转金币 请解释！！！ 5-7 11:58
lei0736(金币-11,VIP+0):请解释 5-9 13:32
lei0736(金币+0,VIP+0):马甲行为 自己给自己转金币 请解释！！！ 5-7 11:58
lei0736(金币-11,VIP+0):请解释 5-9 13:32
| 问题已解决！！ |
6楼2009-04-13 11:36:01
7楼2009-05-07 10:28:32
8楼2009-05-07 14:02:05
bay__gulf
金虫 (著名写手)
刘苏州
- 模拟EPI: 8
- 应助: 9 (幼儿园)
- 贵宾: 4.85
- 金币: 2332.8
- 红花: 1
- 帖子: 1344
- 在线: 271小时
- 虫号: 592012
- 注册: 2008-09-03
- 专业: 理论和计算化学
- 管辖: 分子模拟
9楼2009-05-07 14:11:16
|
请问我是按所给程序运行的，所得结果怎样转化为我在其他文献上看到此程序所得的总能量随步数演化的图形、模拟中动能的分布图以及其他一些物理量的值等，请问应该怎样去操作？ 下面是我运行结果的前面的一部分，后面是从100到3000部分我未写出 MOLECULAR DYNAMICS SIMULATION PROGRAM NUMBER OF PARTICLE IS 256 SIDE LENGTH OF THE BOX IS 6.752840 CUT OFF IS 2.500000 REDUCED TEMPERATURE IS 0.7220000 BASIC TIME STEP IS 6.4000003E-02 TOTAL LINEAR MOMENTUM IN X DIRECTION IS 6.7786448E-02 TOTAL LINEAR MOMENTUM IN Y DIRECTION IS -6.6548891E-02 TOTAL LINEAR MOMENTUM IN Z DIRECTION IS 0.1201641 VELOCITY ADJUSTMENT TEMPERATURE BEFORE SCALING IS 45.41020 SCALE FACTOR IS 0.1260933 2 270.799927 -1699.370728 -1428.570801 4 252.771790 -1681.374146 -1428.602295 6 218.253845 -1646.222656 -1427.968750 8 171.139252 -1598.697998 -1427.558716 10 132.214249 -1559.054932 -1426.840698 12 116.545746 -1543.839722 -1427.293945 14 119.052315 -1546.832153 -1427.779785 16 127.827530 -1555.922119 -1428.094604 18 134.569000 -1562.949707 -1428.380737 20 137.285446 -1565.422729 -1428.137329 22 137.336853 -1565.384644 -1428.047852 24 136.563614 -1564.932617 -1428.369019 26 134.083801 -1562.593628 -1428.509766 28 129.720490 -1558.355103 -1428.634644 30 126.987198 -1555.609863 -1428.622681 32 127.910660 -1556.748169 -1428.837524 34 131.486847 -1560.345459 -1428.858643 36 135.715729 -1564.671997 -1428.956299 38 138.538544 -1567.548828 -1429.010254 40 139.642487 -1568.478882 -1428.836426 42 140.040283 -1568.885742 -1428.845459 44 140.190887 -1569.073975 -1428.883057 46 139.712982 -1568.837280 -1429.124268 48 138.952988 -1568.239868 -1429.286865 VELOCITY ADJUSTMENT TEMPERATURE BEFORE SCALING IS 0.3634446 SCALE FACTOR IS 1.409449 50 276.164978 -1568.357178 -1292.192139 52 271.862701 -1564.260986 -1292.398315 54 259.602783 -1552.116089 -1292.513306 56 242.570480 -1535.268799 -1292.698364 58 225.266388 -1518.148560 -1292.882202 60 211.759369 -1504.485107 -1292.725708 62 205.603088 -1498.357300 -1292.754150 64 209.019333 -1501.927979 -1292.908691 66 213.457825 -1506.205566 -1292.747803 68 215.538666 -1508.412109 -1292.873413 70 217.871109 -1510.706299 -1292.835205 72 223.228973 -1516.197266 -1292.968262 74 229.282425 -1522.213135 -1292.930664 76 230.661316 -1523.383911 -1292.722656 78 227.155609 -1519.937866 -1292.782227 80 221.844513 -1514.623657 -1292.779175 82 215.889771 -1508.402832 -1292.513062 84 211.333252 -1503.536133 -1292.202881 86 211.333344 -1503.780518 -1292.447144 88 216.985687 -1509.352051 -1292.366333 90 227.358459 -1519.550293 -1292.191895 92 235.516296 -1527.578247 -1292.062012 94 236.259933 -1528.119873 -1291.859985 96 232.384796 -1524.102417 -1291.717651 98 227.971985 -1519.898438 -1291.926514 VELOCITY ADJUSTMENT TEMPERATURE BEFORE SCALING IS 0.5878155 SCALE FACTOR IS 1.108276 100 276.164978 -1516.690308 -1240.525391 |
10楼2009-05-07 16:42:21
