Relationship between Structure and Zero-Field Splitting of Octahedral Nickel(II) Complexes with a Low-Symmetric Tetradentate Ligand

: Octahedral nickel(II) complexes are among the simplest systems that exhibit zero-field splitting by having two unpaired electrons. For the purpose of clarifying the relationship between structure and zero-field splitting in a low-symmetric system, distorted octahedral nickel(II) complexes were prepared with a tetradentate ligand, 2-[bis(2-methoxyethyl)aminomethyl]-4-nitrophenolate(1 − ) [(onp) − ]. The complex [Ni(onp)(dmso)(H 2 O)][BPh 4 ] · 2dmso ( 1 ) (dmso = dimethyl sulfoxide) was characterized as a bulk sample by IR, elemental analysis, mass spectrometry, electronic spectra, and magnetic properties. The powder electronic spectral data were analyzed based on the angular overlap model to conclude that the spectra were typical of D 4 -symmetric octahedral coordination geometry with a weak axial ligand field. Simultaneous analysis of the temperature-dependent susceptibility and field-dependent magnetization data yielded the positive axial zero-field splitting parameter D ( H = g u β S u H u + D [ S z 2 − S ( S + 1)/3]), which was consistent with the weak axial ligand field. Single-crystal X-ray analysis revealed the crystal structures of [Ni(onp)(dmso)(H 2 O)][BPh 4 ] · dmso ( 2 ) and [Ni(onp)(dmf) 2 ][BPh 4 ] ( 3 ) (dmf = N , N -dimethylformamide). The density functional theory (DFT) computations based on the crystal structures indicated the D 4 -symmetric octahedral coordination geometries with weak axial ligand fields. This study also showed the importance of considering g -anisotropy in magnetic analysis, even if g -anisotropy is small.


