Molden2AIM

Features

ECP/MCP

In the case of ECP or MCP, a data block of [Core] should be defined in the MOLDEN file. The format is

[Core]
Iatom : Ncore     or
Element: Ncore
...
where Ncore is the number of core electrons replaced by ECP or MCP. Atom/element with Ncore=0 can be ignored.

For example, for a cluster with the atoms N_1, N_2, N_3, Pt_4, and Pt_5, it can be

[Core]
Pt: 60
N : 2
2 : 0

This means that the numbers of core electron are 60 in Pt_4 and Pt_5 and 2 in N_1 and N_3. In N_2 the number of core electron is set to 2 but then reset to 0. It is equivalent to

[Core]
1 : 2
3 : 2
4 : 60
5 : 60

Another way is to define a data block of [Pseudo] in the MOLDEN file, which is supported by Molden. The format is

[Pseudo]
Name1   IAtom1   ZA1-Ncore1
Name2   IAtom2   ZA2-Ncore2
...

Starting from Version 5.1.0, the [Nval] block suggested by MultiWFN may also be used. The format is

[Nval]
Element1 nval1   (here nval is ZA - Ncore)
Element2 nval2
...

Ghost atoms

Ghost atoms in the MOLDEN file may be specified by a prefix "bq-", a prefix "ghost-", a suffix "-bq", a suffix "-ghost", or atomic_number = 0. In the following example, all the five carbon atoms are ghost ones.

[Atoms] AU
C       1   0     0.0000000     2.6361503     0.0000000
C-bq    2   6    -2.2829731     1.3180752     0.0000000
C-ghost 3   6    -2.2829731    -1.3180752     0.0000000
bq-C    4   6     0.0000000    -2.6361503     0.0000000
ghost-C 5   6     2.2829731    -1.3180752     0.0000000
...

About the MOLDEN file

The MOLDEN (or GAB) file can be generated using the following quantum chemistry programs, but some of them have not been supported by Molden2AIM.

1.ACES-II 2.9Use a fixed reorder.F a), and insert a command line
[Program] aces2
into the MOLDEN file or specify PROGRAM=5 in the m2a.ini file.
2.BagelBoth spherical and Cartesian functions are supported.
3.BDF-GNR and SR with GTO. Thanks to Dr. Bingbing Suo for testing.
4.CADPACThe utility program cad2mol.f in Molden is required to generate a MOLDEN file.
5.CFourUse a fixed reorderdf.F b), and insert a command line
[Program] Cfour
into the MOLDEN file or specify PROGRAM=2 in the m2a.ini file.
6.ColumbusNatural orbitals from MCSCF or CI with Cartesian SPDF functions c). Thanks to Dr. Marat Talipov for testing.
7.CP2kMolecular calculation by Quickstep with GTFs. Both spherical and Cartesian functions are supported.
8.CrystalOnly the 0D-molecule calculation with spdf functions is supported; insert a command line
[Program] crystal
into the MOLDEN file or specify PROGRAM=10 in the m2a.ini file.
About how to save a MOLDEN file, see this link.
9.DALTON(> Ver.2013) HF d)/DFT d)/MP2/MCSCF with spherical basis functions.
10.deMon2kSupported.
11.eT(≥ Ver.1.4) Both spherical and Cartesian functions are supported.
12.GabeditGabedit can read geometry, basis functions, and MOs from the output file of Gaussian, Gamess, and Firefly programs, and save a GAB file. However, in the cases of Gamess and Firefly, only s, p, and d basis functions can get correct results. Other QC programs are not supported.
13.JaguarS, P, D, and F basis functions are supported. You have to insert [Molden Format] into the first line.
14.MOLCASBoth spherical and Cartesian functions are supported except their combinations; in the case of Cartesian functions insert a command line
[Program] Molcas
into the MOLDEN file or specify PROGRAM=6 in the m2a.ini file.
15.MoldenMolden can read geometry, basis functions, and MOs from the output file of Gaussian/Gamess/Gamess-UK/Firefly/Q-Chem, and save a standard MOLDEN file via Write/Molden Format. However, in the case of Q-Chem, only s and p basis functions can get correct results.

