VASP Calculation Tutorial: Transition State Calculations of Adsorbed Atom Diffusion on Platinum Surfaces
1. Research Background and Theoretical Foundation
Transition state theory holds significant importance in surface science and catalysis research, helping us understand the diffusion behavior of atoms and molecules on material surfaces. As an important noble metal catalyst, studying the surface atom diffusion behavior of platinum (Pt) is crucial for understanding catalytic reaction mechanisms. This article will detail how to use the Nudged Elastic Band (NEB) method in VASP software to calculate diffusion energy barriers using Pt adsorbed atoms' self-diffusion process on a Pt(001) surface as an example.
The NEB method is an effective algorithm for calculating minimum energy paths (MEP), proposed by Jónsson et al. in 1998 and continuously refined since then. This method inserts several intermediate images between reactant and product configurations, applying elastic forces between adjacent images to converge the entire "elastic band" onto the minimum energy path. In a converged state, the image with the highest energy corresponds to the transition state configuration, with its energy difference from that of the initial state representing the diffusion barrier.
2. Computational System and Model Construction
2.1 Crystal Structure of Platinum and Surface Model
Platinum exhibits a face-centered cubic (fcc) structure at room temperature and pressure, with its (001) surface showing significant reconstruction phenomena. To accurately simulate surface properties, we constructed a supercell model with at least 2×2 reconstructions. Calculations indicate that Pt atoms have their most stable adsorption site at hollow positions marked as h sites on Pt(001). Additionally, there are two metastable adsorption sites: bridge position (b site) and top position (t site).
2.2 Diffusion Path Analysis
Adsorbed atoms diffuse from initial h sites to neighboring h sites via two main pathways: first is hollow-top-hollow (hth), along [1-10] crystal direction; second is hollow-bridge-hollow (hbh), along [100] crystal direction. Studies show that collective hopping mechanism involving two Pt atoms moving cooperatively along [1-10] has the lowest barrier height, serving as this system's primary diffusion mechanism.
3. Computational Methods and Parameter Settings
3.1 Density Functional Theory Computation We used VASP software based on plane-wave basis sets for first-principles calculations; PBE form was chosen for exchange-correlation functional settings: p - Plane wave cutoff energy set at 400 eV, k-point mesh utilized was a Monkhorst-Pack grid of 3×3×1, electronic step convergence criterion set at EDIFF = 1×10⁻⁶ eV, and ionic relaxation convergence criterion set at EDIFFG = -0 .01 eV/Å. n **3 .2 Details about Surface Model Construction ** n Modeling clean Pt(001 ) surfaces requires considering several key factors : First , sufficient layers must be present to avoid quantum size effects ; it’s recommended to use no less than three layers of Pt atoms ; Second , vacuum layer thickness should be large enough (>15 Å )to prevent interactions between upper/lower surfaces ; Finally , reconstruction phenomena must be adequately accounted for—at least basic requirement being a (2 imes{ }{ }{ }{ } { }) o c o n s t r u c t i o n . n ###4.NEBCalculationProcessDetailedExplanation **4 . InitialPreparationStage **Before conducting NEB calculations,two endpoint structures need preparation :initialstate(POSCAR_i )andfinalstate(POSCAR_f ).These files must meet following requirements :no atomic name line(inclusive only within VASP version5 .0 );no“selectivedynamics”keyword;no empty lines;atomic velocity informationmustbe removed.Speciallynotedisthatallatomsinbothinitial/finalstructuresmustremainonthe same sideofsupercelltoavoidinterpolationissueswithatomcrossingunit cell boundaries. * 4 . GeneratingIntermediateImagesUsing interpolatePOS script generates intermediateimages through linear interpolation methodsbetweeninitial/final states.The script specifies total steps number(rep=images+two);usersneedmodifythisparameterbasedonactualnumberdesired.Forinstance,twointermediateimagesrequire settingrep=four.Generatedimageswillbestoredinsubdirectoriesfrom00to0(n+one). * 4 . OptimizingComputationalParametersKeyparametersforNEBcalculationsinclude:imagecount(IMAGES ),springconstant(SPRING ),maximumionicsteps(NSW ),etc.Practiceindicatesusingfewerimagecounts(not exceeding fourusually)yieldsbetterconvergence.Ifmoreprecisedescriptionaroundtransitionstatesneeded,a stagedstrategycanbeemployed:firstconductlowprecisioncalculatingapproximatepositionofthetransitionstatebeforeutilizingadjacentimagestothetransitionstatetoconductsecondroundcalculations ###5.ResultAnalysisandDiscussion *5 . BarrierCalculations&PathOptimization Uponcompletionofcalculation,splineinterpolation(cubictypicallyused)makespossibleobtainingpotentialenergylandscapealongdiffusionpath.Highestpointcorrespondswiththetransitionalstatewhoseenergydifferencefromtheinitialonegivesitsbarrierheight.InPt@Pt(001 )system,hthpathhasabarrierapproximatelyequaltozero pointseventy-fiveeVwhilehbhpathaboutzero pointeighty-fiveeVisconsistentwithexperimentalobservedresults.NoteForfastcomputingsettings(pt_neb_fast.tgz),theseareonlysuitableforlearningbasicprinciplesbehindNEBcomputing,andresultscannotreproduceexperimentalvalues.Toobtainreliableoutcomes,mustusehigheraccuracysettings(pt_neb.tgz):includinggreaterlayernumbers(recommendedfiveormore),denserk-pointmeshes(atleastfivebyfivebyone),highercutoffenergies(over500eVorstricterconvergencecriteria.) ###6.References&FurtherReading Kellogg,G.L.&Feibelman,P.J.SurfaceSelf-DiffusiononPt(001 )byanAtomicExchangeMechanism.PhysicalReviewLetters64,(1990). Mills,G.,Jónsson,H.&Schenter,G.K.ReversibleWorkTransitionStateTheory:ApplicationToDissociativeAdsorptionOfHydrogen.SurfaceScience324,(1995). Jónsson,H.,Mills,G.&Jacobsen,K.W.NudgedElasticBandMethodForFindingMinimumEnergyPathsOfTransitions.InClassicalAndQuantumDynamicsInCondensedPhaseSimulations(ed.Berne,B.J.,Ciccotti,G.&Coker,D.F.)385–404(WorldScientific)(1998). This tutorial provides detailed guidance throughout all processes involvedinperformingNEBcalculationsusingVASPincludingtheoreticalfoundationalaspectsandpracticaloperationalsteps.Readersmayapplymethodstoother surfaceddiffusionsystemsmerelyadjustingmodelparametersaccordingly.Forcomplexmulti-atomcooperative diffusiveprocesses,increasingimagerangesandswitchingtomoreaccurate computationalsettingsmightbenecessary.