Jeg har udført en CASPT2-frekvensberegning ved hjælp af MOLCAS8-programpakken. Jeg har brug for den fulde kartesiske hessiske matrix som input til en efterfølgende programpakke.
MOLCAS8 udfører CASPT2-frekvensberegninger numerisk. Programmet MCKINLEY
kaldes til dette formål, men det er nu min opfattelse, at alle numeriske sekundære afledte beregninger i MOLCAS håndteres ikke af MCKINLEY
, men af det næppe dokumenterede modul supermac . Under alle omstændigheder er Hessian fuldt udskrevet i interne koordinater, når beregningen er afsluttet. For et system med 10 atomer har den interne Hessian størrelse 24 x 24, hvilket antyder ikke-redundante interne koordinater (MOLCAS kalder dem "primitive interne koordinater").
Er der nogen der kender en software, der er i stand til at transformere en hessisk matrix fra interne koordinater til kartesiske koordinater? Jeg er klar over, at dette sandsynligvis vil være kompliceret, da interne koordinater ikke er så entydige og ligetil som kartesianerne. Min første søgningsindsats gav imidlertid ikke en enkelt softwarepakke, der var i stand til at omdanne en beregnet Hessian i interne koordinater til Cartesians. Enhver hjælp ville blive værdsat!
Alternativt kender nogen en kvantekemisk programpakke, hvor jeg muligvis kan få en fuld kartesisk hessian på CASPT2-niveau?
Jeg ved ikke, om dette hjælper, men MOLCAS ser også ud til at udskrive kraftkonstanterne for hvert atom i kartesiske koordinater. Desuden produceres en .molden
-fil, hvor de normale tilstande er inkluderet. Den fulde Hessian ser ikke ud til at være skrevet til denne Molden-fil.