Shubha
Dubey
*a,
Vipin
Kumar
*bcd,
Kumud
Dubey
a,
Kuldeep
Singh
e,
Marcus
Einert
d,
Dipti
Gawande
a,
R. K.
Sharma
f,
Gitanjali
Pagare
g and
N. K.
Gaur
a
aDepartment of Physics, Barkatullah University, Bhopal 462 026, India. E-mail: shubha.dubey4@gmail.com
bDepartment of Physical Electronics, School of Electrical Engineering, Tel Aviv University, Tel Aviv 699 780, Israel. E-mail: vipinkumar@mail.tau.ac.il
cThe Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 699 780, Israel
dSurface Science Laboratory, Department of Materials- and Geosciences, Technical University of Darmstadt, Peter-Grünberg-Straße 4, 64287 Darmstadt, Germany. E-mail: vkumar@surface.tu-darmstadt.de
eSchool of Sciences, Noida International University, Greater Noida, Uttar Pradesh-201308, India
fTechnical Physics Division, Bhabha Atomic Research Centre, Mumbai 400 085, India
gDepartment of Physics, Sarojini Naidu Govt. Girls PG Autonomous College, Bhopal 462016, India
First published on 10th April 2025
The intricate interaction of emergent degrees of freedom, ground state degeneracy, and competing magnetic interactions in oxide-based double perovskite frustrated magnets can give rise to exotic excitations and correlated quantum phenomena with promising technological applications. This study investigates the structural, magnetic, and electronic properties of Dy2CoMnO6 synthesized via solid-state reaction. Rietveld refinement of XRD data confirms a monoclinic P21/n structure, with Raman spectra revealing BO6 polyhedral dynamics. XPS analysis indicates mixed valence states of Co and Mn. Magnetic measurements show ferromagnetic ordering below Tc = 87 K due to Co2+–Mn4+ superexchange, along with a Griffiths phase above Tc (TG = 91 K) and canted antiferromagnetic correlations at low temperature. A coercivity of 6.8 kOe at 10 K further supports FM behavior. Density functional theory (DFT) calculations validate the stability of Dy2CoMnO6 in the ferromagnetic phase, showing that the ferromagnetic state is energetically more favorable than the antiferromagnetic state. The direct band gap of Dy2CoMnO6 was found to be approximately 1.47 eV, while the indirect band gap was around 1.04 eV. Electronic band structure and density of states calculations predict band gaps of 1.60 eV and 3.20 eV for the minority and majority spin orientations, respectively, confirming the semiconductor nature of Dy2CoMnO6. These calculated values closely match the experimentally observed band gap.
Rare-earth-based double perovskites (A2B′B′′O6) with rock-salt ordering have garnered significant attention due to their intriguing phenomena, including magnetoelectric coupling, multiferroicity, spin-polarized conductivity, and superconductivity.9,10 Particularly, systems with B-site cations like Co, Mn, and Ni exhibit insulating behavior and high-temperature ferromagnetism, drawing interest for their potential in multifunctional devices.11–16 The RE2CoMnO6 (RCMO) family, where RE denotes a rare-earth element, offers remarkable tunability in electronic and magnetic properties through modifications of the rare-earth ionic radii (rR3+) and external pressure, making them suitable for spintronic applications.17–20
The magnetic properties of RE2CoMnO6 compounds are dominated by superexchange interactions between Co2+ and Mn4+ ions, governed by the Goodenough–Kanamori rules.21,22 Structural distortions and variations in the ionic radii significantly influence these interactions, impacting the magnetic transition temperatures.23,24 Additionally, antisite disorder, where Co and Mn cations interchange positions, introduces complexities such as spin-glass states and secondary magnetic transitions due to competing magnetic interactions (Co2+–Co2+ or Mn4+–Mn4+). Such disorder, along with antiphase boundaries, can profoundly alter the material's magnetic and electronic behaviour.25–28
Research on analogous rare-earth double perovskites underscores the critical role of antisite disorder in governing structural, electronic, and magnetic properties. Notable studies on La2CoMnO6,29,30 Sm2CoMnO6,11,18,31 Nd2CoMnO611,32 and Tb2CoMnO633 reveal a spectrum of fascinating phenomena, including exchange bias effects, Griffiths phases, and spin–phonon coupling, emphasizing the intricate interplay between cationic ordering and physical properties. Nayak et al.31 examined Sm2CoMnO6 (SCMO) polycrystals and found both ordered (monoclinic P21/n) and disordered (orthorhombic Pnma) phases. They observed an exchange bias effect of 10 kOe at 2 K, a spontaneous exchange bias below 20 K, and a ferromagnetic ground state with a 1.3 eV band gap.
Further advancing the field, Aswathi et al.34 investigated Pr2FeMnO6, uncovering a negative exchange bias below 200 K due to the Dzyaloshinskii–Moriya (DM) interaction and magnetic exchange anisotropy. Their thermomagnetic analysis revealed a high-temperature magnetic transition at Tc = 576 K, with the presence of ferromagnetic short-range correlations above Tc, suggesting a Griffiths-like phase. Das et al.29 explored Nd2CoMnO6, identifying a Griffiths phase, and Bhatti et al.18 investigated Sm2CoMnO6, revealing spin–phonon coupling, a high dielectric constant, and frequency-dependent dielectric response driven by electron hopping.
Dy2CoMnO6 (DCMO), in particular, exhibits complex phenomena such as Griffiths phase, mixed valence states, and short-range magnetic ordering. Xin et al.35 reported that DCMO crystallizes in a monoclinic P21/n structure and hosts a Griffiths phase, unlike Gd2CoMnO6, due to distinct 3d–4f exchange interactions and antisite disorder disrupting long-range ferromagnetic order. Chanda et al.36 demonstrated excellent photocatalytic efficiency of DCMO nanoparticles with a bandgap of 1.93 eV and p-type polaronic conduction. Pedro et al.37 revealed ferrimagnetic ordering accompanied by structural phase coexistence, whereas Das et al.38 observed ferroelectric-paraelectric transitions, polaronic conduction following the OLPT model, and an antiferromagnetic ground state. Muddelwar et al.39 further explored DCMO as a supercapacitor electrode, revealing significant electrochemical performance. Structural studies showed CoO6 and MnO6 octahedral tilting and electron density distortions mapped using GFourier analysis. The micron-sized flat particles provided faster electrochemical reactions and shorter ion/electron diffusion paths. Electrochemical characterization indicated a specific capacitance of 100.9 F g−1 at 1 A g−1 with remarkable cyclic stability of 87% over 10
000 cycles, suggesting its potential as a next-generation supercapacitor electrode material. Further advancing this discourse, Valian et al.40 synthesized urchin-like DCMO nanostructures via a modified Pechini route, employing novel combinations of chelating and cross-linking agents to regulate particle morphology. The resultant nanostructures displayed superior photocatalytic performance, achieving 88% methylene blue dye degradation under ultraviolet irradiation within 120 minutes, underscoring their efficacy as photocatalytic materials. Essaoud et al.41 utilized density functional theory (DFT) to probe the electronic, magnetic, and thermoelectric properties of DCMO, revealing its ferrimagnetic semiconducting behavior with spin-dependent indirect bandgaps of 0.31 eV (minority spin) and 1.65 eV (majority spin). The mixed covalent-ionic bonding characteristic of Dy–O, Co–O, and Mn–O was highlighted by topological analysis using quantum theory of atoms in molecules (QTAIM). At higher temperatures (ZT > 0.5 above 500 K), when p-type hole doping outperformed electron doping, the material demonstrated encouraging thermoelectric performance. Notably, DCMO demonstrated exceptional magnetic robustness, with dominant dysprosium contributions, maintaining structural and magnetic stability under hydrostatic pressures up to 35 GPa. Anirban et al.42 studied Dy2CoMnO6, revealing a monoclinic phase with structural distortions. Charge transport showed activation energies of 0.16 eV (≤75 °C) and 0.37 eV (≥100 °C). Dielectric properties exhibited non-Debye behavior, and the leakage current density followed Ohmic conduction, described by the Schottky barrier model.
In this work, we undertake a comprehensive investigation of the structural, electronic, and magnetic properties of DCMO, with particular emphasis on elucidating the contributions of antisite disorder and 3d–4f exchange interactions to its functional behavior. By bridging the domains of structural modulation and physical phenomena, this study seeks to advance the understanding of structure–property correlations and broaden the technological applicability of DCMO despite significant progress in understanding Dy2CoMnO6, a comprehensive understanding of antisite disorder's impact on its structural, electronic, and magnetic properties remains elusive. Current research often separates experimental observations from theoretical models, resulting in fragmented insights into long-range magnetic ordering, charge transport, and electronic interactions. This study combines experimental techniques and DFT calculations to address these gaps, focusing on the effects of antisite disorder. By linking experimental data with theoretical predictions, it aims to deepen the understanding of rare-earth double perovskites, paving the way for their integration into advanced multifunctional devices.
• Characterization details – Post-synthesis, an exhaustive characterization of the polycrystalline Dy2CoMnO6 double perovskite was undertaken utilizing a suite of sophisticated analytical techniques. A Bruker D8 Advance apparatus was used to perform the powdered X-ray diffraction (XRD) investigation using CuKα radiation (λ = 1.54056 Å). Subsequent Rietveld refinement of the diffraction patterns was meticulously executed using the FullProf Suite, where the pseudo-Voigt function facilitated peak profile fitting, leveraging crystallographic input data sourced from relevant literature. Using the KBr pellet approach, Fourier-transform infrared (FT-IR) absorption spectra were obtained using a PerkinElmer FT-IR spectrometer (Spectrum 1000, Japan). Using a 633 nm He–Ne laser running at full intensity, Raman spectroscopy was carried out using a LaBRAM HR spectrometer (Horiba, France SAS). X-Ray photoelectron spectroscopy (XPS) measurements were carried out at the Indus-2 synchrotron facility, RRCAT Indore, employing photoemission electron spectroscopy (PES) on beamline BL-14. In an ultra-high vacuum environment with a pressure of 5 × 10−9 mbar, the experimental setup included a bending magnet source of 1.5 T, a Si(111) double-crystal monochromator, and a hemispherical analyzer/detector system (Phoibos 225, Specs manufacture). Magnetization measurements were executed using a Quantum Design SQUID VSM, spanning the temperature range of 10 ≤ T ≤ 350 K and subjected to a magnetic field of H = 500 Oe. For low-field magnetization assessments, the reset magnet mode was employed to negate any residual or stray field effects. Ultraviolet-visible (UV-Vis) absorption spectra were recorded with a CARY Varian-300 spectrophotometer, ensuring a comprehensive elucidation of the material's optical characteristics.
• Computational details: The Vienna Ab initio Simulation Package (VASP) was used to conduct geometrical optimization and electronic calculation in this study within the context of density-functional theory (DFT).43 Both electron–electron interactions and ionic core effects were precisely incorporated using the projector-augmented wave (PAW) approach. The calculations for the exchange–correlation potential were carried out using the Perdew–Burke–Ernzerhof (PBE) functional within the generalized gradient approximation (GGA). To better handle the localization of the d-electron states in the Co and Mn atoms, we applied a Hubbard correction (DFT+U), setting the effective U values at 4.80 eV for Co and 4.50 eV for Mn.44 To guarantee accurate geometrical and electrical properties, a kinetic cut-off energy of 500 eV was established. We also conducted band structure calculations that included spin–orbit coupling (SOC), applying the DFT+U method to precisely capture electron interactions and localization. For optimisation and partial density of states analysis, a dense 10 × 8 × 6 k-point mesh created with the Monkhorst–Pack (MP) technique was also used. Structure relaxation is performed using the conjugate gradient algorithm until the total energy converges to within 10−4 eV and the Hellmann–Feynman force on each atom is less than 0.003 eV A−1.
![]() | ||
| Fig. 1 (a) Rietveld refinement of powder X-ray diffraction pattern of Dy2CoMnO6; (b) enlarged view of the selected 2θ range (32.25°–34.50°) highlighting the fitting accuracy in that region. | ||
| Dy2CoMnO6 | |||||
|---|---|---|---|---|---|
| Space group | P21/n (14) | ||||
| Lattice parameters | |||||
| a(Å) | b (Å) | c (Å) | α | β | γ |
| 5.2540 | 5.5755 | 7.4902 | 90° | 90.0126° | 90° |
| Volume (Å3) | 219.415504 | ||||
| Atomic positions | |||||
|---|---|---|---|---|---|
| Atom | Site | x | y | z | Occupancy |
| Dy | 4e | −0.01200 | 0.07276 | 0.25208 | 1.0 |
| Co | 2c | 0.00000 | 0.50000 | 0.00000 | 0.5 |
| Mn | 2d | 0.50000 | 0.00000 | 0.00000 | 0.5 |
| O1 | 4e | 0.08839 | 0.47656 | 0.27729 | 1.0 |
| O2 | 4e | 0.62803 | 0.28493 | 0.06634 | 1.0 |
| O3 | 4e | 0.72874 | 0.26667 | 0.45298 | 1.0 |
| Reliability factors | |||||
|---|---|---|---|---|---|
| R p | R wp | R exp | R Bragg | R F | χ 2 |
| 26.2 | 20.7 | 17.97 | 5.694 | 5.984 | 1.32 |
| Bond distance (Å) | |||||
|---|---|---|---|---|---|
| (Co–O1) × 2 | (Co–O2) × 2 | (Co–O3) × 2 | (Mn–O1) × 2 | (Mn–O2) × 2 | (Mn–O3) × 2 |
| 2.13215 | 2.34620 | 1.94401 | 1.73660 | 1.79527 | 1.96149 |
| Bond angle | |||||
|---|---|---|---|---|---|
| Mn–O1–Co | 150.793° | Mn–O2–Co | 134.888° | Mn–O3–Co | 157.506° |
Fig. 2 illustrates the schematic illustration of the Dy2CoMnO6 unit cell, depicting the distribution of ions across crystallographic positions. As per Table 1, these positions correspond to 4e for Dy3+ ions, 2c for Co2+ ions, 2d for Mn4+ ions, and 4e for O2− ions. Six O2− ions encircle each Co2+ and Mn4+ ion to create CoO6 and MnO6 octahedra, respectively.5 The BO6 octahedra tilting in the direction of the b-axis is described by the Glazer system a−b+a−. The schematic diagram shows tilt in both CoO6 and MnO6 octahedra. To determine the cations’ valence states quantitatively, bond valence sums (BVS) were calculated. The final structural parameters are shown in Table 1, which also includes important bond distances and angles related to CoO6 and MnO6 octahedra, as well as BVS values for Dy2CoMnO6.
![]() | ||
| Fig. 2 (a) The crystal structure with monoclinic space group of Dy2CoMnO6 double perovskite along c-axis, (b) and (c) visual representation of bond angles between CoO6 and MnO6 octahedra. | ||
The refined structural parameters presented in Table 1, along with the monoclinic crystal structure depicted in Fig. 2, clearly illustrate that the bond lengths and bond angles, such as Co–O and Mn–O, and Co/Mn–O–Mn/Co, respectively, do not exhibit identical values, with Co–O/Mn–O bond lengths and Co/Mn–O–Mn/Co bond angles deviating from 180°.29 This finding strongly implies that the Dy2CoMnO6 (DCMO) system exhibits a distorted crystal structure. Furthermore, the Co–O/Mn–O bond lengths and Co/Mn–O–Mn/Co bond angles in DCMO differ from those observed in an ordered double Perovskite (DP) compound at room temperature. In disordered systems, it is well known that disorder and deformed crystal structures in perovskites frequently affect magnetic contests and cause frustration.14,47,48 Thus, a similar situation of magnetic competition and frustration is anticipated in our disordered Dy2CoMnO6 system.
The Raman spectra, stemming from the monoclinic structure identified by space group P21/n (C2h5, no. 14), can be interpreted as a combination of in-phase (+) and antiphase (−) tilts of CoO6 and MnO6 octahedra, occurring along the [001] direction and within the basal plane of the pseudo-cubic cell.52–54 This aligns with the a−a−c+ notation according to Glazer's system. These vibrational modes observed in Raman spectroscopy predominantly reflect the dynamic characteristics of Co/MnO6 octahedra and Co/Mn–O bonds. Additionally, this monoclinic arrangement shows a direct group–subgroup relationship with tetragonal symmetry of I4/m. By employing factor group analysis, we ascertain the allocation of zone-center vibrational modes in relation to the corresponding representations of C2h5 point groups and monoclinic space group.55,56 The outcomes are outlined in Table 2.
| Atom | Site | Symmetry | Distribution |
|---|---|---|---|
| Dy | 4e | 1 | 3Ag + 3Au + 3Bg + 3Bu |
| Co | 2c | −1 | 3Au + 3Bu |
| Mn | 2d | −1 | 3Au + 3Bu |
| O1 | 4e | 1 | 3Ag + 3Au + 3Bg + 3Bu |
| O2 | 4e | 1 | 3Ag + 3Au + 3Bg + 3Bu |
| O3 | 4e | 1 | 3Ag + 3Au + 3Bg + 3Bu |
| Γ = 12A g + 18A u + 12B g + 18B u | |||
| Γ Raman = 12Ag + 12Bg | |||
| Γ Acoustic = Au + 2Bu | |||
| Γ IR = 17Au + 16Bu | |||
Hence, it is anticipated that the monoclinic Dy2CoMnO6 compound will exhibit 24 Raman-active vibrations, i.e.; 12Ag and 12Bg modes. Utilizing the factor-group analyses outlined above, a meticulous fitting of the Raman spectra for the Dy2CoMnO6 double perovskite compound was conducted. The outcomes are illustrated in Fig. 4(a)–(c), wherein the spectrum acquired at 300 K is segmented into three wavenumber ranges: 60–300 cm−1 (Fig. 4(a)), 300–650 cm−1 (Fig. 4(b)), and 650–950 cm−1 (Fig. 4(c)) for enhanced clarity. The observed modes and their corresponding proposed assignments, which will be discussed further below, are documented in Table 3. Out of the anticipated 24 phonons for monoclinic compounds, we successfully detected and fitted a sum of 21 Lorentzian peaks [illustrated by the solid lines in Fig. 4(a)–(c)], as only a few peaks with very low intensity were disregarded during the Lorentzian fitting process. The peak positions and full-width at half-maximum (FWHM) values are detailed in Table 3.
![]() | ||
| Fig. 4 Lorentzian fitted Raman spectra of Dy2CoMnO6 into three wavenumber ranges: (a) 60–300 cm−1, (b) 300–650 cm−1, and (c) 650–950 cm−1. | ||
| Band | Wavenumber (cm −1 ) | FWHM (cm −1 ) |
|---|---|---|
| 1 | 96.4568 | 19.34464 |
| 2 | 105.3659 | 7.90978 |
| 3 | 118.3337 | 27.41265 |
| 4 | 134.258 | 21.6095 |
| 5 | 152.5968 | 6.34473 |
| 6 | 197.5378 | 4.90579 |
| 8 | 261.193 | 14.39484 |
| 9 | 281.8333 | 1.59602 |
| 10 | 349.6054 | 12.33295 |
| 11 | 374.234 | 21.79637 |
| 12 | 432.6127 | 20.40047 |
| 13 | 482.8098 | 7.27032 |
| 14 | 495.4944 | 10.77868 |
| 15 | 522.7165 | 8.58732 |
| 16 | 565.8443 | 19.57891 |
| 17 | 620.1975 | 10.97171 |
| 18 | 619.9716 | 9.44831 |
| 19 | 689.9308 | 11.72895 |
| 20 | 728.7796 | 155.1789 |
| 21 | 927.2323 | 27.29348 |
The Raman representation corresponding to the monoclinic (P21/n) structure of Dy2CoMnO6 can be expressed as given in eqn (i):
| Γ = T(3Ag + 3Bg) + L(3Ag + 3Bg) + ν1(Ag + Bg) + ν2(2Ag + 2Bg) + ν5(3Ag + 3Bg) | (i) |
In this context, L and T denote the librational and translational phonon modes, respectively. The notation ν1 refers to the totally symmetric stretching, ν2 denotes the asymmetric stretching, and ν5 represents the symmetric bending. Within the monoclinic lattice framework, the unit cell accommodates two formula units of Dy2CoMnO6, giving rise to the manifestation of Ag and Bg symmetry components across all vibrational modes. It is reasonable that these symmetry components, particularly of the internal vibrational modes, converge at nearly identical frequencies, necessitating the employment of polarized Raman spectroscopy on single crystals to achieve their resolution.
In contrast, for polycrystalline sample, as examined in the present study, the inherent averaging of crystallographic orientations results in a simplification of the observed spectral features. Consequently, the Raman spectra are anticipated to display a solitary band corresponding to the symmetric stretch (ν1), a doublet associated with the asymmetric stretch (ν2), and a triplet attributed to the bending mode (ν5). The L and T vibrational modes encompass the frequency range of 90–300 cm−1, predominantly arising from the dynamics of the (Co/Mn)O6 octahedra and Dy–O bonds.50,57,58 The translational (T) modes within this interval are attributed to the displacement of the Dy cation, with specific lattice modes identified at 118.33, 134.25, 152.59, 197.53, and 261.19 cm−1 peaks located within the 320–350 cm−1 and 430–440 cm−1 ranges correspond to librational (L) lattice vibrations, which stem from the rotational dynamics of the Dy cation. Notably, librational modes are discernible at 281.83 cm−1, 349.60 cm−1, 374.23 cm−1, and 432.61 cm−1. For polycrystalline samples, a rigorous and exhaustive assignment of all experimentally observed modes via polarization analysis is precluded due to orientation averaging. However, important modes related to the MnO6 octahedra can be clarified by comparing the vibrational modes of the prototypical Fm3m cubic structure with the P21/n monoclinic structure.
The peak observed at 728 cm−1 is assigned to the Ag mode, representing the totally symmetric stretching of the MnO6 octahedra (ν1), wherein the Dy atom remains immobile as the oxygen atoms oscillate along the Co–O–Mn axis. Asymmetric oxygen stretching vibrations, with the Dy atom held stationary, are designated as the ν2 (Bg) mode, occurring at 565.84 cm−1 and 689 cm−1.59 The bending dynamics of oxygen within the octahedra, classified under the ν5 mode, manifest at 482.80 cm−1 and 495.49 cm−1, with their frequencies modulated by the chemical nature of the Co ions.60,61 Additionally, the mode at 927 cm−1 is identified as the ν1 mode, signifying the totally symmetric stretching of oxygen atoms along the Co–O–Mn axis within the MnO6 octahedra.58,62,63
O double bond absorbs infrared light at wavenumbers between approximately 1600–1900 cm−1.66 The band at 3433 cm−1 occurs due to the O–H vibrations.66
X-Ray photoelectron spectroscopy (XPS) emerges as a particularly effective method for scrutinizing the valence states and ligand coordination of constituent elements within a material. In materials containing open-shell ions, the interaction of a core electron vacancy with the open-shell leads to the development of numerous structures evident in XPS analysis. To deduce oxidation states and ligand coordination accurately, one must consider not only the primary photoelectron peak features but also related satellite peaks, chemical shifts, and their corresponding intensities. Understanding various physical characteristics of a material hinges upon grasping its electronic structure. Thus, we have thoroughly examined the XPS spectra of the current material, DCMO, to glean insights into its electronic structure. All peak positions have been assigned based on the National Institute of Standards and Technology (NIST) database. The deconvolution analysis of core-level XPS peaks for relevant ions was performed using a combination of Lorentzian and Gaussian distribution functions. The survey scan XPS spectra recorded at 300 K, depicted in Fig. 6, confirm the presence of Dy, Co, Mn, O, and C elements within the system.35 The absence of extrinsic elements confirms the sample's purity, while the detection of the C 1s peak, typical in such analyses, is attributed to molecules from the surrounding air absorbed onto the surface.
High-resolution XPS spectra of Dy-3d, Co-2p, Mn-2p, and O-1s core levels for the as-prepared sample are depicted in Fig. 7(a)–(d), respectively. Due to the spin–orbit interaction, each spectrum is split into two components. One component has a lower J (total angular momentum) value, corresponding to the (l − s) coupling, and it appears on the higher energy side. The other component has a higher J value, corresponding to the (l + s) coupling, and it appears on the lower energy side. The XPS spectra of the Dy 3d core level, along with the fitting of peaks, are depicted in Fig. 7(a). It's apparent from the figure that the two distinct spin-orbital split peaks, Dy 3d5/2 and Dy 3d3/2, are situated at 1297.77 eV and 1335.5 eV, respectively, with a spin-orbital splitting energy of 37.73 eV. Detailed analysis of the XPS spectrum indicates that the Dy cations are present in the +3 oxidation state, consistent with previous reports in the literature.38Fig. 7(b) displays the Co-2p core level XPS spectrum of the Dy2CoMnO6 double perovskite compound at room temperature. The spectrum exhibits two prominent peaks at 782 and 797 eV corresponding to Co-2p3/2 and Co-2p1/2, respectively, due spin orbital coupling.28 Deconvolution of the 2p3/2 peak reveals two distinct peaks corresponding to Co2+ 2p3/2 and Co3+ 2p3/2 at 785.186 and 781.83 eV, respectively. Similarly, the Co 2p1/2 peak can be deconvoluted into Co2+ 2p1/2 and Co3+ 2p1/2 peaks at 798.744 and 797.019 eV, respectively.48,55 In cobalt-based systems, a change of approximately 1 electron volt (eV) in binding energy signifies a corresponding alteration in the valence state of cobalt. Consequently, when employing XPSPEAK 4.1 software for fitting, two peaks are typically considered for each spin–orbit coupling, along with their corresponding charge transfer satellites. Multiple electron excitation of Co2+ and Co3+ in the cobalt species results in the appearance of two shake-up satellite peaks at 790.025 and 804.240 eV. Therefore, the presence of satellite peaks in the Co 2p XPS spectrum suggests the existence of Co2+ ions in the DCMO compound. It is worth mentioning here that the asymmetry and broadening observed in the peaks indicate the presence of mixed valence states of the Co ions. The peak positions and line-shape observed in the Co 2p XPS spectrum are similar to earlier reports, indicating the presence of mixed valence Co ions. For the present system, the observed Co 2p doublet separation of 15 eV suggests the coexistence of both Co2+ and Co3+ ions. The Co 2p XPS spectrum has been deconvoluted to estimate the concentrations of different Co ions (Co2+/Co3+), as illustrated in Fig. 7(b). While Co3+ ions are predominantly present in the system, the presence of Co2+ ions are also unavoidable.5,33 This, in turn, can lead to competing exchange interactions in Dy2CoMnO6.
![]() | ||
| Fig. 7 High-resolution XPS spectra of (a) Dy-3d, (b) Co-2p, (c) Mn-2p, and (d) O-1s core levels Dy2CoMnO6. | ||
The core level Mn 2p XPS spectrum of the present material is illustrated in Fig. 7(c). This spectrum is broadly divided into two spin–orbit coupling peaks, namely Mn 2p3/2 and Mn 2p1/2, observed at 644 eV and 656 eV, respectively.18 These peaks result from the spin–orbit splitting of 2p orbitals, with a splitting energy of 12 eV. The peaks observed near ∼643.38 and 654.76 eV suggest the presence of Mn3+ ions, whereas peaks near ∼645.75 and 656.79 eV indicate Mn4+ charge states. As a result, XPS analysis demonstrates the presence of Co2+, Co3+, Mn3+, and Mn4+ charge states for the B-site cations.56,63 The observed magnitudes of these charge states are provided in Table 4. Comparing the chemical formula of the investigated compound with the selected formula, Dy2CoMnO6, using the content of Co and Mn verifies the formula. Furthermore, the presence of Mn3+ and Co3+ confirms our hypothesis regarding the Jahn–Teller (JT) distortion in the system. The O 1s core level XPS spectrum is depicted in Fig. 7(d). The peak at 530.467 eV corresponds to the distinctive peak of O2− ions, indicating the bonding between metal and oxygen (M–O–M). Meanwhile, the peak at 531.976 eV is associated with defect sites having low oxygen coordination. The broad signal at 532.625 eV suggests that oxygen ions are ionized and associated with the subsurface of the Dy2CoMnO6 ceramic through physisorption or chemisorption of oxygen molecules.33 This implies that the subsurface is composed of oxygen ions, which may be classified as excess oxygen or O− species, given their lower electron concentration compared to O2− ions. The presence of mixed states of Co and Mn elements suggests the absence of full long-range ordering between Co and Mn elements in the Dy2CoMnO6 double perovskite.68
| Atom | Orbital | Binding energy | FWHM | Area | Peak specification |
|---|---|---|---|---|---|
| Dy | 3d | 1297.77 | 8.364 | 254 824.90 |
3d5/2 |
| 1305.22 | 7.927 | 21 698.55 |
Satellite | ||
| 1335.5 | 6.574 | 217 812.6 |
3d3/2 | ||
| 1339.58 | 3.932 | 10 180.71 |
Satellite | ||
| Co | 2p 3/2 | 781.830 | 3.974 | 7991.513 | Co3+ |
| 785.186 | 5.463 | 7180.572 | Co2+ | ||
| 790.025 | 4.373 | 2422.829 | Satellite | ||
| 2p 1/2 | 797.019 | 3.678 | 2944.458 | Co3+ | |
| 798.744 | 4.095 | 3240.034 | Co2+ | ||
| 804.240 | 5.924 | 4204.529 | Satellite | ||
| Mn | 2p 3/2 | 643.379 | 4.052 | 5074.523 | Mn4+ |
| 645.748 | 5.624 | 6495.841 | Mn3+ | ||
| 2p 1/2 | 654.759 | 4.419 | 2286.352 | Mn4+ | |
| 656.791 | 5.233 | 3580.431 | Mn3+ | ||
| O | 1s | 530.467 | 2.556 | 11 910.12 |
Lattice oxygen |
| 531.976 | 3.113 | 12 448.57 |
O2− vacancy | ||
| 532.625 | 6.058 | 48 062.01 |
Chemisorbed oxygen |
| Cationic charge states observed in XPS analysis | |||
|---|---|---|---|
| Cations/valence state | 2 + | 3 + | 4 + |
| Co (valence state %) | 47.33 | 52.67 | |
| Mn (valence state %) | 56.14 | 43.86 | |
The presence of Co3+ and Mn3+ ions confirm our hypothesis regarding the system's Jahn–Teller (JT) distortion. Recent studies have shown that JT distortion can be induced by Co2+ ions in the high-spin (HS) state within an octahedral (Oh) unit. Additionally, the presence of Mn ions in the 3+ and 4+ states can also contribute to JT distortion in the ceramic. This disordering of cationic charges results in non-collinear magnetic behavior in the ceramics. Further discussion on field and temperature-dependent measurements of Dy2CoMnO6 will be presented in the subsequent section.
To identify the most stable magnetic phase between antiferromagnetic and ferromagnetic configurations, energy–volume curves were plotted (Fig. 8). Dy2CoMnO6 is energetically more stable in the ferromagnetic phase, as the curves clearly demonstrate that the ferromagnetic state exhibits a lower energy compared to the antiferromagnetic state.69
To ascertain the magnetic ground state, we systematically examined the AFM configuration along with multiple spin-alignment variants of ferrimagnetic (FiM) states. The computed relative energies are presented in Table 5. A comprehensive analysis of these calculations revealed that the FM configuration exhibits the lowest energy, in agreement with magnetization measurements. This finding indicates that the predominant exchange interaction between the Co and Mn atoms in DCMO is ferromagnetic in nature.
| Magnetic ordering | Spin alignment | Relative energy (eV) | |||
|---|---|---|---|---|---|
| Co1 | Co2 | Mn1 | Mn2 | ||
| FM | ↑ | ↑ | ↑ | ↑ | 0.00 |
| AFM | ↑ | ↓ | ↑ | ↓ | 0.10 |
| Ferri-1 | ↑ | ↓ | ↑ | ↑ | 1.27 |
| Ferri-2 | ↓ | ↓ | ↑ | ↓ | 0.91 |
| Ferri-3 | ↓ | ↓ | ↓ | ↑ | 0.90 |
| Ferri-4 | ↓ | ↑ | ↑ | ↑ | 1.27 |
| Ferri-5 | ↑ | ↓ | ↓ | ↓ | 0.19 |
| Ferri-6 | ↑ | ↑ | ↓ | ↑ | 1.23 |
Quantum exchange interactions, which are connected to atoms’ magnetic moments, distances from one another, and any applied external magnetic field, are what cause the microscopic magnetic behaviour in materials. The Heisenberg effective spin Hamiltonian can be used to model these interactions:
![]() | (ii) |
represents the spin operator,
is external magnetic field, and Jij is exchange coupling constant, which depends on the distance amongst atoms.
The total and partial magnetic moments of the Dy2CoMnO6 were calculated in its most stable state, as depicted in Fig. 8 and Table 5. The results indicate that the molecule has a significant net magnetic moment of 11.843μB. Dysprosium (Dy) delivers 0.076μB to the unit cell's magnetic moment, according to an analysis of the contributions from each of the atoms. Manganese (Mn) contributes significantly with 6.305μB, followed by cobalt (Co) with 5.324μB, playing a key role in the perovskite structure. This large overall magnetic moment likely arises from the overlap between Co-3d and Mn-3d orbitals. Oxygen atoms contribute 0.138μB, further supporting the magnetic properties. Additionally, superexchange interactions between Mn4+ and Co2+ ions, mediated by oxygen, involve the half-filled Mn-3d-eg and empty Co-3d-eg orbitals, enhancing the magnetic moment.
The dc magnetization of the DCMO ceramic sample was recorded across a temperature range of 10 to 350 K, employing both zero-field-cooled (ZFC) and field-cooled (FC) protocols under an applied field of 500 Oe. As depicted in Fig. 9(a), the magnetization versus temperature (M–T) curves exhibit pronounced irreversibility between the ZFC and FC processes, along with a notable peak in the ZFC curve. The observed peak at approximately 80 K is likely indicative of the localized freezing of magnetic moments, reflecting the characteristics of a conventional spin-glass state, where the magnetic moments freeze randomly at lower temperatures.9,13,34 A typical transition from paramagnetic (PM) to ferromagnetic (FM) behavior, occurring around the transition temperature (Tc) of 87 K as shown in Fig. 9(b), is attributed to the ferromagnetic superexchange interactions between Co2+ (t52g,e2g) and Mn4+ (t32g,e0g) ions,16,32 with Tc being determined from the minimum of the first derivative of the magnetization in the zero-field-cooled (ZFC) mode, as illustrated in the inset of Fig. 9(a).
In DCMO, the magnetic properties are influenced by the contributions from both rare-earth Dy3+ (4f10, 5I8; J = 8) ions and transition metal ions (Co2+ and Mn4+). From the Curie–Weiss (CW) fit, expressed as
[eqn (iii)]; where θ is the Curie–Weiss temperature and C is the Curie constant, the values of C and θ are determined to be 33.37 emu K Oe−1 mol−1 and 18.77 K, respectively. The positive value of θ indicates the presence of ferromagnetic (FM) interactions between Dy3+ moments and the transition metal ion sublattice (Co2+/Mn4+), which may also be associated with the crystal electric field excitations of the 5I8 multiplet of the Dy3+ ion.5
Fig. 10(a) exhibits the variation of inverse susceptibility with temperature for DCMO compound. A down turn can be noticed in the magnetic susceptibility below Tc. There is a broaden peak around ∼30 K, which is relevant to the canted – AFM correlations below Tc. The emergence of a decrease in magnetic susceptibility below Tc is a clear indication of the existence of the Griffiths phase (GP). The Griffiths temperature TG = 91 K is defined as the maximum of the d(1/χ/)/dT vs. temperature plot as shown in the inset of Fig. 10(b).
Between Tc and TG temperatures, it is evident that the 1/χ − (T) curve doesn’t exhibit a linear trend, suggesting a deviation from the Curie–Weiss (C–W) law for the inverse susceptibility within this temperature range. However, the 1/χ data above TG can be accurately modelled using the C–W law, signifying a characteristic indicative of paramagnetic behavior. The effective magnetic moment of the compound was evaluated experimentally in units of Bohr magneton (μB) by using the following equation:70
![]() | (iii) |
![]() | (iv) |
As depicted in Fig. 10(a), the inverse susceptibility follows the Curie–Weiss law at higher temperatures; however, it deviates from linear fitting below ∼100 K. The detected deviation is the characteristic feature of the GP, which originated because of the presence of short-range ferromagnetic clusters in the paramagnetic region. The GP is typically described by the following relation:72,73
| χ−1 = (T − Trdc)1−λ | (v) |
![]() | ||
| Fig. 11 Inverse susceptibility as a function of temperature under 500 Oe applied magnetic field. Solid red line shows the fitting by Griffith model. | ||
Magnetization (M) as a function of applied magnetic field (H) recorded at 10 K and 350 K is presented in Fig. 12(a). As depicted in Fig. 12(a), the magnetization experiences a rapid increase at lower magnetic fields, and at higher fields, saturation is not achieved, even when the field reaches ±50 kOe. The observed hysteresis, along with the absence of saturation at higher fields, may be because of the presence of antiferromagnetic (AFM) correlations and the coexistence of competing magnetic orders in both phases. Hysterias loop at 10 K shows a coercivity of 6.8 kOe also suggests the FM nature of the compound. The hysteresis loop recorded at 350 K exhibits a linear pattern in response to the applied magnetic field, lacking any indication of magnetization saturation. This behavior aligns with expectations for the paramagnetic region of a FM/AFM system.
In Fig. 12(b), we have presented Arrott plot extracted from virgin curve of M–H data shown in Fig. 12(a). Arrott plot (M2vs. H/M) offers an effective tool to understand the nature of magnetic state. Fig. 12(b) displays a negative intercept, indicating the absence of spontaneous magnetization in this material at 10 K. The absence of spontaneous magnetization confirms that the compound is not pure ferromagnetic.
| αhν = A(hν − Eg)n | (vi) |
Here, n is either ½ or 2 for direct or indirect transitions, hν is the photon energy, α denotes the absorption coefficient, A represents a constant that is independent of photon energy, and Eg is the optical band gap. As illustrated in in Fig. 13(b), visualisations of (αhν)2 and (αhν)1/2vs. hν were plotted based on this relationship.
The optical band gap values were determined by extrapolating the linear region of the absorption edge to the point where (αhν)2 = 0 and (αhν)1/2 = 0. The direct band gap was found to be approximately 1.47 eV, while the indirect band gap was approximately 1.04 eV. In comparison, Chanda et al. reported a band gap of 1.93 eV for Dy2CoMnO6 nanoparticles synthesized via the sol–gel method.36 This value, located within the visible spectrum, highlights the material's suitability for visible-light-driven applications.
The observed direct band gap in this study is slightly lower than previously reported values, potentially due to the particle size and polycrystalline nature of the synthesized Dy2CoMnO6. The reduced band gap enhances its potential as a visible-light photocatalyst, emphasizing its applicability in photocatalytic processes.
To further understand the optical properties of Dy2CoMnO6, first-principles density functional theory (DFT) calculations were performed to obtain the absorption coefficient as a function of photon energy. The computed absorption spectrum, presented in Fig. S1 of the ESI,† reveals two distinct optical band gaps at approximately 0.9 eV and 2.61 eV. These values were compared with the experimental results obtained from UV-Vis spectroscopy, which indicated an indirect band gap of 1.04 eV and a direct band gap of 1.47 eV, with a broad absorption centered around 2.21 eV due to charge transfer transitions between Mn(3d) and O(2p) states. The lower-energy DFT band gap (∼0.9 eV) aligns well with the experimentally determined indirect band gap (1.04 eV), suggesting the presence of intermediate electronic states or defect-related transitions. However, the higher-energy DFT band gap (2.61 eV) is larger than the experimental direct band gap (1.47 eV) and the reported literature value of 1.93 eV. This discrepancy may arise due to the limitations of the exchange–correlation functional used in the calculations, as generalized gradient approximation (GGA) functionals are known to overestimate or underestimate band gaps depending on the material. Additionally, factors such as the polycrystalline nature, particle size effects, and synthesis conditions of the experimental sample could contribute to the observed deviation.76,77 Despite these variations, the experimental and theoretical results consistently highlight Dy2CoMnO6 as a promising material for visible-light-driven applications.
Dy2CoMnO6 exhibits semiconductor behavior, characterized by an indirect band gap. For spin-up electrons, the band gap spans from the Z to D points in the Brillouin zone, with a value of 1.6054 eV. In the spin-down channel, the band gap is larger, measured at 3.2029 eV, corresponding to a transition between the Γ and Z points. The overall band gap, considering the spin-up channel, is 1.6054 eV. This spin-dependent band gap structure underlines the compound's complex electronic behavior, with the differing gaps for spin-up and spin-down electrons contributing to its magnetic and electronic properties. Essaoud et al.41 meticulously analyzed the electronic properties, focusing on the total and partial density of states, utilizing the advanced full-potential linearized augmented plane wave (FP-LAPW) method within the WIEN2k framework. Their results revealed a conspicuous absence of electronic states near the Fermi level, extending up to 1.65 eV for the minority spin channel and 0.31 eV for the majority spin channel, highlighting distinct spin-dependent electronic characteristics.
The total and projected density of states (TDOS and PDOS) of Dy2CoMnO6 have been analyzed to elucidate its electronic and magnetic properties. As shown in Fig. 15(a)–(f), the system exhibits a significant lack of states at the Fermi level in both spin channels, confirming its semiconducting nature. The computed energy gap is approximately 1.6054 eV and 3.2029 eV for the minority and majority spin orientations, respectively.
![]() | ||
| Fig. 15 Total and partial density of states (DOS) of Dy2CoMnO6 in spin up and spin down configuration. | ||
To understand the origin of magnetism in Dy2CoMnO6, we further examine the element- and orbital-resolved projected density of states (PDOS) contributions (Fig. 15(c)–(f)). The Co-3d and Mn-3d states exhibit strong spin polarization, with Co contributing 5.324μB and Mn contributing 6.305μB, indicating their dominant role in the magnetic ordering. The Co-3d orbitals show significant exchange splitting, where the spin-up and spin-down states are highly asymmetric, confirming a high-spin configuration for Co2+. Similarly, the Mn-3d states display notable exchange splitting, indicative of Mn4+ in a high-spin state. This spin asymmetry is a key factor contributing to the overall magnetization of the system.
In addition, the Dy-4f states are located deep in the valence band and do not significantly contribute to the overall magnetic moment, apart from minor hybridization effects, leading to a small Dy moment of 0.076μB. The oxygen atoms (O-2p states) display weak spin polarization, contributing 0.138μB, reinforcing the role of superexchange interactions in stabilizing the ferromagnetic state. The strong hybridization between O-2p and the transition metal 3d orbitals further supports the observed magnetic ordering.
The computed total magnetic moment of Dy2CoMnO6 is 11.843μB per formula unit, primarily arising from the cooperative exchange interactions between Mn4+ (6.305μB) and Co2+ (5.324μB) ions, mediated by oxygen ligands. The half-filled Mn-3d-eg and empty Co-3d-eg orbitals participate in a strong superexchange mechanism, stabilizing the ferromagnetic state.
We have employed the GGA+U approach to accurately account for electron–electron correlation in the localized Co-3d and Mn-3d orbitals. The applied Hubbard-U values (Ueff = 4.80 eV for Co and 4.50 eV for Mn) provide an improved description of their electronic and magnetic behavior. The PDOS analysis, along with the computed magnetic moments, confirms that the magnetism in Dy2CoMnO6 is predominantly driven by the exchange splitting of Co-3d and Mn-3d states, reinforced by oxygen-mediated superexchange interactions.
Spin–orbit coupling (SOC) in Dy2CoMnO6 (DCMO) significantly affects its electronic structure, magnetic anisotropy, and exchange interactions. The strong SOC associated with Dy leads to substantial splitting of the Dy 4f states, which influences their hybridization with the Co and Mn 3d orbitals. Additionally, SOC modifies the Co2+–O–Mn4+ superexchange interactions, altering the stability of different magnetic configurations. This effect induces magnetocrystalline anisotropy, impacting spin alignment and potentially stabilizing non-collinear magnetic structures. Furthermore, SOC enhances magnetoelectric coupling, contributing to complex spin–lattice interactions that can lead to spin-driven ferroelectricity in DCMO.
The electronic band structure of Dy2CoMnO6, shown in Fig. S2 of ESI,† is computed with SOC effects and reveals a direct band gap of 1.59 eV. The dispersion curves along high-symmetry points of the Brillouin zone (Γ–Z–D–B–Γ–A) illustrate the electronic transitions and the influence of SOC on band splitting. The separation of the valence band maximum (VBM) and conduction band minimum (CBM) confirms the semiconducting nature of DCMO. The inclusion of SOC modifies band dispersion, impacting the effective mass of charge carriers and the transport properties of the material. The electronic interactions between the Co and Mn 3d states and O 2p orbitals play a critical role in determining charge transfer and hybridization effects in the system.
The total and partial density of states (PDOS) calculations, depicted in Fig. S3 of ESI,† provide further insights into the electronic structure. The valence band is primarily composed of O 2p states hybridized with Mn 3d and Co 3d orbitals, whereas the conduction band is dominated by Dy 4f states. The strong localization of Dy 4f orbitals in the conduction band suggests a minimal contribution to electrical conductivity but plays a crucial role in magnetic exchange interactions. The presence of sharp peaks in the DOS indicates strong electronic correlations, particularly in the Dy 4f states, which are expected due to their localized nature. The hybridization of Dy 4f, Mn 3d, Co 3d, and O 2p states highlights the complex electronic structure of DCMO, where SOC significantly influences both electronic and magnetic properties.
The incorporation of SOC in electronic structure calculations reveals critical modifications in the band structure and density of states of Dy2CoMnO6. The SOC-induced band splitting, 3d–4f hybridization, and enhanced magnetic anisotropy play crucial roles in determining the transport and magnetic behavior of the material. These findings provide a deeper understanding of the fundamental electronic interactions in DCMO and suggest its potential application in spintronic and magnetoelectric devices, where SOC-driven effects can be utilized for advanced functionalities.
Zeta potential is a critical property of particles, reflecting the charges developed at the interface between the particle surface and the surrounding liquid medium.78 It is typically measured indirectly through electrophoretic mobility or particle velocity under an applied electric field. The relationship between zeta potential and surface charges on particles is complex and can be described by Henry's equation:
![]() | (vii) |
Zeta potential is influenced by factors such as particle size, shape, surfactant type and coverage, as well as the solvent's properties, including chain length, polarity, and dielectric constant. By adjusting the electrostatic interactions between particles, the zeta potential can be increased, leading to enhanced stability. A stable dispersion is generally indicated by a zeta potential value greater than ±30 mV. The size of the primary particles, surface coating, shape, material composition, and additives are some of the factors that affect the stability or aggregation of particle suspensions. Aggregation kinetics are controlled by the mobility of aggregates or the charge on particle surfaces; increased aggregation reduces particle mobility in the dispersion medium, which results in a lower zeta potential and, in turn, a less stable colloidal solution.
To examine the zeta potential (ζ) and hydrodynamic diameter (HD), we first optimized the solvent to achieve good stability, choosing deionized water (DI) as the dispersion medium for Dy2CoMnO6. The intensity distribution of HD and ζ, dispersed in DI, is shown in Fig. 16(a). In the first run, two peaks were observed in the HD data at 757 nm and 7288 nm, with the absolute value of ζ recorded at 243 mV and −29.68 mV in DIW (Fig. 16(b)). The HD results indicate that particles exhibit a narrower size distribution in ethanol, suggesting greater stability compared to water. To ensure measurement accuracy, the data were recorded three times. Interestingly, the two peak values of HD at 757 nm and 7288 nm in the first run shifted to three peaks at 352.28 nm, 1851.93 nm, and 10
497.63 nm in the second run. In the third run, the distribution returned to two peaks at 871.523 nm and 6113.69 nm. The average zeta potential was found to be −29.78 mV, with a mobility of −2.325 × 10−4 cm2 V−1 s−1 and an electric field of −16.37 V cm−1.
![]() | ||
| Fig. 16 (a) Intensity distribution of hydrodynamic diameter and (b) zeta potential in deionized water of Dy2CoMnO6. | ||
DC magnetization measurements in zero-field-cooled (ZFC) and field-cooled (FC) modes revealed a spin-glass state at ∼80 K and a ferromagnetic transition at 87 K. Curie–Weiss fitting confirmed ferromagnetic interactions between Dy3+ and Co2+/Mn4+, with crystal field effects from Dy3+. The inverse susceptibility curve indicated a Griffiths phase (GP) between Tc and TG, deviating from the Curie–Weiss law. The experimentally measured magnetic moment (16.33μB) exceeded the theoretical value (16.09μB), likely due to short-range ferromagnetic ordering within the paramagnetic matrix, possibly influenced by magnetic inhomogeneity or oxygen deficiencies. Studies on antisite disorder in Dy2CoMnO6 reveal that Co/Mn mixing on the B-site induces magnetic inhomogeneity and competing interactions, complicating the system's magnetic behavior. Magnetization measurements showed no saturation at high fields, indicating coexisting magnetic orders. The 10 K hysteresis loop exhibited 6.8 kOe coercivity, confirming ferromagnetism, while the 350 K loop was linear, indicating paramagnetism. These results highlight the role of antisite disorder in tuning the material's magnetic properties, making Dy2CoMnO6 promising for spintronic and magnetic applications.
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d5tc00299k |
| This journal is © The Royal Society of Chemistry 2025 |