DFT magnetizabilities. Contributors: M. Olejniczak and Trond Saue.
Limitations: Magnetizabilities and NMR shieldings calculated at the DFT level are so far restricted to C1 symmetry, but we expect to fix this soon.
Enhancements to the frozen density embedding (FDE) functionality
FDE contributions to magnetic properties (NMR shieldings, indirect spin-spin coupling constants, magnetizabilities), see ".FDE" and "*FDE" entries of the manual for details. Contributors: M. Olejniczak, R. Bast, A. S. P. Gomes
FDE interaction energies with CCSD, MP2 and mean-field densities. Contributors: M. Olejniczak, A. Shee, R. Bast, A. S. P. Gomes
Equation of motion coupled cluster
Energies for electronic excitations (EE), electron attachment (EA) and electron detachment (IP), see ".EOMCC", "*EOMCC" and "*CCDIAG" entries of the manual for details. Contributiors: A. Shee, T. Saue, L. Visscher, A. S. P. Gomes
A. Shee, T. Saue, L. Visscher, A. S. P. Gomes Equation-of-motion coupled-cluster theory based on the 4-component Dirac-Coulomb(-Gaunt) Hamiltonian. Energies for single electron detachment, attachment, and electronically excited states.J. Chem. Phys. 149 (2018) 174113
Polarized embedding can be done with xyz-files. Contributor: Trond Saue
Improved quaternion diagonalization Contributor: H. J. Aa. Jensen
Improvements in the visualization module (**VISUAL) Contributors: M. Olejniczak and T. Saue.
the possibility to calculate the NMR shielding tensor in a selected point in space, see ".NICS"
the possibility to visualize various densities on an imported 3D grid, see ".3D_IMP"
the possibility to calculate magnetic properties densities using the imported magnetically-induced current density, see ".READJB"
the possibility to scale densities by Cartesian products xiyjzk, see ".CARPOW"
the possibility to generate and visualize radial distributions, see ".RADIAL"
Change of defaults
New convergence criterium for CC amplitude equation: The convergence criterium for the amplitude equations that determine the CC energy has been revised and made consistent with the criterium used in the lambda equations used for molecular properties. In both cases we now take the norm of the differences between amplitudes of subsequent iterations. In practice this typically means the program will use a few iterations less. For normal calculations this is of no consequence as the default is still to converge very tightly, but if extremely high precision is required one may need to check the achieved convergence.
Change in the reorthonormalization terms in the calculation of magnetic properties with London atomic orbitals: the reorthonormalization and response contributions involve the same orbital pairs, for instance if all rotations between occupied and virtual orbitals are present in response equations, the reorthonormalization terms are also constructed from all orbital blocks; the keywords .DOEPRN and .NOEPRN under *NMR are depreciated;
New features in DIRAC17
Kramers-restricted Polarization Propagator in the ADC framework for electronic excitations, activated with ”.POLPRP”.
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. 14, 1510 (2017).
Reference: M. Denis, M. Nørby, H. J. Aa . Jensen, A. S. P. Gomes, M. K. Nayak, S. Knecht, and T. Fleig. Theoretical study on ThF+, a prospective system in search of time-reversal violation.New J. Phys., 17:043005, 2015
Functional derivatives using automatic differentiation (XCFun)
New visualization options
RKBIMP: MO-coefficients generated using restricted kinetic balance (RKB) can be extended by their unrestricted kinetic balance (UKB) complement, thus providing magnetic balance for response calculations involving external magnetic fields
New and improved 2c Hamiltonian schemes
New build system and infrastructure
New compilation scheme: configure replaced by CMake mechanism
New pam script (python)
Alternative launcher: wrapper.py (python)
New testing framework based on python (runscript)
Many static allocation calls replaced by dynamic allocation; in practice this means that you may need less WORK array memory and/or more space for dynamic allocation compared to DIRAC10.