C2対称性の時IVO計算で出力されるDFPCMONEWの不具合

[minoria-hiroshima]

どういうバグ？

DFPCMONEWのサイズが、DFPCMOよりも小さく、ファイルの最後の終わり方も不自然。(CS2)2 + のC2対称性の時に起こっている。同じ分子でもC2hの時にはうまく流れた。（前も同じ問題があった？C2はIVOが未対応？）

Expected behavior

DFPCMONEWがもとのDFPCMOと同じ量になり、moltra after ivoが流せる。

再現性

再現性がある場合、再現するための手順を記述してください

環境

ビルド環境を選択、記述してください

• gfortran
• Intel Fortran
• mkl

DIRACのバージョンなどの情報は以下に書いてください
(例) DIRAC ver 22.0
parallel build

(optional) inputs

計算時のインプット等について記述してください

mol ファイル

DIRAC
cs2+ S' Muramatsu sensei
***
C 2 1
16. 4
S1 1.801983829 -1.797030661 -1.829205757
S3 -1.801983829 1.797030661 -1.829205757
S2 -0.693902403 -5.708329981 1.815312555
S4 0.693902403 5.708329981 1.815312555
LARGE BASIS aug-cc-pVTZ
6. 2
C1 0.511625239 -3.770967216 0.037084012
C2 -0.511625239 3.770967216 0.037084012
LARGE BASIS aug-cc-pVTZ
FINISH

x2c.inp

**DIRAC
.TITLE
CS2+ S' Muramatsu
.WAVE FUNCTION
!.ANALYZE
.4INDEX
**GENERAL
.PCMOUT
**HAMILTONIAN
.X2C
*AMFI
.MXITER
100
!.AMFICH
!+3
**INTEGRALS
.NUCMOD
2
*READIN
.UNCONTRACT
**WAVE FUNCTIONS
.SCF
*SCF
!.CLOSED SHELL
!316
!.OPEN SHELL
!1
!5/10
.MAXITR
300
.ERGCNV
1.0E-10
**MOLTRA
.ACTIVE
23..314
! up to 250 a.u. orbitals
*END OF

active.ivo.inp

ninact
0
nact
32
nsec
552
nelec
32
nocc
38
nvcut
102
nroot
1
selectroot
1
totsym
1
diracver
22
end

75 electrons system
44 electrons are frozen.
76 - 44 = 32 active spinors for ivo

(optional) Screenshots

If applicable, add screenshots to help explain your problem.

(optional) Additional context

Add any other context about the problem here.

[minoria-hiroshima]

前も野田さんによって同じ問題が指摘されており、確か、UX4のようなTdをC2にする場合には大丈夫だったかが、直線分子をわざとC2にさせるとDIRACのDFPCMOの出力が予想よりも多くなっていた。DIRACの対称性のところが理解できていないので、ここを理解して解決する必要がある。同様にC1分子の時も。

[kohei-noda-qcrg]

@minoria-hiroshima

https://github.com/RQC-HU/dirac_caspt2/tree/dfpcmo-coefs-cmplx
上記ブランチでDFPCMOの出力が多い(cmplxのとき)
の対応をやろうとしたことがあります

ただし現時点では
正しいDFPCMONEWを作ることはできていません

[minoria-hiroshima]

@kohei-noda-qcrg
ブランチ名の記述ありがとうございます。