# Polymer photonic crystal membrane for thermo-regulating textile

Salim Assaf, Mohamed Boutghatin, Yan Pennec, V. Thomy, Alexander Korovin, A. Treizebre, Michèle Carette, Abdellatif Akjouj, Bahram Djafari-Rouhani

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Salim Assaf, Mohamed Boutghatin, Yan Pennec, V. Thomy, Alexander Korovin, et al., Polymer photonic crystal membrane for thermo-regulating textile. Scientific Reports, 2020, 10 (1), pp.9855. 10.1038/s41598-020-66731-1. hal-03142220

HAL Id: hal-03142220

https://hal.science/hal-03142220

Submitted on 11 Jul 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

OPEN

# Polymer photonic crystal membrane for thermo-regulating textile

Salim Assaf, Mohamed Boutghatin, Yan Pennec*, Vincent Thomy, Alexander Korovin, Anthony Treizebre, Michèle Carette, Abdellatif Akjouj & Bahram Djafari-Rouhani.

We study numerically the absorption and scattering properties of a polymer photonic membrane to thermoregulate the human body microclimate which corresponds to the area between the skin and a textile. We first show that the structuration of the absorbing photonic membrane with air holes leads to a modulation of the optical spectrum in the Mid-Infrared range. Indeed, we show that the membrane is able to modulate the transmission amplitude by 28% in benefit or deficit of both the absorption and reflection. We then studied the thermal balance between the human body and the surrounding environment through the photonic membrane. We found that, compared to a regular membrane, the photonic crystal structure behaves as a heating component that offers the possibility to reduce the temperature of the room up to +1 °C. The membrane is flexible, low cost, 3D-printable, free of metallic particles, and can easily be added to usual textiles.

For the past ten years, photonic nanostructures have represented a paradigm for the control of thermal radiations, offering a panel of exciting properties for energy applications¹,². Because of their abilities to control and manage electromagnetic waves at the Mid-Infrared (Mid-IR) wavelength scale, photonic nanostructures demonstrate their ability to manage thermal radiations properties in a way drastically different from conventional thermal emitters. Indeed, photonic crystal performances authorize to overtake the constraints of usual thermal emitters³ toward coherence, narrowband emission, polarization, directionality... The modification of Planck blackbody radiation by photonic band gap materials has been discussed, based on either simple photonic multilayered structures⁴ or more complex 1D⁵,⁶, 2D⁷,⁸ and 3D⁹ ones, thus making the demonstration of the enhancement, suppression and selectivity of thermal emission based on the photonic crystals properties. The fundamental advances in controlling thermal radiation led to different applications in the energy domain, as thermophotovoltaic devices that convert sunlight into thermal emission using hot absorber-emitter¹⁰. Although, the concept of daytime radiative cooling has been introduced with the objective to cool passively terrestrial structures by the use of a broad band photonic mirror for solar light, then emitting strongly in the Mid-IR within the atmospheric transparency window¹¹,¹².

More recently, another field of application has appeared in the thermal radiation control, with the introduction of photonic nanostructures in textiles for personal thermoregulation. The main motivation comes from the reduction of the building energy consumption of Heating, Ventilation and Air-Conditioning (HVAC) systems by locally providing heating or cooling in the human body close environment. To this end, personal thermoregulation photonic textiles have been proposed with the objective to obtain specific spectral properties. For personal heating textiles¹³,¹⁴, integration of metallic nanowires or particles embedded in polyethylene or cotton matrices has been proposed. The heating production is then generated from the human body radiation supplementary reflection that can even being increased by adding Joule heating to complement the passive insulation. For cooling purpose, recent papers have proposed different fibers¹⁵⁻¹⁸ or membranes¹⁹⁻²¹ fabrics to increase IR transparency. In the fiber configuration, Tong et al.¹⁵ proposed an infrared-transparent-visible-opaque-fabric (ITVOF) made of synthetic polymer fibers with intrinsically low IR absorbance and structured fibers to maximize IR transparency and visible opaqueness. For membranes fabrics, Hsu et al.¹⁹ develop experimentally a textile that promotes effective radiative cooling composed of nanoporous polyethylene (PE), transparent to Mid-IR human body radiations. The dual functionalities (cooling and heating) within a same textile has also been proposed in two different ways. First, Hsu et al.²² investigated a textile for human body radiation using a passive bilayer

