<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">WJCMP</journal-id><journal-title-group><journal-title>World Journal of Condensed Matter Physics</journal-title></journal-title-group><issn pub-type="epub">2160-6919</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/wjcmp.2015.51004</article-id><article-id pub-id-type="publisher-id">WJCMP-54131</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Coexistence of Superconductivity and Ferromagnetism in Superconducting HoMo&lt;sub&gt;6&lt;/sub&gt;S&lt;sub&gt;8&lt;/sub&gt;
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>adesse</surname><given-names>Desta</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Gebregziabher</surname><given-names>Kahsay</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Physics, Addis Ababa University, Addis Ababa, Ethiopia</addr-line></aff><aff id="aff2"><addr-line>Department of Physics, College of Science, Bahir Dar University, Bahir Dar, Ethiopia</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>tad4jju@gmail.com(AD)</email>;<email>michige_90@yahoo.com(GK)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>21</day><month>01</month><year>2015</year></pub-date><volume>05</volume><issue>01</issue><fpage>27</fpage><lpage>36</lpage><history><date date-type="received"><day>29</day>	<month>January</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>13</month>	<year>February</year>	</date><date date-type="accepted"><day>16</day>	<month>February</month>	<year>2015</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  This work focuses on the theoretical investigation of the coexistence of superconductivity and ferromagnetism in the superconducting HoMo
  <sub>6</sub>S
  <sub>8</sub>. By developing a model Hamiltonian for the system and using the Green’s function formalism and equation of motion method, we have obtained expressions for superconducting transition temperature (
  T<sub>c</sub>), magnetic order temperature (
  T<sub>m</sub>), superconductivity order parameter (
  D
  ) and magnetic order parameter (η). By employing the experimental and theoretical values of the parameters in the obtained expressions, phase diagrams of energy gap parameter versus transition temperature, superconducting transition temperature versus magnetic order parameter and magnetic order temperature versus magnetic order parameter are plotted separately. By combining the phase diagrams of superconducting transition temperature versus magnetic order parameter and magnetic order temperature versus magnetic order parameter, we have demonstrated the possible coexistence of superconductivity and ferromagnetism in superconducting HoMo<sub>6</sub>S<sub>8</sub>.&lt;
 
