<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.1 20151215//EN" "http://jats.nlm.nih.gov/publishing/1.1/JATS-journalpublishing1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="1.1">
<front>
<journal-meta>
<journal-id journal-id-type="pmc">CMC</journal-id>
<journal-id journal-id-type="nlm-ta">CMC</journal-id>
<journal-id journal-id-type="publisher-id">CMC</journal-id>
<journal-title-group>
<journal-title>Computers, Materials &#x0026; Continua</journal-title>
</journal-title-group>
<issn pub-type="epub">1546-2226</issn>
<issn pub-type="ppub">1546-2218</issn>
<publisher>
<publisher-name>Tech Science Press</publisher-name>
<publisher-loc>USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">25741</article-id>
<article-id pub-id-type="doi">10.32604/cmc.2022.025741</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Article</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Hybrid Chaotic Salp Swarm with Crossover Algorithm for Underground Wireless Sensor Networks</article-title>
<alt-title alt-title-type="left-running-head">Hybrid Chaotic Salp Swarm with Crossover Algorithm for Underground Wireless Sensor Networks</alt-title>
<alt-title alt-title-type="right-running-head">Hybrid Chaotic Salp Swarm with Crossover Algorithm for Underground Wireless Sensor Networks</alt-title>
</title-group>
<contrib-group content-type="authors">
<contrib id="author-1" contrib-type="author" corresp="yes">
<name name-style="western"><surname>Ayedi</surname><given-names>Mariem</given-names></name><xref ref-type="aff" rid="aff-1">1</xref>
<xref ref-type="aff" rid="aff-2">2</xref><email>m.ayedi@psau.edu.sa</email>
</contrib>
<contrib id="author-2" contrib-type="author">
<name name-style="western"><surname>ElAshmawi</surname><given-names>Walaa H.</given-names></name><xref ref-type="aff" rid="aff-3">3</xref>
<xref ref-type="aff" rid="aff-4">4</xref></contrib>
<contrib id="author-3" contrib-type="author">
<name name-style="western"><surname>Eldesouky</surname><given-names>Esraa</given-names></name><xref ref-type="aff" rid="aff-1">1</xref>
<xref ref-type="aff" rid="aff-3">3</xref></contrib>
<aff id="aff-1"><label>1</label><institution>First Department of Computer Science, College of Computer Engineering and Sciences, Prince Sattam Bin Abdulaziz University</institution>, <addr-line>Al-Kharj, 11942</addr-line>, <country>Saudi Arabia</country></aff>
<aff id="aff-2"><label>2</label><institution>MEDIATRON Lab., SUP&#x0027;COM, Carthage University</institution>, <addr-line>Tunis, 2083</addr-line>, <country>Tunisia</country></aff>
<aff id="aff-3"><label>3</label><institution>Department of Computer Science, Faculty of Computers and Informatics, Suez Canal University</institution>, <addr-line>Ismailia, 41522</addr-line>, <country>Egypt</country></aff>
<aff id="aff-4"><label>4</label><institution>Faculty of Computer Science, Misr International University</institution>, <addr-line>Cairo</addr-line>, <country>Egypt</country></aff>
</contrib-group>
<author-notes>
<corresp id="cor1"><label>&#x002A;</label>Corresponding Author: Mariem Ayedi. Email: <email>m.ayedi@psau.edu.sa</email></corresp>
</author-notes>
<pub-date pub-type="epub" date-type="pub" iso-8601-date="2022-03-26"><day>26</day>
<month>03</month>
<year>2022</year></pub-date>
<volume>72</volume>
<issue>2</issue>
<fpage>2963</fpage>
<lpage>2980</lpage>
<history>
<date date-type="received"><day>03</day><month>12</month><year>2021</year></date>
<date date-type="accepted"><day>17</day><month>1</month><year>2022</year></date>
</history>
<permissions>
<copyright-statement>&#x00A9; 2022 Ayedi et al.</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Ayedi et al.</copyright-holder>
<license xlink:href="https://creativecommons.org/licenses/by/4.0/">
<license-p>This work is licensed under a <ext-link ext-link-type="uri" xlink:type="simple" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution 4.0 International License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p>
</license>
</permissions>
<self-uri content-type="pdf" xlink:href="TSP_CMC_25741.pdf"></self-uri>
<abstract>
<p>Resource management in Underground Wireless Sensor Networks (UWSNs) is one of the pillars to extend the network lifetime. An intriguing design goal for such networks is to achieve balanced energy and spectral resource utilization. This paper focuses on optimizing the resource efficiency in UWSNs where underground relay nodes amplify and forward sensed data, received from the buried source nodes through a lossy soil medium, to the aboveground base station. A new algorithm called the Hybrid Chaotic Salp Swarm and Crossover (HCSSC) algorithm is proposed to obtain the optimal source and relay transmission powers to maximize the network resource efficiency. The proposed algorithm improves the standard Salp Swarm Algorithm (SSA) by considering a chaotic map to initialize the population along with performing the crossover technique in the position updates of salps. Through experimental results, the HCSSC algorithm proves its outstanding superiority to the standard SSA for resource efficiency optimization. Hence, the network&#x0027;s lifetime is prolonged. Indeed, the proposed algorithm achieves an improvement performance of 23.6&#x0025; and 20.4&#x0025; for the resource efficiency and average remaining relay battery per transmission, respectively. Furthermore, simulation results demonstrate that the HCSSC algorithm proves its efficacy in the case of both equal and different node battery capacities.</p>
</abstract>
<kwd-group kwd-group-type="author">
<kwd>Underground wireless sensor networks</kwd>
<kwd>resource efficiency</kwd>
<kwd>chaotic theory</kwd>
<kwd>crossover algorithm</kwd>
<kwd>salp swarm algorithm</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1"><label>1</label><title>Introduction</title>
<p>With the precipitous growth of microelectronics and sensing technologies, Wireless Sensor Networks (WSNs) have been categorized as an intense research arena. The outstanding privileges of such networks, including their easy configuration, mobility, and flexibility, lead to their deployment in many environments. UWSNs are an important extension of WSNs applications for use in the underground world. These underground networks have a broad range of applications including tracking coal mining facilities, soil monitoring, and gas/oil pipeline control. In such networks, buried sensors continuously collect sensitive data regarding the sensed environment and forward it to the base station [<xref ref-type="bibr" rid="ref-1">1</xref>]. However, the limited communication range and the energy constraints along with the complex and unpredictable conditions of the underground medium are the primary challenges of UWSNs. Thus, the role of relay nodes is vital in UWSNs because they represent a promising method for achieving high bandwidth and expanding network coverage [<xref ref-type="bibr" rid="ref-2">2</xref>]. In UWSNs, communications among nodes come in three distinct channels: UnderGround-to-AboveGround (UG2AG), UnderGround-to-UnderGround (UG2UG), and AboveGround-to-UnderGround (AG2UG) [<xref ref-type="bibr" rid="ref-3">3</xref>]. The UG2AG channel connection is mainly used to transmit sensed data from buried sensors to relay nodes or aboveground base stations [<xref ref-type="bibr" rid="ref-3">3</xref>&#x2013;<xref ref-type="bibr" rid="ref-5">5</xref>]. Relay node deployment has been studied in UWSNs [<xref ref-type="bibr" rid="ref-5">5</xref>&#x2013;<xref ref-type="bibr" rid="ref-8">8</xref>]. In [<xref ref-type="bibr" rid="ref-5">5</xref>], an underground coal mine was divided into separate regions and addressed optimal relay node placement to support robust coverage of the network. In [<xref ref-type="bibr" rid="ref-6">6</xref>], Wu targeted controlling the amount of energy used by underground sensors to map water pipelines through optimal relay placement. In the same way, the optimum relay node location was debated for the goal of extending the network&#x0027;s duration subjected to reducing the load balance and the number of relays [<xref ref-type="bibr" rid="ref-7">7</xref>]. In [<xref ref-type="bibr" rid="ref-8">8</xref>], two approximation algorithms for relay node deployment and assignment to sensor nodes were introduced to reduce transmission loss among nodes. Since high throughput and capacity are critically constrained by nodes energy consumption, research work in WSN focus recently on studying the trade-off between spectral efficiency and energy efficiency metrics called the resource efficiency [<xref ref-type="bibr" rid="ref-9">9</xref>&#x2013;<xref ref-type="bibr" rid="ref-12">12</xref>]. The primary objective is to jointly evaluate the efficient use of a limited frequency spectrum along with energy consumption. In UWSNs, this problem was first addressed in [<xref ref-type="bibr" rid="ref-13">13</xref>], where optimal powers used by underground sources and relay nodes for data forwarding to an aboveground base station were computed to maximize the energy and spectral efficiency tradeoff. The work [<xref ref-type="bibr" rid="ref-13">13</xref>] proposes a power allocation algorithm that utilizes the Salp Swarm Algorithm (SSA) [<xref ref-type="bibr" rid="ref-14">14</xref>] to solve the considered problem since the swarm intelligence models are interesting for various computer science fields [<xref ref-type="bibr" rid="ref-15">15</xref>,<xref ref-type="bibr" rid="ref-16">16</xref>]. The SSA is suggested in [<xref ref-type="bibr" rid="ref-14">14</xref>] as a recent metaheuristic algorithm which outperforms many other metaheuristic algorithms through tests on 19 different benchmark functions. In nearby research works, SSA proves its efficiency in node localization optimization in WSN [<xref ref-type="bibr" rid="ref-17">17</xref>,<xref ref-type="bibr" rid="ref-18">18</xref>], energy consumption and lifetime optimization in WSN [<xref ref-type="bibr" rid="ref-19">19</xref>]. The work [<xref ref-type="bibr" rid="ref-13">13</xref>] proves that the SSA-based scheme offers a better resource efficiency, given similar bandwidth and battery cost resources, compared with the traditional UWSN scheme. In this paper, we propose to further enhance the resource efficiency of the UWSN considered in [<xref ref-type="bibr" rid="ref-13">13</xref>] by modifying and improving the SSA. A novel algorithm called Hybrid Chaotic Salp Swarm with Crossover (HCSSC) is proposed to determine the optimal powers required by the source and relay nodes that enhance the network resource efficiency considering the initial nodes&#x2019; battery capacities. The proposed algorithm uses chaos theory [<xref ref-type="bibr" rid="ref-20">20</xref>] to generate feasible initial solutions. This can enhance the diversity of solutions due to the randomness and dynamic features of the chaos. Moreover, to compute the final optimal solution, a uniform crossover operator [<xref ref-type="bibr" rid="ref-21">21</xref>,<xref ref-type="bibr" rid="ref-22">22</xref>] is integrated in the exploration phase of the optimization algorithm. This leads to accelerated algorithm convergence due to the wide exploration of the search space. The proposed power optimization scheme is evaluated in terms of its effect on the average relay power and battery remaining per transmission. Since the relay node can have a higher battery capacity than the source node, we propose to study the efficiency of the proposed algorithm under both equal and different node capacities.</p>
<p>The main contributions of this work are listed as follows
<list list-type="order">
<list-item><p>A new algorithm called HCSSC algorithm is proposed to improve the standard SSA by considering a logistic chaotic map to initialize the population and integrating the uniform crossover operator in the update of salps positions.</p></list-item>
<list-item><p>The resource efficiency performance is ameliorated compared to that obtained in the work [<xref ref-type="bibr" rid="ref-13">13</xref>].</p></list-item>
<list-item><p>The average consumed relay power per transmission is minimized and the average remaining relay battery per transmission is maximized.</p></list-item>
<list-item><p>The efficiency of the proposed algorithm is proved in case of both equal and different node batteries capacities.</p></list-item>
</list></p>
<p>This paper is structured as follows. In Section 2, the UWSN system model is presented and the considered problem is formulated. In Section 3, the proposed HCSSC algorithm is detailed. In Section 4, the experimental results and performance analysis are discussed, and finally, Section 5 concludes the paper.</p>
</sec>
<sec id="s2"><label>2</label><title>System Model and Problem Formulation</title>
<p>The considered UWSN model consists of sensor source node <italic>S</italic> that gathers and forwards sensory data to an aboveground base station <italic>B</italic> through a half-duplex Amplify-and-Forward (AF) relay node <inline-formula id="ieqn-1"><mml:math id="mml-ieqn-1"><mml:mspace width="thickmathspace" /><mml:mi>R</mml:mi></mml:math></inline-formula>. The communication link between <inline-formula id="ieqn-2"><mml:math id="mml-ieqn-2"><mml:mspace width="thickmathspace" /><mml:mi>S</mml:mi></mml:math></inline-formula> and <inline-formula id="ieqn-3"><mml:math id="mml-ieqn-3"><mml:mi>R</mml:mi><mml:mspace width="thickmathspace" /></mml:math></inline-formula> is an UG2UG link, whereas the communication link between <inline-formula id="ieqn-4"><mml:math id="mml-ieqn-4"><mml:mrow><mml:mtext>&#x00A0;&#xA0;</mml:mtext></mml:mrow><mml:mi>R</mml:mi><mml:mrow><mml:mtext>&#x00A0; and&#x00A0;&#xA0;</mml:mtext></mml:mrow><mml:mi>B</mml:mi><mml:mspace width="thickmathspace" /></mml:math></inline-formula> is an UG2AG link. The mathematical expressions of the path losses for UG2UG and UG2AG, the communication mechanism among nodes, and the optimization problem under consideration are described in this section.</p>
<sec id="s2_1"><label>2.1</label><title>UG2UG and UG2AG Channel Model</title>
<p>Both source and relay nodes are buried in the soil in UWSNs. Here, the sensor nodes are buried deeper than the relay nodes according to the ground surface. In <xref ref-type="fig" rid="fig-1">Fig. 1</xref>, a two-dimensional plane represents UWSN deployment where the nodes distances in the plane using Cartesian coordinates defined on the <inline-formula id="ieqn-5"><mml:math id="mml-ieqn-5"><mml:mi>x</mml:mi><mml:mspace width="thickmathspace" /></mml:math></inline-formula> and <italic>y</italic> axes with origin <italic>O</italic>. The base station is positioned at the origin of <italic>x</italic>-axis whereas the ground surface is represented as the origin of <italic>y</italic>-axis. <inline-formula id="ieqn-6"><mml:math id="mml-ieqn-6"><mml:mrow><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mi>X</mml:mi><mml:mi>G</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is the node burial depth at positions <inline-formula id="ieqn-7"><mml:math id="mml-ieqn-7"><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mi>X</mml:mi></mml:msub></mml:mrow><mml:mo>,</mml:mo><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mi>Y</mml:mi></mml:msub></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> for node <inline-formula id="ieqn-8"><mml:math id="mml-ieqn-8"><mml:mi>X</mml:mi><mml:mo>&#x2208;</mml:mo><mml:mo fence="false" stretchy="false">{</mml:mo><mml:mrow><mml:mi>S</mml:mi><mml:mo>,</mml:mo><mml:mi>R</mml:mi></mml:mrow><mml:mo fence="false" stretchy="false">}</mml:mo></mml:math></inline-formula>.</p>
<fig id="fig-1"><label>Figure 1</label><caption><title>The topology of UWSN</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_25741-fig-1.png"/></fig>
<p>According to [<xref ref-type="bibr" rid="ref-23">23</xref>], the UG2UG path loss <inline-formula id="ieqn-9"><mml:math id="mml-ieqn-9"><mml:mrow><mml:mi mathvariant="script">L</mml:mi></mml:mrow><mml:mi>u</mml:mi><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>S</mml:mi><mml:mi>R</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is defined as
<disp-formula id="eqn-1"><label>(1)</label><mml:math id="mml-eqn-1" display="block"><mml:mrow><mml:mi mathvariant="script">L</mml:mi></mml:mrow><mml:mi>u</mml:mi><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>S</mml:mi><mml:mi>R</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>6.4</mml:mn><mml:mrow><mml:mtext>&#x00A0;&#xA0;</mml:mtext></mml:mrow><mml:mo>+</mml:mo><mml:mn>20</mml:mn><mml:mrow><mml:msub><mml:mi>log</mml:mi><mml:mrow><mml:mn>10</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mrow><mml:mtext>SR</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo>+</mml:mo><mml:mn>20</mml:mn><mml:mrow><mml:msub><mml:mi>log</mml:mi><mml:mrow><mml:mn>10</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>&#x03B2;</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>+</mml:mo><mml:mrow><mml:mtext>&#x00A0;&#xA0;</mml:mtext></mml:mrow><mml:mn>8.69</mml:mn><mml:mrow><mml:mi>&#x03B1;</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mrow><mml:mtext>SR</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mn>10</mml:mn><mml:mrow><mml:msub><mml:mi>log</mml:mi><mml:mrow><mml:mn>10</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>V</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></disp-formula>where <inline-formula id="ieqn-10"><mml:math id="mml-ieqn-10"><mml:mrow><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mrow><mml:mtext>SR</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mspace width="thickmathspace" /><mml:mo>=</mml:mo><mml:msqrt><mml:mrow><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mi>S</mml:mi></mml:msub></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mi>R</mml:mi></mml:msub></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mi>S</mml:mi></mml:msub></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mi>R</mml:mi></mml:msub></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:msqrt></mml:math></inline-formula> is the distance in metres between <inline-formula id="ieqn-11"><mml:math id="mml-ieqn-11"><mml:mspace width="thickmathspace" /><mml:mi>S</mml:mi></mml:math></inline-formula> and <inline-formula id="ieqn-12"><mml:math id="mml-ieqn-12"><mml:mi>R</mml:mi><mml:mo>.</mml:mo></mml:math></inline-formula> The constants <inline-formula id="ieqn-13"><mml:math id="mml-ieqn-13"><mml:mi>&#x03B1;</mml:mi></mml:math></inline-formula> and <inline-formula id="ieqn-14"><mml:math id="mml-ieqn-14"><mml:mi>&#x03B2;</mml:mi></mml:math></inline-formula> measure the attenuation and phase shifting, respectively. The factor <italic>V</italic> represents the attenuation of the reflection path obtained when the wave is reflected by the ground surface. Hence, the UG2UG communication results from the propagation of the signal in the reflection path and in the direct path between the two sensors [<xref ref-type="bibr" rid="ref-22">22</xref>]. The authors in [<xref ref-type="bibr" rid="ref-24">24</xref>] utilize the electromagnetic propagation model to provide the detailed expressions the constants <inline-formula id="ieqn-15"><mml:math id="mml-ieqn-15"><mml:mspace width="thickmathspace" /><mml:mi>&#x03B1;</mml:mi></mml:math></inline-formula> and <inline-formula id="ieqn-16"><mml:math id="mml-ieqn-16"><mml:mi>&#x03B2;</mml:mi></mml:math></inline-formula>. The soil medium dielectric characteristics along with the system operating frequency <inline-formula id="ieqn-17"><mml:math id="mml-ieqn-17"><mml:mspace width="thickmathspace" /><mml:mi>q</mml:mi></mml:math></inline-formula>, the sand and clay percentages, the bulk density and the Volumetric Water Content constitute the main features that these constants depend on [<xref ref-type="bibr" rid="ref-25">25</xref>]. The soil and air constitute the two media throughout passes the communication between <inline-formula id="ieqn-18"><mml:math id="mml-ieqn-18"><mml:mspace width="thickmathspace" /><mml:mi>R</mml:mi></mml:math></inline-formula> and <italic>B</italic>. Yet, there is no refraction loss from under to aboveground transition due to the perpendicular propagation of signal from higher to lower density medium [<xref ref-type="bibr" rid="ref-3">3</xref>,<xref ref-type="bibr" rid="ref-26">26</xref>]. Then, based on [<xref ref-type="bibr" rid="ref-3">3</xref>,<xref ref-type="bibr" rid="ref-26">26</xref>], the path loss <inline-formula id="ieqn-19"><mml:math id="mml-ieqn-19"><mml:mrow><mml:mi mathvariant="script">L</mml:mi></mml:mrow><mml:mi>u</mml:mi><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>R</mml:mi><mml:mi>B</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> of the channel between <italic>R</italic> and <italic>B</italic> is the result of adding the path losses for both underground and aboveground portions as in <xref ref-type="disp-formula" rid="eqn-2">Eq. (2)</xref>
<disp-formula id="eqn-2"><label>(2)</label><mml:math id="mml-eqn-2" display="block"><mml:mrow><mml:mi mathvariant="script">L</mml:mi></mml:mrow><mml:mi>u</mml:mi><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>R</mml:mi><mml:mi>B</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mi mathvariant="script">L</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>R</mml:mi><mml:mi>G</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:mi mathvariant="script">L</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>G</mml:mi><mml:mi>B</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>,</mml:mo></mml:math></disp-formula>
<disp-formula id="eqn-3"><label>(2.1)</label><mml:math id="mml-eqn-3" display="block"><mml:mrow><mml:mi mathvariant="script">L</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>R</mml:mi><mml:mi>G</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>6.4</mml:mn><mml:mo>+</mml:mo><mml:mn>20</mml:mn><mml:mrow><mml:msub><mml:mi>log</mml:mi><mml:mrow><mml:mn>10</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mrow><mml:mtext>RG</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo>+</mml:mo><mml:mn>20</mml:mn><mml:mrow><mml:msub><mml:mi>log</mml:mi><mml:mrow><mml:mn>10</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>&#x03B2;</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>+</mml:mo><mml:mrow><mml:mtext>&#x00A0;&#xA0;</mml:mtext></mml:mrow><mml:mn>8.69</mml:mn><mml:mrow><mml:mi>&#x03B1;</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mrow><mml:mtext>RG</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>,</mml:mo></mml:math></disp-formula>
<disp-formula id="eqn-4"><label>(2.2)</label><mml:math id="mml-eqn-4" display="block"><mml:mrow><mml:mi mathvariant="script">L</mml:mi></mml:mrow><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>G</mml:mi><mml:mi>B</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>,</mml:mo><mml:mo>=</mml:mo><mml:mo>&#x2212;</mml:mo><mml:mn>147.6</mml:mn><mml:mrow><mml:mtext>&#x00A0;&#xA0;</mml:mtext></mml:mrow><mml:mo>+</mml:mo><mml:mn>10</mml:mn><mml:mi>&#x03B7;</mml:mi><mml:mspace width="thickmathspace" /><mml:mrow><mml:msub><mml:mi>log</mml:mi><mml:mrow><mml:mn>10</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mrow><mml:mtext>BG</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo>+</mml:mo><mml:mn>20</mml:mn><mml:mrow><mml:msub><mml:mi>log</mml:mi><mml:mrow><mml:mn>10</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>q</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>,</mml:mo></mml:math></disp-formula>
<disp-formula id="eqn-5"><label>(2.3)</label><mml:math id="mml-eqn-5" display="block"><mml:mrow><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mrow><mml:mtext>BG</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:msqrt><mml:msubsup><mml:mi>d</mml:mi><mml:mrow><mml:mi>O</mml:mi><mml:mi>R</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup><mml:mo>+</mml:mo><mml:msubsup><mml:mi>d</mml:mi><mml:mrow><mml:mi>O</mml:mi><mml:mi>B</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup></mml:msqrt><mml:mo>,</mml:mo></mml:math></disp-formula>where <inline-formula id="ieqn-20"><mml:math id="mml-ieqn-20"><mml:mrow><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mi>O</mml:mi><mml:mi>R</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> represents the horizontal distance between <italic>R</italic> and the origin <italic>O</italic>, <inline-formula id="ieqn-21"><mml:math id="mml-ieqn-21"><mml:mrow><mml:msub><mml:mi>d</mml:mi><mml:mrow><mml:mi>O</mml:mi><mml:mi>B</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is the height of the aboveground base station <italic>B</italic> and <inline-formula id="ieqn-22"><mml:math id="mml-ieqn-22"><mml:mi>&#x03B7;</mml:mi></mml:math></inline-formula> is the attenuation coefficient on air [<xref ref-type="bibr" rid="ref-26">26</xref>].</p>
</sec>
<sec id="s2_2"><label>2.2</label><title>UWSN System Model</title>
<p>Here, the uplink communication procedure among the trio link is presented. The Time Division Multiple Access (TDMA) scheme is considered to mitigate the signal interference. Each node <inline-formula id="ieqn-23"><mml:math id="mml-ieqn-23"><mml:mi>X</mml:mi><mml:mo>&#x2208;</mml:mo><mml:mo fence="false" stretchy="false">{</mml:mo><mml:mrow><mml:mi>S</mml:mi><mml:mo>,</mml:mo><mml:mi>R</mml:mi></mml:mrow><mml:mo fence="false" stretchy="false">}</mml:mo></mml:math></inline-formula> has a limited battery with power capacity <inline-formula id="ieqn-24"><mml:math id="mml-ieqn-24"><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mi>X</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>. At each packet transmission, <italic>t</italic>, node <italic>X</italic> consumes a power <inline-formula id="ieqn-25"><mml:math id="mml-ieqn-25"><mml:msubsup><mml:mi>P</mml:mi><mml:mi>X</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mspace width="thickmathspace" /></mml:mrow></mml:msubsup><mml:mo>&#x2208;</mml:mo><mml:mo stretchy="false">[</mml:mo><mml:mrow><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>,</mml:mo><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:mrow><mml:mo stretchy="false">]</mml:mo></mml:math></inline-formula>. The communication process at each transmission <inline-formula id="ieqn-26"><mml:math id="mml-ieqn-26"><mml:mi>t</mml:mi><mml:mo>&#x2208;</mml:mo><mml:mo stretchy="false">[</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mi>&#x03C4;</mml:mi></mml:mrow><mml:mo stretchy="false">]</mml:mo></mml:math></inline-formula> needs two phases.</p>
<p>In the first phase, the source node <italic>S</italic> transmits a data packet <inline-formula id="ieqn-27"><mml:math id="mml-ieqn-27"><mml:mrow><mml:msup><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow><mml:mspace width="thickmathspace" /></mml:math></inline-formula> to the relay node <italic>R</italic>. The resulting received signal <inline-formula id="ieqn-28"><mml:math id="mml-ieqn-28"><mml:msubsup><mml:mrow><mml:mtext>y</mml:mtext></mml:mrow><mml:mi>R</mml:mi><mml:mi>t</mml:mi></mml:msubsup></mml:math></inline-formula> at <italic>R</italic> is given by
<disp-formula id="eqn-6"><label>(3)</label><mml:math id="mml-eqn-6" display="block"><mml:msubsup><mml:mrow><mml:mtext>y</mml:mtext></mml:mrow><mml:mi>R</mml:mi><mml:mi>t</mml:mi></mml:msubsup><mml:mo>=</mml:mo><mml:msqrt><mml:msubsup><mml:mi>P</mml:mi><mml:mi>S</mml:mi><mml:mi>t</mml:mi></mml:msubsup></mml:msqrt><mml:msubsup><mml:mi>h</mml:mi><mml:mrow><mml:mi>S</mml:mi><mml:mi>R</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup><mml:mrow><mml:msup><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow><mml:mo>+</mml:mo><mml:msubsup><mml:mrow><mml:mtext>n</mml:mtext></mml:mrow><mml:mi>R</mml:mi><mml:mi>t</mml:mi></mml:msubsup></mml:math></disp-formula>such that <inline-formula id="ieqn-29"><mml:math id="mml-ieqn-29"><mml:mspace width="thickmathspace" /><mml:msubsup><mml:mi>h</mml:mi><mml:mrow><mml:mi>S</mml:mi><mml:mi>R</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup></mml:math></inline-formula> represents the gain of the UG2UG channel between the two nodes <inline-formula id="ieqn-30"><mml:math id="mml-ieqn-30"><mml:mspace width="thickmathspace" /><mml:mi>S</mml:mi></mml:math></inline-formula> and <italic>R</italic>, which obeys to the Rayleigh distribution [<xref ref-type="bibr" rid="ref-23">23</xref>] with the underground path loss <inline-formula id="ieqn-31"><mml:math id="mml-ieqn-31"><mml:mspace width="thickmathspace" /><mml:mrow><mml:mi mathvariant="script">L</mml:mi></mml:mrow><mml:mi>u</mml:mi><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>S</mml:mi><mml:mi>R</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> calculated in <xref ref-type="disp-formula" rid="eqn-1">Eq. (1)</xref> and <inline-formula id="ieqn-32"><mml:math id="mml-ieqn-32"><mml:msubsup><mml:mrow><mml:mtext>n</mml:mtext></mml:mrow><mml:mi>R</mml:mi><mml:mi>t</mml:mi></mml:msubsup><mml:mspace width="thickmathspace" /></mml:math></inline-formula> is the zero-mean complex Additive White Gaussian Noise (AWGN) vector with variance <inline-formula id="ieqn-33"><mml:math id="mml-ieqn-33"><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>. Hence, the Signal-to-Noise Ratio (SNR) between <italic>S</italic> and <italic>R</italic> designated by <inline-formula id="ieqn-34"><mml:math id="mml-ieqn-34"><mml:msubsup><mml:mi>&#x0393;</mml:mi><mml:mrow><mml:mi>S</mml:mi><mml:mi>R</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup></mml:math></inline-formula>, resulting of the communication between <italic>S</italic> and <italic>R</italic> at transmission <italic>t</italic>, is as follows
<disp-formula id="eqn-7"><label>(4)</label><mml:math id="mml-eqn-7" display="block"><mml:msubsup><mml:mi>&#x0393;</mml:mi><mml:mrow><mml:mi>S</mml:mi><mml:mi>R</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msubsup><mml:mi>P</mml:mi><mml:mi>S</mml:mi><mml:mi>t</mml:mi></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo stretchy="false">|</mml:mo></mml:mrow><mml:mrow><mml:msubsup><mml:mi>h</mml:mi><mml:mrow><mml:mi>S</mml:mi><mml:mi>R</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup></mml:mrow><mml:mrow><mml:mo stretchy="false">|</mml:mo></mml:mrow></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mrow><mml:mrow><mml:mi>b</mml:mi><mml:mspace width="thickmathspace" /><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow></mml:mrow></mml:mfrac></mml:math></disp-formula></p>
<p>Where <inline-formula id="ieqn-35"><mml:math id="mml-ieqn-35"><mml:mspace width="thickmathspace" /><mml:mi>b</mml:mi></mml:math></inline-formula> denotes the channel bandwidth in Hz, which is equal to the operating frequency <italic>q</italic> when TDMA is applied. In the second phase, the signal <inline-formula id="ieqn-36"><mml:math id="mml-ieqn-36"><mml:msubsup><mml:mrow><mml:mtext>y</mml:mtext></mml:mrow><mml:mi>R</mml:mi><mml:mi>t</mml:mi></mml:msubsup></mml:math></inline-formula> is amplified and forwarded by <italic>R</italic> to <italic>B</italic>. Consequently, the received signal <inline-formula id="ieqn-37"><mml:math id="mml-ieqn-37"><mml:msubsup><mml:mrow><mml:mtext>y</mml:mtext></mml:mrow><mml:mi>B</mml:mi><mml:mi>t</mml:mi></mml:msubsup></mml:math></inline-formula> at <italic>B</italic> is calculated as follows.
<disp-formula id="eqn-8"><label>(5)</label><mml:math id="mml-eqn-8" display="block"><mml:msubsup><mml:mrow><mml:mtext>y</mml:mtext></mml:mrow><mml:mi>B</mml:mi><mml:mi>t</mml:mi></mml:msubsup><mml:mo>=</mml:mo><mml:msubsup><mml:mi>A</mml:mi><mml:mi>f</mml:mi><mml:mi>t</mml:mi></mml:msubsup><mml:msqrt><mml:msubsup><mml:mi>P</mml:mi><mml:mi>R</mml:mi><mml:mi>t</mml:mi></mml:msubsup></mml:msqrt><mml:msubsup><mml:mi>h</mml:mi><mml:mrow><mml:mi>R</mml:mi><mml:mi>B</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup><mml:msubsup><mml:mrow><mml:mtext>y</mml:mtext></mml:mrow><mml:mi>R</mml:mi><mml:mi>t</mml:mi></mml:msubsup><mml:mo>+</mml:mo><mml:msubsup><mml:mrow><mml:mtext>n</mml:mtext></mml:mrow><mml:mi>B</mml:mi><mml:mi>t</mml:mi></mml:msubsup></mml:math></disp-formula></p>
<p>With <inline-formula id="ieqn-38"><mml:math id="mml-ieqn-38"><mml:msubsup><mml:mi>h</mml:mi><mml:mrow><mml:mi>R</mml:mi><mml:mi>B</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup></mml:math></inline-formula> is the UG2AG Rayleigh distributed channel between <inline-formula id="ieqn-39"><mml:math id="mml-ieqn-39"><mml:mspace width="thickmathspace" /><mml:mi>R</mml:mi></mml:math></inline-formula> and <italic>B</italic> [<xref ref-type="bibr" rid="ref-23">23</xref>] with the underground path loss <inline-formula id="ieqn-40"><mml:math id="mml-ieqn-40"><mml:mrow><mml:mi mathvariant="script">L</mml:mi></mml:mrow><mml:mi>u</mml:mi><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>R</mml:mi><mml:mi>B</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> calculated in <xref ref-type="disp-formula" rid="eqn-2">Eq. (2)</xref>,<inline-formula id="ieqn-41"><mml:math id="mml-ieqn-41"><mml:mspace width="thickmathspace" /><mml:msubsup><mml:mrow><mml:mtext>n</mml:mtext></mml:mrow><mml:mi>B</mml:mi><mml:mi>t</mml:mi></mml:msubsup></mml:math></inline-formula> is the zero-mean complex AWGN vector and <inline-formula id="ieqn-42"><mml:math id="mml-ieqn-42"><mml:msubsup><mml:mi>A</mml:mi><mml:mi>f</mml:mi><mml:mi>t</mml:mi></mml:msubsup><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:msqrt><mml:msubsup><mml:mi>P</mml:mi><mml:mi>S</mml:mi><mml:mi>t</mml:mi></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo stretchy="false">|</mml:mo></mml:mrow><mml:mrow><mml:msubsup><mml:mi>h</mml:mi><mml:mrow><mml:mi>S</mml:mi><mml:mi>R</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup></mml:mrow><mml:mrow><mml:mo stretchy="false">|</mml:mo></mml:mrow></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mo>+</mml:mo><mml:mspace width="thickmathspace" /><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow></mml:msqrt></mml:mrow></mml:mfrac></mml:math></inline-formula> is the amplification factor. Thus, the SNR, denoted by <inline-formula id="ieqn-43"><mml:math id="mml-ieqn-43"><mml:msubsup><mml:mi>&#x0393;</mml:mi><mml:mrow><mml:mi>R</mml:mi><mml:mi>B</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup></mml:math></inline-formula>, of the transmission between <italic>R</italic> and <italic>B</italic> is computed as follows.
<disp-formula id="eqn-9"><label>(6)</label><mml:math id="mml-eqn-9" display="block"><mml:msubsup><mml:mi>&#x0393;</mml:mi><mml:mrow><mml:mi>R</mml:mi><mml:mi>B</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>A</mml:mi><mml:mrow><mml:msup><mml:mrow><mml:msubsup><mml:mi></mml:mi><mml:mi>f</mml:mi><mml:mi>t</mml:mi></mml:msubsup></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:msubsup><mml:mi>P</mml:mi><mml:mi>R</mml:mi><mml:mi>t</mml:mi></mml:msubsup><mml:msubsup><mml:mi>P</mml:mi><mml:mi>S</mml:mi><mml:mi>t</mml:mi></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo stretchy="false">|</mml:mo></mml:mrow><mml:mrow><mml:msubsup><mml:mi>h</mml:mi><mml:mrow><mml:mi>R</mml:mi><mml:mi>B</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup></mml:mrow><mml:mrow><mml:mo stretchy="false">|</mml:mo></mml:mrow></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo stretchy="false">|</mml:mo></mml:mrow><mml:mrow><mml:msubsup><mml:mi>h</mml:mi><mml:mrow><mml:mi>S</mml:mi><mml:mi>R</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup></mml:mrow><mml:mrow><mml:mo stretchy="false">|</mml:mo></mml:mrow></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mrow><mml:mrow><mml:mi>b</mml:mi><mml:mspace width="thickmathspace" /><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>A</mml:mi><mml:mrow><mml:msup><mml:mrow><mml:msubsup><mml:mi></mml:mi><mml:mi>f</mml:mi><mml:mi>t</mml:mi></mml:msubsup></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:msubsup><mml:mi>P</mml:mi><mml:mi>R</mml:mi><mml:mi>t</mml:mi></mml:msubsup><mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo stretchy="false">|</mml:mo></mml:mrow><mml:mrow><mml:msubsup><mml:mi>h</mml:mi><mml:mrow><mml:mi>R</mml:mi><mml:mi>B</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup></mml:mrow><mml:mrow><mml:mo stretchy="false">|</mml:mo></mml:mrow></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mfrac></mml:math></disp-formula></p>
<p>With respect to <xref ref-type="disp-formula" rid="eqn-7">Eqs. (4)</xref> and <xref ref-type="disp-formula" rid="eqn-9">(6)</xref>, the total SNR <inline-formula id="ieqn-44"><mml:math id="mml-ieqn-44"><mml:msubsup><mml:mi>&#x0393;</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mi>o</mml:mi><mml:mi>t</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup></mml:math></inline-formula> and the maximum data rate <inline-formula id="ieqn-45"><mml:math id="mml-ieqn-45"><mml:msubsup><mml:mi>R</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mi>o</mml:mi><mml:mi>t</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup></mml:math></inline-formula> for the trio link are given respectively by
<disp-formula id="eqn-10"><label>(7)</label><mml:math id="mml-eqn-10" display="block"><mml:msubsup><mml:mi>&#x0393;</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mi>o</mml:mi><mml:mi>t</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msubsup><mml:mi>&#x0393;</mml:mi><mml:mrow><mml:mi>S</mml:mi><mml:mi>R</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup><mml:msubsup><mml:mi>&#x0393;</mml:mi><mml:mrow><mml:mi>R</mml:mi><mml:mi>B</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup></mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:msubsup><mml:mi>&#x0393;</mml:mi><mml:mrow><mml:mi>S</mml:mi><mml:mi>R</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup><mml:mo>+</mml:mo><mml:msubsup><mml:mi>&#x0393;</mml:mi><mml:mrow><mml:mi>R</mml:mi><mml:mi>B</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup><mml:mspace width="thickmathspace" /></mml:mrow></mml:mfrac></mml:math></disp-formula>
<disp-formula id="eqn-11"><label>(8)</label><mml:math id="mml-eqn-11" display="block"><mml:msubsup><mml:mi>R</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mi>o</mml:mi><mml:mi>t</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup><mml:mo>=</mml:mo><mml:mi>b</mml:mi><mml:mrow><mml:msub><mml:mi>log</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:msubsup><mml:mi mathvariant="normal">&#x0393;</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mi>o</mml:mi><mml:mi>t</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:math></disp-formula></p>
</sec>
<sec id="s2_3"><label>2.3</label><title>Problem Formulation</title>
<p>The overall paper objective is to determine the power allocation vector <inline-formula id="ieqn-46"><mml:math id="mml-ieqn-46"><mml:mspace width="thickmathspace" /><mml:mi>P</mml:mi><mml:mo>=</mml:mo><mml:mo stretchy="false">[</mml:mo><mml:mrow><mml:msubsup><mml:mi>P</mml:mi><mml:mi>S</mml:mi><mml:mi>t</mml:mi></mml:msubsup><mml:mo>,</mml:mo><mml:msubsup><mml:mi>P</mml:mi><mml:mi>R</mml:mi><mml:mi>t</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy="false">]</mml:mo></mml:math></inline-formula> allowing an efficient balance between two competing metrics: the energy efficiency <inline-formula id="ieqn-47"><mml:math id="mml-ieqn-47"><mml:mi>E</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> and spectral efficiency <inline-formula id="ieqn-48"><mml:math id="mml-ieqn-48"><mml:mi>S</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> at each transmission <italic>t</italic>. Indeed, <inline-formula id="ieqn-49"><mml:math id="mml-ieqn-49"><mml:mi>E</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula>. is the total bits generated per unit energy whereas <inline-formula id="ieqn-50"><mml:math id="mml-ieqn-50"><mml:mi>S</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> represents the total delivered bits per unit bandwidth of link <inline-formula id="ieqn-51"><mml:math id="mml-ieqn-51"><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>S</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>R</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>B</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mspace width="thickmathspace" /></mml:math></inline-formula> at transmission <italic>t</italic>. Their expressions are respectively as follows
<disp-formula id="eqn-12"><label>(9)</label><mml:math id="mml-eqn-12" display="block"><mml:mi>E</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msubsup><mml:mi>R</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mi>o</mml:mi><mml:mi>t</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup></mml:mrow><mml:mrow><mml:msubsup><mml:mi>P</mml:mi><mml:mi>S</mml:mi><mml:mi>t</mml:mi></mml:msubsup><mml:mo>+</mml:mo><mml:msubsup><mml:mi>P</mml:mi><mml:mi>R</mml:mi><mml:mi>t</mml:mi></mml:msubsup></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:math></disp-formula>
<disp-formula id="eqn-13"><label>(10)</label><mml:math id="mml-eqn-13" display="block"><mml:mi>S</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msubsup><mml:mi>R</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mi>o</mml:mi><mml:mi>t</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:msubsup></mml:mrow><mml:mi>b</mml:mi></mml:mfrac></mml:math></disp-formula></p>
<p>As discussed in [<xref ref-type="bibr" rid="ref-9">9</xref>], the resource efficiency metric <inline-formula id="ieqn-52"><mml:math id="mml-ieqn-52"><mml:mspace width="thickmathspace" /><mml:mi>R</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow><mml:mspace width="thickmathspace" /></mml:math></inline-formula> is capable to exploit the trade-off between <inline-formula id="ieqn-53"><mml:math id="mml-ieqn-53"><mml:mi>E</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-54"><mml:math id="mml-ieqn-54"><mml:mi>S</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula>. It is given by
<disp-formula id="eqn-14"><label>(11)</label><mml:math id="mml-eqn-14" display="block"><mml:mi>R</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow><mml:mo>=</mml:mo><mml:mi>E</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow><mml:mo>+</mml:mo><mml:mi>&#x03C9;</mml:mi><mml:mi>S</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></disp-formula>where the weighted factor <inline-formula id="ieqn-55"><mml:math id="mml-ieqn-55"><mml:mi>&#x03C9;</mml:mi></mml:math></inline-formula> is computed as <inline-formula id="ieqn-56"><mml:math id="mml-ieqn-56"><mml:mspace width="thickmathspace" /><mml:mi>&#x03C9;</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mover><mml:mi>&#x03C9;</mml:mi><mml:mo stretchy="false">&#x00AF;</mml:mo></mml:mover></mml:mrow><mml:mfrac><mml:mi>b</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mi>o</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:mfrac></mml:math></inline-formula>, with <inline-formula id="ieqn-57"><mml:math id="mml-ieqn-57"><mml:mrow><mml:mover><mml:mi>&#x03C9;</mml:mi><mml:mo stretchy="false">&#x00AF;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula> is a constant and <inline-formula id="ieqn-58"><mml:math id="mml-ieqn-58"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mi>o</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>+</mml:mo><mml:mspace width="thickmathspace" /><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is the total power budget allocated to the link <inline-formula id="ieqn-59"><mml:math id="mml-ieqn-59"><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>S</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>R</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>B</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mspace width="thickmathspace" /></mml:math></inline-formula> at each transmission. Since, <inline-formula id="ieqn-60"><mml:math id="mml-ieqn-60"><mml:mi>S</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> is practically smaller than <inline-formula id="ieqn-61"><mml:math id="mml-ieqn-61"><mml:mi>E</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula>, this weighted factor <inline-formula id="ieqn-62"><mml:math id="mml-ieqn-62"><mml:mi>&#x03C9;</mml:mi></mml:math></inline-formula> is utilized to achieve the balance between <inline-formula id="ieqn-63"><mml:math id="mml-ieqn-63"><mml:mi>E</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-64"><mml:math id="mml-ieqn-64"><mml:mi>S</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow><mml:mo>.</mml:mo></mml:math></inline-formula> In fact, it solves the inconsistence of adding two metrics with different units since the unit of <inline-formula id="ieqn-65"><mml:math id="mml-ieqn-65"><mml:mi>E</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> is bits/Joule while the unit of <inline-formula id="ieqn-66"><mml:math id="mml-ieqn-66"><mml:mi>S</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> is bits/s/Hz. Then, the unit of <inline-formula id="ieqn-67"><mml:math id="mml-ieqn-67"><mml:mspace width="thickmathspace" /><mml:mi>R</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> is equivalent to <inline-formula id="ieqn-68"><mml:math id="mml-ieqn-68"><mml:mi>E</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> which is still bits/Joule. Besides, maximizing <inline-formula id="ieqn-69"><mml:math id="mml-ieqn-69"><mml:mspace width="thickmathspace" /><mml:mi>R</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> is simply maximizing <inline-formula id="ieqn-70"><mml:math id="mml-ieqn-70"><mml:mi>E</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> if <inline-formula id="ieqn-71"><mml:math id="mml-ieqn-71"><mml:mrow><mml:mover><mml:mi>&#x03C9;</mml:mi><mml:mo stretchy="false">&#x00AF;</mml:mo></mml:mover></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:math></inline-formula> and <inline-formula id="ieqn-72"><mml:math id="mml-ieqn-72"><mml:mi>S</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> if <inline-formula id="ieqn-73"><mml:math id="mml-ieqn-73"><mml:mrow><mml:mover><mml:mi>&#x03C9;</mml:mi><mml:mo stretchy="false">&#x00AF;</mml:mo></mml:mover></mml:mrow><mml:mo>=</mml:mo><mml:mi mathvariant="normal">&#x221E;</mml:mi></mml:math></inline-formula>. The selection of this constant value in the UWSN model is described in [<xref ref-type="bibr" rid="ref-13">13</xref>].</p>
<p>For each node <italic>X</italic>, optimizing the power <inline-formula id="ieqn-74"><mml:math id="mml-ieqn-74"><mml:msubsup><mml:mi>P</mml:mi><mml:mi>X</mml:mi><mml:mi>t</mml:mi></mml:msubsup></mml:math></inline-formula> at each transmission <italic>t</italic> is a mandatory requirement because it depends on the limited battery capacity <inline-formula id="ieqn-75"><mml:math id="mml-ieqn-75"><mml:mrow><mml:msub><mml:mrow><mml:mtext>C</mml:mtext></mml:mrow><mml:mi>X</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> of the buried node, the spent powers for transmissions <inline-formula id="ieqn-76"><mml:math id="mml-ieqn-76"><mml:mo stretchy="false">[</mml:mo><mml:mrow><mml:msubsup><mml:mi>P</mml:mi><mml:mi>X</mml:mi><mml:mn>1</mml:mn></mml:msubsup><mml:mo>,</mml:mo><mml:mspace width="thickmathspace" /><mml:mo>&#x2026;</mml:mo><mml:mo>,</mml:mo><mml:msubsup><mml:mi>P</mml:mi><mml:mi>X</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup></mml:mrow><mml:mo stretchy="false">]</mml:mo><mml:mspace width="thickmathspace" /></mml:math></inline-formula> and the allowed power limitation range at each transmission <inline-formula id="ieqn-77"><mml:math id="mml-ieqn-77"><mml:mo stretchy="false">[</mml:mo><mml:mrow><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>,</mml:mo><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:mrow><mml:mo stretchy="false">]</mml:mo></mml:math></inline-formula>. Then, the proposed optimization problem is stated as
<disp-formula id="eqn-15"><label>(12)</label><mml:math id="mml-eqn-15" display="block"><mml:mrow><mml:mtext>Maximize</mml:mtext></mml:mrow><mml:mspace width="thickmathspace" /><mml:mspace width="thinmathspace" /><mml:mi>R</mml:mi><mml:mi>E</mml:mi><mml:mo>=</mml:mo><mml:munderover><mml:mrow><mml:mo movablelimits="false">&#x2211;</mml:mo></mml:mrow><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>&#x03C4;</mml:mi></mml:munderover><mml:mspace width="thinmathspace" /><mml:mi>R</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></disp-formula></p>
<p>Subject to <inline-formula id="ieqn-78"><mml:math id="mml-ieqn-78"><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:munderover><mml:mrow><mml:mo movablelimits="false">&#x2211;</mml:mo></mml:mrow><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>&#x03C4;</mml:mi></mml:munderover><mml:mo>&#x2061;</mml:mo><mml:msubsup><mml:mi>P</mml:mi><mml:mi>X</mml:mi><mml:mi>t</mml:mi></mml:msubsup></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x2264;</mml:mo><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mi>X</mml:mi></mml:msub></mml:mrow><mml:mo>,</mml:mo></mml:math></inline-formula></p>
<p>where <inline-formula id="ieqn-79"><mml:math id="mml-ieqn-79"><mml:msubsup><mml:mi>P</mml:mi><mml:mi>X</mml:mi><mml:mi>t</mml:mi></mml:msubsup><mml:mo>&#x2208;</mml:mo><mml:mo stretchy="false">[</mml:mo><mml:mrow><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>,</mml:mo><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:mrow><mml:mo stretchy="false">]</mml:mo><mml:mspace width="thickmathspace" /></mml:math></inline-formula> for <inline-formula id="ieqn-80"><mml:math id="mml-ieqn-80"><mml:mi>X</mml:mi><mml:mo>&#x2208;</mml:mo><mml:mo stretchy="false">[</mml:mo><mml:mrow><mml:mi>S</mml:mi><mml:mo>,</mml:mo><mml:mi>R</mml:mi></mml:mrow><mml:mo stretchy="false">]</mml:mo></mml:math></inline-formula>.</p>
<p>In <xref ref-type="disp-formula" rid="eqn-15">Eq. (12)</xref>, the maximization problem presented is considered a NP-hard problem that requires an efficient optimization algorithm to solve it. Therefore, a hybrid meta-heuristic algorithm based on SSA is proposed to obtain optimal values of nodes powers considering the resource efficiency <inline-formula id="ieqn-81"><mml:math id="mml-ieqn-81"><mml:mi>R</mml:mi><mml:mi>E</mml:mi></mml:math></inline-formula> to be maximized.</p>
</sec>
</sec>
<sec id="s3"><label>3</label><title>The Proposed HCSSC for Resource Efficiency</title>
<p>In this section, the main structure of the standard Salp Swarm Algorithm (SSA) is first reviewed. Then, the detailed steps of the proposed Hybrid Chaotic Salp Swarm with Crossover (HCSSC) scheme for maximizing the resource efficiency in the considered UWSN are addressed.</p>
<sec id="s3_1"><label>3.1</label><title>Salp Swarm Algorithm (SSA)</title>
<p>The Salp Swarm Algorithm (SSA) is one of the recent swam algorithms proposed in 2017 [<xref ref-type="bibr" rid="ref-14">14</xref>] and widely used in solving many optimization problems [<xref ref-type="bibr" rid="ref-27">27</xref>,<xref ref-type="bibr" rid="ref-28">28</xref>]. SSA emulates the motion of Salpidae that have a transparent barrel-shaped body and live-in deep oceans [<xref ref-type="bibr" rid="ref-29">29</xref>]. Salps are organized in a form of swarm called salp chain. Mathematically, the salp chain is divided into two groups: leader (i.e., the first salp of the chain and followers (i.e., the remaining salps of the chain which follow the leader). Researchers viewed that their searching for food is an indicator to their behavior [<xref ref-type="bibr" rid="ref-30">30</xref>].</p>
<p>The leader in the swarm updates its position relative to the food source <inline-formula id="ieqn-82"><mml:math id="mml-ieqn-82"><mml:mspace width="thickmathspace" /><mml:mi>F</mml:mi><mml:mspace width="thickmathspace" /></mml:math></inline-formula> according to <xref ref-type="disp-formula" rid="eqn-16">Eq. (13)</xref>
<disp-formula id="eqn-16"><label>(13)</label><mml:math id="mml-eqn-16" display="block"><mml:msubsup><mml:mi>x</mml:mi><mml:mi>j</mml:mi><mml:mn>1</mml:mn></mml:msubsup><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:mrow><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mi>j</mml:mi></mml:msub></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>u</mml:mi><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mi>j</mml:mi></mml:msub></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mi>l</mml:mi><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mi>j</mml:mi></mml:msub></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mo>+</mml:mo><mml:mi>l</mml:mi><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mi>j</mml:mi></mml:msub></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow><mml:mo>&#x2265;</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mi>j</mml:mi></mml:msub></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>u</mml:mi><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mi>j</mml:mi></mml:msub></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mi>l</mml:mi><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mi>j</mml:mi></mml:msub></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow><mml:mo>+</mml:mo><mml:mi>l</mml:mi><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mi>j</mml:mi></mml:msub></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow><mml:mo>&#x003C;</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>Where <inline-formula id="ieqn-83"><mml:math id="mml-ieqn-83"><mml:msubsup><mml:mi>x</mml:mi><mml:mi>j</mml:mi><mml:mn>1</mml:mn></mml:msubsup></mml:math></inline-formula> is the position of the leader in the <italic>j<sup>th</sup></italic> dimension, <inline-formula id="ieqn-84"><mml:math id="mml-ieqn-84"><mml:mi>l</mml:mi><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mi>j</mml:mi></mml:msub></mml:mrow><mml:mspace width="thickmathspace" /></mml:math></inline-formula> and <inline-formula id="ieqn-85"><mml:math id="mml-ieqn-85"><mml:mspace width="thickmathspace" /><mml:mi>u</mml:mi><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mi>j</mml:mi></mml:msub></mml:mrow><mml:mspace width="thickmathspace" /></mml:math></inline-formula> are the lower and the upper bound of the <italic>j<sup>th</sup></italic> dimension, respectively, and <inline-formula id="ieqn-86"><mml:math id="mml-ieqn-86"><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mi>j</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the position of food source in the <italic>j<sup>th</sup></italic> dimension. The parameter <inline-formula id="ieqn-87"><mml:math id="mml-ieqn-87"><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> balances the scale between exploration and exploitation and is computed according to <xref ref-type="disp-formula" rid="eqn-17">Eq. (14)</xref>. The coefficients <inline-formula id="ieqn-88"><mml:math id="mml-ieqn-88"><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-89"><mml:math id="mml-ieqn-89"><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> are random numbers between [&#x2212;1, 1]. Also, <inline-formula id="ieqn-90"><mml:math id="mml-ieqn-90"><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> determines the direction of moving the leader towards a positive infinity or negative infinity.
<disp-formula id="eqn-17"><label>(14)</label><mml:math id="mml-eqn-17" display="block"><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>2</mml:mn><mml:mrow><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mfrac><mml:mrow><mml:mn>4</mml:mn><mml:mi>l</mml:mi></mml:mrow><mml:mi>L</mml:mi></mml:mfrac><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:mrow></mml:msup></mml:mrow></mml:math></disp-formula>where <italic>l</italic> is the current iteration and <italic>L</italic> is the maximum number of iterations.</p>
<p>The remaining of the salps in the chain (i.e., followers) update their positions based on the Newton&#x0027;s law of motion as in <xref ref-type="disp-formula" rid="eqn-18">Eq. (15)</xref>
<disp-formula id="eqn-18"><label>(15)</label><mml:math id="mml-eqn-18" display="block"><mml:msubsup><mml:mi>x</mml:mi><mml:mi>j</mml:mi><mml:mi>i</mml:mi></mml:msubsup><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac><mml:mi>a</mml:mi><mml:mrow><mml:msup><mml:mi>t</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow><mml:mi>c</mml:mi><mml:mo>+</mml:mo><mml:mrow><mml:msub><mml:mi>v</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow><mml:mi>t</mml:mi><mml:mspace width="1em" /><mml:mi mathvariant="normal">&#x2200;</mml:mi><mml:mspace width="thickmathspace" /><mml:mi>i</mml:mi><mml:mo>&#x2265;</mml:mo><mml:mn>2</mml:mn></mml:math></disp-formula>
<disp-formula id="eqn-19"><label>(16)</label><mml:math id="mml-eqn-19" display="block"><mml:mrow><mml:mtext>where</mml:mtext></mml:mrow><mml:mspace width="thinmathspace" /><mml:mi>a</mml:mi><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mrow><mml:msub><mml:mi>v</mml:mi><mml:mrow><mml:mi>f</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi><mml:mi>a</mml:mi><mml:mi>l</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>v</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow></mml:mrow></mml:mfrac><mml:mspace width="thinmathspace" /><mml:mrow><mml:mtext>and</mml:mtext></mml:mrow><mml:mspace width="thinmathspace" /><mml:mi>v</mml:mi><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>x</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow></mml:mrow><mml:mi>t</mml:mi></mml:mfrac></mml:math></disp-formula></p>
<p>With <inline-formula id="ieqn-91"><mml:math id="mml-ieqn-91"><mml:msubsup><mml:mi>x</mml:mi><mml:mi>j</mml:mi><mml:mn>1</mml:mn></mml:msubsup></mml:math></inline-formula> is the <italic>i<sup>th</sup></italic> salp in the <italic>j<sup>th</sup></italic> dimension,<inline-formula id="ieqn-92"><mml:math id="mml-ieqn-92"><mml:mi>t</mml:mi></mml:math></inline-formula> is the time,<inline-formula id="ieqn-93"><mml:math id="mml-ieqn-93"><mml:mrow><mml:msub><mml:mi>v</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> is the initial speed.</p>
<p>In the optimization problems, the time is considered as an iteration where the conflict between iterations is equal to 1 and <inline-formula id="ieqn-94"><mml:math id="mml-ieqn-94"><mml:mrow><mml:msub><mml:mi>v</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> is set to 0. Then, <xref ref-type="disp-formula" rid="eqn-18">Eq. (15)</xref> can be reformulated in <xref ref-type="disp-formula" rid="eqn-20">Eq. (17)</xref> for updating the positions of followers.
<disp-formula id="eqn-20"><label>(17)</label><mml:math id="mml-eqn-20" display="block">
<mml:mrow><mml:msubsup><mml:mi>x</mml:mi><mml:mi>j</mml:mi><mml:mi>i</mml:mi></mml:msubsup><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac><mml:mo stretchy="false">(</mml:mo><mml:msubsup><mml:mi>x</mml:mi><mml:mi>j</mml:mi><mml:mi>i</mml:mi></mml:msubsup><mml:mo>+</mml:mo><mml:msubsup><mml:mi>x</mml:mi><mml:mi>j</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup><mml:mo stretchy="false">)</mml:mo><mml:mo>&#x2200;</mml:mo><mml:mspace width="thickmathspace" /><mml:mi>i</mml:mi><mml:mo>&#x2265;</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:math></disp-formula></p>
<p>While the global optimal solution of any optimization problem is unknown, the best solution can be obtained by moving the leader, followed by the followers, towards the food source. As a result, the salp chain moves towards the global optimum. The overall steps of the SSA are described below.
</p>
<fig id="fig-13">
<graphic mimetype="image" mime-subtype="png" xlink:href="CMC_25741-fig-13.png"/>
</fig>
</sec>
<sec id="s3_2"><label>3.2</label><title><italic>Hybrid Chaotic Salp Swarm with Crossover (HCSSC) Algorithm for</italic> <inline-formula id="ieqn-104"><mml:math id="mml-ieqn-104"><mml:mspace width="thickmathspace" /><mml:mi mathvariant="bold-italic">R</mml:mi><mml:mi mathvariant="bold-italic">E</mml:mi></mml:math></inline-formula></title>
<p>This subsection discusses the main steps of the proposed HCSSC algorithm to find the optimal source and relay nodes powers for maximizing the resource efficiency <inline-formula id="ieqn-105"><mml:math id="mml-ieqn-105"><mml:mi>R</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> at each packet transmission <italic>t</italic> in the UWSN. Assuming that <italic>B</italic> has perfect Channel State Information (CSI) awareness, the proposed HCSSC algorithm is implemented at <italic>B</italic> which will send the obtained powers values to source and relay nodes prior to their packets&#x2019; transmissions. The hybrid algorithm uses a logistic chaotic map to generate a feasible initial population at random. Furthermore, improving the solution during the number of iterations through the uniform crossover operator and chaotic map can guarantee that the optimal solution is the final solution. The steps of the proposed HCSSC are discussed in detail as follows.</p>
<sec id="s3_2_1"><label>3.2.1</label><title>Initial Population</title>
<p>Each salp (i.e.,<inline-formula id="ieqn-106"><mml:math id="mml-ieqn-106"><mml:mrow><mml:msup><mml:mi>x</mml:mi><mml:mi>i</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula>, <inline-formula id="ieqn-107"><mml:math id="mml-ieqn-107"><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mo fence="false" stretchy="false">{</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x2026;</mml:mo><mml:mo>,</mml:mo><mml:mi>n</mml:mi><mml:mo fence="false" stretchy="false">}</mml:mo></mml:math></inline-formula>) in the HCSSC population <italic>N</italic> corresponds to a possible solution for the resource efficiency optimization problem. Moreover, each salp consists of many variables (i.e.,<inline-formula id="ieqn-108"><mml:math id="mml-ieqn-108"><mml:mrow><mml:msub><mml:mrow><mml:mi>var</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>j</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>, <inline-formula id="ieqn-109"><mml:math id="mml-ieqn-109"><mml:mi>j</mml:mi><mml:mo>=</mml:mo><mml:mo fence="false" stretchy="false">{</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x2026;</mml:mo><mml:mo>,</mml:mo><mml:mi>m</mml:mi><mml:mo fence="false" stretchy="false">}</mml:mo></mml:math></inline-formula>) that affect the optimization of the resource efficiency as shown in <xref ref-type="fig" rid="fig-2">Fig. 2</xref>. This paper considers a salp <italic>x</italic> consists of two variables which are <inline-formula id="ieqn-110"><mml:math id="mml-ieqn-110"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi>S</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-111"><mml:math id="mml-ieqn-111"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mi>R</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>. An illustrative example of a salp <italic>x</italic> is shown in <xref ref-type="fig" rid="fig-3">Fig. 3</xref>. The value for each variable can be generated within the lower <inline-formula id="ieqn-112"><mml:math id="mml-ieqn-112"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">min</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> and upper <inline-formula id="ieqn-113"><mml:math id="mml-ieqn-113"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> bounds, respectively.</p>
<fig id="fig-2"><label>Figure 2</label><caption><title>HCSSC population</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_25741-fig-2.png"/></fig>
<fig id="fig-3"><label>Figure 3</label><caption><title>An illustrative example of salp <italic>x</italic></title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_25741-fig-3.png"/></fig>
<p>The diversity of the initial population has a great impact on spreading effectively in the search space. Therefore, for generating an effective initial population, a chaotic map is used in the proposed algorithm. One of the simplest maps is the logistic map that appears in the nonlinear dynamics of a biological population that evidences the chaotic behavior [<xref ref-type="bibr" rid="ref-31">31</xref>] and it is represented mathematically by <xref ref-type="disp-formula" rid="eqn-21">Eq. (18)</xref>.
<disp-formula id="eqn-21"><label>(18)</label><mml:math id="mml-eqn-21" display="block"><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mi>r</mml:mi><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mn>1</mml:mn><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:math></disp-formula>where <inline-formula id="ieqn-114"><mml:math id="mml-ieqn-114"><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the chaotic value at each independent run <italic>k</italic> (i.e.,<inline-formula id="ieqn-115"><mml:math id="mml-ieqn-115"><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow><mml:mo>&#x2208;</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mn>1</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula>), <italic>r</italic> is the growth rate that controls the behavior of chaotic value at a certain time <inline-formula id="ieqn-116"><mml:math id="mml-ieqn-116"><mml:mo stretchy="false">(</mml:mo><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mn>4</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula>. To be more specific, the generated random value of resources at any given time (round) is a function of the growth rate parameter and the previous time step&#x0027;s resource&#x0027;s value. Consequently, the initial population <italic>N</italic> of the proposed HCSSC algorithm is generated according to the following pseudo code.
</p>
<fig id="fig-14">
<graphic mimetype="image" mime-subtype="png" xlink:href="CMC_25741-fig-14.png"/>
</fig>
<p>Each variable is bounded within lower and upper values <inline-formula id="ieqn-123"><mml:math id="mml-ieqn-123"><mml:mo stretchy="false">[</mml:mo><mml:mi>l</mml:mi><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mi>j</mml:mi></mml:msub></mml:mrow><mml:mo>,</mml:mo><mml:mi>u</mml:mi><mml:mrow><mml:msub><mml:mi>b</mml:mi><mml:mi>j</mml:mi></mml:msub></mml:mrow><mml:mo stretchy="false">]</mml:mo></mml:math></inline-formula>, <inline-formula id="ieqn-124"><mml:math id="mml-ieqn-124"><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the chaotic sequence that is generated by logistic chaotic map. The integration of chaotic maps with salp swarm algorithm leads to produce a proper distribution through the characteristic of random and stochasticity of chaos.</p>
</sec>
<sec id="s3_2_2"><label>3.2.2</label><title>Evaluation of Salps</title>
<p>The resource efficiency optimization can be modelled as an optimization problem with a maximized objective function. Therefore, each individual (i.e., salp) in the population is evaluated according to <xref ref-type="disp-formula" rid="eqn-15">Eq. (12)</xref> and the fittest one is assigned to <italic>F</italic>.</p>
</sec>
<sec id="s3_2_3"><label>3.2.3</label><title>Updating Salps</title>
<p>The first salp in the population is called the leader which is responsible for guiding the salp chain and is continuously updating its position towards the direction of the source food. The rest salps in the population are called the followers because they follow the leader in updating their positions. Leader&#x0027;s and followers&#x2019; updates are illustrated below.
<list list-type="bullet">
<list-item><p>Leader&#x0027;s update</p></list-item>
</list></p>
<p>Since the leader updates its position in a positive direction, we ignore the negative direction. Also, the random numbers <inline-formula id="ieqn-125"><mml:math id="mml-ieqn-125"><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-126"><mml:math id="mml-ieqn-126"><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> are replaced by a chaotic sequences <inline-formula id="ieqn-127"><mml:math id="mml-ieqn-127"><mml:msubsup><mml:mi>c</mml:mi><mml:mi>k</mml:mi><mml:mn>1</mml:mn></mml:msubsup></mml:math></inline-formula> and <inline-formula id="ieqn-128"><mml:math id="mml-ieqn-128"><mml:msubsup><mml:mi>c</mml:mi><mml:mi>k</mml:mi><mml:mn>2</mml:mn></mml:msubsup></mml:math></inline-formula> computed according to <xref ref-type="disp-formula" rid="eqn-21">Eq. (18)</xref> while <inline-formula id="ieqn-129"><mml:math id="mml-ieqn-129"><mml:mrow><mml:msub><mml:mi>c</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> is calculated according to <xref ref-type="disp-formula" rid="eqn-17">Eq. (14)</xref>. Thus, <xref ref-type="disp-formula" rid="eqn-16">Eq. (13)</xref> is modified and represented in <xref ref-type="disp-formula" rid="eqn-22">Eq. (19)</xref>
<disp-formula id="eqn-22"><label>(19)</label><mml:math id="mml-eqn-22" display="block"><mml:mrow><mml:msubsup><mml:mi>x</mml:mi><mml:mi>j</mml:mi><mml:mn>1</mml:mn></mml:msubsup><mml:mo>=</mml:mo><mml:msub><mml:mi>F</mml:mi><mml:mi>j</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>c</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mo stretchy="false">(</mml:mo><mml:mi>u</mml:mi><mml:msub><mml:mi>b</mml:mi><mml:mi>j</mml:mi></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mi>l</mml:mi><mml:msub><mml:mi>b</mml:mi><mml:mi>j</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:msubsup><mml:mi>c</mml:mi><mml:mi>k</mml:mi><mml:mn>1</mml:mn></mml:msubsup><mml:mo>+</mml:mo><mml:mi>l</mml:mi><mml:msub><mml:mi>b</mml:mi><mml:mi>j</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo><mml:msubsup><mml:mi>c</mml:mi><mml:mi>k</mml:mi><mml:mn>2</mml:mn></mml:msubsup><mml:mo>&#x2265;</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:math></disp-formula>
<list list-type="bullet">
<list-item><p>Follower&#x0027;s updates</p></list-item>
</list></p>
<p>Each follower <inline-formula id="ieqn-130"><mml:math id="mml-ieqn-130"><mml:msubsup><mml:mi>x</mml:mi><mml:mi>j</mml:mi><mml:mi>i</mml:mi></mml:msubsup></mml:math></inline-formula>; <inline-formula id="ieqn-131"><mml:math id="mml-ieqn-131"><mml:mi>i</mml:mi><mml:mo>&#x2265;</mml:mo><mml:mn>2</mml:mn></mml:math></inline-formula> in the population updates its position based on how far the current position from the best salp&#x0027;s position. To achieve a better follow to the leader, we mate the best individual with the current individual through the uniform crossover [<xref ref-type="bibr" rid="ref-21">21</xref>,<xref ref-type="bibr" rid="ref-32">32</xref>]. The uniform crossover proves its efficiency compared to other crossover operators [<xref ref-type="bibr" rid="ref-20">20</xref>]. An illustrative example of uniform crossover is shown in <xref ref-type="fig" rid="fig-4">Fig. 4</xref>.</p>
<fig id="fig-4"><label>Figure 4</label><caption><title>An illustrative example of uniform crossover</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_25741-fig-4.png"/></fig>
<p>In the process of crossover, two parents (<inline-formula id="ieqn-132"><mml:math id="mml-ieqn-132"><mml:mi>a</mml:mi></mml:math></inline-formula> and <inline-formula id="ieqn-133"><mml:math id="mml-ieqn-133"><mml:mi>b</mml:mi></mml:math></inline-formula>) are selected and a binary mask consists of 1/0 digits is generated randomly with the same length of the individual. The result of this operation is generating two offsprings (<inline-formula id="ieqn-134"><mml:math id="mml-ieqn-134"><mml:mrow><mml:mover><mml:mi>a</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-135"><mml:math id="mml-ieqn-135"><mml:mrow><mml:mover><mml:mi>b</mml:mi><mml:mo>&#x007E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula>) based on the corresponding digit in the binary mask. For the first offspring (<inline-formula id="ieqn-136"><mml:math id="mml-ieqn-136"><mml:mrow><mml:mover><mml:mi>a</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula>), if the digit in a mask is 1, then the digit it is taken from (<inline-formula id="ieqn-137"><mml:math id="mml-ieqn-137"><mml:mi>a</mml:mi></mml:math></inline-formula>) otherwise from (<inline-formula id="ieqn-138"><mml:math id="mml-ieqn-138"><mml:mi>b</mml:mi></mml:math></inline-formula>). Regarding the generation offspring (<inline-formula id="ieqn-139"><mml:math id="mml-ieqn-139"><mml:mrow><mml:mover><mml:mi>b</mml:mi><mml:mo>&#x007E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula>), the complementary of mask is used. The two obtained offsprings <inline-formula id="ieqn-140"><mml:math id="mml-ieqn-140"><mml:mrow><mml:mover><mml:mi>a</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-141"><mml:math id="mml-ieqn-141"><mml:mrow><mml:mover><mml:mi>b</mml:mi><mml:mo>&#x007E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula> are considered as two new salps in the population. Each offspring <inline-formula id="ieqn-142"><mml:math id="mml-ieqn-142"><mml:mrow><mml:mover><mml:mi>x</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mo>&#x2208;</mml:mo><mml:mo fence="false" stretchy="false">{</mml:mo><mml:mrow><mml:mrow><mml:mover><mml:mi>a</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mo>,</mml:mo><mml:mrow><mml:mover><mml:mi>b</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mo fence="false" stretchy="false">}</mml:mo></mml:math></inline-formula> is a vector of two variables <inline-formula id="ieqn-143"><mml:math id="mml-ieqn-143"><mml:mrow><mml:mover><mml:mi>x</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula>&#x003D;<inline-formula id="ieqn-144"><mml:math id="mml-ieqn-144"><mml:mo stretchy="false">[</mml:mo><mml:mrow><mml:msubsup><mml:mi>P</mml:mi><mml:mi>S</mml:mi><mml:mi>t</mml:mi></mml:msubsup><mml:mo>,</mml:mo><mml:msubsup><mml:mi>P</mml:mi><mml:mi>R</mml:mi><mml:mi>t</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy="false">]</mml:mo></mml:math></inline-formula>. Then, the objective function, which is the resource efficiency at each transmission, <italic>t</italic>, <inline-formula id="ieqn-145"><mml:math id="mml-ieqn-145"><mml:mi>R</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> of each offspring <inline-formula id="ieqn-146"><mml:math id="mml-ieqn-146"><mml:mrow><mml:mover><mml:mi>x</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mo>&#x2208;</mml:mo><mml:mo fence="false" stretchy="false">{</mml:mo><mml:mrow><mml:mrow><mml:mover><mml:mi>a</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mo>,</mml:mo><mml:mrow><mml:mover><mml:mi>b</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mo fence="false" stretchy="false">}</mml:mo></mml:math></inline-formula> denoted by <inline-formula id="ieqn-147"><mml:math id="mml-ieqn-147"><mml:mi>R</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:mover><mml:mi>x</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> is evaluated as follows
<disp-formula id="eqn-23"><label>(20)</label><mml:math id="mml-eqn-23" display="block"><mml:mi>R</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:mover><mml:mi>x</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mi>R</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msubsup><mml:mi>P</mml:mi><mml:mi>S</mml:mi><mml:mi>t</mml:mi></mml:msubsup><mml:mo>,</mml:mo><mml:msubsup><mml:mi>P</mml:mi><mml:mi>R</mml:mi><mml:mi>t</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mi>E</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msubsup><mml:mi>P</mml:mi><mml:mi>S</mml:mi><mml:mi>t</mml:mi></mml:msubsup><mml:mo>,</mml:mo><mml:msubsup><mml:mi>P</mml:mi><mml:mi>R</mml:mi><mml:mi>t</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo>+</mml:mo><mml:mrow><mml:mi>&#x03C9;</mml:mi></mml:mrow><mml:mi>S</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msubsup><mml:mi>P</mml:mi><mml:mi>S</mml:mi><mml:mi>t</mml:mi></mml:msubsup><mml:mo>,</mml:mo><mml:msubsup><mml:mi>P</mml:mi><mml:mi>R</mml:mi><mml:mi>t</mml:mi></mml:msubsup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:math></disp-formula></p>
<p>Further, <inline-formula id="ieqn-148"><mml:math id="mml-ieqn-148"><mml:mi>R</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:mover><mml:mi>a</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> and <inline-formula id="ieqn-149"><mml:math id="mml-ieqn-149"><mml:mi>R</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:mover><mml:mi>b</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> are compared. If <inline-formula id="ieqn-150"><mml:math id="mml-ieqn-150"><mml:mi>R</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:mover><mml:mi>a</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mo>&#x003E;</mml:mo><mml:mi>R</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:mover><mml:mi>b</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula>, then, <inline-formula id="ieqn-151"><mml:math id="mml-ieqn-151"><mml:mrow><mml:mover><mml:mi>a</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow><mml:mspace width="thickmathspace" /></mml:math></inline-formula> is chosen to be <inline-formula id="ieqn-152"><mml:math id="mml-ieqn-152"><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>c</mml:mi><mml:mi>r</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>, else, <inline-formula id="ieqn-153"><mml:math id="mml-ieqn-153"><mml:mrow><mml:mover><mml:mi>b</mml:mi><mml:mo stretchy="false">&#x007E;</mml:mo></mml:mover></mml:mrow></mml:math></inline-formula> is chosen to be <inline-formula id="ieqn-154"><mml:math id="mml-ieqn-154"><mml:mspace width="thickmathspace" /><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>c</mml:mi><mml:mi>r</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>. Therefore, each follower in the salp chain updates its position as follows
<disp-formula id="eqn-24"><label>(21)</label><mml:math id="mml-eqn-24" display="block"><mml:msubsup><mml:mi>x</mml:mi><mml:mi>j</mml:mi><mml:mi>i</mml:mi></mml:msubsup><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mrow><mml:mi>u</mml:mi><mml:mi>c</mml:mi><mml:mi>r</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>+</mml:mo><mml:msubsup><mml:mi>x</mml:mi><mml:mi>j</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msubsup></mml:mrow><mml:mo stretchy="false">)</mml:mo><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mspace width="thickmathspace" /><mml:mi mathvariant="normal">&#x2200;</mml:mi><mml:mspace width="thickmathspace" /><mml:mi>i</mml:mi><mml:mo>&#x2265;</mml:mo><mml:mn>2</mml:mn></mml:math></disp-formula></p>
</sec>
<sec id="s3_2_4"><label>3.2.4</label><title>Test the Termination Condition</title>
<p>The leader and followers are updating their positions iteratively until reaching a maximum number of iterations. Once the HCSSC algorithm reached the termination condition, the global best salp is returned as the best solution so far for the resource efficiency <inline-formula id="ieqn-155"><mml:math id="mml-ieqn-155"><mml:mi>R</mml:mi><mml:mrow><mml:msup><mml:mi>E</mml:mi><mml:mi>t</mml:mi></mml:msup></mml:mrow></mml:math></inline-formula> at each packet transmission <italic>t</italic>. The overall steps of the HCSSC algorithm for optimizing the resource efficiency in UWSN are described in the <xref ref-type="fig" rid="fig-5">Fig. 5</xref>.</p>
<fig id="fig-5"><label>Figure 5</label><caption><title>Flowchart of the proposed HCSSC algorithm</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_25741-fig-5.png"/></fig>
</sec>
</sec>
</sec>
<sec id="s4"><label>4</label><title>Experimental Results and Analysis</title>
<p>This section presents the numerical results illustrating the performance of the proposed HCSSC based power optimization scheme. Simulations are done using MATLAB-R2015a running on Windows 7 with 2 GB RAM memory. Simulation results are obtained by averaging over 1000 channel iterations. We assume that source and relay nodes have equal batteries power capacities <inline-formula id="ieqn-156"><mml:math id="mml-ieqn-156"><mml:mi>C</mml:mi></mml:math></inline-formula>(i.e., <inline-formula id="ieqn-157"><mml:math id="mml-ieqn-157"><mml:mi>C</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mi>X</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) and equal maximum allowed powers <inline-formula id="ieqn-158"><mml:math id="mml-ieqn-158"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>. Authors in [<xref ref-type="bibr" rid="ref-13">13</xref>] determine the mathematical computation of UG2UG and UG2AG path losses. <xref ref-type="table" rid="table-1">Tab. 1</xref> summarizes the set of the system parameters and their corresponding values.</p>
<table-wrap id="table-1"><label>Table 1</label><caption><title>Simulation parameters</title></caption><table>
<colgroup>
<col align="left"/>
<col align="left"/>
</colgroup>
<thead>
<tr>
<th align="left">Parameter</th>
<th align="left">Value</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left"><inline-formula id="ieqn-227"><mml:math id="mml-ieqn-227"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">min</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> for <inline-formula id="ieqn-228"><mml:math id="mml-ieqn-228"><mml:mrow><mml:mtext>&#x00A0;&#xA0;</mml:mtext></mml:mrow><mml:mi>X</mml:mi><mml:mrow><mml:mtext>&#x00A0;&#xA0;</mml:mtext></mml:mrow><mml:mspace width="thinmathspace" /><mml:mo>&#x2208;</mml:mo><mml:mspace width="thinmathspace" /><mml:mrow><mml:mtext>&#x00A0;&#xA0;</mml:mtext></mml:mrow><mml:mo fence="false" stretchy="false">{</mml:mo><mml:mrow><mml:mi>S</mml:mi><mml:mo>,</mml:mo><mml:mrow><mml:mtext>&#x00A0;&#xA0;</mml:mtext></mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mo fence="false" stretchy="false">}</mml:mo></mml:math></inline-formula></td>
<td align="left">5&#x2005;mw</td>
</tr>
<tr>
<td align="left"><inline-formula id="ieqn-229"><mml:math id="mml-ieqn-229"><mml:mrow><mml:msub><mml:mi>x</mml:mi><mml:mi>X</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> for <inline-formula id="ieqn-230"><mml:math id="mml-ieqn-230"><mml:mrow><mml:mtext>&#x00A0;&#xA0;</mml:mtext></mml:mrow><mml:mi>X</mml:mi><mml:mspace width="thinmathspace" /><mml:mrow><mml:mtext>&#x00A0;&#xA0;</mml:mtext></mml:mrow><mml:mo>&#x2208;</mml:mo><mml:mspace width="thinmathspace" /><mml:mrow><mml:mtext>&#x00A0;&#xA0;</mml:mtext></mml:mrow><mml:mo fence="false" stretchy="false">{</mml:mo><mml:mrow><mml:mi>S</mml:mi><mml:mo>,</mml:mo><mml:mrow><mml:mtext>&#x00A0;&#xA0;</mml:mtext></mml:mrow><mml:mi>R</mml:mi></mml:mrow><mml:mo fence="false" stretchy="false">}</mml:mo></mml:math></inline-formula></td>
<td align="left">0.2&#x2005;m</td>
</tr>
<tr>
<td align="left"><inline-formula id="ieqn-231"><mml:math id="mml-ieqn-231"><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mi>S</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula></td>
<td align="left">0.6&#x2005;m</td>
</tr>
<tr>
<td align="left"><inline-formula id="ieqn-232"><mml:math id="mml-ieqn-232"><mml:mrow><mml:msub><mml:mi>y</mml:mi><mml:mi>R</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula></td>
<td align="left">0.1&#x2005;m</td>
</tr>
<tr>
<td align="left"><inline-formula id="ieqn-233"><mml:math id="mml-ieqn-233"><mml:mrow><mml:msub><mml:mi>h</mml:mi><mml:mi>B</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula></td>
<td align="left">0.7&#x2005;m</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The improvement of the resource efficiency using the proposed HCSSC algorithm, <inline-formula id="ieqn-159"><mml:math id="mml-ieqn-159"><mml:mi>R</mml:mi><mml:mrow><mml:msub><mml:mi>E</mml:mi><mml:mrow><mml:mi>H</mml:mi><mml:mi>C</mml:mi><mml:mi>S</mml:mi><mml:mi>S</mml:mi><mml:mi>C</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>, against the resource efficiency, which uses the standard SSA, <inline-formula id="ieqn-160"><mml:math id="mml-ieqn-160"><mml:mi>R</mml:mi><mml:mrow><mml:msub><mml:mi>E</mml:mi><mml:mrow><mml:mi>S</mml:mi><mml:mi>S</mml:mi><mml:mi>A</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>, is computed as follows:
<disp-formula id="eqn-25"><label>(22)</label><mml:math id="mml-eqn-25" display="block"><mml:mrow><mml:mi>i</mml:mi><mml:mi>m</mml:mi><mml:mi>p</mml:mi><mml:mi>o</mml:mi><mml:mi>v</mml:mi><mml:mi>e</mml:mi><mml:mi>m</mml:mi><mml:mi>e</mml:mi><mml:mi>n</mml:mi><mml:mi>t</mml:mi><mml:mtext>&#x00A0;</mml:mtext><mml:mi>p</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mi>f</mml:mi><mml:mi>o</mml:mi><mml:mi>r</mml:mi><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>n</mml:mi><mml:mi>c</mml:mi><mml:mi>e</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x0025;</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>R</mml:mi><mml:msub><mml:mi>E</mml:mi><mml:mrow><mml:mtext>HCSSC</mml:mtext></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mi>R</mml:mi><mml:msub><mml:mi>E</mml:mi><mml:mrow><mml:mi>S</mml:mi><mml:mi>S</mml:mi><mml:mi>A</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>R</mml:mi><mml:msub><mml:mi>E</mml:mi><mml:mrow><mml:mi>S</mml:mi><mml:mi>S</mml:mi><mml:mi>A</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac><mml:mo>&#x00D7;</mml:mo><mml:mn>100</mml:mn></mml:mrow></mml:math></disp-formula></p>
<p>In order to demonstrate the effectiveness of chaos theory for generating the initial population, a set of experiments have been conducted with various number of individuals, as shown in <xref ref-type="fig" rid="fig-5">Fig. 5</xref>, for the standard SSA ( i.e., with random initial salp positions) and SSA using chaotic map for initial positions generation with power capacity <inline-formula id="ieqn-161"><mml:math id="mml-ieqn-161"><mml:mi>C</mml:mi><mml:mo>=</mml:mo><mml:mn>3</mml:mn><mml:mi>w</mml:mi></mml:math></inline-formula> and maximum transmission power <inline-formula id="ieqn-162"><mml:math id="mml-ieqn-162"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>50</mml:mn><mml:mi>m</mml:mi><mml:mi>w</mml:mi></mml:math></inline-formula>.</p>
<p><xref ref-type="fig" rid="fig-6">Fig. 6</xref> illustrates the convergence of the resource efficiency <inline-formula id="ieqn-163"><mml:math id="mml-ieqn-163"><mml:mi>R</mml:mi><mml:mi>E</mml:mi></mml:math></inline-formula> obtained using both algorithms for different number of salp. Clearly, the generation of initial population using chaotic map has achieved better results in maximizing the resource efficiency compared to the standard SSA, during the number of iterations, due to the better exploration of the search space. The improvement is increased with the increase of <italic>N</italic>. Indeed, for <inline-formula id="ieqn-164"><mml:math id="mml-ieqn-164"><mml:mi>N</mml:mi><mml:mo>=</mml:mo><mml:mn>5</mml:mn></mml:math></inline-formula>, the SSA with chaotic initial positions gains 4.9&#x0025; better than standard SSA while at <inline-formula id="ieqn-165"><mml:math id="mml-ieqn-165"><mml:mi>N</mml:mi><mml:mo>=</mml:mo><mml:mn>40</mml:mn></mml:math></inline-formula>, it gains 7&#x0025; better.</p>
<fig id="fig-6"><label>Figure 6</label><caption><title>Convergence curves of the standard SSA and SSA using chaotic initial positions</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_25741-fig-6.png"/></fig>
<p>To demonstrate the efficiency of using the chaos theory and the crossover operator to achieve a maximum resource efficiency, the proposed HCSSC algorithm has been tested for different number of salps <inline-formula id="ieqn-166"><mml:math id="mml-ieqn-166"><mml:mi>N</mml:mi><mml:mo>&#x2208;</mml:mo><mml:mo fence="false" stretchy="false">{</mml:mo><mml:mn>5</mml:mn><mml:mo>,</mml:mo><mml:mn>10</mml:mn><mml:mo>,</mml:mo><mml:mn>20</mml:mn><mml:mo>,</mml:mo><mml:mn>40</mml:mn><mml:mo fence="false" stretchy="false">}</mml:mo></mml:math></inline-formula> at <inline-formula id="ieqn-167"><mml:math id="mml-ieqn-167"><mml:mi>C</mml:mi><mml:mo>=</mml:mo><mml:mn>3</mml:mn><mml:mi>w</mml:mi></mml:math></inline-formula> and <inline-formula id="ieqn-168"><mml:math id="mml-ieqn-168"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>50</mml:mn><mml:mi>m</mml:mi><mml:mi>w</mml:mi></mml:math></inline-formula>, as shown in <xref ref-type="fig" rid="fig-7">Fig. 7</xref>. The figure depicts the convergence of resource efficiency as a function of iterations numbers obtained using the proposed HCSSC algorithm in the nodes powers.</p>
<fig id="fig-7"><label>Figure 7</label><caption><title>Convergence curve of the resource efficiency using the HCSSC algorithm</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_25741-fig-7.png"/></fig>
<p>From the figure, the convergence of the proposed algorithm is rapidly obtained for various number of salps to reach the optimal <inline-formula id="ieqn-169"><mml:math id="mml-ieqn-169"><mml:mi>R</mml:mi><mml:mi>E</mml:mi></mml:math></inline-formula> values. Also, the increase in the number of salps enhances the computation accuracy of the optimal <inline-formula id="ieqn-170"><mml:math id="mml-ieqn-170"><mml:mi>R</mml:mi><mml:mi>E</mml:mi></mml:math></inline-formula>. Therefore, the proposed algorithm can be efficiently implemented with acceptable cost of computational complexity.</p>
<p>For a fair comparison, the individuals number involved in the swarm population and the maximum number of iterations are equal for both standard SSA and the proposed HCSSC algorithms. <xref ref-type="table" rid="table-2">Tab. 2</xref> illustrates the maximum (Max.), minimum (Min.), Average (Avg.) and Standard deviation (Std.) of <inline-formula id="ieqn-171"><mml:math id="mml-ieqn-171"><mml:mi>R</mml:mi><mml:mi>E</mml:mi></mml:math></inline-formula> at various <italic>N</italic> and various <inline-formula id="ieqn-172"><mml:math id="mml-ieqn-172"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> at <inline-formula id="ieqn-173"><mml:math id="mml-ieqn-173"><mml:mi>C</mml:mi><mml:mo>=</mml:mo><mml:mn>3</mml:mn><mml:mi>w</mml:mi></mml:math></inline-formula>. The best results are shown in bold.</p>
<table-wrap id="table-2"><label>Table 2</label><caption><title>Statistical results of <inline-formula id="ieqn-234"><mml:math id="mml-ieqn-234"><mml:mi>R</mml:mi><mml:mi>E</mml:mi></mml:math></inline-formula> (Mbits/Joule) for <inline-formula id="ieqn-235"><mml:math id="mml-ieqn-235"><mml:mi>C</mml:mi><mml:mspace width="thinmathspace" /><mml:mo>=</mml:mo><mml:mspace width="thinmathspace" /><mml:mn>3</mml:mn><mml:mi>w</mml:mi></mml:math></inline-formula></title></caption><table>
<colgroup>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
</colgroup>
<thead>
<tr>
<th align="left">N</th>
<th>P<sub>max</sub>(mw)</th>
<th align="center" colspan="4">SSA</th>
<th align="center" colspan="4">HCSSC</th>
</tr>
<tr>
<th/>
<th/>
<th align="left">Avg.</th>
<th align="left">Max.</th>
<th align="left">Min.</th>
<th align="left">Std.</th>
<th align="left">Avg.</th>
<th align="left">Max.</th>
<th align="left">Min.</th>
<th align="left">Std.</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left">5</td>
<td align="left">50</td>
<td align="left">561.98</td>
<td align="left">1297</td>
<td align="left">172.94</td>
<td align="left">138.54</td>
<td align="left">610.11</td>
<td align="left">1353</td>
<td align="left">276.21</td>
<td align="left">154.05</td>
</tr>
<tr>
<td/>
<td align="left">300</td>
<td align="left">146.8785</td>
<td align="left">638.6622</td>
<td align="left">3.0190</td>
<td align="left">80.4059</td>
<td align="left">166.576</td>
<td align="left">642.4605</td>
<td align="left">22.5500</td>
<td align="left">86.0926</td>
</tr>
<tr>
<td/>
<td align="left">600</td>
<td align="left">77.5271</td>
<td align="left">418.7184</td>
<td align="left">4.0183</td>
<td align="left">52.5403</td>
<td align="left">89.9588</td>
<td align="left">427.9407</td>
<td align="left">6.2260</td>
<td align="left">59.1111</td>
</tr>
<tr>
<td align="left">20</td>
<td align="left">50</td>
<td align="left">573.0342</td>
<td align="left">1361</td>
<td align="left">244.7552</td>
<td align="left">145.3766</td>
<td align="left">680.53</td>
<td align="left">1340</td>
<td align="left">238.48</td>
<td align="left">180.43</td>
</tr>
<tr>
<td/>
<td align="left">300</td>
<td align="left">151.4191</td>
<td align="left">658.4308</td>
<td align="left">19.0465</td>
<td align="left">86.4999</td>
<td align="left">168.89</td>
<td align="left">828.68</td>
<td align="left">26.56</td>
<td align="left">99.85</td>
</tr>
<tr>
<td/>
<td align="left">600</td>
<td align="left">78.2166</td>
<td align="left">482.6424</td>
<td align="left">4.4647</td>
<td align="left">54.0314</td>
<td align="left">101.43</td>
<td align="left">600</td>
<td align="left">6.063</td>
<td align="left">72.053</td>
</tr>
<tr>
<td align="left">40</td>
<td align="left">50</td>
<td align="left">592.1117</td>
<td align="left">1308</td>
<td align="left">279.2015</td>
<td align="left">159.5827</td>
<td align="left">681.45</td>
<td align="left">1413</td>
<td align="left">330.57</td>
<td align="left">176.21</td>
</tr>
<tr>
<td/>
<td align="left">300</td>
<td align="left">156.5011</td>
<td align="left">533.7262</td>
<td align="left">20.2041</td>
<td align="left">83.3139</td>
<td align="left">184.65</td>
<td align="left">809.02</td>
<td align="left">23.108</td>
<td align="left">96.072</td>
</tr>
<tr>
<td/>
<td align="left">600</td>
<td align="left">85.7410</td>
<td align="left">432.3858</td>
<td align="left">3.8117</td>
<td align="left">60.1572</td>
<td align="left">102.93</td>
<td align="left">586.96</td>
<td align="left">7.174</td>
<td align="left">7.3319</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>In comparison with the traditional SSA based optimization scheme, the proposed HCSSC outperforms it in all evaluated performance measurements. According to <xref ref-type="table" rid="table-2">Tab. 2</xref>, for <inline-formula id="ieqn-174"><mml:math id="mml-ieqn-174"><mml:mi>N</mml:mi><mml:mo>=</mml:mo><mml:mn>5</mml:mn></mml:math></inline-formula> and <inline-formula id="ieqn-175"><mml:math id="mml-ieqn-175"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>300</mml:mn><mml:mi>m</mml:mi><mml:mi>w</mml:mi></mml:math></inline-formula>, the average value of the <inline-formula id="ieqn-176"><mml:math id="mml-ieqn-176"><mml:mi>R</mml:mi><mml:mi>E</mml:mi></mml:math></inline-formula> based HCSSC reaches 166.576 (Mbits/Joule) and the average value of the <inline-formula id="ieqn-177"><mml:math id="mml-ieqn-177"><mml:mi>R</mml:mi><mml:mi>E</mml:mi></mml:math></inline-formula> based SSA reaches 146.8785 (Mbits/Joule). Hence, the gain is about 13&#x0025;. While for <inline-formula id="ieqn-178"><mml:math id="mml-ieqn-178"><mml:mi>N</mml:mi><mml:mo>=</mml:mo><mml:mn>40</mml:mn></mml:math></inline-formula> and <inline-formula id="ieqn-179"><mml:math id="mml-ieqn-179"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>50</mml:mn><mml:mi>m</mml:mi><mml:mi>w</mml:mi></mml:math></inline-formula>, the average value of <inline-formula id="ieqn-180"><mml:math id="mml-ieqn-180"><mml:mi>R</mml:mi><mml:mi>E</mml:mi></mml:math></inline-formula> based HCSSC is 681.45 (Mbits/Joule) while the average value of the <inline-formula id="ieqn-181"><mml:math id="mml-ieqn-181"><mml:mi>R</mml:mi><mml:mi>E</mml:mi></mml:math></inline-formula> based SSA equals 592.1117 (Mbits/Joule). So, the gain reaches 15&#x0025;.</p>
<p><xref ref-type="fig" rid="fig-8">Fig. 8</xref> shows the convergence behavior of the resource efficiency using the proposed algorithm HCSSC against the standard SSA, in nodes power optimization, with different number of salps <inline-formula id="ieqn-182"><mml:math id="mml-ieqn-182"><mml:mi>N</mml:mi><mml:mo>&#x2208;</mml:mo><mml:mo fence="false" stretchy="false">{</mml:mo><mml:mn>5</mml:mn><mml:mo>,</mml:mo><mml:mn>20</mml:mn><mml:mo>,</mml:mo><mml:mn>40</mml:mn><mml:mo fence="false" stretchy="false">}</mml:mo></mml:math></inline-formula> at <inline-formula id="ieqn-183"><mml:math id="mml-ieqn-183"><mml:mi>C</mml:mi><mml:mo>=</mml:mo><mml:mn>3</mml:mn><mml:mi>w</mml:mi></mml:math></inline-formula> and <inline-formula id="ieqn-184"><mml:math id="mml-ieqn-184"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>50</mml:mn><mml:mi>m</mml:mi><mml:mi>w</mml:mi></mml:math></inline-formula>.</p>
<fig id="fig-8"><label>Figure 8</label><caption><title>Convergence curve of HCSSC against standard SSA</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_25741-fig-8.png"/></fig>
<p>With different number of salps, the proposed HCSSC obtains a better <inline-formula id="ieqn-185"><mml:math id="mml-ieqn-185"><mml:mi>R</mml:mi><mml:mi>E</mml:mi></mml:math></inline-formula> convergence than standard SSA for different number of iterations. The amelioration is raising with the increase of <italic>N</italic>. At <inline-formula id="ieqn-186"><mml:math id="mml-ieqn-186"><mml:mi>N</mml:mi><mml:mo>=</mml:mo><mml:mn>5</mml:mn><mml:mo>,</mml:mo><mml:mn>20</mml:mn></mml:math></inline-formula> and <inline-formula id="ieqn-187"><mml:math id="mml-ieqn-187"><mml:mn>40</mml:mn></mml:math></inline-formula> the <inline-formula id="ieqn-188"><mml:math id="mml-ieqn-188"><mml:mi>R</mml:mi><mml:mi>E</mml:mi></mml:math></inline-formula> of HCSSC is raising up to 9&#x0025;, 15&#x0025; and 16&#x0025;, respectively.</p>
<p>Furthermore, <xref ref-type="fig" rid="fig-9">Fig. 9</xref> illustrates the effect of the maximum power permitted for a single transmission <inline-formula id="ieqn-189"><mml:math id="mml-ieqn-189"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> on resource efficiency for both HCSSC and SSA with respect to different power capacity <inline-formula id="ieqn-190"><mml:math id="mml-ieqn-190"><mml:mi>C</mml:mi><mml:mo>&#x2208;</mml:mo><mml:mo fence="false" stretchy="false">{</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mn>3</mml:mn><mml:mo>,</mml:mo><mml:mn>4</mml:mn><mml:mo fence="false" stretchy="false">}</mml:mo><mml:mi>w</mml:mi></mml:math></inline-formula> and <inline-formula id="ieqn-191"><mml:math id="mml-ieqn-191"><mml:mi>N</mml:mi><mml:mo>=</mml:mo><mml:mn>20</mml:mn></mml:math></inline-formula>.</p>
<fig id="fig-9"><label>Figure 9</label><caption><title>Resource efficiency of HCSSC and SSA <italic>vs.</italic> the maximum allowed power <inline-formula id="ieqn-225"><mml:math id="mml-ieqn-225"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula></title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_25741-fig-9.png"/></fig>
<p>Once more, compared to the power optimized SSA scheme, the HCSSC achieves a higher resource efficiency at the same power cost. Indeed, the combination of chaotic map and the cross over operations in the proposed power algorithm improves the search of optimal nodes powers considering the power physical limitations. Hence, the resource efficiency <inline-formula id="ieqn-192"><mml:math id="mml-ieqn-192"><mml:mi>R</mml:mi><mml:mi>E</mml:mi></mml:math></inline-formula> increases with the increase of nodes batteries capacity. Notably, the gain gap between HCSSC and SSA schemes is higher as <italic>C</italic> increases. In fact, for <inline-formula id="ieqn-193"><mml:math id="mml-ieqn-193"><mml:mi>C</mml:mi><mml:mo>=</mml:mo><mml:mn>2</mml:mn><mml:mi>w</mml:mi></mml:math></inline-formula>, the <inline-formula id="ieqn-194"><mml:math id="mml-ieqn-194"><mml:mi>R</mml:mi><mml:mi>E</mml:mi></mml:math></inline-formula> improvement reaches 15.3&#x0025; at <inline-formula id="ieqn-195"><mml:math id="mml-ieqn-195"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>50</mml:mn><mml:mi>m</mml:mi><mml:mi>w</mml:mi></mml:math></inline-formula> and reaches 19&#x0025; at <inline-formula id="ieqn-196"><mml:math id="mml-ieqn-196"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>400</mml:mn><mml:mi>m</mml:mi><mml:mi>w</mml:mi></mml:math></inline-formula>. While for <inline-formula id="ieqn-197"><mml:math id="mml-ieqn-197"><mml:mi>C</mml:mi><mml:mo>=</mml:mo><mml:mn>4</mml:mn><mml:mi>w</mml:mi></mml:math></inline-formula>, the gain in <inline-formula id="ieqn-198"><mml:math id="mml-ieqn-198"><mml:mi>R</mml:mi><mml:mi>E</mml:mi></mml:math></inline-formula> is 15.5&#x0025; at <inline-formula id="ieqn-199"><mml:math id="mml-ieqn-199"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>50</mml:mn><mml:mi>m</mml:mi><mml:mi>w</mml:mi></mml:math></inline-formula> and is 23.6&#x0025; at <inline-formula id="ieqn-200"><mml:math id="mml-ieqn-200"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>400</mml:mn><mml:mi>m</mml:mi><mml:mi>w</mml:mi></mml:math></inline-formula>.Since the increase in the maximum allowed power <inline-formula id="ieqn-201"><mml:math id="mml-ieqn-201"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> degrades the energy efficiency <inline-formula id="ieqn-202"><mml:math id="mml-ieqn-202"><mml:mi>E</mml:mi><mml:mi>E</mml:mi></mml:math></inline-formula> and the weighted spectral efficiency <inline-formula id="ieqn-203"><mml:math id="mml-ieqn-203"><mml:mi>&#x03C9;</mml:mi><mml:mi>S</mml:mi><mml:mi>E</mml:mi></mml:math></inline-formula> as well, the resource efficiency diminishes when <inline-formula id="ieqn-204"><mml:math id="mml-ieqn-204"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> increases. In addition, network designers can regulate <inline-formula id="ieqn-205"><mml:math id="mml-ieqn-205"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> depending on the accessible power resource of batteries capacities for a given resource efficiency specification.</p>
<p>The average consumed relay power and the average remaining relay battery per transmission for HCSSC and SSA <italic>vs.</italic><inline-formula id="ieqn-206"><mml:math id="mml-ieqn-206"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> with <inline-formula id="ieqn-207"><mml:math id="mml-ieqn-207"><mml:mi>N</mml:mi><mml:mo>=</mml:mo><mml:mn>20</mml:mn></mml:math></inline-formula> are shown in <xref ref-type="fig" rid="fig-10">Figs. 10</xref> and <xref ref-type="fig" rid="fig-11">11</xref>, respectively. In <xref ref-type="fig" rid="fig-10">Fig. 10</xref>, The average consumed relay power is shown for <inline-formula id="ieqn-208"><mml:math id="mml-ieqn-208"><mml:mi>C</mml:mi><mml:mo>=</mml:mo><mml:mn>3</mml:mn><mml:mi>w</mml:mi></mml:math></inline-formula>.</p>
<fig id="fig-10"><label>Figure 10</label><caption><title>Average consumed relay energy per transmission (HCSSC and SSA)</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_25741-fig-10.png"/></fig>
<fig id="fig-11"><label>Figure 11</label><caption><title>Average remaining relay battery per transmission (HCSSC and SSA)</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_25741-fig-11.png"/></fig>
<p>Clearly, the deployment of HCSSC in nodes powers optimization allows a better relay power conservation since the average consumed relay power per transmission is minimized compared with the standard SSA along with <inline-formula id="ieqn-209"><mml:math id="mml-ieqn-209"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> values. In <xref ref-type="fig" rid="fig-11">Fig. 11</xref>, the effect of the HCSSC algorithm on the average remaining relay battery per transmission compared to the SSA is illustrated for <inline-formula id="ieqn-210"><mml:math id="mml-ieqn-210"><mml:mi>C</mml:mi><mml:mo>&#x2208;</mml:mo><mml:mo fence="false" stretchy="false">{</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mn>3</mml:mn><mml:mo>,</mml:mo><mml:mn>4</mml:mn><mml:mo fence="false" stretchy="false">}</mml:mo><mml:mi>w</mml:mi></mml:math></inline-formula>. The proposed HCSSC based scheme significantly ameliorates the average remaining relay battery per transmission for all battery capacity value <italic>C</italic>. At <inline-formula id="ieqn-211"><mml:math id="mml-ieqn-211"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>350</mml:mn><mml:mi>m</mml:mi><mml:mi>w</mml:mi></mml:math></inline-formula>, the obtained gain equals 23&#x0025; at <inline-formula id="ieqn-212"><mml:math id="mml-ieqn-212"><mml:mi>C</mml:mi><mml:mo>=</mml:mo><mml:mn>2</mml:mn><mml:mi>w</mml:mi></mml:math></inline-formula>, it equals 21.6&#x0025; at <inline-formula id="ieqn-213"><mml:math id="mml-ieqn-213"><mml:mi>C</mml:mi><mml:mo>=</mml:mo><mml:mn>3</mml:mn><mml:mi>w</mml:mi></mml:math></inline-formula> and it reaches 20.4&#x0025; at <inline-formula id="ieqn-214"><mml:math id="mml-ieqn-214"><mml:mi>C</mml:mi><mml:mo>=</mml:mo><mml:mn>4</mml:mn><mml:mi>w</mml:mi></mml:math></inline-formula>. Thus, the proposed HCSSC based power optimization extends the nodes batteries lifetime and consequently the whole network lifetime.</p>
<p>In <xref ref-type="fig" rid="fig-12">Fig. 12</xref>, we propose to study the resource efficiency performance in the case where the relay node has a higher battery capacity than the source node and can forward packets considering higher maximum allowed power. In fact, the relay node is expected to consume more power than the source node since it collects data from different sources to forward it to the sink node and, possibly, retransmits lost packets. Interestingly, the proposed algorithm proves its efficiency, not only in the case of equal source and relay batteries capacities, but also in the case of different source and relay batteries capacities as clearly shown in <xref ref-type="fig" rid="fig-12">Fig. 12</xref>. We assume that <inline-formula id="ieqn-215"><mml:math id="mml-ieqn-215"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>2</mml:mn><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-216"><mml:math id="mml-ieqn-216"><mml:mi>N</mml:mi><mml:mo>=</mml:mo><mml:mn>20</mml:mn></mml:math></inline-formula>. The proposed HCSSC based power optimization scheme achieves higher resource efficiency performance than the SSA as <inline-formula id="ieqn-217"><mml:math id="mml-ieqn-217"><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mi>S</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula id="ieqn-218"><mml:math id="mml-ieqn-218"><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mi>R</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> increase. At <inline-formula id="ieqn-219"><mml:math id="mml-ieqn-219"><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mi>S</mml:mi></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>2</mml:mn><mml:mi>w</mml:mi></mml:math></inline-formula> and <inline-formula id="ieqn-220"><mml:math id="mml-ieqn-220"><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mi>R</mml:mi></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>3</mml:mn><mml:mi>w</mml:mi></mml:math></inline-formula>, the HCSSC algorithm achieves a gain of 21&#x0025; better than the SSA at <inline-formula id="ieqn-221"><mml:math id="mml-ieqn-221"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>s</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>50</mml:mn><mml:mi>m</mml:mi><mml:mi>w</mml:mi></mml:math></inline-formula>. While at <inline-formula id="ieqn-222"><mml:math id="mml-ieqn-222"><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mi>S</mml:mi></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>3</mml:mn><mml:mi>w</mml:mi></mml:math></inline-formula> and <inline-formula id="ieqn-223"><mml:math id="mml-ieqn-223"><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mi>R</mml:mi></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>4</mml:mn><mml:mi>w</mml:mi></mml:math></inline-formula>, the HCSSC algorithm achieves a gain of 28&#x0025; better than the SSA at <inline-formula id="ieqn-224"><mml:math id="mml-ieqn-224"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mo>=</mml:mo><mml:mn>200</mml:mn><mml:mi>m</mml:mi><mml:mi>w</mml:mi></mml:math></inline-formula>.</p>
<fig id="fig-12"><label>Figure 12</label><caption><title>Resource efficiency of HCSSC and SSA for <inline-formula id="ieqn-226"><mml:math id="mml-ieqn-226"><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mspace width="thinmathspace" /><mml:mo>=</mml:mo><mml:mspace width="thinmathspace" /><mml:mn>2</mml:mn><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> for different source and relay batteries capacities</title></caption><graphic mimetype="image" mime-subtype="png" xlink:href="CMC_25741-fig-12.png"/></fig>
<p>In some applications, the benefit from of sensor nodes in UWSNs is restricteddue to harsh environmental conditions. For agricultural application for example, authors in [<xref ref-type="bibr" rid="ref-4">4</xref>] mention that if the water volume fraction of the mixture is high (passed 25&#x0025;), the UG2UG communication link is interrupted specially with some particular soil type whose the capacity to hold the bound water is low. Consequently, the UG2UG communication can be interrupted for a long period in case of a rainfall. Also, for underground mine applications, UG2UG communications may be interrupted in many unpredictable situations such as rock falls or explosions. Moreover, we notice that the base station should implement a powerful operating system to support the additional computing complexity of the HCSSC algorithm.</p>
</sec>
<sec id="s5"><label>5</label><title>Conclusions</title>
<p>This paper proposed a Hybrid Chaotic Salp Swarm with Crossover (HCSSC) algorithm for an UWSN to maximize the network resource efficiency. This last is a global metric that jointly considers the energy and the spectral efficiencies to balance the power consumption and the bandwidth usage. The algorithm improves the standard metaheuristic SSA by the use of logistic chaotic map in the generation of the initial population and the deployment of the uniform crossover operator to compute the final solution. At each packet transmission, the HCSSC is applied to provide the optimal source and relay nodes powers considering the remaining nodes batteries capacities constraints. Simulations showed that the integration of the chaotic map in the population initialization and the use of the crossover method in the positions&#x2019; updates improved the resource efficiency compared to the standard SSA for different nodes batteries capacities and different maximum allowed powers. Also, the use of the HCSSC algorithm offered a better relay power conservation proved by the minimization of the average consumed relay power per transmission and the maximization of the average remaining relay battery per transmission. Moreover, the efficiency of the proposed algorithm is demonstrated in the case where the relay node has a higher battery capacity than the source node and can forward packets considering higher maximum power. As future work, the efficiency of the proposed HCSSC algorithm in multi-relay UWSN, where many relay nodes cooperate with source nodes to transmit sensory data to the base station, will be addressed. Thus, the impact of the numbers of variables on HCSSC algorithm performance will be effectively studied. Moreover, the impact of the data packet size on the RE performance will be studied.</p>
</sec>
</body>
<back>
<fn-group>
<fn fn-type="other"><p><bold>Funding Statement:</bold> The authors received no specific funding for this study.</p></fn>
<fn fn-type="conflict"><p><bold>Conflicts of Interest:</bold> The authors declare that they have no conflicts of interest to report regarding the present study.</p></fn>
</fn-group>
<ref-list content-type="authoryear">
<title>References</title>
<ref id="ref-1"><label>[1]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>G.</given-names> <surname>Zhou</surname></string-name></person-group>, &#x201C;<article-title>Recent advances in wireless underground sensor network</article-title>,&#x201D; <source>International Journal of Distributed Sensor Networks</source>, vol. <volume>15</volume>, no. <issue>1</issue>, pp. <fpage>11</fpage>&#x2013;<lpage>21</lpage>, <year>2019</year>.</mixed-citation></ref>
<ref id="ref-2"><label>[2]</label><mixed-citation publication-type="conf-proc"><person-group person-group-type="author"><string-name><given-names>X.</given-names> <surname>Zhang</surname></string-name>, <string-name><given-names>A.</given-names> <surname>Andreyev</surname></string-name>, <string-name><given-names>C.</given-names> <surname>Zumpf</surname></string-name>, <string-name><given-names>M. C.</given-names> <surname>Negri</surname></string-name>, <string-name><given-names>S.</given-names> <surname>Guha</surname></string-name> <etal>et al.,</etal></person-group> &#x201C;<article-title>A Fully-buried wireless underground sensor network in an urban environment</article-title>,&#x201D; in <conf-name>Proc. COMSNETS</conf-name>, <conf-loc>Bengaluru, India</conf-loc>, pp. <fpage>239</fpage>&#x2013;<lpage>250</lpage>, <year>2019</year>.</mixed-citation></ref>
<ref id="ref-3"><label>[3]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>X.</given-names> <surname>Dong</surname></string-name>, <string-name><given-names>M. C.</given-names> <surname>Vuran</surname></string-name> and <string-name><given-names>S.</given-names> <surname>Irmak</surname></string-name></person-group>, &#x201C;<article-title>Autonomous precision agriculture through integration of wireless underground sensor networks with center pivot irrigation systems</article-title>,&#x201D; <source>Ad Hoc Networks</source>, vol. <volume>11</volume>, no. <issue>1</issue>, pp. <fpage>1975</fpage>&#x2013;<lpage>1987</lpage>, <year>2013</year>.</mixed-citation></ref>
<ref id="ref-4"><label>[4]</label><mixed-citation publication-type="conf-proc"><person-group person-group-type="author"><string-name><given-names>K.</given-names> <surname>Lin</surname></string-name>, <string-name><given-names>T.</given-names> <surname>Hao</surname></string-name>, <string-name><given-names>Z.</given-names> <surname>Yu</surname></string-name>, <string-name><given-names>W.</given-names> <surname>Zhen</surname></string-name> and <string-name><given-names>W.</given-names> <surname>He</surname></string-name></person-group>, &#x201C;<article-title>A preliminary study of ug2ag link quality in LoRa-based wireless underground sensor networks</article-title>,&#x201D; in <conf-name>Proc. LCN</conf-name>, <conf-loc>Osnabrueck, Germany</conf-loc>, pp. <fpage>51</fpage>&#x2013;<lpage>59</lpage>, <year>2019</year>.</mixed-citation></ref>
<ref id="ref-5"><label>[5]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>N. T.</given-names> <surname>Tam</surname></string-name>, <string-name><given-names>D. A.</given-names> <surname>Dung</surname></string-name>, <string-name><given-names>T. H.</given-names> <surname>Hung</surname></string-name>, <string-name><given-names>H. T. T.</given-names> <surname>Binh</surname></string-name> and <string-name><given-names>S.</given-names> <surname>Yu</surname></string-name></person-group>, &#x201C;<article-title>Exploiting relay nodes for maximizing wireless underground sensor network lifetime</article-title>,&#x201D; <source>Applied Intelligence</source>, vol. <volume>50</volume>, no. <issue>12</issue>, pp. <fpage>4568</fpage>&#x2013;<lpage>4585</lpage>, <year>2020</year>.</mixed-citation></ref>
<ref id="ref-6"><label>[6]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>D.</given-names> <surname>Wu</surname></string-name>, <string-name><given-names>D.</given-names> <surname>Chatzigeorgiou</surname></string-name>, <string-name><given-names>K.</given-names> <surname>Youcef-Toumi</surname></string-name>, <string-name><given-names>S.</given-names> <surname>Mekid</surname></string-name> and <string-name><given-names>R.</given-names> <surname>Ben-Mansour</surname></string-name></person-group>, &#x201C;<article-title>Channel-aware relay node placement in wireless sensor networks for pipeline inspection</article-title>,&#x201D; <source>IEEE Transactions on Wireless Communications</source>, vol. <volume>13</volume>, no. <issue>7</issue>, pp. <fpage>3510</fpage>&#x2013;<lpage>3523</lpage>, <year>2014</year>.</mixed-citation></ref>
<ref id="ref-7"><label>[7]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>B.</given-names> <surname>Yuan</surname></string-name>, <string-name><given-names>H.</given-names> <surname>Chen</surname></string-name> and <string-name><given-names>X.</given-names> <surname>Yao</surname></string-name></person-group>, &#x201C;<article-title>Optimal relay placement for lifetime maximization in wireless underground sensor networks</article-title>,&#x201D; <source>Information Sciences</source>, vol. <volume>418</volume>, no. <issue>2</issue>, pp. <fpage>463</fpage>&#x2013;<lpage>479</lpage>, <year>2017</year>.</mixed-citation></ref>
<ref id="ref-8"><label>[8]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>R.</given-names> <surname>Sharma</surname></string-name> and <string-name><given-names>S.</given-names> <surname>Prakash</surname></string-name></person-group>, &#x201C;<article-title>Enhancement of relay nodes communication approach in wsn-iot for underground coal mine</article-title>,&#x201D; <source>Journal of Information and Optimization Sciences</source>, vol. <volume>41</volume>, no. <issue>1</issue>, pp. <fpage>521</fpage>&#x2013;<lpage>531</lpage>, <year>2020</year>.</mixed-citation></ref>
<ref id="ref-9"><label>[9]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>J.</given-names> <surname>Tang</surname></string-name>, <string-name><given-names>D. K.</given-names> <surname>So</surname></string-name>, <string-name><given-names>E.</given-names> <surname>Alsusa</surname></string-name> and <string-name><given-names>K. A.</given-names> <surname>Hamdi</surname></string-name></person-group>, &#x201C;<article-title>Resource efficiency: A new paradigm on energy efficiency and spectral efficiency tradeoff</article-title>,&#x201D; <source>IEEE Transactions on Wireless Communications</source>, vol. <volume>13</volume>, no. <issue>8</issue>, pp. <fpage>4656</fpage>&#x2013;<lpage>4669</lpage>, <year>2014</year>.</mixed-citation></ref>
<ref id="ref-10"><label>[10]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>L.</given-names> <surname>Wei</surname></string-name>, <string-name><given-names>R. Q.</given-names>, <surname>Hu</surname></string-name>, <string-name><given-names>Y.</given-names> <surname>Qian</surname></string-name> and <string-name><given-names>G.</given-names> <surname>Wu</surname></string-name></person-group>, &#x201C;<article-title>Energy efficiency and spectrum efficiency of multihop device-to-device communications underlaying cellular networks</article-title>, &#x201C; <source>IEEE Transactions on Vehicular Technology</source>, vol. <volume>65</volume>, no. <issue>1</issue>, pp. <fpage>367</fpage>&#x2013;<lpage>380</lpage>, <year>2015</year>.</mixed-citation></ref>
<ref id="ref-11"><label>[11]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>J.</given-names> <surname>Tang</surname></string-name>, <string-name><given-names>D. K.</given-names> <surname>So</surname></string-name>, <string-name><given-names>E.</given-names> <surname>Alsusa</surname></string-name>, <string-name><given-names>K. A.</given-names> <surname>Hamdi</surname></string-name> and <string-name><given-names>A.</given-names> <surname>Shojaeifard</surname></string-name></person-group>, &#x201C;<article-title>On the energy efficiency&#x2013;spectral efficiency tradeoff in mimo-ofdma broadcast channels</article-title>,&#x201D; <source>IEEE Transactions on Vehicular Technology</source>, vol. <volume>65</volume>, no. <issue>7</issue>, pp. <fpage>5185</fpage>&#x2013;<lpage>5199</lpage>, <year>2015</year>.</mixed-citation></ref>
<ref id="ref-12"><label>[12]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>Y.</given-names> <surname>Wu</surname></string-name>, <string-name><given-names>Y.</given-names> <surname>Chen</surname></string-name>, <string-name><given-names>J.</given-names> <surname>Tang</surname></string-name>, <string-name><given-names>D. K.</given-names> <surname>So</surname></string-name>, <string-name><given-names>Z.</given-names> <surname>Xu</surname></string-name> <etal>et al.,</etal></person-group> &#x201C;<article-title>Green transmission technologies for balancing the energy efficiency and spectrum efficiency trade-off</article-title>,&#x201D; <source>IEEE Communications Magazine</source>, vol. <volume>52</volume>, no. <issue>11</issue>, pp. <fpage>112</fpage>&#x2013;<lpage>120</lpage>, <year>2014</year>.</mixed-citation></ref>
<ref id="ref-13"><label>[13]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>M.</given-names> <surname>Ayedi</surname></string-name>, <string-name><given-names>E.</given-names> <surname>Eldesouky</surname></string-name> and <string-name><given-names>J.</given-names> <surname>Nazeer</surname></string-name></person-group>, &#x201C;<article-title>Energy-spectral efficiency optimization in wireless underground sensor networks using salp swarm algorithm</article-title>,&#x201D; <source>Journal of Sensors</source>, vol. <volume>2021</volume>, no. <issue>1</issue>, pp. <fpage>1</fpage>&#x2013;<lpage>16</lpage>, <year>2021</year>.</mixed-citation></ref>
<ref id="ref-14"><label>[14]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>S.</given-names> <surname>Mirjalili</surname></string-name>, <string-name><given-names>A. H.</given-names> <surname>Gandomi</surname></string-name>, <string-name><given-names>S. Z.</given-names> <surname>Mirjalili</surname></string-name>, <string-name><given-names>S.</given-names> <surname>Saremi</surname></string-name>, <string-name><given-names>H.</given-names> <surname>Faris</surname></string-name> <etal>et al.,</etal></person-group> &#x201C;<article-title>Salp swarm algorithm: A bio-inspired optimizer for engineering design problems</article-title>,&#x201D; <source>Advances in Engineering Software</source>, vol. <volume>114</volume>, no. <issue>1</issue>, pp. <fpage>163</fpage>&#x2013;<lpage>191</lpage>, <year>2017</year>.</mixed-citation></ref>
<ref id="ref-15"><label>[15]</label><mixed-citation publication-type="book"><person-group person-group-type="author"><string-name><given-names>A.</given-names> <surname>Nayyar</surname></string-name>, <string-name><given-names>S.</given-names> <surname>Garg</surname></string-name>, <string-name><given-names>D.</given-names> <surname>Gupta</surname></string-name> and <string-name><given-names>A.</given-names> <surname>Khanna</surname></string-name></person-group>, in <source>Evolutionary Computation: Theory and Algorithms</source>, <publisher-loc>6000 Broken Sound Parkway NW, Suite 300 Boca Raton, FL 33487-2742</publisher-loc>: <publisher-name>Chapman and Hall/CRC Press</publisher-name>, <year>2018</year>. [Online]. Available: <uri xlink:href="https://cutt.ly/IUPfU5Y">https://cutt.ly/IUPfU5Y</uri>.</mixed-citation></ref>
<ref id="ref-16"><label>[16]</label><mixed-citation publication-type="book"><person-group person-group-type="author"><string-name><given-names>A.</given-names> <surname>Nayyar</surname></string-name>, <string-name><given-names>D. N.</given-names> <surname>Le</surname></string-name> and <string-name><given-names>N. G.</given-names> <surname>Nguyen</surname></string-name></person-group>, in <source>Advances in Swarm Intelligence for Optimizing Problems in Computer Science</source>, <publisher-loc>6000 Broken Sound Parkway NW, Suite 300 Boca Raton, FL 33487-2742</publisher-loc>: <publisher-name>Chapman and Hall/CRC Press</publisher-name>, <year>2018</year>. [Online]. Available: <uri xlink:href="https://cutt.ly/nUPft8p">https://cutt.ly/nUPft8p</uri>.</mixed-citation></ref>
<ref id="ref-17"><label>[17]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>H. M.</given-names> <surname>Kanoosh</surname></string-name>, <string-name><given-names>E. H.</given-names> <surname>Houssein</surname></string-name> and <string-name><given-names>M. M.</given-names> <surname>Selim</surname></string-name></person-group>, &#x201C;<article-title>Salp swarm algorithm for node localization in wireless sensor networks</article-title>,&#x201D; <source>Journal of Computer Networks and Communications</source>, vol. <volume>2019</volume>, no. <issue>1</issue>, pp. <fpage>1</fpage>&#x2013;<lpage>12</lpage>, <year>2019</year>.</mixed-citation></ref>
<ref id="ref-18"><label>[18]</label><mixed-citation publication-type="conf-proc"><person-group person-group-type="author"><string-name><given-names>X.</given-names> <surname>Shi</surname></string-name>, <string-name><given-names>J.</given-names> <surname>Su</surname></string-name>, <string-name><given-names>Z.</given-names> <surname>Ye</surname></string-name>, <string-name><given-names>F.</given-names> <surname>Chen</surname></string-name>, <string-name><given-names>P.</given-names> <surname>Zhang</surname></string-name> <etal>et al.,</etal></person-group> &#x201C;<article-title>A wireless sensor network node location method based on salp swarm algorithm</article-title>,&#x201D; in <conf-name>Proc. IDAACS</conf-name>, <conf-loc>Metz, France</conf-loc>, pp. <fpage>357</fpage>&#x2013;<lpage>361</lpage>, <year>2019</year>.</mixed-citation></ref>
<ref id="ref-19"><label>[19]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>M. A.</given-names> <surname>Syed</surname></string-name> and <string-name><given-names>R.</given-names> <surname>Syed</surname></string-name></person-group>, &#x201C;<article-title>Weighted salp swarm algorithm and its applications towards optimal sensor deployment</article-title>,&#x201D; <source>Journal of King Saud University-Computer and Information Sciences</source>, vol. <volume>2019</volume>, pp. <fpage>1</fpage>&#x2013;<lpage>11</lpage>, <year>2019</year>.</mixed-citation></ref>
<ref id="ref-20"><label>[20]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>S. E.</given-names> <surname>Jorgensen</surname></string-name> and <string-name><given-names>B. D.</given-names> <surname>Fath</surname></string-name></person-group>, &#x201C;<article-title>Encyclopedia of ecology</article-title>,&#x201D; <source>Chaos</source>, <volume>1 A-C</volume>, pp. <fpage>550</fpage>&#x2013;<lpage>551</lpage>, <year>2014</year>.</mixed-citation></ref>
<ref id="ref-21"><label>[21]</label><mixed-citation publication-type="conf-proc"><person-group person-group-type="author"><string-name><given-names>G.</given-names> <surname>Syswerda</surname></string-name></person-group>, &#x201C;<article-title>Uniform crossover in genetic algorithms</article-title>,&#x201D; in <conf-name>Proc. ICGA</conf-name>, <conf-loc>San Francisco, CA, United States</conf-loc>, pp. <fpage>2</fpage>&#x2013;<lpage>9</lpage>, <year>1989</year>.</mixed-citation></ref>
<ref id="ref-22"><label>[22]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>A.</given-names> <surname>Hussain</surname></string-name>, <string-name><given-names>Y. S.</given-names> <surname>Muhammad</surname></string-name> and <string-name><given-names>M. N.</given-names> <surname>Sajid</surname></string-name></person-group>, &#x201C;<article-title>An efficient genetic algorithm</article-title>,&#x201D; <source>International Journal of Mathematical Sciences and Computing</source>, vol. <volume>4</volume>, no. <issue>4</issue>, pp. <fpage>41</fpage>&#x2013;<lpage>55</lpage>, <year>2018</year>.</mixed-citation></ref>
<ref id="ref-23"><label>[23]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>M. C.</given-names> <surname>Vuran</surname></string-name> and <string-name><given-names>I. F.</given-names> <surname>Akyildiz</surname></string-name></person-group>, &#x201C;<article-title>Channel model and analysis for wireless underground sensor networks in soil medium</article-title>,&#x201D; <source>Physical Communication</source>, vol. <volume>3</volume>, no. <issue>4</issue>, pp. <fpage>245</fpage>&#x2013;<lpage>254</lpage>, <year>2010</year>.</mixed-citation></ref>
<ref id="ref-24"><label>[24]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>M. C.</given-names> <surname>Dobson</surname></string-name>, <string-name><given-names>F. T.</given-names> <surname>Ulaby</surname></string-name>, <string-name><given-names>M. T.</given-names> <surname>Hallikainen</surname></string-name> and <string-name><given-names>M. A.</given-names> <surname>El-Rayes</surname></string-name></person-group>, &#x201C;<article-title>Microwave dielectric behavior of wet soil-part ii: Dielectric mixing models</article-title>,&#x201D; <source>IEEE Transactions on Geoscience and Remote Sensing</source>, vol. <volume>1</volume>, no. <issue>1985</issue>, pp. <fpage>35</fpage>&#x2013;<lpage>46</lpage>, <year>1985</year>.</mixed-citation></ref>
<ref id="ref-25"><label>[25]</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name><given-names>W. E.</given-names> <surname>Patitz</surname></string-name>, <string-name><given-names>B. C.</given-names> <surname>Brock</surname></string-name> and <string-name><given-names>E. G.</given-names> <surname>Powell</surname></string-name></person-group>, &#x201C;<article-title>Measurement of dielectric and magnetic properties of soil</article-title>,&#x201D; <source>Sandia National Labs</source>, <year>1995</year>.</mixed-citation></ref>
<ref id="ref-26"><label>[26]</label><mixed-citation publication-type="conf-proc"><person-group person-group-type="author"><string-name><given-names>X.</given-names> <surname>Dong</surname></string-name> and <string-name><given-names>M. C.</given-names> <surname>Vuran</surname></string-name></person-group>, &#x201C;<article-title>Impacts of soil moisture on cognitive radio underground networks</article-title>,&#x201D; in <conf-name>Proc. BlackSeaCom</conf-name>, <conf-loc>Batumi, Georgia</conf-loc>, pp. <fpage>222</fpage>&#x2013;<lpage>227</lpage>, <year>2013</year>.</mixed-citation></ref>
<ref id="ref-27"><label>[27]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>L.</given-names> <surname>Abualigah</surname></string-name>, <string-name><given-names>M.</given-names> <surname>Shehab</surname></string-name>, <string-name><given-names>M.</given-names> <surname>Alshinwan</surname></string-name> and <string-name><given-names>H.</given-names> <surname>Alabool</surname></string-name></person-group>, &#x201C;<article-title>Salp swarm algorithm: A comprehensive survey</article-title>,&#x201D; <source>Neural Computing and Applications</source>, vol. <volume>32</volume>, no. <issue>15</issue>, pp. <fpage>11195</fpage>&#x2013;<lpage>11215</lpage>, <year>2019</year>.</mixed-citation></ref>
<ref id="ref-28"><label>[28]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>W. H.</given-names> <surname>El-Ashmawi</surname></string-name> and <string-name><given-names>A. F.</given-names> <surname>Ali</surname></string-name></person-group>, &#x201C;<article-title>A modified salp swarm algorithm for task assignment problem</article-title>,&#x201D; <source>Applied Soft Computing</source>, vol. <volume>94</volume>, no. <issue>2020</issue>, pp. <fpage>106445</fpage>, <year>2020</year>.</mixed-citation></ref>
<ref id="ref-29"><label>[29]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>L. P.</given-names> <surname>Madin</surname></string-name></person-group>, &#x201C;<article-title>Aspects of jet propulsion in salps</article-title>,&#x201D; <source>Canadian Journal of Zoology</source>, vol. <volume>68</volume>, no. <issue>4</issue>, pp. <fpage>765</fpage>&#x2013;<lpage>777</lpage>, <year>1990</year>.</mixed-citation></ref>
<ref id="ref-30"><label>[30]</label><mixed-citation publication-type="conf-proc"><person-group person-group-type="author"><string-name><given-names>P. A. V.</given-names> <surname>Anderson</surname></string-name> and <string-name><given-names>Q.</given-names> <surname>Bone</surname></string-name></person-group>, &#x201C;<article-title>Communication between individuals in salp chains. ii. physiology</article-title>,&#x201D; in <conf-name>Proc. of the Royal Society of London. Series B. Biological Sciences</conf-name>, <conf-loc>London</conf-loc>, vol. <volume>210</volume>, no. <issue>1181</issue>, pp. <fpage>559</fpage>&#x2013;<lpage>574</lpage>, <year>1980</year>.</mixed-citation></ref>
<ref id="ref-31"><label>[31]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>R.</given-names> <surname>Tang</surname></string-name>, <string-name><given-names>S.</given-names> <surname>Fong</surname></string-name> and <string-name><given-names>N.</given-names> <surname>Dey</surname></string-name></person-group>, &#x201C;<article-title>Metaheuristics and chaos theory</article-title>,&#x201D; <source>Chaos Theory, a</source>, vol. <volume>1</volume>, pp. <fpage>182</fpage>&#x2013;<lpage>196</lpage>, <year>2018</year>.</mixed-citation></ref>
<ref id="ref-32"><label>[32]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><given-names>P.</given-names> <surname>Kora</surname></string-name> and <string-name><given-names>P.</given-names> <surname>Yadlapalli</surname></string-name></person-group>, &#x201C;<article-title>Crossover operators in genetic algorithms: A review</article-title>,&#x201D; <source>International Journal of Computer Applications</source>, vol. <volume>6</volume>, no. <issue>1</issue>, pp. <fpage>1</fpage>&#x2013;<lpage>10</lpage>, <year>2017</year>.</mixed-citation></ref>
</ref-list>
</back>
</article>