Physical realization of topological Roman surface by spin-induced … –

Thank you for visiting You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Carousel with three slides shown at a time. Use the Previous and Next buttons to navigate three slides at a time, or the slide dot buttons at the end to jump three slides at a time.
16 June 2021
Marc A. Wilde, Matthias Dodenhöft, … Christian Pfleiderer
17 December 2021
Chiara Devescovi, Mikel García-Díez, … Aitzol García-Etxarri
06 July 2020
Korbinian Geirhos, Boris Gross, … István Kézsmárki
09 February 2022
Dorin Rusu, Jonathan J. P. Peters, … Marin Alexe
25 January 2021
Kosuke Karube, Licong Peng, … Yasujiro Taguchi
06 May 2019
Niels B. M. Schröter, Ding Pei, … Yulin Chen
22 July 2020
Jia-Xin Yin, Wenlong Ma, … M. Zahid Hasan
21 September 2021
Ken Osumi, Tiantian Zhang & Shuichi Murakami
21 January 2022
M. J. Pearce, K. Götze, … P. A. Goddard
Nature Communications volume 13, Article number: 2373 (2022)
1921 Accesses
1 Citations
11 Altmetric
Metrics details
Topology, an important branch of mathematics, is an ideal theoretical tool to describe topological states and phase transitions. Many topological concepts have found their physical entities in real or reciprocal spaces identified by topological invariants, which are usually defined on orientable surfaces, such as torus and sphere. It is natural to investigate the possible physical realization of more intriguing non-orientable surfaces. Herein, we show that the set of spin-induced ferroelectric polarizations in cubic perovskite oxides AMn3Cr4O12 (A = La and Tb) reside on the topological Roman surface—a non-orientable two-dimensional manifold formed by sewing a Möbius strip edge to that of a disc. The induced polarization may travel in a loop along the non-orientable Möbius strip or orientable disc, depending on the spin evolution as controlled by an external magnetic field. Experimentally, the periodicity of polarization can be the same or twice that of the rotating magnetic field, which is consistent with the orientability of the disc and the Möbius strip, respectively. This path-dependent topological magnetoelectric effect presents a way to detect the global geometry of a surface and deepens our understanding of topology in both mathematics and physics.
Topological physics is one of the most successful integrations of the mathematics and physics field. As one of the branches of pure mathematics, topology has recently been employed to describe and understand many intriguing physical phenomena, such as the quantum Hall effect1,2, topological insulators, topological semimetals, and topological superconductors3,4,5. These topological electronic states are classified by various topological invariants, such as the Chern number N, topological classes Z2 and Z4, and Euler class χ6,7,8,9,10,11. It is known that the Chern number is defined as the integral of the vector field, i.e., the Berry curvature over the first Brillouin zone of a two-dimensional (2D) periodic system, which is considered as a torus (T2) and serves as the base manifold. Meanwhile, the magnetic spins and electric dipoles in real space can form topological objects or defects such as skyrmions12,13, polar skyrmions14, and polar merons15. These can be identified with an integer winding number16 based on the field of spin or electrical dipole orientations over the base manifold of a sphere (S2). In this sense, topological invariants are global characteristics, and topology-related physical properties are typically robust against perturbations3,4,8.
In addition to torus and sphere, there are three other fundamental 2D manifolds found in algebraic topology (see Supplementary Table 1): the Möbius strip, Klein bottle, and Roman surface. The torus and sphere are orientable and have two sides. The Berry curvature flux enclosed by them can be a nonzero quantized value and taken as a topological invariant, such as the Chern number. By contrast, the Möbius strip, Klein bottle, and Roman surface are non-orientable surfaces with only one side. The flux of any vector field through them is zero. A well-known example is that a one-sided Klein bottle cannot be filled with water, but a two-sided bottle can. Moreover, the non-orientable surfaces feature a geometrical property in which a body can be mirror-reflected when traveling along a special loop path on the surface, such as along the Möbius strip without crossing the boundary, and recovers after the second loop of travel. It will be intriguing to realize a physical system that can evolve globally on a non-orientable surface to reflect such unique geometrical properties.
In real space, a Möbius strip with a single boundary has been visualized in nanocrystal/DNA bands17,18. Möbius insulators have been proposed and realized in the reciprocal spaces of electronic19,20,21 and acoustic crystals22,23. Under proper nonsymmorphic symmetry or projective translation symmetry, the Bloch wavefunctions of its boundary states have 4(pi) periodicity instead of 2(pi), which mimics characteristics of a Möbius strip. Recently, nonorientability has been regarded as a Z2 gauge charge to describe the topological defect of shear deformation on a Möbius strip24. However, the smooth and global evolution of a physical quantity on a non-orientable surface and the detection of its intrinsic geometry remains a challenge. Furthermore, a non-orientable 2D surface without a boundary must self-intersect in a three-dimensional (3D) space, as shown by the Klein bottle25 (one-sided bottle) and Roman surface26 (RP2, real projective plane). This inhibits the physical realization and exhibition of topological behaviors in a 3D real space. Therefore, it is natural to ask whether the Roman surface, the simplest non-orientable 2D manifold without a boundary, can be realized by some physical entity and whether it can demonstrate exotic topological properties or phenomena.
Jakob Steiner discovered and named the Roman surface during his visit to Rome in early 1844. It is composed of a Möbius strip sewed to the edge of a disc and has a high tetrahedral Td symmetry, which is different from other non-orientable 2D surfaces. Three self-intersecting double lines Di (i = a, b, and c are three mutually perpendicular Cartesian axes) form a right trihedron in a Euclidean 3D phase space. Each double line ends with two pinch points, and they all intersect at the triple point (see Fig. 1 and Supplementary Fig. 1).
a, b Roman surface with four three-fold axes (3), one triple point, six pinch points, and three double lines (Da, Db, Dc). c Three geometrically different fundamental loops on a Roman surface. The red loop effectively winds twice to pass the two sides of the triple point as it rides on the Möbius strip being part of the Roman surface. The blue loop winds around the double line Dc once to have Wc = 1. The green loop does not enclose any double lines with zero Wi (i = a, b, c). Both blue and green loops only wind once to close the path.
The triple point on the Roman surface indicates the intersection of the three surfaces. The analytical equation of the Roman surface RP2 in mathematics is
where δ is a constant and (x, y, z) is a point in the a, b, and c coordinate systems. As a 2D manifold, there are three types of fundamental closed loops on a Roman surface, and other closed loops are formed by their combination. As Fig. 1c shows, the first type of loop is a red loop that passes through at least one point twice (e.g., this point can be the triple point). It winds twice to close the path on this non-orientable surface because it rides on the Möbius strip, which is a component of the Roman surface. The other two types of loops are blue and green. They are topologically different from the red ones because they do not pass any point twice and wind once to close the path. The blue loop can only wind along one of the double lines Di (i = a, b, c) by Wi = 1 time. The green loop does not enclose any double lines and is topologically trivial. The physical entity that matches Eq. (1) may provide an opportunity to realize a Roman surface and exhibit path-dependent properties. Therefore, realizing the Roman surface is intriguing and may extend the territory in the field of topological physics.
Magnetoelectric (ME) multiferroics with spin-induced ferroelectric polarization (P) have received considerable attention owing to their coupled spin and electrical dipole orders27, where the electric polarization is switched using a magnetic field (B), and the magnetization (M) is switched using an electric field (E)27,28, giving rise to significant ME effects. A few physical mechanisms have been proposed for ME multiferroics29,30,31,32, providing a fertile and unique playground for investigating topological properties as a mapping between the spin vector space and electric dipole vector space. Although some structurally and ferroelectrically coupled vortex domain patterns have been observed in hexagonal RMnO3 (R = rare earth or Y, In) and BiFeO327,33, spin-dependent topological multiferroics remain to be elucidated, the topological ME effects with versatile controls of P by B or M by E also need to be studied. Here, in A-site ordered cubic perovskites AMn3Cr4O12 with high lattice symmetry, we discovered that the set of spin-induced ferroelectric polarizations can be described by a non-orientable Roman surface, as deduced through a 3D symmetry analysis. Various angle-dependent P oscillations under different magnetic fields were observed and can be understood based on the changes in the local topological structure and symmetry of the magnetism.
To prove the high cubic symmetry of synthesized TbMn3Cr4O12 (TMCO) at room temperature, we measured its powder X-ray diffraction (XRD) pattern and performed the Rietveld refinement (see Fig. 2a). The Rietveld analysis demonstrates that the compound crystallizes in an A-site ordered perovskite structure with the chemical formula AA3B4O12 (see Fig. 2b)34,35. The space group is cubic Im-3 with central symmetry, where Tb and Mn are 1:3 ordered at special Wyckoff sites 2a (0, 0, 0) and 6b (0, 0.5, 0.5), and Cr and O at special 8c (0.25, 0.25, 0.25) and 24 g (x, y, 0) sites, respectively. Compared with a simple ABO3 perovskite, the introduction of a transition metal into the A′ site in the ordered AA3B4O12 perovskite results in the formation of square-coordinated A′O4 units and heavily tilted BO6 octahedra, as shown in Fig. 2b and described in detail in literatures36,37. Supplementary Table 2 lists the refined structural parameters of TMCO. According to the refined Cr–O and Mn–O bond lengths, the bond valence sum calculations illustrate that the valence states of Mn and Cr are both close to +3, which agrees with the reported isostructural compound LaMn3Cr4O12 (LMCO)38. Moreover, when the temperature was decreased to 5 K, there is no trace of the structural phase transition in TMCO (see Supplementary Fig. 2).
a XRD pattern measured at room temperature and the Rietveld refinement results. The observed (black circles), calculated (red line), and difference (blue line) values are shown. The ticks indicate the allowed Bragg reflections with space group Im-3. b Crystal structure of A-site ordered perovskite AA3B4O12 with Im-3 symmetry. c Temperature dependence of magnetic susceptibility χ measured at 0.1 T and its inverse χ-1. Both zero-field-cooling (ZFC) and field-cooling (FC) modes were adopted to measure χ. The yellow line shows the Curie–Weiss fitting above 175 K. The fitted effective moment μeff (15.52 μB f.u.-1) is comparable to the theoretical value (15.04 μB f.u.-1) by considering the contribution from Tb3+, Mn3+ and Cr3+ ions. The inset shows the derivative of inverse susceptibility versus temperature. d Temperature dependence of specific heat Cp measured at zero field. e, f Temperature dependent polarization P after poling the sample under an electric field (= Epole) and different magnetic fields B, where B was not removed until the temperature increased to 50 K. For B//E and BE configurations, both +Poled and -Poled conditions were measured.
Figure 2c, d shows the temperature dependence of the magnetic susceptibility and specific heat of TMCO, respectively. With the temperature decreasing to TN1 ≈ 136 K, a sharp λ-type anomaly is observed in the specific heat, suggesting the occurrence of a second-order long-range magnetic phase transition. Moreover, the derivative of the inverse susceptibility displays a variation at TN1 (inset of Fig. 2c). Further, upon cooling to TN2 ≈ 36 K, both the susceptibility and specific heat show a clear anomaly, illustrating a second long-range magnetic transition. Below TN2, the specific heat does not show any visible anomaly, indicating the absence of spin ordering of Tb3+ at temperatures down to 2 K. The inverse susceptibility above 175 K follows the Curie–Weiss law (Fig. 2c). The fitted Weiss temperature of −28.0 K is indicative of antiferromagnetic (AFM) interactions. The linear isothermal magnetization behaviors within ±2 T further confirm the AFM orders occurring at TN1 and TN2 (see Supplementary Fig. 3a). However, at higher fields, a metamagnetic transition occurs near 2.5 T at temperatures below TN1 ≈ 136 K, as demonstrated by the field derivative of magnetization, indicating field-induced variation in spin alignment (Supplementary Fig. 3b). Similar to LMCO38, at zero field, the B-site Cr3+ spins are antiferromagnetically ordered in a collinear G-type fashion at TN1; the A′-site Mn3+ sublattice experiences another G-type AFM ordering at TN2 in the current TMCO. The total spin structure, composed of Cr3+ and Mn3+ magnetic sublattices with spin moments along the [111] direction, forms a type-II polar magnetic point group of 31′, which breaks the space inversion symmetry. Therefore, a spin-induced ferroelectric phase transition is expected to occur at TN2.
Corresponding to the AFM phase transition, a frequency-independent dielectric peak is observed at TN2 ≈ 36 K in TMCO, whereas no dielectric anomaly emerges at TN1 (Supplementary Fig. 4a, b). When magnetic field B and electric field E are simultaneously applied to measure the dielectric constant, an apparent anisotropic behavior can be observed. Specifically, a sharp dielectric peak is always observed with B up to 9 T for B//E (Supplementary Fig. 4c). By contrast, B significantly suppresses the dielectric peak in the BE configuration, whereas the dielectric anomaly nearly disappears at 6 T (Supplementary Fig. 4d). These results strongly suggest that spin-induced ferroelectric polarization also shows distinctive responses to external magnetic fields along the parallel and perpendicular configurations in TMCO38. Furthermore, when the pyroelectric current Ip was measured to obtain temperature-dependent polarization, a sharp ferroelectric phase transition was observed at TN2 (Fig. 2e), as expected from the formation of the polar magnetic point group mentioned earlier. The polarization is reversible without changing the magnitude if the sign of the poling electric field (Epole) is reversed, further confirming the spin-induced ferroelectric phase transition in TMCO. Similar to the dielectric measurements, applying a magnetic field with a B//E configuration can significantly enhance the polarization from 33.3 μC/m2 at 0 T to 84.4 μC/m2 at 6 T (Fig. 2e). However, in the BE configuration, the spin-induced polarization is suppressed to 10.7 μC/m2 at B = 3 T, and only 4.51 μC/m2 at B = 6 T (Fig. 2f). All experimental results indicate that TMCO has the same crystal and magnetic structures as LMCO38.
Based on the crystal and spin structure symmetry analyses, the peculiar ME multiferroicity of AMn3Cr4O12 can be described through the topological Roman surface RP2. As mentioned, both TMCO and LMCO have a cubic crystal symmetry of Im-3, where inversion centers are located on each Cr atom, with four three-fold rotation axes along four body diagonals and three mirror planes perpendicular to the three Cartesian axes (Fig. 3a). The total magnetic structure composed of Cr3+ spins (we use Néel vector SCr to trace the coherent motion of spins in the AFM ordering of Cr ions) and Mn3+ spins (Néel vector SMn) forms a polar magnetic point group 31′ with SMn//SCr//[111] (see Fig. 3a). According to the proposed anisotropic symmetric exchange mechanism39 (see Supplementary Section 1) originating from spin–orbit coupling, such a special spin structure induces polarization P along the [111] direction. To describe a general relationship between the spin and polarization directions, we assign SCr = (μCr, θ, φ) and SMn = (μMn, θ, φ). Both are on sphere S2 (Fig. 3a and b), where (theta in left[0,pi right]) and (varphi in [-pi ,pi ]) with fixed magnitudes μCr and μMn. Thus, the polarization P is expressed as a function of spins P(SMn, SCr) = P(θ, φ) mapping from S2 to RP2. Based on strict symmetrical analysis40, as described in Supplementary Sections 2 and 3, we see that P(θ, φ) has the following form:
where the coefficients β and Ω depend on temperature and crystal properties, respectively. The polarization vector P(θ, φ) = (Pa, Pb, Pc) is rewritten using an analytical equation as follows:
a Crystal and spin configurations where Cr and Mn atoms are presented in yellow and blue, respectively, and O atoms are omitted for clarity. The lattice symmetries of AMn3Cr4O12: two space inversions (bar{{{{{{boldsymbol{1}}}}}}}) located at the center and at the Cr atom, one of the mirror planes (see the light-blue plane), and one of the three-fold rotation axes (see the dash line) are indicated. All spins composed of Mn and Cr magnetic sublattices are parallel or antiparallel to the [111] direction. Antiferromagnetic sublattices of Cr and Mn are denoted by Néel vectors SCr and SMn, respectively, and the induced polarization P is shown by the arrow in green. b Orientation of SCr is represented in spherical coordinates (μCr, θ, φ). c Complete set of allowed spin-induced P vectors in three dimensions is a Roman surface when the magnetic structure shown in a coherently rotates on a sphere in b. d Distribution of P vectors along four of the eight diagonal directions. The yellow sphere and the black arrow denote the Cr atom and its spin moment. SMn is always parallel to SCr, and it is not shown for clarity. Thus, AFM configuration composed of Mn and Cr sublattices is indicated by SCr. The black cube represents the crystal unit cell.
Equation (3) is a specific case of Eq. (1) with x = Pa, y = Pb, z = Pc, and δ = 2βμCrμMnΩ. This result indicates that the induced P vectors in AMn3Cr4O12 may constitute a non-orientable Roman surface RP2 (see Fig. 3c).
In this cubic multiferroic family, the collinear G-type AFM locks the Néel vectors SMn and SCr in parallel; they pass through all eight equivalent diagonal directions coherently (see Supplementary Section 3 for details). Only four of the eight diagonal directions are allowed for the polarization P vectors, which form a regular tetrahedron and agrees with the tetrahedral symmetry of the Roman surface (Fig. 3d). Note that there is another AFM configuration with SMn antiparallel to SCr (see Supplementary Fig. 5), which leads to a set of polarization vectors P forming a Roman surface with negative δ and having similar features as the positive one discussed here.
The spin–polarization relationship in cubic perovskites AMn3Cr4O12 leads to topological multiferroicity owing to the topological and symmetry properties of the Roman surface discussed earlier. The mapping from S2 to RP2 covers the Roman surface twice, because SCr and –SCr induce identical P vectors. This is similar to the spin–momentum locking effect on the Fermi surface enclosing Weyl nodes. The Hamiltonian around the Weyl node maps the momentum space (Fermi sphere, S2) to the spin space (another sphere called the Bloch sphere). If the enclosed Weyl node have a topological charge of one or two, the momentum moves around the node once and the spin rotates once or twice, respectively. The relative rotational directions of the momentum and spin are determined by the sign of the charge41. In this sense, a closed evolution loop of Néel vectors SCr on the sphere may map to a non-trivial path of the P vector on the Roman surface, showing topological ME effects.
The red loop in Fig. 1c passes through the triple point of the Roman surface twice, where polarization P = 0. This can be induced by the evolution of SCr on sphere S2, which includes six special points where the sphere intersects the a-, b-, and c-axis (see Fig. 4a). For a longitudinal great circle of the sphere with φ fixed, as Fig. 4a shows, SCr passes through the special points twice (the north and south poles) if φ is not 0, π/2, π, or 3π/2. This results in a double-winding elliptical trajectory of the P vector passing through the triple point twice (see Fig. 4b and Supplementary Video 1). The induced P rotates twice, whereas SCr travels through one loop. When φ = 0, π/2, π, or 3π/2, the trajectories of P collapse to the double lines and pass through special points four times (see Fig. 4b and Supplementary Videos 1 and 2). It is observed that only the magnitude of P vibrates, and the period is twice the loop evolution of SCr.
a Longitudinal (fixed φ) path of the Néel vector rotation mode with a series of selected φ values and continuously varying (theta in left[{{{{mathrm{0,2}}}}}pi right]). b Trajectories of P induced by loops in a with the same color. c Latitudinal (fixed θ) Néel vector rotation mode with a series of selected θ values and continuously varying (varphi in left[-pi ,pi right]). d Trajectories of P induced by loops in c with the same color. The six special points, the intersections between sphere and a-, b– and c-axis, are illustrated in black point in a and c. The thicker lines in a and c represent two special cases where four diagonal directions are passed through; the resultant trajectories of P are indicated by the thicker lines in b and d, respectively. Rotation directions of SCr and P are indicated by black arrows toward opposite directions away from [111] direction. The animations of ad are shown in Supplementary Videos 13. e Two closed loops winding the a and b axes once are shown in orange and purple, respectively. f Two closed trajectories of P on the Roman surface are induced by the loops in e. The loops with Wa = 1 are in orange and those with Wb = 1 are in purple. g, h A band included in the equator of the Néel vector sphere g is mapped to h, a self-crossed Möbius strip as a part of the Roman surface. The residual upper and lower caps are mapped to a disc on the Roman surface simultaneously. All the loops in b, d, and f are fundamental loops defined in Fig. 1c.
When the Néel vectors evolve along the latitudinal parallels with the polar axis being the c-axis and a fixed θ value (θ is not 0 or π/2), as Fig. 4c shows, the induced saddle-like trajectory of P on the Roman surface only winds around the c-axis (Fig. 4d and Supplementary Video 3). This trajectory of P is nontrivial, in the sense that Wc = 1. The projection onto the ab-plane is a circle. The rotation direction, clockwise or counterclockwise, of the projection of the P vector on the circle is opposite to that of the Néel vector. Note that the rotation period of the P vector is the same as that of SCr, instead of double that in the cases shown in Fig. 4b. Similarly, when SCr winds around either the a– or b-axis once, the corresponding P vector loop on the Roman surface is characterized by Wa = 1 or Wb = 1, respectively, as indicated in Fig. 4e and f. Notably, nonzero Wi corresponds to the component of the P vector perpendicular to the i-axis that reverses Wi times during one loop evolution. The green circle in Fig. 1c is topologically trivial with zero values of Wi. It can be mapped from a loop on the SCr sphere that satisfies these conditions: there is no special point being wound, and –SCr will not be passed by for any SCr on the path. To reveal the underlying topological origins of the double-winding and single-winding paths in Figs. 4b, d, and 1c, we mapped a strip composed of the neighborhoods of the equator and itself on the spin sphere (Fig. 4g) to a Möbius strip (Fig. 4h) on the Roman surface. Meanwhile, the residual upper and lower caps on the sphere are mapped to a saddle-like disc. Thus, according to the orientable and non-orientable nature of the disc and the Möbius strip, their paths on them will wind once and twice, respectively, to form a closed loop.
We experimentally present the topological ME effects in TMCO and LMCO polycrystalline samples with identical crystal and spin structures. To define the relative orientations of E and B on the polycrystalline sample, we defined a local x, y, and z coordinate system (see Fig. 5a). First, the samples were cooled down to 10 K under Epole//y at zero B. Subsequently, we removed the electric field Epole and measured the Bdependent P by sweeping B in either the B//y or B//x configuration (Fig. 5a). In both configurations, quadratic ME effects were dominant in TMCO (Fig. 5b). This is consistent with the 31′ point group with independent time-reversal symmetry, where a quadratic rather than a linear ME effect is allowed. To investigate the postulated topological ME effects in TMCO, rotating B was applied to continuously tune the orientation of the Néel vector SCr. The B direction dependence of Py under a selected field-strength by changing ψ (initially ψ = 0 with B//y; Fig. 5c–f) was measured after electrical poling to 10 K. The ψ dependence of Py exhibited double-angle sinusoidal modulation curves at B within ~2 T without changing the sign (Fig. 5c). This is a direct manifestation of the intrinsic quadratic ME behavior in this system. Based on Monte Carlo simulations, we find that under a relatively small rotating B without inducing metamagnetic transition, the Néel vector SCr will slightly tilt away from the diagonal direction, leading to a topologically trivial path on the sphere. However, it is to be noted that SCr itself is double period of rotating B (see Supplementary Section 4 and Supplementary Fig. 6a). Moreover, the resulting loop of the P vector on the Roman surface is a small perturbation to its original position, with zero values of Wi (Supplementary Fig. 6b), as exemplified by the green loop in Fig. 1c. Therefore, SCr travels along a green-type path twice when B rotates 360° (Supplementary Fig. 7). This results in a double winding of the loop on the Roman surface, which is consistent with the experimentally observed double-angle sinusoidal modulation of Py for small B, indicating trivial mapping from SCr to P along the green loop. Additionally, under applied magnetic field, B will change the orientation of antiferromagnetic axis and induce a small ferromagnetic component simultaneously. The Zeeman energy and antiferromagnetic exchange energy will compete with each other. The AFM components of spins tend to align normally with the direction of B while the induced ferromagnetic components point towards the B direction. This is also consistent with our Monte Carlo simulations under a small B in Supplementary Figs. 6a and 7a where B//b (Φ = 90°) leads to a significant weakening of ({S}_{{Cr}}^{y}) in spins.
a Schematic illustration of measurements for magnetoelectric current under rotating B within the xy-plane; ψ denotes the angle between B and the y-direction. b, g Magnetic-field-dependent Py at T = 10 K after +E poling for B//E and BE configurations of TMCO and LMCO, respectively. Py is defined as the component of P along the y-axis. cf ψ-dependent polarization of TMCO at B = 2 T c, 3 T d, 6 T e, and 9 T f. The red and black arrows indicate an increase or decrease in angle, respectively. h, i ψ-dependent polarization of LMCO at B = 1 T h and 6 T i. All polarizations were integrated from the ME current as a function of time.
For B values greater than 3 T, a sinusoidal Py variation appeared, owing to the metamagnetic transition occurring at a critical field of approximately 2.5 T (Supplementary Fig. 3b). At B = 6 and 9 T, the Py value can be reversed with pure sinusoidal behavior (Fig. 5e, f). The energy barrier of magnetic anisotropy in the diagonal directions can be overcome by increasing the Zeeman energy, and a spin-flop transition occurs at a larger rotating B ( > 3 T) to induce a nontrivial loop of the Néel vector on the sphere. The loop of the Néel vector winds its principal axis under a sufficiently high B, which induces a P vector loop with one of Wi = 1 (i = a, b, c) on the Roman surface. This leads to a reversal of the P components perpendicular to the i-axis, as the blue loop in Fig. 1c. Thus, the pure sinusoidal P behavior at 6 T and 9 T in TMCO can be understood.
As Fig. 5d shows, both single- and double-angle sinusoidal Py features coexist at B = 3 T with clear hysteresis behavior. P = 0 values have been reached multiple times in one period. There is a possibility that the loops of the Néel vector pass through some of the crystal axes, and some naught P values occur during a 360° rotation of B, as in the red loop in Fig. 1c. Note that owing to the polycrystalline specimens used in the experiments, it is impossible to determine the detailed grain configurations or P rotating paths of each grain for the TMCO. In isostructural LMCO, all the aforementioned ME behaviors can be reproduced (Fig. 5g–i), indicating the universality of these phenomena in the cubic multiferroic perovskite family of AMn3Cr4O12.
We propose that topological nontrivial rotation of P is realized in the multiferroic domain walls of AMn3Cr4O12 without an external magnetic field. As Supplementary Fig. 8 shows, there are three types of multiferroic domains with SMn parallel to SCr, considering that the Néel vectors prefer eight diagonal directions, i.e., 71°, 109°, and 180°. AFM domains are induced at 109°, 109°, and 360°. In particular, ferroelectric domains are induced at 360°. The ferroelectric domain wall represents a half winding of the Möbius strip since a 180° rotation of SCr leads to one loop of P and a mirror reflected on the Roman surface (Supplementary Fig. 8c). This type of multiferroic domain wall cannot be destroyed by an external electric field because the polarizations in two neighboring ferroelectric domains can be coherently aligned in the direction of the electric field; whereas, the related AFM configurations will always be opposite. This would be a type of topological object unique in AMn3Cr4O12, which deserves further investigation. Moreover, it is worthy to note that the magnetic skyrmion of Cu2OSeO342 may also map the induced polarization as a complete Roman surface owing to the similar polarization distribution to that of AMn3Cr4O12.
In summary, we investigated the unusual spin–polarization relationship in A-site-ordered cubic perovskite oxides AMn3Cr4O12. The uncommon polarization occurring in the cubic system is coupled with the inversion symmetry breaking caused by G-type collinear AFM ordering composed of Mn3+ and Cr3+ magnetic sublattices. We discovered that the ME behavior, namely, the mapping from spin to the induced polarization, can be considered as a physical entity realizing the Roman surface, which is a non-orientable 2D manifold. This is different from other topological objects defined on orientable manifolds, such as Chern insulators defined on a torus and skyrmions defined on a sphere. The non-orientability of the Roman surface leads to topological ME effects depending on the evolution of the path. In our experiments, different types of angle-dependent polarization are observed in polycrystalline TMCO and LMCO under different magnetic fields. The evolution paths of the spin vectors are controlled by a magnetic field, and the induced P evolves along orientable or non-orientable surfaces to show different topological properties. This study provides an interesting fact that topology has a close and deep relationship with the physical world, opening up an avenue to achieve more topological objects from an orientable manifold to a non-orientable manifold.
Polycrystalline TMCO was synthesized via a solid-state reaction at a high temperature and high pressure. Tb4O7, MnO, Cr2O3, and Mn2O3 powders were used as the starting materials at a molar ratio of 1:2:8:5. The finely mixed reactants were packed into a gold capsule 4.0 mm in length and 2.8 mm in diameter and then treated at 1300 K and 7 GPa for 30 min on a cubic-anvil-type high-pressure apparatus.
Powder XRD analysis was performed using a Huber diffractometer with Cu Kα1 radiation at 40 kV and 30 mA. Diffraction data were collected in the 2θ angle range of 10° to 100°. A program from the general structure analysis system was used to refine the XRD data using the Rietveld method43. Raman spectra were collected on a T64000 spectrometer using a liquid-nitrogen-cooled charge-coupled device. Volume Bragg gratings of 600 G and 488 nm laser were used to measure the spectra.
The magnetic susceptibility and magnetization were measured using a superconducting quantum interference device magnetometer (Quantum Design, MPMS-7T). Zero-field-cooling and field-cooling modes were adopted to measure the susceptibility at a magnetic field of 0.1 T. Magnetization was measured between −14 T and 14 T at selected temperatures. Specific heat measurements were performed using a physical property measurement system (Quantum Design, PPMS-9T) at zero field.
A hard disk-shaped piece 2.0 mm in diameter and 0.2 mm in thickness was used to measure the relative dielectric constant on an Agilent-4980A LCR meter. Silver paste was used as the electrodes. Subsequently, the sample was used to measure the pyroelectric current and isothermal polarized current under different magnetic fields using a Keithley 6517B electrometer. To measure the pyroelectric current for TMCO, the sample was poled with an electric field of Epole =±  11.7 kV/cm and a magnetic field from 50 K to 10 K. Then, the electric field was removed. After waiting for 30 min in an electric short circuit, pyroelectric current data were collected at temperatures ranging from 10 K to 50 K at a speed of 2 K/min. Polarization P was obtained by integrating the current data as a function of time. The polycrystalline sample LMCO was also prepared in a disk-shaped piece with 0.23 mm in thickness to measure the polarization. In the measurements of magnetic-field dependent and angle-dependent polarizations, the polycrystalline sample was poled by an electric field of 8.33 kV/cm for TMCO and 8.68 kV/cm for LMCO down to 10 K. After the electrical poling procedure, different magnetic field was applied to the sample to measure the isothermal polarized current; then, the magnetic-field-dependent polarization or angle-dependent polarization was obtained by integrating the current.
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
The data supporting this study’s findings are available from the corresponding authors upon reasonable request.
The codes in this study are available from the corresponding authors upon reasonable request.
Klitzing, K. V., Dorda, G. & Pepper, M. New method for high-accuracy determination of the fine-structure constant based on quantized hall resistance. Phys. Rev. Lett. 45, 494–497 (1980).
Article  ADS  Google Scholar 
Weng, H., Yu, R., Hu, X., Dai, X. & Fang, Z. Quantum anomalous Hall effect and related topological electronic states. Adv. Phys. 64, 227–282 (2015).
Article  ADS  Google Scholar 
Hasan, M. Z. & Kane, C. L. Colloquium: topological insulators. Rev. Mod. Phys. 82, 3045–3067 (2010).
Article  ADS  CAS  Google Scholar 
Qi, X.-L. & Zhang, S.-C. Topological insulators and superconductors. Rev. Mod. Phys. 83, 1057–1110 (2011).
Article  ADS  CAS  Google Scholar 
Armitage, N. P., Mele, E. J. & Vishwanath, A. Weyl and Dirac semimetals in three-dimensional solids. Rev. Mod. Phys. 90, 015001 (2018).
Article  ADS  MathSciNet  CAS  Google Scholar 
Thouless, D. J., Kohmoto, M., Nightingale, M. P. & den Nijs, M. Quantized Hall conductance in a two-dimensional periodic potential. Phys. Rev. Lett. 49, 405–408 (1982).
Article  ADS  CAS  Google Scholar 
Berry, M. V. Quantal phase factors accompanying adiabatic changes. Proc. R. Soc. Lond. Ser. A 392, 45 (1984).
Article  ADS  MathSciNet  Google Scholar 
Kane, C. L. & Mele, E. J. Z2 topological order and the quantum spin Hall effect. Phys. Rev. Lett. 95, 146802 (2005).
Article  ADS  CAS  Google Scholar 
Song, Z., Zhang, T., Fang, Z. & Fang, C. Quantitative mappings between symmetry and topology in solids. Nat. Commun. 9, 3530 (2018).
Article  ADS  Google Scholar 
Ahn, J., Park, S. & Yang, B.-J. Failure of Nielsen-Ninomiya theorem and fragile topology in two-dimensional systems with space-time inversion symmetry: application to twisted bilayer graphene at magic angle. Phys. Rev. X 9, 021013 (2019).
CAS  Google Scholar 
Bouhon, A. et al. Non-Abelian reciprocal braiding of Weyl points and its manifestation in ZrTe. Nat. Phys. 16, 1137–1143 (2020).
Article  CAS  Google Scholar 
Mühlbauer, S. et al. Skyrmion lattice in a chiral magnet. Science 323, 915–919 (2009).
Article  ADS  Google Scholar 
Yu, X. Z. et al. Real-space observation of a two-dimensional skyrmion crystal. Nature 465, 901–904 (2010).
Article  ADS  CAS  Google Scholar 
Das, S. et al. Observation of room-temperature polar skyrmions. Nature 568, 368–372 (2019).
Article  ADS  CAS  Google Scholar 
Wang, Y. J. et al. Polar meron lattice in strained oxide ferroelectrics. Nat. Mater. 19, 881–886 (2020).
Article  ADS  Google Scholar 
Nagaosa, N. & Tokura, Y. Topological properties and dynamics of magnetic skyrmions. Nat. Nanotechnol. 8, 899–911 (2013).
Article  ADS  CAS  Google Scholar 
Tanda, S. et al. A Möbius strip of single crystals. Nature 417, 397–398 (2002).
Article  CAS  Google Scholar 
Han, D., Pal, S., Liu, Y. & Yan, H. Folding and cutting DNA into reconfigurable topological nanostructures. Nat. Nanotechnol. 5, 712–717 (2010).
Article  ADS  CAS  Google Scholar 
Young, S. M. & Kane, C. L. Dirac semimetals in two dimensions. Phys. Rev. Lett. 115, 126803 (2015).
Article  ADS  Google Scholar 
Wang, Z., Alexandradinata, A., Cava, R. J. & Bernevig, B. A. Hourglass fermions. Nature 532, 189–194 (2016).
Article  ADS  CAS  Google Scholar 
Zhang, R. X., Wu, F., & Das Sarma, S. Mӧbius insulator and higher-order topology in MnBi2nTe3n+1. Phys. Rev. Lett. 124, 136407 (2020).
Article  ADS  CAS  Google Scholar 
Li, T. et al. Acoustic Möbius insulators from projective symmetry. Phys. Rev. Lett. 128, 116803 (2022).
Article  ADS  CAS  Google Scholar 
Xue, H. et al. Projectively enriched symmetry and topology in acoustic crystals. Phys. Rev. Lett. 128, 116802 (2022).
Article  ADS  CAS  Google Scholar 
Bartolo, D. & Carpentier, D. Topological elasticity of nonorientable ribbons. Phys. Rev. X 9, 041058 (2019).
CAS  Google Scholar 
Hilbert, D. & Cohn-Vossen, S. Geometry and the Imagination (Chelsea, New York, 1999).
Gray, A. Modern Differential Geometry of Curves and Surfaces with Mathematica. 2nd ed (CRC Press, Boca Raton, 1997).
Dong, S., Liu, J.-M., Cheong, S.-W. & Ren, Z. Multiferroic materials and magnetoelectric physics: Symmetry, entanglement, excitation, and topology. Adv. Phys. 64, 519–626 (2015).
Article  ADS  CAS  Google Scholar 
O′Dell, T. H. The Electrodynamics of Magneto-electric Media. (North-Holland, Amsterdam, 1970).
Google Scholar 
Jia, C., Onoda, S., Nagaosa, N. & Han, J. H. Bond electronic polarization induced by spin. Phys. Rev. B 74, 224444 (2006).
Article  ADS  Google Scholar 
Sergienko, I. A. & Dagotto, E. Role of the Dzyaloshinskii-Moriya interaction in multiferroic perovskites. Phys. Rev. B 73, 094434 (2006).
Article  ADS  Google Scholar 
Sergienko, I. A., Sen, C. & Dagotto, E. Ferroelectricity in the magnetic E-phase of orthorhombic perovskites. Phys. Rev. Lett. 97, 227204 (2006).
Article  ADS  Google Scholar 
Murakawa, H., Onose, Y., Miyahara, S., Furukawa, N. & Tokura, Y. Ferroelectricity induced by spin-dependent metal-ligand hybridization in Ba2CoGe2O7. Phys. Rev. Lett. 105, 137202 (2010).
Article  ADS  CAS  Google Scholar 
Lu, C., Wu, M., Lin, L. & Liu, J.-M. Single-phase multiferroics: New materials, phenomena, and physics. Natl Sci. Rev. 6, 653–668 (2019).
Article  CAS  Google Scholar 
Shimakawa, Y. A-site-ordered perovskites with intriguing physical properties. Inorg. Chem. 47, 8562–8570 (2008).
Article  CAS  Google Scholar 
Long, Y. W. et al. Temperature-induced A-B intersite charge transfer in an A-site-ordered LaCu3Fe4O12 perovskite. Nature 458, 60–63 (2009).
Article  ADS  CAS  Google Scholar 
Long, Y. A-site ordered quadruple perovskite oxides AA’3B4O12. Chin. Phys. B 25, 078108 (2016).
Article  Google Scholar 
Long, Y. & Shimakawa, Y. Intermetallic charge transfer between A-site Cu and B-site Fe in A-site-ordered double perovskites. N. J. Phys. 12, 063029 (2010).
Article  Google Scholar 
Wang, X. et al. Observation of magnetoelectric multiferroicity in a cubic perovskite system: LaMn3Cr4O12. Phys. Rev. Lett. 115, 087601 (2015).
Article  ADS  Google Scholar 
Feng, J. S. & Xiang, H. J. Anisotropic symmetric exchange as a new mechanism for multiferroicity. Phys. Rev. B 93, 174416 (2016).
Article  ADS  Google Scholar 
Chai, Y. S., Chun, S. H., Cong, J. Z. & Kim, K. H. Magnetoelectricity in multiferroic hexaferrites as understood by crystal symmetry analyses. Phys. Rev. B 98, 104416 (2018).
Article  ADS  Google Scholar 
Xu, G., Weng, H., Wang, Z., Dai, X. & Fang, Z. Chern semimetal and the quantized anomalous Hall effect in HgCr2Se4. Phys. Rev. Lett. 107, 186806 (2011).
Article  ADS  Google Scholar 
Seki, S., Ishiwata, S. & Tokura, Y. Magnetoelectric nature of skyrmions in a chiral magnetic insulator Cu2OSeO3. Phys. Rev. B 86, 060403 (2012).
Article  ADS  Google Scholar 
Larson, A. C. & von Dreele, R. B. General Structure Analysis System (GSAS). Report No. LAUR 86-748 (Los Alamos National Laboratory, 1994).
Download references
This study was supported by the National Key R&D Program of China (Grant No. 2021YFA1400300, Y.L.; 2018YFE0103200, Y.L.; and 2018YFA0305700, Y.L.), the Beijing Natural Science Foundation (Grant No. Z200007, Y.L.), the National Natural Science Foundation of China (Grant No. 11934017, Y.L.; 11974065 Y.C.; 11921004, Y.L.; 11904392, L.Zhou; 51772324, Y.L.; 11925408, H.W.; and 12188101, H.W.), the Chinese Academy of Sciences (Grant No. XDB33000000, Y.L.), Fundamental Research Funds for the Central Universities (2020CDJQY-A056, Y.C.) and Projects of President Foundation of Chongqing University (2019CDXZWL002, Y.C.). We acknowledge the fruitful discussion with X. H. Chen and W. M. Liu.
These authors contributed equally: Guangxiu Liu, Maocai Pi.
Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing, China
Guangxiu Liu, Long Zhou, Zhehong Liu, Xudong Shen, Xubin Ye, Shijun Qin, Lin Zhao, Bowen Zhou, Jia Guo, Xiaohui Yu, Hongming Weng & Youwen Long
School of Physical Sciences, University of Chinese Academy of Sciences, Beijing, China
Guangxiu Liu, Zhehong Liu, Xubin Ye, Shijun Qin, Lin Zhao, Bowen Zhou, Jia Guo, Xiaohui Yu, Hongming Weng & Youwen Long
Center of Quantum Materials and Devices, Chongqing University, Chongqing, China
Maocai Pi, Xinrun Mi, Xue Chen & Yisheng Chai
Low Temperature Physics Laboratory and Chongqing Key Laboratory of Soft Condensed Matter Physics and Smart Materials, College of Physics, Chongqing University, Chongqing, China
Maocai Pi, Xinrun Mi, Xue Chen & Yisheng Chai
Songshan Lake Materials Laboratory, Dongguan, Guangdong, China
Xudong Shen, Hongming Weng & Youwen Long
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
Y.L., Y.C., and H.W. conceived and led this project. M.P., X.C., and Y.C. performed the theoretical calculations. G.L., X.Ye, B.Z., and J.G. synthesized the samples. G.L., L.Zhou, X.S., and X.M. performed dielectric constant and polarization measurements and data analysis. Z.L., S.Q., L.Zhao, and X.Yu performed the structural and magnetic measurements. G.L., M.P., H.W., Y.C., and Y.L. wrote the manuscript. All the authors discussed the results and commented on the manuscript.
Correspondence to Yisheng Chai, Hongming Weng or Youwen Long.
The authors declare no competing interests.
Nature Communications thanks the other anonymous reviewer(s) for their contribution to the peer review of this work. Peer review reports are available.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit
Reprints and Permissions
Liu, G., Pi, M., Zhou, L. et al. Physical realization of topological Roman surface by spin-induced ferroelectric polarization in cubic lattice. Nat Commun 13, 2373 (2022).
Download citation
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Nature Communications (Nat Commun) ISSN 2041-1723 (online)
© 2023 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.


Leave a Reply

Your email address will not be published. Required fields are marked *