Institute of Electronic, Microelectronic and Nanotechnology (IEMN), UMR 8520, University of Lille, Villeneuve d'Ascq, France. e-mail: yan.pennec@univ-lille.fr
<figure><img src="image_4.png" /><figcaption>**Figure 1.** (a) Schematic representation of the polymer membrane textile under human body radiation from the skin where (I), (R) and (T) represent respectively the incident, reflection and transmission coefficients of the electromagnetic waves. (b) Elementary unit cell used for the FE calculations with perfect matching layers (PML) and periodic boundary conditions (PBC). (c) 3D and in-plane view of the polymer membrane of thickness *h*, drilled with a triangular array of air holes with period *P* and diameter *D*.</figcaption></figure>

thermal emitter embedded inside an IR-transparent nano-PE that can perform both radiative heating and cooling using the same non-symmetric piece of textile. Second, reversible humidity sensitive clothing for personal thermoregulation was also proposed using shape memory polymer<sup>23</sup>. This smart textile has been designed to reversibly adapt the thermal insulation functionality, thus permitting the air flow and reducing the humidity level and the apparent temperature. Recently, Zhang *et al.*<sup>24</sup> have proposed a dynamic control the Mid-IR in reaction of the relative humidity of the underlaying skin by coating triacetate-cellulose bimorph fibers with a thin layer of carbon nanotubes.

The personal thermoregulation appears as an exciting scientific challenge to investigate, especially in the energy efficiency context. Many subjective parameters are attached to the thermoregulation of the human body linked to the apparent temperature and personal resentment<sup>25</sup>. Therefore, a tremendous effort is still necessary to develop smart wearable thermoregulating textiles to meet user demand for better thermal comfort. Until now, the previous papers have reported on transparent polymers (as PE) for cooling or on the introduction of metallic particles for heating effects. Complex geometry or properties for the dual heating and cooling purposes has been also proposed. In the present paper, we aim to propose a simple micro-structured photonic membrane for the heating functionality. Our purpose is to take advantage of the modulation of the electromagnetic waves at the Mid-Infrared (Mid-IR) to increase the human body temperature. This is done through the structuration of the membrane, as compared to the regular structure. To reach this objective, we propose to take benefit of the absorption properties of the polymer together with the scattering properties of the photonic membrane to thermoregulate the temperature of the membrane. We chose the benzocyclobutene (BCB) for the demonstration but the results can be easily extended to every polymer containing such an absorption in the Mid-IR range.

In the first section of the paper we investigate the optical properties of the polymer membrane, considering the effect of its structuration on the reflection/transmission/absorption coefficients in the Mid-IR range, and analyze the origin of the specific features that occur. In the second part, we define and analyze the thermal balance between the human body and the indoor environment through the photonic membrane, considering the radiation, convection and conduction mechanisms.

## Optical properties of the BCB polymer membrane

**Model and method.** The investigation of the properties of a photonic membrane in the Mid-IR has been done numerically with the help of the Finite Element (FE) method. As seen schematically Fig. 1(a), we consider an incident electromagnetic waves (I) radiating in air from the human body skin and normally interacting with the photonic polymer membrane. All calculations have been done considering the polarization (E<sub>y</sub>, H<sub>x</sub>) of the incident light. Nevertheless, due to the honeycomb structuration of the photonic crystal array, the results can be extended to the other polarization. We record the transmission (T) and reflection (R) coefficients and then deduce the absorption coefficient following the equation A = 1 - R - T. The source and the detector are placed in air, respectively before and after the membrane. The structure of the elementary 3D unit cell used for the FE calculation is shown in Fig. 1(b). The incident source is a plane wave generated in the air medium which propagates along the z-direction, perpendicularly to the periodic structure. Along the x and y direction, periodic boundary conditions (PBC) are applied on each side of the unit cell to build the periodic structure. As the structure is supposed infinite along the z direction, perfect matching layers (PMLs) are applied at the finite boundaries of the box, along z, in order to avoid any reflections of outgoing waves. All calculations have been performed on
<figure><img src="image_5.png" /><figcaption>**Figure 2.** (a) Representation of the refractive index components (real (n) in black and imaginary (k) in blue) in the Mid-IR range [5-15] µm. (b, c) Reflection (black), transmission (blue) and absorption (red) spectra for the (b) non-structured and (c) structured BCB membrane with the geometrical parameters: *P* = 7.0 µm, *D* = 5.5 µm, *h* = 4.0 µm. (d) Zoom of the reflection curve in the wavelength range [5, 8] µm. The blue hatched area represents the black body emissivity of the human body, calculated from the Planck law at 34 °C. (e) (Left) Schematic representation of the unit cell in the (y, z) plane. The plane wave is launched from the bottom of the unit cell. (Right) Snapshots of the modulus of the electric field E, at the wavelength λA = 6.08 µm (A), λB = 6.41 µm (B) and λC = 6.88 µm (C).</figcaption></figure>

