| 查看: 1501 | 回复: 7 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
想请教一个有关二面角的问题！ 已有1人参与
|
|||
|
我现在使用tinker程序中的xyzint程序，将笛卡尔坐标(xyz coor.)转化为内坐标(internal coor.)的形式； 然后再使用其中的intxyz程序将内坐标转化回笛卡尔坐标形式，但是却得不到原来的构型！ 关键的问题出现在转化二面角上，比如：二面角ABCD 程序中是先求出angle(ABC)所在平面的法向量；然后再求出angle(BCD)所在平面的法向量，最后求出两个平面的法向量，最终求出这两个平面的夹角！ 我的问题是： 1. 当这两个平面重合的情况下，是应该没有二面角，还是0度或者180度？ 因为这种情况下，程序中是处理成，二面角直接等于angle(ABD)，这个平面角了，这样的处理不太理解！ 2. 如果ABC三点或者BCD三点在同一条直线上，这样的二面角又该怎么处理？ 3. 同样，如果两个平面重合的话，它们各自对应的法向量应该是什么关系呢？ 我是学化学的，二面角以及向量这块儿，很是欠缺，求大侠帮助！ |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有9人回复
-
找导师
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有11人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教一个 in-plane 和 out-of-plane 问题。
已经有10人回复
- 请教一个分子量分布的问题
已经有13人回复
- 请教一个问题，谢谢解答！
已经有11人回复
- 请教一个关于皂化酸化的问题
已经有7人回复
- 请教一个关于三氟甲基化的问题，懂得请进
已经有11人回复
- 请教一个电化学腐蚀的问题
已经有15人回复
- 请教一个用vesta建模的问题，谢谢！
已经有4人回复
- 请教一个关于香农公式的问题！急求！
已经有4人回复
- 各位师兄师姐，我想向大家请教下关于研究方向的一个选择问题
已经有9人回复
- 又是一个问题，关于二面角的？？？
已经有6人回复
- 晶体E编辑意见，急求怎么解决
已经有13人回复
- 请教一个手性物质的旋光性问题
已经有12人回复
- 请教一个关于origin模拟的问题，紧急，谢谢！
已经有15人回复
- 我是新手想请教一个关于细胞铺板的问题
已经有4人回复
- 请教一个plasmonics投稿的问题
已经有6人回复
- 【求助】有关酸度的问题想请教一下
已经有17人回复
- 【请教】关于二面角
已经有6人回复
- 【求助】请教一个关于解方程的问题
已经有8人回复
- 【求助/交流】请教一个NCBI上注释的问题
已经有5人回复
- 【求助】请问，怎么计算两个环的二面角
已经有7人回复
- 【请教】请教一个六方晶系的问题
已经有12人回复
- 【请教】请教一个关于图片精度的问题
已经有6人回复
polypro
木虫 (正式写手)
- 应助: 197 (高中生)
- 金币: 3921.3
- 散金: 160
- 红花: 22
- 帖子: 512
- 在线: 62.6小时
- 虫号: 1809438
- 注册: 2012-05-11
- 性别: GG
- 专业: 理论和计算化学
|
极坐标转换 !this sub to transfer Cartesian coordinates (x,y,z) with Polar coordinates (rou,theta,phi) (radius) !control number 1 for (x,y,z)=>(rou,theta,phi) and -1 is reverse !x=rou*sin(phi)*cos(theta) !x=rou*sin(phi)*sin(theta) !z=rou*cos(phi) subroutine trans(a,b,c,iflag) parameter(pi=3.1415926,h=0.01745329252) !pi/180 real rou,theta,phi real a,b,c integer iflag !1 from Cartesan to Polar, -1 is reverse (degree) !print*,"a b cini",a,b,c,iflag if(iflag==1)then !Cartesan to Polar rou=sqrt(a*a+b*b+c*c) if(rou.eq.0)then a=0 b=0 c=0 goto 80 endif phi=acos(c/rou) if(abs(rou-c).le.1e-3.or.abs(b).le.1e-3)then theta=0 !point on Z aixs or Y axis, thets==0 fixed! else if(abs(a).le.1e-3)then theta=pi/2.0 else sinphi=sqrt(1-c*c/rou/rou) theta=acos(a/rou/sinphi) endif endif if(b.lt.0)theta=-theta+pi !all is in degree a=rou b=theta/h c=phi/h !in radius else !polar to Cartesan rou=a if(rou.eq.0)then a=0 b=0 c=0 goto 80 endif theta=b*h phi=c*h c=rou*cos(phi) !z b=rou*sin(phi)*sin(theta) !y a=rou*sin(phi)*cos(theta) !x endif 80 continue !print*,"a b c",a,b,c,iflag END subroutine trans |
7楼2014-09-28 15:07:01
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3891
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24217
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
2楼2014-09-26 15:45:37
polypro
木虫 (正式写手)
- 应助: 197 (高中生)
- 金币: 3921.3
- 散金: 160
- 红花: 22
- 帖子: 512
- 在线: 62.6小时
- 虫号: 1809438
- 注册: 2012-05-11
- 性别: GG
- 专业: 理论和计算化学
3楼2014-09-28 08:28:25
4楼2014-09-28 11:16:12
20