Abstract
This article aims to highlight the synthesis of foams from a blend of hydroxyurethane of castor oil and epoxy resin. An epoxidized castor oil of 4% oxirane oxygen was first converted to cyclic carbonate of castor oil at 120°C, 1 atm CO2 pressure and then it was reacted with three different aliphatic diamines to yield amine terminated Polyhydroxyurethane (PHU). Foams were prepared in a metal mould from the blend of PHU, epoxy resin, epoxy hardener and polymethylhydrogensiloxane blowing agent which releases hydrogen gas upon reaction with amine. FTIR and 1H NMR of cyclic carbonate of castor oil and PHU of castor oil were done to confirm their chemical structures. Optical microscopy and scanning electron microscopy of foams was done to assess their cellular morphology along with DSC and TGA to evaluate their thermal properties. Both flexible and rigid type of foams were synthesised in this study. Resilience of flexible foams was inspected using a ball rebound test and compression-recovery test while thermal insulation property was checked by measuring thermal conductivity, thermal diffusivity and R-values of rigid foams from heat transfer study using a heat transfer apparatus.
Introduction
Polymer industry was boosted in the 20th century not just with the advent of world war but also with the catapulting emergence of technology. In 1937, it was Otto Bayer et al. 1 who first synthesised and introduced to the world, polyurethane, from polyols and isocyanates in an attempt to replace natural rubber. Since then, polyurethane have been commercially synthesised by a polyaddition reaction between di/poly-isocyanates and diols/polyols. 2 Polyurethane exhibits a wide range of mechanical properties such as flexibility, rigidity or elasticity, good damping properties and also resistance to weathering, abrasion and impact.3–5 Owing to these peculiar characteristics PU applications cover foams (67%) and elastomers (6%), 2 adhesives and sealants (6%),6,7 coatings (3%) 8 and others (18%). 9 Production of polyurethane is done by various isocyanates such as methylene diphenyl diisocyanate (MDI), 10 toluene diisocyanate (TDI), 11 hexamethylene diisocyanate (HDI),12,13 isophorone diisocyanate14–16 and methylene bis-(4-cyclohexyl-isocyanate). 17 However, it is a disappointing fact that most of these isocyanates contribute to certain acute and chronic ailments which is often overlooked especially when it comes to their commercial use in synthesis of polyurethane. Kapp 18 proposed toxicological studies on various isocyanates and reported that under mild or prolonged exposure to these chemicals a person is likely to acquire symptoms or ailment from any of the following: asthma, dermatitis, eczema, skin and respiratory sensitization, pulmonary oedema, coughing and respiratory distress, immune sensitization, infertility, foeticide, eye irritation and more. Moreover, down the decades researchers have assessed the adverse effects of TDI and MDI on human health19–21 and concluded that their prolonged exposure may become the concern of causing tumours in human organs and eventually cancer. Therefore, the European Council in the Commission Regulation for Chemicals deemed these diisocyanates to be Carcinogenic, Mutagenic and Reprotoxic. In the latest edition, COMMISSION REGULATION (EU) 2020/1149 of 3 August 2020, 22 some serious restrictions were imposed on the use of diisocyanates in terms of concentration to be used, training and safety protocols to be followed and marketing with adequate information.
In order to overcome the toxicity related issue in conventional isocyanate-based polyurethane, a new pathway has emerged that completely discard the use of isocyanates designated as Non-isocyanate Polyurethane (NIPU). This pathway has been practised a lot recently as it utilizes a simple polyaddition reaction between a 5-membered cyclic carbonate and an amine to yield polyurethane with a secondary hydroxyl group in the main chain, hence also called as polyhydroxyurethane.23–36 Although amines are readily available substances, synthesis of cyclic carbonates correspond to key aspect of a methodical approach toward PHU. So far in literature synthesis of cyclic carbonates have been carried out from 1,2-diols,37,38 olefins, 39 halohydrins, 40 propargylic alcohols,41,42 urea43–45 and linear carbonates.46–48 However, the most cheap, convenient and preferred route to obtain cyclic carbonates has been attributed to the chemical fixation of CO2 into epoxides.33–36,49–55
Majority of polyurethane applications are covered by flexible foams for mattresses and cushions, semi-rigid and rigid foams for thermal insulation. Production of polyurethane foams is done by making use of physical blowing agents like hydrocarbons (HC’s), hydrofluorocarbons (HFC’s), hydrochlorofluorocarbons (HCFC’s), and chemical blowing agents such as dialkyl carbonates, azo compounds and peroxides. 56 Chemical blowing agents decompose to produce gas that blows the polymer matrix whereas physical blowing agents expand quickly at foaming temperature as a result of vaporization. Unfortunately, during NIPU synthesis no such blowing agent is produced in any step or reaction. Therefore, requirement here is of a chemical species which is non-toxic, stable and which will produce enough gas to blow the polymer matrix to an appreciable volume. Polymethylhydrogensiloxane has been used as a blowing agent in literature for synthesis of foams especially in the system where amines are primary reactants.57–62 It undergoes a reaction with amines to yield hydrogen gas.
Researchers worldwide have reported the synthesis of non-isocyanate polyurethane foams from cyclic carbonates and amines as starting materials. Cornille et al.61,62 synthesized NIPU foam using cyclic carbonates, diamines and polymethylhydrogensiloxane as blowing agent. While Figovsky et al. 63 for the first time patented the composition of a sprayable NIPU foam, Lauth et al. 64 patented NIPU foams from a carbonate blend, hexamethylene diamine and HFC’s as a blowing agent. Also, Grignard et al. 65 reported NIPU foams of low density by utilizing supercritical CO2 as a foaming medium. Glucose-based NIPU foams from linear carbonates66,67 along with sorbitol and lysine based NIPU foams 68 have also been reported recently in literature.
The biggest limitation of NIPU technology is the low reactivity between cyclic carbonate and amine which imparts inferior properties to the polymer such as low molecular weight and poor performance. In order to circumvent these challenges hybrid NIPU materials have been introduced constituting polyhydroxyurethane, epoxy and/or acrylic as major components.54,69–71 Cornille et al. 72 prepared epoxy-urethane polymers of high molecular weight (1,200–3,000 g mol−1) from oligo polyhydroxyurethane. In this article, cyclic carbonate of castor oil and aliphatic diamines have been used to prepare PHU based foams. Due to limitation of low reactivity and poor conversion reaction of CCO, synthesis of foams from PHU of castor oil alone was far from possible. Therefore, commercial epoxy/hardener system was used and it was blended with PHU as a supporting medium. 73 Amine terminated PHU of castor oil was prepared from reaction between cyclic-carbonate of castor oil (CCO) and various diamines such as ethylene diamine, hexamethylene diamine and diethylene triamine. The blend of ingredients was blown in a metal mould using polymethylhydrogensiloxane as a chemical blowing agent to obtain foams of different formulations.
Experimental
Materials
Epoxidized castor oil (ECO) (4% oxirane oxygen, MW∼980 g mol−1) was purchased from Jayant Agro-Organics Ltd. Mumbai, India. Ethylene diamine (EDA), diethylene triamine (DETA), tetrabutylammonium bromide (TBAB) (catalyst for carbonation of oil) and sodium sulphate (used for washing the carbonated oil) were purchased from Amrut Chemicals Mumbai, India. Hexamethylene diamine (HMDA) was purchased from Thermo Fischer Ltd. India. Commercial epoxy resin (EEW∼190 g eq−1) and its hardener was obtained from Bhor Bond Ltd. India. 1,8-diazabicyclo [5.4.0] undec-7-ene (DBU) (catalyst for synthesis of PHU of castor oil) was purchased from Avra Chemicals Ltd. Bangalore, India. Polymethylhydrogensiloxane (PMHSO) (MW∼1900 g mol−1) blowing agent was purchased from Alfa Aesar Ltd. India. Deionized water was used for product washing purposes. All chemicals were used as such without any further purification.
Synthesis of cyclic-carbonate of castor oil
Epoxidized castor oil and TBAB catalyst (0.035 mol per mole of epoxy group
35
) was charged into a glass reactor and heated to dissolve the catalyst. The carbonation reaction of oil was carried out by passing CO2 gas for 70 h at 120°C and 1 atm. For every 5 h, oil samples were withdrawn and checked for their reduction in oxirane value by HBr in acetic acid titration method as per ASTM D 1652. The reaction was stopped after a constant value was obtained as shown in Figure 1. The product obtained was washed twice with DI water and sodium sulphate solution to remove the catalyst. The reaction pathway for this step is shown in Scheme 1. The yield of carbonation reaction of epoxidized castor oil was estimated from equation (1). Reduction in oxirane value of oil with reaction time. 5-membered cyclic-carbonate from epoxidized castor oil.

