research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoSTRUCTURAL SCIENCE
CRYSTAL ENGINEERING
MATERIALS
ISSN: 2052-5206

On the magnetic and crystal structures of NiO and MnO

crossmark logo

aLaboratory for Neutron Scattering and Imaging, Paul Scherrer Institut (PSI), 5232 Villigen PSI, Switzerland
*Correspondence e-mail: vladimir.pomjakushin@psi.ch

Edited by V. Garlea, Oak Ridge National Laboratory, USA (Received 13 May 2024; accepted 31 July 2024; online 10 September 2024)

This article is part of a focused issue on Magnetic Structures.

The magnetic and crystal structures of manganese and nickel monoxides have been studied using high-resolution neutron diffraction. The known 1k-structures based on the single propagation vector [½ ½ ½] for the parent paramagnetic space group Fm3m are forced to have monoclinic magnetic symmetry and are not possible in rhombohedral symmetry. However, the monoclinic distortions from the rhombohedral crystal metric allowed by symmetry are very small, and the explicit monoclinic splittings of the diffraction peaks have not been experimentally observed. We analyse the magnetic crystallographic models metrically compatible with our experimental data in full detail by using isotropy subgroup representation approach, including rhombohedral solutions based on the propagation vector star {[½ ½ ½], [−½ ½ ½], [½−½ ½], [½ ½ −½]}. Although the full star rhombohedral RI3c structure can equally well fit our diffraction data for NiO, we conclude that the best solution for the crystal and magnetic structures for NiO and MnO is the 1k monoclinic model with the magnetic space group Cc2/c (Belov–Neronova–Smirnova No. 15.90, UNI symbol C2/c.1′c[C2/m]).

1. Introduction

Since the first neutron diffraction studies of manganese and nickel monoxides (Shull et al., 1951[Shull, C. G., Strauser, W. A. & Wollan, E. O. (1951). Phys. Rev. 83, 333-345.]), interest in the magnetic and crystal structures of MnO and NiO has been maintained [see e.g. experimental studies (Cheetham & Hope, 1983[Cheetham, A. K. & Hope, D. A. O. (1983). Phys. Rev. B, 27, 6964-6967.]; Lee et al., 2016[Lee, S., Ishikawa, Y., Miao, P., Torii, S., Ishigaki, T. & Kamiyama, T. (2016). Phys. Rev. B, 93, 064429.]) and theoretical first principal calculations (Schrön et al., 2012[Schrön, A., Rödl, C. & Bechstedt, F. (2012). Phys. Rev. B, 86, 115134.]; Lim et al., 2016[Lim, J. S., Saldana-Greco, D. & Rappe, A. M. (2016). Phys. Rev. B, 94, 165151.]) and references cited therein]. It is well established experimentally that the cubic paramagnetic crystal structure [{Fm\bar{3}m}] transforms metrically to rhombo­hedral [{R\bar{3}m}] at the transition to an antiferromagnetically ordered structure below the Néel temperature TN. The antiferromagnetic structure is well described with one propagation vector k1 = [½  ½  ½] (in cubic crystal metric) and the moment directed perpendicular to k1. This type of structure corresponds to the monoclinic magnetic space group (MSG) Cc2/m [No. 12.63 Belov–Neronova–Smirnova (BNS)], whose properties and relation to the parent cubic group are shown in Table 1[link]. The second high-symmetry monoclinic MSG is Cc2/c (No. 15.90) and it allows spin components both parallel and perpendicular to k1 (the x-component in Cc2/c settings is perpendicular to cubic k1, and the z-component makes an angle ∼35° to the k1 direction). Experimentally, one does see rhombohedral distortion in the title compounds, but the monoclinic distortions from the [{R\bar{3}m}] crystal metric allowed by symmetry are very small and the explicit monoclinic splittings of the diffraction peaks have not been experimentally observed. Thus, a quite logical idea would be try to find the solutions in rhombohedral symmetry that would fit the experimental data. The only rhombohedral MSG permitting non-zero magnetic moments on Ni/Mn for a one-k1 structure is [{R_{I}\bar{3}c}] (1k) with a single Wyckoff position for Mn, but this structure cannot fit experimental diffraction intensities at all because it gives zero intensity for the most intense magnetic peak (½ ½ ½). However, the rhombohedral structures based on several propagation vectors (arms) from the propagation vector star {k1} = {[½ ½ ½], [−½ ½ ½], [½ −½ ½], [½ ½ −½]} could, in principle, be compatible with experimental data. In this paper we study possible rhombohedral and monoclinic solutions for multi-k magnetic structures based on several arms from the star {k1} and compare the models with high-resolution neutron diffraction data.

Table 1
The crystal and magnetic structure parameters in NiO and MnO at T = 2 K for the 1k monoclinic MSG Cc2/c (BNS No. 15.90, UNI C2/c.1′c[C2/m]) and multi-k [R_{I}\bar{3}] (BNS No. 148.20, UNI [R\bar{3}.1^{\prime}_{c}[R\bar{3}]]) structure for NiO