a periodic triangular array of holes in the membrane. The geometrical parameters involved in the study are the lattice parameter *P*, the hole's diameter *D*, and the thickness of the membrane *h* (Fig. 1(c)).

We aim to develop a photonic crystal membrane that can modulate the Mid-IR radiations of the human body
radiations. To reach the objective, we propose the study of a Benzocyclobutene (BCB) membrane, currently used
in microelectronic manufacturing processes. As seen Fig. 2(a), the optical characterization of the BCB mem-
brane, shows a variation of the real (n) and imaginary (k) part of the refractive index in the Mid-IR. Because of
its non-polar chemical structure, its average refractive index of 1.57 makes it a good candidate to ensure a suffi-
cient refractive-index contrast with air together with an absorption in the frequency range of interest (see Supp.
Information 1 for more details).

**Spectral analysis of the photonic membrane in the Mid-IR.** We first present Fig. 2(b), as a reference, the reflection, transmission and absorption coefficient spectra of a non-structured BCB membrane, having a thickness of *h* = 4.0 µm. At wavelengths close to 8.0 µm, 9.5 µm and 12.0 µm, several dips appear in the transmission spectrum (blue line), corresponding to absorption peaks in the spectrum (red line). These peaks are directly linked to the extinction coefficient of the refractive index of the BCB (Fig. 2(a)) and are thus attributed to the absorption of the polymer. One can note that all peaks appear in the human body emittance, calculated from the Planck law at the skin temperature of the human body *T*<sub>s</sub> = 34 °C, represented within the blue hatched area. A small reflection (black curve), less than 20%, is recorded with a maximum at 8.5 µm and a zero at 6.0 µm. This modulation corresponds to Fabry-Perot oscillations through the BCB membrane of finite thickness.

We now consider the structured membrane, following a triangular array of air holes, with the geometrical
parameters *P* = 7.0 µm, *D* = 5.5 µm, and *h* = 4.0 µm. Figure 2(c) shows the corresponding calculation of the reflec-
tion (black), transmission (blue) and absorption (red) spectra. We first find the absorption peaks of the BCB
polymer, close to 8.0 µm, 9.5 µm and 12.0 µm, described previously. Qualitatively, it can be seen that the average
level of transmission increases to the detriment of the reflection and absorption. Then, one can see the occurrence
of new features appearing at low wavelengths, between 6.0 µm and 7.5 µm. To get a higher wavelength resolution,
the reflection spectrum is magnified in the Fig. 2(d). We can remark that the high reflection occurs at two main
peaks, at 6.08 µm and 6.88 µm, and a small one at 6.41 µm in the middle. To understand their origin, we performed
calculations of the modulus *E* of the electric field at the corresponding wavelengths (Fig. 2e). One can see that the
input signal, launched from the bottom, does not transmit through the membrane. In each case, the origin comes
from the excitation of a stationary mode which belongs to the BCB membrane and couple to the incident wave.
The two high reflection peaks (A and C) correspond to modes spread over the BCB membrane and are attributed
to guided modes in the plate: one is antisymmetric (λA = 6.08 µm) and the other is symmetric (λC = 6.88 µm) with
respect to the mid-plane of the membrane. The small peak in the middle (B) corresponds to an antisymmetric
<figure><img src="image_6.png" /><figcaption>**Figure 3.** (a) Evolution of the reflection (R), transmission (T) and absorption (A) coefficients as a function of the wavelength in the Mid-IR for three photonic membranes, designed following the scaling factors $α_1, α_4$ and $α_9$. (b) Evolution of the efficiency coefficient $η$, associated to the reflection (black), transmission (blue), absorption (red) responses of the non-structured (dashed lines) and structured (solid lines) BCB membrane.</figcaption></figure>

mode strongly confined inside the air holes. One can notice that the interaction of this localized mode with the incident wave propagating in the continuum, gives rise to a peak of asymmetric shape in the reflection spectrum. Such wave interaction phenomenon, known as a Fano resonance, is widely reported in the literature$^{26,27}$

