Multipole expansion for magnetic structures: A generation scheme for symmetry-adapted orthonormal basis set in crystallographic point group
M. -T. Suzuki T. Nomoto R. Arita Y. Yanagi S. Hayami H. Kusunose



Center for Computational Materials Science, Institute for Materials Research, Tohoku University, Sendai, Miyagi 980-8577, Japan

Department of Applied Physics, University of Tokyo 7-3-1 Hongo Bunkyo-ku, Tokyo 113-8656, Japan

Department of Applied Physics, University of Tokyo 7-3-1 Hongo Bunkyo-ku, Tokyo 113-8656, Japan
RIKEN Center for Emergent Matter Science (CEMS), Wako, Saitama 351-0198, Japan

Center for Computational Materials Science, Institute for Materials Research, Tohoku University, Sendai, Miyagi 980-8577, Japan

Department of Physics, Hokkaido University, Sapporo 060-0810, Japan

Department of Physics, Meiji University, Kawasaki 214-8571, Japan



We propose a systematic method to generate a complete orthonormal basis set of multipole expansion for magnetic structures in arbitrary crystal structure.
The key idea is the introduction of a virtual atomic cluster of a target crystal, on which we can clearly define the magnetic configurations corresponding to symmetry-adapted multipole moments.
The magnetic configurations are then mapped onto the crystal so as to preserve the magnetic point group of the multipole moments, leading to the magnetic structures classified according to the irreducible representations of crystallographic point group.
We apply the present scheme to pyrhochlore and hexagonal ABO_3 crystal structures, and demonstrate that the multipole expansion is useful to investigate the macroscopic responses of antiferromagnets.


Multipole expansion for magnetic structures: A generation scheme for symmetry-adapted orthonormal basis set in crystallographic point group
    H. Kusunose
    December 30, 2023
===========================================================================================================================================





§ INTRODUCTION

Diversity of physical properties of magnets provides a fascinating playground in condensed matter physics.
When we explore this exciting arena, what is interesting to note is that there are many restrictions imposed by the symmetry of the magnetic structures.
For example, it has been well known that the structure of linear response tensors are determined by the magnetic point group<cit.>.
Furthermore, identifying the order parameter for the magnetic phase is most useful for deeper understanding of physical phenomena.
There has been many studies to investigate the relation between the
order parameters in a particular magnetic structure and macroscopic
phenomena such as anomalous Hall (AH)
effect <cit.>, electromagnetic (EM)
effect <cit.>, and optical responses <cit.>.
Thus we have a significant chance to specify or even design a magnet exhibiting desired physical properties by investigating the order parameters which characterize the magnetic structures.

For this purpose, the multipole expansion of the magnetic structure is an efficient and powerful approach.
Indeed, the multipole moments inherent in the magnetic structure, which we call cluster multipole moments <cit.>, have played a crucial role as a key order parameter in a variety of studies: EM effect has been discussed in terms of the magnetic (M) rank-0 monopole and rank-2 quadrupole as well as the magnetic toroidal (MT) rank-1 dipole as the order parameter to characterize the specific magnetic structures <cit.>.
The relation between parity odd multipoles and electromagnetic effect
has recently been investigated based on generalized forms of the multipole expansions for magnetic distributions <cit.>.
It has also been shown that the rank-3 M multipole (octupole) plays a key role for a large AH effect <cit.>, anomalous Nernst effect <cit.>, and magneto-optical Kerr effect <cit.> in the coplanar antiferromagnets Mn_3Z (Z=Sn, Ge) <cit.>.
In these studies, the interplay between the physical properties and the magnetic structure through cluster multipoles has been investigated extensively.
However, there has been no concrete scheme to make a complete basis set of cluster multipoles for a given crystal structure.

