Strain-mediated ferromagnetism and low-field magnetic reversal in Co doped monolayer \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$WS_2$$\end{document}WS2

Strain-mediated magnetism in 2D materials and dilute magnetic semiconductors hold multi-functional applications for future nano-electronics. Herein, First principles calculations are employed to study the influence of biaxial strain on the magnetic properties of Co-doped monolayer \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$WS_2$$\end{document}WS2. The non-magnetic \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$WS_2$$\end{document}WS2 shows ferromagnetic signature upon Co doping due to spin polarization, which is further improved at low compressive (-2 %) and tensile (+2 %) strains. From the PDOS and spin density analysis, the opposite magnetic ordering is found to be favourable under the application of compressive and tensile strains. The double exchange interaction and p-d hybridization mechanisms make Co-doped \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$WS_2$$\end{document}WS2 a potential host for magnetism. More importantly, the competition between exchange and crystal field splittings, i.e. (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta _{ex}>\Delta _{cfs}$$\end{document}Δex>Δcfs), of the Co-atom play pivotal roles in deciding the values of the magnetic moments under applied strain. Micromagnetic simulation reveals, the ferromagnetic behavior calculated from DFT exhibits low-field magnetic reversal (190 Oe). Moreover, the spins of Co-doped \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$WS_2$$\end{document}WS2 are slightly tilted from the easy axis orientations showing slanted ferromagnetic hysteresis loop. The ferromagnetic nature of Co-doped \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$WS_2$$\end{document}WS2 suppresses beyond \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pm 2~\%$$\end{document}±2% strain, which is reflected in terms of decrease in the coercivity in the micromagnetic simulation. The understanding of low-field magnetic reversal and spin orientations in Co-doped \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$WS_2$$\end{document}WS2 may pave the way for next-generation spintronics and straintronics applications.


Computational method
We perform the first principles spin-polarised DFT calculations using the Vienna ab initio simulation package (VASP) 33 which implements the projector augmented wave (PAW) method to describe electron-ion interaction. The electronic exchange-correlation potential is described by the generalized gradient approximation (GGA) in the Perdew-Burke-Ernzerhof (PBE) parametrization 27 . The correction of a Hubbard term U (GGA+U) is necessary to describe the strongly correlated 3d orbital of the TM impurities. However, the half-metallic doped magnetic TMDC systems are found to be independent of U values (U = 2.5, 3). The value of the magnetic moments remain little affected using different U values. Additionally, the optimization of lattice constants have been performed with various U values under unstrained and strained conditions by energy minimization and are observed to be unaltered with GGA+U (see SI). Henceforth, we find no relevant modifications on our conclusions and the further discussions are based on without considering the on-site interaction. We use kinetic energy cutoff of 500∼ eV for the plane-wave expansion of the wave functions. All the structures are fully relaxed under conjugate gradient algorithm until a total energy convergence and Hellmann Feynman force up to 10 −5 eV and 0.01 eV/ Å are achieved. The Brillouin zone integration for self-consistent and projected-density calculations is approximated by Monkhorst-Pack K point mesh of 5 × 5 × 1 and 12 × 12 × 1 respectively. We use 4 × 4 × 1 and 5 × 5 × 1 supercells for doping and strain-based calculations. The micromagnetic simulations is carried out by open freeware object-oriented-micromagnetic-framework (OOMMF) 34 . By considering the previous experimental reports 35,36 , the parameters such as magnetic anisotropy, saturation magnetization, exchange-length, etc are employed in Co-doped WS 2 by micromagnetic calculations.