**Tuning of the spectral coefficients in the Mid-IR.** The objective is now to tune the spectral coefficients in the Mid-IR range. To reach this objective, we proceed to a systematic calculation of the R, T and A coefficients as a function of the geometrical parameters of the photonic membrane. We showed that the reflection, transmission and absorption curves are affected by the structuration of the membrane, more particularly in the lower part of the emissive spectrum. To see the impact of the geometric parameters on the emissivity of the human body, we repeat the same calculations when all the geometrical dimensions of the membrane are multiplied by a scaling factor $α_i$, where $i = 1$ corresponds to the geometrical parameters chosen in the previous section as a reference. Indeed, such a homothetic variation of the geometrical parameters, following a scale law, induces a similar scaling in the wavelength of the radiation effect of the membrane, thus emerging inside the human body emissivity. The supplementary information 2 reports the set of parameters associated to the scaling factors $α_i$.

The calculated reflection, transmission and absorption coefficients are shown in Fig. 3(a) for three scaling factors $α_1, α_4$ and $α_9$. In the reflection spectrum, we found that, increasing the size of the geometrical parameters leads to the shift of the three reflective features A, B and C, towards higher wavelengths. In the transmission and absorption spectra, one can qualitatively follow the signature of the reflective modes, that fully cross the human body emissivity, and drastically affects both spectra.

To get a quantitative estimation of these behaviors, we defined an efficiency coefficient, $η$, corresponding to the integration of the R, T and A coefficients over a defined wavelength range corresponding to the human body emissivity. Such a coefficient is commonly used for the estimation of photovoltaic solar cells efficiency$^{28}$. This coefficient is explained through the formula:

$$ \eta = \frac{\int_{\lambda_{min}}^{\lambda_{max}} E_{\lambda} \cdot \chi_{\lambda} \cdot d\lambda}{\int_{\lambda_{min}}^{\lambda_{max}} E_{\lambda} \cdot d\lambda} \qquad (1) $$

where $E_λ$ is the human body emissivity and $χ_λ$ is one of the R, T or A coefficients at the wavelength $λ$. The integration will be done over the wavelength interval [$λ_{min} = 7.5 \mu m - λ_{max} = 11.5 \mu m$] that cover 70% of the human body emissivity with a maximum at $λ = 9 \mu m$.

Figure 3(b) reports the evolution of the efficiency coefficient $η$ for the reflection (black solid line), transmission (blue solid line) and absorption (red solid line) as a function of the scaling factors $α_i$. One can see that, increasing $α_0$ to $α_3$, the transmission drops by 28%, i.e. from 94% to 66% in benefit of the reflection and absorption which increase respectively by 16% (from 0.5% to 16%) and 12% (from 5.5% to 17.5%). The efficiency coefficient is then almost constant until the geometrical parameters reach $α_6$ where the transmission increases again. As demonstrated previously, such tunability of the efficiency coefficients is directly related to the occurrence of the reflective peaks A, B and C over the wavelength interval. The tunability has been compared with the one of a non-structured membrane of same thickness (dashed lines). We clearly see that the structuration of the
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Evolution of the reflection (black), transmission (blue), and absorption (red) efficiency coefficients $\eta$ in the range [7.5–11.5] µm, as a function of: (a) the diameter of air holes, keeping constant the period $P = 10$ µm, and the thickness $h = 5.7$ µm of the BCB membrane, (b) the thickness of the membrane for $P = 10$ µm, $D = 7.9$ µm, and (c) the incident angle $\theta$ of the incoming human body radiation for the scaling factor $\alpha_4$. The solid (resp. dashed) lines correspond to the structured (resp. non structured) responses.</figcaption></figure>

membrane made possible the tunability of the transmission coefficient, from transparent to opaque in the Mid-IR, by increasing the geometrical factor $\alpha_i$. We then obtained a photonic membrane with both a transparency ($\alpha_i \le \alpha_1$) and an opacity ($\alpha_i > \alpha_3$), higher by respectively ≈14% and ≈13% compared to an unstructured membrane.

