Hello,
<br>
<br>I am trying to save a transition density matrix from an EOM-CCSD calculation.
<br>If I am not mistaken, according to the manual, this can be done by setting trans=2
<br>and choosing a record, either with the dm directive or the densave keyword.
<br>
<br>The problem is that I can't get any one of these to work.
<br>
<br>With densave I get the following error:
<br>
<br>********************************************************************************
<br>
<br> EOM-CCSD properties program.                 Author: T. Korona (2003)
<br>
<br>********************************************************************************
<br>
<br>
<br> Records on file 2
<br>
<br> IREC   NAME  TYPE        OFFSET    LENGTH   IMPLEMENTATION   EXT   PREV   PARENT  MPP_STATE
<br>   1     500  VAR          4096.   102002.         df          0      0      0      1
<br>   2     610  BASINP     106098.     8192.         df          0      0      0      1
<br>   3     700  GEOM       114290.    27110.         df          0      0      0      1
<br>   4    1000  BASIS      141400.      933.         df          0      0      0      1
<br>   5    2100  RHF        142333.     7552.         df          0      0      0      1
<br>   6    5555             149885.  2527727.         df          0      0      0      1
<br>   7    5040            2677612.      300.         df          0      0      0      1
<br>   8    5700            2677912.    10624.         df          0      0      0      1
<br>
<br> READM:  RECORD=    5700.2  EXTENSION= 0  OFFSET=          0.  ADDRESS=    2677912.  LEN=          20  IMPLEMENTATION=df
<br> CURRENT WRITTEN FILE LENGTH:  2677912.TRYING TO READ UP TO ADDRESS:  2677932.
<br> ? Error
<br> ? Trying to read beyond end of file
<br> ? The problem occurs in readm
<br>
<br> ERROR EXIT
<br> CURRENT STACK:      MAIN
<br>
<br>
<br>If I use dm I get a different one:
<br>
<br> Natural orbitals saved on record    5600.2 (orbital set 2)
<br>
<br> Bra <   1.1|,       Ket |  2.1>
<br>
<br> EOM-CCSD trans. dipole moment          0.00000000     0.00000000     0.36439832   <  1.1|  2.1>
<br>
<br>
<br> ?ERROR IN WRITE_DUMP: OVERWRRITING RECORD    5600.2  ENTRY DENSITY/TRANSITION SET= 2 WITH ILLEGAL LENGTH. NEW=     577 
OLD=     300
<br>
<br> ERROR EXIT
<br> CURRENT STACK:      MAIN
<br>
<br>
<br>A sample input file is below:
<br>
<br>memory,100,m
<br>file,2,w_tdm.wfu
<br>symmetry,nosym
<br>geometry
<br>3
<br>
<br>O         -0.32532        0.22659        0.00000
<br>H          0.62464        0.21785        0.00000
<br>H         -0.65067       -0.66595        0.00000
<br>endg
<br>basis=vdz;hf;{ccsd;eom,2.1,trans=2,densave=5600.2;}
<br>
<br>
<br>Am I missing something in my input or is there a bug in the program?
<br>
<br>The molpro version used is 2010.1 patch level 20 (latest).
<br>
<br>Thank you,
<br>Hugo