In this paper, we propose a scheme to generate cluster multipoles which form a complete orthonormal basis set for arbitrary magnetic structures.
Here, we introduce a virtual atomic cluster, which depends only on the crystallographic point group of the system.
We define unambiguously the magnetic configurations corresponding to the symmetry-adapted multipoles in the atomic cluster.
The obtained magnetic configurations are mapped to the original crystal structure with the magnetic point group symmetry preserved.
The generated complete basis set for the magnetic structure in crystal is useful to measure the symmetry breaking as an order parameter according to the magnetic point group <cit.>.
Although we restrict ourselves to the case with “uniform” magnetic structures characterized by the ordering vector q=0 in this paper, an extension to cases with nonzero q ordering vectors is straightforward.

To demonstrate the efficiency of the present scheme, we apply the cluster multipole expansion of magnetic structures to pyrochlore and hexagonal ABO_3 crystal structures.
For the pyrochlore structure, it is shown that the all-in all-out magnetic structure corresponds to a 
M octupole, and two-in two-out and one-in three-out magnetic structures are expressed by the linear combinations of the M dipole and octupole that 
belong to the same irreducible representation (IREP) of the crystallographic point group.
The two-in two-out and one-in three-out magnetic structures can be transformed continuously to the pure antiferromagnetic structures, indicating that the antiferromagnetic structures without net magnetization yield the AH effect.
For the hexagonal ABO_3 structure, higher rank multipoles such as MT quadrupole and M octupole are necessary to describe the magnetic structures exhibiting the AH and EM effects within the uniform magnetic ordering.










§ MULTIPOLE EXPANSION OF MAGNETIC STRUCTURES IN CRYSTAL



 §.§ Outline of multipole expansion of magnetic structures in crystal

In this paper, we present a method to generate an orthonormal complete basis set of magnetic structures corresponding to the cluster multipoles classified according to the IREPs of the crystallographic point group.
To make free from confusion and ambiguity, we define some words used in this paper.
“Atomic cluster” is used for the atoms transformed to each other only by the rotation operations of a point group and distinguished from “crystal” which assumes a periodicity for an atomic configuration.
“Magnetic configuration” is used for the alignment of magnetic dipole moments on atoms of the atomic cluster and distinguished from “magnetic structure” which is for the alignment of magnetic dipole moments on atoms in periodic crystal.
A magnetic configuration (a magnetic structure) characterized by the M and/or MT multipoles is called multipole magnetic configuration (multipole magnetic structure).
The outline of the generation procedure for the orthonormal multipole
magnetic structures in crystal is following (see also Fig. <ref>):



  * Set a virtual atomic cluster corresponding to the crystallographic point group of a target crystal.

  * Generate magnetic configurations corresponding to the symmetrized M and MT multipoles in the virtual atomic cluster.

  * Map the magnetic dipole moments on the atoms in the virtual atomic cluster to those on the crystallographically equivalent atoms in crystal, which represent the multipole magnetic structures.

  * Orthonormalize the bases of multipole magnetic structures by using
      the Gram-Schmidt orthogonalization procedure.


Details for the procedure will be explained in Sec. <ref>, <ref>, and <ref>.
With this method, all the generated orthonormal magnetic structures are characterized by the symmetry-adapted multipoles.
Thus, the multipoles can be a useful measure of symmetry breaking of the crystallographic point group in the presence of the uniform magnetic orders.
It is a natural extension of the conventional dipole magnetization which measures the symmetry breaking of the ferromagnetic order.
The multipole expansion for magnetic structures in crystal is very efficient to investigate the relation between the magnetic structures and physical properties beyond the symmetry analysis.



 §.§ Multipole magnetic configurations in atomic clusters

The multipole expansion of the vector gauge potential is given under the Coulomb gauge ∇· A( r)=0 as follows:

    A( r)=∑_ℓ m(b_ℓM_ℓ m Y_ℓ m^ℓ(r̂)/r^ℓ+1 + c_ℓ T_ℓ m Y_ℓ m^ℓ+1(r̂)/r^ℓ+2),