</p></abstract><kwd-group><kwd>Superconductivity</kwd><kwd> Ferromagnetism</kwd><kwd> Coexistence</kwd><kwd> Green’s Function</kwd><kwd> HoMo&lt;sub&gt;6&lt;/sub&gt;S&lt;sub&gt;8&lt;/sub&gt;</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Superconductivity was discovered in 1911 by Kamerlingh Onnes [<xref ref-type="bibr" rid="scirp.54131-ref1">1</xref>] when a so-called “Blue Boy” noticed that the resistivity of Hg metal vanished abruptly at a temperature of about 4.2 K. Ferromagnetism is a phenomenon by which a material can exhibit a spontaneous magnetization and is one of the strongest forms of magnetism. It is responsible for most of the magnetic behaviors encountered in everyday life and is the basis for all permanent magnets (as well as the metals that are noticeably attracted to them). In particular, a material is ferromagnetic in narrower sense only if all of its magneticions add a positive contribution to the net magnetization. If some of the magneticions subtract from the net magnetization, that is, if they are partially anti-aligned, then the material is ferrimagnetic. If the ions anti-align completely so as to have zero netmagnetization, despite the magnetic ordering, then it is an antiferromagnet. Thus, ferromagnetic materials exhibit parallel alignment of moments resulting in large net magnetization even in the absence of a magnetic field.</p><p>Superconductivity in Ferromagnetic must result from a different type of electronpairing mechanisms. In these materials, electrons with spins pointing in the same direction team up with each other to form Cooper pairs with one unit of spin resulting in a triplet superconductivity. In contrast, conventional superconductivity also known as s-wave singlet superconductivity occurs when electrons with oppositespins bind together to form Cooper pairs with zero momentum and spin.</p><p>The coexistence of superconductivity and ferromagnetism has been studied theoretically and experimentally. The coexistence of ferromagnetism and superconductivity was first addressed theoretically by Ginzburg in 1957 [<xref ref-type="bibr" rid="scirp.54131-ref2">2</xref>] and experimental investigation was made by Matthias et al. [<xref ref-type="bibr" rid="scirp.54131-ref3">3</xref>] . The interplay between superconducting and ferromagnetic long range order has been recently attracting new interest due to the discovery of superconductivity in ferromagnetic compounds such as UGe<sub>2</sub> [<xref ref-type="bibr" rid="scirp.54131-ref4">4</xref>] , URhGe [<xref ref-type="bibr" rid="scirp.54131-ref5">5</xref>] , ZrZn<sub>2</sub> [<xref ref-type="bibr" rid="scirp.54131-ref6">6</xref>] , and in RuSr<sub>2</sub>RECu<sub>2</sub>O<sub>8</sub> compounds (with RE = Eu or Gd) [<xref ref-type="bibr" rid="scirp.54131-ref7">7</xref>] . The relationship between magnetism and superconductivity has received renewed attention since the discovery of ternary superconducting materials which also achieved long-range magnetic ordering at low temperatures. Ferromagnetic alignment can be expected to be strongly opposed by superconductivity. Such a long-period magnetic ordering was actually found in HoMo<sub>6</sub>S<sub>8</sub> and in ErRh<sub>4</sub>B<sub>4</sub>. In ErRh<sub>4</sub>B<sub>4</sub>, Sinha et al. [<xref ref-type="bibr" rid="scirp.54131-ref8">8</xref>] carried out a detailed study on a single crystal in order to characterize this phase. For HoMo<sub>6</sub>S<sub>8</sub>, the study was done by Lynn et al. [<xref ref-type="bibr" rid="scirp.54131-ref9">9</xref>] only on polycrystalline samples.</p><p>In HoMo<sub>6</sub>S<sub>8</sub>, the ferromagnetic state destroys the superconductivity at sufficiently low temperatures. Recently, an experiment on HoMo<sub>6</sub>S<sub>8</sub> [<xref ref-type="bibr" rid="scirp.54131-ref10">10</xref>] has shown that, the superconducting ordering parameter has a distinct maximum between some critical temperatures T<sub>c</sub><sub>1</sub> and T<sub>c</sub><sub>2</sub> (lower and upper superconducting critical temperatures), respectively and vanishes for T &lt; T<sub>c</sub><sub>1</sub> and T &gt; T<sub>c</sub><sub>2</sub>.</p><p>Among the “Chevrel phases”, HoMo<sub>6</sub>S<sub>8</sub> has been extensively studied in recent years [<xref ref-type="bibr" rid="scirp.54131-ref11">11</xref>] . HoMo<sub>6</sub>S<sub>8</sub> becomes superconducting at T<sub>c</sub><sub>1</sub> ≌ 1.82 K, but at a lower temperature T<sub>c</sub><sub>2</sub> ≌ 0.64 K, it re-enters the normal state at the onset of long range ferromagnetic order. In a narrow temperature range T<sub>c</sub><sub>2</sub> &lt; T &lt; T<sub>m</sub>, superconductivity coexists with a modulated magnetic structure [<xref ref-type="bibr" rid="scirp.54131-ref12">12</xref>] .</p></sec><sec id="s2"><title>2. Model System Hamiltonian</title><p>In order to study the coexistence of ferromagnetism and superconductivity in superconducting HoMo<sub>6</sub>S<sub>8</sub> theoretically in general and to find the expressions for transition temperature and order parameters in particular, a system of conduction and localized electrons have been considered. The exchange interaction acts between the conduction and the localized electrons. Thus, within the frame work of the BCS model [<xref ref-type="bibr" rid="scirp.54131-ref13">13</xref>] , the model system Hamiltonian can be formulated as follows.</p><disp-formula id="scirp.54131-formula94"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x5.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.54131-formula95"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x6.png"  xlink:type="simple"/></disp-formula><p>and is the Hamiltonian or energy of mobile (conduction) electrons and localized electrons respectively.</p><p>Here, the operators <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x7.png" xlink:type="simple"/></inline-formula> <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x8.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x9.png" xlink:type="simple"/></inline-formula> <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x10.png" xlink:type="simple"/></inline-formula> are the creation (annihilation) operators for conduction</p><p>and localized electrons respectively with the wave vector <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x11.png" xlink:type="simple"/></inline-formula> and the spin projection on z-axis σ = ↑ or ↓. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x12.png" xlink:type="simple"/></inline-formula>is the one electron kinetic energy measured relative to the chemical potential. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x13.png" xlink:type="simple"/></inline-formula>is the interaction (electron- electron) through boson (phonon) exchange and is given by,</p><disp-formula id="scirp.54131-formula96"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x14.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x15.png" xlink:type="simple"/></inline-formula> defines the matrix element of the interaction potential. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x16.png" xlink:type="simple"/></inline-formula>is the interaction term between con-</p><p>duction electrons and localized electrons due to some unspecified mechanism with some coupling constant (α) and is expressed as,</p><disp-formula id="scirp.54131-formula97"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x17.png"  xlink:type="simple"/></disp-formula></sec><sec id="s3"><title>3. Equation of Motion for Mobile (Conduction) Electrons</title><p>Now, let us evaluate the following commutation relation,</p><disp-formula id="scirp.54131-formula98"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x18.png"  xlink:type="simple"/></disp-formula><p>From which we obtain,</p><disp-formula id="scirp.54131-formula99"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x19.png"  xlink:type="simple"/></disp-formula><p>Following similar procedure as above, we get,</p><disp-formula id="scirp.54131-formula100"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x20.png"  xlink:type="simple"/></disp-formula><p>and</p><disp-formula id="scirp.54131-formula101"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x21.png"  xlink:type="simple"/></disp-formula><p>Substituting Equations (6), (7) and (8) into the equation of motion given by,</p><disp-formula id="scirp.54131-formula102"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x22.png"  xlink:type="simple"/></disp-formula><p>we obtain,</p><disp-formula id="scirp.54131-formula103"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x23.png"  xlink:type="simple"/></disp-formula><p>In general, we have to write the higher order Green’s function into lower order Green’s function by using Wick’s theorem. Thus, we have,</p><disp-formula id="scirp.54131-formula104"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x24.png"  xlink:type="simple"/></disp-formula><p>Now, substituting Equation (11) into Equation (10), we get,</p><disp-formula id="scirp.54131-formula105"><graphic  xlink:href="http://html.scirp.org/file/4-4800277x25.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.54131-formula106"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x26.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x27.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x28.png" xlink:type="simple"/></inline-formula></p><p>One can also obtain the equation of motion for the expression <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x29.png" xlink:type="simple"/></inline-formula> and obtain,</p><disp-formula id="scirp.54131-formula107"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x30.png"  xlink:type="simple"/></disp-formula><p>For <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x31.png" xlink:type="simple"/></inline-formula> <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x32.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x33.png" xlink:type="simple"/></inline-formula>, we obtain,</p><disp-formula id="scirp.54131-formula108"><label>(14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x34.png"  xlink:type="simple"/></disp-formula><p>Now, using Equations (12) and (14), the equation of motion becomes,</p><disp-formula id="scirp.54131-formula109"><graphic  xlink:href="http://html.scirp.org/file/4-4800277x35.png"  xlink:type="simple"/></disp-formula><p>From which we obtain,</p><disp-formula id="scirp.54131-formula110"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x36.png"  xlink:type="simple"/></disp-formula><p>Using the relation for ∆, given by,</p><disp-formula id="scirp.54131-formula111"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x37.png"  xlink:type="simple"/></disp-formula><p>and by changing the summation into integration and by introducing the density of states at the Fermi level, (N(0)), we get,</p><disp-formula id="scirp.54131-formula112"><label>(17)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x38.png"  xlink:type="simple"/></disp-formula><p>Now, changing<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x39.png" xlink:type="simple"/></inline-formula>, we use the Matsubara frequency [<xref ref-type="bibr" rid="scirp.54131-ref14">14</xref>] ,</p><disp-formula id="scirp.54131-formula113"><label>(18)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x40.png"  xlink:type="simple"/></disp-formula><p>Now, using Equation (18) in Equation (17), we get,</p><disp-formula id="scirp.54131-formula114"><label>(19)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x41.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x42.png" xlink:type="simple"/></inline-formula> Since attraction is effective in the region <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x43.png" xlink:type="simple"/></inline-formula> and taking the density of state to be constant in this region and using the relation,</p><disp-formula id="scirp.54131-formula115"><graphic  xlink:href="http://html.scirp.org/file/4-4800277x44.png"  xlink:type="simple"/></disp-formula><p>We can write Equation (19) as,</p><disp-formula id="scirp.54131-formula116"><label>(20)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x45.png"  xlink:type="simple"/></disp-formula><p>Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x46.png" xlink:type="simple"/></inline-formula></p><disp-formula id="scirp.54131-formula117"><label>(21)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x47.png"  xlink:type="simple"/></disp-formula><sec id="s3_1"><title>3.1. Effect of Temperature on Superconducting Order Parameter (∆) and Magnetic Order Parameter (η)</title><p>Now, let us study equation (21) by considering different cases.</p><p>Case (I): As <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x48.png" xlink:type="simple"/></inline-formula> so that, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x49.png" xlink:type="simple"/></inline-formula></p><p>Hence, Equation (21) becomes,</p><disp-formula id="scirp.54131-formula118"><label>(22)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x50.png"  xlink:type="simple"/></disp-formula><p>Using the integral <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x51.png" xlink:type="simple"/></inline-formula> where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x52.png" xlink:type="simple"/></inline-formula>, Equation (21) becomes,</p><disp-formula id="scirp.54131-formula119"><label>(23)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x53.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.54131-formula120"><label>(24)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x54.png"  xlink:type="simple"/></disp-formula><p>For<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x55.png" xlink:type="simple"/></inline-formula>, Equation (24) reduces to,</p><disp-formula id="scirp.54131-formula121"><label>(25)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x56.png"  xlink:type="simple"/></disp-formula><p>This implies that,</p><disp-formula id="scirp.54131-formula122"><label>(26)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x57.png"  xlink:type="simple"/></disp-formula><p>For <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x58.png" xlink:type="simple"/></inline-formula> Equation (26) reduces to the well-known BCS model.</p><p>The experimental value of HoMo<sub>6</sub>S<sub>8</sub> is, T<sub>c</sub> ≈ 1.82 K.</p><p>Thus, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x59.png" xlink:type="simple"/></inline-formula>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x59.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x60.png" xlink:type="simple"/></inline-formula> (for BCS model).</p><p>Case (II): At<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x61.png" xlink:type="simple"/></inline-formula>,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x62.png" xlink:type="simple"/></inline-formula>. Thus, we get,<sub> </sub></p><disp-formula id="scirp.54131-formula123"><label>(27)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x63.png"  xlink:type="simple"/></disp-formula><p>Now, employing Equation (27) and the experimental value of T<sub>c</sub> for the superconducting HoMo<sub>6</sub>S<sub>8</sub> and plausible approximations for other parameters, we plotted the transition temperature (T<sub>c</sub>) versus magnetic ordering parameter (η) as shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p><p>For η = 0, we get the expression for T<sub>c</sub> to be,</p><disp-formula id="scirp.54131-formula124"><label>(28)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x64.png"  xlink:type="simple"/></disp-formula><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Transition temperature (T<sub>c</sub>) versus magnetic order parameter (η) for superconducting HoMo<sub>6</sub>S<sub>8</sub></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-4800277x65.png"/></fig></sec><sec id="s3_2"><title>3.2. Equation of Motion for Localized Electrons</title><p>Using Green’s function formalism, the equation of motion for the localized electrons is obtained to be,</p><disp-formula id="scirp.54131-formula125"><label>(29)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x66.png"  xlink:type="simple"/></disp-formula><p>Now, using the Hamiltonian given in Equation (1), we evaluated the commutation <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x67.png" xlink:type="simple"/></inline-formula> and obtained,</p><disp-formula id="scirp.54131-formula126"><label>(30)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x68.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.54131-formula127"><label>(31)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x69.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x70.png" xlink:type="simple"/></inline-formula></p><p>Applying similar procedure as above and assuming<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x71.png" xlink:type="simple"/></inline-formula>, we can obtain the expression for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x72.png" xlink:type="simple"/></inline-formula> to be,</p><disp-formula id="scirp.54131-formula128"><label>(32)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x73.png"  xlink:type="simple"/></disp-formula><p>Now, from Equations (31) and (32), we get,</p><disp-formula id="scirp.54131-formula129"><label>(33)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x74.png"  xlink:type="simple"/></disp-formula><p>From which we get,</p><disp-formula id="scirp.54131-formula130"><label>(34)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x75.png"  xlink:type="simple"/></disp-formula></sec><sec id="s3_3"><title>3.3. Correlation between Conduction and Mobile Electrons</title><p>The equation of motion that shows the correlation between the conduction and localized electrons can be demonstrated. Using similar definition as for ∆, we can write the magnetic ordering parameter, η as,</p><disp-formula id="scirp.54131-formula131"><label>(35)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x76.png"  xlink:type="simple"/></disp-formula><p>Changing the summation into integration and by introducing the density of states, N(0), we get,</p><disp-formula id="scirp.54131-formula132"><label>(36)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x77.png"  xlink:type="simple"/></disp-formula><p>Using the Matsubara frequency,</p><disp-formula id="scirp.54131-formula133"><graphic  xlink:href="http://html.scirp.org/file/4-4800277x78.png"  xlink:type="simple"/></disp-formula><p>Equation (36) becomes,</p><disp-formula id="scirp.54131-formula134"><label>(37)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x79.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x80.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x80.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x81.png" xlink:type="simple"/></inline-formula></p><p>Now, let us first solve the following expression.</p><disp-formula id="scirp.54131-formula135"><label>(38)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x82.png"  xlink:type="simple"/></disp-formula><p>Using Laplace’s transform and Matsubara frequency, Equation (38) becomes,</p><disp-formula id="scirp.54131-formula136"><graphic  xlink:href="http://html.scirp.org/file/4-4800277x83.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.54131-formula137"><label>(39)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x84.png"  xlink:type="simple"/></disp-formula><p>and</p><disp-formula id="scirp.54131-formula138"><label>(40)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x85.png"  xlink:type="simple"/></disp-formula><p>Then,</p><disp-formula id="scirp.54131-formula139"><label>(41)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x86.png"  xlink:type="simple"/></disp-formula><p>Since ∆<sub>l</sub> is very small, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x87.png" xlink:type="simple"/></inline-formula>can be neglected and thus Equation (41) becomes,</p><disp-formula id="scirp.54131-formula140"><graphic  xlink:href="http://html.scirp.org/file/4-4800277x88.png"  xlink:type="simple"/></disp-formula><p>From which we get,</p><disp-formula id="scirp.54131-formula141"><label>(42)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x89.png"  xlink:type="simple"/></disp-formula><p>Using Equation (42) and the experimental value, T<sub>m</sub> ≈ 0.67 K for HoMo<sub>6</sub>S<sub>8</sub> and some plausible approximations for other parameters in the equation, we plot the magnetic order temperature (T<sub>m</sub>) versus magnetic order parameter as shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>.</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Magnetic order temperature (T<sub>m</sub>) versus magnetic order parameter (η) for superconducting HoMo<sub>6</sub>S<sub>8</sub></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-4800277x90.png"/></fig></sec><sec id="s3_4"><title>3.4. Equation of Motion for Pure Superconducting System</title><p>For pure superconducting system, that is, when magnetic order cannot appear or magnetic effect is zero, we can ignore η and our previous calculation gives the following results which is similar to the well-known BCS model.</p><p>As T → 0, η → 0 and tanh(βE/2) → 1, Equation (21) reduces to,</p><disp-formula id="scirp.54131-formula142"><label>(43)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x91.png"  xlink:type="simple"/></disp-formula><p>From which we obtain,</p><disp-formula id="scirp.54131-formula143"><label>(44)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x92.png"  xlink:type="simple"/></disp-formula><p>Furthermore, for T → T<sub>c</sub>, η = 0 and for low temperature, i.e.<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x93.png" xlink:type="simple"/></inline-formula>, Equation (21) yields,</p><disp-formula id="scirp.54131-formula144"><graphic  xlink:href="http://html.scirp.org/file/4-4800277x94.png"  xlink:type="simple"/></disp-formula><p>From which we get,</p><disp-formula id="scirp.54131-formula145"><label>(45)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x95.png"  xlink:type="simple"/></disp-formula><p>To obtain the temperature dependency of energy gap in Equation (21), we used the same techniques to solve the integral,</p><disp-formula id="scirp.54131-formula146"><label>(46)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x96.png"  xlink:type="simple"/></disp-formula><p>But from the BCS model<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x97.png" xlink:type="simple"/></inline-formula>, (as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x97.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x98.png" xlink:type="simple"/></inline-formula>).</p><p>For<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x99.png" xlink:type="simple"/></inline-formula>, Equation (46) can be simplified and obtain,</p><disp-formula id="scirp.54131-formula147"><graphic  xlink:href="http://html.scirp.org/file/4-4800277x100.png"  xlink:type="simple"/></disp-formula><p>Using the relation<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-4800277x101.png" xlink:type="simple"/></inline-formula>, we get,</p><disp-formula id="scirp.54131-formula148"><graphic  xlink:href="http://html.scirp.org/file/4-4800277x102.png"  xlink:type="simple"/></disp-formula><p>From which we can get,</p><disp-formula id="scirp.54131-formula149"><label>(47)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/4-4800277x103.png"  xlink:type="simple"/></disp-formula><p>Equation (47) shows how the superconducting order parameter, ∆(T) varies with temperature when η = 0 and is similar to the BCS model.</p><p>Using the experimental value, T<sub>c</sub> ≈ 1.82 K for HoMo<sub>6</sub>S<sub>8</sub> and some plausible approximations, we plot ∆ versus T<sub>c</sub> as shown in <xref ref-type="fig" rid="fig3">Figure 3</xref>.</p><p>Now, by combining <xref ref-type="fig" rid="fig1">Figure 1</xref> and <xref ref-type="fig" rid="fig2">Figure 2</xref>, we demonstrated the possible coexistence of superconductivity and ferromagnetism in HoMo<sub>6</sub>S<sub>8</sub> as shown in <xref ref-type="fig" rid="fig4">Figure 4</xref>.</p></sec></sec><sec id="s4"><title>4. Results and Discussion</title><p>In this section, we describe the results which are obtained using the model Hamiltonian developed. We obtain the expressions for the superconducting ordering parameter (∆) and magnetic order parameter (η) with respect to</p><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Energy gap parameter (∆) versus transition temperature (T<sub>c</sub>) for superconducting HoMo<sub>6</sub>S<sub>8</sub></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-4800277x104.png"/></fig><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Coexistence of superconductivity and ferromagnetism in superconducting HoMo<sub>6</sub>S<sub>8</sub></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/4-4800277x105.png"/></fig><p>superconducting transition temperature (T<sub>c</sub>) and magnetic order temperature (T<sub>m</sub>) respectively. First, using Equation (27) and the experimental value of T<sub>c</sub> for the superconducting HoMo<sub>6</sub>S<sub>8</sub> and plausible approximations for other parameters, we plotted the transition temperature (T<sub>c</sub>) versus magnetic order parameter (η) as shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>. As can be seen from the figure, when the magnetic order parameter increases the superconducting transition temperature decreases. Second, by employing Equation (42) and the experimental value, T<sub>m</sub> ≈ 0.67 K for HoMo<sub>6</sub>S<sub>8</sub> and some suitable approximations for the other parameters in the equation, we plotted the magnetic order temperature versus magnetic order parameter as demonstrated in <xref ref-type="fig" rid="fig2">Figure 2</xref>. From the figure, it is vivid that, as the magnetic order parameter increases the magnetic order temperature also increases. Furthermore, the superconducting order parameter (∆) is expressed as a function of the transition temperature (T<sub>c</sub>) and is plotted in <xref ref-type="fig" rid="fig3">Figure 3</xref>. The expression we obtained for the pure superconductor in Equation (47) is in agreement with the BCS model for η = 0. It is clear that from <xref ref-type="fig" rid="fig3">Figure 3</xref>, the superconducting order parameter, which is the measure of pairing energy decreases with increasing temperature and vanishes at the transition temperature (T<sub>c</sub>). From <xref ref-type="fig" rid="fig4">Figure 4</xref>, we observe that, T<sub>c</sub> decreases with increasing η, whereas T<sub>m</sub> increases with increasing η and there is a small region of temperature where both superconductivity and ferromagnetism coexist in HoMo<sub>6</sub>S<sub>8</sub>. Our finding is in agreement with the experimental observation [<xref ref-type="bibr" rid="scirp.54131-ref15">15</xref>] .</p></sec><sec id="s5"><title>5. Conclusion</title><p>In the present work, we have demonstrated the basic concepts of superconductivity with special emphasis on the BCS model and Cooper pair focusing on the interaction between superconductivity and ferromagnetism which are closely connected to the particular crystal of superconducting HoMo<sub>6</sub>S<sub>8</sub>. Employing the double time temperature dependent retarded Green’s functions formalism, we developed the model Hamiltonian for the system and derived equations of motion for conduction electrons, localized electrons and for pure superconducting system and carried out various correlations by using suitable decoupling procedures. In developing the model Hamiltonian, we considered spin triplet pairing mechanism and obtained expressions for superconducting order parameter, magnetic order parameter, superconducting transition temperature and magnetic order temperature. By using appropriate experimental values and considering suitable approximations, we plotted figures using the equations developed. As is well-known, superconductivity and ferromagnetism are two cooperative phenomena which are mutually antagonistic since superconductivity is associated with the pairing of electron states related to time reversal while in the magnetic states the time reversal symmetry is lost. Because of this, there is a strong competition between the two phases. This competition between superconductivity and magnetism made coexistence unlikely to occur. However, the model we employed in this work, shows that, there is a small region of temperature where both superconductivity and ferromagnetism can coexist in superconducting HoMo<sub>6</sub>S<sub>8</sub>.</p></sec><sec id="s6"><title>Cite this paper</title><p>TadesseDesta,GebregziabherKahsay, (2015) Coexistence of Superconductivity and Ferromagnetism in Superconducting HoMo<sub>6</sub>S<sub>8</sub>. World Journal of Condensed Matter Physics,05,27-36. doi: 10.4236/wjcmp.2015.51004</p></sec></body><back><ref-list><title>References</title><ref id="scirp.54131-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Onnes, H.K. (1911) Akad. van Wetenschappen (Amsterdam), 14, 113.</mixed-citation></ref><ref id="scirp.54131-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Ginzburg, V.L. (1957) Ferromagnetic Superconductors. Soviet Physics—JETP, 4, 153.</mixed-citation></ref><ref id="scirp.54131-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Matthias, B.T., Suhl, H. and Corenzwit, E. (1958) Spin Exchange in Superconductors. Physical Review Letters, 1, 92.http://dx.doi.org/10.1103/PhysRevLett.1.92</mixed-citation></ref><ref id="scirp.54131-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Saxena, S.S., Agarwal, P., Ahilan, K., Grosche, F.M., Haselwimmer, R.K.W., Steiner, M.J., Pugh, E., Walker, I.R., Julian, S.R., Monthoux, P., Lonzarich, G.G., Huxley, A., Sheikin, I., Braithwaite, D. and Flouquet, J. (2000) Superconductivity on the Border of Itinerant-Electron Ferromagnetism in UGe2 Superconductivity. Nature, 406, 587.http://dx.doi.org/10.1038/35020500</mixed-citation></ref><ref id="scirp.54131-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Aoki, D., Huxley, A., Ressouche, E., Braithwaite, D., Flouquet, J., Brison, J.P., Lhote, E. and Paulsen, C. (2001) Coexistence of Superconductivity and Ferromagnetism in URhGe. Nature, 413, 613-616. http://dx.doi.org/10.1038/35098048</mixed-citation></ref><ref id="scirp.54131-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Pfleiderer, C., Uhlarz, M., Hayden, S.M., Vollmer, R., Lohneysen, H.V., Bernhoeft, N.R. and Lonzarich, G.G. (2001) Coexistence of Superconductivity and Ferromagnetism in the d-Band Metal ZrZn2. Nature, 412, 58.http://dx.doi.org/10.1038/35083531</mixed-citation></ref><ref id="scirp.54131-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Cuoco, M., Gentile, P. and Noce, C. (2003) Coexistence of Ferromagnetism and Singlet Superconductivity via Kinetic Exchange. Physical Review Letters, 91, Article ID: 197003. http://dx.doi.org/10.1103/PhysRevLett.91.197003</mixed-citation></ref><ref id="scirp.54131-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Sinha, S.K., Crabtree, G.W., Hinks, D.G., Mook, H.A. and Pringle, O.A. (1983) Neutron Scattering Studies of Magnetic Superconductors. Journal of Magnetism and Magnetic Materials, 31-34, 489.http://dx.doi.org/10.1016/0304-8853(83)90545-0</mixed-citation></ref><ref id="scirp.54131-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Lynn, J.W., Raggazoni, A., Pynn, R. and Joffrin, J. (1981) Observation of Long Range Magnetic Order in the Reentrant Superconductor HoMo6S8. Journal de Physique Lettres, 42, 45-49. http://dx.doi.org/10.1051/jphyslet:0198100420204500</mixed-citation></ref><ref id="scirp.54131-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Morales, F., Escudero, R., Briggs, A., Monceau, P., Horyn, R., Le Berre, F. and Pena, O. (1996) Point Contact Spectroscopy on the Ferromagnetic Superconductor HoMo6S8. Physica, 218B, 193. http://dx.doi.org/10.1016/0921-4526(95)00591-9</mixed-citation></ref><ref id="scirp.54131-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">(1982) Superconductivity in Ternary Compounds II. In: Maple and Fischer, Eds., Springer-Verlag, Berlin.</mixed-citation></ref><ref id="scirp.54131-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Lynn, J.W., Shirane, G., Thomlinson, W., Shelton, R.N. and Moncton, D.E. (1981) Magnetic Properties of the Reentrant Ferromagnetic Superconductor HoMo6S8. Physical Review B, 24, 3817. http://dx.doi.org/10.1103/PhysRevB.24.3817</mixed-citation></ref><ref id="scirp.54131-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Bardeen, J., Cooper, L.N. and Schrieffer, J.R. 1957) Theory of Superconductivity. Physical Review, 108, 1175. http://dx.doi.org/10.1103/PhysRev.108.1175</mixed-citation></ref><ref id="scirp.54131-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Zubarev, D.N. (1960) Double-Time Green Functions in Statistical Physics. Soviet Physics Uspekhi, 3, 320.http://dx.doi.org/10.1070/PU1960v003n03ABEH003275</mixed-citation></ref><ref id="scirp.54131-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Ishikawa, M., Fischer, O. and Muller, J. (1978) Long Range Magnetic Order in the Superconducting State of Heavy Rare Earth Molybdenum Sulfides and Their Pseudo Ternary Compounds. Journal de Physique, 39, C6-1379.</mixed-citation></ref></ref-list></back></article>