<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Phys.</journal-id>
<journal-title>Frontiers in Physics</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Phys.</abbrev-journal-title>
<issn pub-type="epub">2296-424X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1393279</article-id>
<article-id pub-id-type="doi">10.3389/fphy.2024.1393279</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Physics</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Natural near field coupled leaky-mode resonant anti-reflection structures: the setae of <italic>Cataglyphis bombycina</italic>
</article-title>
<alt-title alt-title-type="left-running-head">Schwind et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fphy.2024.1393279">10.3389/fphy.2024.1393279</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Schwind</surname>
<given-names>Bertram</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2624440/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wu</surname>
<given-names>Xia</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Tiemann</surname>
<given-names>Michael</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Fabritius</surname>
<given-names>Helge-Otto</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Department of Chemistry</institution>, <institution>Paderborn University</institution>, <addr-line>Paderborn</addr-line>, <country>Germany</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Bionics and Materials Development</institution>, <institution>Department Lippstadt 1</institution>, <institution>Hamm-Lippstadt University of Applied Sciences</institution>, <addr-line>Hamm</addr-line>, <country>Germany</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Bruker Optik GmbH</institution>, <addr-line>Ettlingen</addr-line>, <country>Germany</country>
</aff>
<aff id="aff4">
<sup>4</sup>
<institution>Department Microstructure Physics and Alloy Design</institution>, <institution>Max-Planck-Institut f&#xfc;r Eisenforschung GmbH</institution>, <addr-line>D&#xfc;sseldorf</addr-line>, <country>Germany</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/301952/overview">S&#xe9;bastien Robert Mouchet</ext-link>, University of Exeter, United Kingdom</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2615579/overview">Branko Kolaric</ext-link>, University of Mons, Belgium</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2670334/overview">Diana Skigin</ext-link>, University of Buenos Aires, Argentina</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Bertram Schwind, <email>bertram.schwind@hshl.de</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>22</day>
<month>05</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>12</volume>
<elocation-id>1393279</elocation-id>
<history>
<date date-type="received">
<day>28</day>
<month>02</month>
<year>2024</year>
</date>
<date date-type="accepted">
<day>23</day>
<month>04</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Schwind, Wu, Tiemann and Fabritius.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Schwind, Wu, Tiemann and Fabritius</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Leaky mode resonances of the setae of <italic>Cataglyphis bombycina</italic> are found to enhance the thermal emission of the animals by near field coupling to the chitinous exoskeleton. This is remarkable, as the setae are also an adaption to enhance the reflectivity in the visible wavelength range. Both effects are dependent on morphology, dimensions and spatial arrangement. These parameters were experimentally characterized and simulated by finite difference time domain simulations to elucidate the optical impact of the setae in the mid infrared range and the contribution of leaky mode resonances. This mode of action and the setae&#x2019;s optical properties in the visible range explain evolutionary strains that led to the actual morphology and size of the setae.</p>
</abstract>
<kwd-group>
<kwd>
<italic>Cataglyphis bombycina</italic>
</kwd>
<kwd>leaky mode resonance</kwd>
<kwd>near field coupling</kwd>
<kwd>thermal emission enhancement</kwd>
<kwd>low refractive index materials</kwd>
<kwd>anti-reflection structures</kwd>
<kwd>FDTD</kwd>
</kwd-group>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Optics and Photonics</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>
<italic>Cataglyphis bombycina</italic> is a highly heat-tolerant ant species living in desert areas with extreme insolation. Since they forage during the hottest hours of the day, they have to sustain temperatures that are close to their critical thermal limit [<xref ref-type="bibr" rid="B1">1</xref>]. For survival, this species has evolved several traits, including long legs, heat shock proteins, visual orientation, and a dense coverage with hair-like setae [<xref ref-type="bibr" rid="B1">1</xref>&#x2013;<xref ref-type="bibr" rid="B3">3</xref>]. These setae have been shown to have amazing optical properties, helping the ant to deal with a wide spectrum of electromagnetic waves from sunlight. Nano-ribs on their upper surface in combination with a central canal efficiently scatter ultraviolet (UV) light [<xref ref-type="bibr" rid="B4">4</xref>]. The setae have cross sections in the shape of an isosceles right triangle with a hypothenuse of (3.4 &#xb1; 0.4) &#x3bc;m [<xref ref-type="bibr" rid="B4">4</xref>]. Because of this cross-sectional shape, they reflect visible light by total internal reflection [<xref ref-type="bibr" rid="B5">5</xref>, <xref ref-type="bibr" rid="B6">6</xref>] and scatter light in the visible range by higher order Mie-resonances [<xref ref-type="bibr" rid="B4">4</xref>, <xref ref-type="bibr" rid="B6">6</xref>]. In the mid infrared (MIR) the setae are suggested to form an antireflection layer, which enhances the radiative cooling of their body [<xref ref-type="bibr" rid="B6">6</xref>]. However, the size of the setae in that wavelength range is large enough to cause scattering due to low-order Mie-resonances [<xref ref-type="bibr" rid="B4">4</xref>]. Since the setae have an elongated shape, these resonances resemble resonances in optical wires. In optical wires these resonances are eigenmodes, that can be perfectly linked with Mie-resonances in dielectric structures [<xref ref-type="bibr" rid="B7">7</xref>]. These eigenmodes are also named leaky mode resonances (LMRs), since their electric field decays in the surrounding of the structure. Leaky modes determine optical properties of dielectric structures including absorption, scattering and emission. These properties can be engineered by changing, e.g., morphological details as well as substrate and superstrate [<xref ref-type="bibr" rid="B7">7</xref>]. In the natural system represented by the ant&#x2019;s integument, the distances between the setae and the cuticle are small compared to the wavelength of the thermal radiation. In this case, near field effects may play a role for radiative cooling [<xref ref-type="bibr" rid="B8">8</xref>, <xref ref-type="bibr" rid="B9">9</xref>]. Reid et al. showed that the radiative heat transfer between patterned and between un-patterned slabs increases by orders of magnitude if they are brought in close proximity, where <inline-formula id="inf1">
<mml:math id="m1">
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>&#x003c;</mml:mo>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mi>T</mml:mi>
</mml:msub>
<mml:mo>&#x003D;</mml:mo>
<mml:mi mathvariant="normal">&#x210f;</mml:mi>
<mml:mi>c</mml:mi>
<mml:mo>/</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> [<xref ref-type="bibr" rid="B9">9</xref>]. For this thermal radiation region, Shi et al. suggested that the reflectivity of the ant body is reduced by a refractive index gradient of an effective medium layer formed by the setae, which enhances the emissivity according to Kirchhoff&#x2019;s law, and thus enhances the radiative cooling of the ant body [<xref ref-type="bibr" rid="B6">6</xref>]. Since the setae, however, are spatially arranged with a defined distance from the underlying cuticle, the question can be raised whether effective medium theory is sufficient to explain the effects observed by [<xref ref-type="bibr" rid="B6">6</xref>]. Therefore, we experimentally characterized the arrangement of the setae and directly simulated the hemispherical emission of the system as a whole with finite difference time domain simulations (FDTD).</p>
</sec>
<sec id="s2" sec-type="materials|methods">
<title>2 Materials and methods</title>
<sec id="s2-1">
<title>2.1 Measurements of the distances between setae and cuticle</title>
<p>The spatial arrangement of the setae of the ants was investigated on cross-sectional samples of the cuticle. For this, small parts of the abdominal tergites were dissected and fixed in a mixture of 2.5% glutaraldehyde and 2% paraformaldehyde in H<sub>2</sub>O (pH 7.8). After washing the cuticle pieces 3 times in bi-distilled water (10&#xa0;min) they were post-fixed for 60&#xa0;min in an aqueous solution of 1% OsO<sub>4</sub> and 0.8% K<sub>4</sub> [Fe(CN)<sub>6</sub>]. Subsequently, the samples were washed a second time in bi-distilled water and dehydrated in a series of isopropanol solutions. After a final washing step in 100% acetone (4&#xa0;min), the cuticle pieces were embedded in EPON-resin. Cross-sectional block faces of the resin embedded cuticle were polished using a Leica RM 2155 Ultramicrotome and a diamond knife (Diatome, Switzerland). The polished cross sections were analyzed using a confocal laser scanning microscope (Olympus LEXT OLS3100) with a MPLAPON &#xd7;10 plan apochromat objective lens.</p>
</sec>
<sec id="s2-2">
<title>2.2 Estimating emission via FDTD reflectance simulations of the setae-covered cuticle</title>
<p>The emissivity can be investigated via simulations of the reflection <italic>R</italic>. Since the absorption <italic>A</italic> and the emission <italic>E</italic> are correlated via Kirchhoff&#x2019;s Law, the following holds for an opaque body:<disp-formula id="e1">
<mml:math id="m2">
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mo>&#x002B;</mml:mo>
<mml:mi>R</mml:mi>
<mml:mo>&#x003D;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
</p>
<p>The reflection of the ant&#x2019;s setae-cuticle system at an angle of incidence of 90&#xb0; was simulated with a classical 2D-FDTD approach using Meep [<xref ref-type="bibr" rid="B10">10</xref>] to investigate the influence of the setae in the MIR. The simulations were carried out with both TE and TM polarization. TE represents the case where the electric field is transversal to the long axis of the setae, while TM represents the case where the magnetic field is transversal to the long axis of the setae. The simulations were reduced to two dimensions with the assumption that the dimension of the setae&#x2019;s cross section is much smaller compared to its length. Therefore, the results represent the case of infinitely long setae.</p>
<p>Three scenarios were compared by simulating a plane wave pulse incident perpendicular to the cuticle surface to obtain the reflection. As first reference, the pulse was simulated in empty space for normalization purposes. As second reference, the pulse was directed towards a plain cuticle surface without setae. The reflection of the cuticle was calculated by the ratio of the results of these two simulations. The third scenario includes a layer of setae on top of the cuticle to simulate the reflection of the entire setae-cuticle system. Comparing these results with the plain cuticle can be used to investigate the influence of the setae on the reflection.</p>
<p>To investigate the influence of the gap between setae and cuticle as accurately as possible, randomized setups with 40 setae at five different mean distances <italic>d</italic>
<sub>mean</sub> to the cuticle were created (0&#xa0;&#x3bc;m, 0.38&#xa0;&#xb5;m, 0.63&#xa0;&#xb5;m, 1.0&#xa0;&#xb5;m, and 2.5&#xa0;&#xb5;m). The simulation results of 32 arrangements for each mean distance were averaged. The simulation parameters are summarized after the next section.</p>
</sec>
<sec id="s2-3">
<title>2.3 Simulation of the angle resolved emission of the setae covered cuticle</title>
<p>The emissive properties were additionally investigated by directly simulating the black body radiation of the ant with FDTD following the approach in [<xref ref-type="bibr" rid="B11">11</xref>]. The thermal radiation was simulated as a random noise <inline-formula id="inf2">
<mml:math id="m3">
<mml:mrow>
<mml:mover accent="true">
<mml:mi>K</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mfenced close=")" open="(" separators="&#x007C;">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> in the polarization (Eq. <xref ref-type="disp-formula" rid="e2">2</xref>) for dispersive materials [<xref ref-type="bibr" rid="B11">11</xref>].<disp-formula id="e2">
<mml:math id="m4">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mi>d</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mover accent="true">
<mml:mi>P</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:msup>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x002B;</mml:mo>
<mml:mi>&#x3b3;</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mover accent="true">
<mml:mi>P</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x002B;</mml:mo>
<mml:msubsup>
<mml:mi>&#x3c9;</mml:mi>
<mml:mn>0</mml:mn>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mover accent="true">
<mml:mi>P</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
<mml:mo>&#x003D;</mml:mo>
<mml:mi>&#x3c3;</mml:mi>
<mml:mover accent="true">
<mml:mi>E</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
<mml:mo>&#x002B;</mml:mo>
<mml:mover accent="true">
<mml:mi>K</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mfenced close=")" open="(" separators="&#x007C;">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
</p>
<p>In Meep, the random noise in the polarization is realized via a noisy Lorentzian susceptibility. Therefore, the complex material dispersion must be considered. For this, the absorption of chitin was experimentally determined and fitted with a Lorentzian dispersion function (see <xref ref-type="sec" rid="s12">Supplementary Material</xref>). The results are summarized in <xref ref-type="table" rid="T1">Table 1</xref> together with additional noise amplitudes simulating the respective emission temperature. These have been weighted by the Planck radiation of a black body at 53&#xb0;C. The electrical fields are created by the random material polarization. Therefore, no other electro-magnetic source was specified in these simulations.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Parameters of the fitted material dispersion of chitin used in the simulations.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">
<inline-formula id="inf3">
<mml:math id="m5">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> [cm<sup>&#x2212;1</sup>]</th>
<th align="center">3,306</th>
<th align="center">2,925</th>
<th align="center">1,650</th>
<th align="center">1,533</th>
<th align="center">1,413</th>
<th align="center">1,245</th>
<th align="center">1,064</th>
<th align="center">795</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">
<inline-formula id="inf4">
<mml:math id="m6">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b3;</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> [cm<sup>&#x2212;1</sup>]</td>
<td align="center">267</td>
<td align="center">131</td>
<td align="center">70</td>
<td align="center">54</td>
<td align="center">164</td>
<td align="center">76</td>
<td align="center">126</td>
<td align="center">97</td>
</tr>
<tr>
<td align="center">
<inline-formula id="inf5">
<mml:math id="m7">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">0.0069</td>
<td align="center">0.0023</td>
<td align="center">0.0123</td>
<td align="center">0.0064</td>
<td align="center">0.0104</td>
<td align="center">0.0061</td>
<td align="center">0.0175</td>
<td align="center">0.0152</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The electric field <inline-formula id="inf6">
<mml:math id="m8">
<mml:mrow>
<mml:mover accent="true">
<mml:mi>E</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mfenced close=")" open="(" separators="&#x007C;">
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> of the thermal radiation emitted by the ant was then monitored in space and time above the setae-cuticle system. To obtain the dependence on wavevector and frequency of the electric field according to [<xref ref-type="bibr" rid="B12">12</xref>], it was transformed by a 2-D Fourier-transformation:<disp-formula id="e3">
<mml:math id="m9">
<mml:mrow>
<mml:mover accent="true">
<mml:mi>E</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mfenced close=")" open="(" separators="&#x007C;">
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
<mml:mo>/</mml:mo>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x003D;</mml:mo>
<mml:mover accent="true">
<mml:mi>E</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mfenced close=")" open="(" separators="&#x007C;">
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>&#x3c9;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x003D;</mml:mo>
<mml:msubsup>
<mml:mo>&#x222b;</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>p</mml:mi>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:msubsup>
<mml:mo>&#x222b;</mml:mo>
<mml:mn>0</mml:mn>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:msub>
</mml:msubsup>
<mml:mover accent="true">
<mml:mi>E</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mfenced close=")" open="(" separators="&#x007C;">
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x22c5;</mml:mo>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>i</mml:mi>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>&#x22c5;</mml:mo>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>i</mml:mi>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>d</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
<p>Where the integration boundaries are defined by the unit cell size <italic>p</italic> and the simulation time <inline-formula id="inf7">
<mml:math id="m10">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. From the electric field in Eq. <xref ref-type="disp-formula" rid="e3">3</xref>, the intensity of the thermal radiation can be calculated in dependence of the wavelength and the angle of emittance <inline-formula id="inf8">
<mml:math id="m11">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mfenced close=")" open="(" separators="&#x007C;">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> using the relations (Eqs <xref ref-type="disp-formula" rid="e4">4</xref>&#x2013;<xref ref-type="disp-formula" rid="e6">6</xref>) following [<xref ref-type="bibr" rid="B12">12</xref>]:<disp-formula id="e4">
<mml:math id="m12">
<mml:mrow>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x003D;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
<disp-formula id="e5">
<mml:math id="m13">
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>&#x003D;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
<disp-formula id="e6">
<mml:math id="m14">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:mo>&#x003D;</mml:mo>
<mml:msup>
<mml:mover accent="true">
<mml:mi>E</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
<mml:mo>&#x2a;</mml:mo>
</mml:msup>
<mml:mo>&#x22c5;</mml:mo>
<mml:mover accent="true">
<mml:mi>E</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
</p>
<p>The unit cell must be large enough to prevent a time correlation between the left and the right boundary of the simulation cell. A schematic of the computation cell and the utilized geometry is given in <xref ref-type="fig" rid="F1">Figure 1</xref>. The simulation size was chosen to cover at least 40 randomly distributed setae (the randomization is summarized in the next section) on top of the cuticle, so the unit cell is large enough to employ periodic boundaries and make use of the equations above. The random thermal noise prevents a time correlation between corresponding points in neighboring unit cells. The incoherent thermal sources are thus advantageous for simulating the angle dependent emission.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Sketch of the computation cell for the simulation of the angle resolved emission in Meep, with grey shaded copys of itself on the left and right to illustrate the periodic boundary conditions.</p>
</caption>
<graphic xlink:href="fphy-12-1393279-g001.tif"/>
</fig>
<p>For the correct emission, a Lorentzian susceptibility is necessary. The absorption of chitin was fitted with the same Lorentz-oscillator model used for the reflection simulations. The thermal noise amplitudes of each Lorentz oscillator in <inline-formula id="inf9">
<mml:math id="m15">
<mml:mrow>
<mml:mover accent="true">
<mml:mi>K</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mfenced close=")" open="(" separators="&#x007C;">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> have been weighted by the Planck distribution of states (Eq. <xref ref-type="disp-formula" rid="e7">7</xref>) at a given temperature T to assign the correct temperature to the ant model [<xref ref-type="bibr" rid="B13">13</xref>]:<disp-formula id="e7">
<mml:math id="m16">
<mml:mrow>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mfenced close=")" open="(" separators="&#x007C;">
<mml:mrow>
<mml:mi>&#x3bd;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>d</mml:mi>
<mml:mi>&#x3bd;</mml:mi>
<mml:mo>&#x003D;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>8</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
<mml:msup>
<mml:mi>&#x3bd;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:msup>
<mml:mi>c</mml:mi>
<mml:mn>3</mml:mn>
</mml:msup>
</mml:mfrac>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mi>exp</mml:mi>
<mml:mrow>
<mml:mfenced close=")" open="(" separators="&#x007C;">
<mml:mrow>
<mml:mfrac bevelled="true">
<mml:mrow>
<mml:mi>h</mml:mi>
<mml:mi>&#x3bd;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>B</mml:mi>
</mml:msub>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>d</mml:mi>
<mml:mi>&#x3bd;</mml:mi>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>Where the most relevant temperature is the maximum temperature, the ant can survive T &#x003D; 326,15&#xa0;K &#x003D; 53&#xb0;C [<xref ref-type="bibr" rid="B3">3</xref>]. An average of 800 simulations with just the plain cuticle were used as reference. Since these simulations are based on thermal noise, longer simulation times are employed and up to 430 different runs with randomized setae arrangements were averaged. The emission enhancement caused by the setae is calculated by dividing the emission with setae by the reference without setae.</p>
<sec id="s2-3-1">
<title>2.3.1 Simulation parameters</title>
<p>The size of the computation cell was 16&#xa0;&#x3bc;m &#xd7; 60&#xa0;&#xb5;m bordered above and below by two perfect matching layers (PML) of 20&#xa0;&#xb5;m thickness. In Y-direction, periodic boundaries were employed. The cuticle has a thickness of 30&#xa0;&#xb5;m and the setae are represented as right-angled triangles with a hypotenuse of 3.5&#xa0;&#xb5;m length and a ribbed surface.</p>
<p>The randomization of the setae positions was realized by multiplication of the ordered coordinates with normal distributed random factors &#x3d5;<sub>x</sub> and &#x3d5;<sub>y</sub>, one for the X and one for the Y coordinate of the setae positions. The X coordinate determines the distance between adjacent setae and is always randomized. The Y coordinate determines the distance of each seta to the cuticle. Its randomization depends on the mean size of the gap <italic>d</italic>
<sub>mean</sub> multiplied with &#x3d5;<sub>y</sub>. The dimensions of the triangular setae cross sections were adapted from experimental measurements [<xref ref-type="bibr" rid="B4">4</xref>]. The complex refractive index of the seta material was experimentally determined [<xref ref-type="bibr" rid="B4">4</xref>] and fitted with a complex absorption function including 8 Lorentz resonance frequencies (see <xref ref-type="sec" rid="s12">Supplementary Material</xref>).</p>
<p>One unit in Meep was defined to represent 10&#xa0;&#xb5;m. For the reflection simulations a broadband planewave pulse and a resolution of 128 was used. For the emission simulations, only the noisy material dispersion was used as source. The computation cell and the number of simulations was increased; therefore, the resolution was reduced to 32. This is still reasonable, since the maximum frequency of an emitting body at 53&#xb0;C is lower than the maximum frequency of the pulse used for the reflectance simulations.</p>
<p>The imaginary part of the refractive index was fitted with a sum of 8 Lorentzian susceptibilities: <disp-formula id="e8">
<mml:math id="m17">
<mml:mrow>
<mml:mi>&#x3f5;</mml:mi>
<mml:mrow>
<mml:mfenced close=")" open="(" separators="&#x007C;">
<mml:mrow>
<mml:mi>&#x3c9;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x003D;</mml:mo>
<mml:mi>&#x3f5;</mml:mi>
<mml:mrow>
<mml:mfenced close=")" open="(" separators="&#x007C;">
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x003D;</mml:mo>
<mml:msub>
<mml:mi>&#x3f5;</mml:mi>
<mml:mi>&#x221e;</mml:mi>
</mml:msub>
<mml:mo>&#x002B;</mml:mo>
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mo>&#x003D;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mn>8</mml:mn>
</mml:munderover>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mi>n</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mi>n</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mi>f</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>i</mml:mi>
<mml:mi>f</mml:mi>
<mml:msub>
<mml:mi>&#x3b3;</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
</p>
<p>Where <inline-formula id="inf11">
<mml:math id="m18">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3f5;</mml:mi>
<mml:mi>&#x221e;</mml:mi>
</mml:msub>
<mml:mo>&#x003D;</mml:mo>
<mml:mn>2.4336</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, which corresponds to a refractive index of <italic>n</italic> &#x003D; 1.56, the refractive index of chitin [<xref ref-type="bibr" rid="B14">14</xref>]. The variables for the Lorentzian susceptibilities in Eq. <xref ref-type="disp-formula" rid="e8">8</xref> are given in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
</sec>
</sec>
</sec>
<sec id="s3" sec-type="results">
<title>3 Results</title>
<sec id="s3-1">
<title>3.1 Distribution of the setae on top of the cuticle</title>
<p>The distances between the cuticle and 42 setae distributed over 250&#xa0;&#xb5;m were measured on confocal laser scan images using the software ImageJ [<xref ref-type="bibr" rid="B15">15</xref>]. An exemplary section is shown in <xref ref-type="fig" rid="F2">Figure 2</xref>. The section shows the layered structure of the cuticle with the triangular setae distributed above in embedding resin. The cuticle and the embedding resin material are separated by a black area above the cuticle, which is a gap that opened up during microtome sectioning of the samples. Since the delamination of the embedding material and the cuticle happened after sectioning of the sample, the distances between individual setae were measured from the center of the hypothenuse to the cuticle surface and the delamination width was not included in the resulting distances.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Confocal laser-scanning microscopic image of the embedded cross sectioned cuticle of <italic>Cataglyphis bombycina</italic> with setae.</p>
</caption>
<graphic xlink:href="fphy-12-1393279-g002.tif"/>
</fig>
<p>The measured distances are summarized in the histogram in <xref ref-type="fig" rid="F3">Figure 3</xref>. Since the distance represents a positive random variable, it can be represented with a log-normal distribution (black curve). The according mean value is <italic>d</italic>
<sub>mean</sub> &#x003D; (0.68 &#xb1; 0.72) &#xb5;m.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Histogram of the measured distances between 42 setae and the cuticle surface. The histogram has three classes per micrometer. The according lognormal distribution is plotted in black and results in a mean value of <italic>d</italic>
<sub>
<italic>mean</italic>
</sub> &#x003D; (0.68 &#xb1; 0.72) &#xb5;m.</p>
</caption>
<graphic xlink:href="fphy-12-1393279-g003.tif"/>
</fig>
</sec>
<sec id="s3-2">
<title>3.2 FDTD reflection simulations using Kirchhoff&#x2019;s law for emission</title>
<p>First insights of the influence of the setae on the emissive properties of the ant are gathered by the reflection simulations under perpendicular incidence. In <xref ref-type="fig" rid="F4">Figures 4A, C</xref> the simulated reflection is plotted in the range from 0&#x2013;7,000&#xa0;cm<sup>&#x2212;1</sup>, in (A) for TM polarization (electric field parallel to the setae&#x2019;s long axis) and in (C) for TE polarization (electric field perpendicular to the setae&#x2019;s long axis). The data is averaged from simulations of 32 randomized setae arrangements for each mean gap to the cuticle.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>
<bold>(A,C)</bold> Simulated reflectance of a cuticle covered with setae [polarization: <bold>(A)</bold> TM; <bold>(C)</bold> TE], with different mean gaps between setae and cuticle. <bold>(B,D)</bold> The according emission enhancement [polarization: <bold>(B)</bold> TM; <bold>(D)</bold> TE]. In all graphs for comparison: simulated chitin reference without setae (blue dashed curve), wavenumber of the maximum of the Planck radiation at 53&#xb0;C (black dash-dot vertical line).</p>
</caption>
<graphic xlink:href="fphy-12-1393279-g004.tif"/>
</fig>
<p>The trends of both polarizations are comparable. They mostly differ in magnitude. The reflectance for TM polarization is generally higher than for TE. At 0&#xa0;cm<sup>&#x2212;1</sup>, the reflection is independent of the gap around 5% reflectance, which is the same reflectance as without presence of setae. With increasing wavenumbers all graphs undergo a minimum which is shifted to smaller wavenumbers the larger the gap between setae and cuticle becomes. At larger wavenumbers the reflection beyond this minimum increases the larger the gap between setae and cuticle is, and the absorption becomes more pronounced, e.g., around 1,600&#x2013;1,700&#xa0;cm<sup>&#x2212;1</sup>. Towards even larger wavenumbers, Fabry-Perot-like fringes form when a gap larger than 0&#xa0;&#xb5;m is present. The transition from low reflectance to the fringes, where all graphs with gaps coincide, takes place at larger wavenumbers the smaller the gap is. Interestingly, when using the formula to determine a film thickness from Fabry-Perot fringes: <inline-formula id="inf12">
<mml:math id="m19">
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>&#x003D;</mml:mo>
<mml:mfrac>
<mml:mi>m</mml:mi>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mrow>
<mml:mfenced close=")" open="(" separators="&#x007C;">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bd;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3bd;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msqrt>
<mml:mrow>
<mml:msup>
<mml:mi>n</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:msup>
<mml:mi>sin</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>, the result is a thickness <italic>d</italic> &#x223c; 1.75&#xa0;&#xb5;m, which corresponds to the height of the setae in the simulation. The parameters in the formula are: <italic>m</italic>, the number of fringes between two selected minima/maxima, with their frequencies <italic>&#x3bd;</italic>
<sub>1</sub>, <italic>&#x3bd;</italic>
<sub>2</sub> and the refractive index <italic>n</italic> under the angle of incidence <italic>&#x3b8;</italic>.</p>
<p>With Eq. <xref ref-type="disp-formula" rid="e1">1</xref> the emission can be calculated from the reflection, by assuming that <italic>E</italic> &#x003D; <italic>A</italic> via the reciprocity of light and Kirchhoff&#x2019;s law. In <xref ref-type="fig" rid="F4">Figures 4B, D</xref> the respective emission curves are plotted. The emission is normalized to the emission of the plain cuticle without setae. This results in an emission enhancement due to the coverage with setae. The minimum in the reflection is now a maximum depending on the gap size. The strongest enhancement can be found for the gap of 0.63&#xa0;&#xb5;m (orange) in the range between 0 and 1,500&#xa0;cm<sup>&#x2212;1</sup> in both polarizations. Interestingly, both polarizations have similar enhancement values in this range.</p>
<p>For the first three resonances found in <xref ref-type="fig" rid="F4">Figures 4A, B</xref> (at wavenumbers of <inline-formula id="inf13">
<mml:math id="m20">
<mml:mrow>
<mml:mover accent="true">
<mml:mi>&#x3bd;</mml:mi>
<mml:mo>&#x223c;</mml:mo>
</mml:mover>
<mml:mo>&#x003D;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> 670&#xa0;cm<sup>&#x2212;1</sup>, 2,900&#xa0;cm<sup>&#x2212;1</sup>, and 4,500&#xa0;cm<sup>&#x2212;1</sup>), the corresponding scattering fields of one seta are illustrated in <xref ref-type="fig" rid="F5">Figure 5</xref> for TM polarization. The resonances are indicated with TM<sub>ml</sub> similarly to [<xref ref-type="bibr" rid="B16">16</xref>], where <italic>m</italic> is the azimuthal mode number on the circumference of the seta and <italic>l</italic> is the number of field maxima within the seta. The scattering fields are shown for five different gap sizes between the cuticle and the seta for each wavenumber. The scattered electric fields were calculated by subtracting the electric field of the reflection of a plain cuticle without setae. Periodic boundaries apply for this illustration, so it represents an array of perfectly arranged neighboring setae. Therefore, the leaky mode resonances of the single setae are influenced by their neighbors and the cuticle below. Nevertheless, the leaky mode resonances can still be easily recognized such as also demonstrated in other literature examples [<xref ref-type="bibr" rid="B16">16</xref>&#x2013;<xref ref-type="bibr" rid="B18">18</xref>].</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Field distribution of the LMRs of the setae. Plotted is the scattering electric field of single setae on top of a cuticle with periodic boundaries in <italic>x</italic>-direction for TM-polarized light: <inline-formula id="inf14">
<mml:math id="m21">
<mml:mrow>
<mml:mover accent="true">
<mml:mi>E</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> parallel to the long axes of the setae (<italic>z</italic>-direction). The electric field of the reflection of the cuticle is subtracted (<inline-formula id="inf15">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi>E</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
<mml:mtext>setae</mml:mtext>
</mml:msub>
<mml:mo>&#x003D;</mml:mo>
<mml:msub>
<mml:mover accent="true">
<mml:mi>E</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
<mml:mtext>total</mml:mtext>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mover accent="true">
<mml:mi>E</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
<mml:mtext>cuticle</mml:mtext>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>). <bold>(A)</bold> TM<sub>01</sub> at 670&#xa0;cm<sup>&#x2212;1</sup> <bold>(B)</bold> TM<sub>11</sub> at 2,900&#xa0;cm<sup>&#x2212;1</sup> <bold>(C)</bold> TM<sub>21</sub> at 4,500&#xa0;cm<sup>&#x2212;1</sup>. Gap sizes (from left to right): 0, 0.38, 0.63, 1.0, and 2.5&#xa0;&#xb5;m for each LMR.</p>
</caption>
<graphic xlink:href="fphy-12-1393279-g005.tif"/>
</fig>
</sec>
<sec id="s3-3">
<title>3.3 FDTD simulations of the thermal emission of the ant</title>
<p>The simulations with a plane wave only represent the case for perpendicular incidence. Since the thermal emission occurs in all directions, the hemispherical emissivity of the ant has to be considered. This would require multiple simulations for many angles and wavelengths with the conventional FDTD simulation of the reflectivity. Therefore, a different and more direct approach was used to simulate the angle-resolved emission of the ant. <xref ref-type="fig" rid="F6">Figure 6</xref> shows three simulations of the angle resolved emission with mean gaps of 0&#xa0;&#x3bc;m, 0.63&#xa0;&#xb5;m, and 1.56&#xa0;&#xb5;m. The data is integrated over the frequency and normalized by reference simulations of a plain cuticle. Therefore, the results can be interpreted as an emission enhancement by the setae. Each graph is the average of 430 randomized setae arrangements and 800 reference runs were done. The according standard deviation is shown as confidence interval by the shaded areas around the graphs. As additional information the angle integrated values of the enhancement are listed in the legend entries, these can be seen as the amplification of the hemispherical emissivity of the ant through the presence of setae.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Simulation of the angle dependent emission enhancement for three different mean gaps between setae and cuticle. Data was integrated over the simulated wavelength range (0&#x2013;4,000&#xa0;cm<sup>&#x2212;1</sup>) and averaged over 430 randomized setae arrangements each. The graphs are normalized by the reference of a plain chitin surface that was averaged over 800 runs to reduce the confidence interval. The emission enhancement over all angles is given in percent within the legend.</p>
</caption>
<graphic xlink:href="fphy-12-1393279-g006.tif"/>
</fig>
<p>The emission enhancement with a gap of zero is meant to compare the natural setup of the ant with a typical anti-reflection coating. The emission enhancement of the surface coating at a gap size of 0&#xa0;&#xb5;m amounts to around 3% over the full angle range with little fluctuation and a slight trend to increase towards large angles. At a gap size of 0.63&#xa0;&#xb5;m, the emission of the ant increases over the whole angle of emission. Towards large exit angles, the enhancement is strongly elevated. Increasing the gap to 1.56&#xa0;&#xb5;m reduces the emission for small angles close to normal incidence compared to the previous two cases. The increase of the emission at large exit angles starts at slightly lower angle values and is even more pronounced than in the other simulations. Integrating the angle-resolved emission enhancement values results in an overall increase of the ant&#x2019;s emission of 3.2% if there is no gap between setae and cuticle, whereas the emission increases to 5.6% for a gap of 0.63&#xa0;&#xb5;m. Increasing distances between setae and cuticle further lowers the overall emission enhancement to 3.3%.</p>
<p>To verify the dependence of emission enhancement from setae-cuticle distance, further simulations were performed in a broader range of gap sizes up to 2.5&#xa0;&#xb5;m. To reduce computational efforts, only 72 randomized arrangements were simulated per gap size. This increases the confidence interval of these simulations, but the data is reproducible, and the trend of the emission enhancement is continuous. The integrated data with the according confidence intervals is plotted in <xref ref-type="fig" rid="F7">Figure 7A</xref> together with the data used for <xref ref-type="fig" rid="F6">Figure 6</xref> (smaller confidence intervals). The emission enhancement reaches a maximum of about 5.6% at a mean gap size of 0.63&#xa0;&#xb5;m. At smaller gap sizes, the emission enhancement drops to a value of around 3% at 0&#xa0;&#xb5;m gap size. Gap sizes over 0.63&#xa0;&#xb5;m reduce emission enhancement continuously down to a value close to 0% at 2.5&#xa0;&#xb5;m. In a second set of simulations, the dependence of the emission enhancement on the gap size was investigated for different setae dimensions.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Hemispherical emission enhancement in dependence of the average gap size between setae and cuticle. Simulated with various randomized setae arrangements (numerically integrated over emission angle and wavenumber) for <bold>(A)</bold>: setae with the natural size of <italic>Cataglyphis bombycina</italic> with confidence intervals and <bold>(B)</bold>: setae of different sizes.</p>
</caption>
<graphic xlink:href="fphy-12-1393279-g007.tif"/>
</fig>
<p>The results are shown in <xref ref-type="fig" rid="F7">Figure 7B</xref>. Each data point is an average of 72 randomized setae arrangements with the according mean gap. The confidence intervals are around &#xb1; 6% for all shown data points. For clarity reasons, they are not displayed as error bars. The data of the setae with natural dimensions (scaling factor 1.0; <xref ref-type="fig" rid="F7">Figure 7A</xref>) is included in the same color (orange; <xref ref-type="fig" rid="F7">Figure 7B</xref>). The emission enhancement of the system with smaller sized setae (scaling factor 0.5; blue dashed) has a maximum which is shifted to larger gaps compared to the original setae. For systems with larger setae (scaling factors 1.5, 2.0 and 3.0) this maximum is shifted to smaller gaps ultimately ending up at 0&#xa0;&#xb5;m and the classical behavior of an anti-reflection coating for scaling factors larger than 2. Comparison of the curves at a mean gap of 1&#xa0;&#xb5;m shows that the enhancement is highest for the original sized system. In the gap size range from 0.5&#xa0;&#xb5;m to 1&#xa0;&#x3bc;m, only the largest of the simulated setae (scaling factor 3.0) shows a better enhancement.</p>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<title>4 Discussion</title>
<p>Since the setae of <italic>C. bombycina</italic> have a very anisotropic aspect ratio, it is still a challenge to simulate them with reasonable computational efforts. Therefore, the setae-cuticle system was assumed to extend infinitely along the setae&#x2019;s long axis. The simulations can then be reduced to a 2D computation cell. Therefore, the emission can only be investigated along angles perpendicular to the setae&#x2019;s long axes. In this study, hemispherical is thus defined only via the angle <inline-formula id="inf16">
<mml:math id="m23">
<mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x003D;</mml:mo>
<mml:mrow>
<mml:mfenced close="]" open="[" separators="&#x007C;">
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mn>90</mml:mn>
<mml:mo>&#xb0;</mml:mo>
</mml:mrow>
<mml:mo>;</mml:mo>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mrow>
<mml:mn>90</mml:mn>
<mml:mo>&#xb0;</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, as visualized in <xref ref-type="fig" rid="F1">Figure 1</xref>. The emission was realized with a noisy material dispersion corresponding to the black body radiation at 53&#xb0;C. This corresponds to a wavenumber range of [0&#x2013;2,000&#xa0;cm<sup>&#x2212;1</sup>]. Thus, the simulations consider the maximum body temperature the ants can tolerate (T &#x003D; 326.15&#xa0;K &#x003D; 53&#xb0;C) during their foraging periods in the natural habitat of this species [<xref ref-type="bibr" rid="B3">3</xref>].</p>
<p>The angle resolved emission was simulated for three different systems with increasing gap sizes between setae and cuticle (<xref ref-type="fig" rid="F6">Figure 6</xref>). To demonstrate the effect of the setae on the emission, the three curves are normalized to a cuticle without setae (emission enhancement &#x003D; 1). In the first simulated case the setae are directly attached to the cuticle without any gap (black dashed, gap &#x003D; 0&#xa0;&#xb5;m). This already leads to an average increase in emission of about 3%. At small angles the emission enhancement is constant and only increases slightly at large angles. This mode of action corresponds to a classical anti-reflection surface [<xref ref-type="bibr" rid="B19">19</xref>]. In the second case the emission was simulated with a mean gap of 0.63&#xa0;&#xb5;m between setae and cuticle (orange solid) which corresponds closely to the geometry present in the ants with a mean distance of (0.68 &#xb1; 0.72) &#xb5;m [determined by experimental measurements (<xref ref-type="fig" rid="F2">Figures 2</xref>, <xref ref-type="fig" rid="F3">3</xref>)]. At small angles, the emission enhancement is almost like in the case without a gap, but at large angles it strongly increases. This results in an average emission enhancement of about 6%, which is almost twice as much as if there was no gap between setae and cuticle. If the gap is enlarged to more than double the size of the natural setup (blue dotted, gap &#x003D; 1.56&#xa0;&#xb5;m), the overall average emission enhancement decreases to about 3%. This is a similar value as in the first case, but the curves show different trajectories in similar angle ranges. At small angles, the emission enhancement is lower and at larger angles [sin(&#x3b8;) &#x003e; 0.7] it is much higher and rises with an even larger slope than in case two. The enhanced emission at large angles in case two and three can be explained by much better outcoupling of evanescent waves than in case one.</p>
<p>The lower emission enhancement at small angles in case three shows that when the distance between setae and the cuticle would be too large, the setae would become a hindrance for the thermal radiation. This is also supported by the reflectance simulations under perpendicular incident light (<xref ref-type="fig" rid="F4">Figures 4A, C</xref>).</p>
<p>The reflectance of the setup without a gap between setae and cuticle (dark blue; gap &#x003D; 0&#xa0;&#xb5;m) is lower than that of the plain cuticle (dashed blue) over the full displayed wavenumber range for both polarizations. Thus, the zero-gap case shows the typical behavior of an anti-reflection coating [<xref ref-type="bibr" rid="B19">19</xref>]. The situation is different if a gap between cuticle and the setae is present. In those cases, resonances occur, especially at large wavenumbers (e.g., between 4,000&#x2013;7,000&#xa0;cm<sup>&#x2212;1</sup>). These resonances increase the reflectance of TM as well as TE polarized light (<xref ref-type="fig" rid="F4">Figures 4A, C</xref>). In the spectral range below 2,000&#xa0;cm<sup>&#x2212;1</sup> however, the presence of setae results in a lower reflectance compared to a plain cuticle. Depending on the gap size between setae and cuticle, there is a spectral region where the reflectance is lower than if the setae were attached to the cuticle without gap. The larger the gap size, the narrower the respective region with low reflectivity. Additionally, the minimum in the reflectance is also shifted to lower wavenumbers. The larger the gap, the stronger is the influence of the absorption bands of chitin. They originate from the material dispersion included in the simulation, which is summarized in <xref ref-type="table" rid="T1">Table 1</xref>. The strongest ones occur between 1,500 and 2,000&#xa0;cm<sup>&#x2212;1</sup>.</p>
<p>An opaque body with a higher reflectance has by Eq. <xref ref-type="disp-formula" rid="e1">1</xref> a lower absorption. By the reciprocity of light, a lower absorption also leads to lower emission. Thus, a higher reflectance would cause a lower emission for the ant. Therefore, all gaps between setae and cuticle reduce the emission enhancement in the spectral range between 4,000&#x2013;7,000&#xa0;cm<sup>&#x2212;1</sup> (<xref ref-type="fig" rid="F4">Figure 4B</xref>). In this range, all curves of the simulations with a gap show a higher reflectance and similar resonances. This indicates that these resonances are independent of gap size and the randomized spatial arrangement of the setae. Therefore, it is safe to assume that the resonances originate solely from the setae and correspond to their leaky modes (LMR) or Mie-resonances. Interestingly, deriving a Fabry-Perot thickness from the periodicity of these resonances leads to a value d &#x223c; 1.75&#xa0;&#xb5;m, which corresponds very well to the height of the setae [<xref ref-type="bibr" rid="B4">4</xref>]. The presence of similar resonances was confirmed before with experimental FTIR transmission measurements of single setae [<xref ref-type="bibr" rid="B4">4</xref>]. Since the Planck radiation of a black body at 53&#xb0;C is almost zero in the wavenumber range of 4,000&#x2013;7,000&#xa0;cm<sup>&#x2212;1</sup>, the resonances in this region do not contribute to the thermal emission of the ant. For the thermal emission and resulting possible radiative cooling effects the range between 0&#x2013;1,500&#xa0;cm<sup>&#x2212;1</sup> is much more relevant. In the range of 0&#x2013;1,500&#xa0;cm<sup>&#x2212;1</sup>, the reflectance for all gap sizes decreases below the reflectance of the hypothetical anti-reflection coating represented by the zero-gap case. In this spectral region the emission is strongly influenced by the distance between cuticle and setae. Increasing these gap sizes shifts the maximum of the emission enhancement towards smaller wavenumbers. This effect can be explained by near field coupling of the LMRs. Three examples of these LMRs of setae in different distances to the cuticle are illustrated in <xref ref-type="fig" rid="F5">Figure 5</xref>.</p>
<p>
<xref ref-type="fig" rid="F5">Figure 5A</xref> shows the LMR TM<sub>01</sub> of a seta. Interestingly, it corresponds very well to the TM<sub>01</sub> of a cylinder, as demonstrated for a thin Germanium wire in [<xref ref-type="bibr" rid="B17">17</xref>]. Even the resonance frequency correlates very well to this literature example when considering the larger dimensions of the seta and the lower refractive index of chitin compared to the Germanium wire. Due to the lower refractive index of the seta, the leaky mode covers a broader wavenumber range (<xref ref-type="fig" rid="F4">Figure 4</xref>) and extends further into the surrounding space. Therefore, the LMR can couple to the cuticle and to neighboring setae even if the gap between seta and cuticle becomes quite large. Similar coupling effects have been measured on Si-nanowires, too [<xref ref-type="bibr" rid="B20">20</xref>].</p>
<p>In terms of optics, enhancing the thermal emission of the ant for its survival is a problem which is comparable to maximizing the absorption in nanowire antennas for harvesting solar radiation, although different wavelengths are involved. Cao et al. showed that there would be better structures than triangular prisms for enhancing the absorption, at least at resonances of higher order than TM<sub>01</sub> [<xref ref-type="bibr" rid="B17">17</xref>]. To maximize the photocurrent inside their nanowires, they optimized the absorption of the TM<sub>21</sub>/TE<sub>11</sub> resonances rather than TM<sub>01</sub>, because these are less sensitive on illumination angle. The TM<sub>01</sub> LMR of silicon nanowires has a low absorption cross section (<italic>Q</italic>
<sub>abs</sub>) [<xref ref-type="bibr" rid="B17">17</xref>], which limits the absorption in the wires, a disadvantage that would lower the photocurrent. Cao et al. theoretically showed that the TM<sub>01</sub> LMR would have the largest potential to enhance the photocurrent compared to unstructured silicon with the same thickness, at least for certain illumination angles [<xref ref-type="bibr" rid="B17">17</xref>]. However, a wire that could use TM<sub>01</sub> LMR would require a diameter of around 10&#xa0;nm, which is experimentally not feasible.</p>
<p>Despite the absorption bands in the MIR, which increase the refractive index of the chitinous setae, the optical nanowire antennas investigated by Cao have a much higher refractive index. Thus, the TM<sub>01</sub> LMR of the setae have an even lower <italic>Q</italic>
<sub>abs</sub> [<xref ref-type="bibr" rid="B4">4</xref>] compared to the nanowires. This turns into an advantage for the thermal emission problem of the ant. The setae cannot be treated as isolated antennas, since the underlying cuticle is an integral part of the body that always plays a role for the optical system of the ant. Thus, a large <italic>Q</italic>
<sub>abs</sub> is not necessary for the setae to enhance the outcoupling of the emitted thermal radiation of the cuticle. In this, the ant differs fundamentally from the nanowire antennas, where the substrate does not play an optical role. Therefore, near field heat transfer and near field coupling effects are much more relevant for the ant [<xref ref-type="bibr" rid="B9">9</xref>, <xref ref-type="bibr" rid="B20">20</xref>]. For such coupling, the TM<sub>01</sub> of the setae is beneficial, because of its larger spatial extension compared to higher order resonances or to TM<sub>01</sub> resonances of high refractive index materials like Ge.</p>
<p>In this study, the TM<sub>01</sub> LMR was only investigated for angles of emission perpendicular to the setae&#x2019;s long axis. However, for angles parallel to the long axis Cao et al. [<xref ref-type="bibr" rid="B17">17</xref>] experimentally characterized the LMRs of an isotropic nanowire and theoretically showed that there are only minor differences in the low order LMRs for nanowires of round, rectangular, hexagonal, and triangular cross-sectional shapes. The results in [<xref ref-type="bibr" rid="B17">17</xref>] show that the absorption cross section <italic>Q</italic>
<sub>abs</sub> is largest at normal incidence and decreases at larger angles of incidence. Especially the TM<sub>01</sub> resonance strongly decreases beyond 40&#xb0; to normal incidence and disappears at parallel incidence. The incidence angle on the setae in the setae-cuticle system is limited by the angle of total internal reflection at the cuticle-air surface, which is about 40&#xb0;; larger angles result in evanescent waves, which lead to parallel incidence on the setae. Since the TM<sub>01</sub> resonance of wire shaped objects (including objects with triangular cross sections like the setae of <italic>C. bombycina</italic>) is not excitable from these angles [<xref ref-type="bibr" rid="B17">17</xref>] the respective electric fields can no longer couple to the LMR. Thus, this scenario plays a minor role for the emission enhancement of the ant cuticle. Therefore, we focus on angles of emission perpendicular to the setae&#x2019;s long axis, thus avoiding expensive 3D-simulations of the complex architecture of the natural setae-cuticle system and still gaining the most important information on its mode of action in the MIR. The LMRs of the triangular setae, illustrated in <xref ref-type="fig" rid="F5">Figure 5A</xref>, are almost cylindrical even near the cuticle. Therefore, the TM<sub>01</sub> LMR is excitable over a broad illumination angle and the LMR can radiate into all the seta&#x2019;s azimuthal angles in free space. Thus, the coupled LMR extends the angle of emission of the ant cuticle. Thereby, the thermal emission of the ant is enhanced by outcoupling of evanescent waves. The TM<sub>01</sub> intercoupling of the randomly arranged setae and its coupling to the cuticle causes the system to act like an array of low-Q optical antennas [<xref ref-type="bibr" rid="B20">20</xref>] that emit in the MIR and offload thermal radiation, which helps to keep the ant cool.</p>
<p>
<xref ref-type="fig" rid="F5">Figures 5B, C</xref> show the higher order resonances (TM<sub>11</sub> and TM<sub>21</sub>). TM<sub>11</sub> corresponds well to simulations of triangular nanowires published by [<xref ref-type="bibr" rid="B17">17</xref>]. Since higher order resonances like TM<sub>21</sub> occur at shorter wavelengths, they are more confined to the seta. Therefore, they cannot couple that efficiently to the cuticle or neighboring setae like lower order resonances. Backscattering effects may occur and therefore the setae become a hinderance for the emitted radiation in these wavelength ranges. This is the reason why the emission depicted in <xref ref-type="fig" rid="F4">Figures 4B, D</xref> decreases in these wavelength ranges compared to the zero-gap situation.</p>
<p>The simulations shown in <xref ref-type="fig" rid="F4">Figures 4A, B</xref> reveal a certain wavenumber range where the near field-coupled TM<sub>01</sub> LMR have a lower reflectance than a classical anti-reflection coating. If the peak wavenumber of the Planck radiation at 53&#xb0;C (<inline-formula id="inf17">
<mml:math id="m24">
<mml:mrow>
<mml:mover accent="true">
<mml:mi>&#x3bd;</mml:mi>
<mml:mo>&#x223c;</mml:mo>
</mml:mover>
<mml:mo>&#x003D;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> 640&#xa0;cm<sup>&#x2212;1</sup>) is plotted in the diagram (black dash-dot vertical line), it becomes clear that the reflectance of the TM<sub>01</sub> is minimal in that wavenumber range. This increases the emission of a black body at an emission angle perpendicular to the cuticle surface. However, the whole hemispherical emission at 53&#xb0;C must be considered. This was done by a numerical integration over all angles within direct emission simulations (<xref ref-type="fig" rid="F7">Figure 7A</xref>), where the resulting emission enhancement through the setae coverage for a series of gaps have been meticulously computed. In accordance with <xref ref-type="fig" rid="F4">Figures 4B, D</xref> the maximum of the emission enhancement occurs at a mean gap size larger than 0&#xa0;&#xb5;m. The graph in <xref ref-type="fig" rid="F7">Figure 7A</xref> reveals a maximum of the emission enhancement at a mean gap of 0.63&#xa0;&#xb5;m. Notably, this gap size aligns very well with the mean of the log-normal distribution derived from experimentally measured distances between setae and cuticle in the prepared cross section (Histogram <xref ref-type="fig" rid="F3">Figure 3</xref>), at (0.68 &#xb1; 0.72) &#xb5;m. The match between the dimensions of the natural system and the gap with highest emission enhancement is a strong indication for an evolutionary adaptation of the setae-cuticle system.</p>
<p>The dependence of the emission on the gap size suggests that it results from the coupling of the setae&#x2019;s LMR to the cuticle. If this is the case, the trajectory of the curve in <xref ref-type="fig" rid="F7">Figure 7A</xref> should be significantly influenced by variations in setae dimensions as this would shift the resonances and influence their coupling to the cuticle and between each other. Since the setae of <italic>C. bombycina</italic> are also adapted to reflect incoming sunlight [<xref ref-type="bibr" rid="B4">4</xref>&#x2013;<xref ref-type="bibr" rid="B6">6</xref>], their dimensions are playing a crucial role in the visible and ultraviolet (UV) range [<xref ref-type="bibr" rid="B4">4</xref>]. Previous investigations have shown that the reflectance in the visible and the UV becomes better if the setae would be larger than their actual size [<xref ref-type="bibr" rid="B4">4</xref>]. The ability of the setae to enhance the emission in the MIR could be a possible evolutionary strain for their size restriction. To verify this, simulations of the emission of the setae-cuticle system were also done with varying setae sizes. This was done with a reduced number of simulations and only for a subset of gaps used for the data in <xref ref-type="fig" rid="F7">Figure 7A</xref> to limit the computational effort. In <xref ref-type="fig" rid="F7">Figure 7B</xref> the results of these simulations are summarized together with the curve from <xref ref-type="fig" rid="F7">Figure 7A</xref>. The results indicate that the maximum enhancement for smaller setae (scaling factor &#x003D; 0.5, blue dotted) would be localized at larger gaps, whereas for larger setae, the local maximum would be achieved at smaller gaps. At setae scaling factors of 2 (purple) and 3 (green) the emission enhancement decreases almost linearly with increasing gap size, similar to an anti-reflection coating, where the maximum of the emission enhancement is located at 0&#xa0;&#xb5;m gap size and leaky mode resonances play a minor role. This compellingly suggests that the setae-cuticle system represents an anti-reflection structure in the MIR through nearfield coupling resonant leaky-modes of the setae to the cuticle.</p>
</sec>
<sec id="s5" sec-type="conclusion">
<title>5 Conclusion</title>
<p>The detailed simulations presented in this paper lead to the conclusion, that the setae enhance the thermal emission of the ant in a more complex way then forming an effective medium. The setae not only form a Mie-resonant structure in the visible [<xref ref-type="bibr" rid="B4">4</xref>, <xref ref-type="bibr" rid="B6">6</xref>], but also in the MIR. Thus, they do not only shield the ant from visible and UV-light, but also enhance their ability for radiative cooling. This is crucial for the survival of the animals in their extreme desert habitat, where the body of the ants is heated up by solar radiation and the surrounding sand and air. This heat radiates from the body to the environment. Without setae, the emission angle is limited by the total internal reflection inside the cuticle. This leads to evanescent waves that are trapped inside the body. Through the presence of the setae, these waves can now couple to the first order LMRs of neighboring setae, thereby radiating the energy to the full hemispherical angle range. This leads to a thermal emission enhancement of about 6% compared to the plain cuticle. Varying the setae&#x2019;s size and their distances to the cuticle in our simulations shows, that the parameters realized by the ant are evolutionary optimized for enhancing the thermal emission. Considering the optical effects of the setae in the visible and UV range shows that dimensions, shape, surface structure, and average distance from the cuticle are evolutionary compromises to fulfill a large variety of functions that help this ant species to survive in its environment.</p>
<p>Although the optical functions described for visible and UV-light would be more efficient with larger setae [<xref ref-type="bibr" rid="B4">4</xref>], their actual size might be evolutionary restricted through their ability to enhance their thermal emission. The MIR effects are resonant-based, where the size is crucial and larger setae would become a hinderance for radiative cooling. LMRs in, e.g., nanowires have been shown to be an efficient way to enhance the absorption and the emission of light [<xref ref-type="bibr" rid="B7">7</xref>]. Radiative cooling does not require energy and is therefore a promising concept for a number of technical applications in the energy sector [<xref ref-type="bibr" rid="B21">21</xref>]. Here we show that a concept of radiative cooling based on near field coupled LMRs is realized in a biological system, which mainly differs from technical applications by consisting of low refractive index materials. Such systems can also combine high reflectance in the visible range with a high emittance in the MIR within a single structure. This type of structure has not been considered yet and might provide inspiration for new ways to use passive radiative cooling as a sustainable energy source.</p>
</sec>
</body>
<back>
<sec id="s6" sec-type="data-availability">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/<xref ref-type="sec" rid="s12">Supplementary Material</xref>, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s7">
<title>Ethics statement</title>
<p>Ethical approval was not required for the studies involving animals in accordance with the local legislation and institutional requirements because the specimens used were obtained from a laboratory that rear this ant species for their studies. The specimens used died of natural causes (age).</p>
</sec>
<sec id="s8">
<title>Author contributions</title>
<p>BS: Conceptualization, Investigation, Methodology, Visualization, Writing&#x2013;original draft, Writing&#x2013;review and editing. XW: Writing&#x2013;review and editing. MT: Writing&#x2013;review and editing. H-OF: Conceptualization, Supervision, Writing&#x2013;review and editing.</p>
</sec>
<sec id="s9" sec-type="funding-information">
<title>Funding</title>
<p>The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was supported by the DFG (Deutsche Forschungsgemeinschaft) Priority Program 1839 &#x201c;Tailored Disorder&#x2014;A science- and engineering-based approach to materials design for advanced photonic applications&#x201d; within project 12 &#x201c;Broadband reflecting fibers with tailored structures inspired by desert ants&#x201d;.</p>
</sec>
<ack>
<p>H-OF thanks Prof. Dr.-Ing. Dierk Raabe and the entire team of the Department Microstructure Physics and Alloy Design at the Max-Planck-Institut f&#xfc;r Eisenforschung GmbH for hosting this research and their continuous support. We thank Dr. Matthias Wittlinger and Prof. Dr. Harald Wolf from the Department Neurobiology of the University of Ulm for providing dead <italic>C. bombycina</italic> workers from their cultures.</p>
</ack>
<sec id="s10" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>Author XW was employed by Bruker Optik GmbH.</p>
<p>The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="s11" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s12">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fphy.2024.1393279/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fphy.2024.1393279/full&#x23;supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="DataSheet1.docx" id="SM1" mimetype="application/docx" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lenoir</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Aron</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Cerd&#xe1;</surname>
<given-names>X</given-names>
</name>
<name>
<surname>Hefetz</surname>
<given-names>A</given-names>
</name>
</person-group>. <article-title>Cataglyphis desert ants: a good model for evolutionary biology in Darwin&#x2019;s anniversary year - a review</article-title>. <source>Isr J Entomol</source> (<year>2009</year>) <volume>39</volume>:<fpage>1</fpage>&#x2013;<lpage>32</lpage>.</citation>
</ref>
<ref id="B2">
<label>2.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pfeffer</surname>
<given-names>SE</given-names>
</name>
<name>
<surname>Wahl</surname>
<given-names>VL</given-names>
</name>
<name>
<surname>Wittlinger</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Wolf</surname>
<given-names>H</given-names>
</name>
</person-group>. <article-title>High-speed locomotion in the Saharan silver ant, Cataglyphis bombycina</article-title>. <source>J Exp Biol</source> (<year>2019</year>) <volume>222</volume>:<fpage>jeb198705</fpage>. <pub-id pub-id-type="doi">10.1242/jeb.198705</pub-id>
</citation>
</ref>
<ref id="B3">
<label>3.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wehner</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Marsh</surname>
<given-names>AC</given-names>
</name>
<name>
<surname>Wehner</surname>
<given-names>S</given-names>
</name>
</person-group>. <article-title>Desert ants on a thermal tightrope</article-title>. <source>Nature</source> (<year>1992</year>) <volume>357</volume>:<fpage>586</fpage>&#x2013;<lpage>7</lpage>. <pub-id pub-id-type="doi">10.1038/357586a0</pub-id>
</citation>
</ref>
<ref id="B4">
<label>4.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schwind</surname>
<given-names>B</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>X</given-names>
</name>
<name>
<surname>Tiemann</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Fabritius</surname>
<given-names>H-O</given-names>
</name>
</person-group>. <article-title>Broadband Mie scattering effects by structural features of setae from the Saharan silver ant Cataglyphis bombycina</article-title>. <source>J Opt Soc Am B</source> (<year>2023</year>) <volume>40</volume>:<fpage>B49</fpage>. <pub-id pub-id-type="doi">10.1364/JOSAB.474899</pub-id>
</citation>
</ref>
<ref id="B5">
<label>5.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Willot</surname>
<given-names>Q</given-names>
</name>
<name>
<surname>Simonis</surname>
<given-names>P</given-names>
</name>
<name>
<surname>Vigneron</surname>
<given-names>J-P</given-names>
</name>
<name>
<surname>Aron</surname>
<given-names>S</given-names>
</name>
</person-group>. <article-title>Total internal reflection accounts for the bright color of the saharan silver ant</article-title>. <source>PLoS ONE</source> (<year>2016</year>) <volume>11</volume>:<fpage>e0152325</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0152325</pub-id>
</citation>
</ref>
<ref id="B6">
<label>6.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shi</surname>
<given-names>NN</given-names>
</name>
<name>
<surname>Tsai</surname>
<given-names>C-C</given-names>
</name>
<name>
<surname>Camino</surname>
<given-names>F</given-names>
</name>
<name>
<surname>Bernard</surname>
<given-names>GD</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>N</given-names>
</name>
<name>
<surname>Wehner</surname>
<given-names>R</given-names>
</name>
</person-group>. <article-title>Keeping cool: enhanced optical reflection and radiative heat dissipation in Saharan silver ants</article-title>. <source>Science</source> (<year>2015</year>) <volume>349</volume>:<fpage>298</fpage>&#x2013;<lpage>301</lpage>. <pub-id pub-id-type="doi">10.1126/science.aab3564</pub-id>
</citation>
</ref>
<ref id="B7">
<label>7.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Powell</surname>
<given-names>DA</given-names>
</name>
<name>
<surname>Padilla</surname>
<given-names>WJ</given-names>
</name>
<name>
<surname>Miroshnichenko</surname>
<given-names>AE</given-names>
</name>
</person-group>. <article-title>Resonant leaky modes in all-dielectric metasystems: fundamentals and applications</article-title>. <source>Phys Rep</source> (<year>2023</year>) <volume>1008</volume>:<fpage>1</fpage>&#x2013;<lpage>66</lpage>. <pub-id pub-id-type="doi">10.1016/j.physrep.2023.01.001</pub-id>
</citation>
</ref>
<ref id="B8">
<label>8.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rodriguez</surname>
<given-names>AW</given-names>
</name>
<name>
<surname>Ilic</surname>
<given-names>O</given-names>
</name>
<name>
<surname>Bermel</surname>
<given-names>P</given-names>
</name>
<name>
<surname>Celanovic</surname>
<given-names>I</given-names>
</name>
<name>
<surname>Joannopoulos</surname>
<given-names>JD</given-names>
</name>
<name>
<surname>Solja&#x10d;i&#x107;</surname>
<given-names>M</given-names>
</name>
<etal/>
</person-group> <article-title>Frequency-selective near-field radiative heat transfer between photonic crystal slabs: a computational approach for arbitrary geometries and materials</article-title>. <source>Phys Rev Lett</source> (<year>2011</year>) <volume>107</volume>:<fpage>114302</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevLett.107.114302</pub-id>
</citation>
</ref>
<ref id="B9">
<label>9.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Reid</surname>
<given-names>MTH</given-names>
</name>
<name>
<surname>Rodriguez</surname>
<given-names>AW</given-names>
</name>
<name>
<surname>Johnson</surname>
<given-names>SG</given-names>
</name>
</person-group>. <article-title>Fluctuation-induced phenomena in nanoscale systems: harnessing the power of noise</article-title>. <source>Proc IEEE</source> (<year>2013</year>) <volume>101</volume>:<fpage>531</fpage>&#x2013;<lpage>45</lpage>. <pub-id pub-id-type="doi">10.1109/JPROC.2012.2191749</pub-id>
</citation>
</ref>
<ref id="B10">
<label>10.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Oskooi</surname>
<given-names>AF</given-names>
</name>
<name>
<surname>Roundy</surname>
<given-names>D</given-names>
</name>
<name>
<surname>Ibanescu</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Bermel</surname>
<given-names>P</given-names>
</name>
<name>
<surname>Joannopoulos</surname>
<given-names>JD</given-names>
</name>
<name>
<surname>Johnson</surname>
<given-names>SG</given-names>
</name>
</person-group>. <article-title>Meep: a flexible free-software package for electromagnetic simulations by the FDTD method</article-title>. <source>Comput Phys Commun</source> (<year>2010</year>) <volume>181</volume>:<fpage>687</fpage>&#x2013;<lpage>702</lpage>. <pub-id pub-id-type="doi">10.1016/j.cpc.2009.11.008</pub-id>
</citation>
</ref>
<ref id="B11">
<label>11.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Luo</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Narayanaswamy</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Joannopoulos</surname>
<given-names>JD</given-names>
</name>
</person-group>. <article-title>Thermal radiation from photonic crystals: a direct calculation</article-title>. <source>Phys Rev Lett</source> (<year>2004</year>) <volume>93</volume>:<fpage>213905</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevLett.93.213905</pub-id>
</citation>
</ref>
<ref id="B12">
<label>12.</label>
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Kolle</surname>
<given-names>M</given-names>
</name>
</person-group> <article-title>Theoretical aspects of photonic structures</article-title>. In: <person-group person-group-type="editor">
<name>
<surname>Kolle</surname>
<given-names>M</given-names>
</name>
</person-group>, editor. <source>Photonic structures inspired by nature springer theses</source>. <publisher-loc>Berlin, Heidelberg</publisher-loc>: <publisher-name>Springer</publisher-name> (<year>2011</year>). p. <fpage>5</fpage>&#x2013;<lpage>27</lpage>. <pub-id pub-id-type="doi">10.1007/978-3-642-15169-9_2</pub-id>
</citation>
</ref>
<ref id="B13">
<label>13.</label>
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Sproull</surname>
<given-names>RL</given-names>
</name>
<name>
<surname>Phillips</surname>
<given-names>WA</given-names>
</name>
</person-group> <source>Modern physics: the quantum physics of atoms, solids, and nuclei</source>. <edition>3d ed.</edition> <publisher-loc>New York</publisher-loc>: <publisher-name>Wiley</publisher-name> (<year>1980</year>).</citation>
</ref>
<ref id="B14">
<label>14.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Leertouwer</surname>
<given-names>HL</given-names>
</name>
<name>
<surname>Wilts</surname>
<given-names>BD</given-names>
</name>
<name>
<surname>Stavenga</surname>
<given-names>DG</given-names>
</name>
</person-group>. <article-title>Refractive index and dispersion of butterfly chitin and bird keratin measured by polarizing interference microscopy</article-title>. <source>Opt Express</source> (<year>2011</year>) <volume>19</volume>:<fpage>24061</fpage>. <pub-id pub-id-type="doi">10.1364/OE.19.024061</pub-id>
</citation>
</ref>
<ref id="B15">
<label>15.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schneider</surname>
<given-names>CA</given-names>
</name>
<name>
<surname>Rasband</surname>
<given-names>WS</given-names>
</name>
<name>
<surname>Eliceiri</surname>
<given-names>KW</given-names>
</name>
</person-group>. <article-title>NIH Image to ImageJ: 25 years of image analysis</article-title>. <source>Nat Methods</source> (<year>2012</year>) <volume>9</volume>(<issue>9</issue>):<fpage>671</fpage>&#x2013;<lpage>5</lpage>. <pub-id pub-id-type="doi">10.1038/nmeth.2089</pub-id>
</citation>
</ref>
<ref id="B16">
<label>16.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cao</surname>
<given-names>L</given-names>
</name>
<name>
<surname>White</surname>
<given-names>JS</given-names>
</name>
<name>
<surname>Park</surname>
<given-names>J-S</given-names>
</name>
<name>
<surname>Schuller</surname>
<given-names>JA</given-names>
</name>
<name>
<surname>Clemens</surname>
<given-names>BM</given-names>
</name>
<name>
<surname>Brongersma</surname>
<given-names>ML</given-names>
</name>
</person-group>. <article-title>Engineering light absorption in semiconductor nanowire devices</article-title>. <source>Nat Mater</source> (<year>2009</year>) <volume>8</volume>:<fpage>643</fpage>&#x2013;<lpage>7</lpage>. <pub-id pub-id-type="doi">10.1038/nmat2477</pub-id>
</citation>
</ref>
<ref id="B17">
<label>17.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cao</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Fan</surname>
<given-names>P</given-names>
</name>
<name>
<surname>Vasudev</surname>
<given-names>AP</given-names>
</name>
<name>
<surname>White</surname>
<given-names>JS</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>Z</given-names>
</name>
<name>
<surname>Cai</surname>
<given-names>W</given-names>
</name>
<etal/>
</person-group> <article-title>Semiconductor nanowire optical antenna solar absorbers</article-title>. <source>Nano Lett</source> (<year>2010</year>) <volume>10</volume>:<fpage>439</fpage>&#x2013;<lpage>45</lpage>. <pub-id pub-id-type="doi">10.1021/nl9036627</pub-id>
</citation>
</ref>
<ref id="B18">
<label>18.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>P</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Z</given-names>
</name>
<name>
<surname>Takahara</surname>
<given-names>J</given-names>
</name>
</person-group>. <article-title>Multipole and multimode engineering in Mie resonance-based metastructures</article-title>. <source>Nanophotonics</source> (<year>2020</year>) <volume>9</volume>:<fpage>1115</fpage>&#x2013;<lpage>37</lpage>. <pub-id pub-id-type="doi">10.1515/nanoph-2019-0505</pub-id>
</citation>
</ref>
<ref id="B19">
<label>19.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Boden</surname>
<given-names>SA</given-names>
</name>
<name>
<surname>Bagnall</surname>
<given-names>DM</given-names>
</name>
</person-group>. <article-title>Optimization of moth&#x2010;eye antireflection schemes for silicon solar cells</article-title>. <source>Prog Photovoltaics</source> (<year>2010</year>) <volume>18</volume>:<fpage>195</fpage>&#x2013;<lpage>203</lpage>. <pub-id pub-id-type="doi">10.1002/pip.951</pub-id>
</citation>
</ref>
<ref id="B20">
<label>20.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cao</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Fan</surname>
<given-names>P</given-names>
</name>
<name>
<surname>Brongersma</surname>
<given-names>ML</given-names>
</name>
</person-group>. <article-title>Optical coupling of deep-subwavelength semiconductor nanowires</article-title>. <source>Nano Lett</source> (<year>2011</year>) <volume>11</volume>:<fpage>1463</fpage>&#x2013;<lpage>8</lpage>. <pub-id pub-id-type="doi">10.1021/nl1040429</pub-id>
</citation>
</ref>
<ref id="B21">
<label>21.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yin</surname>
<given-names>X</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Tan</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Fan</surname>
<given-names>S</given-names>
</name>
</person-group>. <article-title>Terrestrial radiative cooling: using the cold universe as a renewable and sustainable energy source</article-title>. <source>Science</source> (<year>2020</year>) <volume>370</volume>:<fpage>786</fpage>&#x2013;<lpage>91</lpage>. <pub-id pub-id-type="doi">10.1126/science.abb0971</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>