The conventional reliability R factors (%) (Rodríguez-Carvajal, 1993[Rodríguez-Carvajal, J. (1993). Physica B, 192, 55-69.]) are also given after SG symbol for the fits with λ = 1.155 Å. The reliability factors on the separate lines, labelled Rp, Rwp, Rexp and χ2 are for the fits with λ = 1.886 Å, aimed at detailing the magnetic structure, with the magnetic moment components provided in the previous line. The fits in rhombohedral space group [R\bar{3}m] with the propagation vector [{[{0}{0}{3\over 2}]}] for the magnetic structures give an equivalent description of the structures. The metric transformation from cubic to monoclinic cells reads: (½, ½, −1), (½, −½, 0), (−1, −1, 0) and to 1k rhombohedral reads: (−½, ½, 0), (0, −½, ½), (1, 1, 1). The monoclinic unit-cell parameters derived from the rhombohedral metric are shown after the comma. Ni and O occupy 3a (000) and 3b (00½) Wyckoff positions in space group [R\bar{3}m], respectively. The magnitude of the magnetic moment is denoted by m, with mx, my and mz representing its components along their respective axes of their MSGs. The magnetic moment has both components perpendicular and parallel to the propagation vector. In the magnetic space group Cc2/c, Ni and O are in 4c (000 | mx,0,mz) and 4b (0 ½ ¼) positions. The moment directions in monoclinic groups along x, y, and z axes correspond to (½, ½, −1), (1, −1, 0) and (−1, −1, 0) in cubic metric, respectively. The Wyckoff positions in the multi-k [R_{I}\bar{3}] MSG are 18e (½, 0, 0) and 6a (0, 0, 0) for Ni1 and Ni2 and 18d (½, 0, ¼) and 6b (0, 0, ¼) for O1 and O2, respectively. Here m is the size of the magnetic moment; θ and ϕ are spherical angles with c and a axes, respectively in MSGs [R_{I}\bar{3}].

  NiO MnO
[R\bar{3}m] χ2 8.94, 9.59, 5.84, 2.70 8.71, 9.87, 7.18, 1.89
a (Å) 2.95006 (3) 3.14980 (7)
c (Å) 7.21235 (9) 7.59158 (16)
mx, mz, m (μB) 1.840 (9), 0.63 (3), 1.943 (9) 4.75 (1), 0, 4.75 (1)
Cc2/c χ2 8.84, 9.01, 5.87, 2.36 8.46, 9.54, 7.14, 1.79
a (Å) 5.10263 (8), 5.10967 5.45583 (15), 5.45657
b (Å) 2.95045 (6), 2.95007 3.15006 (9), 3.15035
c (Å) 5.89158 (11), 5.89263 6.22984 (15), 6.233451
β (°) 125.1722 (9), 125.3159 125.6659 (13), 125.7027
mx, mz, m (μB) 2.28 (2), 0.77 (4), 1.946 (9) 4.61 (1), 0, 4.61 (1)
Rp,  Rwp, Rexp, χ2 3.67, 5.10, 2.31, 4.86 4.51, 6.04, 3.18, 3.61
  SA search for NiO in MSG [R_{I}\bar{3}]
Magnetic moments Ni1 (18e) Ni2 (6a)
mx, my, mz, m (μB) 0.5582, 1.1093 −1.6018, 1.8678 0, 0, 2.0484, 2.0484
  −0.5386, −1.1026, 0.2224, 0.9805 0, 0, −3.4317, 3.4317
R factors 3.65, 5.27, 2.31, 5.19  
 
  LSQ fit in MSG [R_{I}\bar{3}]
a, c (Å) 5.90046 (6), 14.42656 (18)
m (μB), ϕ, θ (°) 1.976 (9), 90, 35 (5) 1.977 (9), 0, 180
mx, my, mz (μB) −0.65 (8), −1.31 (17), 1.6 (1) −1.977 (9) 0 0
Rp,  Rwp, Rexp, χ2 3.72, 5.37, 2.31, 5.38  

2. Experimental

Neutron powder diffraction measurements were carried out on the high-resolution powder diffractometer for thermal neutrons (HRPT) (Fischer et al., 2000[Fischer, P., Frey, G., Koch, M., Könnecke, M., Pomjakushin, V., Schefer, J., Thut, R., Schlumpf, N., Bürge, R., Greuter, U., Bondt, S. & Berruyer, E. (2000). Physica B, 276-278, 146-147.]) at the SINQ neutron spallation source (Paul Scherrer Institute, Switzerland), using high-resolution mode (δd/d > 10−3) with neutron wavelength λ = 1.155 Å and high-intensity mode with neutron wavelength λ = 1.886 Å. The refinements of the structure parameters were carried out using the FullProf (Rodríguez-Carvajal, 1993[Rodríguez-Carvajal, J. (1993). Physica B, 192, 55-69.]) program, with the use of its internal tables for neutron scattering lengths. The symmetry analysis was performed using the ISODISTORT tool based on the ISOTROPY (Stokes & Hatch, 1988[Stokes, H. T. & Hatch, D. M. (1988). Isotropy Subgroups of the 230 Crystallographic Space Groups. Singapore: World Scientific.]; Campbell et al., 2006[Campbell, B. J., Stokes, H. T., Tanner, D. E. & Hatch, D. M. (2006). J. Appl. Cryst. 39, 607-614.]) software, the software tools of Bilbao Crystallographic Server (Aroyo et al., 2011[Aroyo, M., Perez-Mato, J., Orobengoa, D., Tasci, E., De La Flor, G. & Kirov, A. (2011). Bulg. Chem. Commun. 43, 183.]) and the BasiRep program (Rodríguez-Carvajal, 1993[Rodríguez-Carvajal, J. (1993). Physica B, 192, 55-69.]). The magnetic structure figures were prepared using the VESTA (Momma & Izumi, 2011[Momma, K. & Izumi, F. (2011). J. Appl. Cryst. 44, 1272-1276.]) software program.

