<div dir="ltr">I think the error was coming because I was doing 'lsint' with ECP. If I remove that command then spin-orbit works fine with ECP MWB.<div>Amit</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jan 11, 2016 at 4:58 PM, Gerald Knizia <span dir="ltr"><<a href="mailto:knizia@theochem.uni-stuttgart.de" target="_blank">knizia@theochem.uni-stuttgart.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Amit,<br>
<br>
Regarding the SO error:<br>
note that statements like "basis=def2-TZVP" and "basis=dhf-TZVP" do not<br>
only select the basis set, but also its associated effective core<br>
potential (in this case the older ECPxxMWB for the def2 variant and the<br>
newer and better(tm) ECPxxMDF for the dhf variant; In fact, this is the<br>
main reason for the existence of the dhf basis sets---some of MDF ECPs<br>
were not yet available while Weigend made the def2 sets, so Weigend made<br>
the dhf ones as an update).<br>
<br>
My guess is that in this case the MDF ECP has also ECP parameters for<br>
the spin-orbit interaction assigned, and the MWB ECP does not. This can<br>
be checked by looking into pseudo.libmol. While the error message could<br>
certainly be more helpful, it does not necessarily mean that there is a<br>
problem in the program.<br>
<br>
Regarding the compilation error: Have not seen that before. Sorry, no<br>
idea.<br>
<br>
Best wishes,<br>
Gerald<br>
<div><div class="h5"><br>
<br>
<br>
On Sun, 2016-01-10 at 12:46 -0500, Amit Sharma wrote:<br>
> When running spin-orbit calculation for an atom I get the following<br>
> error for a particular basis set (def2-TZVP).<br>
><br>
> READM: TRYING TO READ BEYOND RECORD BOUNDARY:<br>
><br>
> RECORD     5551.4  EXTENSION= 1  OFFSET=          0.  ADDRESS=<br>
> 4096.  LEN=         243  LENREST=          24<br>
><br>
><br>
>  Records on file 4<br>
><br>
><br>
>  IREC   NAME  TYPE        OFFSET    LENGTH   IMPLEMENTATION   EXT<br>
> PREV   PARENT  MPP_STATE<br>
><br>
>    1    5551               4096.      219.         df          0<br>
> 0      0      0<br>
><br>
>    2    7771               4315.      438.         df          0<br>
> 0      0      0<br>
><br>
>    3    5552               4753.      219.         df          0<br>
> 0      0      0<br>
><br>
>    4    7772               4972.      438.         df          0<br>
> 0      0      0<br>
><br>
>    5    5553               5410.      219.         df          0<br>
> 0      0      0<br>
><br>
>    6    7773               5629.      438.         df          0<br>
> 0      0      0<br>
><br>
><br>
><br>
>  GLOBAL ERROR fehler on processor   0<br>
><br>
><br>
> This error does not appear with dhf-TZVP basis. partial output<br>
> attached.<br>
><br>
><br>
> My thought was that this is a bug so I tried to install the latest<br>
> 2012 version (2012.1 patch level 32) but during compilation I get the<br>
> following error.<br>
><br>
><br>
><br>
> Rebuilding neci/depend<br>
><br>
> make[3]: Entering directory<br>
> `/usr/local/molpro2012.1.32/Molpro/src/neci'<br>
><br>
> $(FC) -c $(FFLAGS) $(F90FLAGS) -O3 proc_ptrs.F90<br>
><br>
> f951: internal compiler error: in gfc_typenode_for_spec, at<br>
> fortran/trans-types.c:907<br>
><br>
> Please submit a full bug report,<br>
><br>
> with preprocessed source if appropriate.<br>
><br>
> See <<a href="http://bugzilla.redhat.com/bugzilla" rel="noreferrer" target="_blank">http://bugzilla.redhat.com/bugzilla</a>> for instructions.<br>
><br>
> make[3]: *** [proc_ptrs.o] Error 1<br>
><br>
> make[3]: Target `targets' not remade because of errors.<br>
><br>
> make[3]: Leaving directory<br>
> `/usr/local/molpro2012.1.32/Molpro/src/neci'<br>
><br>
> make[2]: *** [default] Error 2<br>
><br>
> make[2]: Leaving directory<br>
> `/usr/local/molpro2012.1.32/Molpro/src/neci'<br>
><br>
> make[1]: *** [neci] Error 2<br>
><br>
> make[1]: Leaving directory `/usr/local/molpro2012.1.32/Molpro/src'<br>
><br>
> make: *** [src] Error 2<br>
><br>
><br>
> Can anyone please point to what I might be doing wrong in the<br>
> spin-orbit calculation and how do I fix the compile error.<br>
><br>
> Thanks<br>
><br>
> Amit<br>
><br>
><br>
> --------------------------------- output<br>
><br>
><br>
> Using customized tuning parameters: mindgm=7; mindgv=24; mindgc=1;<br>
> mindgr=1; noblas=0; minvec=7<br>
><br>
>  default implementation of scratch files=df<br>
><br>
><br>
>  ***,Rb<br>
><br>
><br>
>             !memory,5,M;<br>
><br>
>  gprint,orbitals,civector<br>
><br>
><br>
>             !gthresh,energy=1.d-8,coeff=1.d-8;<br>
><br>
><br>
>  geometry={Rb};<br>
><br>
><br>
><br>
>             !basis=sto-3g<br>
><br>
><br>
>             !basis=dhf-TZVP<br>
><br>
>  basis=def2-TZVP<br>
><br>
><br>
>  {hf; occ,2,1,1,0,1,0,0,0; wf,9,1,1}<br>
><br>
><br>
>  {multi;<br>
><br>
>   occ,2,2,2,0,2,0,0,0;<br>
><br>
>   frozen,0,0,0,0,0,0,0,0;<br>
><br>
>   closed,0,0,0,0,0,0,0,0;<br>
><br>
>   wf,9,1,1,state,1<br>
><br>
>   wf,9,2,1,state,1<br>
><br>
>   wf,9,3,1,state,1<br>
><br>
>   wf,9,5,1,state,1<br>
><br>
>  }<br>
><br>
><br>
><br>
>             ! ci with no-excitations to save casscf orbitals<br>
><br>
>  {ci;wf,9,1,1; noexc; save, 3011.1}<br>
><br>
>  {ci;wf,9,2,1; noexc; save, 3021.1}<br>
><br>
>  {ci;wf,9,3,1; noexc; save, 3031.1}<br>
><br>
>  {ci;wf,9,5,1; noexc; save, 3051.1}<br>
><br>
><br>
>  {ci;wf,9,1,1; save, 4011.1}<br>
><br>
>  {ci;wf,9,2,1; save, 4021.1}<br>
><br>
>  {ci;wf,9,3,1; save, 4031.1}<br>
><br>
>  {ci;wf,9,5,1; save, 4051.1}<br>
><br>
><br>
>  lsint<br>
><br>
><br>
><br>
>             ! determine so matrix<br>
><br>
><br>
>                                                            !{ci;<br>
> hlsmat,ls,3011.1,3021.1,3031.1,3051.1; option,matel=1;print,hls=0}<br>
><br>
>  {ci; hlsmat,ls,4011.1,4021.1,4031.1,4051.1;<br>
> option,matel=1;print,hls=0}<br>
><br>
><br>
><br>
>  Variables initialized (778), CPU time= 0.01 sec<br>
><br>
>  Commands  initialized (547), CPU time= 0.01 sec, 516 directives.<br>
><br>
>  Default parameters read. Elapsed time= 0.06 sec<br>
><br>
><br>
><br>
><br>
><br>
> ------------------------------------- spin orbit calculation<br>
><br>
> 1PROGRAM * CI (Multireference internally contracted CI)     Authors:<br>
> H.-J. Werner, P.J. Knowles, 1987<br>
><br>
><br>
><br>
>    ******************************<br>
><br>
>    *** Spin-orbit calculation ***<br>
><br>
>    ******************************<br>
><br>
><br>
><br>
>    Spin-orbit matrix elements<br>
><br>
>    ==========================<br>
><br>
><br>
>  Preparing effective Fock matrices<br>
><br>
>  Total X Y Z Fock matrices evaluated:     3    3    3<br>
><br>
><br>
>  Wavefunction restored from record  4011.1  Symmetry=1  S= 0.5<br>
> NSTATE=1<br>
><br>
>  Wavefunction restored from record  4021.1  Symmetry=2  S= 0.5<br>
> NSTATE=1<br>
><br>
>  Wavefunction restored from record  4031.1  Symmetry=3  S= 0.5<br>
> NSTATE=1<br>
><br>
>  Wavefunction restored from record  4051.1  Symmetry=5  S= 0.5<br>
> NSTATE=1<br>
><br>
><br>
>  Bra-wavefunction restored from record   4031.1<br>
><br>
>  Ket-wavefunction restored from record   4051.1<br>
><br>
><br>
><br>
>  Symmetry of spin-orbit operator:         7<br>
><br>
>  Symmetry of bra wavefunction:            3   S= 0.5  MS=-0.5<br>
><br>
>  Symmetry of ket wavefunction:            5   S= 0.5  MS= 0.5<br>
><br>
><br>
>  READM: TRYING TO READ BEYOND RECORD BOUNDARY:<br>
><br>
>  RECORD     5551.4  EXTENSION= 1  OFFSET=          0.  ADDRESS=<br>
> 4096.  LEN=         243  LENREST=          24<br>
><br>
><br>
>  Records on file 4<br>
><br>
><br>
>  IREC   NAME  TYPE        OFFSET    LENGTH   IMPLEMENTATION   EXT<br>
> PREV   PARENT  MPP_STATE<br>
><br>
>    1    5551               4096.      219.         df          0<br>
> 0      0      0<br>
><br>
>    2    7771               4315.      438.         df          0<br>
> 0      0      0<br>
><br>
>    3    5552               4753.      219.         df          0<br>
> 0      0      0<br>
><br>
>    4    7772               4972.      438.         df          0<br>
> 0      0      0<br>
><br>
>    5    5553               5410.      219.         df          0<br>
> 0      0      0<br>
><br>
>    6    7773               5629.      438.         df          0<br>
> 0      0      0<br>
><br>
><br>
><br>
>  packet_write_wait: Connection to <a href="http://130.108.150.21" rel="noreferrer" target="_blank">130.108.150.21</a>: Broken pipe<br>
><br>
><br>
><br>
><br>
><br>
><br>
</div></div>> _______________________________________________<br>
> Molpro-user mailing list<br>
> <a href="mailto:Molpro-user@molpro.net">Molpro-user@molpro.net</a><br>
> <a href="http://www.molpro.net/mailman/listinfo/molpro-user" rel="noreferrer" target="_blank">http://www.molpro.net/mailman/listinfo/molpro-user</a><br>
<br>
<br>
</blockquote></div><br></div>