where expansion coefficientsli
; aij andbijk correspond to dipole
moment, dipole polarizability and first hyperpolarizability, respectively. Eq.conforms to the Einstein’s summation convention.
The diagonal components of these tensorial quantities were calculated based on the formulas derived using the finite-difference
method: