Kramers-restricted Polarization Propagator in the ADC framework for electronic excitations, activated with ”.POLPRP”.
References:
M. Pernpointner. The relativistic polarization propagator for the calculation of electronic excitations in heavy systems.J. Chem. Phys. 140, 084108 (2014)
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.
Reference: see T. Fleig and M. K. Nayak. Electron electric-dipole-moment interaction constant for HfF^{+} from relativistic correlated all-electron theory. Phys. Rev. A, 88:032514, 2013
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
New ”.MVOFAC” option in *KRMC input section for Modified Virtual Orbitals in MCSCF. Contributor: H. J. Aa. Jensen.
New and numerically stable procedure for elimination/freezing of orbitals at SCF level. Contributor: T. Saue.
Support for use of DIRAC in PyADF and QMFlows workflow engines. Contributors: Lucas Visscher, Andre Gomes and Christoph Jacob
New easier options for point charges in the .mol file: “LARGE POINTCHARGE” or “LARGE NOBASIS” (the two choices are equivalent), see here
Provided memory counter for RelCC calculations, suitable for memory consuming large scale Coupled Cluster calculations, see here for details. Contributor: Miroslav Iliaš
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).
Restart for RELCCSD. Contributor: Andre Gomes. See the keyword .RESTART and the section *CCRESTART
Performance Improvements
Restored integral screening. Contributor: Hans Jørgen Aagaard Jensen
POLPRP module + Davidson diagonalizer now parallel. Contributor: Markus Pernpointner
Corrections
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 error for parallel complex CI or MCSCF with GASCIP
Fixed compilation of XCFun on Mac OS X High Sierra.
Change of defaults
Change of final (open shell) orbital energies + SCF cycle modification. Contributors: Hans Jørgen Aagaard Jensen and Trond Saue
.SKIPEP is now default for KR-MCSCF, new keyword .WITHEP to include e-p rotations
Basis set news
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.
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.
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
New features in DIRAC16
RELCCSD expectation values. For more information, see J. Chem. Phys. 145 (2016) 184107 as well as test/cc_gradient for an example.
Improved start potential for SCF: sum of atomic LDA potentials, generated by GRASP.
Change of defaults
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). And .AOFOCK documented.
Corrections
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 moments 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)
Fixed Carbon basis set (wrong contraction coefficients, see [MOLCAS ANO-RCC](http://www.molcas.org/ANO/).
Modified the 3 Th h-functions by replacing them with the 3 Ac h-functions to Th.
Fixed reading of ANO-RCC and ANO-DK3 basis sets from the included basis set library.
New defaults
For open-shell SCF calculations, .OPENFAC = 0.5 by default, as this seems to improve convergence. Final orbital energies are recalculated with .OPENFAC 1.0, for IP interpretation.
Intrinsic Atomic Orbitals (IAOs), as formulated by Gerald Knizia, have been implemented to eliminate the polarization contribution in projection analysis.
The Polarizable Continuum Model (PCM) is available for the inclusion of solvent effects. For more details, see this paper
As a byproduct of the PCM implementation, molecular electrostatic potential (MEP) maps are available for 4-component electronic-structure calculations, see this paper
+Q corrections (size-consistency corrections) for KR-CI calculations
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.