Density functional restricted-unrestricted approach for nonlinear properties: application to electron paramagnetic resonance parameters of square planar copper complexes

The density functional restricted-unrestricted approach for treatments of spin polarization effects in molecular properties using spin restricted Kohn-Sham theory has been extended from linear to nonlinear properties. It is shown that the spin polarization contribution to a nonlinear property has the form of a quadratic response function that includes the zero-order Kohn-Sham operator, in analogy to the lower order case where the spin polarization correction to an expectation value has the form of a linear response function. The developed approach is used to formulate new schemes for computation of electronic g-tensors and hyperfine coupling constants, which include spin polarization effects within the framework of spin restricted Kohn-Sham theory. The proposed computational schemes are in the present work employed to study the spin polarization effects on electron paramagnetic resonance spin Hamiltonian parameters of square planar copper complexes. The obtained results indicate that spin polarization gives rise to sizable contributions to the hyperfine coupling tensor of copper in all investigated complexes, while the electronic g-tensors of these complexes are only marginally affected by spin polarization and other factors, such as choice of exchange-correlation functional or molecular structures, will have more pronounced impact on the accuracy of the results.