Moreover, these tunabilities are supported by a complementary calculation where we change one geometrical parameter at a time while keeping the two others constant. In Fig. 4(a), the period $P = 10$ µm and thickness $h = 5.7$ µm are kept constant while we investigate the effect of the variation in the diameter D. When the diameter of the holes is small we reproduce the behavior of a non-structured membrane which is 80% transparent to the human body radiation. For a diameter $D = 7.9$ µm, corresponding to the scaling factor $\alpha_4$, we recover a decrease of the transmission down to 66%. In between, the spectral coefficients can still be tuned, for instance we obtain an increased opacity of 15% for $D = 7$ µm. We have also investigated the effect of the thickness $h$ of the membrane, keeping constant $P = 10.0$ µm and $D = 7.9$ µm (Fig. 4(b)). One can see that the minimum of transmission is obtained when the reflection is maximum, i.e. for $h = 5.7$ µm. Finally, the low transmission has been investigated as a function of the angle of incidence of the human body radiations for the scaling factor $\alpha_4$. The numerical results (Fig. 4(c)) reveals a robustness of the transmission opacity with respect to the incidence angle at least until 10 degrees which is acceptable with respect to the hypothesis of an incident plane wave coming from the human body radiation.

## Thermal investigation

The objective of this section is to examine the effect of the optical properties found in section 1 on the thermal flows between the human body and the surrounding environment, through the membrane, and to demonstrate the heating power of the structured membrane. The comparison is done with a non-structured equivalent membrane.

**Thermal balance analysis.** To quantify the impact and benefit of the photonic BCB membrane on personal thermoregulation, we rely on a one-dimensional steady-state heat transfer model analysis including radiative, conductive, and convective heat transfer mechanisms as shown in Fig. 5.

In the following, we provide a summary of the assumptions, input parameters, and analytical formula useful for the thermal calculation. Throughout this work we limit ourselves to a functionality of the membrane for everyday use, assuming that the human body is in a sedentary state, inside a room. As a consequence, uniform skin temperature and heat generation from the human body are considered. The air between the skin and the polymer membrane is assumed to be stationary, meaning that convective heat transfer will be neglected in this region. In the model, both skin and environment are assumed to be respectively an ideal blackbody emitter and absorber. The membrane will then partially reflect, transmit and absorb the incoming optical radiation, in support of the optical results obtained in section 1. We used Kirchhoff's law and define the emissivity $\epsilon_m(\lambda)$ of the membrane throughout its absorption coefficient, as:

$$ \epsilon_m(\lambda) = A = 1 - R - T \qquad (2) $$

where A, R and T are respectively the absorbance, reflectance and transmittance of the polymer membrane calculated in section 1.
<figure><img src="image_8.png" /><figcaption>Figure 5. Schematic representation of the heat transfer model of the clothed human body to the surrounding environment, considering all radiative (rad.), conductive (cond.) and convective (conv.) mechanisms. The indices s, a, mo and mi respectively mean skin, ambient air, membrane inner and membrane outer surfaces. d is the distance between the human body and the membrane.</figcaption></figure>

The equation of the thermal balance can be defined considering two volumes of control, respectively embedding the human body (CV1) and the BCB membrane (CV2) (Fig. 5). By calculating the incoming and outgoing flows passing through these volumes, the energy balance at the CV1 and CV2 can be expressed as:

$$q_{\text{gen}} - q_{\text{rad,s}}(1-R) - q_{\text{cond}} + T \cdot q_{\text{rad,a}} + q_{\text{rad,mi}} = 0 \quad (3)$$

$$q_{\text{cond}} + A \cdot q_{\text{rad,s}} + A \cdot q_{\text{rad,a}} - [q_{\text{rad,mi}} + q_{\text{rad,mo}} + q_{\text{conv}}] = 0 \quad (4)$$

where $q_{gen}=73W.m^{-2}$ is the metabolic body heat production rate per unit area in a middle-aged man standing and relaxed, $q_{rad,s}$ is the radiative heat flux from the skin, $q_{rad,a}$ is the radiative heat flux from the ambient air, $q_{rad,mi}$ is the radiative heat flux from the polymer membrane inner surface, $q_{rad,mo}$ is the radiative heat flux from the polymer membrane outer surface, $q_{cond}$ is the conductive heat flux in the air gap between the skin and the textile, $q_{conv}$ is the convective heat flux from the membrane to the ambient air.

The conductive, convective, and radiative heat flux terms follow respectively the Fourier, Newton, and Stefan-Boltzmann laws and take the following expressions:

$$q_{\text{rad,s}} = \sigma \cdot T_s^4 \quad (5)$$