3. Symmetry analysis

The paramagnetic crystal structure of MnO and NiO is [Fm\bar{3}m] with Mn or Ni in the 4a (0,0,0) position and oxygen in the 4b (½, ½, ½) position. All magnetic peaks are indexed with the propagation vector k1 = [½ ½ ½], which is the L-point of the Brillouin zone. The decomposition of the magnetic representation for the 4a position into irreducible representations (irreps) reads [{\rm mL}_{2}^{+}(\tau_{3})\oplus {\rm mL}_{3}^{+}(\tau_{5})] with the dimensions of irreps being 1D and 2D, respectively. The nomenclature for the irreps is given in accordance with Campbell et al. (2006[Campbell, B. J., Stokes, H. T., Tanner, D. E. & Hatch, D. M. (2006). J. Appl. Cryst. 39, 607-614.]) and Aroyo et al. (2011[Aroyo, M., Perez-Mato, J., Orobengoa, D., Tasci, E., De La Flor, G. & Kirov, A. (2011). Bulg. Chem. Commun. 43, 183.]), with Kovalev (1993[Kovalev, O. V. (1993). Representations of the Crystallographic Space Groups: Irreducible Representations, Induced Representations, and Corepresentations, 2nd ed. New York: Gordon and Breach]) notation in the parentheses. For the single propagation vector 1k structure, there are four subgroups shown in Fig. 1[link]. The one-dimensional irrep [{\rm mL}_{2}^{+}] results in the magnetic structure [R_{I}\bar{3}c] with the moments parallel to k1. The subgroups generated by the irrep [{\rm mL}_{3}^{+}] as the primary order parameter result in the magnetic structure with the moments perpendicular to k1 in MSG Cc2/m, whereas two other groups Cc2/c and [P_{S}\bar{1}] allow the presence of irrep [{\rm mL}_{2}^{+}] as a secondary order parameter, which allows a spin component along the propagation vector.

[Figure 1]
Figure 1
Group–subgroup graph for 4k structures based on the propagation vector star generated by [½ ½ ½] in space group [{Fm\bar{3}m}1^{\prime}] with the magnetic atom in the 4a (0,0,0) position. The possible 1k structures are indicated by label 1k. Only cubic, rhombohedral and monoclinic groups are shown, whose crystal metrics are compatible with the experiment.

We note that for the 1k-case, Cc2/m and [R_{I}\bar{3}c] are maximal magnetic subgroups of the parent grey Shubnikov group [{Fm\bar{3}m}{1^{\prime}}]. The basis transformation from cubic to monoclinic cells for both Cc2/m and Cc2/c reads: (½, ½, −1), (½, −½, 0), (−1, −1, 0). In MSG Cc2/c, Ni and O are in the 4c (000 | mx,0,mz) and 4b [(0\,{1\over 2}\,{1\over 4})] positions. In MSG Cc2/m, Ni and O are in the 4a (000 | 0,my,0) and 4d (0 ½ ¼) positions, respectively. The (ab) monoclinic plane is perpendicular to the propagation vector k1, so in Cc2/m the allowed spin direction is only perpendicular to k1. The propagation vector k1 and the a and c axes lie in the same plane with an ideal pseudocubic monoclinic angle (between a and c) β = arccos[−1/(3)1/2] ≃ 125.264°, and an angle between k1 and the c axis 270 − β ≃ 144.736°.

For the multi-k structures of the propagation vector star {[½ ½ ½], [−½ ½ ½], [½−½ ½], [½ ½ −½]} of [{Fm\bar{3}m}], the number of possible subgroups is significantly larger. In the primitive basis, the arms of the star are {[½ ½ ½], [0 0 ½], [0 ½ 0], [½ 0 0]}. As we will show below from the diffraction experiments, the crystal metric in the AFM state is rhombohedral. For that reason we exclude possible tetragonal and orthorhombic MSGs and their subgroups because the crystal metric from experiment cannot be fitted to the orthorhombic one. The eight MSGs that are left to consider are shown as a tree graph in Fig. 1[link]. Four subgroups labelled with 1k are MSGs allowed for a single propagation vector k1. We note that although the labels for the groups for 1k, 2k, 3k and 4k cases are the same, the magnetic structures are different. Three rhombohedral MSGs are [R_{I}\bar{3}m], ([{\rm mL}_{3}^{+}], P22), [R_{I}\bar{3}c] ([{\rm mL}_{3}^{+}], P21), and [R_{I}\bar{3}] ([{\rm mL}_{3}^{+}], C56), with the active primary irrep and direction of the order parameter shown in the parenthesis in accordance with Campbell et al. (2006[Campbell, B. J., Stokes, H. T., Tanner, D. E. & Hatch, D. M. (2006). J. Appl. Cryst. 39, 607-614.]). The basis transformation is (−1, 0, −1), (0, 1, 1), (2, 2, −2), and the single 4a site in [{Fm\bar{3}m}] is split into 18e and 6a positions with spin components 18e (mx,0,0), 6a (0,0,0); 18e (mx,2mx,mz), 6a (0,0,mz) and 18e (mx,my,mz), 6a (0,0,mz) for the above-listed rhombohedral groups, respectively. The last two MSGs allow the presence of the secondary irrep [{\rm mL}_{2}^{+}]. The first MSG does not allow the non-zero magnetic moment on all magnetic atom sites and hence can be excluded. Interestingly, the MSG [R_{I}\bar{3}c] generated by the primary irrep [{\rm mL}_{3}^{+}] with three active k vectors can be also generated by the irrep [{\rm mL}_{2}^{+}], C7 with four k vectors. The magnetic structures in both cases are the same. The MSG [R_{I}\bar{3}] is a direct subgroup of [R_{I}\bar{3}c] with identity basis transformation and has the same positions for Mn/Ni and O atoms, but due to the lost of the mirror plane, the point group allows a general direction of the magnetic moment in the 18e position.