Gamess/Gamess-UK/Firefly: no special keywords/options are needed

Gaussian: specify the keywords/options GFINPUT IOP(6/7=3)

Q-Chem: use the following keywords/options in $REM
  PRINT_GENERAL_BASIS=TRUE
  PRINT_ORBITALS=TRUE
  SCF_FINAL_PRINT=2

16.MOLPROSupported.
17.MRCCBoth spherical and Cartesian functions are supported. In the case of Cartesian functions insert a command line
[Program] MRCC
into the MOLDEN file or specify PROGRAM=8 in the m2a.ini file.
18.MultiWFNMultiWFN can read fchk file of Gaussian and Q-Chem, and save a MOLDEN file through the options 100/2/6. Until h basis functions are supported.
19.NBO(> May.2014) Both spherical and Cartesian functions are supported. You have to insert a command line
[Program] NBO6
into the MOLDEN file or specify PROGRAM=9 in the m2a.ini file.
20.NWChem(≥ Ver. 6.8) Using MOLDEN_NORM JANPA or NONE to generate a MOLDEN file. See the attached examples.
21.ORCA(Optional) Insert a command line
[Program] Orca
into the MOLDEN file or specify PROGRAM=1 in the m2a.ini file.
* Since Version 5.0.2, Molden2AIM may identify ORCA's MOLDEN file by the title "Molden file created by orca_2mkl for BaseName=...".
22.PrirodaThanks to Dr. Evgeniy Pankratyev for testing.
23.PSI4Both spherical (spdfg) and Cartesian functions (spdf) are supported except their combinations; in the case of Cartesian functions insert a command line
[Program] PSI4
into the MOLDEN file or specify PROGRAM=7 in the m2a.ini file.
24.PySCFSupported.
25.Q-ChemSpherical S, P, D, and F or Cartesian S, P, and D basis functions.
26.StoBeSupported.
27.TeraChemS, P, D, and F basis functions are supported.
28.TurbomoleInsert a command line
[Program] Turbomole
into the MOLDEN file or specify PROGRAM=3 in the m2a.ini file.
29.ADFNot supported.
30.NRLMOLTo be tested.
31.SeqQuestTo be tested.

a) (ACES-II) Put reorder.F into $ACES-II-2.9.0/a2proc, and recompile the ACES-II 2.9 program. To reduce numerical errors, it is also recommended to modify $ACES-II-2.9.0/a2proc/molden_rdorb.F: replace the line 170
  130 format(I4,A,F12.6)
  by
  130 format(I4,A,F18.10)
b) (CFour) For the version 1.x and 2.0beta: put reorderdf.F into $CFour/libr, and recompile the CFour program. Do nothing for 2.1 and higher versions.

c)
(Columbus) Dalton AO integrals have to be used; spherical basis functions cannot pass the examination; occupation numbers of HF and MCSCF MOs are wrong.
d) (Dalton 2013) Occupation numbers of HF and DFT are all zero values by mistake. One has to correct them by hand or by writing a script. No such a problem in Dalton 2015.

Saving a MOLDEN or GAB file using Molden or Gabedit

Limitations

  • In general, only S, P, D, F, and G Gaussian basis functions are supported by Molden2AIM. For the MOLDEN files generated by Multiwfn, ORCA, DALTON, or CFour,  H functions may also be used.
  • The total energy and virial ratio printed at the end of the WFN and WFX files do not make sense. If they are used in your AIM analysis, please modify them manually.
  • The data saved in the *.47 file are very limited.
  • Because of some uncertainties in the MOLDEN format, the atomic indices must be ordered ascendingly.

(2024.03.09)