$$q_{\text{rad,a}} = \sigma \cdot T_a^4 \quad (6)$$

$$q_{\text{rad,mo}} = \varepsilon_m(\lambda) \cdot \sigma \cdot T_2^4 \quad (7)$$

$$q_{\text{rad,mi}} = \varepsilon_m(\lambda) \cdot \sigma \cdot T_1^4 \quad (8)$$

$$q_{\text{cond}} = k_{\text{air}} \frac{T_s - T_1}{d} \quad (9)$$

$$q_{\text{conv}} = h_c(T_2 - T_a) \quad (10)$$

where $T_s$ is the temperature of the skin ($T_s = 34$ °C), $T_a$ is the temperature of the ambient air, $T_1$ and $T_2$ are respectively the temperatures of the inner and outer surfaces of the membrane. $k_{air}$ is the thermal conductivity of air ($k_{air} = 0.027 W.m^{-1}.K^{-1}$), $h_c$ is the convective heat transfer coefficient ($h_c = 3.0 W.m^{-2}.K^{-1}$), and $\sigma$ is the Stefan-Boltzmann constant ($\sigma = 5.6710^{-8} W.m^{-2}.K^{-4}$). Note that the convective heat transfer coefficient $h_c$ can be adapted to simulate different environmental conditions as, for example, the air circulation within the room (not considered here). Additionally, the air gap thickness, $d$, has been fixed to $d = 2$ mm and can be adapted as well.

To these two Eqs. (3) and (4), we have to add a third equation describing the heat transfer (conduction) mechanism through the membrane, $q_{cond,m}$. From the temperature profile within the textile demonstrated in ref. 16 (see SI, Eq. S22), we are able to define a relation between the temperatures $T_1$ and $T_2$ respectively at the inner and outer faces of the membrane as:

$$T_2 = \frac{h}{2k_m} \{ \varepsilon_m \sigma (T_1)^4 + \varepsilon_m \sigma (T_2)^4 - \varepsilon_m \sigma T_s^4 - \varepsilon_m \sigma T_a^4 \} - \frac{k_{\text{air}}h}{k_m d} (T_s - T_1) + T_1 \quad (11)$$
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Evolution of (a) the temperature of the membrane T<sub>BCB</sub> = T<sub>1</sub> = T<sub>2</sub> and (b) the ambient temperature as a function of the scaling factor α<sub>i</sub> for a structured (solid lines) and non-structured (dashed lines) membrane.</figcaption></figure>

where *k*<sub>m</sub> = 0.2 W·m<sup>-1</sup>·K<sup>-1</sup> is the BCB's thermal conductivity<sup>29</sup>, *h* is the thickness of the membrane, and *d* is the thickness of the microclimate. Actually, due to very small thickness of our membrane, it is expected that the temperatures T<sub>1</sub> and T<sub>2</sub> will be almost the same as we shall see numerically in the next section.

**Determination of the BCB and ambient temperatures.** The efficiency of the membrane is demonstrated by fixing the temperature of the skin, typically around T<sub>s</sub> = 34°C and to determine the unknown temperatures T<sub>1</sub>, T<sub>2</sub>, and T<sub>a</sub> from the resolution of Eqs. (3), (4) and (11).

In Fig. 6(a,b), we show the evolution of the temperatures T<sub>1</sub>, T<sub>2</sub> and T<sub>a</sub> as a function of the scaling factor α<sub>i</sub>. To estimate the effect of the structuration of the membrane, we compare the results (solid lines in Fig. 6) with those of the non-structured membrane of equivalent thickness (dashed lines). First, one can notice that the inner and outer temperatures of the membrane are exactly the same, reported as T<sub>BCB</sub> in Fig. 6a. This is due to the extreme thinness of the membrane, typically smaller than 10 µm. Then, when the scaling factor increases, the temperature varies from 32.4 °C to 31.75 °C, reaching a minimum for α<sub>i</sub> = α<sub>5</sub>. This modulation of temperature follows the variation in transmission (see Fig. 3b) realized in benefit of the reflection and absorption of the membrane. Finally, if we consider the geometrical parameters of the membrane associated to α<sub>5</sub>, one can see that the temperature is lower in the structured membrane as compared to the non-structured membrane of equivalent thickness by about 0.25 °C. This means that, to keep the skin at the comfort temperature of 34 °C, a lower temperature is needed in the structured membrane as compared to the unstructured one. Therefore, as we shall see from Fig. 6b, a lower ambient temperature will be sufficient to provide the same comfort.