The full star monoclinic subgroup Cc2/c ([{\rm mL}_{2}^{+}], S2) is constructed in a similar way to the 1k Cc2/c, but the single 4c position is split into three special positions 4c (mx,0,mz), 4d (mx,0,mz) and 8f (mx,my,mz). If the spin components are constrained to be the same for all three positions, then the magnetic 4k structure is equivalent to 1k structure.

4. Neutron diffraction results

4.1. Crystal structure

The diffraction patterns were collected at a base temperature of T = 2 K for all samples and at some selected temperatures below and above the Néel temperature TN = 115 K for MnO. The diffraction pattern for NiO taken at room temperature (∼300 K), is similar to the one at base temperature, because TN = 525 K is significantly higher than room temperature. Fig. 2[link] shows the diffraction patterns of MnO and NiO at 2 K together with the refinement curves. One can see that the rhombohedral distortion in MnO is very large: the group of five peaks around 2θ = 140° is mainly the splitting of the single cubic peak (711)/(551). In space group [{R\bar{3}m}], they are (321), (045), (137), (309) and (2,0,11). For NiO the rhombohedral splitting is about ten times smaller than for MnO, but is clearly seen as an [{R\bar{3}m}] doublet (404), (0,0,12) from the splitting of the (444) cubic peak, around 2θ = 147°. The rhombohedral distortions from the parent cubic metrics can be defined as the ratio δR = [(6)1/2c/a]/(6)1/2 and amount to δR ≃ 1.9 × 10−3 and 16 × 10−3 for NiO and MnO, respectively, where a and c are rhombohedral lattice constants.

[Figure 2]
Figure 2
The Rietveld refinement pattern and the difference plot of the neutron diffraction pattern for MnO (top) and NiO (bottom) at T = 2 K measured on the HRPT with the wavelength λ = 1.1545 Å. The line is the refinement pattern based on the magnetic model Cc2/c shown in Table 1[link]. The row of ticks show the Bragg peak positions. The insets show fragments of diffraction pattern at high angles.

