| 9 | 1/1 | 返回列表 |
| 查看: 2111 | 回復: 8 | ||
yangxuezhang金蟲 (正式寫手)
|
[求助]
請大家?guī)兔纯磗iesta輸入文件有哪些錯誤
|
|
我這是第一次寫輸入文件,以前也沒師兄做這個,也每個人問,希望高手不要嫌低級,幫我看一下,哪里有錯誤,不甚感激 SystemName Test for Ti3SiC2 SystemLabel Ti3SiC2 NumberOfAtoms 11 # Number of atoms NumberOfSpecies 3 # Number of species %block Chemical_Species_Label 1 22 Ti 2 6 C 3 14 Si %endblock Chemical_Species_Label PAO.BasisSize SZ #??? # Lattice, coordinates, k-sampling AtomicCoordinatesFormat ScaledCartesian # Format for coordinates AtomicCoorFormatOut Ang %block AtomicCoordinatesAndAtomicSpecies 0.000000000 0.000000000 0.000000000 1 0.000000000 0.000000000 4.401250000 3 0.000000000 0.000000000 13.20375 3 1.541000000 0.889696765 1.271081005 2 1.541000000 0.889696765 7.531419 2 0.000000000 1.779393530 10.073581 2 0.000000000 1.779393530 16.333919 2 0.000000000 1.779393530 2.385477481 1 0.000000000 1.779393530 6.4170225 1 1.541000000 0.889696765 11.1879775 1 1.541000000 0.889696765 15.2195225 1 %endblock AtomicCoordinatesAndAtomicSpecies LatticeConstant 3.082 Ang %block LatticeVectors 1.000 0.500 0.000 0.000 0.866 0.000 0.000 0.000 1.000 %endblock LatticeVectors kgrid_cutoff 7. Ang #????? # DFT, Grid, SCF XC.functional LDA # Exchange-correlation functional type XC.authors CA # Particular parametrization of xc func SpinPolarized .false. # Spin unpolarized calculation MeshCutoff 200. Ry # Equivalent planewave cutoff for the grid MaxSCFIterations 100 # Maximum number of SCF iterations per step DM.MixingWeight 0.3 # New DM amount for next SCF cycle DM.Tolerance 1.d-4 # Tolerance in maximum difference # between input and output DM DM.NumberPulay 3 # Number of SCF steps between pulay mixing # Eigenvalue problem: order-N or diagonalization SolutionMethod diagon # OrderN or Diagon ElectronicTemperature 5 K # Temp. for Fermi smearing # Molecular dynamics and relaxations MD.TypeOfRun cg # Type of dynamics: # Output options WriteCoorInitial WriteCoorStep WriteForces WriteKpoints .false. WriteEigenvalues .false. WriteKbands .false. WriteBands .false. WriteMullikenPop 1 # Write Mulliken Population Analysis WriteCoorXmol .false. WriteMDCoorXmol .false. WriteMDhistory .false. WriteCoorXmol .false. # Options for saving/reading information DM.UseSaveDM # Use DM Continuation files MD.UseSaveXV .false. # Use stored positions and velocities MD.UseSaveCG .false. # Use stored positions and velocities SaveRho # Write valence pseudocharge at the mesh SaveDeltaRho # Write RHOscf-RHOatm at the mesh SaveElectrostaticPotential .false. # Write the total elect. pot. at the mesh SaveTotalPotential .false. # Write the total pot. at the mesh WriteSiestaDim .false. # Write minimum dim to siesta.h and stop WriteDenchar # Write information for DENCHAR |
金蟲 (正式寫手)
金蟲 (正式寫手)
|
這是提交后的scrip.err文件里面的 ERROR STOP from Node: 0 -------------------------------------------------------------------------- MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD with errorcode 1. NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them. -------------------------------------------------------------------------- -------------------------------------------------------------------------- mpiexec has exited due to process rank 0 with PID 14863 on node dell745-2-1.local exiting without calling "finalize". This may have caused other processes in the application to be terminated by signals sent by mpiexec (as reported here). -------------------------------------------------------------------------- real 0m0.866s user 0m0.041s sys 0m0.054s |
木蟲 (著名寫手)