Results and discussion
To study the effects of Co-doping for the magnetic properties of WS 2 mono layer (ML), we replace one host W-atom by foreign Co-atom in 4 × 4 × 1 [ (Co, W)S 2 ] and extended 5 × 5 × 1 super-cell with doping concentrations 6.25% and 4.16% respectively. The graphical representation of (Co, W)S 2 ML with its hexagonal unit cell (black solid line) is shown in Fig. 1a, where the optimized lattice constant is found to be 3.16 Å. The structural and magnetic parameters extracted from the calculated results are listed in Table 1. The obtained results show that the bond lengths of W-S shrink around the doping sites locally after Co-doping when compared to pristine WS 2 . However, the Co-S bond length remains same with W-S bond length in pristine WS 2 at various doping concentration. Such behavior may be attributed to various factors such as electron affinities, ionic radii differences between the Co and W ions after the geometrical relaxation and the changed electrostatic forces that arise due to the occurrence of multivalent Co ions at the host W site 37 . It is apparent to notice that WS 2 exhibits effective magnetic ordering with net magnetic moments 2.62 µ B and 2.58 µ B at different Co-doping concentration of 4.16% and 6.25%, respectively. The obtained results offer similar/higher magnetic moments than the previously  21 . The observation of such high magnetic moments in Co-doped WS 2 ML suggests its potential candidature for 2D DMSs. In order to investigate the effects of external perturbation such as strain on the magnetic properties, both compressive ( −5% to 0% ) and tensile ( 0% to 5% ) isotropic biaxial strain (along x-and y-axis simultaneously) have been applied to (Co, W)S 2 ML. From Fig. 1a, strained unit cells (dotted orange and violet cells for tensile and compressive strain, respectively) have been shown schematically along with the unstrained unit cell. It has been previously observed that the biaxial strain attributes isotropic changes to the magnitude of the lattice vectors of various TMDCs ML 38 . The relationship between the reciprocal lattice in Brillouin zone and the primitive vectors in real space is expressed as: where a i and b i are the primitive vectors and reciprocal lattice, respectively. Accordingly, under an isotropic/ uniform biaxial tensile (compressive) strain, the lattice constant will be increased (decreased), and the reciprocal lattice will be shrinked (enlarged), which in turn reduces (expands) the Brillouin zone as shown in Fig. 1b 38 . During our calculations, the periodicity is well preserved under the both biaxial strains. The structural stability of (Co, W)S 2 ML under isotropic biaxial strain can be estimated by formation energies, which can be expressed as: , where E pure and E doped represent the total energy of the pristine WS 2 and Co-doped WS 2 ML, respectively. µ W and µ Co are the chemical potential for host W and foreign Co-atom, respectively. n is the number of dopants in the studied supercell. The formation energies at various biaxial strains for W-and S-rich conditions are listed in Table 2. The formation energies of the strained systems are seen to increase monotonically with the increase in both compressive and tensile strains, as displayed in Fig. 2a. Previous experimental results evident that the S-rich condition is more suitable for the practical growth of pristine WS 2 ML than the W-rich condition 39 . For our system, it can be noticed that E form for S-rich is lower when compared to W-rich conditions at each studied strain system, as supported by the experimental condition. From Fig. 2b, it can be worth noting that the magnetic moments are further modified under various applied strains, which explore possible directions for TMDCs based spintronic and straintronic application 9,32 . Among all the applied strains, the maximum magnetic moment of 3.25 µ B and 2.69 µ B are achieved for +2% tensile and -2% compressive strain, respectively. The obtained results ensure a higher magnetic moment at low strain as compared to other TM doped TMDCs 30,31,40 . According to Luo et. al., the enhancement in the magnetic moments was only observed under compressive strain in Al-doped WS 2 ML 30 ; however, the present study reveals the improvement of the magnetic moment under both compressive and tensile strain by Co-doping. The enhanced magnetic properties under strained conditions may be attributed to the overlapping of spin-polarized electronic band alignments originating from both magnetic cations and nearest neighbour anions.       www.nature.com/scientificreports/ The pristine WS 2 ML exhibits a direct bandgap of E g = 1.87 eV (see SI: Fig S2) at the K-point having very close approximation with both experimental ( E g = 1.88 eV) 41 and other theoretical ( E g = 1.68 eV ) results 42 . Figure 3 elucidates the spin-polarized band structures of unstrained and strained (Co, W)S 2 ML at their equilibrium lattice constant. However, it can be noticed from Fig. 3b that some impurity states appear within the bandgap when a single W-atom is replaced by Co in both spin-up and spin-down channels. These impurity states are mainly contributed by Co-atom, while the effects from the nearest neighboring W and S atoms can be neglected, which has been explained in other TM doped TMDCs 21 . The asymmetric behavior of spin-up and spin-down components of Co-doped WS 2 ML indicates WS 2 behaves magnetically active after Co-doping. Moreover, the doped ML systems show half-metallic characteristics due to the suppression of bandgap. The spin-up component in the doped band structure shows dispersion behavior due to crossing of band lines across the Fermi level, as shown in Fig. 3b. Similarly, the spin-down component of the doped system has a bandgap of 1.168 eV with a half-metallic spacing of 0.294 eV for (Co, W)S 2 ML. Furthermore, the tuning of bandgap at various isotropic biaxial strains brings out effective variations in the magnetic properties of the doped system. The increase in the bond-lengths and bond-angles under various isotropic tensile strains are listed in Table 2. This indicates the decrease in the bond-energy between TM elements and S atoms leads to increase the number of spin-polarized electrons at the vicinity of fermi-level. As a result of which, the maximum obtained magnetic moment ( 3.25 µ B ) is observed at +2% tensile strain. Figure 3a and c shows the spin-polarized band structure under the application of −2% compressive and +2% tensile strain, respectively. Moreover, one can infer that in the case of −2%, spin-up impurity states are localized near the Fermi level occupying Co-d states, whereas, for +2%, the spin-down Co-d states are filled to produce the magnetism. This behavior can be correlated to the opposite nature of the magnetic ordering in compressive and tensile strains and may show good performance in transferring one particular spin oriented electrons such as in a spin filtering devices. The highly dispersive impurity d-band lines around the fermi level in case of +2% tensile strain renders more electrons on the edge of Co atom spin-polarised when compared to the less dispersive d-band lines in case of −2% compressive strain. Interestingly, the bandgap is minimum at strain ±2% , which ensures the bandgap tunability plays a significant role in TMDC based DMSs. The variation of half-metallic gaps with two different U values (U = 2.5, 3) under unstrained and strained (both compressive and tensile) conditions have also been incorporated which are listed in supplementary information (see SI: Table S3). However, the bandgap is further increased at higher applied strains (see SI: Fig. S2) resulting in the decrease of the magnetic moments in both strain direction.