Indeed, another way to demonstrate the efficiency of the structuration is to represent the ambient temperature calculated from Eqs. (3), (4) and (11) as a function of the scaling factor (Fig. 6b). This representation allows to compare the necessary temperature of the room to keep the thermal comfort of the human body. For the design α<sub>5</sub>, the calculated ambient temperature can be up to +1 °C lower for the structured membrane as compared to the non-structured one. This property offers the opportunity of maintaining the individual comfort while reducing the external energy supply of the room.

We also calculated the set of temperatures (T<sub>1</sub>, T<sub>2</sub>, T<sub>a</sub>) for different skin temperatures slightly lower than 34°C (see SI, table 3). Figure 7 shows the evolution of the skin temperature as a function of the ambient temperature for the structured and unstructured membrane. One can clearly see that the temperature of the skin is systematically 1°C higher when the human body is clothed with the structured membrane. It demonstrates the heating capacity of the photonic membrane. A variation of -1°C would represent a significant reduction in energy consumption, leading to economic and environmental benefits.

**Conclusion**

We have demonstrated the thermoregulation properties of a photonic crystal membrane based on the modulation of the human body optical radiations in the Mid-IR range. We first study the transmission, reflection and absorption properties of the polymer membrane in the Mid-IR domain. The polymer under consideration is the Benzocyclobutene (BCB) which presents both an absorbance and sufficient refractive-index contrast with air to deal with the photonic crystal properties following a triangular array of air holes. We found that the structuration leads to new features in the spectral reflection, transmission and absorption spectra that have been identified as guided and localized modes of the photonic crystal membrane. We then introduced a quantitative coefficient to estimate the efficiency of the transmission of the human body emissivity through the photonic membrane. We found that, by scaling the geometrical parameters of the crystal, the transmission can be modulated up to 28% of the human body emissivity in the wavelength range [7.5, 11.5] µm.