The lowest oxirane value was 1.3 at the end of carbonation step (Figure 1), hence from equation (1), % yield = (4–1.3/4) × 100 = 67.5%.
Carbonate equivalent weight (CEW) of CCO was determined by integrating the carbonate peak in 1H NMR spectra by using trimethylsilane standard.
AHEW of EDA, HMDA and DETA were calculated from equation (2).
CEW and AHEW were used in setting the molar ratio for synthesis of PHU of castor oil.
Synthesis of amine terminated PHU of castor oil
Reduction in amine value of PHU with reaction time.

Amine terminated PHU from cyclic carbonate of castor oil.
Synthesis of foams
Formulation of foams (quantities in parts by weight).
aWeight percentage of PMHSO with respect to amine.

Reaction of PMHSO with amine to yield hydrogen gas.
Swell index and gel content
Swell index of foam samples was calculated as follows: Foams were immersed in 30 mL THF for 24 h. Weight of the samples measured before (w1) and after (w2) immersion was noted and swell index (%) was obtained from equation 3,
Gel content of foams was determined as follows: Foam samples were dried in oven at 60°C for 24 h. Weight of samples taken before (m1) and after (m2) drying was used to calculate gel content (%) from equation 4,
Density of foams
Density of foams was obtained by measuring dimensions of foams with the help of vernier calliper. Volume was obtained from dimensions and with known mass density was calculated from equation 5,
Nuclear magnetic resonance spectroscopy
Chemical structures of cyclic carbonate of castor oil and PHU of castor oil was studied under 1H NMR spectroscopy using a SA-Varian 400 MHz NMR spectrometer equipped with a OneNMR pulse-field gradient probe at room temperature. CDCL3 was used as solvent and trimethylsilane (TMS) as external reference.
Fourier transform infrared spectroscopy
Chemical structures of products were determined by Infrared Spectroscopy on Bruker instrument in the ATR (attenuated total reflectance) mode from 700 cm−1 to 4,000 cm−1 wavelength range with 100% transmittance.
Cellular morphology of foams
Internal cellular structure, cell size and cell compactness of the foams was studied by optical microscopy and scanning electron microscopy at various magnifications. Olympus BX41 polarized optical microscope (POM) and Hitachi N-3400 scanning electron microscope (SEM) were used.
Thermogravimetric analysis
TGA study was done to determine thermal stability of foams by using a Perkin Elmer Pyris 1 thermogravimetric analyser. Heating was done in the range of 30–600°C with 20 mL/min nitrogen purge at a heating rate of 10°C min−1.
Differential scanning calorimetry
Thermal properties of foams were studied by employing TA Q100 DSC analyzer to determine Tg where heating and cooling cycles were recorded between −40 and 100°C under nitrogen atmosphere at heating and cooling rate of 5°C min−1. Temperature versus heat flow was plotted and Tg of the foams were determined accordingly.
Thermal conductivity
Thermal conductivity of rigid foams synthesized from DETA was determined by Lee’s Disc experiment. The experimental setup consisted of circular brass metal disc, brass steam chamber, steam source and thermometers. Foam sample having 115 mm diameter and 2–3 mm thickness was kept between steam chamber at the top and metal disc at the bottom where steam chamber was equipped with inlet and outlet for steam. Two thermometers were inserted through holes, one in steam chamber and one in metal disc. Steam was then passed through the chamber and sample until the temperature indicated by thermometers became steady. Then steam chamber was removed and the bottom metal disc with sample kept over it was heated to 10°C above the noted steady temperature. Reduction in temperature with time was observed till the temperature of the bottom disc fell 10°C below the steady state temperature. Lastly time-temperature graph was plotted (Figure 10) and thermal conductivity of foam was determined by equation 6,
Where λ is thermal conductivity of foam, m is mass, s is specific heat of metal disc,
Thermal diffusivity is related to thermal conductivity and it was calculated using equation (7).
Compression and recovery study
The mechanical compressions of the foam samples were done on a UTM with compression holders using 20 mm × 12 mm × 10 mm foam slabs. A maximum force of 40 N with application time of 200 s and a maximum of 80% of strain was applied to the sample. Time of recovery of foam sample after removal of load was measured in seconds.
Resilience by ball rebound test
Resilience of foam was measured by the ball-rebound test according to a standard method. In this test, a steel ball weighing 16 g and 10 mm diameter was dropped from 30 cm height on the foam slab of 70 mm × 30 mm × 15 mm and the rebound height was then measured in percentage compared to the height of drop as shown in equation (8). The resilience value of regular Polyurethane (PU) foams is between 25 and 40%.
Determination of R-value
R-value also known as thermal resistance relates to insulation characteristic of foam. For closed cell polyurethane foams, it is around 0.97–1.4 m2KW−1. R-value of the foam was determined by heat transfer study on a heat transfer apparatus. The foam slabs of 120 mm × 90 mm × 18 mm dimensions incorporated between aluminum panels (Figure 13) were inserted in the heat transfer apparatus (Figure 2) between two sections which measured the temperatures TA and TB respectively. TA was gradually raised from room temperature and finally maintained at constant temperature of 40°C while changes in TB with respect to time were noted. R-value was determined at steady state conditions by employing equations (9–11). Heat transfer apparatus for determining R-value.

