# Elaboration and Characterization of Active Films Containing Iron-Montmorillonite Nanocomposites for O2 Scavenging

Erland-Modeste Kombaya-Touckia-Linin, Sebastien Gaucel, Moulay Tahar Sougrati, Lorenzo Stievano, Nathalie Gontard, Valérie Guillard

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

## ► To cite this version:

Erland-Modeste Kombaya-Touckia-Linin, Sebastien Gaucel, Moulay Tahar Sougrati, Lorenzo Stievano, Nathalie Gontard, et al. Elaboration and Characterization of Active Films Containing Iron-Montmorillonite Nanocomposites for O2 Scavenging. Nanomaterials, 2019, 9 (9), 10.3390/nano09091193. [hal-02285618](https://hal.science/hal-02285618)

HAL Id: hal-02285618

https://hal.science/hal-02285618

Submitted on 12 Sep 2019

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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>

Article

# Elaboration and Characterization of Active Films Containing Iron-Montmorillonite Nanocomposites for O₂ Scavenging

Erland-Modeste Kombaya-Touckia-Linin ¹D, Sébastien Gaucel ¹D, Moulay T. Sougrati ²D, Lorenzo Stievano ²D, Nathalie Gontard <fnref n="1" /> and Valérie Guillard <fnref n="1" />,*

Received: 25 June 2019; Accepted: 16 August 2019; Published: 23 August 2019

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

check for
updates

**Abstract:** Iron particles of sizes between 6 and 20 nm forming aggregates of 57 ± 17 nm were synthesized by chemical reduction of iron precursors on the surface of montmorillonite (MMT). This active MMT-Fe powder was then uniformly distributed in a linear low-density polyethylene (LLDPE) matrix by extrusion at atmospheric conditions, as confirmed by wide-angle X-ray scattering (WAXS), which also detected a partial exfoliation of the nanoclays. Thermogravimetric analysis (TGA) did not detect any significant modification of the degradation temperature between nanocomposites and active nanocomposites. ⁵⁷Fe Mössbauer spectroscopy evidenced the formation of a majority of iron boride in MMT-Fe as well as in the active film containing it. The LLDPE.Fu15.MMT-Fe3.75 and LLDPE.Fu15.MMT-Fe6.25 films had oxygen-scavenging capacities of 0.031 ± 0.002 and 0.055 ± 0.009 g(O₂)/g(Fe), respectively, while the neat powder had an adsorption capacity of 0.122 g(O₂)/g(Fe). This result confirms that the fresh film samples were partially oxidized shortly after thermomechanical processing (60% of oxidized species according to Mössbauer spectroscopy). No significant difference in oxygen permeability was observed when MMT-Fe was added. This was related to the relatively small film surface used for measuring the permeability. The reaction-diffusion model proposed here was able to reproduce the observed data of O₂ adsorption in an active nanocomposite, which validated the O₂ adsorption model previously developed for dried MMT-Fe powder.

**Keywords:** iron nanoparticles; iron boride; montmorillonite; nanocomposites; oxygen adsorption kinetics; modeling; ⁵⁷Fe Mössbauer spectroscopy

## 1. Introduction

Much attention is currently being paid to the development of O₂ barrier monolayer films that are as performing as the multilayer complex materials currently used for commercial applications and can still maintain the end of life characteristics of a monolayer pure matrix (recyclability, biodegradability, etc.). Indeed, multilayer packaging materials that nowadays offer the required protection against O₂ and oxidation are difficult to recycle due to their complex formulations and thus have very high environmental impact [1–3]. Moreover, they are expensive and represent an economic burden for the food industry.

To obtain monolayer films with good O₂ barrier properties, one strategy is to develop materials based on nanocomposites such as, for instance, clay-polymer nanocomposites. Thanks to the dispersion of nanoscale clays in the bulk of the polymer, a significant decrease of O₂ permeability can be achieved

<footnote n="1">¹ UMR, Ingénierie des Agropolymères et Technologies Emergentes, INRA, Univ. Montpellier, Montpellier SupAgro, CIRAD, 34060 Montpellier, France</footnote>

<footnote n="2">² Institut Charles Gerhardt Montpellier, Univ. Montpellier, CNRS, 34090 Montpellier, France</footnote>

<footnote>* Correspondence: valerie.guillard@umontpellier.fr</footnote>
compared with that of a pure polymer matrix [4–6]. This is observed when a highly exfoliated structure is obtained, that is, when the polymer penetrates into the layered structure of the clay platelets, taking them apart and leading to individual clay foils randomly dispersed in the polymer. The observed improvement in barrier properties can be ascribed to the extension of the length and tortuosity of the diffusing gas flow paths, which leads to a reduction of the diffusivity of the gas through the polymer nanocomposite [6–8]. In practice, however, highly exfoliated structures are difficult to achieve, and the expected decrease in O₂ permeability is not always as important as expected, as recently explained by Wolf et al. [6].

To improve the O₂ barrier performance of monolayer films, another strategy relies on the introduction of active agents reacting with the O₂ that permeates from the surroundings [7–9]. The most widely used O₂ scavengers are based on powders containing polyamide, sulfite, and iron species. Among iron-based systems, zerovalent iron nanoparticles (ZVINPs) have recently gained attention for their application in active packaging. ZVINPs are used for their very high surface reactivity in environmental applications (remediation of inorganic contaminants in water and soil, e.g., U, Cd, Zn, NO₃⁻, etc.) [10–13] and medical applications (detoxification of biological fluids, tissue healing, etc.) [14]. Their use in packaging applications to produce materials with active oxygen-scavenging properties is much more recent and still less studied. A previous work from Kombaya-Touckia-Linin et al. (2019) [15] recently demonstrated the feasibility of synthesizing via a liquid route low-size spherical ZVINPs of about 50 nm in diameter directly on a nanoclay support. These authors provided a thorough characterization of their iron-based nanoscale oxygen scavengers, including full characterization of the O₂ adsorption kinetics. The as-synthetized active powder was found to have an adsorption capacity of 122 mg of O₂ per gram of absorber (91.7 mL of O₂ per gram of iron). The main advantage of using ZVINPs synthetized on nanoclay supports is that they can be easily included in the bulk polymer using classical thermomechanical processes [15,16]. In such composites, both a passive barrier (provided by the tortuosity induced by nanoclays) and an active effect (due to O₂-scavenging compounds) are expected with enhancement of the global performance of the material. Moreover, ZVINPs in such composites are entrapped in the bulk polymer and thus do not have a potential mobility to migrate into food preventing human exposure when used as food contact materials. It is generally assumed that consumers will not be exposed to nanoparticles from food contact polymers as long as their size exceeds 4 nm in at least one of their dimensions and the contact surface is not altered by mechanical surface stress during application [17].

Only three articles could be found in the scientific literature dealing with active nanocomposites with ZVINPs as active oxygen scavengers: that of Khalaj et al. [18], who prepared nanocomposite polypropylene films containing ZVINPs anchored on montmorillonite (MMT) via a thermomechanical process; the work of Foltynowicz et al. [7], who synthesized ZVINPs and dispersed them in a silicon matrix; and the paper of Busolo and Lagaron [9], who evaluated the oxygen-scavenging potential of a synthetic ZVINPs-containing kaolinite for food packaging plastics. In all these studies, ZVINPs are able to intercept and trap oxygen with adsorption capacities varying from 2.3 to 4.3 mL O₂ per gram of active film composite. However, in all these previous studies, the lack of full structural characterization of these ZVINPs systems did not permit the authors to carry out deep, mechanistic analyses of the adsorption of O₂ and to provide a precise relationship between iron speciation and the observed adsorption properties.

The current work had the original ambition of describing, from a mechanistic point of view, the phenomenon of O₂ diffusion–reaction within an active film containing nanoscale iron-based oxygen scavengers. This goal was sought via the production of active nanocomposites containing well-defined and characterized ZVINPs supported on nanosized MMT-Fe according to the MMT-Fe synthesis previously set up by Kombaya-Touckia-Linin et al. (2019) [15]. These nanoparticles were dispersed in a linear low-density polyethylene (LLDPE) matrix via an extrusion process. LLDPE was chosen as a model polymer since it is widely used in the packaging industry and its thermal, mechanical, and barrier properties are well known. Its processing with nanoclays using classical thermomechanical
processes is well characterized too, and the exfoliation of MMT is well mastered, especially when polyethylene-grafted maleic anhydride (PE-g-MA) is added to improve the compatibility between the clay and polymer matrix [19,20]. The synthesized composites were characterized by transmission electron microscopy (TEM), wide-angle X-ray scattering (WAXS), differential scanning calorimetry (DSC), thermogravimetric analysis (TGA), and $^{57}$Fe Mössbauer spectroscopy. The oxygen adsorption rate and maximal adsorption capacity were experimentally measured and correlated with the iron speciation in the film and its evolution with time. The oxygen adsorption kinetic was successfully modeled using a two-species kinetic reaction. Finally, the O$_{2}$ permeability of the final materials was also evaluated.

## 2. Materials and Methods

### 2.1. Materials

LLDPE (LL 1002YB melt flow index 2.0 g·10 min$^{-1}$, density 0.918 g·cm$^{-3}$) was supplied by Exxon Mobil Chemicals (Germany). MMT naturally modified by alkyl ammonium (Cloisite 20A: CEC of 95 meq·100 g$^{-1}$ clay mineral and a basal distance value $d_{001}$ = 2.65 nm) was obtained from Southern Clay Products Inc.

LDPE-g-MA (Fusabond E226: maleic-anhydride-modified polyethylene with a melt flow index of 1.75 g·10 min$^{-1}$ and density of 0.93 g·cm$^{-3}$), used to enhance the compatibility between the nanoclay and polymer matrix, was kindly provided by Dupont™ (Germany).

High-purity reagents, including ethanol (99%), NaOH (32% aqueous solution), FeCl$_{3}$·6H$_{2}$O (98%), and sodium borohydride (NaBH$_{4}$), were purchased from Sigma-Aldrich (Germany) and used as received. Distilled and demineralized water was used throughout these experiments.

#### 2.1.1. Synthesis of MMT-Fe

Montmorillonite-supported zerovalent iron particles (MMT-Fe) were obtained through the liquid route by the reduction of ferric iron with sodium borohydride according to a previous work [15]. Briefly, for the preparation of each batch, 10.1 g of FeCl$_{3}$·6H$_{2}$O was dissolved in 500 mL of a solution containing 1:4 v/v of deionized water and absolute ethanol. MMT (10.5 g) was added to the FeCl$_{3}$ solution and the suspension was stirred for 40 h at room temperature.

The iron-reducing solution (7.5 g of NaBH$_{4}$ in 500 mL of a basic aqueous solution made of 475 mL of deionized water and 25 mL of 32% NaOH) was stirred for 16 h. Use of a basic solution to dissociate sodium borohydride prevents the production of hydrogen [21]. The FeCl$_{3}$-MMT suspension and the NaBH$_{4}$ solution were both degassed in flowing argon for 2 h to remove the dissolved oxygen.

Dropwise, 250 mL of the NaBH$_{4}$ solution were added to the FeCl$_{3}$-MMT suspension under vigorous stirring (400 rpm). The molar ratio of borohydride to Fe was 4:1. Borohydride is expected to reduce Fe(III) to Fe(0) according to Equation (1):

$$ \mathrm{Fe}(\mathrm{H}_{2}\mathrm{O})_{6}^{3+} + 3\mathrm{BH}_{4}^{-} + 3\mathrm{H}_{2}\mathrm{O} \rightarrow \mathrm{Fe}^{0} + 3\mathrm{B}(\mathrm{OH})_{3} + 10.5\mathrm{H}_{2}. \quad (1) $$

After the addition of NaBH$_{4}$, the solution turned black, indicating the reduction of ferric iron. Synthesis was done at room temperature. A continuous argon flow was used in the reactor to remove the oxygen. A constant stirring speed was maintained throughout the synthesis.

The solid recovered by filtration was washed three times in absolute ethanol (99%), dried in flowing argon for 24 h at room temperature, and then stored in anoxic conditions before use for a maximum of one month before disposal. The nanostructured MMT-Fe final product was used for active nanocomposite design.
### 2.1.2. Processing of MMT-Fe/LLDPE Films

The active polyolefin composites were prepared by mixing a masterbatch of LLDPE containing either 9 or 15 wt % of LLDPE-g-MA with 3.75 or 6.25 wt % of MMT-Fe, respectively, using a corotating twin-screw extruder (Thermo Scientific™ EuroLab 16, Germany) in atmospheric conditions with a L/D ratio of 40 and a screw diameter of 16 mm at a screw speed of 200 rpm and feed rate of 1.0 kg·h⁻¹. The temperature profile of the extruder from the feeder to the die increased regularly from 160 to 180 °C, as described in detail elsewhere [19,20]. The maximum temperature was set at 180 °C to prevent decomposition of the modified MMT. Control samples containing pure LLDPE as well as samples containing iron-free MMT were also prepared in the same conditions. No specific precautions were taken to protect the active material in order to mimic the processing conditions of the industrial operation environment. The extruder was filled regularly and continuously with a very small quantity of active components to avoid its overexposure to atmospheric O₂ in the feeder (<1 min in the presence of atmospheric O₂). Uniform thin films were produced using a flat die at the end of the extruder, which was followed by calendaring. The as-produced active films took less than 10 min to be cold, cut, and stored in anoxic conditions. Composition and labeling of the prepared active films are summarized in Table 1.

<table><caption>**Table 1.** Labeling and composition of the active films prepared.</caption><thead><tr><th rowspan="2">Sample *</th><th colspan="4">Composition</th></tr><tr><th>LLDPE wt %</th><th>PE-g-MA wt %</th><th>MMT wt %</th><th>MMT-Fe wt %</th></tr></thead><tbody><tr><td>LLDPE</td><td>100</td><td></td><td></td><td></td></tr><tr><td>LLDPE.Fu15</td><td>85</td><td>15</td><td></td><td></td></tr><tr><td>LLDPE.Fu15.MMT3</td><td>82</td><td>15</td><td>3</td><td></td></tr><tr><td>LLDPE.Fu5.MMT5</td><td>80</td><td>15</td><td>5</td><td></td></tr><tr><td>LLDPE.Fu15.MMT-Fe3.75</td><td>81.25</td><td>15</td><td></td><td>3.75</td></tr><tr><td>LLDPE.Fu15.MMT-Fe6.25</td><td>78.75</td><td>15</td><td></td><td>6.25</td></tr></tbody></table>

Abbreviations: LLDPE—linear low-density polyethylene; PE-g-MA—polyethylene-grafted maleic anhydride; MMT—montmorillonite. * The number after "Fu" stands for the wt % of PE-g-MA used in the formulation, the number after "MMT" for the wt % of MMT used in the preparation of iron-free reference samples, and the number after "MMT-Fe" for the wt % of MMT-Fe used in active film formulation.

## 2.2. Methods

### 2.2.1. Structural and Morphological Analysis

For TEM analyses, thin sections (100 nm) of the prepared films obtained by ultramicrotomy under cryogenic conditions (−100 °C) were deposited on a holey carbon film supported on a 300-mesh copper TEM grid and analyzed using a Jeol 1200EX2 transmission electron microscope with an acceleration voltage of 100 kV.

WAXS experiments were performed at Laboratoire Charles Coulomb (Univ. Montpellier, France) using an in-house setup consisting of a high brightness and low power Cu Kα X-ray tube coupled with a spheric multilayer optic (GeniX³D from Xenocs) delivering an ultralow divergent beam (0.5 mrad). Scatterless slits were used to produce a clean 0.8 mm diameter beam (35 Mphotons/s) at the sample stage. A transmission configuration was used and the scattered intensity was measured using a Schneider 2D image plate detector placed 0.2 m from the sample. The films were mounted on a custom-made sample holder having Kapton windows to protect the samples from ambient air.

All intensities were corrected by transmission and the empty holder contribution was subtracted. Interplanar distance (d) was calculated using Bragg's equation: $nλ = 2d \sin θ$, where n is a positive integer, λ is the wavelength of the X-ray beam, and θ is the diffraction angle.
### 2.2.2. Thermal Properties

TGA was performed using a TGA Q50 thermobalance (TA Instruments) with a heating ramp of 10 °C·min⁻¹ from room temperature to 1000 °C under flowing air (50 mL·min⁻¹). Approximately 10 mg of the sample was analyzed in an open sapphire crucible. Thermal degradation was characterized by two temperatures: T onset (beginning of the degradation process) and T offset (end of the degradation process).

DSC measurements were carried out using a thermomodulated calorimeter (Q200 modulated DSC, TA Instruments). The sample films (around 7–10 mg) were placed in an aluminum pan and heated from 30 to 200 °C at a heating rate of 10 °C·min⁻¹ and held at this temperature for 5 min to erase the thermal history. The samples were then cooled to room temperature with a cooling rate of 10 °C·min⁻¹. Crystallization and melting thermograms of the samples were obtained. From the DSC thermograms, the melting temperature (T m) and heat of fusion (ΔH m) were obtained. The degree of crystallinity x c was then calculated as follows:

$$x_c = \frac{\Delta H_m + \Delta H_c}{\Delta H_m^0} 100 \quad (2)$$

where ΔH m and ΔH c are the enthalpies of fusion and crystallization, respectively, in J·g⁻¹, and ΔH m 0 is the enthalpy of fusion of fully crystalline polyethylene (293 J·g⁻¹) [22].

### 2.2.3. O₂ and Vapor Barrier Properties

Oxygen permeability was determined according to the standard procedure DIN 53380-5 using an oxygen permeation cell based on an optical luminescence quenching method (OTR-Pst6, PresSens-GmbH) as described elsewhere [23,24]. Oxygen permeability, P O₂, in mol·m⁻¹·s⁻¹·Pa⁻¹, was calculated using the following equation:

$$P_{O_2} = \frac{\dot{P} 1}{A P_{\text{atm}}} \quad (3)$$

where P (mol·s⁻¹) is the slope of the increase of the oxygen partial pressure in the upper test chamber, l (m) is the average thickness of the film, A (m²) is the film area, and P atm (Pa) is the standard atmospheric pressure of oxygen. The measurements were replicated four times. The thickness of the films was determined using a precision micrometer (Mitutoyo).

Water vapor permeability (WVP) (mol·m⁻¹·s⁻¹·Pa⁻¹) measurements were performed by a gravimetric method [25]. Film discs were hermetically sealed (with Teflon seals) in a glass permeation cell containing distilled water and placed at 20 °C inside a desiccator with silica gel. The weight evolution was measured every 24 h for 7 days and plotted against time. The water vapor that transferred through the films was calculated by linear regression:

$$WVP = \frac{w L}{A \Delta_p} \quad (4)$$

where w (mol·s⁻¹) is the slope of the weight loss versus time, L (m) is the film thickness at equilibrium measured at the end of the experiment, A (m²) is the area of exposed film, and Δp the differential vapor pressure of water at 20°C (2337 Pa).

### 2.2.4. Oxygen-Scavenging Capacity

The oxygen-scavenging capacity of the nanocomposite films was obtained by measuring the time–O₂ depletion curve in the headspace of a tightly closed reactor (at least five replications). To do that, 25–30 g of the sample was placed in the reactor, into which 10 mL of pure water was also poured at the bottom to achieve a saturated relative humidity (RH). The reactor was then closed and placed at 20 ± 1 °C in a temperature-controlled oven.
The O₂ partial pressure in the reactor was monitored using a noninvasive optical oxygen sensor (Presens Precision Sensing, Germany) connected to an oxygen measuring device. The limit of detection was 0.03% and the accuracy of 0.4% O₂ at 20.9% of O₂. O₂ partial pressure was recorded until scavenger saturation, which was achieved when no further O₂ consumption was detected (constant O₂ partial pressure in headspace). The mass of absorbed O₂ (in grams per gram of iron) was calculated from the O₂-headspace depletion curve using the perfect gas law. The procedure is well described in detail in the previous work of Kombaya-Touckia-Linin et al. [15].

### 2.2.5. Iron Content by Inductively Coupled Plasma Optical Emission Spectrometry (ICP-OES)

The elemental iron content in MMT-Fe was obtained by ICP-OES (Agilent Technologies) on mineralized solutions of the samples previously obtained by fluoro-nitro perchloric etching according to standard procedures NF X31-147 and ISO 14869-1 [26,27].

### 2.2.6. Mössbauer Spectroscopy

Transmission ⁵⁷Fe Mössbauer spectra were collected with a ⁵⁷Co:Rh source at room temperature. The samples were cooled at 5 K in a helium-flow cryostat (SHI-850 Series, Janis). The velocity driver was operated with a triangular velocity waveform. The velocity scale was calibrated at room temperature with the magnetically split sextet of a high-purity α-Fe foil. The spectra were fitted by least-squares methods to appropriate combinations of Lorentzian profiles representing the quadrupole doublets, magnetic sextets, or octets using the program PC-Mos II [28]. In this way, spectral parameters such as the hyperfine magnetic field (H) quadrupole splitting (Δ), the isomer shift (δ), and the line width at half-maximum (Γ) of the different spectral components were determined for the whole series of spectra. Isomer shifts are given relative to α-Fe at room temperature. To quantify the amount of iron within the films, 16 layers of film samples were mounted together with an iron reference of known iron content on the cryostat.

## 2.3. Modeling of Oxygen Adsorption in the Active Nanocomposite

### 2.3.1. Modeling of O₂ Adsorption for Iron Nanoparticles

An iron-based O₂ scavenger was considered to be composed of iron nanoparticles supported on MMT clay nanoplatelets. The adsorption properties of this scavenger were characterized and formalized using a predictive model, as in the previous work of Kombaya-Touckia-Linin et al. [15]. The oxidation of both Fe⁰ and Fe(OH)₂ species were simultaneously considered, each one following an independent second-order kinetic law. The oxidation model is summarized in Equation (5):

$$ 
\begin{aligned}
\frac{dC_{O_2}}{dt} &= -k_1 n_1 C_{O_2} n_{Fe} - k_2 n_2 C_{O_2} C_{Fe(OH)_2} \\
\frac{dC_{Fe}}{dt} &= -k_1 C_{O_2} C_{Fe} \\
\frac{dC_{Fe(OH)_2}}{dt} &= -k_2 C_{O_2} C_{Fe(OH)_2}
\end{aligned}
\qquad (5) $$

where C₀, CFe, and CFe(OH)₂ are the concentrations (mol·m⁻³) of O₂, Fe, and Fe(OH)₂; n₁ and n₂ are the apparent stoichiometric coefficients for oxidation of Fe and Fe(OH)₂ by O₂; and k₁ and k₂ are the kinetic coefficients (m³·s⁻¹·mol⁻¹) for oxidation of Fe and Fe(OH)₂, respectively.

### 2.3.2. Modeling of O₂ Adsorption for an Active Nanocomposite

The phenomena occurring at the nanocomposite level can be described by a reaction-diffusion system. Fick's law of diffusion was applied to describe the diffusion of O₂. The reaction part was modeled according to Equation (5).

MMT nanoplatelets and iron nanoparticles were assumed to be immobile in the polymer matrix, and the nanocomposite was considered as a homogenous material, that is, with a single apparent
$O_2$ diffusivity $D_{O_2}(\text{m}^2\cdot\text{s}^{-1})$ and the iron species uniformly distributed in the material. Finally, the mathematical model for a plane film geometry reduced to the one-dimensional reaction-diffusion system is given in Equation (6), for $x \in ]-L/2, L/2[$, where L is the thickness of the film:

$$ 
\begin{aligned}
\frac{\partial C_{O_2}(t,x)}{\partial t} &= D_{O_2} \frac{\partial^2 C_{O_2}(t,x)}{\partial x^2} - k_1 n_1 C_{O_2}(t,x) C_{Fe}(t,x) - k_2 n_2 C_{O_2}(t,x) C_{Fe(OH)_2}(t,x) \\
\frac{\partial C_{Fe}(t,x)}{\partial t} &= -k_1 C_{O_2}(t,x) C_{Fe}(t,x) \\
\frac{\partial C_{Fe(OH)_2}(t,x)}{\partial t} &= -k_2 C_{O_2}(t,x) C_{Fe(OH)_2}(t,x)
\end{aligned}
\qquad (6) $$

with the initial (Equation (7)) and boundary (Equation (8)) conditions given hereafter.

Initial condition: at initial time $t_0$, uniform concentrations $C_{Fe}(t_0)$, and $C_{Fe(OH)_2}(t_0)$ were considered and computed according to Equation (7):

$$ 
\begin{aligned}
C_{Fe}(t_0) &= \frac{\rho_f x_{Fe}^f x_{Fe}}{M_{Fe}} \\
C_{Fe(OH)_2}(t_0) &= \frac{\rho_f x_{Fe}^f x_{Fe(OH)_2}}{M_{Fe}}
\end{aligned}
\qquad (7) $$

where $x_{Fe}^f$ (kilograms of iron per kilogram of sample) is the mass fraction of iron in the sample, $x_{Fe}$ and $x_{Fe(OH)_2}$ (kilograms per kilogram of total iron) are the mass fraction of zerovalent and divalent iron, respectively, $M_{Fe}$ (kg·mol⁻¹) is the molar mass of iron, and $\rho_f$ (kg·m⁻³) is the apparent density of the nanocomposite. A null oxygen concentration, $C_{O_2}(t_0) = 0$, was assumed.

Boundary condition: the external conditions at the upper and lower surfaces of the active nanocomposite film surrounded by headspace were equal. For instance, the flow that went in the film at the upper surfaces $\Gamma_{L/2}$ is given by

$$ \varphi_{L/2} = k A (C_{O_2, HS} - C_{O_2, HS, L/2}) \qquad (8) $$

where $k$ (m·s⁻¹) is the external mass transfer coefficient; $A$ (m²) is the surface area of the headspace-film interface, which was assumed to be identical for the upper and lower surfaces; and $C_{O_2, HS}$ (mol·m⁻³) and $C_{O_2, HS, i}$ (mol·m⁻³) are, respectively, the O₂ concentration in the headspace and at the vicinity of the composite surface. Converting concentration into partial pressure thanks to the ideal gas law and Henry's law (more details in Appendix A) finally leads to the two following boundary conditions:

$$ 
\begin{aligned}
D_{O_2} \frac{\partial C_{O_2}(t,x)}{\partial x} &= \frac{\varphi_{L/2}}{A} = \frac{k}{RT} \left( p_{O_2, HS} - \frac{C_{O_2}(t,x)}{k_H} \right) \quad \text{at } x = \frac{L}{2} \text{ and } \forall t \ge 0 \\
D_{O_2} \frac{\partial C_{O_2}(t,x)}{\partial x} &= \frac{\varphi_{-L/2}}{A} = \frac{k}{RT} \left( p_{O_2, HS} - \frac{C_{O_2}(t,x)}{k_H} \right) \quad \text{at } x = -\frac{L}{2} \text{ and } \forall t \ge 0.
\end{aligned}
\qquad (9) $$

The external mass transfer coefficient $k$ can be expressed through the dimensionless Biot number $Bi$, as presented in Equation (10):

$$ k = \frac{2 Bi D_{O_2}}{L}. \qquad (10) $$

Mass balance of $O_2$ in the headspace: the composite film was assumed to be inside an airtight container, so that the headspace volume $V_{HS}$ (m³) was constant and the gas flow through the container was negligible. Consequently, the variation of the headspace $O_2$ partial pressure could be computed by integrating the flows $\varphi_{L/2}$ and $\varphi_{-L/2}$ on the film surface area and again using the ideal gas law, which leads to Equation (11):

$$ \frac{\partial p_{O_2, HS}}{\partial t} = \frac{RT}{V_{HS}} \left( \iint_{\Gamma_{L/2}} \varphi_{L/2} dS + \iint_{\Gamma_{-L/2}} \varphi_{-L/2} dS \right). \qquad (11) $$
As the flows are homogeneous on the surface, Equation (11) is rewritten as

$$ \frac{\partial p_{0, \text{HS}}}{\partial t} = k \frac{A}{V_{\text{HS}}} \left( 2 p_{0, \text{HS}} - \frac{C_{\text{O}_2}(t, x = L/2)}{k_{\text{H}}} - \frac{C_{\text{O}_2}(t, x = L/2)}{k_{\text{H}}} \right). \quad (12) $$

### 2.3.3. Numerical Simulations

Numerical simulations were performed with a high Biot number ($Bi = 10^4$), characteristic of this diffusion-limited case. The partial differential equation system (Equations (2), (3), (8), and (11)) was transformed into an ordinary differential equation (ODE) system by spatial discretization with a second-order central difference method and a mesh of 100 nodes. The resulting ODE system was numerically solved using MATLAB (MathWorks).

### 2.4. Statistics

Statistical tests were performed using R software for statistical computing (R, 2014). One-way ANOVA tests were performed to verify if the variances of the datasets were statistically different or not. Comparisons between the compositions were performed by pairwise comparisons using Tukey's test. Different letters were used to denote significant differences between datasets (level of significance $\alpha = 0.05$, unless stated).

## 3. Results

### 3.1. Structural and Morphological Characterization of the MMT-Fe and MMT-Fe/LLDPE Films

The as-synthesized MMT-Fe powder contained about 0.22 ± 0.01 g of Fe per gram of powder, as determined by ICP-OES. This powder was embedded, after a short storage period in anoxic conditions (15 days max.), in a matrix of LLDPE containing PE-g-MA (15 wt %) as an additive in quantities previously tested and validated [19,20,29].

TEM images of the MMT-Fe particles are shown in Figure 1. Particles of sizes between 6 and 20 nm forming aggregates of 57 ± 17 nm [15] were observed on the surface of montmorillonite.

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Transmission electron microscopy (TEM) images of MMT-Fe at different magnifications, showing the dispersion of the aggregates (A) and their microstructure (B).</figcaption></figure>
TEM analyses were also performed to evaluate the dispersion of MMT and MMT-Fe in the LLDPE matrix. The micrographs shown in Figure 2 show a quite good dispersion of the MMT in all the composites, even though a few small agglomerates of clay platelets were observed in LLDPE/MMT and LLDPE/MMT-Fe. A preferred orientation of the nanoclay platelets was observed for all nanocomposites, with the exception of the nanocomposite with high MMT loading. This effect may have been related to the shear forces induced during extrusion and calendaring, which aligned the nanoplatelets in the direction of the flow. The presence of iron nanoparticles was not clearly visible from TEM pictures except at a high MMT-Fe content, where some iron aggregates are visible in the TEM picture (Figure 2B).

<figure><img src="image_6.png" /><figcaption>**Figure 2.** TEM micrographs of (A) LLDPE.Fu15.MMT5 film, (B) LLDPE.Fu15.MMT-Fe6.25 film, (C) LLDPE.Fu15.MMT3 film, and (D) LLDPE.Fu15.MMT-Fe3.75 film.</figcaption></figure>

The WAXS results, shown in Figure 3, confirmed that the nanoclays were partially exfoliated in the polymer independent of the formulation and of the presence of the active compounds. The WAXS pattern of MMT showed an intense diffraction peak, which corresponded to a spacing distance ($d_{001}$) equal to 3.34 nm [30]. This peak was no longer observed on patterns of MMT-Fe and the active nanocomposites, indicating a virtually complete delamination of the clay in these materials.

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Wide-angle X-ray scattering (WAXS) patterns for nanoclay powder MMT, MMT-Fe, LLDPE, and LLDPE matrix compatibilized with PE-g-MA and active nanocomposite films (the curves are vertically offset).</figcaption></figure>
## 3.2. Effect of MMT-Fe Content on the Thermal Stability of LLDPE Nanocomposite Films

The results of DSC and TGA, which were used to study the thermal stability of MMT-Fe and the active films, are shown in Figure 4. The corresponding data (heat of fusion, melting point, and crystallinity determined by DSC and temperature of degradation determined by TGA) are compiled in Table 2. For MMT and MMT-Fe, a first mass loss was observed around 84 °C, corresponding to the desorption of volatile species (ethanol), while a second mass loss between 200 and 500 °C, characteristic of Cloisite 20A, corresponded to the departure of organic substances (alkyl ammonium). No significant difference was observed by TGA between the different polymer nanocomposites: the presence of the PE-g-MA, MMT, and/or MMT-Fe had no significant influence on the thermal degradation of the polymer matrix, which decomposed at 450 °C in all of studied materials.

<figure><img src="image_8.png" /><figcaption>**Figure 4.** (a) Thermogravimetric (TG) curves in flowing air of LLDPE, MMT, MMT-Fe, and nanocomposites with different MMT and MMT-Fe contents. (b) Differential scanning calorimetry (DSC) thermograms of LLDPE and LLDPE nanocomposites.</figcaption></figure>

<table><caption>**Table 2.** Thermal properties of LLDPE and LLDPE nanocomposites.</caption><thead><tr><th>Sample</th><th>Heat of Fusion (J/g)</th><th>Tm (°C)</th><th>Crystallinity</th><th>Temp. Onset Degradation (°C)</th><th>Temp. Endset Degradation (°C)</th></tr></thead><tbody><tr><td>LLDPE</td><td>104 ± 10</td><td>122.1 ± 0.5</td><td>35.8 ± 3.6</td><td>416 ± 9</td><td>475 ± 4</td></tr><tr><td>LLDPE.Fu15</td><td>104 ± 8</td><td>121.3 ± 1.2</td><td>35.7 ± 3.1</td><td>416 ± 2</td><td>457 ± 2</td></tr><tr><td>LLDPE.Fu15.MMT3</td><td>104 ± 5</td><td>121.9 ± 0.7</td><td>34.8 ± 0.9</td><td>416 ± 22</td><td>478 ± 19</td></tr><tr><td>LLDPE.Fu15.MMT5</td><td>99 ± 10</td><td>122.2 ± 1.0</td><td>31.3 ± 0.6</td><td>429 ± 9</td><td>480 ± 4</td></tr><tr><td>LLDPE.Fu15.MMT-Fe3.75</td><td>107 ± 2</td><td>121.6 ± 0.4</td><td>36.3 ± 0.7</td><td>419 ± 5</td><td>468 ± 3</td></tr><tr><td>LLDPE.Fu15.MMT-Fe6.25</td><td>102 ± 4</td><td>122.3 ± 1.2</td><td>34.5 ± 1.2</td><td>418 ± 7</td><td>466 ± 2</td></tr></tbody></table>

The data reported here were not significantly different (ANOVA test, *p*-value ≤ 0.05). Raw data corresponding to this figure are uploaded at https://doi.org/10.5281/zenodo.3253103. See supplementary material.

DSC was used to study the influence of MMT and MMT-Fe on melting temperature (Tm), specific melting enthalpy (ΔHm), and the degree of crystallinity of LLDPE (Table 2). No significant differences in the melting temperatures were observed for LLDPE, LLDPE.Fu15, and its active nanocomposite samples (LLDPE.Fu.MMT3, LLDPE.Fu15.MMT-Fe3.75, and LLDPE.Fu15.MMT-Fe6.25). Similar results were obtained for the heat of fusion and crystallinity, except for LLDPE.Fu15.MMT5, which displayed significantly lower crystallinity. We can confirm from these results that the MMT and MMT-Fe powder did not act as nucleating agents for the LLDPE pure matrix.

## 3.3. Mössbauer Spectroscopy Study of the Active Nanocomposite

The ⁵⁷Fe Mössbauer spectra at 5 K of the MMT-Fe powder, the fresh film obtained directly after extrusion (LLDPE.Fu15.MMT-Fe6.25), and the same sample exposed to ambient air for four months are shown in Figure 5. The spectra were fitted to several spectral components as described in
Kombaya-Touckia-Linin et al. [15]. The same analysis was performed on similar samples based on LLDPE.Fu15.MMT-Fe3.75.

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Examples of ⁵⁷Fe Mössbauer spectra at 5 K for (A) MMT-Fe, (B) fresh LLDPE.Fu15%.MMT-Fe6.25%, and (C) LLDPE.Fu15%.MMT-Fe6.25% exposed to ambient air for four months.</figcaption></figure>

In the spectrum of MMT-Fe, four different iron environments can be detected: magnetically ordered Fe²⁺ (blue octet), paramagnetic Fe²⁺ (cyan doublet), Fe³⁺ oxide (red sextet), and the blue sextet characterized by an isomer shift of 0.27 mm/s and a hyperfine field of 30.7 T indicating the presence of Fe₁⁻ₓBₓ, as already observed previously [15,31]. These parameters, significantly different from those of α-Fe metal, confirm the formation of amorphous iron boride from the reduction of iron by the borohydride anion [32,33]. Fe₁⁻ₓBₓ was also detected on the fresh films LLDPE.Fu15.MMT-Fe6.25 and LLDPE.Fu15.MMT-Fe3.75. In the films, however, this sextet was no longer detected and was replaced by a large paramagnetic singlet, probably resulting from the film processing treatment, which could have induced structural disorder in Fe₁⁻ₓBₓ. None of the observed Fe₁⁻ₓBₓ components were stable in air and were rapidly converted to iron oxides upon air exposure.

For the fresh LLDPE.Fu15.MMT-Fe6.25 film, the two dominant components were typical of Fe³⁺ in iron oxides or oxyhydroxides (Table 3), summing up to 75% of the total resonance area. This confirms that MMT-Fe was already partially oxidized during the film processing step. The presence of two large sextets indicated the presence of a distribution of hyperfine magnetic fields, typical of relatively ill-defined oxide or oxyhydroxide species. For LLDPE.Fu15.MMT-Fe6.25 exposed to air for four months, an increase of the oxide resonance area to 86% was observed, confirming the gradual oxidation of the sample with time.
<table><caption>**Table 3.** $^{57}$Fe Mössbauer parameters at 5 K for samples MMT-Fe, fresh LLDPE.Fu15.MMT-Fe6.25, and aged LLDPE.Fu15.MMT-Fe6.25.</caption><thead><tr><th>Sample</th><th>Component</th><th>δ (mm/s)</th><th>Δ (mm/s)</th><th>Γ (mm/s)</th><th>B<sub>HF</sub> (T)</th><th>Area (%)</th></tr></thead><tbody><tr><td rowspan="4">MMT-Fe</td><td>Fe<sub>1-x</sub>B<sub>x</sub></td><td>0.08</td><td>-0.03</td><td>0.47</td><td>30.8</td><td>53</td></tr><tr><td>Fe<sup>2+</sup> para</td><td>1.25</td><td>3.03</td><td>0.66</td><td>-</td><td>6</td></tr><tr><td>Fe<sup>2+</sup> magn</td><td>1.21</td><td>-2.78</td><td>1.72</td><td>15.0</td><td>6</td></tr><tr><td>Fe<sup>3+</sup> oxide</td><td>0.37</td><td>-0.12</td><td>0.60</td><td>49.0</td><td>35</td></tr><tr><td rowspan="6">Fresh LLDPE.Fu15.MMT-Fe6.25</td><td>Fe<sub>1-x</sub>B<sub>x</sub> type1</td><td>0.00</td><td>0.00</td><td>0.60</td><td>25.00</td><td>3</td></tr><tr><td>Fe<sub>1-x</sub>B<sub>x</sub> type2</td><td>0.24</td><td>0.00</td><td>0.86</td><td>-</td><td>8</td></tr><tr><td>Fe<sup>2+</sup> para</td><td>1.28</td><td>3.09</td><td>0.64</td><td>-</td><td>4</td></tr><tr><td>Fe<sup>2+</sup> magn</td><td>1.21</td><td>-2.07</td><td>1.14</td><td>19.45</td><td>10</td></tr><tr><td>Fe<sup>3+</sup> oxide type1</td><td>0.34</td><td>0.00</td><td>0.85</td><td>51.62</td><td>75</td></tr><tr><td>Fe<sup>3+</sup> oxide type 2</td><td>0.41</td><td>0.00</td><td>2.5</td><td>47.17</td><td></td></tr><tr><td rowspan="6">LLDPE.Fu15.MMT-6.25 after four months of oxidation</td><td>Fe<sub>1-x</sub>B<sub>x</sub> type1</td><td>0.00</td><td>0.00</td><td>0.60</td><td>24.55</td><td>3</td></tr><tr><td>Fe<sub>1-x</sub>B<sub>x</sub> type2</td><td>0.30</td><td>0.00</td><td>1.11</td><td>-</td><td>3</td></tr><tr><td>Fe<sup>2+</sup> para</td><td>1.30</td><td>3.04</td><td>0.56</td><td>-</td><td>3</td></tr><tr><td>Fe<sup>2+</sup> magn</td><td>1.21</td><td>-2</td><td>1.81</td><td>19.52</td><td>5</td></tr><tr><td>Fe<sup>3+</sup> oxide type1</td><td>0.37</td><td>0.00</td><td>1.27</td><td>50.76</td><td>86</td></tr><tr><td>Fe<sup>3+</sup> oxide type 2</td><td>0.41</td><td>0.00</td><td>2.5</td><td>46.16</td><td></td></tr></tbody></table>

### 3.4. Oxygen-Scavenging Capabilities for MMT-Fe and Active Nanocomposites

The oxygen adsorption capacity of the synthesized MMT-Fe powder and the films was measured in static conditions in a saturated water vapor atmosphere (100% relative humidity), as the oxidation of iron by oxygen requires the presence of water molecules (0.5 g of water·g⁻¹ of iron) [34,35]. At least five replications for each sample were carried out. Figure 6a shows the headspace–O₂ depletion curve and the derived O₂ adsorption curve (calculated by using the perfect gas law approximation) for the active film LLDPE.Fu15.MMT-Fe6.25.

The adsorption capacity of the films (LLDPE.Fu15.MMT-Fe6.25 and LLDPE.Fu15.MMT-Fe3.75) was compared with that of MMT-Fe (Figure 6b). The film adsorption capacity remained low compared with that of the powder used alone. In fact, oxygen-scavenging capacities of 31 ± 2 and 55 ± 9 mg O₂·g⁻¹ iron were obtained for LLDPE.Fu15.MMT-Fe3.75 and LLDPE.Fu15.MMT-Fe6.25, respectively, while the adsorption capacity of MMT-Fe was on average 122 mg O₂·g⁻¹ iron. This result confirms the substantial oxidation of the iron in the pristine film samples after the thermomechanical process. More precisely, iron in fresh LLDPE.Fu15.MMT-Fe6.25 and LLDPE.Fu15.MMT-Fe3.75 was approximately twice and three times more oxidized, respectively, than in MMT-Fe, which agrees with Mössbauer spectroscopy analysis.

Moreover, the iron oxidation kinetics was slower in LLDPE.Fu15.MMT-Fe3.75 than in LLDPE.Fu15.MMT-Fe6.25 and did not reach equilibrium at the end of the experiment (two weeks). This might have been due to the larger thickness of LLDPE.Fu15.MMT-Fe3.75, owing to extrusion constraints at lab scale, which limited O₂ diffusion to the core of the film and led to a lower apparent O₂ adsorption capacity.
<figure><img src="image_10.png" /><figcaption>**Figure 6.** (a) Example of experimental O₂-headspace depletion curves and corresponding O₂ adsorption curves for LLDPE.Fu15.MMT-Fe6.25 and (b) comparison of oxygen adsorption of dried MMT-Fe powder, LLDPE.Fu15.MMT-Fe3.75, and LLDPE.Fu15.MMT-Fe6.25 films.</figcaption></figure>

The adsorption capacity of the active films prepared here was in the same order of magnitude of those of previously published systems (Table 4). A calculation of the benefit of using these active films for food packaging applications is provided in Appendix B.

<table><caption>**Table 4.** Comparison with results reported in the literature.</caption><thead><tr><th>Film</th><th>O₂ Adsorption Capacity<br>mg O₂/g of Absorber</th><th>Ref.</th></tr></thead><tbody><tr><td>PE.Fu15.MMT-Fe6.25 film</td><td>55 ± 2</td><td>This work</td></tr><tr><td>PE.Fu15.MMT-Fe3.75 film</td><td>31 ± 2</td><td>This work</td></tr><tr><td>Polymer with iron powder and additives ("SHELFPLUS"); masterbatch</td><td>26–48</td><td>[36]</td></tr><tr><td>O₂Block1 (NanoBioMatters S.L., Paterna, Spain); masterbatch</td><td>>10–25</td><td>[37]</td></tr><tr><td>Cyclo-olefin bonded to a silicate backbone "ORMOCER1"; lacquer</td><td>90</td><td>[38]</td></tr><tr><td>Nanoscale iron in silicon matrix</td><td>>32 ± 80</td><td>[7]</td></tr><tr><td>Polyolefin films containing 10 wt % of iron modified kaolinite</td><td>5.7</td><td>[9]</td></tr><tr><td>Nanoscale iron powder</td><td>>68 ± 2</td><td>[7]</td></tr><tr><td>Nanoscale iron-MMT powder</td><td>122</td><td>[15]</td></tr><tr><td>Stoichiometric capacity of iron</td><td>394</td><td>[15]</td></tr></tbody></table>

### 3.5. Oxygen and Water Vapor Barrier Properties of LLDPE/MMT-Fe Active Film Composite

The oxygen permeability of the composite films at 23 °C and 0% RH conditions is shown in Figure 7a. The addition of PE-g-MA to the LLDPE reduced the permeability by about 15%. The addition of 3% or 5% MMT, producing an intercalated nanocomposite structure, further reduced the permeability by 20% compared with LLDPE. No significant difference in oxygen permeability was observed when MMT-Fe was added instead of MMT. Even if the active films displayed significant O₂-scavenging properties (see previous paragraph), this scavenging effect did not impact the O₂ permeability property in the conditions of measurement of the present work.
<figure><img src="image_11.png" /><figcaption>**Figure 7.** (a) Oxygen permeability and (b) water vapor permeability of (a) LLDPE film and LLDPE with compatibilizer and active films containing MMT-Fe with different contents (3.75% or 6.25%). The data are reported as mean ± standard deviation of four replicates. Bars with different letters represent significant difference (ANOVA test, *p*-value ≤ 0.05). Raw data corresponding to this figure has been uploaded at https://doi.org/10.5281/zenodo.3253103.</figcaption></figure>

Also, the water vapor permeability was checked for all the films. No significant differences were detected among the different films (Figure 7b), also considering the relatively large experimental error bars. Further experiments are thus needed to confirm this result.

### 3.6. Modeling the Oxygen-Scavenging Activity

The reaction model developed in the previous work of Kombaya-Touckia-Linin et al. was here combined with a model describing O₂ solubilization and diffusion in the polymer matrix. The parameters used for modeling O₂ adsorption by LLDPE.Fu15.MMT-Fe6.25, calculated according to Equations (6), (7), (9), and (12), are listed in Table 5, while the comparison between experimental and simulated O₂ depletion is shown in Figure 8. Good agreement was found between the experimental and modeled curves. It is worth noting that no additional parameters were required. Even with the simplified assumptions, the reaction-diffusion model was able to reproduce the observed data of O₂ adsorption in an active nanocomposite, validating the O₂ adsorption model previously developed for MMT-Fe.

The O₂-scavenging properties of the produced films were strongly impacted by the extrusion process, which induced substantial iron oxidation. As a comparison, a simulation was performed assuming that no iron oxidation occurs during extrusion (Figure 9). This “theoretical” case obviously produces higher O₂ adsorption capacities, twice as much as actually measured. The production of such highly active O₂-scavenging film will require a thorough optimization of both the synthesis of MMT-Fe active powder and the extrusion process leading to the nanocomposite film.
<table><caption>**Table 5.** Parameter descriptions and values for numerical simulation of O₂ adsorption for the PE.Fu15.MMT-Fe6.25 sample.</caption><thead><tr><th>Symbol</th><th>Definition</th><th>Unit</th><th>Value</th><th>Source</th></tr></thead><tbody><tr><td>L</td><td>Thickness of the nanocomposite sample</td><td>m</td><td>296 × 10<sup>-6</sup></td><td>This work</td></tr><tr><td>A</td><td>Surface area of a single side of the nanocomposite film sample</td><td>m<sup>2</sup></td><td>0.1376</td><td>This work</td></tr><tr><td>V<sub>HS</sub></td><td>Volume of the headspace</td><td>m<sup>3</sup></td><td>458 × 10<sup>-6</sup></td><td>This work</td></tr><tr><td>x<sub>F<sub>e</sub></sub><sup>f</sup></td><td>Mass fraction of iron in the sample</td><td>kg(iron)·kg(sample)<sup>-1</sup></td><td>0.22</td><td>This work</td></tr><tr><td>x<sub>F<sub>e</sub></sub></td><td>Mass fraction of zerovalent iron</td><td>kg·kg(iron)<sup>-1</sup></td><td>0.11</td><td>This work</td></tr><tr><td>x<sub>Fe(OH)<sub>2</sub></sub></td><td>Mass fraction of bivalent iron</td><td>kg·kg(iron)<sup>-1</sup></td><td>0.14</td><td>This work</td></tr><tr><td>M<sub>F<sub>e</sub></sub></td><td>Molar mass of iron</td><td>kg·mol<sup>-1</sup></td><td>55.8 × 10<sup>-3</sup></td><td></td></tr><tr><td>ρ<sub>f</sub></td><td>Apparent density of the nanocomposite, taken as the corresponding value for the pure polymer (LLDPE)</td><td>kg·m<sup>-3</sup></td><td>918</td><td>[15]</td></tr><tr><td>k<sub>1</sub></td><td>Kinetic coefficient of Fe oxidation</td><td>m<sup>3</sup>·s<sup>-1</sup>·mol<sup>-1</sup></td><td>2.4 × 10<sup>-6</sup></td><td>[15]</td></tr><tr><td>k<sub>2</sub></td><td>Kinetic coefficient of Fe(OH)<sub>2</sub> oxidation</td><td>m<sup>3</sup>·s<sup>-1</sup>·mol<sup>-1</sup></td><td>3.2 × 10<sup>-5</sup></td><td>[15]</td></tr><tr><td>n<sub>1</sub></td><td>Apparent stoichiometric coefficient of Fe oxidation</td><td>mol·mol<sup>-1</sup></td><td>0.19</td><td>[15]</td></tr><tr><td>n<sub>2</sub></td><td>Apparent stoichiometric coefficient of Fe(OH)<sub>2</sub> oxidation</td><td>mol·mol<sup>-1</sup></td><td>0.40</td><td>[15]</td></tr><tr><td>D<sub>O<sub>2</sub></sub></td><td>Apparent O<sub>2</sub> diffusivity in the nanocomposite, taken as the corresponding value for the pure polymer (LLDPE)</td><td>m<sup>2</sup>·s<sup>-1</sup></td><td>1.68 × 10<sup>-11</sup></td><td>[39]</td></tr><tr><td>k<sub>H</sub></td><td>O<sub>2</sub> solubility in the nanocomposite material</td><td>molPa·m<sup>-3</sup></td><td>3.54 × 10<sup>-6</sup></td><td>[40]</td></tr><tr><td>Bi</td><td>Biot number</td><td>dimensionless</td><td>1 × 10<sup>5</sup></td><td>This work</td></tr><tr><td>R</td><td>Ideal gas constant</td><td>J·mol<sup>-1</sup>·K<sup>-1</sup></td><td>8.314</td><td></td></tr><tr><td>T</td><td>Temperature of the experiment</td><td>K</td><td>293.15</td><td>This work</td></tr></tbody></table>

<figure><img src="image_12.png" /><figcaption>**Figure 8.** Comparison between simulated (red line) and experimental (grey dots) headspace–O₂ depletion curves.</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>**Figure 9.** Simulated headspace O₂ partial pressure: actual case, i.e., with oxidation during extrusion (red line), and theoretical case, i.e., with no oxidation during extrusion (black dashed line).</figcaption></figure>

**4. Conclusions**

The oxygen-scavenging capacity of an MMT-Fe composite, designed as the active species for monolayer film applications, as well as its LLDPE nanocomposites was evaluated. The MMT-Fe contained Fe₁⁻ₓBₓ nanoparticles with an average aggregate diameter of ~57 nm deposited on the surface of MMT clay platelets. Active nanocomposites containing 6.25 or 3.75 wt % of MMT-Fe were manufactured by extrusion with LLDPE, showing adsorption capacities two times lower than that of pristine MMT-Fe. This decrease in adsorption capacity was caused by the oxidation of the Fe₁⁻ₓBₓ nanoparticles during the extrusion process, as assessed by Mössbauer spectroscopy. The active films containing 3.75 and 6.25 wt % were able to uptake between 31 ± 2 and 55 ± 9 mg O₂·g⁻¹ iron, respectively, while MMT-Fe had an adsorption capacity of 122 mg O₂·g⁻¹ iron at 20 °C and 100% RH. The reaction–diffusion model was able to reproduce the observed data of O₂ adsorption. This study suggests that there is significant potential for the use of this new oxygen-scavenging composite, MMT-Fe, which, if included in packaging under the correct anoxic conditions, might help to extend the shelf life of oxygen-sensitive food products.

**Supplementary Materials:** Raw data of Table 2 and Figure 7 are available online at http://www.mdpi.com/2079-4991/9/9/1193/s1; and https://doi.org/10.5281/zenodo.3253103.

**Author Contributions:** E.-M.K.-T.-L. performed the experiments and contributed to the writing of the draft; S.G. did the modeling and parameter optimization and contributed to the writing and proofreading of the manuscript; M.S. and L.S. performed the Mössbauer experiments and their analysis and contributed to the writing and proofreading of the manuscript; N.G. contributed to the work supervising and manuscript proofreading; and V.G. contributed to the writing and proofreading of the manuscript as well as the general supervision of the work.

**Funding:** This work has been carried out in the framework of the Labex SERENADE (ANR-11-LABX-0064) funded by the Investissements d'Avenir French Government program managed by the French National Research Agency (ANR) and in the framework of the French program Mutations alimentaires 2019, cofinanced by INRA and CNRS.

**Conflicts of Interest:** The authors declare no conflict of interest.

**Appendix A**

More details about the initial and boundary conditions of the diffusion–reaction model developed in the present work are provided hereafter.

**Initial and boundary conditions.** At initial time $t_0$, uniform concentrations $C_{\text{Fe}}(t_0)$ and $C_{\text{Fe(OH)}_2}(t_0)$ were considered and computed according to Equation (7), renamed hereafter as Equation (A1):
$$ 
\begin{aligned}
C_{\text{Fe}}(t_0) &= \frac{\rho_f x_{\text{Fe}}^f x_{\text{Fe}}}{M_{\text{Fe}}} \\
C_{\text{Fe(OH)}_2}(t_0) &= \frac{\rho_f x_{\text{Fe}}^f x_{\text{Fe(OH)}_2}}{M_{\text{Fe}}}
\end{aligned}
\qquad (\text{A1}) $$

where $x_{Fe}^f$ (kilograms of iron per kilogram of sample) is the mass fraction of iron in the sample, $x_{Fe}$ and $x_{Fe(OH)2}$ (kilograms per kilogram of total iron) are the mass fraction of zerovalent and divalent iron, respectively, $M_{Fe}$ (kg·mol⁻¹) is the molar mass of iron, and $\rho_f$ (kg·m⁻³) is the apparent density of the nanocomposite. A null oxygen concentration, $C_{O_2}(t_0) = 0$, was assumed.

Focusing on the $O_2$ adsorption by an active nanocomposite film, the material was surrounded by a unique headspace, and the external conditions at upper and lower surfaces were equal. The $O_2$ molar flow $\varphi_{L/2}$ (mol·s⁻¹) at the upper headspace-film interface $\Gamma_{L/2}$ is given in Equation (A2):

$$ \varphi_{L/2} = k A (C_{O_2, \text{HS}} - C_{O_2, \text{HS,L/2}}) \qquad (\text{A2}) $$

where $k$ (m·s⁻¹) is the external mass transfer coefficient at the upper and lower surfaces $\Gamma_{L/2}$ and $\Gamma_{-L/2}$; $A$ (m²) is the surface area of the headspace-film interface, which was assumed to be identical for the upper and lower surfaces; and $C_{O_2, \text{HS}}$ (mol·m⁻³) and $C_{O_2, \text{HS,i}}$ (mol·m⁻³) are, respectively, the $O_2$ concentration in the headspace and at the vicinity of the composite surface.

From the ideal gas law, Equation (A2) can be rewritten as

$$ \varphi_{L/2} = k \frac{A}{RT} (P_{O_2, \text{HS}} - P_{O_2, \text{HS,L/2}}) \qquad (\text{A3}) $$

where $T$ (K) is the temperature, $R$ is the ideal gas constant, $P_{O_2, \text{HS}}$ (Pa) and $P_{O_2, \text{HS,i}}$ (Pa) are the $O_2$ partial pressure in the headspace and at the vicinity of the composite surface, respectively.

According to Henry's law, assuming that the $O_2$ is at a thermodynamic equilibrium between the composite and the headspace, the concentration $C_{O_2}(t, x = L/2)$ (mol·m⁻³) of dissolved $O_2$ at the composite surface is proportional to $P_{O_2, \text{HS,L/2}}$, the $O_2$ partial pressure in the vicinity of the composite surface. The $O_2$ flow $\varphi_{L/2}$ at the upper headspace-film interface $\Gamma_{L/2}$ becomes

$$ \varphi_{L/2} = k \frac{A}{RT} \left( P_{O_2, \text{HS}} - \frac{C_{O_2}(t, x = L/2)}{k_H} \right) \qquad (\text{A4}) $$

where $k_H$ (mol·Pa⁻¹·m⁻³) is the $O_2$ solubility in the composite material. In the same way, the $O_2$ flow $\varphi_{-L/2}$ at the lower headspace-film interface $\Gamma_{-L/2}$ is described by Equation (A5):

$$ \varphi_{-L/2} = k \frac{A}{RT} \left( P_{O_2, \text{HS}} - \frac{C_{O_2}(t, x = -L/2)}{k_H} \right). \qquad (\text{A5}) $$

Finally, the boundary conditions can be written as

$$ 
\begin{aligned}
D_{O_2} \frac{\partial C_{O_2}(t,x)}{\partial x} &= \frac{\varphi_{\frac{L}{2}}}{A} = \frac{k}{RT} \left( P_{O_2, \text{HS}} - \frac{C_{O_2}(t,x)}{k_H} \right) \quad \text{at } x = \frac{L}{2} \text{ and } \forall t \ge 0 \\
D_{O_2} \frac{\partial C_{O_2}(t,x)}{\partial x} &= \frac{\varphi_{-\frac{L}{2}}}{A} = \frac{k}{RT} \left( P_{O_2, \text{HS}} - \frac{C_{O_2}(t,x)}{k_H} \right) \quad \text{at } x = -\frac{L}{2} \text{ and } \forall t \ge 0
\end{aligned}
\qquad (\text{A6}) $$

## Appendix B

In order to verify the possible advantage of using active films prepared in the present work as $O_2$ scavengers for food packaging applications, their theoretical adsorption capacity in a real situation was estimated. We considered a food/packaging system stored at ambient temperature, with a packaging headspace of 100 mL. We considered the adsorption of $O_2$ by the food as well as its permeation through the packaging to be negligible. In that condition, including about 1 g of MMT-Fe powder in the
headspace polymer would permit a decrease of the volume of O₂ from 21 to 1.5 mL. For the same amount of MMT-Fe but embedded in the bulk polymer (active film), considering its level of oxidation after processing and that the film, in the ideal case, would be protected from the outside atmosphere by a barrier layer before use, the decrease of O₂ in the headspace would be about half of the initial content (10.5 mL) (Figure A1). In such conditions, the active packaging would be much less efficient than MMT-Fe used alone in the headspace and insufficient for acceptable removal of O₂ from the packaging. However, such a film could have a positive impact when used in combination with other methods, such as modified-atmosphere packaging (e.g., N₂ flush) for the adsorption of residual amounts of O₂ in the packaging atmosphere.

<figure><img src="image_14.png" /><figcaption>**Figure A1.** Schematic representation of adsorption capacity compatibility with food packaging.</figcaption></figure>

References

1. Singh, N.; Hui, D.; Singh, R.; Ahuja, I.P.S.; Feo, L.; Fraternali, F. Recycling of plastic solid waste: A state of art review and future applications. *Compos. Part B Eng.* **2017**, *115*, 409–422. [CrossRef]

2. Ragaert, K.; Delva, L.; Van Geem, K. Mechanical and chemical recycling of solid plastic waste. *Waste Manag.* **2017**, *69*, 24–58. [CrossRef] [PubMed]

3. Kaiser, K.; Schmid, M.; Schlummer, M. Recycling of polymer-based multilayer packaging: A Review. *Recycling* **2018**, *3*, 1. [CrossRef]

4. Cui, Y.; Kundalwal, S.I.; Kumar, S. Gas barrier performance of graphene/polymer nanocomposites. *Carbon* **2016**, *98*, 313–333. [CrossRef]

5. Bharadwaj, R.K. Modeling the barrier properties of polymer-layered silicate nanocomposites. *Macromolecules* **2001**, *34*, 9189–9192. [CrossRef]

6. Wolf, C.; Angellier-Coussy, H.; Gontard, N.; Doghieri, F.; Guillard, V. How the shape of fillers affects the barrier properties of polymer/non-porous particles nanocomposites: A review. *J. Memb. Sci.* **2018**, *556*, 394–418. [CrossRef]

7. Foltynowicz, Z.; Bardenshtein, A.; Sängerlaub, S.; Antvorskov, H.; Kozak, W. Nanoscale, zero valent iron particles for application as oxygen scavenger in food packaging. *Food Packag. Shelf Life* **2017**, *11*, 74–83. [CrossRef]

8. Pant, A.F.; Sängerlaub, S.; Müller, K. Gallic acid as an oxygen scavenger in bio-based multilayer packaging films. *Materials* **2017**, *10*, 489. [CrossRef]

9. Busolo, M.A.; Lagaron, J.M. Oxygen scavenging polyolefin nanocomposite films containing an iron modified kaolinite of interest in active food packaging applications. *Innov. Food Sci. Emerg. Technol.* **2012**, *16*, 211–217. [CrossRef]
10. Wang, C.; Luo, H.; Zhang, Z.; Wu, Y.; Zhang, J.; Chen, S. Removal of As(III) and As(V) from aqueous solutions using nanoscale zero valent iron-reduced graphite oxide modified composites. *J. Hazard. Mater.* **2014**, *268*, 124–131. [CrossRef]

11. Shahwan, T.; Üzüm, C.; Eroğlu, A.E.; Lieberwirth, I. Synthesis and characterization of bentonite/iron nanoparticles and their application as adsorbent of cobalt ions. *Appl. Clay Sci.* **2010**, *47*, 257–262. [CrossRef]

12. Frost, R.L.; Xi, Y.; He, H. Synthesis, characterization of palygorskite supported zero-valent iron and its application for methylene blue adsorption. *J. Colloid Interface Sci.* **2010**, *341*, 153–161. [CrossRef]

13. Uzum, C.; Shahwan, T.; Eroglu, A.; Hallam, K.; Scott, T.; Lieberwirth, I. Synthesis and characterization of kaolinite-supported zero-valent iron nanoparticles and their application for the removal of aqueous Cu²⁺ and Co²⁺ ions. *Appl. Clay Sci.* **2009**, *43*, 172–181. [CrossRef]

14. Gupta, A.K.; Gupta, M. Synthesis and surface engineering of iron oxide nanoparticles for biomedical applications. *Biomaterials* **2005**, *26*, 3995–4021. [CrossRef] [PubMed]

15. Kombaya-Touckia-Linin, E.-M.; Gaucel, S.; Sougrati, M.T.; Khederlou, K.; Pen, N.; Stievano, L.; Gontard, N.; Guillard, V. Hybrid iron montmorillonite nano-particles as an oxygen scavenger. *Chem. Eng. J.* **2019**, *357*, 750–760. [CrossRef]

16. Fan, M.; Yuan, P.; Zhu, J.; Chen, T.; Yuan, A.; He, H.; Chen, K.; Liu, D. Core-shell structured iron nanoparticles well dispersed on montmorillonite. *J. Magn. Magn. Mater.* **2009**, *321*, 3515–3519. [CrossRef]

17. Störmer, A.; Bott, J.; Kemmer, D.; Franz, R. Critical review of the migration potential of nanoparticles in food contact plastics. *Trends Food Sci. Technol.* **2017**, *63*, 39–50.

18. Khalaj, M.; Ahmadi, H.; Lesankhosh, R. Trends in Food Science & Technology Study of physical and mechanical properties of polypropylene nanocomposites for food packaging application: Nano-clay modified with iron nanoparticles. *Trends Food Sci. Technol.* **2016**, *51*, 41–48.

19. Nasiri, A.; Gontard, N.; Gastaldi, E.; Peyron, S. Contribution of nanoclay to the additive partitioning in polymers. *Appl. Clay Sci.* **2017**, *146*, 27–34. [CrossRef]

20. Nasiri, A.; Peyron, S.; Gastaldi, E.; Gontard, N. Effect of nanoclay on the transfer properties of immanent additives in food packages. *J. Mater. Sci.* **2016**, *51*, 9732–9748. [CrossRef]

21. Schlesinger, H.I.; Brown, H.C.; Finholt, A.E.; Gilbreath, J.R.; Hoekstra, H.R.; Hyde, E.K. Sodium borohydride, its hydrolysis and its use as a reducing agent and in the generation of hydrogen. *J. Am. Chem. Soc.* **1953**, *75*, 215–219. [CrossRef]

22. Kim, H.; Kobayashi, S.; Abdurrahim, M.A.; Zhang, M.J.; Khusainova, A.; Hillmyer, M.A.; Abdala, A.A.; Macosko, C.W. Graphene/polyethylene nanocomposites: Effect of polyethylene functionalization and blending methods. *Polymer* **2011**, *52*, 1837–1846. [CrossRef]

23. Turan, D.; Sängerlaub, S.; Stramm, C.; Gunes, G. Gas permeabilities of polyurethane films for fresh produce packaging: Response of O₂ permeability to temperature and relative humidity. *Polym. Test.* **2017**, *59*, 237–244. [CrossRef]

24. Müller, K.; Scheuerer, Z.; Florian, V.; Skutschik, T. Comparison of test methods for oxygen permeability: Optical method versus carrier gas method. *Polym. Test.* **2017**, *63*, 126–132. [CrossRef]

25. Angellier-Coussy, H.; Gastaldi, E.; Gontard, N.; Guillard, V. Influence of processing temperature on the water vapour transport properties of wheat gluten based agromaterials. *Ind. Crops Prod.* **2011**, *33*, 457–461. [CrossRef]

26. NF X31-147 *Qualité des Sols-Sols, Sédiments—Mise en Solution Totale par Attaque Acide*; AFNOR: Paris, France, 1996.

27. NF ISO 14869-1 *Qualité du Sol—Mise en Solution Pour la Détermination des Teneurs Elémentaires Totales-Partie 1: Mise en Solution par L'acide Fluorhydrique et L'acide Perchlorique*; AFNOR: Paris, France, 2001.

28. Grosse, G. PC-Mos, 2nd ed.; Technische Universität München: Munich, Germany, 1993.

29. Cerclé, C.; Sarazin, P.; Favis, B.D. High performance polyethylene/thermoplastic starch blends through controlled emulsification phenomena. *Carbohydr. Polym.* **2013**, *92*, 138–148. [CrossRef]

30. Wang, K.; Wang, C.; Li, J.; Su, J.; Zhang, Q.; Du, R.; Fu, Q. Effects of clay on phase morphology and mechanical properties in polyamide 6/EPDM-g-MA/organoclay ternary nanocomposites. *Polymer* **2007**, *48*, 2144–2154. [CrossRef]

31. De Resende, V.G.; De Grave, E.; da Costa, G.M.; Janssens, J. Influence of the borohydride concentration on the composition of the amorphous Fe-B alloy produced by chemical reduction of synthetic, nano-sized iron-oxide particles. *J. Alloy. Compd.* **2007**, *440*, 236–247. [CrossRef]
32. Shen, J.; Li, Z.; Yan, Q.; Chen, Y. Reactions of bivalent metal ions with borohydride in aqueous solution for the preparation of ultrafine amorphous alloy particles. *J. Phys. Chem.* **1993**, *97*, 8504–8511. [CrossRef]

33. Nurmi, J.T.; Tratnyek, P.G.; Amonette, J.E.; Pecher, K.; Wang, C.; Linehan, J.C.; Matson, D.W.; Penn, R.L.E.E.; Driessen, M.D. Characterization and properties of metallic iron nanoparticles: Spectroscopy, electrochemistry, and kinetics. *Environ. Sci. Technol.* **2005**, *39*, 1221–1230. [CrossRef]

34. Kumar, N.; Auffan, M.; Gattacceca, J.; Rose, J.; Olivi, L.; Borschneck, D.; Kvapil, P.; Jublot, M.; Kaifas, D.; Malleret, L.; et al. Molecular insights of oxidation process of iron nanoparticles: Spectroscopic, magnetic, and microscopic evidence. *Environ. Sci. Technol.* **2014**, *48*, 13888–13894. [CrossRef] [PubMed]

35. Miltz, J.; Perry, M. Evaluation of the performance of iron-based oxygen scavengers, with comments on their optimal applications. *Packag. Technol. Sci.* **2005**, *18*, 21–27. [CrossRef]

36. Lehner, M.; Schlemmer, D.; Sängerlaub, S. Recycling of blends made of polypropylene and an iron-based oxygen scavenger—Influence of multiple extrusions on the polymer stability and the oxygen absorption capacity. *Polym. Degrad. Stab.* **2015**, *122*, 122–132. [CrossRef]

37. Lagarón, J.-M.; Busolo, M.A. Active nanocomposites for food and beverage packaging. In *Emerging Food Packaging Technologies*; Yam, K.L., Lee, D.S., Eds.; Elsevier: Philadelphia, PA, USA, 2012; pp. 55–65.

38. Amberg-Schwab, S.; Weber, U.; Burger, A.; Nique, S.; Xalter, R. Development of passive and active barrier coatings on the basis of inorganic-organic polymers. *Monatshefte Chem.* **2006**, *137*, 657–666. [CrossRef]

39. Poley, L.H.; Ana, P.; Siqueira, L.; Silva, M.G.; Vargas, H. Photothermal characterization of low density polyethylene food packages. *Polímeros Ciência Tecnol.* **2004**, *14*, 8–12. [CrossRef]

40. Mrkić, S.; Galić, K.; Ivanković, M.; Hamin, S.; Ciković, N. Gas transport and thermal characterization of mono- and di-polyethylene films used for food packaging. *J. Appl. Polym. Sci.* **2006**, *99*, 1590–1599. [CrossRef]

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

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license ([http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)).