The structure can be further refined in Cc2/c, leading to additional monoclinic splitting of the above peaks. The fit to the monoclinic group is converged with the metrics shown in Table 1[link]. We note that all the atoms preserve their special positions in the Cc2/c symmetry, so there are no extra structure parameters to refine, allowing a convergent fit. The monoclinic distortion can be refined as well. Despite the fact that the structures can be refined in Cc2/c, there is no apparent peak splitting and the improvement of the goodness-of-fit χ2 is marginal. The refined monoclinic distortion from rhombo­hedral metric corresponds to the splitting of the Bragg peaks (shown in the inset of Fig. 2[link]) δd/d ≃ 3 × 10−4 and 6 × 10−4 for the NiO doublet around 2θ = 147° and MnO around 2θ = 139.2°, respectively. The metric transformation from [{R\bar{3}m}] to Cc2/c reads: (−1, −2, 0), (−1, 0, 0), ([{2\over 3}, {4\over 3},-{2\over 3}]). This transformation results in the following monoclinic metric that preserves rhombohedral symmetry: A = (3)1/2a, B = b, [C = {2\over 3}\left( {3a^{2}+c^{3}} \right)^{1/2}], β = 180 − [\arccos[a/\left( {a^{2}+c^{2}/3} \right)^{1/2}]], where capital letters stand for the monoclinic metric. This `ideal' value of β, derived by the above formula, is also listed in Table 1[link]. We think that the refined monoclinic deviations from the rhombohedral symmetry may suffer from systematic errors related to profile parameters. From the viewpoint of high-resolution neutron diffraction, the crystal metric is indistinguishable from the rhombohedral one.

4.2. One-k1 magnetic structure

All magnetic Bragg peaks can be indexed with a single propagation vector k1 in [{Fm\bar{3}m}] settings. More precisely, due to very strong rhombohedral distortions, the propagation vector that fits the experimental peak positions is [{[{0}\,{0}\,{3\over 2}]}] in [{R\bar{3}m}] settings with the basis transformation from the cubic group shown in Table 1[link]. The one-k rhombohedral magnetic subgroup [R_{I}\bar{3}c] has the magnetic ion in the 6b Wyckoff position, which allows a spin component only parallel to [{[{0}\,{0}\,{3\over 2}]}], and gives the structure with ferromagnetic (ab) layers with antiparallel orientation of neighbouring layers along c. This type of structure does not fit the experiment because it does not explain the first magnetic peak (½ ½ ½) (in cubic setting) or (0 0 3) in the above rhombohedral group.

The data are very well fitted in the Cc2/c MSG (Fig. 3[link]) with the structure parameters shown in Table 1[link]. The directions of moments along specific monoclinic axes are dictated only by the symmetry of the magnetic space group and are conditional because the crystal metric does not show any explicit deviations from the rhombohedral one, according to our experimental data. There was no convergence of the fit with the mz component released for MnO, so it was fixed to mz = 0. Since the experimental data do not show monoclinic distortions, the absence of a spin component along the k1 vector raises the symmetry of MSG to Cc/2m, where the allowed spin direction is only my, which is perpendicular to k1, similar to mx. The fit in MSG Cc/2m gives the same goodness-of-fit and the same magnetic moment for MnO. The description of the magnetic structure for both samples is given in Cc/2c in Table 1[link]. The formal description of the magnetic structures in the mcif format is given in supporting information. Figs. 3[link] and 4[link] show the neutron diffraction patterns together with the fits to the magnetic structure models. For MnO, the figure shows the diffraction patterns below and above the Néel temperature TN = 115 K, where one can observe the progressive development of rhombohedral distortions as the temperature decreases to 100 K and 2 K. This effect is particularly noticeable at diffraction angles around 2θ ≃ 90° and 140°, evidenced by an increased splitting of the single cubic peaks that were observed above the Néel temperature at 125 K. The presence of the short-range magnetic correlations just above TN is also clearly seen as a broad (½ ½ ½) peak at 2θ ≃ 20°. The magnetic structures are illustrated in Fig. 5[link]. The magnetic moment sizes are close to the expected spin-only values 2 μB and 5 μB for NiO and MnO, respectively.

[Figure 3]
Figure 3
The Rietveld refinement pattern and the difference plot of the neutron diffraction pattern for NiO at T = 2 K measured at HRPT with the wavelength λ = 1.886 Å. The lines are the refinement patterns based on the magnetic model [{R_{I}\bar{3}c}] (top) and Cc2/c (bottom) shown in Table 1[link]. The magnetic contributions are shown by the green lines. The rows of ticks show the Bragg peak positions for crystal (upper) and magnetic (lower) structures. The difference between the observed and calculated intensities is shown by the dotted blue line. Since both refinements are visually identical, we zoom the bottom Cc2/c figure to better show the magnetic contribution. The minor unindexed peaks are due to λ/2, as discussed in the text and illustrated in Fig. 8[link].
[Figure 4]
Figure 4
The Rietveld refinement pattern and the difference plot of the neutron diffraction pattern for MnO at T = 2 K (top), 100 K (middle) and 125 K (bottom) measured at HRPT with the wavelength λ = 1.886 Å. The lines are the refinement patterns for MSG Cc2/c below the Néel temperature for T = 100 K and paramagnetic SG [{Fm\bar{3}m}] for 125 K. The magnetic contributions are shown by the green lines. The rows of ticks show the Bragg peak positions for crystal (upper) and magnetic (lower) structures. The difference between the observed and calculated intensities is shown by the dotted blue line.
[Figure 5]
Figure 5
Magnetic structure of NiO (left) and MnO (right) in MSG Cc2/c. The details of the magnetic models are given in Table 1[link].

The structure description in [{R\bar{3}m}] given in Table 1[link] does not correspond to any rhombohedral MSG but is equivalent to the magnetic structure with MSG Cc2/c, with the crystal metric fixed as explained in Section 4.1[link]. In this description, the crystal structure is refined in space group [{R\bar{3}m}], whereas the magnetic structure is defined in P1 with rhombohedral centring translations (R1). In the rhombohedral metric, the moment components are along and perpendicular to the c axis. For NiO, there is a spin component parallel to the propagation vector, whereas it is absent for MnO. This type of representation is useful for comparison with the literature data and the multi-k model considered in Section 4.3[link].

4.3. Multi-k magnetic structures in NiO

Two rhombohedral groups [R_{I}\bar{3}c] [No. 167.108 [R\bar{3}c.1^{\prime}_{c} \,(R\bar{3}m)] (UNI symbol)] and its subgroup [R_{I}\bar{3}] (No. 148.20 [R\bar{3}]) allow a spin component on both magnetic sites 18e and 6a. We performed a simulated annealing (SA) search (Kirkpatrick et al., 1983[Kirkpatrick, S., Gelatt, C. D. & Vecchi, M. P. (1983). Science, 220, 671-680.]; Rodríguez-Carvajal, 1993[Rodríguez-Carvajal, J. (1993). Physica B, 192, 55-69.]) of the full diffraction profile for possible solutions in most general rhombohedral MSG [R_{I}\bar{3}] for both NiO and MnO samples. The SA search starts from random values of the free parameters, and we have repeated the search more than several hundred times, so we are confident that we have not lost any solutions.

For the MnO, there is no reasonable multi-k rhombohedral solution in [R_{I}\bar{3}] that would be compatible with the 1k structure based on the goodness-of-fit. For NiO, we have found two structures that give a goodness-of-fit similar to the one-k structure in the Cc2/c 1k model. The results of the SA search are shown in Table 1[link]. One of the structures has `good' sizes of the magnetic moments, about 2 μB, which is shown on the first line for the [R_{I}\bar{3}] model SA search in Table 1[link]. One can see that the components along x and y are related by a factor of two, suggesting the [R_{I}\bar{3}c] symmetry. Using these initial values of the parameters, we have performed a standard least squares (LSQ) fit, but the fit did not converge for all parameters released. We attempted a restricted fit under the assumption of the same moment sizes on both Ni sites using spherical coordinates with a fixed ratio of mx/my = 0.5, but there was no convergence for the spherical θ angle to the z axis. So, in the final LSQ fit, the angle was fixed to −35° (close to the value from the SA search). The error bars correspond to the θ variation between −30° and −40°. The formal description of the structure is given in supporting information. The illustration of the fit quality can be seen in Fig. 3[link]. The magnetic structure in this multi-k model in MSG [R_{I}\bar{3}c] is shown in Fig. 6[link]. The magnetic structure description with respect to the atomic positions and the crystal metrics is the same in both MSG [R_{I}\bar{3}c] and its subgroup [R_{I}\bar{3}], provided that we correlate my = 2mx in the latter. In supporting information, we give the description of the structure in [R_{I}\bar{3}] and, for completeness, we also provide an mcif file with the equivalent description of this structure in MSG [R_{I}\bar{3}c].

