| 4 | 1/1 | 返回列表 |
| 查看: 643 | 回復: 3 | |||
[交流]
【求助】怎么將transmission在origin里作圖 已有2人參與
|
| 我想問一下各位前輩,怎么將ATK 中的tarnsmission在origion 里面作圖 |
金蟲 (著名寫手)
實習版主

|
更高層次的是修改腳本,程序自動生成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()) 這是我的輸入文件,但是我想導出一個數(shù)據(jù)文件,你能告訴我怎么修改嗎 |
金蟲 (著名寫手)
實習版主
|
你在程序最后加上這段腳本,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ā)表 | |
|---|---|---|---|---|
|
[考研] 343求調(diào)劑 +4 | 愛羈絆 2026-03-28 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 本科新能源科學與工程,一志愿華理能動285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +6 | 芋泥小鈴鐺 2026-03-28 | 6/300 |
|
|
[考研] 一志愿太原理工安全工程300分,求調(diào)劑 +4 | 0857求調(diào)劑. 2026-03-24 | 5/250 |
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-24 | 4/200 |
|
|
[考研] 266分,求材料冶金能源化工等調(diào)劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
[考研] 339求調(diào)劑,想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
[考研] 0703化學/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +9 | 丹青奶蓋 2026-03-26 | 10/500 |
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學-環(huán)境工程專碩-298 +4 | 手機用戶 2026-03-26 | 4/200 |
|
|
[考研] 328求調(diào)劑 +7 | 嗯滴的基本都 2026-03-27 | 7/350 |
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
[考研] 070300化學求調(diào)劑 +4 | 起個名咋這么難 2026-03-27 | 4/200 |
|
|
[考研] 一志愿華東理工大學081700,初試分數(shù)271 +6 | kotoko_ik 2026-03-23 | 7/350 |
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
[考研] 考研調(diào)劑 +10 | 呼呼?~+123456 2026-03-24 | 10/500 |
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
[考研] 求調(diào)劑 +8 | Auroracx 2026-03-22 | 8/400 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
[考研] 材料調(diào)劑 +3 | iwinso 2026-03-23 | 3/150 |
|