where Y_ℓ m^ℓ'(r̂) (ℓ≥ 1, -ℓ≤ m ≤ℓ, ℓ'=ℓ-1, ℓ,ℓ+1) is the vector spherical harmonics that transforms as conventional scalar spherical harmonics Y_ℓ m(r̂) for rotation operation with its orbital angular momentum ℓ' <cit.>.
The expansion coefficients in Eq. (<ref>) 
are the so-called M multipole, M_ℓ m, and MT multipole, T_ℓ m, respectively (b_ℓ and c_ℓ are introduced for convenience).
The M and MT multipoles around a single magnetic ion in the unit of Bohr magneton are defined as

    M_ℓ m   =√(4π/2ℓ+1)∑_j(2ℓ_j/ℓ +1+σ_j)·O_lm(r_j),
    
    
    T_ℓ m   =√(4π/2ℓ+1)∑_j{ r_j/ℓ +1×(2ℓ_j/ℓ +2+σ_j)}·O_lm(r_j),

with

    O_lm(r)≡∇[ r^ℓY^*_ℓ m(r̂) ],

where ℓ_j and σ_j are the orbital and spin angular momentum of an electron at r_j.
Here, r=|r| and r̂=r/r.
Since the vector spherical harmonics Y_ℓ m^ℓ' are an orthonormal complete basis set of a vector function on a sphere, the M and MT multipoles M_ℓ m and T_ℓ m represent the arbitrary angular dependence of the magnetization distribution on a single magnetic ion.
Focusing on the spin part of Eqs. (<ref>) and (<ref>), the M and MT multipoles can be extended straightforwardly to characterize the classical magnetic configurations { m_i } in an atomic cluster whose atoms are transformed to each other through the point-group symmetry operations with respect to the symmetry center of the atomic cluster.
The explicit definitions are given as

    M_ℓ m   ≡√(4π/2ℓ+1)∑_i=1^N_ atom m_i·O_lm(R_i), 
    
    T_ℓ m   ≡1/ℓ+1√(4π/2ℓ+1)∑_i=1^N_ atom( R_i× m_i)·O_lm(R_i),

where R_i is the position vector of i-th atom and N_ atom is the number of atoms in the atomic cluster.
These multipoles have the same transformation property with that of the spherical harmonics Y_ℓ m for rotation operation of the point group, and hence they are classified according to the IREPs of the point group <cit.>.
The spherical harmonics are usually symmetrized according to the IREPs as

    𝒴_ℓγ(r̂) =∑_m c^γ_ℓ m Y_ℓ m(r̂),

where γ runs from 1 to 2ℓ+1 in order to distinguish the IREP and its component including multiplicity which is necessary when the same IREPs are multiply appeared in the same rank ℓ.
The coefficients c^γ_ℓ m of the symmetrized spherical harmonics 𝒴_ℓγ(r̂) are tabulated in Ref. Kusunose2008 for instance, where c^γ_ℓ m are chosen so that 𝒴_ℓγ(r̂) is real.
This is always possible in the presence of the time-reversal symmetry.
The symmetry-adapted multipoles are thus reexpressed as follows:

    M_ℓγ   =∑_i=1^N_ atom u_ℓγ i^M· m_i,
    
    
    T_ℓγ   =∑_i=1^N_ atom u_ℓγ i^T· m_i,

where

    u^M_ℓγ i   =√(4π/2ℓ+1)𝒪_lγ(R_i),
    
     u^T_ℓγ i   =1/ℓ+1√(4π/2ℓ+1)(𝒪_lγ(R_i)× R_i),

and 𝒪_lγ is given by replacing Y_lm^* in Eq. (<ref>) with 𝒴_lγ.

In generating an orthogonal basis set of the magnetic structures based on the symmetry, the conventional projection operator method requires nontrivial trials to find out the suitable trial functions, and become complicated when the crystal contains many magnetic atoms.
Moreover, that procedure has large ambiguity for low-symmetry crystals, and for the complicated magnetic structures, it would be difficult to find out the correspondence between the magnetic structures and symmetry-adapted multipoles.
On the contrary, our new method is highly advantageous since it is possible to automatically generate a complete orthonormal basis set of magnetic structures, and it is, by definition, based on the symmetry-adapted multipoles.

