Sun Hua,
Feng Chen,
Zheng-ying Liu,
Wei Yang and
Ming-bo Yang*
College of Polymer Science and Engineering, State Key Laboratory of Polymer Materials Engineering, Sichuan University, Chengdu 610065, Sichuan, People's Republic of China. E-mail: yangmb@scu.edu.cn
First published on 21st December 2015
In order to improve the melt strength of poly(lactic acid) (PLA), microcrystalline cellulose-graft-polylactic acid (MCC-g-PLA) copolymer was prepared and introduced into the PLA matrix. The MCC-g-PLA copolymers were synthesized by melt copolycondensation of lactic acid (LA) with microcrystalline cellulose (MCC) which was pretreated to improve accessibility. The MCC-g-PLA copolymers with a molar substitution (MS) of PLA in the range of 1.67–5.97 were synthesized by adjusting reaction temperature, molar ratio of LA monomer to MCC and washing times in the pretreatment process. Compared with MCC, the crystalline structures of MCC-g-PLA copolymers were not perfect due to the existence of PLA side chains. A glass transition temperature (Tg) appeared in the copolymers, which had never been observed in unmodified MCC, and decreased with the increase of MS. When MS was above 4.41, the thermal degradation of PLA side chains emerged. When MCC-g-PLA copolymer was introduced into the PLA matrix, good dispersion of MCC-g-PLA was verified with SEM results. The extension rheology results showed that the melt strength of PLA can be effectively enhanced with the addition of MCC-g-PLA, especially at low elongation rate. Meanwhile, MCC-g-PLA also improved the crystallization ability of PLA in the non-isothermal crystallization process.
Although PLA possesses many advantages, the drawbacks of the poor viscoelastic behavior, low melt strength, the inherent brittleness, and slow crystallization rate, have limited its large-scale commercialization especially as film materials.14–16 As is known, blow molding, by which a bi-axial stress is imposed on the film bubble, is the most convenient and high-performance-price-ratio method to manufacture films. From such a process, impact properties of the film are excellent and mechanical and optical properties can be controlled.17 However, polymer melt undergoes a complex three-dimensional deformation process during blow molding, so to prepare blown film steadily and successfully, high melt strength is the basic requirement for materials, such as the case for branched low density polyethylene (LDPE). Unfortunately, the blow molding of PLA is quite difficult owing to its poor viscoelastic behavior and low melt strength. The brittleness of PLA also prevents it from being used as plastic bags and plastic film. Therefore, in order to improve the processability of PLA to be used as blown films and to prepare PLA films with good properties, PLA needs to be modified to improve the melt strength, viscoelastic behavior and toughness.16,18
Increasing the molecular weight and introducing long branched chains have been proved to be effective methods to improve the melt strength of PLA, which can increase chain entanglement and extend melt relaxation time of PLA chains. Although high molecular weight (Mw) PLA can be synthesized by ring-opening polymerization of lactide and chain extension in laboratory studies,19,20 Mw of commercial PLA is still not high enough to ensure the melt strength needed for blow molding. Introducing the long branched chains is another effective method, which also has a great influence on the viscoelastic behavior of polymer melt, and has been proved to be able to increase the melt strength and chain entanglement.21–23 So, if we can synthesized copolymer with long branched chains of PLA, the existence of PLA branched chains can not only increase entanglement of PLA chains and melt strength of PLA, but also improve the dispersion state of the copolymer in PLA, thus making PLA easy to process.
Cellulose is the most abundant, biodegradable, and renewable natural polymer in the earth, and holds many attractive properties, such as low cost, nontoxic, good mechanical properties and low density.24,25 If we can graft long branched chains of PLA onto cellulose, the biodegradable property of products is guaranteed, simultaneously, the applications of PLA and cellulose are both widened.
Graft copolymerization is a common way to introduce other polymers onto cellulose or cellulose derivatives.26–32 Unfortunately, cellulose has strong intramolecular and intermolecular hydrogen-bonding network and high crystallinity, which make it non-thermoplastic and insoluble in common solvent. Most reaction reagents penetrate only into the amorphous regions of cellulose,33 while a lot of hydroxyl groups are sealed in crystalline region, then it is difficult to access to reagents and participate in reaction. These restrict the reaction activity of cellulose. In recent years, new solvents of cellulose have been discovered, such as lithium chloride/N,N-dimethylacetamide and ionic liquid.34–36 They had been used as reaction media to graft polymer chains onto cellulose,37–42 which could improve the accessibility and reactivity of cellulose. But the use of solvents could bring some problems, such as toxicity, complicated preparation process, or degradability.
How can we improve the accessibility of cellulose in non-toxic solvent? According to previous studies,43–46 treating cellulose with aqueous sodium hydroxide (NaOH) solutions had a considerable effect on the hydrogen-bonding network and crystalline structure of cellulose, in this case the crystalline structure of cellulose can transform from cellulose I to cellulose II. In this process, the alkali can destroy the hydrogen bonds and infiltrate into the crystalline region, which could be utilized to improve the accessibility and activity of cellulose. Then the solvent can be exchanged. Thus, the reaction monomer could get inside of cellulose.
Because lactide is easy to decompose, lactic acid was chosen as reaction monomer. In this work, firstly we improved the accessibility of cellulose by solvent exchanged method. Then, melt copolycondensation of lactic acid with microcrystalline cellulose (MCC) was carried out to obtain microcrystalline cellulose-graft-polylactic acid copolymers using tin(II) chloride dehydrate catalyst. Finally, the synthesized copolymer was added into PLA matrix, and the properties of PLA composites were studied.
As a contrast, using the same process, we used water to precipitate the dissolved cellulose and wash the precipitate, then the product was dried in an oven at 60 °C, and was coded as MCC-W. The weight of MCC-W was 0.9 g.
Sample | Washing times | LAa/AHG (mol mol−1) | T (°C) | DPs | DS | MS |
---|---|---|---|---|---|---|
a Indicates the amount added in melt copolycondensation process, which does not include the LA remained in MCC-L.b Specific reaction temperature and time: 120 °C, 2 h; 130 °C, 6 h.c Specific reaction temperature and time: 120 °C, 2 h; 130 °C, 2 h; 140 °C, 4 h.d Specific reaction temperature and time: 130 °C, 2 h; 140 °C, 2 h; 150 °C, 4 h. | ||||||
MP-C0W4-T130 | 4 | 0 | 130b | 2.79 | 0.60 | 1.67 |
MP-C0W4-T140 | 4 | 0 | 140c | 3.40 | 0.60 | 2.04 |
MP-C0W4-T150 | 4 | 0 | 150d | 3.47 | 0.64 | 2.22 |
MP-C0W4-T160 | 4 | 0 | 160 | 4.67 | 0.67 | 3.13 |
MP-C1W4-T160 | 4 | 18/1 | 160 | 5.17 | 0.71 | 3.67 |
MP-C2W4-T160 | 4 | 36/1 | 160 | 4.60 | 0.64 | 2.94 |
MP-C2W3-T160 | 3 | 36/1 | 160 | 4.70 | 0.88 | 4.14 |
MP-C2W2-T160 | 2 | 36/1 | 160 | 5.90 | 0.88 | 5.19 |
MP-C2W1-T160 | 1 | 36/1 | 160 | 6.22 | 0.96 | 5.97 |
The thermal analysis of PLA/MCC-g-PLA composites was also performed by DSC. The samples of 5–7 mg were dried at 60 °C for 24 h before measurement, and were heated from 25 °C to 200 °C at a heating rate of 10 °C min−1 and kept at 200 °C for 3 min to eliminate the thermal history, then cooled to 25 °C at a rate of 5 °C min−1, and finally were reheated to 200 °C at a rate of 10 °C min−1 under nitrogen atmosphere. The cooling and second heating curves were recorded.
In cellulose, the hydrogen-bonding network makes molecular chains arrange regularly and form crystal. Different from cellulose I with parallel chain structure, cellulose II was comprised of an array of antiparallel chain molecules, and the unit cell was monoclinic according to Kolpak and Blackwell.44 When MCC was treated in NaOH solution, alkali reached to amorphous and crystalline regions of MCC, and destroyed the hydrogen bonds in the crystalline regions, then MCC was dissolved in NaOH solution.47 When dissolved MCC was regenerated in water, the conformation of hydroxymethyl group changed,43–45 and cellulose I transformed to cellulose II with an antiparallel chain structure. So MCC-W was cellulose II. But for cellulose II, the hydrogen-bonding network and structure were also very dense, and the reaction reagents were still difficult to penetrate into the inside of cellulose and react with cellulose. So, simple alkali treatment could not improve the accessibility of cellulose. Therefore, if we wanted to introduce reaction reagents into cellulose, we had to use the stage when hydrogen bonds were destroyed.
According to the above analysis, after MCC was dissolved in NaOH solution, hydrogen bonds in cellulose were destroyed by alkali. At this time, if the reaction reagents were added into the mixture, reaction reagents may enter into the cellulose before the hydrogen bonds regenerated. For MCC-L, perfect crystalline structure was not observed, indicating that the crystal structure and hydrogen-bonding network had been disrupted. This may be owing to that adding lactic acid into the dissolved MCC solution did make lactic acid enter into cellulose, and the lactic acid remained in cellulose after drying at 60 °C. The boiling point of lactic acid was about 122 °C,48 and it was not easy to evaporate at 60 °C. This was also confirmed by increase in the weight of MCC-L compared with MCC-W. The remained LA molecules formed steric hindrance to hinder the configuration change of C–O bond at C6 position and the formation of hydrogen bonds, and broke the regularity of hydrogen-bonding network in cellulose, thus MCC-L cannot crystallize. That is, lactic acid had reached into the crystalline region of cellulose, and had a possibility of reacting with hydroxyl groups both in amorphous and crystalline region, and the accessibility of cellulose to lactic acid was improved.
The structure of MCC-g-PLA copolymers was characterized by FTIR and the results were presented in Fig. 2. Because the byproduct PLA homopolymer produced in melt copolycondensation process had been removed by purification process, the FTIR spectra should reflect the chemical structure of the grafted copolymers. As shown in Fig. 2, the strong peak at around 3410 cm−1 was assigned to the absorption of hydroxyl groups in cellulose. Compared with FTIR result of MCC, a new band at about 1750 cm−1 was observed in each grafted copolymer spectrum, which was assigned to the carbonyl group. The methyl asymmetric deformation in PLA side chains at around 1475 cm−1, and the symmetric C–O–C stretching modes in ester groups at around 1200 cm−1 also appeared in the spectra of the copolymers.49 The results clearly proved that PLA chains had been grafted onto cellulose backbone by melt copolycondensation of LA with hydroxyl groups in cellulose and MCC-g-PLA copolymers were synthesized successfully, as shown in Scheme 1.
Fig. 2 FTIR spectra of (a) MCC, (b) MP-C0W4-T130, (c) MP-C0W4-T140, (d) MP-C0W4-T150, and (e) MP-C0W4-T160. |
The structure of the copolymers was also investigated by 1H NMR. An 1H NMR spectrum obtained for one MCC-g-PLA copolymer (MP-C0W4-T160) was displayed in Fig. 3 and the signal peaks were assigned as follows: δ = 1.32 ppm (methyl protons of terminal lactyl); δ = 1.50 ppm (methyl protons of internal lactyl); δ = 4.23 ppm (methine protons of internal lactyl); δ = 5.07–5.29 ppm (methine protons of terminal lactyl); δ = 3–5.4 ppm (protons of anhydroglucose). The appearance of proton signals of lactyl groups in this spectrum, further confirmed that MCC-g-PLA copolymers were synthesized successfully.
The contents of grafted PLA in MCC-g-PLA copolymers were calculated according to 1H NMR spectra. As shown in Fig. 3, we designated the resonance peak area derived from methine protons of internal lactyls in PLA side chains as A, and the resonance peak area from methine protons of terminal lactyls as B. The resonance peak area derived from methyl protons of internal lactyls was labeled as C, and the area from methyl protons of terminal lactyls was labeled as D. The integral of all protons of anhydroglucose (AHG) unit was designated as E. Then the average degree of polymerization of a PLA side chain (DPs), the degree of lactyl substitution (DS, the average number of hydroxyl groups substituted for lactyls per AHG of cellulose), and the molar substitution (MS, the average number of introduced lactyl units per AHG of cellulose) were estimated by eqn (1)–(3),38,40 respectively, and the results were summarized in Table 1.
DPs = (IC + ID)/ID | (1) |
DS = 7ID/3[IE − (IC + ID)/3] | (2) |
MS = DPs × DS | (3) |
According to Table 1, the amount and length of PLA side chains can be controlled by adjusting reaction conditions and molar ratios of LA monomer to MCC, and the relationship between them will be discussed in detail subsequently.
From Table 1, the difference of the reaction conditions of samples MP-C0W4-TX (X = 130, 140, 150, 160) was reaction temperature, where the average numbers of hydroxyl groups substituted for lactyls per AHG (DS values) were almost constant, ranging from 0.60 to 0.67 with increasing reaction temperature, meaning that the accessibility of LA monomer to cellulose kept the same. Based on previous analysis, the accessibility of cellulose was mainly affected by solvent exchange pretreatment process. So the influence of reaction temperature on DS was small. However DPs values of the four samples increased with the increasing of reaction temperature. Lactic acid melting polycondensation reaction was an equilibrium reaction, producing PLA and water, so removing water was beneficial to the conversion of lactic acid into PLA.50,51 Under the same vacuum condition, high temperature could accelerate the removal of water, including free water and that formed in reaction process. Besides, elevated temperature could also increase the thermal degradation rate. Therefore, the influence of temperature on reaction was a combination of these factors. In the temperature range of our experiments, the DPs value increased with the temperature rising.
Samples MP-CXW4-T160 (X = 0, 1, 2) were synthesized by adjusting the molar ratios of LA monomer to MCC. According to Table 1, with increasing molar ratios of LA/AHG, DS and DPs values of the copolymers showed no obvious changes, indicating that increasing the content of LA monomer had almost no effect on the copolymers. Theoretically, increasing the content of reactant was helpful to produce the product for an equilibrium reaction. But in our work, the melt copolycondensation happened in a heterogeneous condition, and the steric hindrance influenced the reaction, which made the added monomer not participate in the reaction completely, so the influence of monomer content on graft copolymers was not so significant.
Except for the reaction process, the pretreatment process, in which we used solvent exchanged method to change the accessibility of cellulose, also influenced the formation of copolymers. The solvent exchanged process was related to the washing step with lactic acid in the pretreatment process, then we wondered whether the washing times with lactic acid had influence on the graft copolymers. So we synthesized samples MP-C2WX-T160 (X = 1, 2, 3, 4) through the same reaction process, only changing the washing times with lactic acid in pretreatment process.
From Table 1, we observed the dependence of DS and DPs on washing times. Obviously, DS and DPs increased gradually when the number of washing times reduced, indicating that the consumption of lactic acid can be reduced. During the process that the dissolved cellulose was regenerated and washed with lactic acid, firstly, lactic acid entered into cellulose, and combined with alkali, then the alkali was brought out of cellulose, at this time hydroxyl groups surrounded by alkali were freed out. These free hydroxyl groups probably formed hydrogen bonds with other free hydroxyl groups, or probably were surrounded by lactic acid again. If hydroxyl group contacted with lactic acid, it had the possibility of reacting with LA. But, once it formed hydrogen bond with other free hydrogen group, it was hard to take part in a reaction because the hydrogen bonds were difficult to be broken again. So in each washing process, some hydroxyl groups which were surrounded with alkali or lactic acid were freed out again and then formed hydrogen bonds, thus losing the possibility of reaction. So with the increase of washing times, DS value decreased gradually. The hydrogen bonds also limited the diffusion of lactic acid in cellulose, which affected the amount of lactic acid that could participate in reaction, causing the decrease of DPs.
The crystalline structures of MCC and MP-C2W1-T160 were studied by WAXD and the results were shown in Fig. 4 (Spectra of other copolymers were similar to MP-C2W1-T160 and not shown here for brevity). MCC was a typical cellulose I, as mentioned above. Different from WAXD spectrum of MCC, MP-C2W1-T160 showed two broad and small peaks at 19.8°and 22.0°, which were the characteristic diffraction peaks of cellulose II. The intensity of these two peaks was very weak, and the characteristic peak at 12.1° did not appear, indicating that MCC-g-PLA copolymers presented an imperfect crystalline structure of cellulose II. According to Kolpak and Blackwell,44 for cellulose II intramolecular hydrogen bonds (O2′–H⋯O6) and intermolecular hydrogen bonds (O6–H⋯O3, O6–H⋯O2) were in the 020 plane, and intermolecular hydrogen bonds (O2–H⋯O2′) were in the 110 plane. These hydrogen bonds made the molecular chains arrange regularly and form a perfect crystal. Reaction between LA and hydroxyl groups reduced the amount of hydroxyl groups capable of forming hydrogen bonds. On the other hand, the existence of PLA side chains generated steric hindrance, also affecting the formation of hydrogen bonds and regular arrangement of molecular chains. So the hydrogen-bonding network was weakened and MCC-g-PLA copolymers could not form perfect crystals. These results also proved that copolycondensation reaction occurred both in amorphous and crystalline region. Therefore we had improved the reactivity of cellulose by this method and synthesized MCC-g-PLA successfully.
Sample | DPs | DS | MS | Tg (°C) | Tonset (°C) | Tmax (°C) |
---|---|---|---|---|---|---|
MP-C0W4-T130 | 2.79 | 0.60 | 1.67 | 77.50 | 248.1 | 317.8 |
MP-C0W4-T140 | 3.40 | 0.60 | 2.04 | 70.50 | 244.2 | 315.8 |
MP-C0W4-T150 | 3.47 | 0.64 | 2.22 | 70.50 | 253.7 | 312.9 |
MP-C0W4-T160 | 4.67 | 0.67 | 3.13 | 60.93 | 264.7 | 316.7 |
MP-C1W4-T160 | 5.17 | 0.71 | 3.67 | 58.56 | 257.7 | 321.5 |
MP-C2W4-T160 | 4.60 | 0.64 | 2.94 | 61.94 | 255.1 | 321.9 |
MP-C2W3-T160 | 4.70 | 0.88 | 4.14 | 62.79 | 242.3 | 285.7, 321.6 |
MP-C2W2-T160 | 5.90 | 0.88 | 5.19 | 59.99 | 239.4 | 277.1, 313.2 |
MP-C2W1-T160 | 6.22 | 0.96 | 5.97 | 57.42 | 216.4 | 256.7, 301.7 |
MCC | 320.1 | 347.6 |
TG measurement was used to evaluate the thermal stability of copolymers and the results were displayed in Fig. 6 and Table 2. Obviously, the onset decomposition temperature (Tonset) and the maximum decomposition temperature (Tmax) of MCC-g-PLA copolymers were lower than those of MCC. It was owing to that the introduction of PLA side chains on the cellulose backbone destroyed the hydrogen-bonding network and crystalline structure of MCC to some extent, resulting in the decrease of thermal stability. This was also confirmed by DSC and WAXD results. On the derivative thermogravimetry (DTG) curves, there was only one peak for copolymers when MS < 4.14, demonstrating that the thermal degradation of copolymers proceeded in one step. Because the content of PLA side chains was low, the decomposition peak of PLA chains did not emerge, and the peak in DTG curve was mainly caused by decomposition of cellulose backbone. However a shoulder peak appeared near the main peak in DTC curves of MP-C1W4-T160 and MP-C2W4-T160, showing that thermal decomposition process changed gradually. When MS ≥ 4.14, two derivate thermogravimetry peaks appeared, indicating that thermal degradation of these samples proceeded in two steps. This should be attributed to the different decomposition temperatures of MCC backbone and PLA side chains, and the degradation occurring at low temperature was caused by PLA side chains, and that at high temperature was caused by cellulose backbone, which was also reported in literature.40
Fig. 7 SEM micrographs of (a) PLA, (b) PLA-MP-2, (c) PLA-MP-3, (d) PLA-MP-4, (e) PLA-MCC-2, (f) PLA-MCC-3, and (g) PLA-MCC-4. |
In order to investigate the influence of MCC-g-PLA on the melt strength of PLA, the extensional rheological properties of PLA and PLA/MCC-g-PLA composites were characterized by the extensional rheology at the extensional strain rates of 0.1 s−1 and 0.5 s−1 in steady uniaxial extension, and the results were shown in Fig. 8. It was obvious that the elongational viscosity of neat PLA kept almost constant at the full time of stretching, and the strain-hardening behavior could not be observed, which was the typical characteristic of linear polymer. The elongational viscosity was improved with the addition of MCC-g-PLA at low extensional strain rate (0.1 s−1), and the increment increased with increasing content of MCC-g-PLA, which was attributed to the entanglement between MCC-g-PLA and PLA matrix. According to the 1H NMR analysis of MCC-g-PLA, the average degree of polymerization of PLA side chain was 6.22, that is the average molecular weight was 425 g mol−1, however the critical molecular weight for entanglement of PLA was near 9000 g mol−1,55 the length of PLA side chains in MCC-g-PLA was not long enough to form strong entanglement with PLA matrix. So the entanglement should be between MCC-g-PLA filler particles and PLA matrix, but not between the PLA side chain and PLA matrix. However, at high extensional strain rate (0.5 s−1), when the content of filler was low (sample PLA-MP-2 and PLA-MP-3), the elongational viscosity could not be improved and was approximately equal to that of neat PLA; only PLA-MP-4 showed higher elongational viscosity than neat PLA. According to the SEM results, MCC-g-PLA was well dispersed in the matrix, so the higher the content of MCC-g-PLA, the more the entanglement in the composites. For PLA-MP-2, due to the low entanglement, the elongational viscosity was almost equal to neat PLA even at low extensional strain rate. For PLA-MP-3, the entanglement was increased, so the elongational viscosity increased at low extensional strain rate. But at high extensional strain rate, the elongational viscosity decreased because of disentanglement and was approximately equal to neat PLA indicating hardly any entanglement between filler and matrix. For PLA-MP-4, the elongational viscosity decreased with the increase of extensional strain rate, but was still higher than neat PLA, showing that there was entanglement remained in spite of the existence of disentanglement.
Fig. 8 Variation of elongation viscosity as a function of time for PLA and PLA/MCC-g-PLA composites at 158.5 °C with elongation rates of (a) 0.1 s−1 and (b) 0.5 s−1. |
The effects of MCC-g-PLA on the non-isothermal crystallization and melting behaviors of PLA matrix were also investigated by DSC. Cooling curves of different samples were shown in Fig. 9(a). For pure PLA, there was no exothermic peak representing the crystallization of PLA in the cooling curve due to the relatively rigid molecular chain and slow crystallization rate. However, an exothermic peak emerged when MCC-g-PLA was introduced, and the peak became stronger with increasing content of MCC-g-PLA, except for sample PLA-MP-4. In the second heating process (Fig. 9(b)), there was no cold crystallization peak, which was also observed in literature,56,57 and melting peak for pure PLA. For PLA/MCC-g-PLA composites, we observed the appearance of cold crystallization peak, and the cold crystallization temperature (TCC) shifted to the lower temperature with increasing content of MCC-g-PLA, and returned to higher temperature for PLA-MP-4. These results indicated that the addition of MCC-g-PLA acted as nucleation agent and promoted the crystallization ability of PLA. At first, the nucleation ability of MCC-g-PLA was enhanced with the increasing content of nucleation agent, thus making crystallization peak become stronger and TCC shift to lower temperature. Then the content of MCC-g-PLA increased further, the entanglement between PLA and MCC-g-PLA also increased, resulting in the hindered mobility of PLA molecular chains. This was in agreement with the extensional rheology test. In addition, the shape of melting peak also changed with the content of MCC-g-PLA. As shown in Fig. 9(b), there were double melting peaks in the second heating curves of PLA/MCC-g-PLA composites, which was attributed to the melt/recrystallization behavior of PLA.58,59 The melting peak at lower temperature was caused by the melting of imperfect crystals and the melting peak at higher temperature was caused by the melting of the relatively perfect crystal and recrystallized crystals. Taking PLA-MP-3 as an example, the crystallization ability was the strongest, and the amount of relatively perfect crystal formed in crystallization process was more, so the higher temperature melting peak was intensified.
Fig. 9 DSC thermograms of pure PLA and PLA/MCC-g-PLA composites (a) for the first cooling process, (b) for the second heating process. |
MCC-g-PLA copolymer can be well dispersed in PLA matrix, and the melt strength of PLA can be improved by the addition of MCC-g-PLA at low extensional strain rate, but the improvement decreased at high extensional strain rate, due to the weak entanglement. MCC-g-PLA also acted as nucleation agent to enhance the crystallization ability of PLA. The results in this work proved that the performance of PLA can be improved by this method, but the length of PLA side chains was short, if we improve its length over the critical entanglement length, the PLA side chain can also form entanglement with matrix, thus making melt strength increase significantly, so the studies on further improving the length of PLA side chains are now in progress.
This journal is © The Royal Society of Chemistry 2016 |