An up-to-date overview of the CFOUR program system is given. After providing a brief outline of the evolution of the program since its inception in 1989, a comprehensive presentation is given of its well-known capabilities for high-level coupled-cluster theory and its application to molecular properties. Subsequent to this generally well-known background information, much of the remaining content focuses on lesser-known capabilities of CFOUR, most of which have become available to the public only recently or will become available in the near future. Each of these new features is illustrated by a representative example, with additional discussion targeted to educating users as to classes of applications that are now enabled by these capabilities. Finally, some speculation about future directions is given, and the mode of distribution and support for CFOUR are outlined in the appendix.
We present the self-consistent implementation of current-dependent (hybrid) meta-generalized gradient approximation (mGGA) density functionals using London atomic orbitals. A previously proposed generalized kinetic energy density is utilized to implement mGGAs in the framework of Kohn-Sham current density functional theory (KS-CDFT). A unique feature of the nonperturbative implementation of these functionals is the ability to seamlessly explore a wide range of magnetic fields up to 1 au (∼235 kT) in strength. CDFT functionals based on the TPSS and B98 forms are investigated, and their performance is assessed by comparison with accurate coupled-cluster singles, doubles, and perturbative triples (CCSD(T)) data. In the weak field regime, magnetic properties such as magnetizabilities and nuclear magnetic resonance shielding constants show modest but systematic improvements over generalized gradient approximations (GGA). However, in the strong field regime, the mGGA-based forms lead to a significantly improved description of the recently proposed perpendicular paramagnetic bonding mechanism, comparing well with CCSD(T) data. In contrast to functionals based on the vorticity, these forms are found to be numerically stable, and their accuracy at high field suggests that the extension of mGGAs to CDFT via the generalized kinetic energy density should provide a useful starting point for further development of CDFT approximations.
We present a formalism and an implementation for calculating spin-orbit couplings (SOCs) within the EOM-CCSD (equation-of-motion coupled-cluster with single and double substitutions) approach. The following variants of EOM-CCSD are considered: EOM-CCSD for excitation energies (EOM-EE-CCSD), EOM-CCSD with spin-flip (EOM-SF-CCSD), EOM-CCSD for ionization potentials (EOM-IP-CCSD) and electron attachment (EOM-EA-CCSD). We employ a perturbative approach in which the SOCs are computed as matrix elements of the respective part of the Breit-Pauli Hamiltonian using zeroth-order non-relativistic wave functions. We follow the expectation-value approach rather than the response-theory formulation for property calculations. Both the full two-electron treatment and the mean-field approximation (a partial account of the two-electron contributions) have been implemented and benchmarked using several small molecules containing elements up to the fourth row of the periodic table. The benchmark results show the excellent performance of the perturbative treatment and the mean-field approximation. When used with an appropriate basis set, the errors with respect to experiment are below 5% for the considered examples. The findings regarding basis-set requirements are in agreement with previous studies. The impact of different correlation treatment in zeroth-order wave functions is analyzed. Overall, the EOM-IP-CCSD, EOM-EA-CCSD, EOM-EE-CCSD, and EOM-SF-CCSD wave functions yield SOCs that agree well with each other (and with the experimental values when available). Using an EOM-CCSD approach that provides a more balanced description of the target states yields more accurate results.
An implementation of coupled-cluster (CC) theory to treat atoms and molecules in finite magnetic fields is presented. The main challenges for the implementation stem from the magnetic-field dependence in the Hamiltonian, or, more precisely, the appearance of the angular momentum operator, due to which the wave function becomes complex and which introduces a gauge-origin dependence. For this reason, an implementation of a complex CC code is required together with the use of gauge-including atomic orbitals to ensure gauge-origin independence. Results of coupled-cluster singles-doubles-perturbative-triples (CCSD(T)) calculations are presented for atoms and molecules with a focus on the dependence of correlation and binding energies on the magnetic field.
A program for the direct calculation of excitation energies of atoms and molecules in strong magnetic fields is presented. The implementation includes the equation-of-motion coupled-cluster singles-doubles (EOM-CCSD) method for electronically excited states as well as its spin-flip variant. Differences to regular EOM-CCSD implementations are due to the appearance of the canonical angular-momentum operator in the Hamiltonian causing the wave function to become complex. The gauge-origin problem is treated by the use of gauge-including atomic orbitals. Therefore, a modified Davidson method for diagonalizing complex non-Hermitian matrices is used. Excitation energies for selected atoms and molecules that are of importance in the astrochemical context are presented and their dependence on the magnetic field is discussed.
Trygve (2015) The importance of current contributions to shielding constants in density-functional theory. Physical Chemistry Chemical Physics, 17 (28 A note on versions:The version presented here may differ from the published version or from the version of record. If you wish to cite this item you are advised to consult the publisher's version. Please see the repository url above for details on accessing the published version and note that access may require a subscription.For more information, please contact eprints@nottingham.ac.uk J o u r n a l Na me The sources of error in the calculation of nuclear-magnetic-resonance shielding constants determined by density-functional theory are examined. Highly accurate Kohn-Sham wave functions are obtained from coupled-cluster electron density functions and used to define accurate-but current independent-density-functional shielding constants. These new reference values, in tandem with high-accuracy coupled-cluster shielding constants, provide a benchmark for the assessment of errors in common density-functional approximations. In particular the role of errors arising in the diamagnetic and paramagnetic terms is investigated, with particular emphasis on the role of current-dependence in the latter. For carbon and nitrogen the current correction is found to be, in some cases, larger than 10 ppm. This indicates that the absence of this correction in general purpose exchange-correlation functionals is one of the main sources of errors in shielding calculations using density functional theory. It is shown that the current correction improves the shielding performance of many popular approximate DFT functionals.
In this review, we discuss the current status of analytic derivative theory in relativistic quantum chemistry. A brief overview of the basic theory for the available relativistic quantum chemical methods as well as the state-of-the-art development of their analytic energy derivatives is given. Among the various relativistic quantum chemical methods, cost-effective approaches based on spin separation and/or on the matrix representation of two-component theory have been proven particularly promising for the accurate and efficient treatment of relativistic effects in electron correlation calculations. We highlight analytic derivative techniques for these cost-effective relativistic quantum chemical approaches including direct perturbation theory, the spin-free Dirac-Coulomb approach, and the exact two-component theory in its one-electron variant. An outlook is given on future developments of analytic energy derivative techniques for relativistic quantum chemical methods, again with an emphasis on cost-effective schemes.
Finite-field EOM-CCSDT: a highly accurate method for the theoretical prediction of excitation energies and electronic spectra in strong magnetic fields.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.