Here, we comment on the so-called M monopole or magnetic flux configuration that corresponds to a magnetic configuration in which all the magnetic dipole moments point to the center of the atomic cluster.

 The monopole magnetic configuration is parity odd and invariant
under all the rotational symmetry operations of the point group. 
Such a monopole magnetic configuration can be defined as <cit.>:

    u^M_01i≡ R_i/R_i^2.

In this paper, however, we do not include this type of M monopole in the cluster multipole expansion since the ordinary multipole expansion in Eq. (<ref>) does not contain the monopole term and the corresponding magnetic configuration always appears as a higher-rank multipole of the parity odd fully symmetric 
IREP, as will be shown in the case of hexagonal ABO_3 in Sec. <ref>.




 §.§ Virtual atomic clusters for crystallographic point groups





In crystal, crystallographically equivalent atoms are transformed to
each other by combinations of point group symmetry
operation, ℛ_i, non-primitive translation, τ_j,
and primitive translation, T_j, which constitute
symmetry operations of space group.
According to Neumann's principle, the appearance of macroscopic phenomena such as AH and EM effects are determined by crystallographic point group, whose symmetry operations consist only of the rotation part ℛ_i of the space group.
Therefore, in order to discuss the macroscopic phenomena induced by the antiferromagnetic order, it is necessary to define appropriate order parameters for antiferromagnetic order reflecting the symmetry breaking of the crystallographic point group.

The crystallographic point group corresponding to the space group of crystal is more definitely defined as follows.
The space group 𝒢 is decomposed as

    𝒢 = ∑_i=1^N_ coset{ℛ_i|τ_i}ℋT,

where the subgroup ℋ consists only of the pure rotation symmetry operations {h_1, h_2, …,h_N_0}, T is the group consisting of lattice translations T_j, and {ℛ_i|τ_i} are representative elements of 𝒢 with rotation operation ℛ_i (ℛ_1=E (identity operation)) and non-primitive translation operation τ_i (τ_1=0).
N_ coset is the number of cosets.
Note that the terms with i≥ 2 in Eq. (<ref>) exist only in nonsymmorphic space groups.
The crystallographic point group 𝒫 corresponding to the space group 𝒢 is then defined as

    𝒫 = ∑_i=1^N_ coset{ℛ_i| 0}ℋ.


Magnetic configurations in 
an atomic cluster is unambiguously defined by multipoles through Eqs. (<ref>) and (<ref>).
It is highly contrast to a direct generation of magnetic structures in crystal, which is not straightforward since the atoms in crystal are related not only by rotation operations but also by translation operations as we discussed in Ref. Suzuki2017.
To avoid difficulties in the direct generation of magnetic structures in crystal, our strategy to generate the symmetry-adapted magnetic structures in crystal is to generate the multipole magnetic configurations at first in a virtual atomic cluster defined under the point group 𝒫, and then map the magnetic configurations onto the magnetic structure in crystal so as to preserve the magnetic point group symmetry of the multipole configurations.
The virtual atomic cluster is defined as an atomic cluster consisting of the same number of atoms as the symmetry operations in the crystallographic point group 𝒫,
in which their atomic positions are given as the general Wyckoff positions of the corresponding symmorphic space group,

as listed in Table <ref> for representative point groups.
In the virtual atomic cluster, the magnetic configurations corresponding to the M and MT multipoles, Eqs. (<ref>) and (<ref>), are obtained by Eqs. (<ref>) and (<ref>), respectively, with respect to the origin of the virtual atomic cluster.
The examples of the virtual atomic clusters are depicted in Fig. <ref> for O_h, C_4h, D_6h, C_6v point groups.
Note that the virtual atomic cluster has ambiguity due to the choice of the parameter x, y, z, as seen in the Table <ref>.
This leads to an arbitrariness of the magnetic configurations in the virtual atomic cluster.
This arbitrariness is largely reduced when the magnetic configuration is mapped onto the atoms at high symmetry sites in crystal, as explained in Sec. <ref>.




 




 §.§ Mapping of multipoles from virtual atomic cluster to crystal