Results and discussion
Spectroscopy by NMR and FTIR
Figure 3a shows 1H NMR spectra of cyclic-carbonate of castor oil in which a peak at 4.5 ppm shows carbonate functionality obtained from carbonation of castor oil. Similarly Figure 3b shows NMR spectra of PHU of castor oil. In this spectra, disappearance of carbonate peak at 4.5 ppm and appearance of new peaks at 2.8 and 3.2 ppm corresponded to formation of urethane linkages. Also, peak at 4.6 ppm showed formation of hydroxyl group in the PHU of castor oil. Carbonate equivalent weight of CCO was determined from integration of peak at 4.5 ppm in 1H NMR titration using TMS standard, and it was found to be 340. 1H NMR spectra of (a) cyclic carbonate of castor oil and (b) PHU of castor oil.
Figure 4 shows FTIR spectra of ECO, CCO, Castor oil PHU and PHU foam. It was observed that epoxy peak at 1240 cm−1 of ECO disappeared and carbonyl peak at 1790 cm−1 of CCO appeared after carbonation reaction of the oil. Also, during synthesis of castor oil based PHU, carbonyl peak at 1790 cm−1 disappeared while CN stretch at 1070 cm−1, CO band at 1690 cm−1 and NH stretch at 3328 cm−1 could be seen which concluded the presence of urethane functionality. The appearance of SiH elongation at 2164 cm−1 showed the use of PMHSO blowing agent in the synthesis of PHU foam. CH2 symmetric and asymmetric stretching vibrations of linear chain of carbon atoms at 2925 cm−1 and 2857 cm−1 along with C=O at 1728 cm−1 of ester linkage could be seen in all spectra thereby confirming the use of castor oil in the synthesis of foams. FTIR Spectra of ECO, CCO, PHU of castor oil and PHU foam.
Density, swell index and gel content of foams
Figure 5 shows all the foams made in this study. Table 3 shows density, swell index and gel content of the foams. The highest density reported was 235 kg m−3 of PHU-DETA-F4 and lowest density was 155 kg m−3 of PHU-EDA-F1. The density values of each type of foam went on increasing with the decrease in quantity of castor oil based PHU in the formulation.
74
The foams made from EDA and HMDA were flexible due to long chain aliphatic ethylene and hexamethylene structures. Perhaps their density values were in the range 155–215 kg m−3. However, density values of foams made from DETA were comparatively higher 188–235 kg m−3. This is due to the cross-linking nature of DETA that imparts toughness and rigidity to the foam. Rigid foams also had low swell index values between 45 and 92% while the flexible foams had higher swell index between 162 and 325%. Furthermore, gel content of the foams made from DETA was highest at 99% compared to flexible foam at 91%. Gel content typifies the extent of cross-linking between the molecular chains
75
and it can be evidently said that for foams made from DETA it must be highest due to is cross linking nature. Other factors apart from structure of amine such as curing time and foaming time does play a significant role in fluctuation of gel content values. Polyhyrdoxyurethane foams. Physical properties of foams
Cellular morphology of foams
Figure 6 shows the cell structure of foams observed under optical microscope and Figure 7 shows cell structure as seen under scanning electron microscope. Foams made from linear aliphatic diamines, EDA and HMDA, were flexible and had cells (Figure 6(a) and (b)) of 142 μm diameter compared to cells of foam from DETA (Figure 6c) which were 115 μm in diameter. Optical microscopic images of foams. (a) PHU-EDA-F2, (b) PHU-HMDA-F1, (c) PHU-DETA-F2. SEM images of PHU foam at varying magnifications. (a) 20X, (b) 30X, (c) 100X, (d) 200X.