[Figure 6]
Figure 6
Magnetic structure of NiO multi-k MSG [R_{I}\bar{3}c]. Details of the magnetic model are given in Table 1[link].

5. Discussion

The magnetic structure in MnO can be described using two different one-k MSGs, Cc2/c and Cc2/m, which give different magnetic moment directions with respect to the same monoclinic metrics in these groups. Since our diffraction data do not show significant deviations from the rhombohedral crystal metrics, we cannot unambiguously assign the MSG symmetry. We note that the above MSGs are not group–subgroup related. Therefore, we use the less restricted monoclinic MSG Cc2/c for the magnetic structure description for both samples. Furthermore, as we have shown, our data for both MnO and NiO samples are not sensitive to the orientation of the magnetic moment in the plane perpendicular to the rhombohedral c axis in [R\bar{3}m] SG (direction of k1 = [½ ½ ½] in parent SG [{Fm\bar{3}m}]) – we can only determine the components of the magnetic moment along and perpendicular to k1.

The neutron diffraction patterns allow two very different magnetic structure models for NiO, which fit the data equally well. The directions of the magnetic moments in the monoclinic one-k structure Cc2/c and the rhombohedral 4k structure [R_{I}\bar{3}c] are practically orthogonal to each other. In the one-k structure, the Ni moments form ferromagnetic layers with the moment direction approximately perpendicular to the propagation vector k1 and stacked antiferromagnetically (AFM) along k1. In the multi-k magnetic structure, the moments of the neighbouring Ni atoms are approximately AFM coupled. These two structures are completely incompatible with each other and we think that the 4k structure [R_{I}\bar{3}c] is a false solution and is only possible because of the very small rhombohedral distortions in NiO, as well as the small moment value. We cannot formally disregard this magnetic structure because it gives the same goodness-of-fit, but we like to show the reasons for excluding this model. Fig. 7[link] shows a comparison between the models. The lines show the difference between calculated diffraction profiles for the models [{R_{I}\bar{3}c}] and Cc2/c. The allowed magnetic Bragg reflections are localized at approximately the same 2θ values, but due to rhombohedral distortion, all peaks are split in the larger unit cell in the [{R_{I}\bar{3}c}] model. For instance, the first magnetic peak at 22.5°, which is a single (001) peak in Cc2/c, becomes a doublet in [{R_{I}\bar{3}c}] (101)/(003) with (101) at slightly lower 2θ. Thus, the (003) peak has zero intensity in [{R_{I}\bar{3}c}], but all calculated intensity belongs to the (101) peak. Conversely, the single (001) peak in Cc2/c, metrically corresponding to (003), describes all the experimental intensity. This repartitioning of the peak intensities results in the up-down shape of the calculated intensity difference clearly seen at all magnetic peak positions in Fig. 7[link]. Although the models in Figs. 5[link] and 6[link] may be indistinguishable by neutron diffraction, they almost certainly have very different total energies. Ab initio calculations can and should be used to select the correct model in cases like NiO.

[Figure 7]
Figure 7
The difference between the calculated neutron diffraction patterns shown in Fig. 3[link] for the models [{R_{I}\bar{3}c}] and Cc2/c with refinement parameters shown in Table 1[link] for NiO at T = 2 K. The line around zero level is the overall difference between the calculated patterns. The line shifted up to the count level 2500 is the difference of pure magnetic contributions between the two models. The rows of ticks show the Bragg peak positions for crystal (upper) and magnetic (lower) structures.

