Spray coating of cellulose aerogel particles in a miniaturized spouted bed

Aim of this work is to apply protective and homogeneous shellac coating layers on the surface of hydrophilic open-pore cellulose aerogel particles with low densities ≤ 0.1 g/cm3 and high specific surface areas in the range of ~ 400–450 m2/g while keeping the aerogels’ microstructure intact during processing. For this purpose, an innovative miniaturized spouted bed setup was used. Successful process settings for application of enclosed films on aerogel surfaces without intrusion of coating material into the pores were determined. Precise control of coating layer thickness in the range of 10–50 µm was achieved due to variation of coating solution amount without agglomeration and clogging events occurring during processing. Comparison of bulk densities and specific surface areas before and after coating proved the intactness of the porous structure. Coating of particles loaded with vanillin led to controlled release, enhancing release half-life times from 20 to 1600 min. Overall, a successful strategy for coating of organic low-density aerogels was developed.


Introduction
Aerogels are solid materials with a unique combination of properties such as high porosities, low densities and high specific surface areas (S V ). Their synthesis usually involves gelation, solvent exchange and supercritical drying steps (García-González et al. 2011). Besides many other applications, aerogels are suitable as tailorable carrier matrices for flavors and drugs, since they enable the targeted release of the Supplementary Information The online version contains supplementary material available at https://doi.org/10.1007/ s10570-021-04032-0. components from the porous structure and offer high loading capacities (Smirnova and Gurikov 2017;Bian et al. 2020). Food-and pharma sectors demand properties such as biodegradability and biocompatibility for in vivo use, which are provided by biopolymeric materials like starch-, alginate-or cellulosebased aerogels which can e.g. be loaded with various active pharmaceutical ingredients (Ulker and Erkey 2014;Veres et al. 2018;Pantić et al. 2016;Mehling et al. 2009;Zeng et al. 2020). The particular suitability of cellulose aerogel matrices as drug delivery systems is rooted in the possibility to induce controlled and hierarchical porosities, their relatively high stability when brought into aqueous environment, high specific surface areas, and the flexibility of macroscopic shape control, which allows the production of spherical particles with defined sizes (Budtova 2019;Volkert et al. 2009;Liebner et al. 2009). Besides, cellulose is a cheap and abundant raw material.
However, some challenges for in vivo use of biopolymer aerogels as carriers arise if a long-term release of the active agent is targeted: due to the hydrophilic surface and the open pore structure of aerogels the penetration of moisture or reactive molecules into the inner part of the materials may cause damage of the pores during storage and lead to aging of encapsulated material. Moreover, the pores may undergo fast collapse when the materials are exposed to an aqueous environment. Finally, depending on the loading process and interaction between loading substance and aerogel surface, active substances may be loaded in an amorphous state in form of a thin layer into the matrix (García-González et al. 2011). The given points may result in a poor kinetic control and contribute to a burst release of substances from the inner part, which is mainly controlled by the collapse of the 3D-structure instead of desired slower diffusion processes (Ulker and Erkey 2014).
One possibility to address these issues is a hydrophobic post-modification of aerogel inner surfaces, which prevents or considerably slows down the penetration of water into the pores on the one hand, but changes the original surface reactivity on the other hand. Application of a protective coating to the aerogel outer surfaces is an alternative approach, which leaves the inner part unmodified and provides a physical barrier as a protection of the inner porous structure and encapsulants. In addition, external coatings provide the possibility to achieve enhanced control over the release profile by tailoring the coating properties, like layer thickness, pH-and temperature sensitivity. Therefore, the possibility to coat cellulose aerogel particles with protective polymer layers was investigated in this work.
Application of homogeneous coating to the surface of particles is usually carried out in spray coating processes, where moving particles are brought in contact with dispersed coating solutions. Spray coating is a complex process involving many simultaneous steps and process parameters to be optimized (Jones 1985). Specifically, the light weight and low density of aerogel particles have to be considered, since they result in a non-classical fluidization behavior, which may not be described adequately by classical Geldart-classification (Akgün and Erkey 2019). The open porous nature of aerogels may allow significant penetration of the solvent into the aerogel core, resulting in changes or destruction of the porous network, due to the high capillary forces, which occur during subsequent evaporative drying. To prevent this, solidification of the polymer layer has to be faster than penetration of the solution into the pores (Alnaief et al. 2012).
Spouted bed systems are promising setups in this regard since they enable extreme fast heat and mass transfer between gas and particle phases. Aerogel particle fluidization dynamics and spray coating in spouted bed systems have already been described in work, which focused on simulation and mechanical aspects of the process, determining the important process parameters to keep particle macroscopic shape intact during processing (Antonyuk et al. , 2013. Besides particle intactness, two additional requirements can be formulated for spray coating of aerogels in spouted bed systems: (1) The coating materials should be applied in a homogenous manner and defined layer thickness. (2) The inner pore structure has to be preserved to keep the aerogel functionality intact. Especially the second point is not easy to quantify via characterization methods. The shrinkage behavior of silica and starch/alginate hybrid aerogels during coating processes with an aqueous suspension was investigated by Antonyuk et al. (2015) and Alnaief et al. (2012). They also several measures for the preservation of the porous structure from shrinkage were suggested. Since the direct coating with aqueous Eudragit Ò solution led to pore collapse, a double layer approach was necessary, wherein the particles were firstly coated with a melt, before being exposed to water born coating material (Alnaief et al. 2012). In this regard, measurements of specific surface area can provide a sensitive assessment of aerogel quality under varying operating conditions (Subrahmanyam et al. 2015). Goslinska et al. (2019) found a direct correlation between the decrease of specific surface area and amount of applied coating material, showing the intactness of the pores and quality of the coating layer. In the same work, ethanolic shellac solutions were used to coat micrometer-sized whey protein aerogel particles with an initial surface area of 207 m 2 g -1 , firstly creating coated, fully edible aerogels with intact porosity.
Shellac is a natural resin and is used for instance in commercially available modified atmosphere packaging systems (Sandhya 2010). The use of edible shellac coatings as barrier system has been documented since the beginning of the last century, e.g. for pharmaceuticals, confections, fruits and vegetables. A combination of hydrophobicity and biodegradability makes shellac a promising renewable material in formulations with cellulose. In addition, cellulose/shellac systems fulfill the regulations of American food and drug administration legislations on components for food contact materials (Hult et al. 2010).
In this work we aim to investigate the possibility to coat cellulose aerogels with shellac to protect the cellulose matrix from moisture and to enable an (edible) product suitable for long term release of active substances in aqueous surrounding.
Our strategy towards shellac coated cellulose aerogel particles involves the following steps ( Fig. 1): (1) Synthesis of highly spherical cellulose aerogel particles with different particle sizes, densities and narrow size distributions. Sphericity of particles is an important factor with respect to the target application, since it brings benefits in terms of flowability, homogeneous heat-and mass transfer as well as for controlled release of encapsulants (Lee et al. 2013). (2) Spray coating of cellulose aerogel particles with ethanolic shellac solutions in a miniaturized spouted bed setup under variation of process parameters like shellac content in coating solution and the ratio of a spray solution to the bed mass in order to determine successful process conditions. (3) Characterization of coating properties, coating-pore interaction and aerogel microstructure.
(4) Quantification of drug delivery properties by determination of release kinetics from particles loaded in supercritical carbon dioxide.

Preparation of coagulation bath solutions
Different coagulation bath solutions consisting of diluted (33 wt.%) aqueous sulfuric acid (H 2 SO 4 , Merck Millipore, Emsure Ò ), and a mixture of ethyl acetate (EtAc Carl Roth GmbH & Co. KG), trichloromethane (CHCl 3 , Carl Roth GmbH & Co. KG) and acetic acid (CH 3 COOH, Carl Roth GmbH & Co. KG) were used to induce the coagulation of the cellulose solutions. EtAc and CHCl 3 ratios were kept constant at equal mass ratios, CH 3 COOH ratio was 10 wt.% regarding the overall mass.

Preparation of cellulose aerogel particles
A quantity of 10 or 30 g microcrystalline cellulose type II (JRS Pharma GmbH & Co. KG, Vivapure Ò , 101) was dispersed in 200 g demineralized water and left for swelling at 5°C for at least 30 min. 38 g of NaOH (Merck KGaA) were dissolved in 232 or 252 g demineralized water and cooled to -9°C. The swollen cellulose/water dispersion was added to the NaOH solution and stirred (500 rpm) at -9°C for two hours in a double jacket glass reactor with a volume of 500 mL to achieve a homogeneous cellulose solution with a cellulose content of 2 or 6 wt.%. The solutions were stored at 5°C for a maximum of 24 h and used for further processing. Cellulose particles with particle diameter d particle \ 2 mm were produced via jet cutting, following the procedure described by Schroeter et al. (2021). For production of particles with d particle [ 2 mm, dripping method was applied. Cellulose solutions were extruded from tubes and syringes (diameters = 0.4-4.1 mm) into the coagulation bath solutions (3 L, 20°C). The distance between the syringe tip and the coagulation bath was kept small (7 cm) and coagulation solutions where gently stirred via magnetic stirrer during processing. Cellulose solutions were kept at a temperature of 5°C in order to avoid thermal gelation of cellulose during the runs prior to getting in contact with the coagulation bath. Two compositions of the coagulation bath were employed: with aqueous sulfuric acid and with mixtures of EtAc/CHCl 3 /CH 3 COOH. A maximum of 1 L cellulose solution per 3 L of coagulation solution was used in each run. Hydrogel particles were collected with a mesh. Hydrogel particles produced in the aqueous sulfuric acid bath were washed with demineralized water until a pH value of 2 was reached in the water phase, while hydrogel particles coagulated in EtAc/CHCl 3 /CH 3 COOH were directly transferred to the solvent exchange step. Solvent exchange was carried out by immersion of the resulting hydrogel particles in 99.9 wt.% ethanol until a minimum final concentration of 98.0 wt.% ethanol inside the hydrogel particles was reached (controlled by density measurements, Anton Paar, DMA 4500 M). Approximately 7 L EtOH per liter hydrogel were used. After the solvent exchange, the particles were sealed into a filter paper bag and placed into a high pressure autoclave with an overall volume of 3.9 L for supercritical drying with CO 2 . The supercritical drying was performed at a temperature of 60°C and pressure of 120 bar. A continuous flow of CO 2 (flow rate = 120-140 g/min) was set through the autoclave until complete extraction of ethanol was achieved after 4 h. The dried particles were collected after slow depressurization (1 bar/min) of the autoclave. Dried particles were stored in sealed tubes (25 mL, Sarstedt) in a desiccator over calcined silica gel prior to analysis.
Specific surface area and pore size distribution Low temperature N 2 adsorption-desorption analysis was used to investigate microstructural properties of the aerogels (Nova 3000e Surface Area Analyzer, Quantachrome Instruments, Boynton Beach, USA). An overall sample mass of 20-30 mg was used for each analysis, which corresponds to 6-20 intact particles per measurement, depending on the overall particle size. The specific surface area was determined using the BET (Brunauer-Emmett-Teller) method (single determinations). The standard error of 5.5% for the individual BET-results was estimated from the average error of four double determinations. The pore volume and mean pore diameter were estimated by the BJH (Barrett-Joyner-Halendia) method and applying approximation of cylindrical pores, respectively. All

Bulk density
A container with a volume of 12.5 cm 3 was completely filled with dry aerogel particles. The full as well as the empty container were weighted using an analytical scale (Excellence Plus XP, Mettler Toledo) and the bulk density (q bulk ) was calculated as a ratio of the particles' mass and the known volume of the container. All measurements were carried out in triplicate.

Preparation of shellac solutions
Coating solutions with three different shellac concentrations c shellac were prepared as follows: 37.5-75.0 g of shellac flakes (SSB56 Pharma FL) were mixed continuously with 175.0-212.5 g of pure ethanol (Table 1). After dissolution of the shellac flakes, 1.2 g of a red or blue food dye (Ruth GmbH, strawberry coloured, patent blue V) was added to the solution to enable the visual evaluation of the coating procedure during the process.

Coating of cellulose aerogel particles in a miniaturized spouted bed
A miniaturized conical-cylindrical spouted bed was used for the coating process of the aerogel particles with the shellac-ethanol solution. The process chamber, the expansion zone, and the lower part of the relaxation zone of the spouted bed are conical and the upper part of the relaxation zone is cylindrical (Fig. 2).
A charge of 4.0-18 g (depending on the bulk densities for particles of different composition and size) of aerogel particles were introduced into the process chamber and the fluidization was started. The process conditions for all experiments are shown in Table 2. The particles were fluidized with heated air at gas volume flow rates _ V fluidisation between 1.2 and 5.1 m 3 /h and at a temperature of 80°C at underpressure conditions. The high gas velocities in the annular gap result in a turbulent jet in the middle of the process chamber taking along the particles. The gasparticle jet passes the conical expansion zone and enters the cylindrical relaxation zone, and as a result of the decrease of gas velocity due to the increasing diameter of the chamber, a fountain is formed. The particles drift toward the walls of the apparatus and slide back into the process chamber due to gravity and are then entrained by the jet again. Thus, the typical, circulating particle movement inside a spouted bed is formed (Goslinska et al. 2019;Rothberg et al. 2020).
The shellac/ethanol solutions (c shellac = 15-30 wt.%) were filled into a syringe and different quantities of the spraying solution (8-50 mL) were injected into the apparatus during the experiments via a syringe pump (Legato 270, World Precision Instruments Germany GmbH, Friedberg, Germany) with a spraying rate of 0.5 mL/min through a two-fluid nozzle (Model 970/3 S1, Düsen-Schlick GmbH, Germany, nozzle tip size: 0.5 mm). The nozzle was installed in a bottom-spray configuration. In order to obtain different ratios r v/m between the overall particle mass and amount of coating solution, the mass of particles m bed was kept constant at given reaction conditions and the quantity of the injected solution v spray was varied. To avoid the risk of elutriation of the particles the atomization air was kept at low flow rates ( _ V atomization = 0.5-0.9 m 3 /h). At these spray air volume flows the formation of small droplets was still ensured. The droplets are spreading on the particle surfaces, ethanol evaporates and a shellac layer solidifies on the particle surfaces. Before the completion of the experiments, an additional 10 min post drying step (spouting without spraying) was applied to ensure the complete evaporation of the solvent.

3D laser scanning confocal microscopy
The surface roughness is described by the arithmetical mean height (Sa) of the surface and was measured via a 3D laser scanning confocal microscope VK-X160K (Keyence). Six measurements at 109 and six measurements at 1009 magnification were performed for each sample. The taken images were evaluated using the surface area tool of the MultiFileAnalyzer software (Keyence). To smooth the curvature caused by the spherical shape of the particles as well as signal noise and thus avoid aberrations in the roughness calculations, a sphere correction tool and filters were used in the software that neglect depressions and heights below 8 and above 80 lm. Through image analysis (ImageJ) of the laser microscope pictures the  Particle size and composition Spray coating process parameters coating layer thickness was determined and afterwards compared with the image analysis results of the scanning electron microscope images. For this purpose, respectively three coated particles were cut in half and six pictures at 209 and/or 1009 magnifications of the two opposing sides of one half were taken. No differences resulting from the magnification lenses were noticed in the evaluation of the pictures.

Scanning electron microscopy
The surface properties and inner structure of the aerogels as well as the layer thickness of the shellac coating were characterized via scanning electron microscopy (SEM, Zeiss Supra VP55, Jena, Germany). The intact and cut particles were sputtered with a thin layer of gold (Sputter Coater SCD 050, BAL-TEC) before analysis was started. The measurements were carried out under high vacuum at an accelerating voltage of 2-5 kV.

Image analysis of particles
Size and shape of dry aerogel particles were determined with a Camsizer XT system (Retsch Technology) in free fall mode. All characterizations were carried out according to DIN 66141. As particle diameter d particle the calculated Sauter mean diameter was adopted: where x min is the particle size below which there are no particles in a given size distribution, x max is the particle size above which there are no particles in a given size distribution, q 3 (x) is the volume (mass) distribution density, M -1,3 is the complete first moment of a q 3 (x)sample distribution. Particles' sphericity SPH is given according to ISO 9276-6 as: where P is the measured perimeter/circumference of a particle projection and A is the measured area covered by a particle projection. For an ideal sphere SPH is equal to 1 and for real samples lower than 1. The polydispersity of the samples is given as the Span value based on volume: with xQ 3 = Q 3 values, where a given particle diameter x is reached, based on volume.

Supercritical impregnation
Vanillin was used as an UV-Vis active and sublimable flavor encapsulant, which can be supercritically impregnated into the porous aerogel structure due to its good solubility in supercritical CO 2 (Rojas-Á vila et al. 2016). For supercritical impregnation, aerogel particles were sealed in filter bags and transferred to an autoclave with an overall volume of 250 mL, vanillin powder was placed in an open metal cup at the bottom of the autoclave. Particles with an initial diameter of d particle = 3.29 mm and a cellulose content of c cell-= 6 wt.% were used. A weight ratio of aerogel:vanillin = 1.0:1.5 was used. Impregnation was carried out at a pressure of 125 bar and a temperature of 60°C for 16 h. The impregnated particles were collected

Release experiments
Impregnated aerogel particles were placed in a beaker filled with 500 g demineralized water kept at a temperature of 30°C and stirred at 200 rpm with an overhead stirrer. Six to eight particles were used in each experiment. Samples of the water phase were taken periodically (2 mL each) and transferred to cuvettes (Hellma Analytics, 117-QS). The vanillin concentration was determined via UV-Vis-spectroscopy (Thermofisher Evolution 300) at a wavelength of 230 nm. Sample solutions were returned to the beaker after determination of vanillin concentration.

Results and discussion
Production and properties of cellulose aerogel particles Spherical cellulose aerogel particles with different densities were obtained by variation of coagulation baths compositions. Aqueous sulfuric acid based coagulation has been reported to be a successful and green way towards particles at higher cellulose content with high SPH and was therefore also used in this work for production for particles with c cell = 6 wt.% (Schroeter et al. 2021). However, this method is not suitable at low cellulose contents of 1-3 wt.%. For the production of particles with c cell = 2 wt.% the method reported by Bian et al. was adapted and a coagulation bath consisting of EtAc/CHCl 3 /CH 3 COOH was applied. At c cell = 6 wt.%, particles with five different diameters in the size range of 1.01-4.37 mm were produced, and at c cell = 2 wt.% one particle size of 3.18 mm was set. Narrow particle size distributions (SPAN = 0.12-0.58) and highly spherical shapes (SPH = 0.91-0.96) were obtained for particles with d particle [ 2.00 mm produced by dripping method.
Slightly lower values of SPH = 0.86 were obtained for particles with c cell = 2 wt.%, which can be attributed to partial breakage of some particles during shape analysis in the Camsizer XT system. All aerogels showed high specific surface areas (428 ± 35 m 2 /g) with no notable trend with respect to the coagulation conditions. Obtained surface areas are in the higher range for cellulose II aerogels and in good agreement with the previously reported results (Budtova 2019;Schroeter et al. 2021). In case of aerogels with c cell = 2 wt.% the achieved surface area of 405 m 2 /g exceeds the values reported by Bian et al. (2020) by 39-61%, who used freeze drying in the final processing step. This difference can be explained by the difference of applied drying techniques, since supercritical drying preserves the mesopores, while they are mainly lost during the freeze-drying leading to macroporous materials with significantly lower specific surface areas (Baudron et al. 2019).
Microstructure, density and surface roughness of the cellulose aerogel particles were found to be significantly dependent on the processing route and cellulose content (Table 3).
While differences in the bulk densities can be related to the different overall cellulose contents, structural and textural properties are most likely affected by the coagulation methods. The results of characterization of both the pore structure and the particle surface via SEM support the results given in Table 3 and allow more detailed insights into the structural properties of the materials. Aerogels with a cellulose content of 2 wt.% show an open macroporous and wavy surface (Fig. 4a, b). In addition, the porous network contains macro channels (Fig. 4c) with diameters of approx. 200-600 nm which reach up to the particles surface. In contrast, particles with cellulose content of 6 wt.% show a homogenous and smooth surface which contains no macropores (Fig. 4d, e), whereas homogeneous distributed mesopores with diameters B 40 nm can be identified. The inner porous network appears denser compared to particles with lower cellulose content (Fig. 4f). In addition, a thin, dense layer in form of a skin is formed on the surface of particles coagulated in sulfuric acid, as also reported in previous works (Schroeter et al. 2021), while such layer could not be observed for particles coagulated in EtAc/CHCl 3 /CH 3 COOH.
Even though the mechanisms of acidic induced coagulation of sodium hydroxide based cellulose solutions are not completely understood up to now, differences of microstructure and surface properties may be attributed to solvent effects occurring during the hydrogel formation. While in aqueous sulfuric acid baths coagulation is controlled by one process (diffusion of water/acid mixture into the inner part of the particle), solvent mixing/demixing effects occurring in presence of acetic acid/ethyl acetate/chloroform mixtures may lead to solvent gradient influenced network formation which contributes to the formation of less dense pore structures and less homogenous surfaces.
Summarizing we can say, that we adapted already existing procedures for production of cellulose particles at different cellulose contents in this work, showing the differences in the properties and enhancing the accessible particle sizes. In case of coagulation in EtAc/CHCl 3 /CH 3 COOH-bath (c cell = 2 wt.%) we could show, that high surface areas of 400 m 2 /g are achievable due to supercritical drying. Therefore, we have demonstrated the possibility to induce variable pore structures and surface properties in cellulose aerogel particles on the one hand while keeping the specific area on a comparable level on the other hand.
In the next step the coating of the produced aerogel particles was investigated.

Coating layer properties
One of the most important characteristics of coating quality is the uniformity of coating material on the particle surface, since non-homogeneous coating results in incomplete protection of the inner porous part and may lead to decreased release control in case of loaded particles. Image analysis via laser microscopy and SEM showed, that the porous structure was completely covered by the coating material and that no cracks and fissures were obtained when working with aerogels with c cell = 6 wt.% and a shellac content of c shellac = 20-30 wt.% in the coating solution (Online Resource 1, Figs. S1-S8). A slight increase in surface Coagulation bath (-) Sa (lm)* d pore (nm) v pore (cc/g) q bulk (g/cm 3 )* 2.0 EtAc/CHCl 3 /CH 3 COOH 3.89 ± 1.81 8.65 ± 0.05 2.92 ± 0.20 4.10 9 10 -2 ± 1.56 9 10 -3 6.0 Aqueous H 2 SO 4 1.72 ± 0.08 13.7 ± 0.07 3.99 ± 0.28 1.07 9 10 -1 ± 9.04 9 10 -3 d pore is the average pore diameter, v pore is the pore volume *At c cell = 6 wt.%: results for all different particle sizes used to calculate the average and standard deviation roughness from Sa = 1.7 ± 0.1 lm to Sa = 3.6 ± 0.8 lm was determined for coated particles of different particle sizes compared to the non-coated species (r v/ m = 2.0 mL/g, c shellac = 20 wt.%), whereas no specific influence of particle size was detected. This increase in surface roughness can be explained by the rapid drying of the solution after coming in contact with the particles, which causes some individual droplets to solidify on the surface before they spread. In contrast, the surface roughness at lower cellulose content does not change significantly after application of coating material (r v/m = 2.0 mL/g, c shellac-= 20 wt.%): on the one hand, the solidified shellac encloses gaps in the surface, and on the other hand the coating material is applied in a non-uniform manner, as shown by microscopic pictures (Online Resource 1, Figs. S9-S11). In addition, parts of the porous structure can be identified which remain uncovered, therefore particles with c cell = 2 wt.% were only partly coated. The results indicate that initially smoother and non-macroporous aerogel surfaces are important for the formation of uniform coating layers. Therefore, further optimization and quantification of process conditions was carried out at c cell = 6 wt.%. The influence of the shellac concentration on the surface morphology was investigated at two different settings (c shellac = 20 wt.% and 30 wt.%) and four different ratios between coating solution and bed mass (r v/m = 1.0-3.5 mL/g).
The surface roughness increases linearly (R 2 = 0.985-0.997) with the amount of applied coating solution (Fig. 5, left), which is a result of the increasing solidification of droplets on particles surface (Fig. 5, right). In addition, the surface roughness is more effected at higher shellac concentrations. This can be explained by the higher viscosity of the spray solution at c shellac = 30 wt.% which leads to slower spreading of the droplets after hitting the particle surface. In consequence, the liquid evaporates before the droplet has spread completely on the particle surface, and the solid material left behind forms an irregular surface. Due to improved homogeneity of the coating obtained at c shellac = 20 wt.%, the given shellac concentration was identified as preferable one. Variation of coating layer thicknesses was therefore achieved by changing r v/m while leaving all other process parameters constant (c cell = 6 wt.%, c shellac = 20 wt.%). Different layer thicknesses were set for particles loaded with vanillin (d particle-= 3.06 mm) and non-loaded particles (d particle-= 4.37 mm) as determined from images analysis of cut open particles (Fig. 6). A clear separation is observed at the interface between the coating material and the inner porous part, which indicates that no coating material and solvent has penetrated into the pore structure during processing (Fig. 6a, b).
The layer thickness increases linearly (R 2 = 0.975-0.997) as function of applied coating  (Fig. 6c). This corresponds well to the difference in the overall particle volume (factor 2.1), which was calculated from the particle diameter. In order to directly compare the influence of r v/m on the total amount of applied coating per particle for both particle sizes, the coating volume per particle was calculated from the difference of the overall particle volume prior coating and that volume with added layer thickness. All calculations were carried out under the assumption of a perfectly spherical shape for all particles. A similar influence of r v/m on the layer volume per particle was determined for both particle species (Fig. 6d), which shows that the spray coating process runs under steady state conditions in the investigated range of particle sizes and densities.

Particle properties
In addition to changes of surface properties, the coating process can lead to changes in particles size and shape, influencing the bulk density q bulk , sphericity SPH and the particle diameter d particle . Changes of the latter in terms of shrinkage is a measure for the pore collapse caused by the coating process. A large value of the shrinkage corresponds to a loss of aerogel specific properties. A comparative study of given properties was carried out for of non-loaded aerogel particles at constant cellulose content (c cell = 6 wt.%) and constant amount of coating solution (r v/m = 2.0).
Coating of particles results in an increase of the bulk density (? 22-45%), whereas the change correlates mainly with the initial particle diameter and results in relatively low overall values of q bulk-B 0.18 g/cm 3 for coated particles (Fig. 7). Only small changes of SPH and d particle were observed, Fig. 6 Exemplary image analysis of a coated particle (c cell-= 6 wt.%, d particle = 4.37 mm, c Shellac = 20 wt.%, r v/m = 1): a scanning electron microscope; b laser microscope. Influence of coating amount on layer thickness; c layer volume per particle. d Straight lines correspond to linear fits, error bars to the standard deviation of image analysis, dashed lines to the 95% confidence interval with exception of particles with the smallest initial values. The results clearly show, that no significant shrinkage occurred during processing, in most cases even slightly higher particle diameters were determined for coated particles. Based on values of layer thicknesses determined by microscopy, an increase in particle diameters in the range of approx. 0.5-4.0% would have to be expected in case of complete lack of aerogel-shrinkage. A direct quantitative comparison of results regarding changes of particle size obtained via microscopy and particle shape analysis with the Camsizer XT system is however not possible, since the results obtained via Camsizer XT are also influenced by interactions of obtained particle size distributions with particles sphericities, as reported previously (Schroeter et al. 2021). For particles at lower cellulose content of 2 wt.% and initial particle diameter of 3.18 mm a significant linear shrinkage of 14.6% was found under the same process conditions. Obtained shrinkage is however still relatively low when compared to values reported in the literature: spray coating of silica aerogels with waterborne solutions led to an overall linear of shrinkage of up to 75%, and exposure of alginate/starch based aerogels to air humidity to volumetric shrinkage up to approx. 85% (Alnaief et al. 2012;Antonyuk et al. 2015). Therefore we can conclude, that only partial pore collapse occurred due to ethanol entering the comparatively open surface in case of particles with c cell = 2wt.%. Changes of SPH (-4.2) were in the same range as obtained for particles at higher cellulose content, while the increase in q bulk was comparatively high (? 55%). Overall, almost no negative influence of the coating process on the aerogels's characteristics was found, especially for particles with high initial sphericity [ 0.9 and a cellulose content of c cell = 6wt.%.
Additionally, the influence of the amount of spray solution on particle properties was investigated for particles previously impregnated with vanillin and at a constant initial particle size (d particle = 3.06 mm) in the range of r v/m = 1.0-4.0 mL/g. Only small changes of overall particle sizes (average of -3.4 ? -4.4%) were found at all values of r v/m . The bulk density increases as a linear function of r v/m (R 2 = 0.983) as well as of the calculated layer volume per particle (R 2 = 0.982) (Fig. 8). The particle sphericity was not affected significantly.
If the bulk density of coated and loaded particles is considered as the sum of the contributions of the pure aerogel, vanillin introduced by the impregnation and the coating material, the contributions of the individual components can be calculated for different values of r v/m (Eq. 4). q bulk;total ¼ q bulk;aerogel þ q bulk;vanillin þ q bulk;coating ð4Þ with q bulk, total = total bulk density of impregnated and coated particles, q bulk, aerogel = bulk density of pure aerogel particles in relation to total bulk density, q bulk, vanillin = difference in the bulk density of vanillin loaded particles and pure aerogel particles in relation to the total bulk density, q bulk, coating = difference in the bulk density of coated and non-coated, loaded particles in relation to the total bulk density.
The contribution of the coating material to the total bulk density ranges from 11-24% and correlates linearly with r v/m (R 2 0.981) (Table 4). Consequently, Fig. 7 Influence of spray coating on bulk density and d particle .
Process conditions: c cell = 6wt.%, r v/m = 2.0, c shellac = 20wt.%. Straight lines correspond to the linear fitting, error bars to the standard deviation of bulk density measurements. Bracketed lines symbolize the size change of a given particle species, values in % correspond to the linear shrinkage decrease (-), or increase (?) in particle size after coating relative to the noncoated species. Numbers without unit correspond to the SPH value of the according particle species bulk density can be used as a direct measure of the overall weight increase of particles during to the coating process.

Specific surface area
Analysis of cutted aerogel particles via microscopic methods suggested that no change of pore structure and no penetration of coating material into the pores occurred during the processing. At the same time, quantification of the influence of the spray coating process on the microstructure is possible by determination of the specific surface area before and after coating. Since determination of the specific surface area via BET method is based on gravimetric analysis of adsorbed nitrogen on the pores' surfaces, the question arises whether nitrogen can diffuse through the applied shellac layers of whole particles during the measurement process. In order to determine the influence of the outer coating layer, BET results from whole and cut open coated particles were compared at different ratios of r v/m . The same results for whole and cut particles were obtained for all samples (Online Resource 1, Fig. S12), therefore the shellac layer had no effect on the measurement results at all ratios of applied coating material.
A systematic study on the effect of r v/m and c shellac was carried out at constant c cell = 6 wt.%. Particles showed high surface areas after the coating process which correlated linearly (R non-loaded 2 0.939, R loaded 2 0.969) with the amount of applied shellac solution (Fig. 9, left). The loss in specific surface area of loaded particles compared to non-loaded ones is independent of r v/m (average of 50.0 ± 0.8%) and can be related to the additional contribution to the overall weight of the sample from vanillin. Furthermore, a general linear relation between the increase in bulk density due to coating and the loss of specific surface area after coating process shows, that changes of the specific surface area result from the weight of coating material (Fig. 9, right). In contrast, loss of specific surface area was found to be 19-32% higher compared to increase in bulk density for particles with lower cellulose content (c cell = 2 wt.%), which results partly from higher shrinkage under these conditions. Overall, we Fig. 8 Bulk density and sphericity auf vanillin loaded particles as function of r v/m , at constant initial d particle (4.37 mm) and c shellac (20 wt.%) (left) and layer volume per particle as function of the bulk density (right). Error bars correspond to the standard deviation of bulk density measurements, and determination of layer thickness via microscopy. Straight lines correspond to linear fitting, dashed lines to the 95% confidence interval can conclude that the differences in initial surface properties of the particles at different cellulose content are decisive for the quality of the process: The smoother, but open porous surface in form of a skin at c cell = 6 wt.% led to homogeneous distribution of the coating material on the particles surface without intrusion of coating solution into the pore structure, resulting in a complete preservation of the pore structure.

Release kinetics
In order to describe the vanillin release from different coated and non-coated aerogel particles, the release data was fitted via a flexible Weibull model approach.
The Weibull model has proofed to be suitable for highly correlated non-linear kinetic description of chemical reactions (Schroeter et al. 2019;Manso et al. 2001;Kokkinidou et al. 2014) and follows with c t being the concentration at time t, c i is the initial concentration, a is the scale constant and b is the shape constant. The quantity a -1 is analogous to a rate constant. The shape constant b describes mathematically the shape of the reaction profile and is temperature independent, if the reaction mechanism is not related to temperature (Kokkinidou et al. 2014). If b = 1 equation Eq. (5) corresponds to the first-order kinetic model. For description of growth processes the Weibull model can be adapted (Manso et al. 2001). Equation (5) describes the concentration decrease in the aerogel particles. Since we measured the concentration increase in the aqueous phase, Eq. (5) has to be rewritten by setting c i to c t ¼ 1 ð Þ¼c end and calculating the quantity c end À cðtÞ. This yields for the concentration of vanillin concentration in the water phase, cðtÞ: with c end being the vanillin concentration in the water phase at the end of the measurement. For individual measurements c end was set to 1, i.e. individual measurements were normalized to the concentration at the last measured point. In contrary, in Fig. 10 each data set is compared considering the total amount of vanillin, i.e. c end corresponds to the absolute amount of vanillin in the aerogel. Fitting was applied with high correlation at all reaction conditions (Table 5).
Without coating the loaded cellulose particles released 50% of the vanillin in 20.4 min and 95% after 110 min (Fig. 10, left). The release profile of non-coated particles shows, that fast release occurs up to * 60% release, followed by a slower release phase, which indicates that diffusion control of the process becomes increasingly dominant.
Application of coating led to sustained release control in dependence of the coating layer thickness and leads to increasingly homogenous release profiles. The enhancement of release times is also reflected by Fig. 9 Specific surface area in dependence of r v/m for loaded and non-loaded particles at constant c shellac = 20 wt.%, c cell-= 6 wt.% and d particle = 3.3 mm (left). Error bars correspond to the standard error of BET-results. Decrease of specific surface area due to coating as function of increase in bulk density, c cell = 6 wt.%, c shellac = 15-30 wt.%, d particle = 2.0-4.4 mm. Error bars correspond to the standard error of BET-results (yaxis) and standard deviation of bulk density measurements (xaxis). Straight line corresponds to linear fitting, dashed lines to the 95% confidence interval the rate constant a -1 , which decreases by 74% at r v/ m = 1.0 mL/g compared to the release from the uncoated species and correlates linearly (R 2 0.985) with the coating amount (Fig. 10, right). While the half-life times of the release t 1/2 could be determined directly for r v/m B 2 mL/g, half-life times for r v/ m [ 2 mL/g were calculated by linear extrapolation of the raw data under assumption of intactness of the coating over the extrapolated time period. A broad range of total release times is accessible within the investigated range for vanillin, which is reflected by the exponential increase of t 1/2 as function of added shellac solution (R 2 0.999) (Fig. 10, right). Mechanistically, drug release in aerogels depends on several factors: the rate determining step in the release of vanillin from non-coated cellulose aerogel particles can expected to be the transport of the dissolved drug from the pores and particles external surface to the water phase, which tend to be fast, since they depend on the hydration and dissolution of a hydrophilic drug from a hydrophilic aerogel (García-González et al. 2021). In addition to wetting and hydration, overall drug release rates from aerogels can also be influenced by erosion of the gel matrix, whereas biopolymer aerogels swell and collapse in water, which results in the formation of hydrogel like particles (García-González et al. 2021;Kéri et al. 2020). Since the shellac layer is practically insoluble in non-alkaline aqeous media, the rate of vanillin permeation in the given case should be proportional to both: to the extent of the hydration and the thickness of the coating layer, the latter being reflected by the rate constants (Table 5; Fig. 10). Furthermore, coating as well as particles were (optically) intact after the measurements, which indicates that the shellac coating may also inhibit possible erosion of aerogel particles, which are less stable then cellulose based materials. Fig. 10 Release patterns of vanillin from impregnated cellulose aerogel particles (c cell = 6 wt.%, d particle = 3.06 mm) with different coating amounts (c shellac = 20 wt.%) (left). All data is normalized to 100% release. The solid lines correspond to Weibull model fitting. Release rate constants and half-life times t 1/2 in dependence of r v/m (right). Solid lines correspond to data fitting (linear = black, exponential = blue), dashed lines to the 95% confidence interval. Error bars correspond to the error from individual fitting Table 5 Fitting of kinetic parameters and rate constants at different ratios of r v/m, c cell = 6 wt.%, d particle = 3.06 mm, c shellac = 20 wt.% 0.0 0.989 3.2 9 10 -2 1 9 10 -3 0.71 0.03 1.0 0.999 8.2 9 10 -3 4 9 10 -5 1.42 0.02 2.0 0.994 6.1 9 10 -3 1 9 10 -4 1.26 0.04 2.5 0.998 5.4 9 10 -3 4 9 10 -5 2.02 0.03 3.0 0.995 4.7 9 10 -3 5 9 10 -5 2.14 0.06 4.0 0.921 4.0 9 10 -3 2 9 10 -4 1.4 0.1

Conclusions
It was firstly shown that the spouted bed technology is a suitable method for coating of low density, hydrophilic and open porous cellulose aerogel materials with ethanolic shellac solutions. On contrast to silica aerogels, no additional protective layer is required for cellulose aerogels. Suitable routes for production of cellulose aerogels with different densities and high sphericities were identified, including different coagulation strategies to tune aerogel particles' surface and textural properties. Specific surface areas in the range of * 400-450 m 2 /g were obtained using supercritical drying technology. By analysis of specific surface area and bulk densities before and after coating, it was unambiguously demonstrated that the inner porous structure of aerogels remains intact during the coating process. A clear border between coating layer material and pore structure could be identified via microscopic methods. In consequence it can be concluded that the drying of the sprayed solution on the particle surface was faster than its penetration into the pores, which was one of the targets of this work. It was shown that increase of the surface roughness resulted from the fast drying of droplets on the particles surface, with higher changes of roughness being obtained at higher amounts of shellac solution and shellac concentrations. The precise setting of coating layer thickness was possible due to variation of the ratio between sprayed coating solution and number of processed particles. Low initial roughness and homogeneities are decisive for the quality of the coating process: particles with a cellulose content of 6 wt.% and a smooth surface showed complete enclosed coating layers, while surfaces with initial defects in from of gaps and higher roughness at lower cellulose content (2 wt.%) led to incomplete coating layers, with parts of the open porous structure remaining un-covered. Stable fluidization without breakage or agglomeration of the particles was achieved for particles of different densities and a broad range of particle sizes from * 1.1 to 4.4 mm, as reflected by the analytics of the size and shape prior and after coating.
Impregnation of the particles with vanillin in presence of supercritical CO 2 resulted in loaded particles with intact pore structure. It was demonstrated that coating of vanillin-loaded particles with different layer thicknesses allows to control the release of vanillin in aqueous environment in a broad range of release times.
While the shellac coating enables controlled release from the inner porous part on the one hand, it serves also as a protection layer for the pore structure against external influences (e.g. humidity) on the other hand, therefore enhancing the aerogels shelf life. The combination of protection and controlled release enhances the applicability and tuneability of biopolymer aerogels, especially as drug delivery excipients in the pharmaceutical field and as part of digestible or injectable formulations.
While the given showcase demonstrates the applicability of spouted bed spray coating for the creation of biocompatible and edible aerogel systems, further studies are necessary to show the transferability of this method to other hydrophilic biopolymer aerogels (e.g. alginate) with low density and variable surface properties. Thus, in our future work different aspects of shellac coating in regards to its drying kinetics, film formation at different process conditions and coating uniformity on other biopolymer aerogels, will be investigated.
Author contributions The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript.
Funding Open Access funding enabled and organized by Projekt DEAL. The authors wish to acknowledge support for this research received from the German Research Foundation (DFG) under the project SM 82/18-1 and HE 4526/24-1.
Data availability Additional Data (Table S1, Figs. S1-S12) is reported in the electronic supplementary data.

Declaration
Conflict of interest All authors declare that they have no conflict of interest.
Human and animal rights This manuscript does not contain any studies with human participants or animals performed by any of the authors.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.