At the instance of foaming, flexible foams exhibited significant expansion which encouraged their cellular growth and therefore most of the cells were abruptly larger than their neighbouring cells (Figure 6(a) and (b)). Whereas in the case of rigid foam the cross-linking reaction due to DETA was predominant and since the foams became rigid while curing, its cell growth was restricted and cells were compacted (Figure 6c). More importantly, the foams made in this study had close cell morphology (Figure 7(a) and (b)) with cell size ranging from 115 to 142 μm in diameter. The cellular morphology of polyurethane foams is dependent on a number of factors that assist in forming the proper cell structure such as type of raw material, temperature of foaming, amount of blowing agent, curing temperature and time.
TGA analysis of foams
Degradation temperatures of foams.

TG (solid lines) and DTG (dot lines) curves of foams.
In this study epoxy resin was blended in castor oil based PHU. For epoxy polymers, the maximum temperature of degradation is between at 410–500°C. 77 While degradation for a general PU may last to a maximum of 350°C or higher, for foams made in this study, the maximum degradation temperature was about 395°C due to presence of epoxy resin in the foam formulation which obviously had an influence in increasing the temperature of degradation. From Table 4, it can be seen that at maximum degradation temperature flexible foams made from EDA and HMDA degraded more (78%, 74% at 395°C, 392°C) than rigid foam from DETA (70% at 385°C) which means that rigid foams had higher thermal stability than flexible foams which is a universal trend in PU foams. This can open avenues for their use as potential insulation materials for applications where higher thermal resistance is a prerequisite.
DSC of foams
Figure 9 shows the Tg values of foams made from various amines. The highest Tg reported was 35.3°C for rigid foam made from DETA while it was lower for HMDA and EDA at 31.5°C and 28.6°C respectively. It is pertinent that Tg increases as molecular weight of any of the component in the formulation increases. Also, foam having cross-linking component (DETA) in the formulation will have higher Tg as opposed to linear aliphatic component (EDA, HMDA).
78
Heat flow vs temperature curves of foams.
Moreover, foams contained epoxy resin in the formulation which caused an increase in Tg of the foams. 58 Therefore, unlike conventional flexible PU foams where Tg can go as low as −20°C, Tg values of PHU foams made in this study were between 28–36°C.
Thermal conductivity of foams
Rigid PU foams have insulative property due to their low thermal conductivity values between 48 and 50 mWm−1K−1.
79
Thermal conductivity of PU foam is dependent on cellular morphology, density and sometimes gas entrapped in pores of the foam which is released from blowing agent. Total heat transmission through closed cell foam is the sum of heat conducted through gas contained in the cells, heat conducted through the solid polymer of the cells and heat radiated across the cell walls.
56
Table 5 shows thermal conductivity and thermal diffusivity values of rigid foams made from PHU-DETA which were calculated from equations (9–11). Temperature vs time graph for thermal conductivity measurement. Thermal conductivity, thermal diffusivity and R-value of foams.
All foams had thermal conductivity in the range 104–130 mWm−1K−1 and thermal diffusivity in the range 0.2–0.3 mm2s−1. Foams made from DETA were rigid and densely compact in cellular structure. Cornille et al. 62 synthesized NIPU foams and stated that density and thermal conductivity values of the foams were closely related to each other. As density of the foams increased, the pore became smaller in size and larger in numbers which led to decrease in radiation of heat through the cells of the foam. Also, he stated that low thermal diffusivity values of NIPU foams (0.2–0.3 mm2s−1) justified that temperature flow within these materials is slow as opposed to PU foam (0.93 mm2s−1).
Compression and recovery analysis
Recovery of flexible foams made from EDA and HMDA after mechanical compression under UTM was tested and results were obtained as shown in Figure 11. All the foams showed full recovery after complete removal of load which can be attributed to excellent flexibility of foams made from castor oil, linear aliphatic diamines and epoxy resin. Foams with high density such as PHU-EDA-F3, PHU-EDA-F4, PHU-HMDA-F3 and PHU-HMDA-F4 took longer time to recover than foams with low density. Moreover, decrease in quantity of EDA and HMDA based PHU and increase in epoxy resin in the formulation led to slight increase in toughness of foams thereby imparting prolonged recovery times after compression of foams. However, complete recovery in all foams was observed due to influence of epoxy resin. Compression and recovery of flexible foams. (a) PHU-EDA, (b) PHU-HMDA.
Fauzi et al. 80 explained that complete shape recovery in epoxy foams is seen after removal of load. This can be due to optimum interconnection between cells, porous cell wall and homogeneous cell size. All these factors played a significant role in the PHU foams made in this study.
Resilience by ball rebound test
Resilience of foam is the property by virtue of which it can regain its original shape after being compressed by an object. Resilience of foams is checked by ball rebound test. Figure 12 shows resilience values of PHU foams obtained by ball rebound test. Flexible PU foams have resilience values between 25 and 40%. Resilience values of all flexible PHU foams made in this study were low compared to reference PU foam. A similar trend was observed between compression-recovery and resilience values of PHU foams. The foam with quickest recovery, 14 s, viz. PHU-EDA-F1 showed highest resilience value 16, followed by PHU-HMDA-F1 with resilience value 15. Similarly, the foam that took longest time to recover, 38 s, was PHU-EDA-F4 with the lowest resilience value 3. All the foams showed resilience values between 3 and 16. Resilience and flexibility are interdependent properties by virtue of which PU foams exhibit high compressibility along with complete and irreversible recovery. Resilience of flexible foams. PHU-DETA rigid foams incorporated between aluminium panels.