The purpose of this section is to obtain the magnetic structures whose
transformation property for the magnetic point group operations is
the same with that of multipole configurations in the virtual atomic cluster.
For this purpose, we first choose one atom in the virtual atomic cluster and one in crystal and set the same M dipole moments on these atoms by setting a mapping from the atoms in virtual atomic cluster to that in crystal.
The whole mapping between the atoms in the virtual atomic cluster and
the symmetrically equivalent atoms in crystal is obtained by identifying
the atoms transformed by the point group symmetry operations, R_i, in
the virtual atomic cluster with the atoms transformed by the symmetry
operations of space group, {E|T_j}{R_i|τ_i}h_k, in the crystal for the initially chosen atoms in both systems.
The mapping from the atoms in the virtual atomic cluster to the symmetrically equivalent atoms in crystal obviously does not have one-to-one correspondence.

The M dipole moment on an atom in crystal is obtained by summing up the M dipole moments on the all of the corresponding atoms in the virtual atomic cluster, as will be discussed in Sec. <ref>.
The generated magnetic structures in crystal are fully characterized by the symmetry-adapted multipoles in the virtual atomic cluster through the mapping, which we call the multipole magnetic structures.
We note that the correspondence between generated magnetic structures and multipole configurations depends on the choice of the first mapping between the atoms in virtual atomic cluster and that in crystal.
It occurs especially for low-symmetry crystallographic point groups having the multiple symmetrized bases within the same IREPs and rank.
In this case, we have to specify the representative atoms in the virtual atomic cluster and in crystal to identify the correspondence between the magnetic structures and multipole configurations.
The other arbitrariness for the multipole structures arises from a choice of the parameter (x, y, z) of virtual atomic cluster as listed in Table <ref>, as already mentioned in Sec. <ref>.

The advantage to introduce the virtual atomic cluster is that once we generate the multipole magnetic configurations in the virtual atomic cluster, we can systematically generate complete orthonormal basis sets for arbitrary crystal structures by appropriate mappings, in which they share the common 
point group between the atomic cluster and crystal.
This is in a high contrast with our previous scheme proposed in Ref. Suzuki2017 to identify the multipoles in crystal.
In the previous scheme, the atomic cluster is defined not in a virtual one but directly for the atoms in the crystal.
 We identified N_ coset clusters in which one consist of the atoms transformed to each other by the point group operations of ℋ in Eq. (<ref>) and the others are determined from transformation of {ℛ_i|τ_i}ℋ for the cluster in the crystal unit cell.
The macroscopic multipole moment was obtained by summing up the multipole moments in each cluster

in accordance with a nonsymmorphic space group of Mn_3Z (Z=Sn, Ge).
However, there are some cases that the higher rank multipoles are ill defined, especially when the magnetic atoms are located at high symmetry sites.
 The present scheme does not cause this problem and is advantageous in its efficiency since the arbitrariness of the generated orthonormal magnetic structures is largely reduced as compared with the conventional method using projection operator <cit.>, whose generated magnetic structures highly depend on the choice of the trial functions to be operated.
Based on the present procedure, we can systematically and automatically generate magnetic structures corresponding to the multipoles according to the IREPs of the crystallographic point group of the focusing crystal.




 §.§ Orthonormalization of basis set for magnetic structure in crystal



The multipole expansion in Eq. (<ref>) requires infinite number of components to express magnetization distribution in continuous space.
In contrast, the multipole expansion in magnetic structures that does not break the crystal periodicity is represented by a linear combination of 3N_ atom orthogonal basis set, where N_ atom is the number of atoms in a crystal unit cell.
In this subsection, we explain how to obtain the 3N_ atom orthonormal basis set which is sufficient to express a uniform magnetic order.

