Luca
Pellegrino
a,
Gunjan
Tyagi
a,
Eric S. J.
Robles
b and
João T.
Cabral
*a
aDepartment of Chemical Engineering, Imperial College London, London, UK. E-mail: j.cabral@imperial.ac.uk
bProcter & Gamble, Newcastle Innovation Centre, Newcastle upon Tyne, UK
First published on 16th November 2022
We investigate the phase behavior of model ternary triacylglycerol blends, comprising triolein (C57H104O6, OOO), tripalmitin (C51H98O6, PPP) and tristearin (C57H110O6, SSS), building upon extensive characterisation of single and binary mixtures, in order to rigorously map the thermal transitions of model natural ‘fats’. A combination of calorimetry, X-ray diffraction, and FTIR spectroscopy is employed to determine crystallisation and melting temperatures and identify the corresponding phases in the complex ternary system. We recover the eutectic behaviour of SSS-PPP blends and the invariability of OOO neat transitions, and resolve the complex β′ + β ternary surface, reflecting the roles of unsaturation and polymorphism of its constituents. Our results provide a representation of the OOO:PPP:SSS:temperature phase behaviour into a triangular prism, consistent with binary pair-wise data, which can inform a range of food science, cosmetic, pharmaceutical and cleaning applications that depend strongly on the physical-chemistry of such multicomponent ‘triglycerides’.
Despite their importance, considerably less is known about ternary TAG blends, and no complete phase diagrams have been reported. A classical study of the OOO:PPP:SSS system from 191236,37 provides reasonable estimates of the melting point, but does not consider the polymorphic nature of TAGs, which was not known at the time (in fact, X-ray crystallography was reported by von Laue a few weeks prior to its publication in the same year). In this study, we describe the ternary phase behaviour of a model OOO–PPP–SSS ternary blends. We take into account the polymorphism of the single TAGs and their binary blends, combining experimental and previous literature data as the foundation to provide a detailed description of the ternary phase behaviour, focusing on the effect of the unsaturated component OOO on the PPP–SSS saturated system by means of DSC, XRD and FTIR. The investigated ternary phase space is reported in Fig. 1(a), including the pure TAGs (yellow circles), the binary blends (blue circles) and the experimental ternary blend explored in this study (red circles). To benchmark our results with the previous literature data we repeat the characterisation of pure TAGs and selected OOO–SSS, OOO–PPP and PPP–SSS binary blends, numbered in the phase diagram. In the phase diagram are also indicated the approximate compositions of a range of natural fats described in terms of their OOO–PPP–SSS content with increasing unsaturation fraction towards pure OOO. Evidently, polymorphism, intersolubilty and phase behavior of a blend will be affected by different isomeric permutations of TAGs26,37–40 and, for the current system, 27 possible permutations are conceivable. Therefore, in this simplified, yet comprehensive description, we take OOO, PPP, and SSS as coarse-grained representations of the homotriglycerides and all respective O, P and S isomeric permutations with a majority O, P or S fatty acid chains. Due to the sensitivity to thermal processing, which impacts the crystallisation and melting processes, and to ensure reproducibility and comparability of experimental data, all TAGs and blends were characterised following the same thermal program depicted in Fig. 1(b). To further understand the complexity of natural fats, the respective FFA content is reported in Fig. 1(c) and the molecular structures of the OOO, PPP and SSS triglycerides employed in this work are shown in Fig. 1(d). Our findings, incorporate the prior literature on single triglycerides and their binary blends, and provide a rigorous description of the ternary blend phase behaviour into a triangular prism, including the effect of polymorphism on melting and crystallisation behaviours.
![]() | ||
Fig. 1 (a) Ternary phase space for the investigated triglyceride blends (full circles), literature data (hollow circles) indicating the mass fractions of triolein (OOO), tripalmitin (PPP), and tristearin (SSS); pure components are shown in yellow, binary blends (1–6) in blue, and ternary blends in red. Typical compositions of naturally occurring fats are indicated alongside the ternary blends. Blend 8 is formulated as a synthetic analogue of natural ‘beef tallow’ (0.5![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
![]() | ||
Fig. 2 Characterisation of neat TAGs. (a) DSC thermograms reporting the second heating of OOO, PPP and SSS at 5 °C min−1, of the thermal profile in Fig. 1(b). (b) FTIR spectra of the ‘fingerprint’ region (1500–650 cm−1) following the temperature profile, indicated in Fig. 1(b). (c) XRD profiles, measured at room temperature, showing the three principal polymorphs α, β′ and β. The top XRD patterns in each panel were obtained following the thermal profile (*), cooling at 5 °C min−1, resulting in a β phase for PPP and SSS; OOO is in the liquid phase at this temperature. The second XRD profile (cooling at 25 °C min−1), shown for PPP and SSS, resulted in β′; while the third XRD profile (cooling at 50 °C min−1), shown only for SSS, yields the α phase. (d) Ternary composition diagram showing the three neat components and schematics of the α (H, hexagonal), β (T, triclinic) and β′ (O, orthorombic) polymorphs. |
Fig. 2(b) depicts the fundamental vibrational modes of neat TAGs in the ‘fingerprint’ region at 25 °C (full-range FTIR analysis in the ESI,† Fig. S1). The FTIR spectrum of OOO exhibits broad infrared absorption bands in comparison to PPP and SSS, which exhibit sharp, intense peaks in the region. Mainly the frequency of the CH2 bending or “scissoring” band is dependent upon the methylene chain packing and conformation. Disordered liquid-like chains exhibit a considerably broadened and relatively lower intensity scissoring band in the region of 1468–1466 cm−1, while its shift and sharpening suggest ordering of the methylene chains approaching that of a crystalline geometry. Further changes in the CH2 ‘rocking’ mode of the triglyceride chains, which depends upon the inter chain interactions are indicative of structural alterations. For OOO, the absorption peak appears at 720 cm−1, indicating a hexagonal unit cell, whereas a sharper absorption peak appears at 717 cm−1, characteristic of a triclinic parallel unit cell. The saturated PPP and SSS on the other hand, display a larger number of sharp peaks i.e. from 1400 to 1100 cm−1 (C–H in plane bending, wagging ω and rocking ρ mode) whose number is equal to half of the number of CH2 of the fatty acid chains.43 The three basic polymorphs α, β′ and β, have distinct wide angle X-ray diffraction patterns as displayed in Fig. 2(c) and the schematic of the three crystal cages in Fig. 2(d). In general, a TAG α phase is characterized by a single strong reflection at 4.15 Å [(2θ) = 21.4°] of the hexagonal packing.17 The β′ phases are characterized by an orthorhombic packing, which results in two strong reflections around 4.2 and 3.8 Å, [respectively (2θ) = 21.1° and (2θ) = 23.3°]. The β-phase (triclinic) is identified by a very strong reflection at 4.6 Å, ((2θ) = 19.4°) and two other strong reflections at 3.9 and 3.7 Å ((2θ) = 23.1° and (2θ) = 24.1°). The three reflections at 4.6, 3.9 and 3.7 Å correspond, respectively, to the (010), (110) and (100) planes of the subcell. The d-spacing is identified for PPP and SSS by a single d003 reflection at 15.3 Å ((2θ) = 5.78°).
![]() | ||
Fig. 3 Binary TAG blend characterisation. The binary blends compositions are described in Fig. 1(a) and in Table 1 and numbered from 1 to 6. (a) DSC thermograms of the second heating (5 °C min−1) for different binary blends. (b) Triangular diagram describing the composition in OOO, PPP and SSS for different binary blends. (c) Room temperature XRD patterns of blend 1 to 6, showing the characteristic reflection for the different crystal structures identified. (d) Room temperature FTIR of the binary blends with assignments for the characteristic absorption bands. |
Blend | OOO | PPP | SSS |
---|---|---|---|
1 | 0.5 | 0.5 | — |
2 | 0.7 | 0.3 | — |
3 | 0.5 | — | 0.5 |
4 | 0.8 | — | 0.2 |
5 | — | 0.3 | 0.7 |
6 | — | 0.8 | 0.2 |
7 | 0.3(3) | 0.3(3) | 0.3(3) |
8 | 0.5 | 0.3 | 0.2 |
Isopleth 3PPP:2SSS | x | (1 − x)3/5 | (1 − x)2/5 |
The room temperature XRD patterns after thermal processing (indicated by a star in Fig. 1(b)) are reported in Fig. 3(c). For the OOO blends with either PPP or SSS, the reflections associated with mixed β′ + β phases are identified, following the characterisation in Fig. 2(c), and distinct from the liquid OOO phase, appearing as a ‘halo’ centered at (2θ) ∼ 20°. The presence of OOO does not provide additional reflections at low angle, meaning that the crystalline structure of the saturated components remains largely unchanged, even as the characteristic melting temperatures change. In addition, OOO provides a medium for crystalline PPP or SSS, where (β′ + β) crystals are dispersed. For blend 5 and 6, (0.7:
0.3 PPP
:
SSS) and (0.8
:
0.2 PPP
:
SSS) respectively, the presence of two reflections at 2θ ≃ 5.5° and 6.1° (d003) indicates that PPP and SSS crystallise in a β′ + β phase rich in β′. From the RT FTIR spectrum in Fig. 3(d), the two-component system can be identified by the presence of a shoulder emerging in the C
O stretching area and, in OOO blends, the absorption bands in the 717–719 cm−1 (rocking CH2, β′ + β triclinic and orthorombic) region are broadened and shifted to higher wavenumber (720–723 cm−1).43 Binary blends melting temperatures are plotted in the phase diagrams reported in Fig. 4, respectively for the OOO–PPP (a), PPP–SSS (b) and SSS–OOO systems (c). The experimental blends investigated in this work, extracted from DSC thermograms reported in Fig. 3(a), are labelled as blue full circles and triangles and are combined with DSC and dilatometry measurements form previous literature data1,26,30,33,34,44 (hollow circles and triangles). In OOO–PPP blends, the increase in OOO w/w reduces the onset temperature of the endothermic transitions of PPP, decreasing from the initial value of 65 °C to about 40 °C. In turn, PPP also affects the onset of the β′ phase of OOO increasing it from −10 °C to −3 °C. In agreement with Hale et al.,34 a plateau region for both triglycerides can be observed between 0.2 and 0.7 w/w OOO. Higher OOO content further reduces the onset the β PPP fusion temperature, and concurrently slightly increases the β′ OOO temperature. The β form of OOO, on the other hand, is largely unaffected by the presence of PPP at all concentrations, likely due to the greater stability and thermal reversibility of this polymorph. The presence of distinct melting points for the two triglycerides is indicative of phase separation, with possible partial miscibility at extreme concentrations, OOO-rich or PPP-rich, where a small upturn in
has been reported.14 Overall OOO–PPP is thus classified as a monotectic system (with partial solid solution at extremely high and low OOO fractions). By contrast, the PPP–SSS binary system is characterised by a eutectic temperature close to 64 ± 3°, with the minimum in the liquidus curve occurring at a composition between 16 and 25%, depending on the different techniques and resolution employed in previous reports.26,33,45 Another feature of this system is the presence of non-eutectic, therefore continuous solid solution of the α and β′ phases, first reported by Lutton et al.33 based on XRD and dilatometry measurements. The tendency for the lower melting point component to be dissolved to a greater extent, and for the eutectic point to be displaced towards a composition richest in the lower-melting component is commonly occurring in triglyceride mixtures. For the SSS–OOO system a similar behaviour to that of OOO–PPP system is observed, although the OOO melting point of both β′ and β forms is not lowered by the presence of SSS, denoting vanishing low solubility of SSS in solid OOO. The OOO is therefore included in the crystal lattice of SSS in weaker proportions (less than 2% w/w in the β phase) and the melting behaviour is similar to that of pure SSS.44
![]() | ||
Fig. 4 Phase diagrams for the investigated binary blends (a) OOO–PPP, (b) PPP–SSS and (c) SSS–OOO. Melting temperatures 1–6 (full blue circles, triangles) Tm extracted from DSC measurements reported in Fig. 3(a). Experimental data in this work are combined with data from previous works (hollow blue circle, triangles) from ref. 30, 33, 34 and 44. Yellow symbols, triangles refer to the pure components. For OOO–PPP and SSS–OOO all the reported α and β phases are measurable upon cooling below 45 °C and reheating at 5 °C min−1. The β′ phase reported in the (b) diagram can be formed cooling/annealing at 55 °C. In panels (a and c) there is no evidence of appreciable β and β′ transition temperatures in OOO induced by PPP and SSS. |
![]() | ||
Fig. 5 Ternary blend characterisation. (a) Ternary phase diagram, reporting the composition of the investigated OOO:PPP:SSS blends (red full circles) alongside the compositional data for binary blends reported in Fig. 4 as full blue circle (experimental) and hollow blue circle (literature). The pure components are reported as full yellow circles. (b) Room temperature XRD patterns of blend 7 (0.3![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
DSC characterisation of the explored ternary blends is shown in Fig. 5(c)–(e). Specifically, we focus on the high-melting region of the ternary bends Fig. 5(c), given by the phase transitions of the PPP–SSS pair because OOO transitions were not affected, as reported in Fig. 5(e) and already discussed for the binary blends. Fig. 5(d) reports the extracted melting temperature Tm for each phase identified in Fig. 5(c). Upon heating, at low OOO weight fractions (0.001–0.005 w/w) the phase behaviour is similar to that observed for binary blend 6 (PPP rich) with a common α(PPP+SSS) phase quickly re-crystallising in a mixture of β and β′ forms expressed as (β′ + β), melting in two different phases, a lower melting at 60 °C ascribed to the PPP fraction ((β′ + β)PPP) and a higher melting at 65 °C ascribed to the SSS fraction ((β′ + β)SSS). Overall, the Tm profile of the (β′ + β) phase appears lower than the values recorded for the pure components, as already discussed for binary blends. In this case, however, the decrease is observed at much lower OOO fractions. As the OOO fraction increases (0.005 to 0.01), the metastable α(PPP+SSS) melting peak is disappearing, re-crystallising and contributing to the (β′ + β) phase. This is observed in the DSC as a growing exothermic peak between the α and (β′ + β), until eventually disappearing at 0.06 w/w OOO. The Tm of PPP and SSS phases are initially reduced, then approaching a plateau region up to 0.1 w/w OOO. Here, until the full conversion of the α phase, is it possible to still identify the two (β′ + β) melting peaks for PPP and SSS, although progressively approaching and eventually merging into a single broad endothermic peak for 0.3 w/w OOO. At this stage, a single phase transition ((β′ + β)(PPP+SSS)) is observed. This is interpreted to a solid solution of PPP and SSS crystals in an OOO liquid medium that destabilises and favours the rearrangement of the acyl chains into a more favourable packing. A further increase in the OOO weight fraction leads to the same melting behaviour progressively moving towards lower Tm, until it is eventually no longer detectable for 0.99 w/w OOO. Due to the complex behaviour of the ternary mixture, assignment of the different phase transitions is not trivial, though we can summarise the addition of OOO to the PPP + SSS blend as follows:
• At low OOO w/w (0.001–0.005), the system behaves as a phase separated blend of PPP and SSS crystals incorporating small fractions of OOO;
• At intermediate OOO w/w (0.007 and 0.2), the PPP–SSS blend is not able to contain the liquid OOO, and the system rearranges towards a solid solution PPP + SSS with liquid OOO;
• For high OOO w/w (0.3–0.99), a single melting transition corresponding to the (β′ + β)(PPP+SSS) solid solution is observed progressively at lower temperatures. The system behaves as a dispersion of PPP + SSS co-crystals in a continuous liquid OOO medium.
In the reference study of Kremann and Schoulz in 1912,36 the melting behaviour of the OOO–PPP–SSS ternary system was described for the first time, including the diagram reproduced in the ESI,† Fig. S3, showing a maximum melting point of 64.4 °C corresponding to a mixture containing 0.303:
0.07
:
0.627 OOO
:
PPP
:
SSS. These results, though extremely insightful, deviate from our results in Fig. 5(c), most likely due to the experimental method employed or the chemical purity of components; further, the polymorphic behaviour of TAGs was not considered, and instead a single, average melting temperature was recorded. In accordance with that observed for the PPP–SSS binary system, the authors report a binary eutectic in the ternary phase space at 0.42 w/w SSS (0.58 PPP) and melting at 54 °C. The ternary system is represented by contours defined by the melting points, with the eutectic forming a valley descending towards the pure OOO corner and vanishing at 0.9 w/w OOO. Although this representation could be a simpler guideline especially in understanding the liquidus phase boundary, the solidus lines and solid solutions formed cannot be inferred due to the lack of information and characterisation of the different polymorphs. Here, subsuming the knowledge gained from the binary phase diagram and the experimental results of ternary OOO:PPP:SSS blends, we report a comprehensive description of the phase behaviour of the ternary blends in the triangular prism, composition vs. temperature, shown in Fig. 6.
![]() | ||
Fig. 6 3D ternary phase diagram built combining the experimental and literature data for pure components, binary and ternary blends reported respectively in Fig. 2, 4 and 5. The pure components are identified by yellow spheres, the binary and ternary blends by blue and red spheres respectively. The upper phase space identified the higher melting β form of the ternary blends (red areas) (β′ + β)SSS, (β′ + β)PPP+SSS and for the binary blends (blue lines). The surface beneath shows the ternary phase transitions for (β′ + β)PPP and α(PPP+SSS) alongside binary PPP–SSS phase transitions β′ and α. The lower phase space reports the phase behaviour and thermal transition for the alpha and β forms of OOO. |
The boundaries of the phase space are defined by the β phases of the high-melting components compiled from the three binary blends reported in Fig. 4, OOO–PPP, PPP–SSS and SSS–OOO (blue spheres and lines). The ternary phase space (red) is built from the melting temperature data from Fig. 5(c), where the outer region is defined by the (β′ + β)SSS and (β′ + β)(PPP+SSS) transitions. The inner region is defined by the (β′ + β)PPP and α(PPP+SSS) phases. The two phase region for low OOO weight fractions defines a depression followed by a plateau region in the β′ + β space, eventually leading to the PPP + SSS co-phase at high OOO loading. This behaviour defines a ‘peritectic’ region, usually observed only in mixed saturated/unsaturated systems where at least one triglyceride has two unsaturated acids as SOS/SOO or POP/POO.1 In our representation, the PPP–SSS eutectic follows the (accurate) description made by Lutton33 and we project its evolution in the ternary space with two horizontal red lines. We are able to describe, albeit in a simplified way, the different polymorphs for PPP and SSS and locate a new phase boundary, where most likely the eutectic extinguishes, at 0.4:
0.4
:
0.2 w/w OOO
:
PPP
:
SSS and the co-phase (β′ + β)(PPP+SSS) arises with a monotonic decrease in melting temperature as a function of the increase of the OOO content. All melting temperatures are found to be lower compared to the binary mixtures and pure components, therefore a depression can be identified, although distinct to that described by Kremann and Schoulz. As discussed in the binary and ternary blend characterisation, the OOO phase does not directly interact with the PPP + SSS system due to the distance in Tm. Therefore, in the ternary phase diagram the phase behaviour of the OOO β′ and β phases is described by means of triangular cuts defined by the respective melting temperature.
Natural fats are complex mixtures of triglycerides, fatty acids, and other minor components. Moving across the OOO tieline described in the ternary phase diagram allows us to explore different mixtures which, in OOO:PPP:SSS content, is approximate to the physical behaviour of several natural fats. In order to spatially resolve the structure and chemical interactions of the ternary TAG ternary constituents, we next consider a practically relevant, synthetic ‘beef tallow’ analogue. Natural beef tallow at ambient temperature exists as a mixture of β′ and β phases. Upon cooling from the melt, a β′(sub-α) phase is formed which transforms first in α and then to a mixture of α and β. At 16 °C, the α + β mixture is observed to transform into the normally observed β′ + β mixture noted above. Finally at 43 °C, the β phase disappeared leaving the β′ phase to melt at 48 °C. Thus beef tallow is a system characterised by a β phase which is not the highest melting phase.1,37,46 According to our ternary space representation beef tallow can be simplified as a mixture of OOO:
PPP
:
SSS in a 0.5
:
0.3
:
0.2 mass fraction. Therefore we produced a synthetic analogue of natural tallow, here referred to as blend 8, whose XRD structural characterisation and DSC are shown in Fig. 5. Additionally, we explored the melting behaviour of TAGs and in particular blend 8 to the molecular level. For this we took advantage of a Variable-Temperature Fourier-Transform Infrared spectroscopy (VT-FTIR) and FTIR-imaging. In VT-FTIR, a temperature ramp is applied to the sample by means of a controlled temperature stage where the sample is loaded, in reflection mode, between the torque screw and the diamond. The results are reported in Fig. 7, corresponding to spectra collected in continuous mode, heating from 25 °C to 75 °C at 5 °C min−1, probing the full range of PPP and SSS melting temperatures.
From Fig. 7(a), in the three spectral regions the temperature increase is visualised by a transition from sharp and well defied absorption peaks to broader bands, typical of a solid to liquid phase transition and associated reduction in the conformational order. Specifically in the ‘fingerprint’ region, the C–H in-plane bending mode structure (1350 to 1180 cm−1) typical of the solid β phase is lost and the spectra display a shift to lower wavenumbers and absorption bands similar to the liquid state of neat OOO at RT. The crystalline band centered at 718 cm−1 (inset), indicative of the crystal cage shifts from 717 cm−1, associated with the β′ + β phase, to a broad absorption band centered at 722 cm−1, indicating a mixed phase constituted by an increased fraction of hexagonal phases (720 cm−1) in the liquid fraction. For the CO and C–H regions, a shift to higher wavenumbers is observed. To track the evolution of the melting process, representative absorption shifts for each region were plotted as function of temperature, reported in Fig. 7(b). Specifically, for the “fingerprint region” we took the symmetric stretching of the ester group, νsCOOR, formed by one of the glycerol OH functional group and a COOH group from a fatty acid chain; for the “C
O region”, we considered the symmetric stretching νsC
O of the fatty acid carbonyl group; for the “C–H region” we took the symmetric stretching νsC–H of the fatty acids methylene groups. In all regions, a discontinuity is observed from 45 °C, indicating the onset of the melting process. From 50 to 55 °C and beyond, the shifts align to a new constant value. The red lines are guide to the eye, depicting a first-order solid to liquid phase transition. The simple nature of this transition implies the existence of no intermediate phases, since triglycerides cannot be brought into an amorphous solid state. The amorphous portion could be assigned to a bulk liquid state or it could be described in terms of a disordered portion of the crystallographic unit cell. To further explore the nature of the bulk liquid phase, we performed FTIR imaging of the CH2-CH3, reported in Fig. 7(c). This stretching band is associated with the CH2 group adjacent to the terminal CH3 in the fatty acid chains. Because TAGs pairing in the different crystal cages is determined by the alignment and orientation of the CH2–CH3 end groups, this absorption is ascribed to an ordered crystal phase and therefore not present in the liquid state. In Fig. 7(c), 200 μm × 200 μm absorbance maps of 2960 cm−1 at 25, 40 and 75 °C, with a resolution of 2 μm, report the evolution of the absorption band with temperature. The red area is associated with the high intensity peak that progressively disappears at high temperature, corresponding to the transition from an ordinate solid to an homogeneous bulk liquid (green) with crystal inclusions. These conclusions are quantified in Fig. 7(d), where a color deconvolution was performed on the FTIR maps extracting the red channel, associated with the solid portion and the green channel, associated with the liquid portion, with the areal coverage of each phase plotted as a function of temperature. The behaviour is again reminiscent of first order solid to liquid transition, although with a discontinuity not as sharp as that observed from the absorption shifts. This behaviour therefore represents not an order–disorder transition, but a gradual phase conversion with coexistence of both solid and liquid portions, even at temperature where the blend appears to be in a ‘liquid’ state. Further evidence of this melting behaviour is provided by polarised optical microscopy reported in Fig. S4 (ESI†), where changes in crystal morphology can be observed above the melting point of the blend, expected for the complex polymorphic behaviour of its components.
We expect our work to establish a multicomponent framework to describe the complex phase behaviour of TAG mixtures, in terms of the OOO–PPP–SSS ratio, enabling further studies, including the addition of fatty acids, or the role of specific isomers, which are intrinsic to naturally occurring TAGs. Moreover, the methodology can be exploited in the analysis of the melting behaviour of specific OOO:PPP:SSS combinations, as synthetic analogues of natural fats, such the one reported for “beef tallow”, providing a route towards defining simpler descriptive models that are relevant to a range of food science, cosmetics, pharmaceutical and cleaning/removal applications, which depend strongly on the physical-chemical properties and thermal history of these blends.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d2cp02395d |
This journal is © the Owner Societies 2022 |