Molden2AIM
Features
- It converts the data format from MOLDEN to AIM's WFN. The latter format can be read by AIMPAC, AIMPAC2, AIM2000, AIMALL, AIM-UC, Critic2, DensToolKit, DGrid, IGMPlot, MORPHY98, Multiwfn, ORBKIT, PAMoC, ProMolden, TopChem, TopMoD, Xaim, and so on. The GAB file of Gabedit is compatible.
- It saves NBO's *.47 data file. One can do NBO analysis using the stand-alone GENNBO program.
- After the *.47 file being generated, it can calculate the generalized Wiberg bond order indices (GWBO) in MO (see I. Mayer, C.P.L. 97, 270, 1983). In the case of closed-shell system, they are the Mayer bond orders (MBO) in MO.
- It saves AIM's WFX data file, which can be read by AIMALL, Critic2, DensToolKit, GPView, IGMPlot, Multiwfn, or ORBKIT. An atomic EDF library (Z=3-120) has been included.
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 : 0This 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 : 60Another 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.9 Use 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. Bagel Both spherical and Cartesian functions are supported. 3. BDF-G NR and SR with GTO. Thanks to Dr. Bingbing Suo for testing. 4. CADPAC The utility program cad2mol.f in Molden is required to generate a MOLDEN file. 5. CFour Use 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. Columbus Natural orbitals from MCSCF or CI with Cartesian SPDF functions c). Thanks to Dr. Marat Talipov for testing. 7. CP2k Molecular calculation by Quickstep with GTFs. Both spherical and Cartesian functions are supported. 8. Crystal Only 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. deMon2k Supported. 11. eT (≥ Ver.1.4) Both spherical and Cartesian functions are supported. 12. Gabedit Gabedit 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. Jaguar S, P, D, and F basis functions are supported. You have to insert [Molden Format] into the first line. 14. MOLCAS Both 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. Molden Molden 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=216. MOLPRO Supported. 17. MRCC Both 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. MultiWFN MultiWFN 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. Priroda Thanks to Dr. Evgeniy Pankratyev for testing. 23. PSI4 Both 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. PySCF Supported. 25. Q-Chem Spherical S, P, D, and F or Cartesian S, P, and D basis functions. 26. StoBe Supported. 27. TeraChem S, P, D, and F basis functions are supported. 28. Turbomole Insert a command line
[Program] Turbomole
into the MOLDEN file or specify PROGRAM=3 in the m2a.ini file.29. ADF Not supported. 30. NRLMOL To be tested. 31. SeqQuest To 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)