For notational convenience, we introduce the vector notation for uniform magnetic structures, { a}=( a_1, a_2,⋯, a_N_ atom), where a_i represents the three component vector on the i-th atom in the crystal unit cell.
Usually, the magnetic structure { m_i} is decomposed into the ferromagnetic part { m_i^ FM}, where m_i^ FM=∑_i m_i/N_ atom, and the rest antiferromagnetic part { m_i^ AFM}, where m_i^ AFM= m_i- m_i^ FM<cit.>.
Such a decomposition for the magnetic structure is generalized to obtain the 3N_ atom orthonormal complete basis set by using Gram-Schmidt orthonormalization procedure as follows:

    { e_1γ^1}   ≡{ u^1_1γ}/√(({ u^1_1γ}·{ u^1_1γ})),
    
    { v_ℓγ^μ}   ={ u_ℓγ^μ}-∑_ℓ'=1^ℓ∑_μ'=1^μ∑_γ'=1^γ-1({ u_ℓγ^μ'}·{ e_ℓ'γ'^μ'}){ e_ℓ'γ'^μ'},
    
    { e_ℓγ^μ}   ={ v_ℓγ^μ}/√(({ v_ℓγ^μ}·{ v_ℓγ^μ})),

where μ=1 and 2 represent the M and MT multipoles, respectively.
The initial u^1_1γ in Eq. (<ref>) is set as the M dipole moments, M_x, M_y, and M_z for γ=1,2,3, respectively.
The iterated calculation of the Gram-Schmidt orthonormalization procedure, Eqs. (<ref>) and (<ref>), starting from the lower-rank multipoles to higher ones as shown in Fig. <ref>, automatically generates orthonormal complete basis set of the uniform magnetic structures classified according to the IREPs of crystallographic point group.
After orthonormalization procedure, the magnetic structures corresponding to the higher-rank multipole may not be the pure multipole with definite rank as in Eq. (<ref>).
However, since the subtraction in Eq. (<ref>) is to eliminate the overlap between the highest-rank multipole and the lower ones, and { e_ℓγ^μ} in Eq. (<ref>) always contains { u_ℓγ^μ}, it can be regarded as the rank-ℓ multipole.

The generated finite norm 
of vectors { e_ℓγ^μ} are stored sequentially as { e^i} with the sequential indices i=1,...,3N_ atom.
With this procedure, the ferromagnetic structures { e^i} with i≤ 3 are orthogonal to the antiferromagnetic structures with i>3.
The obtained magnetic structure basis set are orthonormal, i.e., ({ e^i}·{ e^j})=δ_ij.
Since the obtained basis set is complete for the uniform magnetic structures, arbitrary uniform magnetic structure can be expressed as { m}=∑_i=1^3N_ atomc_i{ e^i} with c_i=({ m}·{ e^i}).
Note that the relation ∑_i^3N_ atom|c_i|^2=∑_i^N_ atom| m_i|^2 holds.




§ EXAMPLES OF MULTIPOLE EXPANSION FOR CRYSTALS



 §.§ Simple examples



We here illustrate the correspondence between the symmetry-adapted multipole magnetic configurations in the virtual atomic cluster and magnetic structures in crystal.
We take two simple examples, i.e., the atoms placed at 8l and 2e
Wyckoff sites in symmorphic space group P4/m (No. 83), and 8k and
2a Wyckoff sites in nonsymmorphic space group P4_2/m (No. 84).
The crystallographic point group of both space groups is C_4h, which corresponds to a virtual atomic cluster consisting of eight atoms transformed to each other through the symmetry operations 
as shown in Fig. <ref> (b).
In Fig. <ref>, two fold rotation C_2 transforms the atom 1 to atom 2, four fold rotation C_4 to atom 3, C_4^-1 to atom 4, space inversion I to atom 5, IC_2 to atom 6, IC_4 to atom  7, and IC_4^-1 to atom 8.
These transformation relations provide mapping from the atoms in the virtual atomic cluster to symmetrically equivalent atoms in the crystal as discussed in Sec. <ref>.
Figures <ref> and <ref> show the relation between the atoms in the virtual atomic cluster and the atoms at the Wyckoff sites in crystal.