金蟲 (正式寫手)
|
下面是*.out文件里面的,為了找錯誤,和前面的輸入文件比,我注釋掉了很多項,不過還是不行,不知道是什么問題 Siesta Version: siesta-3.1 Architecture : x86_64-unknown-linux-gnu--Intel Compiler flags: /opt/openmpi-intel11/bin/mpif90 -g -O2 PARALLEL version * Running on 2 nodes in parallel >> Start of run: 24-SEP-2012 15:08:55 *********************** * WELCOME TO SIESTA * *********************** reinit: Reading from standard input ************************** Dump of input data file **************************** SystemName Test for Ti3SiC2 SystemLabel Ti3SiC2 NumberOfAtoms 11 # Number of atoms NumberOfSpecies 3 # Number of species %block Chemical_Species_Label 1 22 Ti 2 6 C 3 14 Si %endblock Chemical_Species_Label #PAO.BasisSize SZ #??? # Lattice, coordinates, k-sampling AtomicCoordinatesFormat Ang # Format for coordinates AtomicCoorFormatOut Ang %block AtomicCoordinatesAndAtomicSpecies 0.000000000 0.000000000 0.000000000 1 0.000000000 0.000000000 4.401250000 3 0.000000000 0.000000000 13.20375 3 1.541000000 0.889696765 1.271081005 2 1.541000000 0.889696765 7.531419 2 0.000000000 1.779393530 10.073581 2 0.000000000 1.779393530 16.333919 2 0.000000000 1.779393530 2.385477481 1 0.000000000 1.779393530 6.4170225 1 1.541000000 0.889696765 11.1879775 1 1.541000000 0.889696765 15.2195225 1 %endblock AtomicCoordinatesAndAtomicSpecies LatticeConstant 3.082 Ang %block LatticeVectors 1.000 0.500 0.000 0.000 0.866 0.000 0.000 0.000 1.000 %endblock LatticeVectors #kgrid_cutoff 7. Ang #????? # DFT, Grid, SCF #XC.functional LDA # Exchange-correlation functional type #XC.authors CA # Particular parametrization of xc func #SpinPolarized .false. # Spin unpolarized calculation #MeshCutoff 200. Ry # Equivalent planewave cutoff for the grid #MaxSCFIterations 100 # Maximum number of SCF iterations per step #DM.MixingWeight 0.3 # New DM amount for next SCF cycle #DM.Tolerance 1.d-4 # Tolerance in maximum difference # between input and output DM #DM.NumberPulay 3 # Number of SCF steps between pulay mixing # Eigenvalue problem: order-N or diagonalization SolutionMethod diagon # OrderN or Diagon #ElectronicTemperature 5 K # Temp. for Fermi smearing # Molecular dynamics and relaxations #MD.TypeOfRun cg # Type of dynamics: # Output options #WriteCoorInitial #WriteCoorStep #WriteForces #WriteKpoints .false. #WriteEigenvalues .false. #WriteKbands .false. #WriteBands .false. #WriteMullikenPop 1 # Write Mulliken Population Analysis #WriteCoorXmol .false. #WriteMDCoorXmol .false. #WriteMDhistory .false. #WriteCoorXmol .false. # Options for saving/reading information #DM.UseSaveDM # Use DM Continuation files #MD.UseSaveXV .false. # Use stored positions and velocities #MD.UseSaveCG .false. # Use stored positions and velocities #SaveRho # Write valence pseudocharge at the mesh #SaveDeltaRho # Write RHOscf-RHOatm at the mesh #SaveElectrostaticPotential .false. # Write the total elect. pot. at the mesh #SaveTotalPotential .false. # Write the total pot. at the mesh #WriteSiestaDim .false. # Write minimum dim to siesta.h and stop #WriteDenchar # Write information for DENCHAR ************************** End of input data file ***************************** reinit: ----------------------------------------------------------------------- reinit: System Name: Test for Ti3SiC2 reinit: ----------------------------------------------------------------------- reinit: System Label: Ti3SiC2 reinit: ----------------------------------------------------------------------- initatom: Reading input for the pseudopotentials and atomic orbitals ---------- Species number: 1 Label: Ti Atomic number: 22 Species number: 2 Label: C Atomic number: 6 Species number: 3 Label: Si Atomic number: 14 Ground state valence configuration: 4s02 3d02 Reading pseudopotential information in formatted form from Ti.psf Pseudopotential generated from a relativistic atomic calculation There are spin-orbit pseudopotentials available Spin-orbit interaction is not included in this calculation Valence configuration for pseudopotential generation: 3s( 2.00) rc: 1.30 3p( 6.00) rc: 1.30 3d( 2.00) rc: 1.30 4f( 0.00) rc: 1.98 Ground state valence configuration: 2s02 2p02 Reading pseudopotential information in formatted form from C.psf Valence configuration for pseudopotential generation: 2s( 2.00) rc: 1.49 2p( 2.00) rc: 1.52 3d( 0.00) rc: 1.58 Ground state valence configuration: 3s02 3p02 Reading pseudopotential information in formatted form from Si.psf Valence configuration for pseudopotential generation: 3s( 2.00) rc: 1.89 3p( 2.00) rc: 1.89 3d( 0.00) rc: 1.89 4f( 0.00) rc: 1.89 Semicore shell(s) with 8 electrons included in the valence for Ti Ti must be in PAO.Basis (it has semicore states) ERROR STOP from Node: 0 希望大俠多多幫忙 |
木蟲 (著名寫手)
|
從輸出文件中看出來,你計算中的Ti元素的贗勢將8個本應該屬于core的電子放到了價電子的考慮范圍,這樣可以更好的描述不同環(huán)境中的Ti的電子結構。但是你需要對Ti的基函數(shù)進行手動設置。而不能用簡單的PAO.BasisSize SZ關鍵字,使用這個關鍵字默認只給價電子產(chǎn)生基函數(shù)。所以你需要使用 %block PAO.Basis %endblock PAO.Basis 來設置,比如 PAO.BasisSize DZP PAO.EnergyShift 100 meV 對應的C元素的基函數(shù)為: %block PAO.Basis C 2 # Species label, number of l-shells n=2 0 2 # n, l, Nzeta 4.633 3.432 1.000 1.000 n=2 1 2 P 1 # n, l, Nzeta, Polarization, NzetaPol 5.659 3.653 1.000 1.000 %endblock PAO.Basis 當然,這個與具體的贗勢還有點關系。 你可以比照著這個對Ti進行一下設置,將那8個電子的基函數(shù)明確的寫出來。我沒有算過Ti,所以不會設置。 還有一個辦法就是你可以使用一個沒有semicore態(tài)電子的Ti贗勢。你從官網(wǎng)上下載一個就是了。我看了一下官網(wǎng)上的LDA和GGA的Ti贗勢就都沒有semicore態(tài)。 http://icmab.cat/leem/siesta/Dat ... init/Tl_html/Tl.psf http://icmab.cat/leem/siesta/Dat ... init/Tl_html/Tl.psf |