W-rich S-rich
To further elucidate the electronic and magnetic properties in strained (Co, W)S 2 ML, the total density of states (TDOS) and partial density of states (PDOS) are plotted in Fig. 4. The pristine WS 2 behaves as a nonmagnetic semiconductor, which can be inferred from the symmetric nature of spin-up and spin-down states (see SI: Fig. S3), which is well consistent with previously studied results 42 . Unlike pristine WS 2 , in (Co, W)S 2 ML, the splitting of spin states is observed, giving rise to net magnetic moment ( 2.58 µ B ), as depicted in Fig. 4b. The origin of magnetization in (Co, W)S 2 ML is mainly contributed from the additional three unpaired electrons of Co − 3d 7 4S 2 than W − 5d 4 4s 2 with 1.67 µ B per Co-atom giving rise to n-type doping. However, the contribution from the nearest neighbor W ( 0.03 µ B per W atom) and S ( 0.11 µ B per S atom) atoms to the net magnetic moment is less than the contribution arising from Co. Similarly, the contribution from the Co-atom and the nearest neighbour W, S atoms at various studied compressive and tensile strains are listed (see SI: Table S1). The TDOS and PDOS for extended 5 × 5 supercell (see SI: Fig. S3) show no significant modification in magnetization after Co-doping. To understand the magnetic exchange behavior in (Co, W)S 2 ML, the interaction between foreign Co-d with neighbouring W-d and S-p are need to be considered. The values of magnetic moments depend on the hybridization among the Co-d, W-d and S-p. The occupied states near the Fermi level mainly arise from Co-3d and W-5d orbitals in the majority spin channel, as evident from Fig. 4b. As the Fermi level lies within the partially occupied majority band of the impurity states, expecting a double exchange coupling between Co and W 43,44 . Similarly, from Fig. 4a, it can be inferred that the broadening of the majority bands along with the Fermi level passing through them indicates a stronger double exchange mechanism under compressive strain 40 . However, in tensile strain, the broadening of the minority charge carriers can be noticed across the Fermi level forming highly delocalized minority spin channels which allows maximum magnetic moment under strained (Co, W)S 2 ML. The interaction of Co-d with nearest adjacent S-p can be explained in terms of p-d hybridization mechanism. The Co atoms are strongly coupled with their neighboring S atoms around the doping sites due to the hybridization of out-of-plane Co-3d and S-3p orbitals creating considerable unbalanced spin populations in spin-split impurity bands near the fermi level. It is worth noting from Table 2 that the Co-S bond length is lower than the W-S bond under compressive strain, which is found to be reversed in the case of applied tensile strain. This change in bond lengths leads to different hybridization mechanisms for varied strains. From the above analysis, it can be believed that the competition between the Co-S and W-S bond lengths may be considered as the prime factor in modulating the magnetic properties under strain engineering.
Next, we consider the spin density distribution to support our understanding of the exchange coupling from the PDOS near the Fermi level, shown in Fig. 5. Figure 5b elucidates the spin polarization between the Co atom and its nearest W/S in unstrained (Co, W)S 2 ML. The coupling between the foreign Co-atom and host W-atom results in parallel spin alignment indicating a double exchange interaction as evident from PDOS in Fig. 4b. However, the interaction between Co and nearest neighbor S results in p-d hybridization from the out-of-plane orbitals as depicted from Fig. 5b. The magnetic coupling behavior in a TM doped TMDCs systems can also be tuned by the introduction of cation and anion vacancies 45    www.nature.com/scientificreports/ delocalization in Co-doped WS 2 monolayer under both unstrained and strained condition has weak effect to the ferromagnetism. From the spin density distribution of (Co, W)S 2 ML, it can be perceived that Co-doping induces a long-range magnetic interaction with nearest neighboring W/S atoms. Moreover, the spin distribution is more localized around its magnetic centres for strained systems when compared the spin distribution around the magnetic centres of unstrained system. In case of unstrained (Co, W)S 2 ML, a continuous network of the magnetically coupled TM impurities over local clusters establishes a long-range magnetic interaction leading to lower percolation threshold. However, under the application of compressive and tensile strains, the local spin clusters around the impurity atoms establishes a short-range magnetic interaction, which limits to the nearest neighbour S atoms. In such a scenario, a higher doping concentration is required to reach the percolation threshold 47 . Additionally, it can be concluded that the induced spin density at the dopant site is maximum under +2% tensile strain, which reflects the ultimate magnetic moment in this case. For the extended 5 × 5 × 1 Co-doped WS 2 supercell, the magnetic coupling between dopant and host atoms shows similar behavior (see SI: Fig. S4). However, with the increase in either compressive or tensile strain, the magnetic moment decreases due to lower spin polarization. The origin of FM behavior in WS 2 ML after doping and strain engineering can be further explained based on orbital decomposed PDOS analysis of the Co-atom, as shown in Fig. 6. According to ligand field theory, the 3d states of Co atom can be split into single [ a 1 ( d Z 2 )] and two two-fold degenerate [ e 1 (d xy,x 2 −y 2 ) , e 2 (d xz,yz ] states. Intra-atomic Hund's exchange splitting ( ex ) is determined by the energy difference of e 1 orbital between the spin-up and spin-down states, whereas the energy difference between e 1 and e 2 orbitals is referred to as crystal field splitting ( cf ) 48,49 . The spin-splitting in Co-doped and strain engineered WS 2 near the Fermi level mainly results from the exchange splitting and crystal field splitting. As suggested by Pan et. al., the FM behavior in TMDCs due to TM doping arises from the competition between the exchange splitting and crystal field splitting 50 . Figure 6d represents the schematic for exchange and crystal field splitting under − 2% compressive and +2% tensile strain compared with the unstrained condition. Moreover, the exchange splitting dominates over crystal field splitting in our studied system due to n-type Co-dopant. From Fig. 6d, it can be observed that the difference between cf and ex increases under the application of strains, which in turn reflects the increased magnetic moments, as listed from Table 2. This ex and cf at each studied strain is plotted in Fig. 6e. The detailed analysis confirms that the larger the separation between cf and ex , the greater will be magnetic moments, as evident from Fig. 6e 50 .
Our detailed observation from the DFT calculation reveals that the Co-doped WS 2 ML behaves as FM in nature. In order to better implement an application point of view, the magnetization reversal of WS 2 under Co-doping engineering needs to be understood. After substituting Co at W-site, there are two possibilities of formation of anisotropy, one is uniaxial, and another one is biaxial anisotropy. As suggested from the previous reports, the uniaxial anisotropy value is high due to the formation of a larger coercive field ( H c ) with the easy axis measurement 51 . In contrast, the biaxial anisotropy strength is lower, resulting in low H c values. Experimental results show that the Co-doped WS 2 bulk, nanosheets, ML exhibit hysteresis at very low H c (few hundreds of Oe) at room temperature 35,36,52 . A similar FM signature is also reported in the Co-doped WSe 2 system 36 . However, the nature of the hysteresis is found to be slanted in all these cases. This intrigues us to understand the behavior of magnetization in our system. However, the FM behavior in Co-doped in WS 2 has not been thoroughly understood in terms of magnetization reversals. The quest for slanted hysteresis in all the reported cases is still elusive in the scientific community. Here we address the behavior of magnetization reversal and the magnetic effects of the system under strain. We approach micromagnetic modeling by using an open freeware object-orientedmicromagnetic-framework (OOMMF) package to perform qualitative analysis for ferromagnetism 34 . This is a to understand the intrinsic magnetic properties of the FM systems, which has very limited reports available on TMDCs materials. For practical applications, we need a sample dimension of nm-in-range. For micromagnetic modeling, we used a sample dimension of 50 × 50 × 1 nm 3 , and cell size is 1 × 1 × 0.5 nm 3 to compute the simulation. This micromagnetic simulation governs by Landau-Lifshitz-Gilbert (LLG) equation, can be written as 53 : dM dt provides the information of the M over time, first term represents precision of moments, while the second term is responsible for damping. γ denotes the gyromagnetic ratio, α stands for damping factor, these values are kept constant throughout the process of simulation. M s represents saturation magnetization, and H eff is the effective field of demagnetization and external magnetic field. M s is correlated with the total energy and H eff of the system, H eff = −(E total /M)/(µ 0 M s ), where E total is the total energy density of the system. E total would be the sum of all micromagnetics energies, which can be written as: E exch is the exchange energy, E anis is the magnetocrystalline anisotropy energy, E demag is the demagnetizing or stray field energy, E Zeeman is due to an external field, and E me is the magnetoelastic energy.
In the micromagnetic simulation, magnetic parameters are set the same throughout the cells. A stable state is achieved from a random energetic state. The simulation is performed implementing finite difference methods. The parameters are optimized by looking at various Co-doped systems and the magnetic properties associated with atomic layers of Co films 51 . The optimized parameters like exchange length: 2.1 × 10 −6 erg/cm , saturation magnetization: 30 emu.cm −3 , anisotropy values: (4 − 6) × 10 −3 erg/cm 3 are considered for micromagnetics simulation. When Co is doped initially at W-site, the magnetic moments are randomly dispersed in the cell matrix. On applying a sufficient magnetic field, the Co moments get oriented in the direction of external field. The snapshots of the various simulated states for the orientation of the spin moments while achieving from a    S5). A single magnetic domain state is observed throughout the hysteresis. Figure 7a represents two kinds of hysteresis taken at two different easy axes. The hysteresis for (1, 0, 0) is quite square-in-nature, where the nucleation of domains occurred near to the remanence. However, this signature of the loop contradicts the loop obtained for Co-doped in WSe 2 , and WS 2 ML systems 35,36 . To achieve the experimental signature of the loop, which is slightly slanted and nucleation occurs at a distance from the remanence, we consider magnetization orientation tilted by 20% off from the original orientation, i.e, (1, 0.2, 0) direction. As a result, the obtained hysteresis is slightly slanted from the previous orientation (1, 0, 0), which is nearly similar to the hysteresis obtained in the WS 2 system 35 . From the previously studied experimental results, it has been observed that WSe 2 shows relatively strong ferromagnetism after doping with TM Co 36 . The experimental observation shows tilting magnetization hysteresis curve at all studied temperatures ranging from very low 5 K to 300 K 36 . A similar FM response was reported in a wide range of temperature including room temperature of WS 2 35 . Although the magnetic parameters of WS 2 system are tuned for various factors such as size dependent (power and nanoribbons), thermal effect, doping, defect engineering, etc.; but, the tilting behavior is persist with all conditions 25,35,36,52,54 . Hence, the behavior of tilting M-H behavior has unaffected under these conditions, which ensures the anisotropic energy of the system attributes the tilting behavior in magnetization hysteresis curve of TMDCs. The coercive field for (1, 0, 0) orientation is around 345 Oe, whereas the same has been reduced to 190 Oe for (1, 0.2, 0). This coercive field is quite agreeing with the result reported in the WS 2 system 35 . Therefore, after Co doping into the system, it is worthy of mentioning that the effective magnetization orientation is slightly off from the easy anisotropic axis.
From DFT calculations it is evident that the FM behavior can be tuned under the application of strain. Here, the micromagnetic simulation provides the effect of strain on the larger scale in Co-doped WS 2 ML. In this www.nature.com/scientificreports/ case, if strain is applied to the system, we expect an alteration in the anisotropic values. In the literature also, it is reported that strain can certainly control the anisotropy in thin films 55,56 . Here, we vary anisotropy values to understand the effect of strain in the Co-doped WS 2 system. In this case, we tune the anisotropy values from −4 × 10 −3 erg/cm 3 to −7.5 × 10 −3 erg/cm 3 to observe the changes in the magnetic properties. The coercive field is gradually decreased with the increase in the anisotropy value, which is represented in Fig. 7b. On the other way, FM properties are getting affected by enhancing anisotropy values. This behavior is quite analogous with the results obtained from DFT calculations, where the increase in the percentages in strain values leads to a decrease in FM nature. Our results reinvigorate the FM coupling behavior in Co-doped WS 2 ML under various strain conditions; a further understanding of magnetic reversal in this system may pave the way for next-generation spintronics and straintronics applications.