For the C_4h point group, the rank-2 spherical harmonics are symmetrized as

    𝒴_21   = Y_20   (A_g), 𝒴_22   = 1/√(2)(Y_2-1-Y_21)   (E_g), 𝒴_23   = 1/√(2)i(Y_2-1+Y_21)   (E_g),
    
    𝒴_24   = 1/√(2)(Y_2-2+Y_22)   (B_g), 𝒴_25   = 1/√(2)i(Y_2-2-Y_22)   (B_g).

The multipole configurations are calculated from Eq. (<ref>) and (<ref>) on the virtual atomic cluster as shown in Fig. <ref>(b), where the atom 1 is placed at (0,0.3,0,3).
Figure <ref> shows the magnetic configuration generated from Eq. (<ref>) with the symmetrized spherical harmonics 𝒴_25 in Eq. (<ref>), corresponding to the MT quadrupole, T_xy, classified to B_g IREP of C_4h point group <cit.>.
This magnetic configuration in the virtual atomic cluster is mapped to symmetrically equivalent atoms, denoted by the signatures of Wyckoff positions, to obtain the corresponding magnetic structures in crystals as shown in Fig. <ref>.
For the general Wyckoff positions in each crystal (8l site of
P4/m and 8k site of P4_2/m), the magnetic structures are straightforwardly obtained through the mapping (Fig. <ref>).
Meanwhile, for 2e site of P4/m and 2a site of P4_2/m, multiple atoms in the virtual cluster are mapped to the same atoms in crystal, leading to the cancellation of the M dipole moment that is consistent with the symmetry in these Wyckoff sites.



 §.§ Pyrochlore structure








Pyrochlore crystal structure belongs to the space group Fd3̅m (No. 227), whose crystallographic point group O_h leads to the virtual atomic cluster as shown in Fig. <ref>(a).
The generation procedure of orthonormal magnetic structures gives twelve orthonormal bases characterized by the M dipole and octupole, and MT quadrupole as shown in Table <ref>.
In the pyrochlore structure, the inversion symmetry breaking is accompanied by the breaking of the translation symmetries since the two atoms related by the space inversion are also transformed by the commensurate translation of the crystal.
Therefore, a uniform magnetic structure 
must preserve the space inversion symmetry, restricting the parity-even multipole structures without EM effects.
The generated magnetic structure coincides with the basis set obtained in earlier work. <cit.>.





The pyrochlore compounds are known to show a variety of magnetic order depending on the atomic constitutions, such as all-in all-out, two-in two-out, and one-in three-out magnetic orders.
The relation between the observed magnetic structures and the M multipoles are shown in Table <ref>.
The corresponding uniform magnetic orders are shown in Fig. <ref>.
All-in all-out magnetic structures, reported in Cd_2Os_2O_7 and Er_2Ti_2O_7, correspond to M_xyz-octupole in Fig. <ref>.
The magnetic point group of the all-in all-out structure prohibits to induce AH effect.
Meanwhile, two-in two-out magnetic structure is expressed by a linear combination of the M dipole and octupole belonging to T_1g IREP as shown in Fig. <ref>.
The two-in two-out magnetic structure is allowed to induce the AH effect by the magnetic point group symmetry.
Since arbitrary linear combination does not change the symmetry, pure T_1g octupole structure M_z^α without net magnetization can also induce the AH effect.
This fact may be relevant to the mechanism of the AH effect of Pr_2Ir_2O_7 without net magnetization <cit.>.
The order parameter for the magnetic structures of Yb_2Ti_2O_7, Yb_2Sn_2O_7, and Tm_2Mn_2O_7 are also characterized by the T_1g multipoles, but in these cases, the M dipole moment is dominant <cit.>.









 §.§ Hexagonal ABO_3