The presence of the extra Bragg peaks in larger unit cell of [{R_{I}\bar{3}c}] could potentially confirm the multi-k magnetic structure, although these peaks have only nuclear contributions in [{R_{I}\bar{3}c}]. Close inspection of the diffraction pattern shows that a few peaks are indeed present (Fig. 8[link]), but they originate from a very small admixture of the λ/2 component in the neutron beam and their intensities are ideally refined in the one-k Cc2/c model. Ironically, we did not realize the presence of 0.05% of λ/2 in the neutron beam, which theoretically should not be present for the (511) reflection of the Ge monochromator used for λ = 1.886 Å. Thus, the presence of these extra peaks was initially attributed to the multi-k structure, and we attempted to find the solutions in all possible multi-k structures shown in the subgroup tree (Fig. 1[link]).

[Figure 8]
Figure 8
A portion of the Rietveld refinement pattern for NiO at T = 2 K measured on the HRPT with wavelength λ = 1.886 Å, demonstrating the presence of a very small but visible contribution (about 0.05%) of λ/2 at 2θ about 26° and 37° corresponding the strongest Bragg peaks (111) and (311) in cubic metric. The black line is the refinement pattern Cc2/c (bottom) shown in Table 1[link]. The magnetic contribution is shown by the green line. The row of ticks shows the Bragg peak positions.

In MnO, due to significantly larger rhombohedral distortions, one can clearly see that the intensity distribution over the doublets like (101)/(003) definitively allows us to reject the rhombohedral model. We believe that the exchange inter­actions in both compounds are similar, and also on these grounds, the monoclinic model is the only acceptable solution for the magnetic structure in NiO. An interesting previous study (Goodwin et al., 2006[Goodwin, A. L., Tucker, M. G., Dove, M. T. & Keen, D. A. (2006). Phys. Rev. Lett. 96, 047209. ]) reports that the monoclinic, collinear structure of MnO is favored when considering the total (Bragg + diffuse) neutron scattering, and it also resolves the direction of the moment in the (111) planes. The solution proposed is similar to the one we have found, but with a different monoclinic symmetry, C2. This group cannot be a magnetic subgroup of the parent grey MSG, but it can, of course, be modelled using the lowest MSG [P_{S}{\bar 1}].

We would like to note that distinguishing between multi-k and one-k structures from diffraction experiments is difficult, if not impossible, unless there is a magnetostructural distortion associated with the different symmetries of the multi-k versus one-k cases. The point is that different k vectors from the propagation vector star do not interfere with each other, and their partial contributions are simply added to the Bragg peaks located at the same Q-value in powder diffraction. The one-k structure contributes to the Bragg peaks located at the same Q, as in the multi-k case, and we think that in most cases one can find one-k and multi-k models that provide the same integral intensity at the same Q. However, this is not always possible if we keep the same magnetic space group symmetry. For instance, in our case, the 1k and 4k magnetic structures with very different arrangements of the magnetic moments in the same MSG [R{\bar 3}c] (mL2+) result in different magnetic Bragg peak intensities, and the 1k structure can be definitely rejected. However, it is true that one can find the one-k structure in monoclinic MSG Cc2/c, which gives the same powder diffraction intensities in the absence of crystal metric distortions. The situation is not much easier in single-crystal diffraction where, instead of powder averaging, we always have magnetic domains for the one-k structure that mimic the multi-k case.

In some cases, one can distinguish between the models using single-crystal neutron diffraction in an applied magnetic field because multi-k and one-k structures will be modified differently [see, for instance, Puphal et al. (2020[Puphal, P., Pomjakushin, V., Kanazawa, N., Ukleev, V., Gawryluk, D. J., Ma, J., Naamneh, M., Plumb, N. C., Keller, L., Cubitt, R., Pomjakushina, E. & White, J. S. (2020). Phys. Rev. Lett. 124, 017202.])]. Alternatively, one could potentially affect the population of the magnetic domains by applying uniaxial pressure, thus identifying the one-k model. Additionally, there is an example of inelastic neutron scattering on powder samples where the authors claim that they can distinguish multi-q from single-q states (Paddison et al., 2024[Paddison, J. A. M., Zhang, H., Yan, J., Cliffe, M. J., McGuire, M. A., Do, S., Gao, S., Stone, M. B., Dahlbom, D., Barros, K., Batista, C. D. & Christianson, A. D. (2024). npj Quantum Mater. 9, 48.]) from the analysis of the magnetic excitations in linear spin-wave theory.

6. Summary

We have performed experiments using high-resolution powder neutron diffraction on MnO and NiO monoxides with Néel temperatures (TN) of 115 K and 525 K, respectively, in order to refine their crystal and magnetic structures. The structures below TN are generated by a propagation vector k1 = [½ ½ ½] in the parent paramagnetic cubic SG [{Fm\bar{3}m}]. The rhombohedral distortions from the paramagnetic parent cubic metrics are clearly seen in our diffraction data and were refined in SG [R\bar{3}m], amounting to 1.9 × 10−3 and 16 × 10−3 for NiO and MnO, respectively. In accordance with the rhombohedral symmetry, there are eight possible magnetic space subgroups for the one-k and multi-k structures. Among them, there are the three most symmetric MSGs: [{R_{I}\bar{3}c}], Cc2/m and Cc2/c. The refined monoclinic distortions from the rhombohedral metric are very small, about a few 10−4, and can be disregarded, implying that the orientation of the magnetic moment in the plane perpendicular to the rhombohedral c axis cannot be inferred from our data.

