Dear&nbsp;apsi:<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Thanks&nbsp;for&nbsp;your&nbsp;reply,&nbsp;but&nbsp;the&nbsp;problem&nbsp;is&nbsp;still&nbsp;exist!<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;I'm&nbsp;troubled&nbsp;in&nbsp;getting&nbsp;the&nbsp;cohesive&nbsp;energy&nbsp;of&nbsp;Al.&nbsp;The&nbsp;cohesive&nbsp;energy&nbsp;calcu<br>lation&nbsp;are&nbsp;performed&nbsp;by&nbsp;subtracting&nbsp;the&nbsp;total&nbsp;energy(E1)&nbsp;of&nbsp;the&nbsp;spin&nbsp;ploarized&nbsp;Al<br>&nbsp;atom&nbsp;from&nbsp;the&nbsp;energy(E2)&nbsp;per&nbsp;atom&nbsp;in&nbsp;the&nbsp;bulk.&nbsp;So&nbsp;I&nbsp;must&nbsp;get&nbsp;the&nbsp;different&nbsp;energ<br>ys&nbsp;E1&nbsp;and&nbsp;E2&nbsp;firstly.<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;In&nbsp;your&nbsp;last&nbsp;letter,&nbsp;you&nbsp;said&nbsp;&quot;The&nbsp;easiest&nbsp;way&nbsp;to&nbsp;calculate&nbsp;the&nbsp;energy&nbsp;of&nbsp;th<br>e&nbsp;single&nbsp;atom&nbsp;is&nbsp;to&nbsp;use&nbsp;MULTIPLICITY&nbsp;of&nbsp;4&nbsp;(three&nbsp;unpaired&nbsp;electrons)&nbsp;and&nbsp;no&nbsp;FREE&nbsp;<br>ENERGY&nbsp;FUNCTIONAL&nbsp;(and&nbsp;you&nbsp;can&nbsp;erase&nbsp;keyword&nbsp;'STATES').&nbsp;Please&nbsp;notice&nbsp;that&nbsp;you&nbsp;ha<br>ve&nbsp;a&nbsp;typo&nbsp;in&nbsp;'TEMPERATURE&nbsp;ELCECTRON',&nbsp;however&nbsp;in&nbsp;this&nbsp;case&nbsp;it&nbsp;wouldn't&nbsp;matter&nbsp;bec<br>ause&nbsp;the&nbsp;default&nbsp;value&nbsp;is&nbsp;the&nbsp;same,&nbsp;but&nbsp;in&nbsp;case&nbsp;of&nbsp;molecular&nbsp;dynamics&nbsp;this&nbsp;would&nbsp;<br>be&nbsp;interpreted&nbsp;as&nbsp;_ionic_&nbsp;temperature...&quot;,&nbsp;so&nbsp;need&nbsp;I&nbsp;use&nbsp;the&nbsp;FREE&nbsp;ENERGY&nbsp;FUNCTION<br>AL&nbsp;if&nbsp;I&nbsp;want&nbsp;to&nbsp;get&nbsp;the&nbsp;energy(E2)&nbsp;per&nbsp;atom&nbsp;in&nbsp;the&nbsp;bulk?&nbsp;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;The&nbsp;energy&nbsp;of&nbsp;spin&nbsp;ploarized&nbsp;Al&nbsp;atom&nbsp;I&nbsp;get&nbsp;is&nbsp;-1.95126380&nbsp;A.U.&nbsp;and&nbsp;-1.823879<br>24&nbsp;A.U.&nbsp;&nbsp;when&nbsp;I&nbsp;specified&nbsp;MULTIPLICITY&nbsp;2&nbsp;and&nbsp;4&nbsp;respectively.(In&nbsp;my&nbsp;mind,&nbsp;the&nbsp;elec<br>tron&nbsp;configuration&nbsp;of&nbsp;Al&nbsp;is&nbsp;1s(2)2s(2)2p(6)3s(2)3p(1)&nbsp;and&nbsp;one&nbsp;unpair&nbsp;electron,&nbsp;so<br>&nbsp;I&nbsp;think&nbsp;the&nbsp;multiplicity&nbsp;is&nbsp;2*(1/2)+1=2.But&nbsp;you&nbsp;specified&nbsp;MULTIPLICITY&nbsp;4&nbsp;in&nbsp;your<br>&nbsp;letter,&nbsp;does&nbsp;it&nbsp;means&nbsp;the&nbsp;electron&nbsp;configuration&nbsp;is&nbsp;1s(2)2s(2)2p(6)3s(1)3P(2)&nbsp;an<br>d&nbsp;three&nbsp;unpair&nbsp;electrons,&nbsp;so&nbsp;what&nbsp;are&nbsp;the&nbsp;difference&nbsp;and&nbsp;which&nbsp;multiplicity&nbsp;value<br>&nbsp;should&nbsp;I&nbsp;specified&nbsp;when&nbsp;I&nbsp;get&nbsp;the&nbsp;energy&nbsp;of&nbsp;spin&nbsp;ploarized&nbsp;Al&nbsp;atom?)&nbsp;The&nbsp;input&nbsp;f<br>ile&nbsp;for&nbsp;the&nbsp;spin&nbsp;Al&nbsp;atom&nbsp;calculation&nbsp;is&nbsp;following:<br>*******************************************************<br>&amp;CPMD<br>&nbsp;&nbsp;OPTIMIZE&nbsp;WAVEFUNCTION<br>&nbsp;&nbsp;CONVERGENCE&nbsp;ORBITALS<br>&nbsp;&nbsp;1.0D-6<br>&nbsp;&nbsp;PCG&nbsp;MINIMIZE<br>&nbsp;&nbsp;TIMESTEP<br>&nbsp;&nbsp;20<br>&nbsp;&nbsp;LSD<br>&nbsp;&nbsp;MIRROR<br>&amp;END<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br>&amp;SYSTEM<br>&nbsp;&nbsp;SYMMETRY<br>&nbsp;&nbsp;0<br>&nbsp;&nbsp;MULTIPLICITY<br>&nbsp;&nbsp;2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;(4&nbsp;or&nbsp;2?)<br>&nbsp;&nbsp;POISSON&nbsp;SOLVER&nbsp;TUCKERMAN<br>&nbsp;&nbsp;ANGSTROM<br>&nbsp;&nbsp;CELL<br>&nbsp;&nbsp;10.0&nbsp;&nbsp;1.0&nbsp;&nbsp;1.0&nbsp;&nbsp;0.0&nbsp;&nbsp;0.0&nbsp;&nbsp;0.0<br>&nbsp;&nbsp;CUTOFF<br>&nbsp;&nbsp;30<br>&amp;END<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br>&amp;ATOMS<br>*Al_SGS_LDA.psp&nbsp;KLEINMAN-BYLANDER<br>&nbsp;LMAX=D<br>&nbsp;1<br>&nbsp;5.0&nbsp;&nbsp;5.0&nbsp;&nbsp;5.0<br>&amp;END<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br>&amp;DFT<br>&nbsp;&nbsp;FUNCTIONAL&nbsp;LDA<br>&amp;END<br>*******************************************************&nbsp;<br><br>Then&nbsp;I&nbsp;calculate&nbsp;the&nbsp;energy&nbsp;per&nbsp;atom&nbsp;in&nbsp;bulk&nbsp;with&nbsp;FREE&nbsp;ENERGY&nbsp;FUNCTIONAL&nbsp;and&nbsp;with<br>out&nbsp;it,&nbsp;the&nbsp;E2&nbsp;are&nbsp;-2.10337731&nbsp;A.U.&nbsp;and&nbsp;-2.01913733&nbsp;A.U.&nbsp;respectively,&nbsp;the&nbsp;follow<br>ing&nbsp;are&nbsp;the&nbsp;input&nbsp;files&nbsp;with&nbsp;and&nbsp;without&nbsp;free&nbsp;energy&nbsp;funcitonal:<br><br>******************************************************<br>&amp;CPMD<br>&nbsp;&nbsp;OPTIMIZE&nbsp;WAVEFUNCTION<br>&nbsp;&nbsp;FREE&nbsp;ENERGY&nbsp;FUNCTIONAL<br>&nbsp;&nbsp;CONVERGENCE&nbsp;ORBITALS<br>&nbsp;&nbsp;1.0D-6<br>&nbsp;&nbsp;MIRROR<br>&amp;END<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br>&amp;SYSTEM<br>&nbsp;&nbsp;SYMMETRY<br>&nbsp;&nbsp;FCC<br>&nbsp;&nbsp;ANGSTROM<br>&nbsp;&nbsp;CELL<br>&nbsp;&nbsp;3.95&nbsp;&nbsp;1.0&nbsp;&nbsp;1.0&nbsp;&nbsp;0&nbsp;&nbsp;0&nbsp;&nbsp;0&nbsp;&nbsp;&nbsp;(#I&nbsp;get&nbsp;the&nbsp;equilibration&nbsp;lattice&nbsp;with&nbsp;the&nbsp;energy&nbsp;vs&nbsp;<br>variety&nbsp;lattice)<br>&nbsp;&nbsp;CUTOFF<br>&nbsp;&nbsp;60<br>&nbsp;&nbsp;TESR<br>&nbsp;&nbsp;3<br>&nbsp;&nbsp;STATES<br>&nbsp;&nbsp;4<br>&nbsp;&nbsp;KPOINTS&nbsp;MONKHORST-PACK<br>&nbsp;&nbsp;8&nbsp;8&nbsp;8<br>&amp;END<br><br>&amp;ATOMS<br>*Al_SGS_LDA.psp&nbsp;KLEINMAN-BYLANDER<br>&nbsp;&nbsp;LMAX=D<br>&nbsp;&nbsp;1<br>&nbsp;&nbsp;0&nbsp;&nbsp;0&nbsp;&nbsp;0<br>&amp;END<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br>&amp;DFT<br>&nbsp;&nbsp;FUNCTIONAL&nbsp;LDA<br>&amp;END<br>*****************************************************<br><br>******************************************************<br>&amp;CPMD<br>&nbsp;&nbsp;OPTIMIZE&nbsp;WAVEFUNCTION<br>&nbsp;&nbsp;CONVERGENCE&nbsp;ORBITALS<br>&nbsp;&nbsp;1.0D-6<br>&nbsp;&nbsp;MIRROR<br>&amp;END<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br>&amp;SYSTEM<br>&nbsp;&nbsp;SYMMETRY<br>&nbsp;&nbsp;FCC<br>&nbsp;&nbsp;ANGSTROM<br>&nbsp;&nbsp;CELL<br>&nbsp;&nbsp;3.95&nbsp;&nbsp;1.0&nbsp;&nbsp;1.0&nbsp;&nbsp;0&nbsp;&nbsp;0&nbsp;&nbsp;0<br>&nbsp;&nbsp;CUTOFF<br>&nbsp;&nbsp;60<br>&nbsp;&nbsp;TESR<br>&nbsp;&nbsp;3<br>&nbsp;&nbsp;KPOINT&nbsp;MONKHORST-PACK<br>&nbsp;&nbsp;8&nbsp;8&nbsp;8<br>&amp;END<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br>&amp;ATOMS<br>*Al_SGS_LDA.psp&nbsp;KLEINMAN-BYLANDER<br>&nbsp;&nbsp;LMAX=D<br>&nbsp;&nbsp;1<br>&nbsp;&nbsp;0&nbsp;&nbsp;0&nbsp;&nbsp;0<br>&amp;END<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br>&amp;DFT<br>&nbsp;&nbsp;FUNCTIONAL&nbsp;LDA<br>&amp;END<br>**********************************************************<br><br>the&nbsp;cohesive&nbsp;energy&nbsp;I&nbsp;get&nbsp;with&nbsp;the&nbsp;above&nbsp;the&nbsp;values&nbsp;are&nbsp;far&nbsp;away&nbsp;the&nbsp;experiment&nbsp;v<br>alues&nbsp;3.45&nbsp;ev.&nbsp;Could&nbsp;you&nbsp;help&nbsp;me&nbsp;solve&nbsp;the&nbsp;problem&nbsp;and&nbsp;get&nbsp;the&nbsp;calculation&nbsp;cohesi<br>ve&nbsp;value&nbsp;approach&nbsp;the&nbsp;experiment&nbsp;values?<br>&nbsp;&nbsp;&nbsp;Another&nbsp;problem&nbsp;about&nbsp;project-dos&nbsp;is&nbsp;that&nbsp;I&nbsp;want&nbsp;to&nbsp;get&nbsp;the&nbsp;information&nbsp;that&nbsp;t<br>he&nbsp;dos&nbsp;vs&nbsp;energy&nbsp;on&nbsp;s,p,d&nbsp;...orbitals,&nbsp;I&nbsp;check&nbsp;the&nbsp;LDOS&nbsp;and&nbsp;PROJECT&nbsp;ANALYSIS&nbsp;and&nbsp;<br>don't&nbsp;know&nbsp;how&nbsp;to&nbsp;get&nbsp;the&nbsp;information&nbsp;I&nbsp;need,&nbsp;could&nbsp;you&nbsp;explain&nbsp;it&nbsp;to&nbsp;me&nbsp;in&nbsp;detai<br>ls?<br>&nbsp;&nbsp;&nbsp;I'm&nbsp;poor&nbsp;on&nbsp;the&nbsp;sphere&nbsp;of&nbsp;ab&nbsp;initio&nbsp;calculation&nbsp;and&nbsp;maybe&nbsp;some&nbsp;questions&nbsp;is&nbsp;st<br>upid&nbsp;and&nbsp;not&nbsp;professional,&nbsp;but&nbsp;I&nbsp;really&nbsp;need&nbsp;your&nbsp;help&nbsp;thirstly!<br>&nbsp;&nbsp;&nbsp;Thanks&nbsp;in&nbsp;advace!<br>&nbsp;&nbsp;&nbsp;Best&nbsp;regards!<br>&nbsp;&nbsp;&nbsp;shyma&nbsp;&nbsp;<br><br>