Hexagonal ABO_3 compounds belong to P6_3cm (No. 185) space group, whose crystallographic point group is C_6v, and the corresponding virtual cluster is 
shown in Fig. <ref>(d).
The magnetic structure generation gives eighteen magnetic structures as shown in Fig. <ref>.
The present scheme generates the real basis set, while it contains the complex expressions for two dimensional IREPs, E_1 and E_2, in the previous work <cit.>.
In Table <ref>, we summarize the relation between the multipole magnetic structures and possible AH and EM effects induced under uniform magnetic order.
Table <ref> shows that the higher-rank multipoles are necessary to fully describe the EM effects.

The magnetic order characterized by M_u quadrupole is recognized in LuFeO_3, and in a temperature range of  75K-129K in ScMnO_3 <cit.>.
The M_u quadrupole induces the diagonal components in the EM tensor, α_xx=α_yy and α_zz.
The M quadrupole belongs to A_2 IREP, which is the same as that of M_z dipole in the C_6v crystallographic point group, and hence it can also induce the AH conductivity σ_xy as shown in Table <ref>.
We note that the M_u quadrupole structure is also regarded as the monopole magnetic structure characterized by M_0 since the same magnetic structure is obtained by using Eq. (<ref>) with the appropriate mapping from the virtual cluster to the crystal unit cell.

The magnetic structures observed in HoMnO_3 <cit.> and YbMnO_3 <cit.> are characterized by M_3a and M_3b octupoles, which do not induce the EM effect (Table <ref>).
The magnetic structures in YMnO_3 and the low-temperature phase in ScMnO_3 in Table <ref> are characterized by linear combinations of the multipole structures with different IREPs, which lowers the magnetic point group symmetry, and may lead to additional finite components in the AH conductivity and EM tensors.
For instance, the (M_u, T_z) multipole structure of low-temperature phase of ScMnO_3 only preserves 6 fold rotation symmetry along z-axis and can have finite σ_xy for the AH conductivity <cit.>, and finite α_xx=α_yy, α_zz, and α_xy=-α_yx for the EM coefficients <cit.>.
Meanwhile, the (M_3b, M_3a) multipole structure of YMnO_3 preserves the magnetic point group symmetry operation of 6 fold rotation along z-axis combined with the time reversal operation, and no AH and EM effects are expected.










§ SUMMARY

We have proposed a scheme to efficiently generate the symmetry-adapted orthonormal magnetic structures in crystallographic point group by introducing a virtual atomic cluster to perform the multipole expansion.
With this method, we can obtain magnetic structure that is fully characterized by the M and MT multipole as the suitable order parameters.
We have introduced a virtual atomic cluster to obtain the magnetic structures preserving the magnetic point group symmetry of the multipoles.
We have applied the proposed method to pyrochlore and hexagonal ABO_3 crystal structures.
For the pyrochlore crystal structure, we have investigated all-in all-out, two-in two-out, and one-in three-out magnetic structures and found that the two-in two-out and one-in three-out magnetic structures are able to be transformed continuously to the pure antiferromagnetic (octupole) structures without net magnetization, leading to the AH effect.
For the hexagonal ABO_3 crystal structure, the expression of the EM effect is fully identified by the multipole magnetic structures.
The proposed scheme paves a way to generate the multipole magnetic structures in compatible with the crystallographic point group, which is essential for macroscopic phenomena following Neumann's principle, and is useful to search for desired functional magnetic materials.



§ ACKNOWLEDGMENTS

This work was supported by JSPS KAKENHI Grant Numbers JP15K17713 (MTS),
JP15H05883 (J-Physics) (MTS), JP18H04230 (MTS), JP16H04021 (MTS), JP15K05176(HK),
JP15H05885 (J-Physics) (HK), JP16H06590(SH), JP18H04296
(J-Physics) (SH), JP18K13488 (SH), JP16H06345 (RA), JST
PREST (MTS), CREST Grant Number JPMJCR15Q5 (RA& MTS).
T. N. is supported by RIKEN Special Postdoctoral Researchers Program.


apsrev


