<div style="line-height:1.7;color:#000000;font-size:14px;font-family:Arial"><p class="MsoNormal" style="font-family: Helvetica, "Microsoft Yahei", verdana; line-height: 28px;"><span style="font-family: "Times New Roman", serif; font-size: 12pt;">Dear experts,</span></p><div style="font-family: Helvetica, "Microsoft Yahei", verdana; line-height: 28px;"><span lang="EN-US" class="" style="font-size: 12pt; line-height: 32px; font-family: "Times New Roman", serif;">This is Yuan Fang. My research focuses on first principles simulation of infrared and Raman spectra of adsorbates at liquid/solid electrochemical interfaces. </span><span class="" style="font-family: "Times New Roman", serif; font-size: 12pt;">Recently, I plan to accurately calculate the adsorption energies of the Pt(111)(2¡Á2)-CO</span><span style="font-family: "Times New Roman", serif; font-size: 10px;">L</span><span class="" style="font-family: "Times New Roman", serif; font-size: 12pt;"> structure (shown in Fig.1) by LDA+U calculation with reading maximally localized Wannier function of  specified bands (pmw.x tool and "U_projection_type='file'" tag). Since the method can¡¯t automatically obtain the atomic forces in Quantum Espresso, I need manually calculate the atomic forces by the finite difference method to decide the next atomic movement and to optimize the structure finally.</span></div><div style="font-family: Helvetica, "Microsoft Yahei", verdana; line-height: 28px;"><span class="" style="font-family: "Times New Roman", serif; font-size: 12pt;"><br></span></div><div style="font-family: Helvetica, "Microsoft Yahei", verdana; line-height: 28px;"><span class="" style="font-family: "Times New Roman", serif; font-size: 12pt;"><img src="http://mail.163.com/js6/s?func=mbox:getMessageData&mid=50:1tbiMg7hBlWBu2pQrgAAsg&part=3" orgwidth="1338" orgheight="203" data-image="1" style="border: none; width: 1203px; height: 182px;"></span></div><div style="font-family: Helvetica, "Microsoft Yahei", verdana; line-height: 28px;"><span class="" style="font-family: "Times New Roman", serif; font-size: 12pt;">Fig.1 the strcture of Pt(111)(2¡Á2)-CO</span><span style="font-family: "Times New Roman", serif; font-size: 10px;">L</span><span class="" style="font-family: "Times New Roman", serif; font-size: 12pt;"> adsorption.</span></div><div style="font-family: Helvetica, "Microsoft Yahei", verdana; line-height: 28px;"><span class="" style="font-family: "Times New Roman", serif; font-size: 12pt;"><br></span></div><div style="font-family: Helvetica, "Microsoft Yahei", verdana; line-height: 28px;"><span class="" style="font-family: "Times New Roman", serif; font-size: 12pt;">However, I find that <u>the energy accuracy of LDA+U calculations with reading MLWF of specified bands is so low</u> that the differential value (atomic force in red line of Fig.2) deviates the scf calculation without LDA+U (the black line) and fails to calculate atomic force. Although I tried increasing the kpoint density (from 4*4*1 to 6*6*1), the scf convergence threshold (from 1E-7 Ry to 1E-8 even 1E-10 Ry), cutoff energy (from 40 Ry to 120 Ry), the spreading of Wannier function (from 0.02 to 0.002) or replacing Methfessel-Paxton first-order spreading with Marzari-Vanderbilt cold smearing, all of them can't improve the energy accuracy obviously in Fig.2. <u>Could you please teach me how to improve the </u></span><u style="font-family: "Times New Roman", serif; font-size: 16px;">the energy accuracy of LDA+U calculations with reading MLWF of specified bands? </u><u> </u><span style="font-family: "Times New Roman"; font-size: 16px;">Many thanks!</span></div><div style="font-family: Helvetica, "Microsoft Yahei", verdana; line-height: 28px;"><span class="" style="font-family: "Times New Roman", serif; font-size: 12pt;"></span></div><div style="font-family: Helvetica, "Microsoft Yahei", verdana; line-height: 28px;"><span class="" style="font-family: "Times New Roman", serif; font-size: 12pt;"><img src="cid:61330a4b$1$17b0547cedd$Coremail$yfang330$163.com" orgwidth="1416" orgheight="1047" data-image="1" style="width: 979.2px; height: 724px; border: 1px solid rgb(153, 153, 153);"></span></div><div style="font-family: Helvetica, "Microsoft Yahei", verdana; line-height: 28px;"><span class="" style="font-family: "Times New Roman", serif; font-size: 12pt;">Fig.2 the force of carbon atom belonging to adsorbed CO molecule along with z axis as the fuction of differential displacement.</span></div><div style="font-family: Helvetica, "Microsoft Yahei", verdana; line-height: 28px;"><span class="" style="font-family: "Times New Roman", serif; font-size: 12pt;"><br></span></div><div style="font-family: Helvetica, "Microsoft Yahei", verdana; line-height: 28px;"><span class="" style="font-family: "Times New Roman", serif; font-size: 12pt;">Best Regards,</span></div><div style="font-family: Helvetica, "Microsoft Yahei", verdana; line-height: 28px;"><span class="" style="font-family: "Times New Roman", serif; font-size: 12pt;"><br></span></div><div style="font-family: Helvetica, "Microsoft Yahei", verdana; line-height: 28px;"><span class="" style="font-family: "Times New Roman", serif; font-size: 12pt;">Yuan Fang</span></div></div><br><br><span title="neteasefooter"><p> </p></span>