| 24小時熱門版塊排行榜 |
| 4 | 1/1 | 返回列表 |
| 查看: 644 | 回復(fù): 3 | |||
zhangying09銅蟲 (小有名氣)
|
[交流]
【求助】怎么將transmission在origin里作圖 已有2人參與
|
| 我想問一下各位前輩,怎么將ATK 中的tarnsmission在origion 里面作圖 |
金蟲 (著名寫手)
實習(xí)版主

銅蟲 (小有名氣)
|
更高層次的是修改腳本,程序自動生成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()) 這是我的輸入文件,但是我想導(dǎo)出一個數(shù)據(jù)文件,你能告訴我怎么修改嗎 |
金蟲 (著名寫手)
實習(xí)版主
|
你在程序最后加上這段腳本,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 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +7 | 芋泥小鈴鐺 2026-03-28 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 332求調(diào)劑 +5 | 蕉蕉123 2026-03-28 | 5/250 |
|
|
[考研] 一志愿北化085600材料專碩275|有文章專利|求調(diào)劑 +7 | Micky11223 2026-03-25 | 7/350 |
|
|
[考研] 321求調(diào)劑 +6 | 璞玉~~ 2026-03-25 | 7/350 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
[考研] 0703一志愿9,初試成績:338,四六級已過,有科研經(jīng)歷,求調(diào)劑! +4 | Zuhui0306 2026-03-25 | 4/200 |
|
|
[考研] 277跪求調(diào)劑 +5 | 1915668 2026-03-27 | 9/450 |
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
[考研] 291求調(diào)劑 +7 | 孅華 2026-03-22 | 7/350 |
|
|
[考研] 考研調(diào)劑 +4 | Sanmu-124 2026-03-26 | 4/200 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
[考研] 一志愿河工大 081700 276求調(diào)劑 +4 | 地球繞著太陽轉(zhuǎn) 2026-03-23 | 4/200 |
|
|
[考研] 296求調(diào)劑 +4 | 汪?! 2026-03-25 | 7/350 |
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|