| 查看: 585 | 回复: 3 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】怎么将transmission在origin里作图已有2人参与
|
|||
| 我想问一下各位前辈，怎么将ATK 中的tarnsmission在origion 里面作图 |
回复此楼» 猜你喜欢
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有9人回复
-
心脉受损
已经有6人回复
-
博士申请都是内定的吗？
已经有8人回复
-
博士读完未来一定会好吗
已经有32人回复
-
读博
已经有5人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
Springer期刊投稿求助
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
求助!急。。。急。。origin作图后上下各出现一条线
已经有14人回复- origin 作图坐标轴问题请教
已经有6人回复
- 请问TXT交流阻抗数据导出之后怎么在origin作图？
已经有6人回复
- 只有一列的数据，怎么用origin作图？
已经有4人回复
- 请教origin作图问题
已经有11人回复
- 【求助】XRD数据如何在origin里作图？
已经有17人回复
- origin 作图请教
已经有7人回复
- origin作图问题
已经有9人回复
- 【求助】怎么在origin里面作图知道曲线那点的斜率最大？
已经有5人回复
163.com
金虫 (著名写手)
实习版主
- 应助: 2 (幼儿园)
- 贵宾: 0.016
- 金币: 1489.8
- 散金: 1347
- 红花: 32
- 帖子: 1215
- 在线: 87小时
- 虫号: 923817
- 注册: 2009-12-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
gavinliu7390(金币+2):谢谢交流！ 2010-06-04 22:19:00
小木虫(金币+0.5):给个红包，谢谢回帖交流
gavinliu7390(金币+2):谢谢交流！ 2010-06-04 22:19:00
|
你在程序最后加上这段脚本，T就会写入到 “trans.dat” 中： f = open('trans.dat','w') import numpy transmission_spectrum = calculateTransmissionSpectrum( self_consistent_calculation = scf, energies = numpy.arange(-5.0, 5.0+0.002, 0.02)*electronVolt, brillouin_zone_integration_parameters = brillouinZoneIntegrationParameters((1, 1)), green_function_infinitesimal = 1.0e-5*electronVolt ) energies = transmission_spectrum.energies() coefficients = transmission_spectrum.coefficients() if len(coefficients.shape)==2: # spin-polarized print ' Transmission' print 'Energy (eV) Spin-up Spin-down' print '------------------------------------------' for i in range(len(energies)): print "%g\t%g\t%g\n" % ( energies.inUnitsOf(eV), coefficients[0,i],coefficients[1,i] ) s = '%g\t%g\t%g\n' % ( energies.inUnitsOf(eV), coefficients[0,i],coefficients[1,i] ) f.write(s) else: print 'Energy (eV) Transmission' print '-----------------------------------' for i in range(len(energies)): print "%g\t%g\n" % ( energies.inUnitsOf(eV),coefficients ) s = '%g\t%g\n' % ( energies.inUnitsOf(eV),coefficients ) f.write(s) f.close() |
4楼2010-06-04 16:37:10
163.com
金虫 (著名写手)
实习版主
- 应助: 2 (幼儿园)
- 贵宾: 0.016
- 金币: 1489.8
- 散金: 1347
- 红花: 32
- 帖子: 1215
- 在线: 87小时
- 虫号: 923817
- 注册: 2009-12-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2010-06-03 23:49:54
|
更高层次的是修改脚本，程序自动生成Transmission,我想问一下什么命令可以呢？ from ATK.TwoProbe import * # Restoring the SCF from the previous calculation scf = restoreSelfConsistentCalculation('si-bias-1.0.nc') # Create a list of energies from -5 to 5 eV, with 0.1 eV spacing import numpy energy_list = numpy.arange(-5.0, 5.0, 0.1)*electronVolt # Set k-points for transmission bz_int_param = brillouinZoneIntegrationParameters( (5,5,100) ) spectrum = calculateTransmissionSpectrum( self_consistent_calculation=scf, energies = energy_list, brillouin_zone_integration_parameters=bz_int_param) # Save the transmission spectrum to VNL file vnlfile = VNLFile("transmission.vnl"  vnlfile.addToSample(spectrum,'transmission') # Print at screen the Transmission Spectrum print 'Transmission Spectrum' print '-----------------------------------------------------------------------' print 'Energies(eV) Coefficients' for i in range(len(spectrum.energies())): print '%17.3f %17.6f' %( spectrum.energies().inUnitsOf(Units.eV), spectrum.coefficients()) 这是我的输入文件，但是我想导出一个数据文件，你能告诉我怎么修改吗 |
3楼2010-06-04 10:14:41