1. Introduction
A correct interpretation of electronic spectra for transition metal complexes (d-d transitions), magnetometric data (magnetic susceptibility and magnetization), and spectra of electron spin resonance requires appropriate theoretical support. A traditional approach is represented by the crystal field theory, which is well elaborated for octahedral complexes (O
h symmetry), even with tetragonal (trigonal) distortion (D
4h, D
3d) [
1,
2,
3,
4,
5]. Analogously, tetrahedral patterns (T
d) and their distortion daughters to prolate and/or oblate bispheoids (D
2d) are also known. However, one is rather helpless when dealing with pentacoordination in its limiting cases represented by trigonal bipyramids (D
3h) and tetragonal pyramids (C
4v) and especially for intermediate geometries on the Berry rotation path (C
2v).
A target of the present work is to elucidate a comprehensive view of the crystal field terms and crystal field multiplets in the case of pentacoordinate d5 to d9 complexes. Whereas multielectron crystal field terms are labelled according to Mulliken notation (A, B, E, T), the involvement of the spin–orbit interaction requires a passage from common symmetry point groups to double groups; therefore, crystal field multiplets are labelled according to Bethe notation (Γ1 to Γ8).
Geometries belonging to point groups O
h, T
d, D
4h, or D
2d are mostly omitted hereafter; numerical computer-assisted treatment is necessary when ligands occupy arbitrary positions. This approach is slightly more complicated, involving algebra of complex numbers due to the occurrence of complex spherical harmonic functions fixing the ligand positions. The treatment used below is termed the Generalized Crystal Field Theory, as outlined elsewhere [
6].
The eigenvalues of the model Hamiltonian refer to the energy levels at a given approximation. The eigenvectors bear all information about the symmetry of the wave function; therefore, they can be utilized to assign irreducible representations (IRs) either of the crystal field terms or the crystal field multiplets . The irreducible representation within point group G is Γ, γ (its component when IR is degenerate), and a (the branching (repetition) number). The same holds true for double group G′.
2. Results
The generalized crystal field theory (GCFT) applied below is fully described elsewhere, along with the closed formulae for the matrix elements of the involved operators in the basis set of the electronic atomic terms
, where the apparatus of the irreducible tensor operators has been utilized [
6,
7]. (Here, the seniority number (
ν) for the terms can distinguish between terms possessing the same set
; the quantum numbers adopt their usual meaning [
8]). These matrix elements refer to five operators:
- (a)
Interelectronic repulsion parametrized by the Racah parameters BM and CM;
- (b)
Crystal field potential depending upon crystal field poles (strengths) Fk(RL) of the k-th order (k = 4, 2) for each ligand (L) and its position;
- (c)
Spin–orbit interaction depending upon the spin–orbit coupling constant (ξM);
- (d)
The orbital Zeeman term , which eventually involves the orbital reduction factors; and
- (e)
The spin Zeeman term , which contains the spin-only magnetogyric (ge) factor.
The position of ligands (L) is arbitrary and fixed by the polar coordinates {}. The model Hamiltonian involves three important cases:
- 1.
Diagonalization of (a) + (b) yields the energies of crystal field terms , which span the IRs of point group G;
- 2.
Diagonalization of (a) + (b) + (c) produces energies of the crystal field multiplets in the zero magnetic field , which span the IRs of double group G′;
- 3.
Diagonalization of (a) + (b) + (c) + (d) + (e) gives the magnetic energy levels in the applied magnetic field.
The energy levels of crystal field multiplets for the half-integral spin (S = 1/2, 3/2, 5/2) appear as Kramers doublets and remain doubly degenerate in the absence of a magnetic field. This is the case of high-spin Fe(III), Mn(II), Co(II), and Cu(II) complexes.
Traditional crystal field theory operates with a set of collective parameters, such as 10Dq = Δ, Ds, Dt, etc., and is useful for cases certain symmetry, such as Oh, Td, and D4h, all of which are derived from the crystal field poles (F4(L) and eventually F2(L)), e.g.,
For Oh/D4h: 10Dq = (10/6)F4(xy);
Dt = (2/21)[F4(xy) − F4(z)];
Ds = (2/7) [F2(xy) − F2(z)]; and
For Td: 10Dq = (20/27)F4(xy).
The crystal field poles originate in the partitioning of the matrix elements of the crystal field potential into radial (
R) and angular (
A) parts in the polar coordinates
.
The integration of the angular part yields some values (manually calculated for some cases, such as O
h symmetry). This part contains the spherical harmonic functions
for the positions of ligand
K, and in general, it is a complex number. The radial part contains the metal–ligand distance (
RK) and defines the crystal field poles:
(
k = 0, 2, 4), where the integration runs over the electronic coordinates. The matrix elements of the crystal field operators can be expressed as:
where for the reduced matrix elements
,
and the 3j symbols, closed formulae exist [
6,
7] and can be evaluated with a desktop computer.
In practice, the crystal field poles are not subject to evaluation; they are taken as parameters of the theory and depend on the quality of the ligand (halide, amine, phosphine, cyanide, carbonyl, etc.), as well as the quality and oxidation state of the central atom. For practical applications, the spectroscopic series is used according to the Δ-value [
4]. The values of Δ can be deduced from the transitions observed in the electronic d-d spectra. Moreover, the Δ value can be estimated based on the empirically determined increments
fL for the ligands and
gM for the central atoms
However, the same ligand can produce different crystal field strengths depending on the actual metal–ligand distance (cf. Equation (2)). For instance, the -NCS– group can be attached at distance R(Ni–N) = 2.2 or 2.0 Å. In the second case, it produces a much stronger crystal field.
For the hexacoordinate complexes, value of F4 = 5000 cm−1 refers to Δ(Oh) = 8300 cm−1, which is a weak crystal field (appropriate for the halido ligand). Then, F4 = 15,000 cm−1 is equivalent to Δ(Oh) = 25,000 cm−1, which refers to the strong crystal field (appropriate for cyanido or carbonyl ligands). For tetrahedral complexes, F4 = 5000 (15,000) cm−1 refers to Δ(Td) = 3700 (11,100) cm−1.
2.1. Crystal Field Terms
Figure 1 displays the relative energies of the crystal field terms (not to scale) for individual d
n configurations. These result from the GCFT calculations using the weak crystal field characterized by the crystal field poles
F4(L) = 5000 cm
−1 for each ligand. For the tetragonal pyramid (C
4v), the angle L
a-M-L
e = 104 deg was maintained. The passage from the fully rotation group R
3 of a free atom to point group D
3h or C
4v is shown as the splitting of the atomic terms by the crystal field. The literature outlines the branching rules for such a reduction process [
9].
The character tables for the point groups usually assign the dipole moment components to the IRs; these are useful in determining the selection rules for the excitation energies. For instance, within group D3h, the direct product of IRs is , meaning that the z component of the dipole moment is active in transition , yielding the non-zero transition moment (orbitally allowed transition). On the contrary, and thus transition are forbidden.
In addition to the energy levels,
Figure 1 also shows the allowed/forbidden polarized electronic dipole transitions, which are displayed as solid/dashed arrows. These data can be compared with the observations of the electronic d-d spectra [
10].
2.2. Crystal Field Multiplets
The crystal field terms represent a starting point for the further precision of the energy levels: upon introduction of the spin–orbit interaction, the crystal field terms are further spit into a set of crystal field multiplets (energy levels in the zero magnetic field) [
11,
12]. The basis set and the resulting multiplets contain 256, 210, 120, 45 and 10 members for electron configurations d
5 through d
9, respectively. In this case, the energy levels are labelled using the Bethe notation for the IRs within double group G’. These symbols involve Γ
1 through Γ
8, and their degeneracy is shown parentheses, e.g., Γ
4(2). (The IR tables for the double groups are useful for practical reasons).
The spin and the orbital parts of the wave function are assessed independently. For instance, in D3h the level, 6A1′(6) transforms its spin according to {2Γ4 + (Γ5 + Γ6)}. The orbital part matches A1 = Γ1. Finally, the spin–orbit wavefunction transforms according to the direct product {2Γ4 + (Γ5 + Γ6)}Γ1, and the result is {2Γ4 + (Γ5 + Γ6)}. In this special case, the levels (Γ5 + Γ6) form a complex conjugate pair that can be abbreviated as Γ5,6(2) or simply Γ5(2). To this end, upon passage from the D3h to double group D’3, the crystal field term 6A1′(6) is split into a set of {2Γ4(2) + Γ5,6(2)} multiplets. However, this part of the theory says nothing about the relative energies of the final three Kramers doublets; these result from numerical calculations by GCFT.
The principal result of the CGTF calculations with spin–orbit coupling in the complete space spanned by d
n configurations is the spectrum of the crystal field multiplets. The lowest zero-field energy gaps are abbreviated as
δ1,
δ2, …, provided that the energy of the ground multiplet (
δ0) is set to zero (
Table 1). For the non-degenerate ground state (A or B type), the lowest multiplet gaps relate to the axial zero-field splitting parameter (
D). For d
5-Fe(II) (and, analogously, d
5-Mn(III)), the sequence of the spin–orbit multiplet does not strictly follow
D and 4
D (there is a small difference (
δa) around 4
D). For Cu(II), the ground electronic term is not split by the spin–orbit interaction; however, the spin–orbit multiplets are slightly influenced by the spin–orbit coupling. The concept of the
D parameter is strictly related to the spin–Hamiltonian theory.
The effect of the spin-orbit interaction leading to the passage from the crystal-field terms to the crystal-field multiplets is depicted in
Figure 2.
2.3. Zero-Field Splitting
The concept of the spin Hamiltonian is a popular and very useful tool for interpretation of the spectra of electron paramagnetic resonance, as well as for analysis of DC magnetometric data. The key formulae of the spin Hamiltonian are based on consideration of only the spin kets
of non-degenerate ground term A or B. The second-order perturbation theory offers the Λ tensor in the following form:
where
K runs over all excited electronic terms, and the magnetic tensors are expressed as follows:
the
κ tensor (reduced, temperature−independent paramagnetic susceptibility tensor):
the
g tensor (magnetogyric ratio tensor):
the
D tensor (spin–spin interaction tensor):
This approximation fails in the case of orbital (pseudo) degeneracy. The matrix elements of the angular momentum can be assessed by exploiting the symmetry of the ground and excited crystal field terms; the matrix element is non-zero only if the direct product () contains the irreducible representation of at least one component of the angular momentum. For instance, within group D3h, and the common character tables indicate that the result contains the irreducible representation of Lx and Ly.
The spin Hamiltonian parameters calculated via the GCFT are listed in
Table 2. For Fe(III) and Mn(II), the ground electronic term
6A does not allow transitions to excited terms with different spin multiplicities. Therefore,
D = 0,
gi =
ge in this approximation. In this case, the spin Hamiltonian formalism is insufficient, so
6A
1 +
4T
1 terms must be considered for the O
h symmetry [
13].
The spin Hamiltonian is often presented in the following form:
where the
D tensor is considered diagonal and traceless, yielding only two independent parameters: the axial zero-field splitting parameter (
D) and the rhombic zero-field splitting parameter (
E). This form is widely used for analysis of magnetometric and EPR data. According to convention, the rhombic part is minor: |
D| > 3
E > 0.
D serves as a measure of zero-field splitting. This energy gap can also be measured also by FAR-infrared spectroscopy (FIRMS and FDMRS techniques), inelastic neutron scattering, calorimetry, etc. [
14].
As mentioned above, the case of d
6-Fe(II) or d
6-Mn(III) is specific, as for C
4v geometry, the sequence of the spin–orbit multiplets differs depending on the exact multiplet splitting {0,
δ1(2),
δ22′(1 + 1) and the spin Hamiltonian formalism {0,
D(2), 4
D(2)}; the number in parentheses corresponds to the multiplicity. The ground crystal field term is
5B
2; the orbital and spin parts transform as B
2 → Γ
4,
S = 2 → Γ
1 + Γ
3 + Γ
4 + Γ
5, and their direct product is Γ
4(Γ
1 + Γ
3 + Γ
4 + Γ
5) = Γ
1(1) + Γ
2(1) + Γ
4(1) + Γ
5(2). Only Γ
5(2) is doubly degenerate, whereas the remaining multiplets are nondegenerate: Γ
1(1), Γ
2(1), and Γ
4(1). The GCFT calculations for d
6-Fe(II) in the complete basis set of 210 kets obtains Γ
4 as the ground multiplet and the multiplet splitting
E(Γ
5) −
E(Γ
4) =
δ1 = 0.31 cm
−1;
E(Γ
1) −
E(Γ
4) =
δ2 = 1.64 cm
−1;
E(Γ
2) −
E(Γ
4) =
δ2′ = 1.87 cm
−1. This feature is reflected in the spectrum of electron paramagnetic resonance. Details about the symmetry rules are listed in
Supplementary Information.
Table 3 shows a comparison of the d
n configurations from the viewpoint of the spin Hamiltonian formalism. This table is also enriched by data for d
1 to d
4 configurations, as well as data for the intermediate geometry with C
2v symmetry and
τ5 = 0.47.
Table 4 analogously summarizes data for the hexacoordinate complexes.
2.4. DC Magnetic Functions
The magnetic energy levels
result from the diagonalization of the interaction matrix ((a) + (b) + (c) + (d) + (e)), which includes interelectronic repulsion, crystal field potential, spin–orbit coupling, and orbital and Zeeman terms in the applied magnetic field. Statistical thermodynamics offers formulae for magnetization and magnetic susceptibility when the partition function is evaluated for three reference fields:
Bm =
B0 −
δ,
B0,
B0 +
δ (allowing numerical derivatives):
Hence, the molar magnetization is:
The molar magnetic susceptibility is expressed as:
where the physical constants adopt their usual meaning. The index
a refers either to the Cartesian coordinates {
x,
y,
z} or to the grid point over a sphere along which the magnetic field is aligned, which is used to obtain the powder-sample average. Therefore, the magnetic susceptibility and magnetization are functions of discrete parameters (atomic parameters
BM,
CM, and
ξM; ligand positions
θL and
φL; crystal field poles
F4(L); and eventually
F2(L)), as well as the continuous parameters, such as the reference field (
Bm) and temperature (
T).
The modelling of the magnetization and susceptibility for pentacoordinate d
n systems is presented in
Figure 3 and
Figure 4. A counterpart of these graphs for the tetragonally distorted octahedral systems can be found elsewhere [
16]. In the case of zero-field splitting with an orbitally non-degenerate ground term, the effective magnetic moment in the high-temperature limit of 300 K remains almost linear with zero slope; at low temperature, it is reduced. This is the case of d
5-D
3h, d
5-C
4v, d
6-C
4v, d
7-D
3h, and d
8-C
4v. For d
9-D
3h and d
9-C
4v, zero-field splitting is absent, so these systems follow the Curie law. The magnetization saturates to the value of
M1 =
Mmol/(
NAμB) =
gavS when the zero-field splitting is small. This is the case of d
5-D
3h, d
5-C
4v, d
6-C
4v, d
9-D
3h, and d
9-C
4v; exceptions are d
7-D
3h and d
8-C
4v, with large zero-field splitting
D parameters.
Systems with E-type orbitally doubly degenerate ground terms, such as d6-D3h, d7-C4v, and d8-D3h behave differently. The effective magnetic moment is enlarged, and it passes through a round maximum. The magnetization is also suppressed and does not reach saturation until B = 10 T.
A positive slope of the effective magnetic moment reflects the effect of the low-lying excited electronic terms mixed considerably with the ground term via the spin–orbit interaction. This results in temperature-independent paramagnetism,
χTIP > 0. This term, along with the underlying diamagnetism (
χdia < 0), need be subtracted from the measured temperature dependence of the magnetic susceptibility. With respect to the underlying diamagnetism, a method of additive Pascal constants is useful and frequently utilized. However, for temperature-independent paramagnetism, the amount of information is considerably limited [
6,
7].
2.5. AC Magnetic Susceptibility
In the oscillating magnetic fields (usually with a low amplitude of
BAC = 0.3 mT and a frequency range of
f = 10
−2–10
5 Hz), the measured magnetic moment of the specimen has two components: in-phase and out-of-phase. This is easily transformed into two components of AC susceptibility:
χ′ (dispersion) and
χ″ (absorption). The absorption component is a measure of the resistivity of the sample used to alter its magnetization; it provides information about the relaxation time, which is a function of temperature, frequency (
f), and the external applied field (
BDC). The relaxation time can be inferred from the position of the maximum at the out-of-phase susceptibility (
f″
max) with the following formula:
τ = 1/(2π
f″
max). It has been reported that the sample can exhibit two or more relaxation channels and that their absorption curves can overlap or merge to form a shoulder. The whole AC susceptibility can be fitted by exploiting the generalized Debye equation [
17,
18]:
where
K is the number of relaxation channels,
χS is the common adiabatic susceptibility (high-frequency limit),
χk is the thermal susceptibilities,
αk is the distribution parameters,
τk is the relaxation times, and the circular frequency is
ω = 2π
f. This complex equation can be decomposed into a real and imaginary part.
The slow magnetic relaxation includes several mechanisms that can be collected to a single equation for the reciprocal relaxation time:
The first term describes the thermally activated Orbach process, which is associated with the height of the barrier to spin reversal (
Ueff); the second is the Raman term, with the temperature exponent typically
n = 5–9; next is the phonon bottleneck term, with
l ~ 2; the fourth term describes the direct relaxation process, with
m = 2–4; the last term refers to the quantum tunnelling of magnetization throughout the barrier to spin reversal. The reciprocating thermal behavior was recently registered with a term analogous to the phonon bottleneck but a negative temperature exponent (
l ~ −1) [
18].
The effectiveness of the slow magnetic relaxation is, as a rule, evaluated by the value of
Ueff (when the Orbach process applies). It is assumed that it is related to the axial zero-field splitting parameter (
D), which must be negative, and the molecular spin (
S) [
19]:
which holds true for Kramers systems with half-integral spin (e.g.,
S = 3/2 for Co
II); for non-Kramers systems with an integer spin, the factor ¼ is dropped (e.g.,
S = 1, for Ni
II). It is common practice for the
Ueff and the pre-exponential factor (
τ0) to be subtracted using the Arrhenius-like plot ln(
τ) vs. 1/
T (
Figure 5-left): a few high-temperature points are fitted by the straight line, tangential of which refers to
Ueff. However, “high-temperature points” refer to the highest temperature among the data considered in our analysis, so there still could be points yielding a higher tangential and thus
Ueff. A preferred approach involves plotting ln(
τ) vs. ln(
T), where the temperature exponent recovering the high-temperature data refers to the slope (
Figure 5, right). When the temperature coefficient is
n > 9, instead of the Raman process the Orbach process is applied.
For high-spin Co(II) complexes with
S = 3/2, eqn (15) implies a relationship of
U = 2|
D|. A collection of experimental data for a series of tetracoordinate Co
II complexes is shown in
Figure 6 based on the analysis of higher-temperature, high-frequency relaxation data in terms of the Orbach process. Evidently, a correlation of
U vs. 2|
D| fails.
D is a field-independent quantity, whereas the extracted value of
U depends upon the applied magnetic field. A positive value of
D contradicts the
D-
U paradigm; however, SIMs behavior can occur (the Raman mechanism is likely the leading term). With increased barrier to spin reversal (
U), the extrapolated relaxation time (
τ0) is shortened, irrespective of the sign of the
D parameter. A violation of the
D-
U paradigm has been discussed elsewhere with consideration of anharmonicity contributions [
21].
3. Discussion
The GCFT approach enables fast and “continuous” mapping of the energy levels, such as electronic terms and the spin–orbit multiplets: one is free to changing the ligand positions {θL, φL} from regular coordination polyhedra to distorted polyhedra and to alter the crystal field poles F4(L) and, eventually, F2(L). On the contrary, the modern ab initio calculations provide high-quality data on energy levels but only for the unique geometry of the complex under investigation. Therefore, it is interesting to utilize and compare both approaches.
Ab initio calculations have been performed using ORCA software [
22] with respect to the experimental geometry of the complexes resulting from X-ray structural analysis (the corresponding cif files are deposited in the Cambridge Crystallographic Data Centre). The relativistic effects were included in the calculations with a second-order Douglas–Kroll–Hess (DKH) procedure. An extended basis set TZVP of Gaussian functions was used, e.g., BS1 = [17s11p7d1f] and BS2 = [17s12p7d2f1g] for Ni(II). The calculations were based on state-average complete active-space self-consistent field (SA-CASSCF) wave functions. The active space of the CASSCF calculations comprised eight electrons in five metal-based d-orbitals. The state-averaged approach was used, whereby all 10 triplet and 15 singlet states were equally weighted. The spin–orbit effects were included according to quasi-degenerate perturbation theory, whereby the spin–orbit coupling operator (SOMF) was approximated according to the Breit–Pauli form. The electronic terms were evaluated at the CASSCF + NEVPT2 level, and the multiplets by considering the spin–orbit interaction (
Table 5). Effective Hamiltonian was used to evaluate the spin Hamiltonian parameters.
The ab initio calculations refer to in silico state, i.e., intermolecular interactions and other solid-state effects are ignored. This is not the case for experimental magnetometric or spectroscopic data, which could be influenced by the environment. The ab initio data, in general, are consistent with those obtained by experimental techniques.
When comparing the CGTF calculations with ab initio calculations, calculated transition energies can be assessed. With a proper set of crystal field poles, the CGTF can reproduce first allowed transitions; however, the electronic spectrum, has a smaller width with respect to ab initio data.
An extended set of similar pentacoordinate Ni(II) complexes based on the fixed skeleton of a pentadentate Schiff base (
Figure 7) was investigated by magnetometry and
ab initio calculations with respect to the experimental geometry; these are listed in
Table 6.
The experimentally reported and calculated
D values cover a broad interval of positive and negative values over a wide range of the
τ5 parameters. These were used to plot
D vs.
τ5, which can be termed the
second magnetostructural D-correlation for Ni(II) complexes (MSDC). (The first magnetostructural
D correlation for hexacoordinate Ni(II) complexes is outlined elsewhere [
30].) The MSDC can be approximated by a straight line (
Figure 8) when the
τ5 parameter guarantees that the ground electronic term is not orbitally quasi-degenerate (the energy gap Δ > 2000 cm
−1). In the opposite case, the calculated
D values tend to diverge. The value of the D parameter switches between positive and negative values at
τ5 ~ 0.2–0.3. Furthermore, the E parameter plays a role that has not be considered so far.
4. Conclusions
Experimental data on magnetic susceptibility, magnetization, and electron paramagnetic resonance require an appropriate model in order be analyzed correctly. For some shapes of coordination polyhedra, such as octahedron O
h, tetragonal bipyramid D
4h, trigonal antiprism D
3d, tetrahedron T
d, and bispehoid D
2d, the crystal field theory offers such a support, and the spin Hamiltonian formalism defines relationships for the set of magnetic parameters (
D,
E,
gx,
gy,
gz,
χTIP). A dearth in the literature with respect to pentacoordinate systems, such as the trigonal bipyramid D
3h and tetragonal pyramid C
4v symmetry, is filled by this publication. The working tool is the generalized crystal field theory in the form of its fully numerical, computer-assisted tool [
31]. The advantage of this approach is that the positions of the ligands can be arbitrary, making it applicable to any geometry of the chromophore and any ligands. Only the set of Racah parameters of the interelectronic repulsion (
BM and
CM), the spin–orbit coupling constant (
ξM), polar angles (or Cartesian coordinates) of each ligand {
θL,
φL}, the crystal field poles
F4(L) and, eventually,
F2(L) are required. This method enables evaluation of the energies of the multielectron crystal field terms, spin–orbit crystal field multiplets, and the magnetic energy levels at the applied magnetic field. Then, the magnetic susceptibility and magnetization can be evaluated as functions of the temperature field via derivatives of the partition function. The eigenvectors provide complete information about the symmetry and can be used to automatically label terms/multiplets.