24小时热门版块排行榜

>论坛更新日志 (3223)
>虫友互识 (264)
>文献求助 (209)
>导师招生 (144)
>硕博家园 (100)
>休闲灌水 (68)
>考博 (67)
>论文道贺祈福 (53)
>考研 (45)
>博后之家 (44)
>教师之家 (42)
>公派出国 (42)
>论文投稿 (40)
>基金申请 (35)
>招聘信息布告栏 (33)
>找工作 (28)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xzz1007

金虫 (小有名气)

应助: 8 (幼儿园)
金币: 300.4
红花: 2
帖子: 189
在线: 76.3小时
虫号: 2275671
注册: 2013-02-06
性别: MM
专业: 理论和计算化学

[求助] 想请教一个有关二面角的问题！已有1人参与

我现在使用tinker程序中的xyzint程序，将笛卡尔坐标(xyz coor.)转化为内坐标(internal coor.)的形式；
然后再使用其中的intxyz程序将内坐标转化回笛卡尔坐标形式，但是却得不到原来的构型！

关键的问题出现在转化二面角上，比如：二面角ABCD

程序中是先求出angle(ABC)所在平面的法向量；然后再求出angle(BCD)所在平面的法向量，最后求出两个平面的法向量，最终求出这两个平面的夹角！

我的问题是：
1. 当这两个平面重合的情况下，是应该没有二面角，还是0度或者180度？
因为这种情况下，程序中是处理成，二面角直接等于angle(ABD)，这个平面角了，这样的处理不太理解！

2. 如果ABC三点或者BCD三点在同一条直线上，这样的二面角又该怎么处理？

3. 同样，如果两个平面重合的话，它们各自对应的法向量应该是什么关系呢？

我是学化学的，二面角以及向量这块儿，很是欠缺，求大侠帮助！

回复此楼

» 猜你喜欢

论文投稿，期刊推荐已经有6人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有3人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有6人回复
孩子确诊有中度注意力缺陷已经有14人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

请教一个 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人回复

阳光总在风雨后，再坚持坚持吧.....

1楼2014-09-26 09:55:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

polypro

木虫 (正式写手)

应助: 197 (高中生)
金币: 3921.3
散金: 160
红花: 22
帖子: 512
在线: 62.6小时
虫号: 1809438
注册: 2012-05-11
性别: GG
专业: 理论和计算化学

引用回帖:

4楼: Originally posted by xzz1007 at 2014-09-28 11:16:12
可以发给我一份吗?我想对比一下！
我总觉得tinker中的程序在转化的时候有点问题！
留个邮箱吧！...

极坐标转换
!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