We have found that the best antiferromagnetic structure, which fits both MnO and NiO, is the 1k structure with MSG Cc2/c (BNS No. 15.90, UNI symbol C2/c.1′c [C2/m]). In this structure, the moments form ferromagnetic layers with the moment direction approximately perpendicular to the propagation vector k1 and stacked antiferromagnetically along k1 with magnetic moment sizes m = 1.94 μB and 4.61 μB for NiO and MnO, respectively. The same magnetic configuration can be described equally well in MSG Cc2/m for MnO. For NiO, the multi-k rhombohedral structure in MSG [{R_{I}\bar{3}c}], which is very different from the Cc2/c structure, can equally well fit our diffraction data with similar value of m. However, we conclude that this is a false solution and is only possible because of the very small rhombohedral distortions and the small magnetic moment value in NiO.

In conclusion, we would like to note that this is not an uncommon case where the crystal distortions due to magnetoelastic coupling below the magnetic transition are small or even undetectable in diffraction experiments. Since a paradigm such as `one propagation vector at the magnetic transition is enough' is still strong in the neutron scattering community, some high-symmetry multi-k magnetic structures are even not considered as possible structure solutions, which could, however, fit the experimental data equally well.

Supporting information


Acknowledgements

The work was performed at the Swiss neutron spallation source SINQ at Paul Scherrer Institute.

References

First citationAroyo, M., Perez-Mato, J., Orobengoa, D., Tasci, E., De La Flor, G. & Kirov, A. (2011). Bulg. Chem. Commun. 43, 183.  Google Scholar
First citationCampbell, B. J., Stokes, H. T., Tanner, D. E. & Hatch, D. M. (2006). J. Appl. Cryst. 39, 607–614.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationCheetham, A. K. & Hope, D. A. O. (1983). Phys. Rev. B, 27, 6964–6967.  CrossRef CAS Web of Science Google Scholar
First citationFischer, P., Frey, G., Koch, M., Könnecke, M., Pomjakushin, V., Schefer, J., Thut, R., Schlumpf, N., Bürge, R., Greuter, U., Bondt, S. & Berruyer, E. (2000). Physica B, 276–278, 146–147.  Web of Science CrossRef CAS Google Scholar
First citationGoodwin, A. L., Tucker, M. G., Dove, M. T. & Keen, D. A. (2006). Phys. Rev. Lett. 96, 047209.   Google Scholar
First citationKirkpatrick, S., Gelatt, C. D. & Vecchi, M. P. (1983). Science, 220, 671–680.  CrossRef PubMed CAS Web of Science Google Scholar
First citationKovalev, O. V. (1993). Representations of the Crystallographic Space Groups: Irreducible Representations, Induced Representations, and Corepresentations, 2nd ed. New York: Gordon and Breach  Google Scholar
First citationLee, S., Ishikawa, Y., Miao, P., Torii, S., Ishigaki, T. & Kamiyama, T. (2016). Phys. Rev. B, 93, 064429.  Google Scholar
First citationLim, J. S., Saldana-Greco, D. & Rappe, A. M. (2016). Phys. Rev. B, 94, 165151.  Google Scholar
First citationMomma, K. & Izumi, F. (2011). J. Appl. Cryst. 44, 1272–1276.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationPaddison, J. A. M., Zhang, H., Yan, J., Cliffe, M. J., McGuire, M. A., Do, S., Gao, S., Stone, M. B., Dahlbom, D., Barros, K., Batista, C. D. & Christianson, A. D. (2024). npj Quantum Mater. 9, 48.  Google Scholar
First citationPuphal, P., Pomjakushin, V., Kanazawa, N., Ukleev, V., Gawryluk, D. J., Ma, J., Naamneh, M., Plumb, N. C., Keller, L., Cubitt, R., Pomjakushina, E. & White, J. S. (2020). Phys. Rev. Lett. 124, 017202.  Web of Science CrossRef PubMed Google Scholar
First citationRodríguez-Carvajal, J. (1993). Physica B, 192, 55–69.  CrossRef Web of Science Google Scholar
First citationSchrön, A., Rödl, C. & Bechstedt, F. (2012). Phys. Rev. B, 86, 115134.  Google Scholar
First citationShull, C. G., Strauser, W. A. & Wollan, E. O. (1951). Phys. Rev. 83, 333–345.  CrossRef CAS Web of Science Google Scholar
First citationStokes, H. T. & Hatch, D. M. (1988). Isotropy Subgroups of the 230 Crystallographic Space Groups. Singapore: World Scientific.  Google Scholar

This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.

Journal logoSTRUCTURAL SCIENCE
CRYSTAL ENGINEERING
MATERIALS
ISSN: 2052-5206
Follow Acta Cryst. B
Sign up for e-alerts
Follow Acta Cryst. on Twitter
Follow us on facebook
Sign up for RSS feeds