<div dir="ltr">Hi, XIangyan,<div><br></div><div>I changed the functional type for atoms, from <span style="background-color:rgb(252,248,255);color:rgb(0,0,0);font-family:Menlo;font-size:13px;font-variant-ligatures:no-common-ligatures">GTH-PBE-q4 => </span><span style="background-color:rgb(252,248,255);color:rgb(0,0,0);font-family:Menlo;font-size:13px;font-variant-ligatures:no-common-ligatures">GTH-BLYP-q4, since I am using BLYP. (Sorry, I just copied your basis set and functionals and did not realize the difference.)</span></div><div><span style="background-color:rgb(252,248,255);color:rgb(0,0,0);font-family:Menlo;font-size:13px;font-variant-ligatures:no-common-ligatures"><br></span></div><div><span style="background-color:rgb(252,248,255);color:rgb(0,0,0);font-family:Menlo;font-size:13px;font-variant-ligatures:no-common-ligatures">Here is the result, and is in good agreement with Gaussian 09.</span></div><div><span style="background-color:rgb(252,248,255);color:rgb(0,0,0);font-family:Menlo;font-size:13px;font-variant-ligatures:no-common-ligatures"><br></span></div><div>
<p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>State<span class="gmail-Apple-converted-space"> </span>Excitation<span class="gmail-Apple-converted-space"> </span>Transition dipole (a.u.)<span class="gmail-Apple-converted-space"> </span>Oscillator</span></p>
<p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>number <span class="gmail-Apple-converted-space"> </span>energy (eV) <span class="gmail-Apple-converted-space"> </span>x <span class="gmail-Apple-converted-space"> </span>y <span class="gmail-Apple-converted-space"> </span>z <span class="gmail-Apple-converted-space"> </span>strength (a.u.)</span></p>
<p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>------------------------------------------------------------------------</span></p>
<p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>TDDFPT|<span class="gmail-Apple-converted-space"> </span>1 <span class="gmail-Apple-converted-space"> </span>2.08325<span class="gmail-Apple-converted-space"> </span>-2.4674E+00<span class="gmail-Apple-converted-space"> </span>8.6733E-03 -6.5509E-03 <span class="gmail-Apple-converted-space"> </span>3.10723E-01</span></p>
<p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>TDDFPT|<span class="gmail-Apple-converted-space"> </span>2 <span class="gmail-Apple-converted-space"> </span>2.16407 <span class="gmail-Apple-converted-space"> </span>1.4338E-02<span class="gmail-Apple-converted-space"> </span>2.9818E+00<span class="gmail-Apple-converted-space"> </span>1.0596E-02 <span class="gmail-Apple-converted-space"> </span>4.71397E-01</span></p>
<p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>TDDFPT|<span class="gmail-Apple-converted-space"> </span>3 <span class="gmail-Apple-converted-space"> </span>2.52106 <span class="gmail-Apple-converted-space"> </span>1.3816E-01<span class="gmail-Apple-converted-space"> </span>1.4082E-03 -1.0486E-04 <span class="gmail-Apple-converted-space"> </span>1.17904E-03</span></p>
<p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>TDDFPT|<span class="gmail-Apple-converted-space"> </span>4 <span class="gmail-Apple-converted-space"> </span>2.55291<span class="gmail-Apple-converted-space"> </span>-8.4802E-01<span class="gmail-Apple-converted-space"> </span>4.4219E-03 -1.5266E-04 <span class="gmail-Apple-converted-space"> </span>4.49799E-02</span></p>
<p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>TDDFPT|<span class="gmail-Apple-converted-space"> </span>5 <span class="gmail-Apple-converted-space"> </span>2.58935 <span class="gmail-Apple-converted-space"> </span>1.4471E-04 -4.8114E-02<span class="gmail-Apple-converted-space"> </span>4.4750E-02 <span class="gmail-Apple-converted-space"> </span>2.73894E-04</span></p><p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><br></span></p><p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures">I will test your functional, since it is hybrid, it will take some time. </span></p><p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><br></span></p><p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures">Regards,</span></p><p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><br></span></p><p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures">Fangyong</span></p></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jan 14, 2021 at 6:47 PM Fangyong Yan <<a href="mailto:fyya...@gmail.com">fyya...@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi, XIangyan,<div><br></div><div>I used you structure, and tested both Gaussian and cp2k, and they have good agreement. </div><div><br></div><div>So Gaussian and cp2k can both give reasonable results for the excitation energies of your molecule. </div><div><br></div><div>Gaussian 09, blyp/6-31G*, </div><div>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures">Excited State <span> </span>1:<span> </span>Singlet-A<span> </span>1.9414 eV<span> </span>638.63 nm<span> </span>f=0.2542<span> </span><S**2>=0.000</span></p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures">
</span></p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures">Excited State <span> </span>2:<span> </span>Singlet-A<span> </span>1.9724 eV<span> </span>628.58 nm<span> </span>f=0.3319<span> </span><S**2>=0.000</span></p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures">
</span></p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures">Excited State <span> </span>3:<span> </span>Singlet-A<span> </span>2.5113 eV<span> </span>493.71 nm<span> </span>f=0.0038<span> </span><S**2>=0.000</span></p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures">
</span></p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures">Excited State <span> </span>4:<span> </span>Singlet-A<span> </span>2.5205 eV<span> </span>491.90 nm<span> </span>f=0.0341<span> </span><S**2>=0.000</span></p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures">
</span></p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures">Excited State <span> </span>5:<span> </span>Singlet-A<span> </span>2.5764 eV<span> </span>481.23 nm<span> </span>f=0.0001<span> </span><S**2>=0.000</span></p></div><div><br></div><div>cp2k, blyp/TZVP, with GTH pseudopotential, </div><div><br></div><div>I get</div><div><br></div><div>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>State<span> </span>Excitation<span> </span>Transition dipole (a.u.)<span> </span>Oscillator</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>number <span> </span>energy (eV) <span> </span>x <span> </span>y <span> </span>z <span> </span>strength (a.u.)</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>------------------------------------------------------------------------</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>TDDFPT|<span> </span>1 <span> </span>2.08760 <span> </span>2.4868E+00 -8.6129E-03<span> </span>6.4653E-03 <span> </span>3.16297E-01</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>TDDFPT|<span> </span>2 <span> </span>2.16739<span> </span>-1.4224E-02 -2.9924E+00 -1.0393E-02 <span> </span>4.75501E-01</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>TDDFPT|<span> </span>3 <span> </span>2.52956<span> </span>-1.3916E-01 -1.4062E-03<span> </span>1.0520E-04 <span> </span>1.20031E-03</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>TDDFPT|<span> </span>4 <span> </span>2.56128 <span> </span>8.4944E-01 -4.4270E-03<span> </span>1.5135E-04 <span> </span>4.52790E-02</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>TDDFPT|<span> </span>5 <span> </span>2.59590<span> </span>-1.3595E-04<span> </span>4.4904E-02 -4.4546E-02 <span> </span>2.54439E-04</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>TDDFPT|<span> </span>6 <span> </span>2.60688<span> </span>-3.1839E-03<span> </span>4.1336E-02<span> </span>2.5674E-02 <span> </span>1.51876E-04</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>TDDFPT|<span> </span>7 <span> </span>2.76485<span> </span>-3.0533E-03 -2.4437E-01<span> </span>1.0913E-03 <span> </span>4.04581E-03</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>TDDFPT|<span> </span>8 <span> </span>2.80552<span> </span>-2.7107E-01 -4.8867E-03 -9.2586E-04 <span> </span>5.05222E-03</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>TDDFPT|<span> </span>9 <span> </span>2.88205<span> </span>-1.2938E-02 -1.3849E-04<span> </span>1.2856E-03 <span> </span>1.19366E-05</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>TDDFPT| <span> </span>10 <span> </span>2.92954 <span> </span>2.8231E-05<span> </span>6.4011E-03 -8.3886E-02 <span> </span>5.07991E-04</span></p></div><div><br></div><div>****************** Below is the input. </div><div>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&SUBSYS</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&CELL</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>ABC<span> </span>30. 30. 6.</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>PERIODIC none</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&END CELL</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&COORD</span></p></div><div><br></div><div>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&KIND C</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>BASIS_SET TZVP-GTH</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>POTENTIAL GTH-PBE-q4</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&END KIND</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&KIND N</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>BASIS_SET TZVP-GTH</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>POTENTIAL GTH-PBE-q5</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&END KIND</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&KIND H</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>BASIS_SET TZVP-GTH</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>POTENTIAL GTH-PBE-q1</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&END KIND</span></p></div><div><br></div><div>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&QS</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>METHOD GPW</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>MAP_CONSISTENT<span> </span>T</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>EPS_DEFAULT 1.0E-9</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&END QS</span></p></div><div><br></div><div>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&XC</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&XC_FUNCTIONAL BLYP</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&END XC_FUNCTIONAL</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255);min-height:15px"><span style="font-variant-ligatures:no-common-ligatures"></span><br></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&XC_GRID</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>XC_DERIV SPLINE2_SMOOTH<span> </span># this is needed for the 2nd derivatives of the XC functional</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&END XC_GRID</span></p></div><div>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&END XC</span></p></div><div>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&SCF</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>MAX_SCF <span> </span>250</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>EPS_SCF <span> </span>1e-8</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>SCF_GUESS atomic</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&END SCF</span></p></div><div><br></div><div>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&MGRID</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>CUTOFF 600</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>NGRIDS 4</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>REL_CUTOFF 60</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&END MGRID</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255);min-height:15px"><span style="font-variant-ligatures:no-common-ligatures"></span><br></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&</span><span style="font-variant-ligatures:no-common-ligatures;color:rgb(252,248,255);background-color:rgb(0,0,0)">POI</span><span style="font-variant-ligatures:no-common-ligatures">SSON</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>PERIODIC none</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>POISSON_SOLVER MT</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&END POISSON</span></p></div><div><br></div><div>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&PROPERTIES</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&TDDFPT</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>NSTATES <span> </span>10 <span> </span># number of excited states</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>MAX_ITER<span> </span>100 <span> </span># maximum number of Davidson iterations</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>CONVERGENCE [eV] 1.0e-3<span> </span># convergence on maximum energy change between iterations</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255);min-height:15px"><span style="font-variant-ligatures:no-common-ligatures"></span><br></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&MGRID</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>CUTOFF 300<span> </span># separate cutoff for TDDFPT calc</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&END</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&END TDDFPT</span></p>
<p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><span> </span>&END PROPERTIES</span></p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><br></span></p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures">Regares,</span></p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures"><br></span></p><p style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(252,248,255)"><span style="font-variant-ligatures:no-common-ligatures">Fangyong</span></p></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jan 12, 2021 at 6:50 PM Lucas Lodeiro <<a href="mailto:eluni...@gmail.com" target="_blank">eluni...@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi Liu,<div><br></div><div>As you mention it is a little weird if other molecules work fine with the same input. I guess the next test is to tighten the convergence criteria in CP2K if it is possible.</div><div>Another thing, remember that for non periodic calculations, you have to set periodic none in the poisson and cell section.</div><div><br></div><div><div><b> &CELL</b></div><div><b> ABC 30.0 30.0 30.0</b></div><div><b> PERIODIC NONE</b></div><div><b> &END CELL</b></div></div><div><br></div><div>I see the transition energies change a little bit, it is at least comfortable, showing that one or more changes you did affect the result. Maybe as the last chance there will be useful to know which change was the one that caused the variation a the energies.</div><div><br></div><div>Sorry, it does not seem a simple problem.</div><div>Regards</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El mar, 12 ene 2021 a las 1:49, liu xiangyang (<<a href="mailto:lxyl...@gmail.com" target="_blank">lxyl...@gmail.com</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi
Vladimir & Lucas，<div><br></div><div>Thanks a lot for your kindly responses. Actually, I have also tested small molecule such as H2CO using present settings and the results obtained by CP2K and G09 are quanlitatively in good agreement with each other and THAT IS THE FACT THAT PUZZLES ME MOST.</div><div><br></div><div>According to your suggestions, I have made several changes to my input file, including removal of the periodic conditions, removal of the ADMM approximations to avoid possible mistakes, and lower the convergence in G09 since the computational efforts of CP2K is much larger. The basis sets are not changed since the TDDFPT in CP2K is implemented only in GPW method, which is unable to do calculations with all electron basis sets such as Pople basis sets. However, I believe the difference of basis set (DZVP-GTH vs. 6-31G*) should not be so large. Unfortunately, the results of CP2K7.1 is still about 1 eV lower than that obtained in G09 (1.19 and 1.25 eV in CP2K vs. 2.18 and 2.26 eV in G09). </div><div><br></div><div>PS. I have also run GAPW calculations in CP2K using the WB97XD/6-31G** and the resulted ground state energy are nearly the same as that obtained in G09 (differenceless than 0.002 hartree). </div><div><br></div><div>The modified DFT parts of CP2K input file is attatched also:</div><div><br></div><div><div>&DFT</div><div> BASIS_SET_FILE_NAME BASIS_MOLOPT</div><div> POTENTIAL_FILE_NAME POTENTIAL</div><div><br></div><div> CHARGE 0 </div><div> &MGRID</div><div> CUTOFF 400</div><div> &END MGRID</div><div> &QS</div><div> METHOD gpw</div><div> EPS_PGF_ORB 1e-12</div><div> &END QS</div><div> &SCF</div><div> MAX_SCF 100</div><div> EPS_SCF 1e-5</div><div> SCF_GUESS atomic</div><div><br></div><div> &DIAGONALIZATION</div><div> ALGORITHM STANDARD</div><div> &END DIAGONALIZATION</div><div><br></div><div> &MIXING T</div><div> ALPHA 0.5</div><div> METHOD PULAY_MIXING</div><div> NPULAY 5</div><div> &END MIXING</div><div> &END SCF</div><div> &POISSON</div><div> PERIODIC NONE</div><div> PSOLVER MT</div><div> &END POISSON</div><div><br></div><div> &XC</div><div> &XC_FUNCTIONAL</div><div> &LIBXC</div><div> FUNCTIONAL HYB_GGA_XC_WB97X_D</div><div> &END</div><div> &END XC_FUNCTIONAL</div><div> &HF</div><div> &SCREENING</div><div> EPS_SCHWARZ 1.0E-6</div><div> &END</div><div> &MEMORY</div><div> MAX_MEMORY 100</div><div> &END</div><div> &INTERACTION_POTENTIAL</div><div> POTENTIAL_TYPE MIX_CL</div><div> OMEGA 0.2</div><div> SCALE_LONGRANGE 0.777964</div><div> SCALE_COULOMB 0.222036</div><div> &END<br></div><div> &END</div><div> &END XC</div><div><br></div><div> &END DFT</div></div><div><br></div><div>Best wishes,</div><div>Xiang-Yang Liu</div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, January 12, 2021 at 7:07:34 AM UTC+8 <a href="mailto:ry...@gmail.com" target="_blank">ry...@gmail.com</a> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Sorry, I really read the differences incorrectly. Thank you, Lukas, for correcting me. <br><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">понедельник, 11 января 2021 г. в 21:02:30 UTC+1, Lucas Lodeiro: <br></div></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi Liu,<div>I did not run TDDFT calculations, but I did some tests between CP2K and other programs like G09. As Vladimir says, your basis sets are not the same, and some difference could appear due to this reason. But in your case the differences are big, 1.6 eV approx. I found that some default settings of convergences criterium are differents, for example the <b>EPS_SCF</b> which is 1E-8 in G09, you could tight your convergence criterion, <b>EPS_SCF, </b><b>EPS_DEFAULT, </b><b>EPS_SCHWARZ</b> to -8, -12 and -8 to get results with similar convergences in both programs. Also, you are using a PBC calculation in a big cell, but maybe it is no sufficient to mimic the isolated molecule as in G09... and yout cutoff radius for HF is a little bit short, if you run a non-periodic calculation, you can use just the long range potential without the truncation.</div><div>Finally, just to speed up, you can use OT instead of diagonalization method, with it you can use <b>ADMM_PURIFICATION_METHOD </b><span style="font-size:medium;color:rgb(0,0,0);font-family:monospace;text-transform:uppercase">MO_DIAG</span>.</div><div><br></div><div>In order to have the same basis sets, as vladimir says, you could explore to use the same basis sets in both programs, you can get basis sets to both programs from: <a href="https://www.basissetexchange.org/" rel="nofollow" target="_blank">https://www.basissetexchange.org/</a></div><div>And the Auxiliary basis set is the minimum one, you could explore FIT and pFIT basis set to check if the result is sensitive.</div><div><br></div><div>Regards</div><div>Lucas Lodeiro</div></div><br><div class="gmail_quote"></div><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El lun, 11 ene 2021 a las 13:01, Vladimir Rybkin (<<a rel="nofollow">ry...@gmail.com</a>>) escribió:<br></div></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Xiang-Yang Liu,<div><br></div><div>most importantly: you are using different basis sets in Gaussian and CP2K (CP2K also use pseudopotentials). With this difference in mind you differences are within reasonable. Generally, difference below 0.1 eV for TDDFT implementations is not "great".</div><div><br></div><div>Yours,</div><div><br></div><div>Vladimir<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">понедельник, 11 января 2021 г. в 14:49:47 UTC+1, <a rel="nofollow">lx...@gmail.com</a>: <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div>Dear All,</div><div><br></div><div>I have tried to use the TDDFPT method implemented in CP2K7.1 to do excited state calculations with the range-separated functionals such as wB97XD. </div><div>However, after several tests with a small molecule, namely H2Pc, I found that the first two excitation energies are greatly underestimated in comparison with the LR-TDDFT results obtained in GAUSSIAN09 (ca. 0.51 and 0.60 eV (CP2K7.1) vs. 2.18 and 2.26 eV (GAUSSIAN09)).</div><div>I wonder whether there are some mistakes with my input file or there are some problem of TDDFPT for such calculations?</div><div><br></div><div>The input file used in my calculations is written as follows:</div><div><br></div><div><div><b>&GLOBAL</b></div><div><b> PROJECT tddfpt</b></div><div><b> RUN_TYPE energy</b></div><div><b> PRINT_LEVEL medium</b></div><div><b>&END GLOBAL</b></div><div><b><br></b></div><div><b>&FORCE_EVAL</b></div><div><b> METHOD Quickstep</b></div><div><b><br></b></div><div><b> &PROPERTIES</b></div><div><b> &TDDFPT</b></div><div><b> NSTATES 5 # number of excited states</b></div><div><b> &END TDDFPT</b></div><div><b> &END PROPERTIES</b></div><div><b><br></b></div><div><b>&DFT</b></div><div><b> BASIS_SET_FILE_NAME GTH_BASIS_SETS</b></div><div><b> POTENTIAL_FILE_NAME POTENTIAL</b></div><div><b> BASIS_SET_FILE_NAME BASIS_ADMM_MOLOPT</b></div><div><b> BASIS_SET_FILE_NAME BASIS_ADMM</b></div><div><b><br></b></div><div><b> &AUXILIARY_DENSITY_MATRIX_METHOD</b></div><div><b> METHOD BASIS_PROJECTION</b></div><div><b> ADMM_PURIFICATION_METHOD NONE </b></div><div><b> EXCH_CORRECTION_FUNC BECKE88X</b></div><div><b> &END AUXILIARY_DENSITY_MATRIX_METHOD</b></div><div><b><br></b></div><div><b> CHARGE 0 </b></div><div><b> &MGRID</b></div><div><b> CUTOFF 400</b></div><div><b> &END MGRID</b></div><div><b> &QS</b></div><div><b> METHOD gpw</b></div><div><b> EPS_PGF_ORB 1e-12</b></div><div><b> &END QS</b></div><div><b> &SCF</b></div><div><b> MAX_SCF 100</b></div><div><b> EPS_SCF 1e-5</b></div><div><b> SCF_GUESS atomic</b></div><div><b><br></b></div><div><b> &DIAGONALIZATION</b></div><div><b> ALGORITHM STANDARD</b></div><div><b> &END DIAGONALIZATION</b></div><div><b><br></b></div><div><b> &MIXING T</b></div><div><b> ALPHA 0.5</b></div><div><b> METHOD PULAY_MIXING</b></div><div><b> NPULAY 5</b></div><div><b> &END MIXING</b></div><div><b> &END SCF</b></div><div><b> &POISSON</b></div><div><b> PERIODIC XYZ</b></div><div><b> &END POISSON</b></div><div><b><br></b></div><div><b> &XC</b></div><div><b> &XC_FUNCTIONAL</b></div><div><b> &LIBXC</b></div><div><b> FUNCTIONAL HYB_GGA_XC_WB97X_D</b></div><div><b> &END</b></div><div><b> &END XC_FUNCTIONAL</b></div><div><b> &HF</b></div><div><b> &SCREENING</b></div><div><b> EPS_SCHWARZ 1.0E-6</b></div><div><b> &END</b></div><div><b> &MEMORY</b></div><div><b> MAX_MEMORY 100</b></div><div><b> &END</b></div><div><b> &INTERACTION_POTENTIAL</b></div><div><b> POTENTIAL_TYPE MIX_CL_TRUNC</b></div><div><b> OMEGA 0.2</b></div><div><b> SCALE_LONGRANGE 0.777964</b></div><div><b> SCALE_COULOMB 0.222036</b></div><div><b> CUTOFF_RADIUS 5.0</b></div><div><b> T_C_G_DATA t_c_g.dat</b></div><div><b> &END</b></div><div><b> &END</b></div><div><b> &XC_GRID</b></div><div><b> XC_DERIV SPLINE2_SMOOTH # this is needed for the 2nd derivatives of the XC functional</b></div><div><b> &END XC_GRID</b></div><div><b> &END XC</b></div><div><b><br></b></div><div><b> &END DFT</b></div><div><b> &SUBSYS</b></div><div><b> &TOPOLOGY</b></div><div><b> &CENTER_COORDINATES</b></div><div><b> &END CENTER_COORDINATES</b></div><div><b> &END TOPOLOGY</b></div><div><b> &CELL</b></div><div><b> ABC 30.0 30.0 30.0</b></div><div><b> PERIODIC XYZ</b></div><div><b> &END CELL</b></div><div><b> &COORD</b></div><div><b> N 0.01468800 2.00544200 -0.29683900</b></div><div><b> N -2.37562900 2.41155000 -0.13186300</b></div><div><b> N -2.41200400 -2.37377100 -0.17491600</b></div><div><b> N 2.41112400 2.37299900 -0.12658700</b></div><div><b> N -0.01307800 -2.00591400 -0.34375500</b></div><div><b> N 2.37788800 -2.41418900 -0.17065800</b></div><div><b> C -4.17222800 0.73645800 -0.00309900</b></div><div><b> C -0.67412900 4.18239500 0.01934000</b></div><div><b> C -4.18380100 -0.67562700 -0.02014000</b></div><div><b> C 0.73808000 4.17070300 0.02178600</b></div><div><b> C -2.78102300 1.14403300 -0.18692600</b></div><div><b> C -1.10192300 2.79691600 -0.16605600</b></div><div><b> C -2.79906300 -1.10173400 -0.21220500</b></div><div><b> C 1.14375300 2.77826000 -0.16260500</b></div><div><b> C -0.73890800 -4.16630900 -0.00019300</b></div><div><b> C 4.18311400 0.67096100 0.00115100</b></div><div><b> C 0.67355900 -4.17905700 0.00226800</b></div><div><b> C 4.17314700 -0.74140800 -0.01277100</b></div><div><b> C -1.14329800 -2.77792600 -0.20621200</b></div><div><b> C 2.79688300 1.09912600 -0.17512500</b></div><div><b> C 1.10326800 -2.79777700 -0.20275300</b></div><div><b> C 2.78178900 -1.14729500 -0.19912900</b></div><div><b> C -5.34524900 1.46617200 0.20504700</b></div><div><b> C -1.38645100 5.36877700 0.20954200</b></div><div><b> C -5.36863500 -1.39063500 0.17089500</b></div><div><b> C 1.47037200 5.34443400 0.21437700</b></div><div><b> C -6.53049700 0.75063900 0.38667600</b></div><div><b> C -0.65342400 6.54324100 0.39406800</b></div><div><b> C -6.54209900 -0.66018000 0.36971000</b></div><div><b> C 0.75719100 6.53149500 0.39645600</b></div><div><b> C -1.47251300 -5.33425600 0.22345700</b></div><div><b> C 5.36943600 1.38205000 0.19903500</b></div><div><b> C 1.38473700 -5.36039200 0.22884400</b></div><div><b> C 5.34798200 -1.47459500 0.17084800</b></div><div><b> C -0.76071300 -6.51537100 0.44153000</b></div><div><b> C 6.54461400 0.64821000 0.37381100</b></div><div><b> C 0.65044900 -6.52827600 0.44412700</b></div><div><b> C 6.53422100 -0.76263800 0.36007100</b></div><div><b> H -5.31982100 2.55648400 0.22741200</b></div><div><b> H -2.47726000 5.36199000 0.21499900</b></div><div><b> H -5.36047700 -2.48141400 0.16783300</b></div><div><b> H 2.56090800 5.31878300 0.22274500</b></div><div><b> H -7.46690300 1.28881400 0.54987000</b></div><div><b> H -1.17840000 7.48942700 0.54335800</b></div><div><b> H -7.48727000 -1.18662900 0.52012100</b></div><div><b> H 1.29740200 7.46875700 0.54746900</b></div><div><b> H -2.56294400 -5.30668600 0.23110100</b></div><div><b> H 5.36158400 2.47267500 0.21857200</b></div><div><b> H 2.47541900 -5.35281000 0.24119000</b></div><div><b> H 5.32277600 -2.56515200 0.16738300</b></div><div><b> H -1.30192600 -7.44735900 0.61932200</b></div><div><b> H 7.49034300 1.17231500 0.52877000</b></div><div><b> H 1.17373400 -7.47001200 0.62407400</b></div><div><b> H 7.47217800 -1.30348900 0.50426500</b></div><div><b> H -0.95410508 0.01210298 -0.52371821</b></div><div><b> H 0.94993522 -0.00596112 -0.48792422</b></div><div><b> N -2.00779500 0.01729500 -0.33767900</b></div><div><b> N 2.00585900 -0.01636500 -0.31524900</b></div><div><b> &END COORD</b></div><div><b> &KIND N</b></div><div><b> BASIS_SET DZVP-GTH</b></div><div><b> POTENTIAL GTH-PBE</b></div><div><b> BASIS_SET AUX_FIT cFIT3</b></div><div><b> &END KIND</b></div><div><b> &KIND C</b></div><div><b> BASIS_SET DZVP-GTH</b></div><div><b> POTENTIAL GTH-PBE</b></div><div><b> BASIS_SET AUX_FIT cFIT3</b></div><div><b> &END KIND</b></div><div><b> &KIND H</b></div><div><b> BASIS_SET DZVP-GTH</b></div><div><b> POTENTIAL GTH-PBE</b></div><div><b> BASIS_SET AUX_FIT cFIT3</b></div><div><b> &END KIND</b></div><div><b> &END SUBSYS</b></div><div><b>&END FORCE_EVAL</b></div><div><b><br></b></div><div><b>Best wishes, </b></div><div><b>Xiang-Yang Liu</b></div></div></blockquote></div>
<p></p></blockquote></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a rel="nofollow">cp...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/4897d4d1-a074-44f6-89ff-3b2c613eab83n%40googlegroups.com?utm_medium=email&utm_source=footer" rel="nofollow" target="_blank">https://groups.google.com/d/msgid/cp2k/4897d4d1-a074-44f6-89ff-3b2c613eab83n%40googlegroups.com</a>.<br>
</blockquote></div>
</blockquote></div></blockquote></div>
<p></p>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/64940b8a-9bef-4f42-9833-55b0ddaa7326n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/64940b8a-9bef-4f42-9833-55b0ddaa7326n%40googlegroups.com</a>.<br>
</blockquote></div>
<p></p>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAOFT4PJ1bAsK-iK5Lnd7YikGQjJYgDxb47dNh5bNH_FeM1U%2BrQ%40mail.gmail.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/CAOFT4PJ1bAsK-iK5Lnd7YikGQjJYgDxb47dNh5bNH_FeM1U%2BrQ%40mail.gmail.com</a>.<br>
</blockquote></div>
</blockquote></div>