- Calculation of Parity-Violating Nuclear Spin-Rotation Tensors. Activated with “.PVCSR”
Contributor: I. Agustín Aucar.
Reference: Ignacio Agustín Aucar and Anastasia Borschevsky, Relativistic study of parity-violating nuclear spin-rotation tensors. J. Chem. Phys. 155, 134307 (2021)
- Calculation of oscillator strenghts beyond the electric dipole approximation in the generalized velocity representation
Contributors: Martin van Horn, Trond Saue, and Nanna Holmgaard List.
Reference: `M. van Horn, T. Saue, N. H. List. Probing Chirality across the Electromagnetic Spectrum with the Full Semi-Classical Light-Matter Interaction. J. Chem. Phys. 156, 054113 (2022) < https://doi.org/10.1063/5.0077502>`__
- (extended) atomic-mean-field two-electron scalar and spin-orbit picture-change corrections for two-component Hamiltonian calculations
Contributors: Stefan Knecht, Michal Repsiky, Hans Joergen Aa. Jensen and Trond Saue.
Reference: Stefan Knecht, Michal Repisky, Hans Jørgen Aagaard Jensen, and Trond Saue, Exact two-component Hamiltonians for relativistic quantum chemistry: Two-electron picture-change corrections made simple, J. Chem. Phys. 157, 114106 (2022)
- ExaCorr module works now also with contracted basis functions.
Contributor: Lucas Visscher.
- Basis set library has been extended with basis sets with diffuse functions for s and d blocks, dyall.aXNz basis sets.
Infrastructure changes and fixes:
- HDF5 checkpoint file is now default for data handling
Contributor: Lucas Visscher
User manual: Data storage in DIRAC
Developer manual: The CHECKPOINT.h5 file
- Added a CI gradient convergence threshold option ‘.THRGCI’ to ‘*KRCICALC’ (could not be changed by user previously).
Contributor: H. J. Aa. Jensen.
- Improved code in luciarel CI program. Note that KRMC and KRCI by default uses luciarel.
Contributor: Andreas Nyvang.
numerical stability has been improved in the complex parallel diagonalization routine
convergence algorithm has also been changed, with considerable improvement in convergence rate. In particular, the truncation of the reduced space in the Davidson iterations is now done in a much better way, leading to the improvement in convergence rate. Both for real and complex parallel cases.
Default for .MVO has been corrected, so it now behaves as expected: modified virtual orbitals based on the Fock potential from the doubly-occupied molecular orbitals after an open-shell SCF calculation (H. J. Aa. Jensen).
New licence: GNU Lesser General Public License v2.1
- Electronic circular dichroism (ECD) beyond lowest-order – full light-matter interaction
Contributors: Martin van Horn, Trond Saue, and Nanna Holmgaard List
Reference: M. van Horn, T. Saue, N. H. List., Probing Chirality across the Electromagnetic Spectrum with the Full Semiclassical Light-Matter Interaction. J. Chem. Phys. (in press) (2022), ChemRxiv
- MP2 frozen virtual natural orbitals via the ExaCorr module
Contributors: Xiang Yuan, Lucas Visscher, André Severo Pereira Gomes.
Reference: X. Yuan, L. Visscher, A. S. P. Gomes., Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations.
- Polarizable Embedding with Complex Polarization Propagator
Contributors: Joel Creutzberg, Erik D. Hedegård
Reference: Polarizable embedding complex polarization propagator in four- and two-component frameworks.
Infrastructure changes and fixes:
- HDF5 checkpoint file and DIRAC data scheme, python utilities to extract data
Contributor: Lucas Visscher
Manual: The CHECKPOINT.h5 file
- DIIS replaced by CROP algorithm in ExaCorr module.
Contributors: Chima Chibueze and Lucas Visscher.
Reference: Pototschnig et al., J. Chem. Theory Comput. 17 (2021) 5509−5529.
There was an error in the construction of the projection operator for secondary orbitals using density matrices for the very special case of 4-component relatistic average-of-configuration HF calculations on systems with inversion symmetry and no occupation in ungerade symmetry.
Enhanced pam script: now it also looks in job directory for .diracrc, and it looks there first. pam looks for .diracrc three places: 1) job directory (“local .diracrc”), 2) home directory (“global .diracrc”), 3) pam script directory (“system .diracrc). The first .diracrc found is used. (H. J. Aa. Jensen). - added ENABLE_EXATENSOR to the options you can change with ‘ccmake –’ (H. J. Aa. Jensen). - enable talsh only with unsupported mpi libraries (Andre Gomes) - improved exatensor MPI identification (Johann Pototschnig)
Fixed pam script so it can find dirac.x and basis set directories in an installed version (i.e. after make install or e.g. “cmake –install build –prefix <installation root dir>”). H.J.Aa.Jensen.
Fixed memory management in LUCIA (Lucas Visscher)
atomic supersymmetry now works with X2Cmmf (Ayaki Sunaga, Stefan Knecht, Trond Saue)
- New parallel coupled cluster implementation ExaCorr. Contributors: J.V. Pototschnig, A. Papadopoulos,
D. I. Lyakh, M. Repisky, L.Halbert, A.S.P. Gomes, H.J.Aa. Jensen, and L. Visscher. To be published.
- Improved start potential with far fewer functions and better accuracy. Contributor: Lucas Visscher.
S. Lehtola, L. Visscher and E. Engel. Efficient implementation of the superposition of atomic potentials initial guess for electronic structure calculations in Gaussian basis sets. J. Chem. Phys. 152, 144105 (2020); http://dx.doi.org/10.1063/5.0004046 doi: 10.1063/5.0004046
- Calculation of Rotational Molecular g-tensors. Activated with “.ROTG”. Contributor: I. Agustín Aucar.
I. A. Aucar, S. S. Gómez, C. G. Giribet and M. C. Ruiz de Azúa. Theoretical study of the relativistic molecular rotational g-tensor. J. Chem. Phys. 141, 194103 (2014); http://dx.doi.org/10.1063/1.4901422 doi: 10.1063/1.4901422
Improved reading of basis sets to be able to read all EMSL basis sets of type DALTON. Contributor: Hans Joergen Aa. Jensen.
“make latex” for making a latex version of the documenation in “build/latex”, used e.g. pdflatex to make a pdf version of the manual. Contributor: Hans Joergen Aa. Jensen
Efficient parallel diagonalization of quaternion Fock matrices using openMP. Contributor: Hans Joergen Aa. Jensen.
Utility (dirac_data.py) to export DIRAC data into a Python dictionary. Contributor: Lucas Visscher.
Gauge origin (.GAUGEO and .GO ANG) now can only be set under **HAMILTONIAN.
Phase and dipole origins (.PHASEO and .DIPORG, respectively) now can only be set under **HAMILTONIAN.
EOM-CC excitation energies collected for all symmetries and summarized at the end
One-electron operator ANGMOM’s origin was moved from the gauge-origin to the molecular center of mass.
The CODATA2018 set of physical constants is now used as default (taken from http://physics.nist.gov/constants). Before this version, CODATA1998 was default.
Non-relativistic and relativistic effective core-potential calculations (“.NONREL” and “.ECP”) now have nuclear point nucleus as default. For non-relativistic calculations this makes it easier to compare to other codes, and for RECP this is the intended behavior which by mistake was not activated.
- Calculation of Nuclear Spin-Rotation tensors. Activated with “.SPINRO”. Contributors: Trond Saue and I. Agustín Aucar.
I. A. Aucar, S. S. Gómez, M. C. Ruiz de Azúa, and C. G. Giribet. Theoretical study of the nuclear spin-molecular rotation coupling for relativistic electrons and non-relativistic nuclei. J. Chem. Phys. 136, 204119 (2012); http://dx.doi.org/10.1063/1.4721627 doi: 10.1063/1.4721627
- Interface with the Quantum Computing package OpenFermion. https://github.com/bsenjean/Openfermion-Dirac. Contributors: Bruno Senjean and Luuk Visscher.
T. E. O’Brien, B. Senjean, R. Sagastizabal, X. Bonet-Monroig, A. Dutkiewicz, F. Buda, L. DiCarlo and L. Visscher Calculating energy derivatives for quantum chemistry on a quantum computer. Accepted to npj: Quantum Information. https://arxiv.org/abs/1905.03742
- Expectation value of Nuclear Magnetic-Quadrupole-Moment interaction constant using “.OPERATOR” and “.NMQM” keyword in KRCI. Contributor: Malaya K. Nayak
T. Fleig, M. K. Nayak and M. G. Kozlov. TaN, a molecular system for probing P,T-violating hadron physics. Phys. Rev. A 93, 012505, 2016; http://dx.doi.org/10.1103/PhysRevA.93.012505 doi: 10.1103/PhysRevA.93.012505
Maximum number of basis functions increased from 30k to 40k.
Fixed use of point charges (do not try to obtain e.g. mass of a point charge).
Fixed spelling error in keyword .GAUNTSCALE (was .GAUTSCALE).
Fixed T1EQN1 function call in src/relccsd/ccgrad.F (thanks to Martin Siegert).
Fixed wrong memory allocation for spin-free or non-relativistic KRMCSCF.
Fixed runtime error for XPROPMAT in src/krmc/krmccan.F property modules.
Fixed runtime error for calling XPRIND in src/prp/pamprp.F property modules.
Basis set corrections:
Fixed Be basis set in cc-pV5Z.
Fixed Al basis set in aug-cc-pV(D+d)Z, aug-cc-pV(T+d)Z and aug-cc-pV(Q+d)Z.
Fixed Si basis set in aug-cc-pV(D+d)Z and aug-cc-pV(Q+d)Z.
Reduced static memory requirement from common blocks.
Restore local build of Sphinx documentation.
The visualization module has been extended to generate various radial distributions using the keyword .RADIAL. Contributor: Trond Saue.
Fix f correlating exponents for 3d in dyall.cv2z and dyall.ae2z basis sets. Contributor: Ken G. Dyall.
Fix f polarizing exponents for 3d in dyall.cv3z and dyall.ae3z basis sets. Contributor: Ken G. Dyall.
Bugfix in KRCI module (subroutine STRTYP_GAS_REL), reported by Shigeyoshi Yamamoto.
Bugfix in polarizable embedding (the printed electronic part of the energy was not corrected with PE contributions).
Fixed the formatting of the carbon potential in the file basis_ecp/ECPCE2SO.
New much more efficient quaternion diagonalization routine using openMP for performance. Contributor: Hans Joergen Aa. Jensen.
New much more efficient complex diagonalization routine using openMP for performance. Contributor: Hans Joergen Aa. Jensen.
Use MPI_REDUCE for adding Fock matrix contributions from MPI worker nodes on MPI master node. Contributor: Hans Joergen Aa. Jensen.
Save memory in direct Fock matrix constructions, now >13000 basis functions can be run in 32GB RAM. Contributor: Hans Joergen Aa. Jensen.
(With these performance improvements, a 4c HF test calculation with 13000 basis functions, 3300 postive energy orbitals (6600 in total) has been run on Titan at ORNL with 800 nodes, each with 8 openMP threads, at 1h59m-2h06m per HF iteration.)
Allow to compile with OpenMP support.
Update in CMake framework: pass compilers as cmake variables, not as environment variables.
- Kramers restricted Polarization Propagator in the ADC framework for electronic excitations, activated with “.POLPRP”.
M. Pernpointner. The relativistic polarization propagator for the calculation of electronic excitations in heavy systems. J. Chem. Phys. 140, 084108 (2014); http://dx.doi.org/10.1063/1.4865964
M. Pernpointner, L. Visscher and A. B. Trofimov. Four-component Polarization Propagator Calculations of Electron Excitations: Spectroscopic Implications of Spin-Orbit Coupling Effects. J. Chem. Theory Comput. (2017) submitted.
Polarizable embedding using PElib (https://gitlab.com/pe-software/pelib-public). Activated with “.PEQM”, additional options under *PEQM.
Reference: E. D. Hedegård, R. Bast, J. Kongsted, J. M. H. Olsen, and H. J. Aa. Jensen: Relativistic Polarizable Embedding., J. Chem. Theory Comput. 13, 2870-2880 (2017); http://doi.org/10.1021/acs.jctc.7b00162
New and numerically stable procedure for elimination/freezing of orbitals at SCF level. Author: T. Saue.
New “.MVOFAC” option in *KRMC input section. Author: H. J. Aa. Jensen.
New easier options for point charges in the .mol file: “LARGE POINTCHARGE” or “LARGE NOBASIS” (the two choices are equivalent).
Added the RPF-4Z and aug-RPF-4Z basis sets for f-elements to the already existing files with sets for s, p and d elements. Deleted the aug-RPF-3Z set as that was not an official set.
Write out effective Hamiltonian in Fock space coupled cluster to a file for post processing. Can be used with external code of Andrei Zaitsevskii (St. Petersburg).
Provided memory counter for RelCC calculations, what is suitable for memory consuming large scale Coupled Cluster calculations. See the web-documentation.
Updated basis_dalton/ with basis set updates in the Dalton distribution:
fix of errors in Ahlrichs-pVDZ (several diffuse exponents were a factor 10 too big)
fix of errors for 2. row atoms in aug-cc-pCV5Z
added many atoms to aug-cc-PVTZ_J
added many Frank Jensen “pc” type basis sets
added Turbomole “def2” type basis sets
Fixed errors for quaternion symmetries in 2-electron MO integrals used in CI calculations with GASCIP. It is now possible to do CI calculations with GASCIP for C1 symmetry (i.e. no symmetry).
Fixed the p exponents for Na in the dyall 4z basis sets to match the archive. The changes are small so should not significantly affect results.
Fixed compilation of XCFun on Mac OS X High Sierra.
Fixed error for parallel complex CI or MCSCF with GASCIP
.SKIPEP is now default for KR-MCSCF, new keyword .WITHEP to include e-p rotations
restored integral screening behavior from DIRAC15
RELCCSD expectation values aka first-order unrelaxed analytic energy derivative. For more information, see J. Chem. Phys. 145 (2016) 184107 as well as test/cc_gradient for an example.
Calculation of approximate dipole intensities together with the excitation energies in Fock space CC. See test/fscc_intensities for an example.
Improved start potential for SCF: sum of atomic LDA potentials, generated by GRASP.
Negative denominators (e.g. appearing in core ionized systems) accepted in RELCCSD
AOFOCK is now default if at least 25 MPI nodes (parallelizes better than SOFOCK). .AOFOCK documented.
- Error corrections and updates in isotope properties for the following atoms:
Br isotope 2: quadrupole moment .2620 -> .2615
Ag isotope 2: magnetic moment .130563 -> -.130691 (note sign change)
In isotope 2: quadrupole moment .790 -> .799
Nd magnetic moments of isotopes 4 and 5 were interchanged: -0.065 -> -1.065 and -1.065 -> -0.065
Gd: quadrupole meoments of isotopes 4 and 5 updated: 1.36 -> 1.35 and 1.30 -> 1.27
Ho isotope 1: quadrupole moment updated 3.49 -> 3.58
Lu isotope 2: quadrupole moment updtaed 4.92 -> 4.97
Hf isotope 1: mass was real*4, not real*8, thus 7 digits instead of 179.9465457D0 (i.e. approx 179.9465)
Ta isotope 1: quadrupole moment added 0.00 -> 3,17
Tl isotope 1: nuclear moment 1.63831461D0 -> 1.63821461D0 (typo, error 1.d-4)
Pb isotope 3: nuclear moment 0.582583D0 -> 0.592583D0 (typo, error 1.d-2)
Po isotope 1: nuclear moment added: 0.000 -> 0.688
Faster Dirac compilation due to implemented OBJECT libraries (requires CMake of version at least 2.8.10).
Run benchmark tests and only benchmark tests with “ctest -C benchmarks -L benchmark”.
Change any tab character in coordinate specifications in .mol file to a blank. (tab characters have been reported to cause errors in reading coordinates).
Added OPENBLAS in search list for blas and lapack libraries.
Made keyword .X2Cmmf case insensitive.
Cache compiler flags.
Bugfix: MCSCF and GASCIP was not correct for odd number of electrons.
Bugfix: now .CASSCF works again for KRMSCF input.
Bugfix: KRMSCF was not converging is some cases (problem reported by Miro Ilias for Os atom).
Bugfix: “all” keyword for “.GAS SHELLS” was not recognized (problem reported by Miro Ilias).
Bugfix for numerical gradient when automatic symmetry detection: do no move or rotate molecule for finding highest possible symmetry because then the numerical gradient will not be correct (problem reported by Miro Ilias).
Bugfix for spinfree X2Cmmf and linear symmetry.
Bugfix: restored “make install” target.
Bugfix: Test reladc_fano does not report failure if Stieltjes is deactivated.
FanoADC-Stieltjes: Calculation of decay widths of electronic decay processes.
Better convergence for open-shell calculations, using .OPENFAC 0.5 as default. Final orbital energies are recalculated with .OPENFAC 1.0, for IP interpretation.
Performance improvement for MP2 NO: only transform (G O|G O) integrals instead of (G G|G G).
Geometry optimization with xyz input.
Performance improvements for determinant generation in GASCIP.
Use Kramers conjugation on CI vectors (cuts time for CI in half for ESR doublets).
ANO-RCC basis: Fixed Carbon basis set (wrong contraction coefficients, see [MOLCAS ANO-RCC](http://www.molcas.org/ANO/).
ANO-RCC basis: Modified the 3 Th h-functions by replacing them with the 3 Ac h-functions to Th. (A mistake was made in the original work when the 3 Th h-functions were made, and this has never been redone. They are too diffuse, exponents (0.3140887600, 0.1256355100, 0.0502542000) for Th, Z=90, compared to (0.7947153600, 0.3149038200, 0.1259615200) for Ac, Z=89, and (0.8411791300, 0.3310795400, 0.1324318200) for Pa, Z=91. We have selected to just replace the 3 Th h-functions with those from the Ac basis set, because the Ac g-functions are quite close tot he Th g-functions, closer than Ac g-functions, and therefore differences in results compared to optimized Th h-functions should be minimal.) Thanks to Kirk Peterson for pointing out the Th problem on http://daltonforum.org.
Fixed reading of ANO-RCC and ANO-DK3 basis sets from the included basis set library.
Update PCMSolver module to its [1.0.4 version](https://github.com/PCMSolver/pcmsolver/releases).
Do not write Hartree-Fock in output if it is a DFT calculation, write DFT.
Fix a bunch of problems where the number of arguments of the call does not match the subroutine definition (submitted by Martin Siegert).
Fixed plotting of the electrostatic potential (electronic part was factor 2 too large).
Configuration framework uses [Autocmake](http://autocmake.org).
Added OpenBLAS to default blas search list.
Fixed misleading error message for .OPEN SHELL input.
Workaround for CMake warning about nonexistent run_pamadm target.
Intel compilers xHost flag automatic detection only upon request (-D ENABLE_XHOST_FLAG_DETECTION=ON).
Fixed minor issues to allow compilation using IBM XL Fortran compiler (thanks to Bob Walkup, IBM).
Update PCMSolver module to its [1.0.3 version](https://github.com/PCMSolver/pcmsolver/releases).
DIRAC14 release, see doc/release/release-statement.txt