Alistair W. T.
King
*a,
Jarno
Jalomäki
a,
Mari
Granström
b,
Dimitris S.
Argyropoulos
c,
Sami
Heikkinen
a and
Ilkka
Kilpeläinen
*a
aLaboratory of Organic Chemistry, Department of Chemistry, University of Helsinki, P.O. Box 55 (A. I. Virtasen Aukio 1), FIN-00014, Finland. E-mail: alistair.king@helsinki.fi; ilkka.kilpeläinen@helsinki.fi; Tel: +358 50527 9446
bBASF SE, GCI/R, M311 67056, Ludwigshafen, Germany
cDepartment of Forest Biomaterials, North Carolina State University, Raleigh, NC 27695-8005, USA
First published on 31st August 2010
Chloroform-soluble palmitic and decanoic acid esters of cellulose were synthesized from the reaction of MCC with acid chlorides in LiCl/DMA and the ionic liquid [amim]Cl, as novel cellulose solvents. A process of derivatization of the remaining hydroxyl groups, as phosphite esters and subsequent 31P NMR analysis, allowed for simultaneous degree of substitution (DS) determination and quantification of the aqueous-quench acid by-product impurity, after the appropriate calculation. The full mathematical treatment for DS determination is presented, including scripts for the Python and Java programming languages for rapid interpretation of results. This method for DS determination was validated against traditional analyses of the palmitoyl cellulose and the fully substituted p-nitrobenzoyl palmitoyl diester product, from additional reaction with p-nitrobenzoyl chloride. DOSY NMR and SCORE analyses were also employed to demonstrate the utility of this rapid 31P derivatization and analytical process over traditional 1D NMR analyses.
O stretch of the conjugate base–acid pair is shifted out of the region of the ester allowing increased accuracy in acid impurity determination. The IL 1-ethyl-3-methylimidazolium acetate ([emim][OAc]) was used for dissolution of their samples, as one of a range of accessible ILs capable of dissolving wood and cellulose,5 allowing for analysis of esters of cellulose in the low to medium range (0–1.5) in addition to those which are fully organically soluble.
![]() | ||
| Fig. 1 Typical fatty acid esterification of cellulose to give the ester and acid by-product, from [amim]Cl 1, as reaction media. | ||
Despite the high novelty and utility of this method, IR analysis is still limited in its accuracy and in the information that it provides. However, the use of cellulose-dissolving media, such as ILs, in the analysis procedure is an important concept as it may allow for the analysis of normally insoluble products such as low DS cellulose esters.
In the present publication, it is our objective to validate and demonstrate the utility of 31P derivatization and NMR analyses for the optimisation of cellulose functionalization reactions, in comparison to traditional methods. For this, we have initially chosen a short set of long-chain cellulose esters in the organically soluble DS range of ∼1.5 to 3 and of variable purity. Based upon our existing and growing knowledge of wood and wood biopolymer dissolution and functionalization in ILs, aided by 31P labelling and NMR analysis,6,7 it is our goal ultimately to develop this technique for the study of a wider range of wood bio-polymer chemical modification reactions.
000 Hz spectral width, 1 s acquisition time, and 10 s relaxation delay. The quantitative 13C spectrum was collected with 20
000 transients utilizing a 90° pulse flip angle, 36
000 Hz spectral width, 1 s acquisition time and 60 s relaxation delay. Quantitative 1H spectra and DOSY (diffusion-ordered spectroscopy) data were recorded at 27 °C on a Varian Unity Inova 500 spectrometer (500 MHz proton frequency) equipped with 5 mm triple-resonance (1H, 13C, and 15N) gradient probe-head. For 1H measurements, 60 transients were collected with 45° pulse flip angle, 1.9 s acquisition time, and 13 s relaxation delay. Spectral width for the 1H spectra was 8000 Hz. DOSY datasets were recorded using Bipolar Pulse Pair Stimulated Echo pulse sequence (BPPSTE)8 using 1 ms gradient duration, 0.2 ms gradient recovery delay, 150 ms diffusion time, 224 transients, 1 s acquisition time, 2 s relaxation delay, and spectral width of 8000 Hz. In DOSY datasets, 60 diffusion gradient amplitudes were used (ranging from 0.5 to 20 G cm−1). DOSY data were analyzed with the Speedy COmponent REsolution (SCORE)-algorithm9 incorporated into the DosyToolBox 0.53-software10 running on Matlab 7.5.0 (MathWorks, Natick, MA, USA). SCORE analysis was carried out with standard parameters over the region of 0–6 ppm for a 2-component system.
:
pyridine (1
:
1, 2 ml). p-Nitrobenzoylchloride (150 mg, 0.809 mmol) was added in one portion and the mixture was allowed to stir at room temperature for 18 h under inert atmosphere. Methanol was added until the product precipitated and the solution became clear. The solution was removed with a pipette and this precipitation process repeated 2 times. The solid was then finally precipitated from chloroform
:
hexane (1
:
1) with methanol addition. The solution was decanted and the product dried under vacuum at 55 °C for 18 h to give high purity p-nitrobenzoyl-palmitoyl cellulose 3 (128 mg, 28% WPG, 97% yield at a combined DS of 3), as a white solid (found: C, 65.64, H, 8.7, N, 2.29. Calc. for C42.03H64.09N1.33O10.66 (DS 1.67): C, 72.6, H, 9.3, N, 2.7%); νmax (ATR)/cm−1 2922, 2853 and 1741; δH (500 MHz, CDCl3, Me4Si) 0.68–1.09 (3H, br t, J 6.8, CH3), 1.09–1.98 (26H, br m, COCH2(CH2)13CH3), 1.98–2.70 (2H, br m, COCH2(CH2)13CH3), 2.70–5.62 (7H, br m, AGU). δC (125 MHz, CDCl3) 14.46, 23.05, 25.01, 29.28–30.39, 32.29, 34.04, 61.84, 70.75–75.81 (br m), 100.83, 123.91, 130.62, 134.25, 151.18, 162.65–164.37 (br m), 171.96–173.46 (br m); DS1H: 1.67 (calculated from 1H NMR analysis); DS13C: 1.67 (calculated from 13C NMR analysis); DSEA: 1.53 (calculated from C-content determined from the EA); DSEA: 1.75 (calculated from N-content determined from the EA).
:
CDCl3/3
:
2, 0.0152 mmol) was added in one portion and the solution was vortexed (∼10 s). 31P NMR spectra (243 MHz for 31P) were recorded with 700 µl samples, in a 5 mm o.d. NMR tube.
A description for the calculation of DS values and impurities from cellulose functionalization reactions follows; the experimental variables required for DS determination from phosphite ester derivatization and 31P NMR analysis are: (1) MWS (substituent molecular weight). This is the molecular weight (g mol−1) of the substituent, not including the linking oxygen atom (between the cellulose backbone and the substituent); (2) ISvol (internal standard volume, µl); (3) ISmol (internal standard molarity, mmol); (4) IR (integration ratio of remaining functionalized cellulose hydroxyls against internal standard); and (5) WS (sample weight, mg). In addition, the constants DSmax (maximum achievable DS value of 3 for unsubstituted cellulose) and OHC (the free hydroxyls per weight unit of cellulose, 3 × 162 = 0.01852 mol g−1) are also required. The eqn (1) for calculation of the free hydroxyls per weight unit of substrate (OHS, mol g−1) is as follows:
![]() | (1) |
The final eqn (2) required to calculate the DS from 31P labelling and NMR analysis (DS31P), is as follows:
![]() | (2) |
In addition to calculating the DS of pure products, the integration ratio (against internal standard) for any impurity (II), e.g. carboxylic acids, in combination with the molecular weight of the impurity (MWI, g mol−1) can be used to determine the weight of the impurity (WI, mg) in the original sample. This is determined using eqn (3):
![]() | (3) |
If impurity exists in a sample a corrected DS value can be eventually attained by simply subtracting WI from the original sample weight (WS) and recalculating OHS, as shown in eqn (4):
![]() | (4) |
The above equations may be adapted to starting polymers other than unsubstituted cellulose, for any given reaction. A good example of this is if an additional substituent is added to an existing substituted cellulose product. In this case, if you have predetermined the DS value of the initial substituent, by 31P labelling and NMR analysis, one can analyse the product of the second substitution reaction in a similar manner. To determine the DS of the second substituent from these results you can replace OHC in the second calculation of eqn (2) with OHS from the first calculation of eqn (1), and for the second calculation of eqn (2), DSmax is replaced with the DS of the first substituent. The result will be DS31P values for both substituents from analysis of both products. Other non-cellulosic polymers or polymer mixtures, such as xylan, lignin, starch or wood, may also be analysed by replacing OHC and DSmax values from their experimental or theoretically calculated values. As some polymers may have a non-integer number of hydroxyls per monomer unit, or non-regular monomer units, DS31P may be represented as a percentage substitution of the total sample by replacing DSmax with 100.
Eqn (1) and (2) have been incorporated in to scripts for the Python and Java programming languages13 in order to allow for rapid processing of data (see ESI†).
![]() | ||
| Fig. 2 Synthesis of palmitoyl cellulose 2 and p-nitrobenzoyl-palmitoyl cellulose 3. | ||
As a secondary method, outside 31P labelling and analysis, the remaining hydroxyls in the palmitoyl cellulose product 2 were further esterified as p-nitrobenzoyl esters. This was achieved by quantitatively reacting the chloroform soluble starting material with the p-nitrobenzoyl chloride (Fig. 2) to give, after purification, fully substituted p-nitrobenzoyl-palmitoyl cellulose 3. The reaction was determined to be complete (almost no remaining free hydroxyl groups) by the absence of any OH-stretch in the IR spectra (see ESI†) or absence of any aliphatic phosphite ester resonances, after 31P derivatization and NMR analysis. In addition to providing an alternative DS1H value, this functionality also provided the opportunity to analyse the product by EA, observing both carbon and nitrogen percentages, and quantitative 13C NMR, with the carbonyl resonances providing a suitable region for integration and DS13C determination (see ESI†). The inclusion of the chromophore also would allow for GPC with UV detection, although this is not a suitable method for determining DS values. Although it was possible to determine DS13C by straightforward integration of the carbonyl resonances, with sufficient exponential line broadening, the signal to noise ratios and resolution, after 20
000 transients (∼66 h collection time), were not adequate enough to allow for accurate determination of regioselectivity, as has been previously observed for cellulose acetates.14
31P derivatization and NMR analysis of palmitoyl cellulose 2, shown in Fig. 3, were performed by dissolution of the substrate 2 into CDCl3 and reaction with 2-Cl-TMDP under basic (pyridine) conditions. The internal standard e-HNDI was also included in the phosphitylating mixture, which also reacted with the phosphite acid chloride (2-Cl-TMDP) or anhydride (TMDP-anhydride), in solution, to form its corresponding phosphite ester (e-HNDI-TMDP). The spectrum was calibrated with TMDP-anhydride at 132.2 ppm. e-HNDI-TMDP resonates at 152.0 ppm under these solvent conditions. Alkoxy phosphite esters (alkoxy-TMDP) are located in the region between 151.5 and 143.5 ppm and both decanoic acid and palmitic acid–phosphite mixed anhydrides (decanoate or palmitate–TMDP) resonate at 134.9 ppm. The reaction was performed in triplicate and the original sample weights were corrected, based upon the presence of 0.94% w/w palmitic acid impurity (as determined by the 31P NMR integrations of palmitate–TMDP). DS31P for this compound 2 was determined to be 1.67 with a standard deviation for the full procedure of 0.012. It is expected that as you approach DS values close to 0 or have very large substituents, the standard deviation will increase although for the chloroform soluble regime (DS ca. 1–3) the error will remain approximately the same.
![]() | ||
| Fig. 3 Derivatization mixture and 31P NMR analysis of TMDP-phosphitylated palmitoyl cellulose 4. | ||
From the initial palmitoyl cellulose substrate 2 and its subsequent characterization, the complete set of determined DS values are listed in Table 1, including short self-explanatory descriptions of the methods used to obtain them. From a comparison of these methods, we can see that EA predictably shows the most deviation from the mean value. This is very much dependent on both purity of materials and skill of the analyst. Due to the simplicity of integration, and the ability to resolve low molecular weight hydroxylated impurities from the polymeric material in solution, 31P derivatization and NMR analysis offer a more accurate and thorough characterization method, in combination with the traditional methods.
| Analyte | Method | Integral regions/element | Integral regions (ppm) | Calculation | DS |
|---|---|---|---|---|---|
| 2 | DS1H | COCH2(CH2)13CH3vs. AGU region | 0.5–1.9 vs. 2.5–5.5 | ((CH2)/29)/(AGU/7) | 1.65 |
| 2 | DSEA | Calculated from C-content | — | — | 1.62 |
| 3 | DS1H | ArH vs. AGU regions | 8.5–7.5 vs. 2.7–6.0 | 3-((ArH)/4)/(AGU/7) | 1.67 |
| 3 | DS13C | COAr region vs. COAlk | 162.7–164.4 vs. 172.0–173.5 | COAlk × 3/(COAr + COAlk) | 1.67 |
| 3 | DSEA | Calculated from C-content | — | — | 1.53 |
| 3 | DSEA | Calculated from N-content | — | — | 1.75 |
| 4 | DS31P | e-HNDI-TMDP vs. alkoxy-TMDP | 151.5–152.5 vs. 151.0–142.0 | See eqn (1) and (2) | 1.67 |
| 2 | DSmean | Mean analysis value | — | — | 1.65 |
As ILs, such as [amim]Cl, are predicted to be potential environmentally benign media for cellulose chemistry, based upon their recyclability and low vapour pressure, [amim]Cl was chosen as media for this reaction. A typical workup procedure for cellulose reactions, from ILs, involves stirring with a solvent such as water to quench the reactants and precipitate the product. Additionally, heating with the ‘quenching’ solvent is often required to remove traces of IL from the product. This procedure will typically preserve the yield and polydispersity of the product by precipitating most of the functionalised polymer provided the functionality is sufficiently hydrophobic, as is the case with most acyl species. This will also co-precipitate a large portion of the quenched reactants, with the desired product. These co-precipitated species are often impossible to resolve using traditional 1H NMR analysis, preventing rapid, consistent and reliable optimisation of these reactions.
The crude decanoylated cellulose product 5 was subjected to DOSY NMR to attempt to resolve any small organic species, including the decanoic acid quenching by-product, from the functionalized biopolymer. Inspection of the SCORE analysis results (Fig. 4), for a 2-component system, clearly shows resolution of the decanoic acid quenching by-product from the desired product 5. The SCORE analysis estimates this impurity at 20.2% mol/mol (of hydrogen) which, when converted into % w/w, gives a value of 18.6% w/w. From the extracted product 5 spectra (Fig. 4), no accurate DS determination is possible.
![]() | ||
| Fig. 4 SCORE analysis of the crude decanoyl cellulose 5 DOSY dataset, for a 2-component system. The separate spectra are from the two major components with widely differing diffusion rates, indicating low and high molecular weight species. | ||
When the same mixture was phosphitylated with 2-Cl-TMDP and analysed by 31P NMR, according to the above method, accurate values of 25.2% w/w for decanoic acid impurity and 2.29 for DS31P of the impure product were rapidly determined (see ESI†). A proper comparison between the % w/w impurity values from the 31P and DOSY experiments is not meaningful as the DOSY procedure is unable to produce quantitative data. When the DOSY spectrum for the palmitoylated cellulose sample 2 was collected, it was not possible to process the data, using the SCORE algorithm, with the purpose of obtaining any meaningful results that would allow for quantification of the palmitic acid impurity. This was due to the very low quantity of acid by-product impurity present in the sample. On the other hand, however, the 31P NMR analysis method allowed for determination of this impurity to be 0.94 ± 0.078% w/w. The significant advantage of DOSY over 1H NMR or 31P NMR seems to be only in a qualitative discrimination between two abundant components in a mixture, where one component is a high molecular weight polymer. This lack of quantitivity with DOSY is mostly dictated by differences in T1 and T2 relaxation that occur during the DOSY pulse sequence. T1 in particular can be greatly affected by the molecular weight of the resonating species.
This 31P-based analysis method limits the optimisation of these reactions to a aqueous quench, drying of the crude product and one in situ derivatization and NMR analysis procedure, to simultaneously determine DS and purity of the product. This is unrivalled by any other NMR procedure, such as DOSY or 1H NMR, which must serve as complementary qualitative techniques, due to long collection times and poor resolution, in the analysis of these crude reaction mixtures.
Footnote |
| † Electronic supplementary information (ESI) available: Python and Java scripts for determination of DS values, including spectroscopic data for the cellulose ester products, for a visual measure of sample quality. See DOI: 10.1039/c0ay00336k |
| This journal is © The Royal Society of Chemistry 2010 |