Introduction
Octahedral nickel(II) complexes are among the simplest systems that exhibit zero-field splitting by having two unpaired electrons.(Please note that octahedral complexes do not necessarily have ideal octahedral symmetry [1], and zero-field splitting is affected by the slight symmetry reduction in the ligand field [2]).However, since the nickel(II) complex is an S = 1 non-Kramers system with an even number of unpaired electrons, highfrequency and high-field electron paramagnetic resonance (HFEPR) is required to observe the electron-paramagnetic-resonance (epr) signals to precisely determine the zero-field splitting parameters, D and E [3].Here, the axial zero-field-splitting parameter D and the rhombic zero-field-splitting parameter E are defined as the phenomenological Hamiltonian H = g u βS u H u + D[S z 2 − S(S + 1)/3] + E(S x 2 − S y 2 ) [4], which can also be derived from the formal spin-Hamiltonian [2].Then, in the S = 1 system, the three originally degenerate microstates are split into −2D/3, D/3 − E, and D/3 + E energy levels without an external magnetic field.On the other hand, attempts to determine the zero-field-splitting parameters from magnetic data analysis have been successful [2,5], although not as accurate as HFEPR.The relationship between the coordination geometry around the nickel(II) ion and the sign of the axial zero-field splitting parameter D has been studied [2], and the result indicates that D is negative when the axial ligand field is strong (e z > e x,y ) and positive when the axial ligand field is weak (e z < e x,y ), where e u represents the metal-ligand interaction in the u direction.Since the sign of the D parameter is related to magnetic anisotropy and the conditions for the single molecule magnets [2,[6][7][8][9], even if there are exceptions [10], it is important to clarify the relationship between the sign of D and structure and to control the sign of D.
To determine the zero-field splitting parameters from the magnetic analysis, simultaneous fitting of the temperature dependence of magnetic susceptibility and the field dependence of magnetization is preferable [2].For this purpose, it is, of course, essential to simulate magnetic susceptibility and magnetization.First, in the days when the cryomagnetic susceptibility was measured only down to about 5 K, it was not necessary to consider the effect of saturation by the applied magnetic field.However, now that it is common to measure magnetic susceptibility down to 2 K, a magnetic-field-dependent magnetic susceptibility equation is required to simulate the effect that the magnetic susceptibility is not proportional to the magnetic field at low temperatures.In addition, the eigenvalues are often series expanded with respect to the magnetic field to simplify the equation [2,4]; however, this should be avoided whenever possible, as it may reduce the accuracy of the calculation when |D| is small.The field-dependent magnetic susceptibility equation without using series expansion was reported for the S = 1 system [11].In recent years, magnetization measurements have also become more common, and magnetization in the principal directions can be expressed in the same way as the magnetic field-dependent susceptibility equation.However, with powder samples, although the methodology is well established, it is not easy to express the powder average magnetization algebraically.In this study, we derive an algebraic expression for the exact powder average magnetization for the S = 1 system.
On the other hand, it is not easy to predict the sign of D when the symmetry of the ligand field around the central metal is not clear.(Please note that the symmetry of the ligand field [1] includes not only the position of the donor atoms but also their ligand-field strengths and the orientation of the donor π-orbitals.)In this study, nickel(II) complexes were prepared with a ligand, 2-[bis(2-methoxyethyl)aminomethyl]-4-nitrophenolate(1−) [(onp) − ] (Figure 1), giving low-symmetric ligand fields to investigate the relationship between the structure and zero-field splitting parameters of low-symmetric nickel(II) complexes.
Magnetochemistry 2024, 10, x FOR PEER REVIEW 2 of 17 the electron-paramagnetic-resonance (epr) signals to precisely determine the zero-field splitting parameters, D and E [3].Here, the axial zero-field-splitting parameter D and the rhombic zero-field-splitting parameter E are defined as the phenomenological Hamiltonian H = guβSuHu + D[Sz 2 − S(S + 1)/3] + E(Sx 2 − Sy 2 ) [4], which can also be derived from the formal spin-Hamiltonian [2].Then, in the S = 1 system, the three originally degenerate microstates are split into −2D/3, D/3 − E, and D/3 + E energy levels without an external magnetic field.On the other hand, attempts to determine the zero-field-splitting parameters from magnetic data analysis have been successful [2,5], although not as accurate as HFEPR.The relationship between the coordination geometry around the nickel(II) ion and the sign of the axial zero-field splitting parameter D has been studied [2], and the result indicates that D is negative when the axial ligand field is strong (ez > ex,y) and positive when the axial ligand field is weak (ez < ex,y), where eu represents the metal-ligand interaction in the u direction.Since the sign of the D parameter is related to magnetic anisotropy and the conditions for the single molecule magnets [2,[6][7][8][9], even if there are exceptions [10], it is important to clarify the relationship between the sign of D and structure and to control the sign of D.
To determine the zero-field splitting parameters from the magnetic analysis, simultaneous fitting of the temperature dependence of magnetic susceptibility and the field dependence of magnetization is preferable [2].For this purpose, it is, of course, essential to simulate magnetic susceptibility and magnetization.First, in the days when the cryomagnetic susceptibility was measured only down to about 5 K, it was not necessary to consider the effect of saturation by the applied magnetic field.However, now that it is common to measure magnetic susceptibility down to 2 K, a magnetic-field-dependent magnetic susceptibility equation is required to simulate the effect that the magnetic susceptibility is not proportional to the magnetic field at low temperatures.In addition, the eigenvalues are often series expanded with respect to the magnetic field to simplify the equation [2,4]; however, this should be avoided whenever possible, as it may reduce the accuracy of the calculation when |D| is small.The field-dependent magnetic susceptibility equation without using series expansion was reported for the S = 1 system [11].In recent years, magnetization measurements have also become more common, and magnetization in the principal directions can be expressed in the same way as the magnetic field-dependent susceptibility equation.However, with powder samples, although the methodology is well established, it is not easy to express the powder average magnetization algebraically.In this study, we derive an algebraic expression for the exact powder average magnetization for the S = 1 system.
On the other hand, it is not easy to predict the sign of D when the symmetry of the ligand field around the central metal is not clear.(Please note that the symmetry of the ligand field [1] includes not only the position of the donor atoms but also their ligand-field strengths and the orientation of the donor π-orbitals.)In this study, nickel(II) complexes were prepared with a ligand, 2-[bis(2-methoxyethyl)aminomethyl]-4-nitrophenolate(1−) [(onp) − ] (Figure 1), giving low-symmetric ligand fields to investigate the relationship between the structure and zerofield splitting parameters of low-symmetric nickel(II) complexes.In the ligand (onp) − , the bis(2-methoxyethyl)aminomethyl motif is introduced at the 2-position of 4-nitrophenol, while the 4-chloro derivative was reported [12].To our knowledge, the introduction of bis(2-methoxyethyl)aminomethyl motif into a phenol-based ligand was first conducted by Sakiyama to obtain 2,6-bis[bis(2-methoxyethyl)aminomethyl]-4-methylphenolate(1−) [(bomp) − ] [13].Subsequently, 4-chlorophenol and 4-nitrophenol derivatives were also synthesized [14].In a series of syntheses, the more electron-withdrawing the 4-substituent is, the slower the reaction proceeds, but the cleaner the synthesis.Especially with 4-nitrophenol, which was used as the ligand in this study with a large electron-withdrawing property, the 2-substituted (onp) − can be made easily and cleanly without 2,6-substituted (bonp) − contamination.In synthesis, it is desirable to follow the raw materials and products in the reaction mixture chromatographically (e.g., by thin-layer chromatography) because slight differences in conditions can cause the reaction to proceed differently.Since formaldehyde, which has been used in past papers, is carcinogenic and volatile, it is better to use paraformaldehyde, which is not volatile.Great care should be taken when handling reactants.

Measurements
Elemental analyses (C, H, and N) were performed at the Elemental Analysis Service Centre of Kyushu University.Powder IR spectra were recorded on a Nicolet iS50 FTIR spectrometer and liquid IR spectra on a Jasco FT/IR-4000 using NaCl plates.Electrospray ionization (ESI) mass spectra were measured on a Waters Quattro micro-API mass spectrometer in mixed solvents (acetonitrile:dimethyl sulfoxide:H 2 O = 4:1:1).Powder electronic spectra were measured on a Shimadzu UV-2600.Magnetic susceptibility measurements were performed with a Quantum Design MPMS-XL7 SQUID magnetometer.Temperature dependence of the susceptibility was measured in the temperature range from 1.9 to 300 K with a static field of 5 kOe.Field dependence of the magnetization was measured in the field range from 0 to 70 kOe at 2 K and 4 K.All data were corrected for the diamagnetism of the capsule and the samples by means of Pascal's constants [15].

Preparations
[Ni(onp)(dmso)(H 2 O)][BPh 4 ]•2dmso 1.To a methanolic solution (20 mL) of Ni(NO 3 ) 2 •6H 2 O (0.31 g, 1.1 mmol) was added H(onp) (0.28 g, 0.98 mmol) and refluxed for 15 min.To this was added sodium tetraphenylborate (0.69 g, 2.0 mmol).After cooling, green precipitates were collected by filtration to obtain crude product (0.24 g).The crude product is considered to be mainly In the ligand (onp) − , the bis(2-methoxyethyl)aminomethyl motif is introduced at the 2-position of 4-nitrophenol, while the 4-chloro derivative was reported [12].To our knowledge, the introduction of bis(2-methoxyethyl)aminomethyl motif into a phenolbased ligand was first conducted by Sakiyama to obtain 2,6-bis[bis(2-methoxyethyl)aminomethyl]-4-methylphenolate(1−) [(bomp) − ] [13].Subsequently, 4-chlorophenol and 4-nitrophenol derivatives were also synthesized [14].In a series of syntheses, the more electron-withdrawing the 4-substituent is, the slower the reaction proceeds, but the cleaner the synthesis.Especially with 4-nitrophenol, which was used as the ligand in this study with a large electron-withdrawing property, the 2-substituted (onp) − can be made easily and cleanly without 2,6-substituted (bonp) − contamination.In synthesis, it is desirable to follow the raw materials and products in the reaction mixture chromatographically (e.g., by thin-layer chromatography) because slight differences in conditions can cause the reaction to proceed differently.Since formaldehyde, which has been used in past papers, is carcinogenic and volatile, it is better to use paraformaldehyde, which is not volatile.Great care should be taken when handling reactants.

Measurements
Elemental analyses (C, H, and N) were performed at the Elemental Analysis Service Centre of Kyushu University.Powder IR spectra were recorded on a Nicolet iS50 FTIR spectrometer and liquid IR spectra on a Jasco FT/IR-4000 using NaCl plates.Electrospray ionization (ESI) mass spectra were measured on a Waters Quattro micro-API mass spectrometer in mixed solvents (acetonitrile:dimethyl sulfoxide:H2O = 4:1:1).Powder electronic spectra were measured on a Shimadzu UV-2600.Magnetic susceptibility measurements were performed with a Quantum Design MPMS-XL7 SQUID magnetometer.Temperature dependence of the susceptibility was measured in the temperature range from 1.9 to 300 K with a static field of 5 kOe.Field dependence of the magnetization was measured in the field range from 0 to 70 kOe at 2 K and 4 K.All data were corrected for the diamagnetism of the capsule and the samples by means of Pascal's constants [15].

Crystallography
Single-crystal X-ray diffraction data were obtained with a Rigaku XtaLAB AFC11 diffractometer with graphite-monochromated Mo Kα radiation (λ = 0.71073 Å).A single crystal was mounted with a cryoloop and flash-cooled with a cold N 2 gas stream.Data were processed using the CrysAlisPro 1.171.39.43a (Rigaku OD, 2018) software packages.The structure was solved by intrinsic phasing methods using the SHELXT 2018/2 (Sheldrick, 2018) [17] software packages and refined on F 2 (with all independent reflections) using the SHELXL 2018/3 (Sheldrick, 2018) [18] software packages.The non-hydrogen atoms were refined anisotropically, and hydrogen atoms, except for those of water molecules, were refined using the riding model.Hydrogen atoms for water molecules were located by Fourier difference map and fixed at the appropriate positions for hydrogen-bonding interactions during the refinement.

Results and Discussion
In this study, the nickel( , obtained as a green powder, was characterized as a bulk sample by IR, elemental analysis, mass spectrometry, electronic spectra, and magnetic properties.The single-crystal sample 2), prepared similarly to 1, but slowly, was obtained only in small quantities and was used only for single-crystal X-ray structural analysis.The related single-crystal sample [Ni(onp)(dmf) 2 ][BPh 4 ] (3) was also used only for X-ray structural analysis.
Crystals of 2 consists of [Ni(onp)(dmso)(H 2 O)] + complex cations and tetraphenylborate anions ([BPh 4 ] − ) in a 1:1 molar ratio.The crystal structure of [Ni(onp)(dmso)(H 2 O)] + complex cation is depicted in Figure 2. The selected bond distances and angles are summarized in Tables 2 and 3.In the complex cation, the (onp) − ligand worked as a tetradentate ligand, and the dmso and aqua ligands worked as monodentate ligands.The dmso ligand coordinated to the central nickel(II) ion through the oxygen atom, and the sulfur atom was disordered [occupancies: S(1), 0.7757 (19); S(1B), 0.2243 (19)].This is a typical phenomenon for an oxygen-coordinating dmso ligand due to the sulfur inversion [26,27].The coordination geometry around the central nickel(II) ion was octahedral with NO 5 donor atoms.The bond distances around the nickel(II) ion are normal for the related nickel(II) complexes [28,29], but the coordination geometry was slightly elongated along the O(3)-Ni(1)-O(6) direction.When the (onp) − ligand acts as a tetradentate ligand in an octahedral nickel(II) complex, the remaining two vacant sites should be in the cis-positions due to the steric requirement of the (onp) − ligand.Furthermore, the space of the vacant sites is not enough for the coordination of two dmso ligands.Therefore, the coordination of the aqua ligand with the dmso ligand is considered inevitable.
Crystals of 2 consists of [Ni(onp)(dmso)(H2O)] + complex cations and tetraphenylborate anions ([BPh4] − ) in a 1:1 molar ratio.The crystal structure of [Ni(onp)(dmso)(H2O)] + complex cation is depicted in Figure 2. The selected bond distances and angles are summarized in Tables 2 and 3.In the complex cation, the (onp) − ligand worked as a tetradentate ligand, and the dmso and aqua ligands worked as monodentate ligands.The dmso ligand coordinated to the central nickel(II) ion through the oxygen atom, and the sulfur atom was disordered [occupancies: S(1), 0.7757 (19); S(1B), 0.2243 (19)].This is a typical phenomenon for an oxygen-coordinating dmso ligand due to the sulfur inversion [26,27].The coordination geometry around the central nickel(II) ion was octahedral with NO5 donor atoms.The bond distances around the nickel(II) ion are normal for the related nickel(II) complexes [28,29], but the coordination geometry was slightly elongated along the O(3)-Ni(1)-O(6) direction.When the (onp) − ligand acts as a tetradentate ligand in an octahedral nickel(II) complex, the remaining two vacant sites should be in the cis-positions due to the steric requirement of the (onp) − ligand.Furthermore, the space of the vacant sites is not enough for the coordination of two dmso ligands.Therefore, the coordination of the aqua ligand with the dmso ligand is considered inevitable.] − anions in a 1:1 molar ratio.The crystal structure of [Ni(onp)(dmf) 2 ] + complex cation is shown in Figure 3.The selected bond distances and angles are summarized in Tables 4 and 5.In the complex cation, the (onp) − ligand worked as a tetradentate ligand, and the two dmf ligands worked as monodentate ligands.The coordination geometry around the central nickel(II) ion was octahedral with NO 5 donor atoms.The bond distances around the nickel(II) ion are normal for the related nickel(II) complexes [28,29], but the coordination geometry is slightly elongated along the O(5)-Ni(1)-O(6) direction.In both 2 and 3, the distortion patterns for the octahedral geometry are very similar to each other, although the exogenous ligands are different.In the case of the dmf ligand, the space required for coordination is smaller than that of the dmso ligand, allowing two dmf ligands to bind to the central nickel(II) ion.
Crystals of 3 consist of [Ni(onp)(dmf)2] + complex cations and [BPh4] − anions in a 1:1 molar ratio.The crystal structure of [Ni(onp)(dmf)2] + complex cation is shown in Figure 3.The selected bond distances and angles are summarized in Tables 4 and 5.In the complex cation, the (onp) − ligand worked as a tetradentate ligand, and the two dmf ligands worked as monodentate ligands.The coordination geometry around the central nickel(II) ion was octahedral with NO5 donor atoms.The bond distances around the nickel(II) ion are normal for the related nickel(II) complexes [28,29], but the coordination geometry is slightly elongated along the O(5)-Ni(1)-O(6) direction.In both 2 and 3, the distortion patterns for the octahedral geometry are very similar to each other, although the exogenous ligands are different.In the case of the dmf ligand, the space required for coordination is smaller than that of the dmso ligand, allowing two dmf ligands to bind to the central nickel(II) ion.To gain a better understanding of the ground state, density functional theory (DFT) computations were conducted for the crystallographically identified complex cations in 2 and 3.The energy diagrams and d-orbital-related molecular orbitals for 2 and 3 are shown in Figures 4 and 5, respectively.The results for the two complex cations were very similar, and we will first discuss the [Ni(onp)(dmso)(H 2 O)] + complex cation in 2. Judging from the splitting pattern of the five d-orbital-related molecular orbitals, the coordination geometry can be well approximated as the D 4 rotation group.(Please note that the D 4 rotation group [30] is a subgroup of the D 4h point group but is defined as a group considering only rotation operations.The D 4 rotation group is the simplest group suitable for reproducing the splitting pattern and the phase of the orbitals based on the coordination geometry).The five orbitals can be approximately assigned to b 2 , e, a 1 , and b 1 from the lower energy, where the e orbital set is doubly degenerate.The three lower orbitals (b 2 and e) are each filled with two electrons, but the two upper orbitals (a 1 and b 1 ) are each singly occupied and are magnetic orbitals.One important point is that the b 1 (d x 2 −y 2 ) orbital is above the a 1 (d z 2 ) orbital in the magnetic orbitals, indicating that the ligand field in the z direction is weaker than those in the x and y directions, which is consistent with the elongation along the O(3)-Ni(1)-O(6) direction (z direction) in the crystal structure.The other important point is the degenerate e orbital set.Although not completely degenerate due to the mixing of the ligand orbitals, the spread of the two orbitals shows linear combinations of the d xz and d yz orbitals.This indicates that the magnitudes of the ligand fields in the x and y directions are approximately the same.The characteristics of the energy diagram and molecular orbitals for 3 are basically the same as those of 2. The ligand field in the z direction is weaker than those in the x and y directions, but the ligand fields in the x and y directions are approximately the same.

Electronic Spectra of 1
The powder electronic spectra were measured for 1 by the diffuse reflection method at room temperature, as shown in Figure 6a.The spectra were typical of an octahedral nickel(II) complex with two spin-allowed bands from near-infrared to visible range.Assuming the O symmetry, the first spin-allowed band at around 10,000 cm −1 was assigned to 3 A2( 3 F) → 3 T2( 3 F), a small shoulder band at around 13,000 cm −1 to 3 A2 → 1 E( 1 D), the second spin-allowed band at around 16,000 cm −1 to 3 A2 → 3 T1( 3 F).The expected third spinallowed band [ 3 A2 → 3 T1( 3 P)] was hidden in the bands originating from the ligand.Focusing on the band shapes of the first and the second spin-allowed bands, the first band is slightly inclined toward the lower energy (the lower-energy side is steeper than the other), but the second band is slightly inclined toward the higher energy (higher energy side is steeper than the other).This feature was clearly reaffirmed by Gaussian curve fitting, assuming three spectral components of equal height and width.(Please note that while these

Electronic Spectra of 1
The powder electronic spectra were measured for 1 by the diffuse reflection method at room temperature, as shown in Figure 6a.The spectra were typical of an octahedral nickel(II) complex with two spin-allowed bands from near-infrared to visible range.Assuming the O symmetry, the first spin-allowed band at around 10,000 cm −1 was assigned to 3 A2( 3 F) → 3 T2( 3 F), a small shoulder band at around 13,000 cm −1 to 3 A2 → 1 E( 1 D), the second spin-allowed band at around 16,000 cm −1 to 3 A2 → 3 T1( 3 F).The expected third spinallowed band [ 3 A2 → 3 T1( 3 P)] was hidden in the bands originating from the ligand.Focusing on the band shapes of the first and the second spin-allowed bands, the first band is slightly inclined toward the lower energy (the lower-energy side is steeper than the other), but the second band is slightly inclined toward the higher energy (higher energy side is steeper than the other).This feature was clearly reaffirmed by Gaussian curve fitting, assuming three spectral components of equal height and width.(Please note that while these

Electronic Spectra of 1
The powder electronic spectra were measured for 1 by the diffuse reflection method at room temperature, as shown in Figure 6a.The spectra were typical of an octahedral nickel(II) complex with two spin-allowed bands from near-infrared to visible range.Assuming the O symmetry, the first spin-allowed band at around 10,000 cm −1 was assigned to 3 A 2 ( 3 F) → 3 T 2 ( 3 F), a small shoulder band at around 13,000 cm −1 to 3 A 2 → 1 E( 1 D), the second spin-allowed band at around 16,000 cm −1 to 3 A 2 → 3 T 1 ( 3 F).The expected third spin-allowed band [ 3 A 2 → 3 T 1 ( 3 P)] was hidden in the bands originating from the ligand.Focusing on the band shapes of the first and the second spin-allowed bands, the first band is slightly inclined toward the lower energy (the lower-energy side is steeper than the other), but the second band is slightly inclined toward the higher energy (higher energy side is steeper than the other).This feature was clearly reaffirmed by Gaussian curve fitting, assuming three spectral components of equal height and width.(Please note that while these restrictions may not be strictly correct, the ability to uniquely determine the solution in the initial process of fitting is an advantage.)The obtained spectral components are also shown in Figure 6a, and the sum of the spectral components fits the observed data very well.This splitting pattern is typical of tetragonally elongated octahedral geometry with a weak axial ligand field [31].
the Racah parameters, and eσ,ax and eσ,eq are the axial (z direction) and the equatorial (x, y direction) AOM parameters, corresponding to the metal-ligand σ interactions [1].The π interactions were not considered in the calculation.The resulting energy levels of the states are shown in Figure 6b, together with the assignment of states in the D4 symmetry.The AOM parameters obtained from the spectral analysis show that the ligand field of complex 1 is D4 symmetric and axially weak (eσ,ax < eσ,eq), which is consistent with the DFT results for 2 and 3.The D4 symmetry and the axially weak ligand field around the central nickel(II) ion may be due to the nature of the (omp) − ligand.Furthermore, it is reasonable to assume that the structure of the complex cation in 1 is the same as that in 2 since the steric requirements of the (omp) − ligand prevent the two dmso ligands from binding to the nickel(II) ion, as discussed in Section 3.

Theory and Magnetic Equations
Prior to magnetic analysis, the equations necessary for magnetic analysis are organized and summarized here.In particular, closed-form algebraic equations describing the magnetization in any field direction in the S = 1 system will be derived in this study.This allows the powder magnetization to be expressed in terms of algebraic equations.
In molecular magnetism, the following procedure is generally used to simulate magnetization and magnetic susceptibility [4,32].First, the Hamiltonian is defined, and then the energy matrix is obtained from the Hamiltonian and the appropriate set of basis functions.Solving the secular equation yields the energy of the microstates as eigenvalues.Furthermore, partial differentiation of the obtained eigenvalues with respect to the magnetic field yields the microscopic magnetization of the microstates.Then, by considering the Boltzmann distribution, the equations for magnetization and magnetic susceptibility can be obtained.In the process of solving the secular equation to obtain eigenvalues, the dimension of the matrix is often too large to obtain an exact solution, but in the present S = 1 system, the dimension of the matrix is only 3 × 3, so the exact solution can be obtained algebraically.Next, angular overlap model (AOM) calculations were performed assuming D 4 symmetry to reproduce the obtained spectral components.The best-fitting parameter set was obtained as (B, C, e σ,ax , e σ,eq ) = (920 cm −1 , 3200 cm −1 , 2690 cm −1 , 3550 cm −1 ), where B and C are the Racah parameters, and e σ,ax and e σ,eq are the axial (z direction) and the equatorial (x, y direction) AOM parameters, corresponding to the metal-ligand σ interactions [1].The π interactions were not considered in the calculation.The resulting energy levels of the states are shown in Figure 6b, together with the assignment of states in the D 4 symmetry.The AOM parameters obtained from the spectral analysis show that the ligand field of complex 1 is D 4 symmetric and axially weak (e σ,ax < e σ,eq ), which is consistent with the DFT results for 2 and 3.The D 4 symmetry and the axially weak ligand field around the central nickel(II) ion may be due to the nature of the (omp) − ligand.Furthermore, it is reasonable to assume that the structure of the complex cation in 1 is the same as that in 2 since the steric requirements of the (omp) − ligand prevent the two dmso ligands from binding to the nickel(II) ion, as discussed in Section 3.1.1.

Theory and Magnetic Equations
Prior to magnetic analysis, the equations necessary for magnetic analysis are organized and summarized here.In particular, closed-form algebraic equations describing the magnetization in any field direction in the S = 1 system will be derived in this study.This allows the powder magnetization to be expressed in terms of algebraic equations.
In molecular magnetism, the following procedure is generally used to simulate magnetization and magnetic susceptibility [4,32].First, the Hamiltonian is defined, and then the energy matrix is obtained from the Hamiltonian and the appropriate set of basis functions.Solving the secular equation yields the energy of the microstates as eigenvalues.Furthermore, partial differentiation of the obtained eigenvalues with respect to the magnetic field yields the microscopic magnetization of the microstates.Then, by considering the Boltzmann distribution, the equations for magnetization and magnetic susceptibility can be obtained.In the process of solving the secular equation to obtain eigenvalues, the dimension of the matrix is often too large to obtain an exact solution, but in the present S = 1 system, the dimension of the matrix is only 3 × 3, so the exact solution can be obtained algebraically.
Since powder samples were used for magnetic measurements, the powder average in Equation (1) [2,4] must be used for simulations in magnetic susceptibility data analysis for magnetically anisotropic compounds.
However, since the magnetic susceptibility in weak magnetic fields is additive, the powder average is the same as the arithmetic average.When the symmetry is axial, the susceptibility is written as in Equation ( 2), which is the same as Equation (A1) in Appendix A.
To calculate the principal magnetic susceptibilities χ z and χ x in Equation ( 2), previously derived Equations (A2) and (A3) [11] in Appendix A are useful for the S = 1 system.Considering the axial anisotropic g-factors (g z and g x ), the axial zero-field splitting (D), temperature-independent paramagnetism (χ tip,z and χ tip,x ), and intermolecular interaction (zJ), the atomic magnetic susceptibility, χ A , and its components are expressed as Equations (A4)-(A7) [11] in Appendix A.
The problem here is the magnetization with anisotropy.While the principal magnetization (M x , M y , and M z ) can be easily expressed, the powder average can be expressed as the arithmetic average of the principal magnetization only in weak magnetic fields where the magnetization is proportional to the field.The principal magnetization M z and M x in axial symmetry can be expressed as Equations (A8) and (A9) in Appendix B. The powder average magnetization required for magnetization simulation in high magnetic fields can be written as in Equation ( 3).In Equation ( 3), the function M(θ, φ) represents the magnetization in the u(θ, φ) direction of the polar coordinates when the magnetic field is applied along the direction u(θ, φ).To calculate the exact powder average of the magnetization, we first need to express the function M(θ, φ).
Fortunately, this matrix can be solved algebraically to obtain three eigenvalues, E u,1 , E u,2 , and E u,3 , as closed-form expressions.Hence, the eigenvalues are partially differentiable with respect to the magnetic field, and the microscopic magnetization [µ u,n = −∂E u,n ⁄∂H u (n = 1, 2, 3)] can also be expressed algebraically.Thus, in axial symmetry, the magnetization in any u(θ) direction, M θ (θ), can be expressed as a closed-form algebraic expression.The obtained expressions are given as Equations (A10)-(A27) in Appendix C, i.e., the magnetic anisotropy is expressed as the θ-dependent (orientation-dependent) magnetization affected by the zero-field splitting (D) and g-anisotropy.
Now the function M(θ, φ) in Equation ( 3) has been obtained as the function M θ (θ) in the axial symmetry (Equation (A27)), and the powder average can be approximated, using numerical integration [11,34].In Equation ( 4), exact integration is conducted with respect to φ, and numerical integration is conducted with respect to θ in 1.0 • steps from 0.5 • to 89.5 • .This equation gives a highly accurate approximation of the powder average.However, if more precision is needed, Equation ( 5) can be used instead of Equation (4), for example.

Magnetic Properties of Complex 1
As discussed in Section 3.3, the structure of the complex cation in 1 is considered to be the same as that in 2. Temperature dependence of magnetic susceptibility (χ A ) and field dependence of magnetization (M) were measured for 1.The χ A T versus T and the M versus H plots for 1 is shown in Figure 7.The χ A T value at 300 K was 1.43 cm 3 K mol −1 , which was slightly larger than the spin-only value (1.00 cm 3 K mol −1 ) for the S = 1 system assuming g = 2.When decreasing the temperature, the χ A T value was almost constant, but below 15 K, the value dropped to show the minimum (0.62 cm 3 K mol -1 ) at 1.9 K.This dropping behavior can be explained by very small antiferromagnetic intermolecular interaction, but also by the zero-field splitting of the ground state.As the magnetic field increases, the M/Nβ value is going to saturate to ~2.0, which is consistent with the S = 1 system.equation gives a highly accurate approximation of the powder average.However, if more precision is needed, Equation ( 5) can be used instead of Equation ( 4 (5)

Magnetic Properties of Complex 1
As discussed in Section 3.3, the structure of the complex cation in 1 is considered to be the same as that in 2. Temperature dependence of magnetic susceptibility (χA) and field dependence of magnetization (M) were measured for 1.The χAT versus T and the M versus H plots for 1 is shown in Figure 7.The χAT value at 300 K was 1.43 cm 3 K mol −1 , which was slightly larger than the spin-only value (1.00 cm 3 K mol −1 ) for the S = 1 system assuming g = 2.When decreasing the temperature, the χAT value was almost constant, but below 15 K, the value dropped to show the minimum (0.62 cm 3 K mol -1 ) at 1.9 K.This dropping behavior can be explained by very small antiferromagnetic intermolecular interaction, but also by the zero-field splitting of the ground state.As the magnetic field increases, the M/Nβ value is going to saturate to ~2.0, which is consistent with the S = 1 system.
In this study, the field-dependent magnetic susceptibility equations (Equations (A1)-(A7) in Appendix A) [11] were used for simulating magnetic susceptibility, and magnetization equations (Equations (A10)-(A27) in Appendix C and Equation ( 4)) for simulating magnetization, considering the axial anisotropic g-factors (gz and gx), the axial zero-field splitting (D), temperature-independent paramagnetism (χtip), and intermolecular interaction (zJ).In the magnetic analysis, the χtip-anisotropy was not considered.As the result of simultaneous analysis, the best-fitting parameter set was obtained as (gz, gx, D, zJ, χtip) = (2.13,2.31, +5.9 cm −1 , −0.70 cm −1 , 534 × 10 −6 cm 3 mol −1 ) at 5 kOe.When χAT decreases at very low temperatures (below 15 K), separating zJ and D is not easy, but it should be possible to separate them because their temperature dependencies are generally different.This time, we were able to find the best combination of zJ and D by simultaneous fitting of magnetic susceptibility and magnetization due to the fact that the difference between the theoretical magnetization curves for 2 K and 4 K is larger for larger In this study, the field-dependent magnetic susceptibility equations (Equations (A1)-(A7) in Appendix A) [11] were used for simulating magnetic susceptibility, and magnetization equations (Equations (A10)-(A27) in Appendix C and Equation ( 4)) for simulating magnetization, considering the axial anisotropic g-factors (g z and g x ), the axial zero-field splitting (D), temperature-independent paramagnetism (χ tip ), and intermolecular interaction (zJ).In the magnetic analysis, the χ tip -anisotropy was not considered.
As the result of simultaneous analysis, the best-fitting parameter set was obtained as (g z , g x , D, zJ, χ tip ) = (2.13, 2.31, +5.9 cm −1 , −0.70 cm −1 , 534 × 10 −6 cm 3 mol −1 ) at 5 kOe.When χ A T decreases at very low temperatures (below 15 K), separating zJ and D is not easy, but it should be possible to separate them because their temperature dependencies are generally different.This time, we were able to find the best combination of zJ and D by simultaneous fitting of magnetic susceptibility and magnetization due to the fact that the difference between the theoretical magnetization curves for 2 K and 4 K is larger for larger −zJ.Determining the sign of D is also difficult, but reversing the sign of D worsens the fitting of the theoretical magnetization curve to the measured values.The positive D value suggests that the axial ligand field is weaker than the equatorial ligand field [2], which is concordant with the result of powder electronic spectra for 1 (Section 3.3) and also with the DFT result for 2, whose cation structure is considered to be the same (Section 3.2).Although it is difficult to obtain precise g-values from magnetic data, considering g-anisotropy improves the fitting of magnetization in the high-field region.The resulting anisotropic g-factor relationship (g z < g x ) agreed with the positive D value [2].The magnetization curves are highly dependent on g-anisotropy, and this study shows the importance of considering g-anisotropy in the analysis, even if g-anisotropy is small.
Based on the obtained magnetic parameters, the three-dimensional magnetic anisotropy can be represented in the image using Equation (A27).First, the spatial distribution of magnetization in the xz cross section when the magnetic field is increased from 5 kOe to 100 kOe is shown in Figure 8a, where the z direction corresponds to θ = 0 • , and the x-direction to θ = 90 • .When the magnetic field is weak (5-20 kOe), magnetization occurs in the xy-plane (easy magnetization plane) direction, and when the magnetic field becomes stronger (60-100 kOe), magnetization increases also in the z-axis (hard magnetization axis) direction.This is consistent with the relatively stronger equatorial ligand field originating from the tetradentate ligand.Next, a three-dimensional representation of magnetization at 5 kOe, the magnetic field in magnetic susceptibility measurement, is depicted in Figure 8b.The shape is a biconcave disk with a depressed center in the z direction, which is different from the shape of the thermal ellipsoid in crystallography.Thus, the observed magnetic anisotropy is characterized as easy-plane anisotropy.
−zJ.Determining the sign of D is also difficult, but reversing the sign of D worsens the fitting of the theoretical magnetization curve to the measured values.The positive D value suggests that the axial ligand field is weaker than the equatorial ligand field [2], which is concordant with the result of powder electronic spectra for 1 (Section 3.3) and also with the DFT result for 2, whose cation structure is considered to be the same (Section 3.2).Although it is difficult to obtain precise g-values from magnetic data, considering g-anisotropy improves the fitting of magnetization in the high-field region.The resulting anisotropic g-factor relationship (gz < gx) agreed with the positive D value [2].The magnetization curves are highly dependent on g-anisotropy, and this study shows the importance of considering g-anisotropy in the analysis, even if g-anisotropy is small.
Based on the obtained magnetic parameters, the three-dimensional magnetic anisotropy can be represented in the image using Equation (A27).First, the spatial distribution of magnetization in the xz cross section when the magnetic field is increased from 5 kOe to 100 kOe is shown in Figure 8a, where the z direction corresponds to θ = 0°, and the xdirection to θ = 90°.When the magnetic field is weak (5-20 kOe), magnetization occurs in the xy-plane (easy magnetization plane) direction, and when the magnetic field becomes stronger (60-100 kOe), magnetization increases also in the z-axis (hard magnetization axis) direction.This is consistent with the relatively stronger equatorial ligand field originating from the tetradentate ligand.Next, a three-dimensional representation of magnetization at 5 kOe, the magnetic field in magnetic susceptibility measurement, is depicted in Figure 8b.The shape is a biconcave disk with a depressed center in the z direction, which is different from the shape of the thermal ellipsoid in crystallography.Thus, the observed magnetic anisotropy is characterized as easy-plane anisotropy.

Powder Average Magnetization and Arithmetic Average Magnetization
For future reference, the present study compares the powder average magnetization Mav,p and the arithmetic average magnetization Mav,a for the S = 1 system.The theoretical magnetization curves are drawn in Figures 9 and 10, using the equations in Appendices B and C and Equation (4).In both cases in Figure 9, gz = gx = 2.00 is assumed in axial symmetry, and Mz and Mx are drawn together.In Figure 9a, D is assumed to be +10 cm −1 , and in Figure 9b, D is assumed to be −10 cm −1 .(Please note that at D = 0 cm −1 , Mz, Mx, Mav,a, and Mav,p are all equal to the Brillouin function).In both Figure 9a,b, each Mav,a curve is close to each Mav,p curve in the low field region.This is concordant with the fact that the arithmetic average magnetic susceptibility is equal to the powder average as long as the magnetic field is weak enough to maintain additivity.

Powder Average Magnetization and Arithmetic Average Magnetization
For future reference, the present study compares the powder average magnetization M av,p and the arithmetic average magnetization M av,a for the S = 1 system.The theoretical magnetization curves are drawn in Figures 9 and 10, using the equations in Appendices B and C and Equation (4).In both cases in Figure 9, g z = g x = 2.00 is assumed in axial symmetry, and M z and M x are drawn together.In Figure 9a, D is assumed to be +10 cm −1 , and in Figure 9b, D is assumed to be −10 cm −1 .(Please note that at D = 0 cm −1 , M z , M x , M av,a , and M av,p are all equal to the Brillouin function).In both Figure 9a,b, each M av,a curve is close to each M av,p curve in the low field region.This is concordant with the fact that the arithmetic average magnetic susceptibility is equal to the powder average as long as the magnetic field is weak enough to maintain additivity.
meaning that the arithmetic average and the powder average are equal.However, since saturation magnetization is proportional to the g-value, it is no longer additive, and the arithmetic average is not a good approximation for a strong magnetic field.It is important to consider g-anisotropy in the simulation because it affects the powder average, especially at high magnetic fields.

Conclusions
In this study, octahedral nickel(II) complexes were prepared with (omp) − ligand.Complex 1 was characterized as a bulk sample, and its powder electronic spectra were analyzed based on the AOM calculation, which concluded that the spectra were typical of meaning that the arithmetic average and the powder average are equal.However, since saturation magnetization is proportional to the g-value, it is no longer additive, and the arithmetic average is not a good approximation for a strong magnetic field.It is important to consider g-anisotropy in the simulation because it affects the powder average, especially at high magnetic fields.

Conclusions
In this study, octahedral nickel(II) complexes were prepared with (omp) − ligand.Complex 1 was characterized as a bulk sample, and its powder electronic spectra were analyzed based on the AOM calculation, which concluded that the spectra were typical of The significant g-anisotropy cases are simulated, as shown in Figure 10.In the weak magnetic field region, the arithmetic average magnetization is equal to the powder average, as in the g-isotropic case.The difference from the g-isotropic cases is the behavior in the high-field region where the magnetization approaches saturation, i.e., the difference between the arithmetic average and the powder average is larger at high magnetic fields when the g-anisotropy is large.The effect of the g-factor on the magnetization curves can be understood by the additivity of the square of the g-value.In terms of the g-value, the g-value is not additive, but the square of the g-value is additive.The slope of magnetization (magnetic susceptibility) in a weak magnetic field, where magnetization is proportional to the magnetic field, is proportional to the square of the g-value, and is additive, meaning that the arithmetic average and the powder average are equal.However, since saturation magnetization is proportional to the g-value, it is no longer additive, and the arithmetic average is not a good approximation for a strong magnetic field.It is important to consider g-anisotropy in the simulation because it affects the powder average, especially at high magnetic fields.

Figure 6 .
Figure 6.Powder electronic spectra and energy levels of 1: (a) observed spectra (-), spectral components (•••••), and sum of the components (---); (b) energy levels of the states in the D4 symmetry obtained from the spectral analysis.

Figure 6 .
Figure 6.Powder electronic spectra and energy levels of 1: (a) observed spectra (-), spectral components (•••••), and sum of the components (---); (b) energy levels of the states in the D 4 symmetry obtained from the spectral analysis.

Figure 9 .
Figure 9.Comparison of powder average magnetization (M av,p ) and arithmetic average magnetization (M av,a ) in the g-isotropic cases (g = 2.00) at 2 K. M z and M x curves are also drawn: (a) when D = +10 cm −1 ; (b) when D = −10 cm −1 .

Figure 10 .
Figure 10.Comparison of powder average magnetization (M av,p ) and arithmetic average magnetization (M av,a ) when g-anisotropy is significant at 2 K. M z and M x curves are also drawn: (a) when D = +10 cm −1 , g z = 1.41, and g x = 2.24; (b) when D = −10 cm −1 , g z = 2.83, and g x = 1.41.

Table 1 .
Crystallographic data and refinement parameters of 2 and 3.