ghost atoms for cartesian coordinates

Hello,

I'm running the following input file for O2 dimer (quintet state) with midbond functions (MRCC standalone):
# CCSD(T) calculation for the O2 dimer with aug-cc-pVDZ basis, with ghost He atom
basis=aug-cc-pVDZ
calc=CCSD(T)
scftype=ROHF
mult=5
mem=1500MB
symm=0

unit=bohr
geom=xyz
5

O 0.991599087 0.000000000 -2.177500000
O -0.991599087 0.000000000 -3.322500000
O 0.991599087 0.000000000 3.322500000
O -0.991599087 0.000000000 2.177500000
He 0.00000 0.00000 0.0000

ghost=serialno
5

---

But it gives the following :

---

Cartesian coordinates [bohr] Atomic mass [AMU] Z
1 O 0.99159909 0.00000000 -2.17750000 15.994915 8
2 O -0.99159909 0.00000000 -3.32250000 15.994915 8
3 O 0.99159909 0.00000000 3.32250000 15.994915 8
4 O -0.99159909 0.00000000 2.17750000 15.994915 8
5 He 0.00000000 0.00000000 0.00000000 0.000000 0

This molecule is an asymmetric top.
This molecule is planar.

Point group: D2h
Computational point group: C1

Cartesian coordinates in standard orientation [bohr]
1 O NaN NaN 0.00000000
2 O NaN NaN 0.00000000
3 O NaN NaN 0.00000000
4 O NaN NaN 0.00000000
5 He NaN NaN 0.00000000

Nuclear repulsion energy [au]: NaN

---

What is wrong with my input? The test job for HF dimer is working fine, but there Z-matrix is used.

And I have one more question: is it possible to do ROHF/UCCSD(T) calculations with MRCC?
I would be appreciated for help.

best regards, Yulia

Thank you very much for reporting this bug. It has been fixed and a new version is available.

If you run a CCSD(T) calculation with ROHF reference, an ROHF-UCCSD(T) calculation is performed by default.

Mihaly Kallay

Thank you very much!