We then coupled the optical properties with thermal balance applied between the human body and the surrounding environment, through the structured membrane. We found that for a set of geometrical parameters (P
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Evolution of the skin temperature as a function of the ambient temperature for the structured (solid line) and unstructured (dashed line) membrane (from SI, table 3).</figcaption></figure>

The membrane behaves as a thermoregulator of the human body. We found that, compared to a non-structured regular membrane of equivalent thickness, the temperature of the room needs to be systematically 1 °C less when the human body is clothed with the structured membrane in order to get the same thermal comfort ($T_s$ = 34 °C). We believe that such a simple photonic membrane offers a new opportunity for the thermal management of textile for everyday use, and targets one of the major domains in energy consumption. Additionally, the presence of pores drilling the membrane enables air permeability and can promote water-wicking.

Received: 19 December 2019; Accepted: 18 May 2020;

Published online: 17 June 2020

References

1. Sawhney, A. P. S. et al. Modern Applications of Nanotechnology in Textiles. *Text. Res. J.* **78**, 731–739 (2008).

2. Li, W. & Fan, S. Nanophotonic control of thermal radiation for energy applications [Invited]. *Opt. Express* **26**, 15995 (2018).

3. Kittel, C. & Kroemer, H. Thermal physics. (W. H. Freeman, 1980).

4. Cornelius, C. M. & Dowling, J. P. Modification of Planck blackbody radiation by photonic band-gap structures. *Phys. Rev. A* **59**, 4736–4746 (1999).

5. Foley, J. M., Itsuno, A. M., Das, T., Velicu, S. & Phillips, J. D. Broadband long-wavelength infrared Si/SiO$_2$ subwavelength grating reflector. *Opt. Lett.* **37**, 1523 (2012).

6. Foley, J. M., Young, S. M. & Phillips, J. D. Narrowband mid-infrared transmission filtering of a single layer dielectric grating. *Appl. Phys. Lett.* **103**, 071107 (2013).

7. Liu, X. et al. Taming the Blackbody with Infrared Metamaterials as Selective Thermal Emitters. *Phys. Rev. Lett.* **107**, (2011).

8. Viallon, M. et al. Modulation of the refractive properties of 1D and 2D photonic crystal polycrystalline silicon-based membranes in the MIR frequency range. *J. Phys. Appl. Phys.* **52**, 205101 (2019).

9. Lin, S.-Y. et al. Enhancement and suppression of thermal emission by a three-dimensional photonic crystal. *Phys. Rev. B* **62**, R2243–R2246 (2000).

10. Lenert, A. et al. A nanophotonic solar thermophotovoltaic device. *Nat. Nanotechnol.* **9**, 126–130 (2014).

11. Rephaeli, E., Raman, A. & Fan, S. Ultrabroadband Photonic Structures To Achieve High-Performance Daytime Radiative Cooling. *Nano Lett.* **13**, 1457–1461 (2013).

12. Raman, A. P., Anoma, M. A., Zhu, L., Rephaeli, E. & Fan, S. Passive radiative cooling below ambient air temperature under direct sunlight. *Nature* **515**, 540–544 (2014).

13. Hsu, P.-C. et al. Personal Thermal Management by Metallic Nanowire-Coated Textile. *Nano Lett.* **15**, 365–371 (2015).

14. Yu, Z., Gao, Y., Di, X. & Luo, H. Cotton modified with silver-nanowires/polydopamine for a wearable thermal management device. *RSC Adv.* **6**, 67771–67777 (2016).

15. Tong, J. K. et al. Infrared-Transparent Visible-Opaque Fabrics for Wearable Personal Thermal Management. *ACS Photonics* **2**, 769–778 (2015).

16. Catrysse, P. B., Song, A. Y. & Fan, S. Photonic Structure Textile Design for Localized Thermal Cooling Based on a Fiber Blending Scheme. *ACS Photonics* **3**, 2420–2426 (2016).

17. Jafar-Zanjani, S., Salary, M. M. & Mosallaei, H. Metafabrics for Thermoregulation and Energy-Harvesting Applications. *ACS Photonics* **4**, 915–927 (2017).

18. Gao, T. et al. Three-Dimensional Printed Thermal Regulation Textiles. *ACS Nano* **11**, 11513–11520 (2017).

19. Hsu, P.-C. et al. Radiative human body cooling by nanoporous polyethylene textile. *Science* **353**, 1019–1023 (2016).

20. Song, Y.-N. et al. Wearable Polyethylene/Polyamide Composite Fabric for Passive Human Body Cooling. *ACS Appl. Mater. Interfaces* **10**, 41637–41644 (2018).

21. Cai, L. et al. Spectrally Selective Nanocomposite Textile for Outdoor Personal Cooling. *Adv. Mater.* **30**, 1802152 (2018).

22. Hsu, P.-C. et al. A dual-mode textile for human body radiative heating and cooling. *Sci. Adv.* **3**, e1700895 (2017).

23. Zhong, Y. et al. Reversible Humidity Sensitive Clothing for Personal Thermoregulation. *Sci. Rep.* **7**, (2017).

24. Zhang, X. A. et al. Dynamic gating of infrared radiation in a textile. *Science* **363**, 619–623 (2019).
25. steadman.pdf.

26. Limonov, M. F., Rybin, M. V., Poddubny, A. N. & Kivshar, Y. S. Fano resonances in photonics. *Nat. Photonics* **11**, 543–554 (2017).

27. Amoudache, S. et al. Optical and acoustic sensing using Fano-like resonances in dual phononic and photonic crystal plate. *J. Appl. Phys.* **119**, 114502 (2016).

28. Zhou, D. et al. Optimization of the optical properties of nanostructured silicon surfaces for solar cell applications. *J. Appl. Phys.* **115**, 134304 (2014).

29. Pinel, S. et al. Thermal modeling and management in ultrathin chip stack technology. *IEEE Trans. Compon. Packag. Technol.* **25**, 244–253 (2002).

## Acknowledgements

This work was supported by the European Interreg project 'Photonitex'. All authors thank the company SEMILAB for technical assistance in the characterization of the BCB refractive index. One author (S.A.) thanks the company DAMART for partial PhD financial support.

## Author contributions

M.B. and S.A. are PhD students involved in the calculation. V.T., A.T. and M.C. are collaborator for experimental discussions about the potential fabrication and characterization of the proposed structure, including the nature of the material or the geometrical stresses. The objective of the paper is to propose a theoretical work that is realistic experimentally. A.K., A.A. and B.D.R. are involved in the analytical and numerical thermal and optical calculations corresponding to the main issue of the ms. Y.P. managed the whole work. All authors reviewed the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41598-020-66731-1.

Correspondence and requests for materials should be addressed to Y.P.

Reprints and permissions information is available at www.nature.com/reprints.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020