Conclusion
We explore the strain-induced ferromagnetism in transition metal Co-doped WS 2 monolayer by using firstprinciples DFT calculations and micromagnetic simulation. Co-doping marks a significant change in magnetic properties with an impressive magnetic moment of 2.58 µ B . The magnetic exchange interaction is found to be double exchange coupling between Co and W and strong p-d hybridization between Co and nearest S, which is further verified from spin density distribution. We find that the resultant impurity bands of the Co-doped WS 2 plays a role of seed to drive novel electronic and magnetic properties under applied strain. Among several biaxial strains, the magnetic moment is found to be a maximum of 3.25 µ B at 2% tensile and 2.69 µ B at −2% compressive strain due to strong double exchange coupling and p-d hybridization among foreign Co and W/S. Further magnetic moments at higher applied strain is decreased due to reduced spin polarization. In addition, the competition between exchange splitting and crystal field splitting of Co d-orbital plays a significant role to determine these values of magnetic moments under the application of strain. From the micromagnetic simulation, it is confirmed that the Co-doped WS 2 monolayer shows slanted ferromagnetic hysteresis with a low coercive field. The effect of higher strain suppresses the ferromagnetic nature, which has a good agreement with the results obtained from DFT calculations. Our findings indicate that induced magnetism in WS 2 monolayer under Co-doping promotes the application of 2D TMDCs for the nano-scale spintronics, and especially, the strain-mediated magnetism can be a promising candidate for future straintronics applications.