金蟲 (正式寫手)

金蟲 (正式寫手)
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 能源材料化學課題組招收碩士研究生8-10名 +3 | 脫穎而出 2026-03-16 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +5 | yy要上岸呀 2026-03-17 | 5/250 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
[考研] 一志愿南京大學,080500材料科學與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 274求調(diào)劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] [導師推薦]西南科技大學國防/材料導師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 286求調(diào)劑 +3 | lemonzzn 2026-03-16 | 5/250 |
|
|
[考研] 0703化學調(diào)劑 +6 | 妮妮ninicgb 2026-03-15 | 9/450 |
|
|
[考研] 283求調(diào)劑 +10 | 小樓。 2026-03-12 | 14/700 |
|
|
[考研] 331求調(diào)劑(0703有機化學 +5 | ZY-05 2026-03-13 | 6/300 |
|
|
[考研] 279求調(diào)劑 +3 | 抓著星星的女孩 2026-03-10 | 3/150 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] 290求調(diào)劑 +9 | ADT 2026-03-11 | 9/450 |
|
|
[考研] 工科,求調(diào)劑 +3 | 我887 2026-03-11 | 3/150 |
|
|
[考研] 270求調(diào)劑 085600材料與化工專碩 +3 | YXCT 2026-03-11 | 3/150 |
|
|
[考研] 321求調(diào)劑(食品/專碩) +3 | xc321 2026-03-12 | 6/300 |
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|