24小时热门版块排行榜    

查看: 5136  |  回复: 11

g_xq96

银虫 (小有名气)

[交流] 【求助】MS晶体建模里的分数坐标如何确定? 已有9人参与

我想自己手动输坐标建立模型!
在直角坐标系中确定了超胞中的几个原的子坐标,但不会转换成分数坐标。
感觉分数坐标应该很有用,想请各位高手介绍一下笛卡尔坐标如何换成分数坐标?
对不同的晶格常数(a,b,c,α,β,γ)分数坐标转换有何异同?
回复此楼

» 猜你喜欢

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

山路十八弯http://blog.sina.com.cn/u/1589327057
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzy870720z

荣誉版主 (文坛精英)

优秀版主优秀版主优秀版主优秀版主

★ ★
g_xq96(金币+1):谢谢参与
gavinliu7390(金币+1):谢谢交流! 2010-05-24 19:28:02
我记不太清楚公式
记得陈舜麟的《计算材料科学》的附录上有公式
用一个过渡矩阵联系直角坐标和内坐标
找那本书看看
博学、审问、慎思、明辨、笃学
2楼2010-05-24 17:00:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liufengbao100

木虫 (正式写手)

foreverhappy


g_xq96(金币+1):谢谢参与
百度一下就知道了
成长总是痛苦的,加油!
3楼2010-05-27 08:59:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主


g_xq96(金币+1):谢谢参与
自己手算一下也可以,中学的立体几何就可以了.
4楼2010-05-27 11:27:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flashpoint

金虫 (小有名气)

★ ★ ★
g_xq96(金币+1):谢谢参与
zzy870720z(金币+2):谢谢交流 2010-05-27 12:54:07
g_xq96(金币+3):QQ:705357207 2010-05-27 18:41:50
通常都是分数坐标,很少直接给直角坐标的,如果是三斜晶系,转换很复杂的,
你可以尝试这个公式:

%Define the crystal parameters
a=5;b=5;c=10;
alpha=90*pi/180;beta=90*pi/180;gamma=120*pi/180;
n2=(cos(alpha)-cos(gamma)*cos(beta))/sin(gamma);
n3=(sin(beta)*sin(beta)-n2^2)^0.5;
A=[a,   b*cos(gamma),    c*cos(beta);
   0,     b*sin(gamma),    c*n2;   
   0,    0,                c*n3];

这个矩阵是把分数坐标转变为直角坐标,反过来就可达到你的要求
引用回帖:
Originally posted by g_xq96 at 2010-05-24 16:35:23:
我想自己手动输坐标建立模型!
在直角坐标系中确定了超胞中的几个原的子坐标,但不会转换成分数坐标。
感觉分数坐标应该很有用,想请各位高手介绍一下笛卡尔坐标如何换成分数坐标?
对不同的晶格常数(a,b,c,α ...

5楼2010-05-27 11:48:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuanfang4123

铁虫 (初入文坛)


g_xq96(金币+1):谢谢参与
引用回帖:
Originally posted by flashpoint at 2010-05-27 11:48:05:
通常都是分数坐标,很少直接给直角坐标的,如果是三斜晶系,转换很复杂的,
你可以尝试这个公式:

%Define the crystal parameters
a=5;b=5;c=10;
alpha=90*pi/180;beta=90*pi/180;gamma=120*pi/180;
n2=( ...

这个是在input里面输入吗?
6楼2010-05-27 16:26:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

g_xq96

银虫 (小有名气)

这个不对呀!

引用回帖:
Originally posted by flashpoint at 2010-05-27 11:48:05:
通常都是分数坐标,很少直接给直角坐标的,如果是三斜晶系,转换很复杂的,
你可以尝试这个公式:

%Define the crystal parameters
a=5;b=5;c=10;
alpha=90*pi/180;beta=90*pi/180;gamma=120*pi/180;
n2=( ...

这里哪里是分数呀,分数坐标得到的是分数。
山路十八弯http://blog.sina.com.cn/u/1589327057
7楼2010-05-27 18:45:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)

★ ★ ★
g_xq96(金币+1):谢谢参与
cenwanglai(金币+2):谢谢,热情解答! 2010-05-28 08:29:22
勤劳一点,自己推导一下还是不难的。
哈哈,好就没做数学题了,貌似前两天有个贴子问过一个类似的情况,当时做过类似的推导:

设直角坐标系三个基矢为x,y和z,设点X在直角坐标系坐标为(x0,y0,z0)。
设晶胞的三个基矢为A,B,C,可以用直角坐标系表示为
A=a1*x+a2*y+a3*z, B=b1*x+b2*y+b3*z, C=c1*x+c2*y+c3*z(1)
以上所有这些量都是已知的,未知的是点X在A,B,C基矢坐标系的坐标。设为(x1,y1,z1),则列方程
x1*A+y1*B+z1*C=x0*x+y0*y+z0*z,
将(1)式代入上式,并解三元一次方程组,(x1,y1,z1)可解

另外5楼给出的是转化矩阵,作用应该是利用分数坐标(x1,y1,z1)乘矩阵,可以得到直角坐标。不过应该有适用条件吧,应该保证A,B,C的方向与直角坐标系方向有某种关系,只定义晶格常数是不够的。同样的一个晶格,我把它相对直角坐标进行胡乱旋转、平移,那么直角坐标系内点所对应的分数坐标肯定也是会跟着发生改变的。因此必须保证A,B,C与直角坐标系方向有某种关系。
8楼2010-05-27 23:56:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flashpoint

金虫 (小有名气)

我终于找到了原文献
你可以看看,

http://scripts.iucr.org/cgi-bin/paper?S0021889800006579
引用回帖:
Originally posted by g_xq96 at 2010-05-27 18:45:36:

这里哪里是分数呀,分数坐标得到的是分数。

9楼2010-05-28 03:39:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

引用回帖:
Originally posted by acridine at 2010-05-27 23:56:06:
勤劳一点,自己推导一下还是不难的。
哈哈,好就没做数学题了,貌似前两天有个贴子问过一个类似的情况,当时做过类似的推导:

设直角坐标系三个基矢为x,y和z,设点X在直角坐标系坐标为(x0,y0,z0)。
设晶胞的 ...

应该不用考虑特别的关系.
两套线性无关的三个三维矢量,是等价的.
10楼2010-05-28 08:28:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 g_xq96 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13 晴天加油 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见