R-value of rigid foams
Figure 14 shows temperature gradient with passage of time inside heat transfer apparatus. A maximum of 5°C difference between section A and section B of heat transfer apparatus was observed for all the foams at steady state conditions. This difference ( Temperature gradient in heat transfer apparatus with time.
It was observed that R-values of the foams increased with the decrease in thermal conductivity of the foams and decrease in thermal diffusivity of the foams. Cornille et al. 62 explained this phenomenon on the basis of density and cellular morphology. That is, with the increase in density cell size decreases. The cluster of small sized cells resists the diffusion of heat through the foam thereby decreasing the thermal conductivity of the foams. From equation (9), it is clear that decrease in thermal conductivity will increase R-value of the foam. Therefore PHU-DETA-F4 had highest R-value of 0.17 mm2KW−1 followed by PHU-DETA-F3 having R-value of 0.16 mm2KW−1, PHU-DETA-F2 with R-value of 0.15 mm2KW−1 and PHU-DETA-F1 with R-value of 0.14 mm2KW−1. R-value is an important factor in determining the potential insulative properties of the rigid PU foams for their applications such as in building for insulation from solar heat, refrigerators, and thermal storage devices.
Conclusion
PHU foams were synthesised in this study by employing cyclic carbonate of castor oil, different aliphatic diamines and epoxy resin. Foams showed closed cell morphology with homogeneous cells. Thermal stability of foams was quite good due to presence of epoxy resin in the blend but lower than PU due to low reactivity of CCO and amine. Foams made from EDA and HMDA were flexible with low resilience values however they showed complete recovery after compression due to excellent shape recovery characteristics of epoxy foam. Foams made from DETA were rigid and had slightly high thermal conductivity than PU foams and low thermal diffusivity and low R-value. However, a maximum of 5°C of temperature difference was observed in heat transfer experiment. Foams with high R-values can be potential materials for insulation applications. Isocyanate free polyurethane route has been researched for decades now and it seems this route will soon be adapted provided toxic isocyanates are completely eradicated from chemical industries.
Footnotes
Acknowledgements
The authors would like to thank Department of Polymer and Surface Engineering, Institute of Chemical Technology Mumbai, for providing instrumental facility for successful completion of this research.
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
