
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

38589622
58838
10.1038/s41598-024-58838-6
Article
Exploring the effects of structure and melting on sweetness in additively manufactured chocolate
Burkard Johannes johannes.burkard@hest.ethz.ch

12
Kohler Lucas 1
Caciagli Sophia 1
Herren Nicolas 1
Kozamernik Mark 2
Mantovani Saskia 2
Windhab Erich J. 1
Denkel Christoph christoph.denkel@bfh.ch

2
1 https://ror.org/05a28rw58 grid.5801.c 0000 0001 2156 2780 ETH Zurich, Institute of Food, Nutrition and Health, 8092 Zürich, Switzerland
2 https://ror.org/02bnkt322 grid.424060.4 0000 0001 0688 6779 School of Agricultural, Forest and Food Sciences HAFL, Food Science and Management, Bern University of Applied Sciences, 3052 Zollikofen, Switzerland
9 4 2024
9 4 2024
2024
14 826119 12 2023
3 4 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ 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/.
In view of the health concerns associated with high sugar intake, this study investigates methods to enhance sweetness perception in chocolate without increasing its sugar content. Using additive manufacturing, chocolate structures were created from masses with varying sugar and fat compositions, where hazelnut oil served as a partial cocoa butter replacement. The study found that while variations in sugar content minimally affected the physical properties of the chocolate masses, hazelnut oil significantly modified melting behavior and consumption time. Chocolate masses with higher hazelnut oil content but similar sugar content exhibited a 24% increase in sweetness perception, likely due to accelerated tastant (i.e., sucrose) release into saliva. Multiphase structures, designated as layered, cube-in-cube, and sandwich structures, exhibited less sensory differences compared to the homogeneous control. Nonetheless, structures with hazelnut oil-rich outer layers resulted in an 11% increase in sweetness perception, even without sugar gradients. This suggests that tastant release plays a more critical role than structural complexity in modifying sweetness perception. This research highlights the efficacy of simpler multiphase structures, such as sandwich designs, which offer sensory enhancements comparable to those of more complex designs but with reduced manufacturing effort, thus providing viable options for industrial-scale production.

Subject terms

Health occupations
Engineering
Materials science
http://dx.doi.org/10.13039/501100001711 Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung 20B2-1-180971/1 20B2-1-180971/1 20B2-1-180971/1 20B2-1-180971/1 20B2-1-180971/1 20B2-1-180971/1 20B2-1-180971/1 20B2-1-180971/1 Burkard Johannes Kohler Lucas Caciagli Sophia Herren Nicolas Kozamernik Mark Mantovani Saskia Windhab Erich J. Denkel Christoph issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The marked increase in life expectancy, coupled with lifestyle changes in Western countries, have led to significant increases in obesity and non-communicable diseases (NCDs)1. Although added sugar consumption has declined in recent years2–5, current intake still exceeds the recommended 5 to 10% of daily energy consumption6. This discrepancy underscores the need for low-sugar food alternatives and compels food manufacturers to explore effective sugar reduction strategies7,8.

However, sugar replacement in processed foods is challenging because it contributes to both taste and texture, as demonstrated in products such as ice cream9, cookies10 or chocolate11. The predominant strategy for reducing dietary sugar involves replacing sugar with non-nutritive sweeteners (such as sucralose, aspartame, acesulfame potassium, saccharin, stevia, thaumatin or monk fruit)8, although limited by their bitter or metallic aftertaste12–15. An emerging alternative which dispenses with non-nutritive sweeteners focuses on the modification of food structure8,16–18. It seeks to alter taste perception by applying a concept known as “pulsatile stimulation”. This technique uses tastant (i.e., sucrose) gradients within samples to discontinuously stimulate taste receptors and enhance overall sweetness perception compared to a homogeneous control. This enhancement has been validated in a variety of matrices, including sweet solutions19–21, multi-layered gels16–18,22, jet-printed chocolate surfaces23, or 3D-printed chocolate24. For example, Kistler et al.22 found that hydrocolloid structures with a sucrose gradient, created using fused deposition modeling (FDM), were perceived as up to 30% sweeter than samples with homogeneously distributed sugar.

In this context, we used FDM to manufacture chocolate structures with inhomogeneous sucrose distribution, aiming to increase the sweetness perception. Critical physical properties, such as flow behavior and particle size distribution, were maintained uniform in both the low-sugar and the high-sugar chocolate masses, which served as masses for the multiphase structures. This methodological choice allowed us to focus how the spatial sugar arrangement influences taste perception in complex fat-based structures, while reducing the influence of confounding factors. In a first set of experiments, we prepared five multiphase structures from the two chocolate masses, varying from layered to cube-in-cube designs, which were then sensorially evaluated by a trained panel (N = 8) against a homogeneous control of similar total sugar content. A follow-up set of experiments comprised the design of structures from chocolate masses that differed in both sugar content and melting speed, achieved by partially substituting cocoa butter with hazelnut oil. Extensive analytical and sensory analyses were carried out on all the chocolate masses. The final stage involved the sensory evaluation of six sandwich structures, designs consisting of either one or three inner layers with distinct sugar and fat compositions and two identical top and bottom layers. The sandwich structures were evaluated by the same panel, using comparative and time-intensity methods (N = 8).

Our experiments explored the influence of mass composition, structural design, and manufacturing effort on the perceived sweetness of chocolate structures, all while maintaining controlled physical properties. This study aimed to improve our understanding of how sweetness perception in chocolate can be adjusted by additive manufacturing and relate it to pulsatile stimulation hypotheses.

Results and discussion

Chocolate mass analysis

Instrumental analysis

To evaluate the effect of inhomogeneous sucrose distribution on sweetness perception within multiphase structures, it was critical to ensure that all chocolate masses had comparable techno-functional properties. This uniformity should minimize the influence of somatosensory factors on taste perception25. To achieve this uniformity, especially in samples with different sugar content, sugar was partially replaced with a mixture of inulin and polydextrose, both known for their bulking properties26,27.Figure 1 Instrumental characterization of chocolate masses with different sugar (Sc) and hazelnut oil (Ho) contents. (a) Particle size distribution, (b) rheological flow curves, (c) differential scanning calorimetry (DSC) curves: the dashed lines delimit the distinct melting domains associated with different cocoa butter polymorphs28. In all three panels (a), (b) and (c) the standard deviation of each sample is depicted as the shaded area around the mean.

After seven hours of grinding, all chocolate masses exhibited a unimodal size distribution pattern with a d90 value below 30μm, which is essential to avoid a gritty texture29,30 (Fig. 1a). The particle size distribution was not significantly influenced by the sugar (F(2,7)[d90] = 0.34, NS) or hazelnut oil (Ho) (F(2,7)[d90] = 0.24, NS) content.

The flow curves in Fig. 1b demonstrate the non-Newtonian pseudoplastic behavior typical of chocolate31–33. Analysis of the viscosity at shear rates of 10, 20, and 50 s-1 showed a linear relationship with the sugar content. Notably, significant differences were found only at a shear rate of 10 s-1, as detailed in Supplementary Information 7. This is in contrast to previous studies by Aidoo et al.27,34, which reported an increased bulk viscosity for chocolates rich in polydextrose/inulin, due to their lower density and the resulting larger surface area that had to be coated with cocoa butter. In our study, the contrasting results may be explained by the irregular fracture patterns of crystalline sugar and amorphous polydextrose/inulin particles resulting from our small-scale grinding process. While this process predominantly affected particle shape, hazelnut oil (e.g., F(2,7)[η10] = 2.07, NS) and sugar content only had minimal effects on particle size distribution and flow behavior.

A consistent endothermic transition between 25.2∘C and 34.9∘C was observed in all the samples detailed in Fig. 1c, with melting peaks around 32∘C to 33∘C. The triglyceride composition of hazelnut oil closely mirrors that of cocoa butter, particularly in the chain length35. This similarity, combined with the moderate hazelnut oil levels used, indicates a uniform crystallization behavior in all chocolate masses, regardless of hazelnut oil content36. While cocoa butter substitution significantly affected several calorimetric characteristics, the most pronounced effect was found for melting enthalpy, which decreased from about 44-45 J/g in pure chocolate to 35-37 J/g in blends with 18 %w hazelnut oil/total fat. However, the sugar content (e.g., F(2,7)[melting enthalpy] = 1.93, NS) showed no significant effect on the melting properties, as detailed in Supplementary Information 7.

Sensory analysis

Figure 2 Sensory evaluation of chocolate masses with different contents of hazelnut oil (Ho) and sucrose (Sc). Panels (a), (b) and (c) show the average relative difference from control (ARD) for sweetness, smoothness, and bitterness perception, respectively. Panel (d) illustrates the interplay between sugar (Sc) and hazelnut oil (Ho) content across the three sensory attributes. Here, m represents the slope of the relationship between Ho content and perceptual ARD for each attribute and sugar level, where a positive m indicates a positive relationship and vice versa. ARD values with matching subscripts are not significantly different (p > 0.05). Panel (e) provides a correlation matrix comparing sensory ratings with instrumental characteristics of the chocolate masses, highlighting both significant (bold and framed) and relevant correlations for clarity.

The expected correlation between the sugar content and sensory sweetness was confirmed (F(2,98)[Sc] = 79, p < 2e-16). However, a more intriguing finding was the significant effect of hazelnut oil content on sweetness perception (F(2,98)[Ho] = 4.2, p = 0.01). Notably, the average relative difference of sweetness from the homogeneous control (hereafter referred to as ARDs) for the samples 42Sc-0Ho and 23Sc-18Ho were comparable, despite the former having nearly twice the sugar content (see Fig. 2a). Blending cocoa butter with hazelnut oil lowered the melting enthalpy as indicated in the previous chapter, reduced the consumption time (see Fig. 2e, R2 = -0.7, p = 0.002), and potentially increased the availability of tastants in saliva.

To isolate the physicochemical effect of hazelnut oil on sucrose release, two matrices with identical sugar content (30 %w/w) but different fat compositions were prepared. In one, the fat component was all cocoa butter (70 %w/w), while in the other, 18 %w/w of the cocoa butter was replaced by hazelnut oil, representing the ratio used in the 23Sc-18Ho sample. The panelists licked both samples for 45 seconds, after which the samples were expectorated and their soluble content was quantified by refractometry (see Supplementary Information 9). Our results showed that salivary secretion (g saliva/g sample) was not significantly different between the two samples, nor did it contribute to the differences in soluble solids content (F(1,33) = 3.58, NS) (see Supplementary Fig. S5c). While previous research has shown that taste influences salivary secretion37, it is likely that the compositional differences between the two samples were too subtle. Despite individual differences in Brix measurements among the panelists (see Supplementary Fig. S5a), a clear trend emerged: samples containing hazelnut oil had higher soluble solids contents (F(1,30) = 15.31, p < 0.001), on average more than 2 Brix higher than those made with pure cocoa butter. This enhancement is probably due to the increased availability of sucrose, a consequence of the differences in melting enthalpy between hazelnut oil and cocoa butter, which allows a more efficient release of sucrose from the fat matrix. These results suggest that the incorporation of hazelnut oil enhances sweetness perception.

Furthermore, the strong positive correlation observed between sensory sweetness and smoothness in our chocolate samples, as illustrated in Fig. 2b+e (R2 = 0.9, p < 10e-6), highlights the importance of understanding the physical factors contributing to sensory smoothness. Andablo-Reyes et al.38 suggested that the perception of smoothness is rooted in oral-tribological interactions, shifting the focus from viscosity to friction as a key determinant. Our tribology data, presented in Supplementary Information 10, shows that chocolate masses with higher sugar content exhibited increased coefficients of friction. This finding challenges the traditional belief of an inverse relationship between smoothness and friction39 and may be explained by the aggregation of sugar particles at low shear rates. The lack of a clear inverse correlation between both the friction coefficient and viscosity, and sensory smoothness, may indicate that the perceived textural differences in our chocolate samples may be influenced by perceptual factors rather than their physicochemical properties.

Sensory bitterness (ARDs) was most pronounced in the low sugar content and high cocoa solids samples (see Fig 2c). This increased bitterness can be attributed to alkaloids such as theobromine and caffeine, as well as flavan-3-ols found in cocoa solids, which are known for their bitter taste40. An inverse relationship was observed between both sugar (F(2,98) = 112.72, p < 2e-16) and hazelnut oil content (F(2,98) = 0.23, p < 0.02) and bitterness perception. This suggests that hazelnut oil can effectively mask the sensory bitterness41. This effect was particularly noticeable in chocolate masses with intermediate levels of sugar, as shown by a pronounced slope between hazelnut oil and sugar in Fig. 2d. At the extremes of sugar content, either the cocoa solids or sugar levels were too high for hazelnut oil to effectively suppress bitterness.

These findings indicate the possibility to modify multimodal perception in chocolate by up to 24% without changing the sugar content. As depicted in Fig. 2d, incorporating hazelnut oil not only increases sweetness and smoothness but also decreases sensory bitterness. This effect is most notable in samples of medium sugar content, as evidenced by the steep slopes (m) in Fig. 2d. Building on these results, it was of interest to explore to what extent these sensory profiles, especially sensory sweetness, are retained or altered when structuring multiple chocolate masses together by FDM.

Multiphase structure analysis and pulsatile stimulation

Figure 3 (a) Evaluation of the sensory attributes of multiphase chocolate structures during experimental set 2 (A = Sw23Sc-0HoI3, B = Sw23Sc-18HoI3, C = Sw42Sc-0HoI3, D = Sw42Sc-18HoI3, E = Sw42Sc-0HoI5, F = Sw42Sc-18HoI5), with results presented as average relative differences in perceived sweetness, smoothness, and bitterness from the control (ARD). ARD values with the same subscript are statistically indistinguishable (p > 0.05). Panel (b) shows the correlation between the overall and maximum sweetness ARD, including a linear fit. Panel (c) summarizes the temporal evolution of sweetness perception as assessed by time-intensity measurements. The arrows at the beginning (Istart,Ref), at the maximum (Imax,Ref), and at the end (Iend,Ref) of the plot indicate the dimension of the deviations from the homogeneous reference at different time points. The time regime is divided into three stages: onset, maximum plateau, and decay. It is shown for the reference sample (light gray line below the sample time-intensity curves). Furthermore, the horizontal arrow between the two vertical hatched areas indicates the delay of the peak intensity onset caused by the substitution of Ho in the tongue proximal layer. The standard deviation of the reference sample is shown (shaded area) to illustrate the variability in sensory perception.

The first set of experiments focused on modifying sweetness perception by inhomogeneous sucrose distribution. Exposure to sucrose was critical for eliciting a sweetness response, which was particularly pronounced in multiphase structures with high sucrose content in the outer layers, which facilitated the immediate release of sucrose upon melting. In particular, the structure CiC42Sc-0Ho I 2 was perceived as 31% sweeter than CiC04Sc-0Ho I 2, as indicated by their ARD at tstart, seen in Supplementary Fig. S4. This observation was further supported by the time-intensity data from the second set of experiments, detailed in Fig. 3c. While this disparity, particularly between Sw42Sc-0Ho I 3 and Sw42Sc-18Ho I 3, was pronounced at the onset of consumption (tstart), it gradually diminished as consumption progressed.

In contrast to chocolate masses, multiphase structures did not show a strong correlation between sweetness perception and smoothness or bitterness (Fig. 3a). This observation extends to both maximum and overall sensory perception, where a consistent correlation is shown across all attributes, with sweetness perception shown as an example in Fig. 3b. In comparison, Khemcheevakul et al.24, who also used FDM to create multilayer chocolate rings, managed to reduce sugar content by 19% in their chocolate multiphase structures without compromising sweetness or consumer preference. However, they did not extensively control the textural differences between high-sugar and low-sugar chocolate masses, which are critical in modulating taste perception at a physicochemical level. This aspect of the influence of texture on taste perception is further emphasized by Mosca et al.17, who showed that altering fracture properties in multi-layered hydrogels can manipulate taste perception by changing the total surface area, a critical factor in taste perception.

Further, it is well documented that sweetness perception can be influenced not only multimodally, but also by the pattern in which taste stimuli are presented16–22,24,42,43. Specifically, a discontinuous presentation of sweet tastants has been shown to modify the perceived intensity. However, our research suggests that changes in fat composition alone, even without changing sweetness levels, can enhance sweetness perception. For example, the ARD in sweetness perception for the structure L23Sc-18Ho I 3 was 15% higher than that of L23Sc-0Ho I 3. This finding suggests the involvement of multiple mechanisms, colloquially referred to as “pulsatile stimulation”.

One explanation is inspired by the “Brück-Bartley” effect observed in vision science, where a flickering light is perceived as brighter than a steady light44. This analogy is applied to taste, where rodents exposed to high-frequency pulsed taste stimuli have been reported to experience a greater taste magnitude45,46, but is unlikely to apply to semi-solid foods19. An alternative hypothesis is that alternating sucrose concentrations could prevent taste adaptation47–49. However, our experiments comparing the structures (L42Sc-0Ho I 2 and L42Sc-0Ho I 12) with different numbers of layers (the transfer from one layer to another being supposed to account for a pulse) did not show increased sweetness perception with more layers. This suggests that adaptation occurred regardless of pulse frequency, or because layers of different composition merged into a single perception at later stages. More compellingly, research suggests that expectation, shaped by textual or direct taste cues, significantly influences taste perception50–52. A study by Wilton et al.53 found that subjects rated low-sweet solutions as much sweeter when they were presented with conflicting cues, highlighting the power of expectation. Our study supports this, suggesting that the initial sweet stimulus sets a contextual framework that influences overall taste judgments. This aligns with the concept of cognitive bias, where the first impression of sweetness is preconsciously stored and integrated into overall perception54. In summary, our findings indicate that the variation in sweetness perception in semi-solid foods through pulsatile stimulation is largely due to this initial contextual shift.

Synthesis and industrial relevance

Manufacturing effort in additive manufacturing (AM), defined here as the relative difference in time required to create a multiphase structure by AM compared to a single-phase homogeneous control, is a significant barrier to the adoption of AM in the food industry55. Our study highlights this challenge, with production time already exceeding four minutes for a single-phase cube (referred to as 0% on the x-axis in Fig. 4a). This increase is primarily due to the complexities of handling chocolate: as a non-Newtonian printing medium it requires careful management of crystallization kinetics and printing speed to achieve precise 3D structures. These factors are critical to ensure sufficient yield stress for layer stacking and to prevent shape distortion56.Figure 4 (a) Summary of the average relative difference from the control (ARD) for sweetness perception in both experiments (y-axis). The x-axis represents the manufacturing effort, calculated as the relative difference in manufacturing time between each multiphase structure (tmp) and a single-phase homogeneous control (tc) when printed by FDM. For reference, the chocolate mass 23Sc-18Ho is added and positioned at 0 on the x-axis for comparison. Although this sample was manufactured by molding (subtractive manufacturing), its position represents a scenario where the sample is made by FDM. A box on the graph summarizes the factors influencing sweetness perception: positive features are marked in green, negative ones in red. The “relevant zone” (green area) indicates an empirically defined region of reduced manufacturing effort and increased taste enhancement, in contrast to the “ineffective zone” (red area). Arrows (i) and (ii) emphasize the importance of the initial phase in sweetness perception, showing similar trends across different structures: (i) for layered samples and (ii) for cube-in-cube structures. (b) Steps of the FDM printing process, divided into pre-printing (steps 1-2) and post-printing (steps 3-4), which are crucial for precise material deposition. (c) Two possible industrial setups for fabricating multi-phase structures: slit-die strand extrusion and single-shot ejection.

This is particularly challenging for multiphase structures, such as cube-in-cube designs, where print tool transitions require additional retraction and extrusion steps to prevent dripping and nozzle blockage57–59, as shown in Fig. 4a+b. In contrast, layer-by-layer deposition techniques, which are suitable for creating sandwich structures, could significantly reduce the manufacturing effort. These structures, especially those with outer layers rich in sugar and/or hazelnut oil, ensure a sweetness perception comparable to more complex designs, but are simpler to manufacture, making them more viable for industrial-scale production, as shown in Fig. 4a. Specifically, sandwich structures, such as Sw42Sc-18Ho I 3, not only have increased tastant availability due to sugar and hazelnut oil enriched contact areas, but also can be efficiently produced using parallel slit-die extruders (Fig. 4c). While cube-in-cube structures could be manufactured by the same technology, they would require more advanced slit-die configurations. However, these would not be able to match the speed of single-shot processes (Fig. 4c), which are the optimal scenario for products containing only spherical objects with a sweet outer layer.

Despite these speed limitations, AM processes offer significantly greater production flexibility and the ability to create texturally more complex structures. This flexibility is a key advantage of AM, allowing for a wider range of product designs and textures that can enhance the consumer experience.

Conclusions

This study investigated the interplay between sweetness perception and the spatial distribution of tastants (sucrose) and fat phase (hazelnut oil fraction) in chocolate structures. Our results reveal that, while variations in sugar content minimally affected physical properties, partially substituting cocoa butter with hazelnut oil markedly decreased both melting enthalpy and consumption time. This change likely accelerated the release of sucrose into the saliva, a theory corroborated by in vivo expectoration tests. These tests demonstrated that the increased tastant release could account for the observed increase in sweetness perception, which reached up to 24% in samples containing hazelnut oil. Interestingly, our study could not confirm the generally recognized inverse correlation between the rheological and tribological properties of chocolate and its perceived smoothness, suggesting that smoothness may be integrated with sweetness perception at a cognitive level.

For multiphase structures, sensory differences were less pronounced. Experiments indicated that sweetness perception was enhanced in cube-in-cube structures with high sugar content in the outer layers, underscoring the importance of initial tastant release in setting a perceptual baseline for overall sweetness perception. However, increasing the structural complexity, such as adding more layers with different sugar content, did not enhance sweetness perception, suggesting that structural intricacy may have a limited effect on altering sweetness perception. Sandwich designs with identical top and bottom layers achieved sweetness perception levels comparable to the cube-in-cube structures. Remarkably, sandwich structures featuring only hazelnut oil gradients were perceived as equally sweet as those with sugar gradients, challenging the prevailing assumption about the necessity of sugar gradients for sweetness enhancement.

Our research suggests that in solid, fat-continuous systems like chocolate, the effect of sugar gradients on sweetness perception is small, particularly when the analytical textural differences within the samples are small. These subtle effects are further diminished in complex designs or when initial contact surfaces have contrasting sugar gradients that neutralize each other. In addition, the manufacturing effort is a significant barrier, especially for cube-in-cube designs, due to the numerous printer-tool transitions that require additional effort to prevent dripping or nozzle clogging. Therefore, future research should prioritize the development of simpler, yet effective, multiphase structures. These should be based on formulations that optimize tastant release and exhibit pronounced analytical and sensory texture gradients. Such advances will change the sensory perception of food products like chocolate, both physicochemically and perceptually.

Methods

Chocolate mass manufacturing

All chocolate samples listed in Table 1 were standardized with respect to solids content (%w/w), and total fat content (%w/w). The sugar content in the samples ranged from 4 %w/w to 42 %w/w. This was achieved by replacing sugar (sourced from Zuckerfabriken Aarberg AG, Switzerland) with varying amounts of polydextrose (Litesse®Two Powder (IP Grade), Danisco AG, Netherlands) and inulin (OraftI®HSI, BENEO-Orafti S.A., Belgium) in a constant ratio of 2.17:1 (polydextrose to inulin). All non-cocoa ingredients, i.e.  polydextrose, inulin, and sucrose, were pre-mixed with melted cocoa mass (“Rondo”, Max Felchlin AG, Switzerland) in a blender (Titanium KMT056, Kenwood Swiss AG, Switzerland) before being transferred to a small melanger (ECGC-12SLTA, Cocoatown LLC, USA). Additional melted cocoa butter or hazelnut oil (Haselnussöl, Demeter, Switzerland), the latter added in the range of 0 to 18 %w/total fat, was added to the mixture after 90 minutes of refining. A detailed description of the preparation and the subsequent molding process (for single mass sensory and instrumental analysis) can be found in Burkard et al.23. For molding, the seeded chocolate masses were transferred into silicone molds (Silicone Culinair®, Siliconesandmore, Netherlands) of defined geometric design (16 mm x 9.6 mm x 12.8 mm). A schematic overview of the mass production can be found in Fig. 5a and the following analytical and sensory characterizations are shown in Fig. 5b.Table 1 Composition of chocolate masses 04Sc-0Ho, 04Sc-18Ho, 23Sc-0Ho, 23Sc-9Ho, 23Sc-18Ho, 42Sc-0Ho, 42Sc-18Ho. The chocolate masses were labeled using the syntax Sugar (%w/w)Sc-Hazelnut oil (%w/w)Ho content, where Sc represents the sucrose content of the chocolate masses, while Ho denotes the proportion of hazelnut oil relative to the total fat content of the samples.

Composition	Chocolate masses	
04Sc-0Ho	04Sc-18Ho	23Sc-0Ho	23Sc-9Ho	23Sc-18Ho	42Sc-0Ho	42Sc-18Ho	
Sugar content (%w/w)	04.15	04.15	23.00	23.00	23.00	41.71	41.71	
Hazelnut content/ total fat (%w/w)	00.00	17.96	00.00	08.98	17.96	00.00	17.96	
Total fat (%w/w)	35.09	35.09	35.09	35.09	35.09	35.09	35.09	
Cocoa mass (%w/w)	50.70	50.70	50.70	50.70	50.70	50.70	50.70	
Sugar (%w/w)	04.15	04.15	23.00	23.00	23.00	41.71	41.71	
Inulin (%w/w)	11.83	11.83	05.55	05.55	05.55	00.00	00.00	
Polydextrose (%w/w)	25.73	25.73	13.16	13.16	13.16	00.00	00.00	
Cocoa butter (melted) (%w/w)	06.30	00.00	06.30	03.15	00.00	06.30	00.00	
Cocoa butter (crystallized) (%w/w)	00.90	00.90	00.90	00.90	00.90	00.90	00.90	
Hazelnut oil (%w/w)	00.00	06.30	00.00	03.15	06.30	00.00	06.30	
Soy lecithin (%w/w)	00.39	00.39	00.39	00.39	00.39	00.39	00.39	

Single mass characterization

Particle size analysis

The particle size distribution of the chocolate masses was measured using a Beckman laser diffractometer (13 320 XR, Beckman Coulter Inc. , United States). For sample preparation, the melted chocolate masses were suspended in Hydriol®SOD.24 (Hydrior AG, Wettingen, Switzerland) and the instrument was operated with Hydriol®SOD.24 as solvent. Prior to measurement, the suspended sample was sonicated for two minutes at 65 Hz in an ultrasonic bath (Elmasonic P 30 Hat, EMAG Nänikon, Switzerland) at 50∘C. The diffraction pattern was interpreted using the Fraunhofer model60. All samples were measured in triplicate (N = 3).

Rheology

Rheological properties were measured using a rheometer (MCR 302, Anton Paar, Austria). Approximately 35 g of chocolate was melted at 50∘C, thoroughly mixed with 0.36 g of crystallized cocoa butter, and transferred to a couette flow cell (CC27, Anton Paar, Austria) at 40∘C. The sample was presheared for two minutes at a shear rate of 10 s-1 at 40∘C. For a logarithmic shear rate ramp ranging from 0.1 s-1 to 100 s-1, ten data points per decade were collected. All data points were collected under steady state conditions. All measurements were performed in triplicate (N = 3) and analyzed in Origin (OriginPro 2021, 9.8.0.200).

Differential scanning calorimetry

The melting profiles of the chocolate masses were characterized by differential scanning calorimetry (DSC 3+/500, Mettler Toledo GmbH, Switzerland). An indium and a water sample were used for calibration. 5 ± 0.2 mg of tempered chocolate was weighed into a 30μL aluminum crucible (Mettler Toledo GmbH, Switzerland). For measurement, the sample was equilibrated at 20∘C for two minutes, followed by a linear heating ramp to 50∘C at a heating rate of 4 K/min. All samples were measured in triplicate (N = 3). Curve characteristics such as onset, offset, peak melting temperatures and melting enthalpy were extracted using Origin (OriginPro 2021, 9.8.0.200) for analysis.Figure 5 (a)–(c). A detailed overview of the preparation of the chocolate masses and the additively manufactured multiphase structures can be found in section (d). Multiphase structures were subjected to two different sensory sessions, referred to as experimental sets 1 and  2. All structures were prepared using one (for control) or two sensorially characterized chocolate masses. In addition, the chocolate masses were analytically evaluated by calorimetry, rheology, tribology and laser diffraction, as described in section (b). In the first set of experiments, sugar was inhomogeneously distributed within cube-in-cube (CiC) and layered (L) structures, varying the manufacturing effort and complexity of the latter by changing the number of layers. The second set of experiments focused on sandwich designs, where the number of layers and the composition of each layer in terms of sugar (Sc) and hazelnut oil (Ho) content were varied. Importantly, the total content of sugar and hazelnut oil was maintained constant across these designs.

Multiphase structure manufacturing

The multiphase structures were labeled with respect to the experimental design, referring to the chocolate mass initially placed on the tongue, and the number of layers (designchocolatemassontongue|numberoflayers). The design notation consisted of abbreviations for layered (abbr.  L) or cube-in-cube (abbr.  CiC) designs, followed by the specific chocolate mass and the number of layers (e.g., Sw23Sc-0Ho|3 was a 3-layer sandwich structure with a bottom and top layer composed of chocolate with 23 %w/w sugar content without hazelnut oil).

In the first set of experiments, three layered structures, two cube-in-cube structures, and a homogeneous control with no sugar gradient and a sugar content of 23 %w/w, were prepared. All multiphase structures in the first set of experiments were made from chocolate with high (42Sc-0Ho) and low (04Sc-0Ho) sugar content. The number of layers of the layered structures varied, allowing to account for structural complexity: either two (L42Sc-0Ho|2), six (L42Sc-0Ho|6), or twelve (L42Sc-0Ho|12), whereby the sweet layer was always placed in contact with the participants’ tongues. For the cube-in-cube structures, high-sugar (CiC42Sc-0Ho|2) or low-sugar chocolate (CiC04Sc-0Ho|2) completely enveloped the other chocolate mass.

In the second set of experiments, sandwich structures were designed, in which the top and bottom layers were composed of the same masses, a distinct difference from the layered configurations. In addition, a fraction of the cocoa butter was replaced with hazelnut oil (Ho). As a result, six sandwich designs were prepared that differed in the number of layers (either three or five) and the composition of the top/bottom chocolate mass, resulting in the multiphase structures Sw42Sc-18Ho|5, Sw42Sc-18Ho|3, Sw42Sc-0Ho|3, Sw42Sc-0Ho|5, Sw23Sc-0Ho|3, and Sw23Sc-18Ho|3. Analogous to the first set of experiments, all multiphase structures were evaluated against a homogeneous control that maintained the same sugar content and, due to the introduction of hazelnut oil, the same hazelnut oil content (9 %w/total fat).

The dimensions of the multiphase structures were set at 16 mm x 9.6 mm x 12.8 mm. The total sugar (23 %w/w) and hazelnut oil (0 %w/total fat for the first, and 9 %w/total fat for the second experiment) contents were held constant with respect to the homogeneous control (notated as 23Sc-0Ho for experimental set 1 and 23Sc-9Ho for experimental set 2). A comprehensive illustration of the experimental setup and the multiphase structures is provided in Fig. 5d.

All chocolate masses were seeded with 1 %w/w of βV/VI cocoa butter prior to being transferred to stainless steel cartridges, equipped with a 15 mm long nozzle with an inner diameter of 0.8 mm (Sigrist & Partner AG, Matzingen, Switzerland). To maintain a constant product temperature, both the cartridge and nozzle temperature were set at 32.5 ± 0.3 ∘C using an aluminum jacket connected to a tempered water bath (Julabo 9153618 FP35, JULABO GmbH, Germany). FDM printing was performed using a custom-built three-axis Cartesian printer, designed by the Institute for Print Technology (Bern University of Applied Sciences)<f. The printer setup is shown in Fig. 5c. It was placed in an environmental chamber (KK-1000 CHLT, Kambic, Slovenia) to maintain the ambient temperature at 10∘C. Printer settings, including print speed, were defined in the Slic3r software (version 1.3.1). Printer operation was controlled by the Repetier Host software (version 2.1.6, 2019). After printing, all samples were stored at 5∘C prior to sensory testing. All printer settings are detailed in Supplementary Information 1.

Sensory analysis

All participants (ten women aged between 35 and 55 years) were recruited by the HAFL Sensory Evaluation Department and participated in 21 sensory sessions. In accordance with the constraints imposed by the COVID-19 pandemic, experimental set 1 was conducted at home, while all other tests were conducted in the HAFL sensory laboratory. The sensory analysis consisted of an in-depth study of the individual chocolate masses, and of a sensory characterization of the 3D-printed multiphase structures during two experimental sets. All sensory experiments were performed by the same group of participants (N = 8), while two participants were replaced during experimental set 2.

A comprehensive summary of all sensory experiments is presented in Table 2. All test sessions/methods were conducted in accordance with the relevant ISO standards and guidelines (ISO 6658 217 (E)), were approved by the University (BFH-HAFL) and guided by a sensory expert from HAFL. To ensure standardized consumption, participants were instructed to place the sample in their mouth and press it lightly against their palate. They were allowed to swallow if necessary and to move their tongue periodically during consumption. Consumption was considered complete once the entire sample had been swallowed. Still water (at 25∘C) and unsalted crackers (M-Classic Microc, Migros SA, Switzerland) were provided to clean the palate between samples. Prior to enrollment, all participants signed an informed consent document indicating their understanding and acceptance of the study procedures. To assess the panel’s effectiveness, we first evaluated individual performance in terms of discrimination, agreement, and repeatability. The panel’s collective performance was then evaluated using the same criteria and can be found in Supplementary Information 11.Table 2 Overview of all sensory tests conducted, both on single-mass and on multiphase structures: methodology used, time sequence of the analysis, experiment details, and attributes investigated. The sensory tests were always preceded by training sessions.

Analyzed sample	Methodology	Sequence	Training	Evaluations	Remarks	Attributes	
Single mass	CATA	1	0	1	-	Flavor (texture, taste and aroma)	
Comparative profiling	2	3	2	-	Sweetness, smoothness, bitterness (maximum)

and time of melting

	
Multiphase structures	Progressive profiling	3	4	2	Experiment 1	Sweetness (start, max and end)	
Comparative profiling	4	2	2	Experiment 2

–

	Sweetness, smoothness, bitterness (maximum

and overall) and time of melting

	
TI	5	3	2	Experiment 2

–

	Sweetness (temporal)	

Single mass characterization

Check-all-that-apply: First, participants performed a check-all-that-apply (CATA) test to describe the sensory profiles of four defined chocolate masses (04Sc-0Ho, 04Sc-18Ho, 42Sc-0Ho and 42Sc-18Ho). These masses were chosen with the expectation that they would cover the full range of sensory perception. Panelists were presented with a list of ten predefined attributes (see Supplementary Information 2), and asked to select the attributes that best described each sample.

Comparative profiling: From the CATA results detailed in Supplementary Information 6, the attributes bitterness, smoothness, and total consumption time were identified as significant in differentiating the chocolate masses. Over three sessions, participants were trained to ensure consistent ratings for the three attributes in addition to sweetness. They practiced rating the maximum intensity of each attribute on a sensation scale ranging from 0 (not bitter/sweet/smooth) to 100 (extremely bitter/sweet/smooth). Participants consumed the control sample 23Sc-9Ho and were instructed to rate its maximum sensation as 50 for all attributes. The rating of sweetness and smoothness was practiced with defined test samples (see Supplementary Information 3). The samples were tempered to 20 ± 1 ∘C and served on plastic trays with random three-digit codes. Each session was conducted by a sensory panel leader, a control warm-up sample was presented at the beginning of each session and the samples were presented in a balanced order. A blind control was also included in the experimental design. Data were collected using EyeQuestion (licensed to HAFL, Zollikofen, Switzerland), including melting time, defined as the time between the start of oral processing and the moment when the chocolate had been completely swallowed.

Multiphase structure analysis

Progressive profiling: Multiphase structures in the first experimental set ( L42Sc-0Ho|2, L42Sc-0Ho|6, L42Sc-0Ho|12, CiC42Sc-0Ho|2 and CiC04Sc-0Ho|2), and a homogeneous control (23Sc-0Ho) were evaluated at three time points tstart, tmax, and tend. Similar to the chocolate masses comparative profiling method, participants were asked to rate the sweetness sensation on a scale from 0 (not sweet) to 100 (extremely sweet). Each session began with a warm-up sample of a defined maximum sweetness intensity of 50. In addition to rating the maximum sweetness intensity Imax, participants were instructed to rate the initial (Istart) and final sweetness perception (Iend), the latter being defined as the sensation before the last piece of sample was completely swallowed. Due to the pandemic, the training and testing sessions were carried out online: the participants were at home and guided remotely by a sensory panel leader. Data from both sessions were recorded using EyeQuestion (licensed to HAFL, Zollikofen, Switzerland).

Comparative profiling: All multiphase structures in the second experimental set (Sw42Sc-18Ho|3, Sw42Sc-18Ho|5, Sw42Sc-0Ho|3, Sw42Sc-0Ho|5, Sw23Sc-0Ho|3, Sw23Sc-18Ho|3) were evaluated by comparative profiling on a unipolar scale from 0 to 100. In addition to rating the maximum perception of the attributes sweetness, smoothness, and bitterness, participants also rated the overall perception of these attributes. After three training sessions, the panel rated the samples twice in two test sessions.

Time-intensity: During three training sessions, participants were instructed to rate the temporal intensity of the sweetness. Specifically, the participants were instructed to begin rating the sweetness intensity after 10 seconds of consumption and continue at 0.1 Hz intervals for a total of three ratings. After this initial period, participants were allowed to continue rating the samples at their own pace. To provide further guidance, participants were asked to rate the intensity when the sample had completely melted (Imelt), as well as the final perceived sweetness sensation (Iend) upon swallowing. Further details can be found in Supplementary Information 4.

Descriptive statistics

A comprehensive breakdown of the statistical analyses is provided in Supplementary Information 5. The results from the Check-All-That-Apply (CATA) approach were analyzed using Cochran’s Q tests to identify significant differences among the attributes. For all other sensory evaluation data, we employed linear mixed-effects models in RStudio to account for the fixed and random effects within our experimental design. Significant findings (p < 0.05) were further explored with a Tukey Honestly Significant Difference (HSD) post-hoc test to pinpoint specific group differences. For an in-depth review, of the specific models applied to each sensory evaluation, please refer to the details in Supplementary Information 5.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-58838-6.

Acknowledgements

The authors thank the Swiss National Science Foundation and Innosuisse for funding the Bridge Project (No. 20B2-1-180971/1). We thank Daniel Kiechl for his help with installing the setup and Catherine Fischer for reviewing the English.

Author Contributions

J.B., L.K., E.W. and C.D. conceived the experiments, J.B., L.K., S.C., N.H., S.M., M.K. conducted the experiments, J.B. analyzed the results. All authors reviewed the manuscript.

Data Availibility

The datasets generated during and/or analyzed during the current study are available from the corresponding author upon reasonable request. The Python script for calculating the time-intensity curves can be found at https://github.com/burkardj/TI-Processing/.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. BVL. Erklärung von Mailand 2019–2024. Tech. Rep. 27.08.2019 (2019).
2. Australian Bureau of Statistics. Australian health survey : Consumption of Added Sugars. Tech. Rep. (2016).
3. Roberts, C. et al. National Diet and Nutrition Survey : Results from Years 7 and 8 (combined) of the Rolling Programme (2014/2015 - 2015/2016). Public Health England 8 (2018).
4. Langlois K Garriguet D Gonzalez A Sinclair S Colapinto CK Change in total sugars consumption among Canadian children and adults Health Rep. 2019 30 10 19 30649778
Langlois, K., Garriguet, D., Gonzalez, A., Sinclair, S. & Colapinto, C. K. Change in total sugars consumption among Canadian children and adults. Health Rep. 30, 10–19 (2019).30649778
5. IUNA. Irish Universities Nutrition Alliance (IUNA) National Teens’ Food Survey II NTFS II (2019-2020). Main Survey Rep. Tech. Rep. (2022).
6. Walton J Bell H Re R Nugent AP Current perspectives on global sugar consumption: Definitions, recommendations, population intakes, challenges and future direction Nutrit. Res. Rev. 2023 36 1 22 10.1017/S095442242100024X 34369326
Walton, J., Bell, H., Re, R. & Nugent, A. P. Current perspectives on global sugar consumption: Definitions, recommendations, population intakes, challenges and future direction. Nutrit. Res. Rev. 36, 1–22. 10.1017/S095442242100024X (2023).34369326 10.1017/S095442242100024X
7. Piernas C Ng SW Popkin B Trends in purchases and intake of foods and beverages containing caloric and low-calorie sweeteners over the last decade in the United States Pediatr. Obes. 2013 8 294 306 10.1111/j.2047-6310.2013.00153.x 23529974
Piernas, C., Ng, S. W. & Popkin, B. Trends in purchases and intake of foods and beverages containing caloric and low-calorie sweeteners over the last decade in the United States. Pediatr. Obes. 8, 294–306. 10.1111/j.2047-6310.2013.00153.x (2013).23529974 10.1111/j.2047-6310.2013.00153.x
8. Hutchings SC Low JY Keast RS Sugar reduction without compromising sensory perception. An impossible dream? Crit. Rev. Food Sci. Nutr. 2019 59 2287 2307 10.1080/10408398.2018.1450214 29561168
Hutchings, S. C., Low, J. Y. & Keast, R. S. Sugar reduction without compromising sensory perception. An impossible dream?. Crit. Rev. Food Sci. Nutr. 59, 2287–2307. 10.1080/10408398.2018.1450214 (2019).29561168 10.1080/10408398.2018.1450214
9. Koeferli CR Piccinali P Sigrist S The influence of fat, sugar and non-fat milk solids on selected taste, flavor and texture parameters of a vanilla ice-cream Food Qual. Prefer. 1996 7 69 79 10.1016/0950-3293(95)00038-0
Koeferli, C. R., Piccinali, P. & Sigrist, S. The influence of fat, sugar and non-fat milk solids on selected taste, flavor and texture parameters of a vanilla ice-cream. Food Qual. Prefer. 7, 69–79. 10.1016/0950-3293(95)00038-0 (1996).10.1016/0950-3293(95)00038-0
10. Pareyt B The role of sugar and fat in sugar-snap cookies: Structural and textural properties J. Food Eng. 2009 90 400 408 10.1016/j.jfoodeng.2008.07.010
Pareyt, B. et al. The role of sugar and fat in sugar-snap cookies: Structural and textural properties. J. Food Eng. 90, 400–408. 10.1016/j.jfoodeng.2008.07.010 (2009).10.1016/j.jfoodeng.2008.07.010
11. Afoakwa EO Paterson A Fowler M Factors influencing rheological and textural qualities in chocolate—a review Trends Food Sci. Technol. 2007 18 290 298 10.1016/j.tifs.2007.02.002
Afoakwa, E. O., Paterson, A. & Fowler, M. Factors influencing rheological and textural qualities in chocolate—a review. Trends Food Sci. Technol. 18, 290–298. 10.1016/j.tifs.2007.02.002 (2007).10.1016/j.tifs.2007.02.002
12. Zorn S Alcaire F Vidal L Giménez A Ares G Application of multiple-sip temporal dominance of sensations to the evaluation of sweeteners Food Qual. Prefer. 2014 36 135 143 10.1016/j.foodqual.2014.04.003
Zorn, S., Alcaire, F., Vidal, L., Giménez, A. & Ares, G. Application of multiple-sip temporal dominance of sensations to the evaluation of sweeteners. Food Qual. Prefer. 36, 135–143. 10.1016/j.foodqual.2014.04.003 (2014).10.1016/j.foodqual.2014.04.003
13. Du Bois GE Prakash I Non-caloric sweeteners, sweetness modulators, and sweetener enhancers Annu. Rev. Food Sci. Technol. 2012 3 353 380 10.1146/annurev-food-022811-101236 22224551
Du Bois, G. E. & Prakash, I. Non-caloric sweeteners, sweetness modulators, and sweetener enhancers. Annu. Rev. Food Sci. Technol. 3, 353–380. 10.1146/annurev-food-022811-101236 (2012).22224551 10.1146/annurev-food-022811-101236
14. Reyes MM Castura JC Hayes JE Characterizing dynamic sensory properties of nutritive and nonnutritive sweeteners with temporal check-all-that-apply J. Sensory Stud. 2017 10.1111/joss.12270
Reyes, M. M., Castura, J. C. & Hayes, J. E. Characterizing dynamic sensory properties of nutritive and nonnutritive sweeteners with temporal check-all-that-apply. J. Sensory Stud.10.1111/joss.12270 (2017).10.1111/joss.12270
15. Wee M Tan V Forde C A comparison of psychophysical dose-response behaviour across 16 sweeteners Nutrients 2018 10 1 16 10.3390/nu10111632
Wee, M., Tan, V. & Forde, C. A comparison of psychophysical dose-response behaviour across 16 sweeteners. Nutrients 10, 1–16. 10.3390/nu10111632 (2018).10.3390/nu10111632
16. Mosca, A. C., Velde, F. v. d., Bult, J. H., van Boekel, M. A. & Stieger, M. Enhancement of sweetness intensity in gels by inhomogeneous distribution of sucrose. Food Qual. Pref. 21, 837–842, 10.1016/j.foodqual.2010.04.010 (2010).
17. Mosca AC van de Velde F Bult JH van Boekel MA Stieger M Taste enhancement in food gels: Effect of fracture properties on oral breakdown, bolus formation and sweetness intensity Food Hydrocolloids 2015 43 794 802 10.1016/j.foodhyd.2014.08.009
Mosca, A. C., van de Velde, F., Bult, J. H., van Boekel, M. A. & Stieger, M. Taste enhancement in food gels: Effect of fracture properties on oral breakdown, bolus formation and sweetness intensity. Food Hydrocolloids 43, 794–802. 10.1016/j.foodhyd.2014.08.009 (2015).10.1016/j.foodhyd.2014.08.009
18. Holm K Wendin K Hermansson AM Sweetness and texture perceptions in structured gelatin gels with embedded sugar rich domains Food Hydrocolloids 2009 23 2388 2393 10.1016/j.foodhyd.2009.06.016
Holm, K., Wendin, K. & Hermansson, A. M. Sweetness and texture perceptions in structured gelatin gels with embedded sugar rich domains. Food Hydrocolloids 23, 2388–2393. 10.1016/j.foodhyd.2009.06.016 (2009).10.1016/j.foodhyd.2009.06.016
19. Burseg KMM Brattinga C de Kok PMT Bult JHF Sweet taste enhancement through pulsatile stimulation depends on pulsation period not on conscious pulse perception Physiol. Behav. 2010 100 327 331 10.1016/j.physbeh.2010.03.007 20307558
Burseg, K. M. M., Brattinga, C., de Kok, P. M. T. & Bult, J. H. F. Sweet taste enhancement through pulsatile stimulation depends on pulsation period not on conscious pulse perception. Physiol. Behav. 100, 327–331. 10.1016/j.physbeh.2010.03.007 (2010).20307558 10.1016/j.physbeh.2010.03.007
20. Burseg KMM Camacho S Knoop J Bult JHF Sweet taste intensity is enhanced by temporal fluctuation of aroma and taste, and depends on phase shift Physiol. Behav. 2010 101 726 730 10.1016/j.physbeh.2010.08.014 20800075
Burseg, K. M. M., Camacho, S., Knoop, J. & Bult, J. H. F. Sweet taste intensity is enhanced by temporal fluctuation of aroma and taste, and depends on phase shift. Physiol. Behav. 101, 726–730. 10.1016/j.physbeh.2010.08.014 (2010).20800075 10.1016/j.physbeh.2010.08.014
21. Burseg KMM Camacho S Franciscus JH Effects of pulsation rate and viscosity on pulsation-induced taste enhancement: New insights into texture—Taste interactions J. Agric. Food Chem. 2011 59 5548 5553 10.1021/jf2002848 21510707
Burseg, K. M. M., Camacho, S. & Franciscus, J. H. Effects of pulsation rate and viscosity on pulsation-induced taste enhancement: New insights into texture—Taste interactions. J. Agric. Food Chem. 59, 5548–5553. 10.1021/jf2002848 (2011).21510707 10.1021/jf2002848
22. Kistler, T., Pridal, A., Bourcet, C. & Denkel, C. Modulation of sweetness perception in confectionary applications. Food Qual. Pref.10.1016/j.foodqual.2020.104087 (2021).
23. Burkard, J. Inkjet-based surface structuring : amplifying sweetness perception through additive manufacturing in foods. Sci. Food, 10.1038/s41538-023-00218-x (2023).
24. Khemacheevakul K Wolodko J Nguyen H Wismer W Temporal sensory perceptions of sugar-reduced 3D printed chocolates Foods 2021 10 1 14 10.3390/foods10092082
Khemacheevakul, K., Wolodko, J., Nguyen, H. & Wismer, W. Temporal sensory perceptions of sugar-reduced 3D printed chocolates. Foods 10, 1–14. 10.3390/foods10092082 (2021).10.3390/foods10092082
25. De Celis Alonso, B. et al. Functional magnetic resonance imaging assessment of the cortical representation of oral viscosity. J. Texture Stud. 38, 725–737, 10.1111/j.1745-4603.2007.00122.x (2007).
26. Farzanmehr H Abbasi S Effects of inulin and bulking agents on some physicochemical, textural and sensory properties of milk chocolate J. Texture Stud. 2009 40 536 553 10.1111/j.1745-4603.2009.00196.x
Farzanmehr, H. & Abbasi, S. Effects of inulin and bulking agents on some physicochemical, textural and sensory properties of milk chocolate. J. Texture Stud. 40, 536–553. 10.1111/j.1745-4603.2009.00196.x (2009).10.1111/j.1745-4603.2009.00196.x
27. Aidoo RP Appah E Van Dewalle D Afoakwa EO Dewettinck K Functionality of inulin and polydextrose as sucrose replacers in sugar-free dark chocolate manufacture—effect of fat content and bulk mixture concentration on rheological, mechanical and melting properties Int. J. Food Sci. Technol. 2017 52 282 290 10.1111/ijfs.13281
Aidoo, R. P., Appah, E., Van Dewalle, D., Afoakwa, E. O. & Dewettinck, K. Functionality of inulin and polydextrose as sucrose replacers in sugar-free dark chocolate manufacture—effect of fat content and bulk mixture concentration on rheological, mechanical and melting properties. Int. J. Food Sci. Technol. 52, 282–290. 10.1111/ijfs.13281 (2017).10.1111/ijfs.13281
28. Wille RL Lutton ES Polymorphism of cocoa butter J. Am. Oil Chem. Soc. 1966 43 491 496 10.1007/BF02641273 5945032
Wille, R. L. & Lutton, E. S. Polymorphism of cocoa butter. J. Am. Oil Chem. Soc. 43, 491–496. 10.1007/BF02641273 (1966).5945032 10.1007/BF02641273
29. Breen SP Etter NM Ziegler GR Hayes JE Oral somatosensatory acuity is related to particle size perception in chocolate Sci. Rep. 2019 9 1 10 10.1038/s41598-019-43944-7 30626917
Breen, S. P., Etter, N. M., Ziegler, G. R. & Hayes, J. E. Oral somatosensatory acuity is related to particle size perception in chocolate. Sci. Rep. 9, 1–10. 10.1038/s41598-019-43944-7 (2019).30626917 10.1038/s41598-019-43944-7
30. Rostagno, W. Chocolate particles size and its organoleptic influence. Manuf. Confect. 81–85 (1969).
31. Beckett, S., Fowler, M. & Ziegler, G. R. Beckett ’ s Industrial Chocolate Manufacture and Use (2017).
32. Kaneko R Kitabatake N Structure-sweetness relationship in thaumatin: Importance of lysine residues Chem. Senses 2001 26 167 177 10.1093/chemse/26.2.167 11238247
Kaneko, R. & Kitabatake, N. Structure-sweetness relationship in thaumatin: Importance of lysine residues. Chem. Senses 26, 167–177. 10.1093/chemse/26.2.167 (2001).11238247 10.1093/chemse/26.2.167
33. Mishra K Controlling lipid crystallization across multiple length scales by directed shear flow J. Colloid Interface Sci. 2023 630 731 741 10.1016/j.jcis.2022.10.005 36274408
Mishra, K. et al. Controlling lipid crystallization across multiple length scales by directed shear flow. J. Colloid Interface Sci. 630, 731–741. 10.1016/j.jcis.2022.10.005 (2023).36274408 10.1016/j.jcis.2022.10.005
34. Aidoo RP Clercq ND Afoakwa EO Dewettinck K Optimisation of processing conditions and rheological properties using Stephan mixer as Conche in small-scale chocolate processing Int. J. Food Sci. Technol. 2014 49 740 746 10.1111/ijfs.12360
Aidoo, R. P., Clercq, N. D., Afoakwa, E. O. & Dewettinck, K. Optimisation of processing conditions and rheological properties using Stephan mixer as Conche in small-scale chocolate processing. Int. J. Food Sci. Technol. 49, 740–746. 10.1111/ijfs.12360 (2014).10.1111/ijfs.12360
35. Savage GP McNeil DL Dutta PC Lipid composition and oxidative stability of oils in hazelnuts (Corylus avellana L.) grown in New Zealand JAOCS, J. Am. Oil Chemists’ Soc. 1997 74 755 759 10.1007/s11746-997-0214-x
Savage, G. P., McNeil, D. L. & Dutta, P. C. Lipid composition and oxidative stability of oils in hazelnuts (Corylus avellana L.) grown in New Zealand. JAOCS, J. Am. Oil Chemists’ Soc. 74, 755–759. 10.1007/s11746-997-0214-x (1997).10.1007/s11746-997-0214-x
36. Rothkopf I Danzl W Changes in chocolate crystallization are influenced by type and amount of introduced filling lipids Eur. J. Lipid Sci. Technol. 2015 117 1714 1721 10.1002/ejlt.201400552
Rothkopf, I. & Danzl, W. Changes in chocolate crystallization are influenced by type and amount of introduced filling lipids. Eur. J. Lipid Sci. Technol. 117, 1714–1721. 10.1002/ejlt.201400552 (2015).10.1002/ejlt.201400552
37. Gavião MBD Engelen L Van Der Bilt A Chewing behavior and salivary secretion Eur. J. Oral Sci. 2004 112 19 24 10.1111/j.0909-8836.2004.00105.x 14871189
Gavião, M. B. D., Engelen, L. & Van Der Bilt, A. Chewing behavior and salivary secretion. Eur. J. Oral Sci. 112, 19–24. 10.1111/j.0909-8836.2004.00105.x (2004).14871189 10.1111/j.0909-8836.2004.00105.x
38. Andablo-Reyes E 3D biomimetic tongue-emulating surfaces for tribological applications ACS Appl. Mater. Interfaces. 2020 12 49371 49385 10.1021/acsami.0c12925 33105986
Andablo-Reyes, E. et al. 3D biomimetic tongue-emulating surfaces for tribological applications. ACS Appl. Mater. Interfaces. 12, 49371–49385. 10.1021/acsami.0c12925 (2020).33105986 10.1021/acsami.0c12925
39. Kokini JL The physical basis of liquid food texture and texture-taste interactions J. Food Eng. 1987 6 51 81 10.1016/0260-8774(87)90021-5
Kokini, J. L. The physical basis of liquid food texture and texture-taste interactions. J. Food Eng. 6, 51–81. 10.1016/0260-8774(87)90021-5 (1987).10.1016/0260-8774(87)90021-5
40. Stark T Bareuther S Hofmann T Molecular definition of the taste of roasted cocoa nibs (Theobroma cacao) by means of quantitative studies and sensory experiments J. Agric. Food Chem. 2006 54 5530 5539 10.1021/jf0608726 16848542
Stark, T., Bareuther, S. & Hofmann, T. Molecular definition of the taste of roasted cocoa nibs (Theobroma cacao) by means of quantitative studies and sensory experiments. J. Agric. Food Chem. 54, 5530–5539. 10.1021/jf0608726 (2006).16848542 10.1021/jf0608726
41. Guinard JX Mazzucchelli R Effects of sugar and fat on the sensory properties of milk chocolate: Descriptive analysis and instrumental measurements J. Sci. Food Agric. 1999 79 1331 1339 10.1002/(SICI)1097-0010(199908)79:11<1331::AID-JSFA365>3.0.CO;2-4
Guinard, J. X. & Mazzucchelli, R. Effects of sugar and fat on the sensory properties of milk chocolate: Descriptive analysis and instrumental measurements. J. Sci. Food Agric. 79, 1331–1339. 10.1002/(SICI)1097-0010(199908)79:11<1331::AID-JSFA365>3.0.CO;2-4 (1999).10.1002/(SICI)1097-0010(199908)79:11<1331::AID-JSFA365>3.0.CO;2-4
42. Mosca AC Rocha JA Sala G van de Velde F Stieger M Inhomogeneous distribution of fat enhances the perception of fat-related sensory attributes in gelled foods Food Hydrocolloids 2012 27 448 455 10.1016/j.foodhyd.2011.11.002
Mosca, A. C., Rocha, J. A., Sala, G., van de Velde, F. & Stieger, M. Inhomogeneous distribution of fat enhances the perception of fat-related sensory attributes in gelled foods. Food Hydrocolloids 27, 448–455. 10.1016/j.foodhyd.2011.11.002 (2012).10.1016/j.foodhyd.2011.11.002
43. Burkard J Nain Shah A Harms E Denkel C Impact of spatial distribution on the sensory properties of multiphase 3D-printed food configurations Food Qual. Prefer. 2023 108 104850 10.1016/j.foodqual.2023.104850
Burkard, J., Nain Shah, A., Harms, E. & Denkel, C. Impact of spatial distribution on the sensory properties of multiphase 3D-printed food configurations. Food Qual. Prefer. 108, 104850. 10.1016/j.foodqual.2023.104850 (2023).10.1016/j.foodqual.2023.104850
44. Wu S Burns SA Reeves A Elsner AE Flicker brightness enhancement and visual nonlinearity Vision. Res. 1996 36 1573 1583 10.1016/0042-6989(95)00226-X 8759460
Wu, S., Burns, S. A., Reeves, A. & Elsner, A. E. Flicker brightness enhancement and visual nonlinearity. Vision. Res. 36, 1573–1583. 10.1016/0042-6989(95)00226-X (1996).8759460 10.1016/0042-6989(95)00226-X
45. Hallock RM Di Lorenzo PM Temporal coding in the gustatory system Neurosci. Biobehav. Rev. 2006 30 1145 1160 10.1016/j.neubiorev.2006.07.005 16979239
Hallock, R. M. & Di Lorenzo, P. M. Temporal coding in the gustatory system. Neurosci. Biobehav. Rev. 30, 1145–1160. 10.1016/j.neubiorev.2006.07.005 (2006).16979239 10.1016/j.neubiorev.2006.07.005
46. Halpern BP Marowitz LA Taste responses to lick-duration stimuli Brain Res. 1973 57 473 478 10.1016/0006-8993(73)90152-2 4722065
Halpern, B. P. & Marowitz, L. A. Taste responses to lick-duration stimuli. Brain Res. 57, 473–478 (1973).4722065 10.1016/0006-8993(73)90152-2
47. Halpern BP Kelling ST Meiselman HL An analysis of the role of stimulus removal in taste adaptation by means of simulated drinking Physiol. Behav. 1986 36 925 928 10.1016/0031-9384(86)90454-3 3714870
Halpern, B. P., Kelling, S. T. & Meiselman, H. L. An analysis of the role of stimulus removal in taste adaptation by means of simulated drinking. Physiol. Behav. 36, 925–928. 10.1016/0031-9384(86)90454-3 (1986).3714870 10.1016/0031-9384(86)90454-3
48. Theunissen MJ Kroeze JH Mouth movements diminish taste adaptation, but rate of mouth movement does not affect adaptation Chem. Senses 1996 21 545 551 10.1093/chemse/21.5.545 8902283
Theunissen, M. J. & Kroeze, J. H. Mouth movements diminish taste adaptation, but rate of mouth movement does not affect adaptation. Chem. Senses 21, 545–551. 10.1093/chemse/21.5.545 (1996).8902283 10.1093/chemse/21.5.545
49. Busch JL Tournier C Knoop JE Kooyman G Smit G Temporal contrast of salt delivery in mouth increases salt perception Chem. Senses 2009 34 341 348 10.1093/chemse/bjp007 19273461
Busch, J. L., Tournier, C., Knoop, J. E., Kooyman, G. & Smit, G. Temporal contrast of salt delivery in mouth increases salt perception. Chem. Senses 34, 341–348. 10.1093/chemse/bjp007 (2009).19273461 10.1093/chemse/bjp007
50. Grabenhorst F Rolls ET Bilderbeck A How cognition modulates affective responses to taste and flavor: Top-down influences on the orbitofrontal and pregenual cingulate cortices Cereb. Cortex 2008 18 1549 1559 10.1093/cercor/bhm185 18056086
Grabenhorst, F., Rolls, E. T. & Bilderbeck, A. How cognition modulates affective responses to taste and flavor: Top-down influences on the orbitofrontal and pregenual cingulate cortices. Cereb. Cortex 18, 1549–1559. 10.1093/cercor/bhm185 (2008).18056086 10.1093/cercor/bhm185
51. Nitschke JB Altering expectancy dampens neural response to aversive taste in primary taste cortex Nat. Neurosci. 2006 9 435 442 10.1038/nn1645 16462735
Nitschke, J. B. et al. Altering expectancy dampens neural response to aversive taste in primary taste cortex. Nat. Neurosci. 9, 435–442. 10.1038/nn1645 (2006).16462735 10.1038/nn1645
52. Woods AT Expected taste intensity affects response to sweet drinks in primary taste cortex NeuroReport 2011 22 365 369 10.1097/WNR.0b013e3283469581 21512419
Woods, A. T. et al. Expected taste intensity affects response to sweet drinks in primary taste cortex. NeuroReport 22, 365–369. 10.1097/WNR.0b013e3283469581 (2011).21512419 10.1097/WNR.0b013e3283469581
53. Wilton M Stancak A Giesbrecht T Thomas A Kirkham T Intensity expectation modifies gustatory evoked potentials to sweet taste: Evidence of bidirectional assimilation in early perceptual processing Psychophysiology 2019 56 1 12 10.1111/psyp.13299
Wilton, M., Stancak, A., Giesbrecht, T., Thomas, A. & Kirkham, T. Intensity expectation modifies gustatory evoked potentials to sweet taste: Evidence of bidirectional assimilation in early perceptual processing. Psychophysiology 56, 1–12. 10.1111/psyp.13299 (2019).10.1111/psyp.13299
54. Mellers B Birnbaum M Loci of contextual effects in judgement J. Exp. Psychol. 1982 4 582 601 10.1037/0096-1523.8.4.582
Mellers, B. & Birnbaum, M. Loci of contextual effects in judgement. J. Exp. Psychol. 4, 582–601. 10.1037/0096-1523.8.4.582 (1982).10.1037/0096-1523.8.4.582
55. Le-Bail A Maniglia BC Le-Bail P Recent advances and future perspective in additive manufacturing of foods based on 3D printing Curr. Opin. Food Sci. 2020 35 54 64 10.1016/j.cofs.2020.01.009
Le-Bail, A., Maniglia, B. C. & Le-Bail, P. Recent advances and future perspective in additive manufacturing of foods based on 3D printing. Curr. Opin. Food Sci. 35, 54–64. 10.1016/j.cofs.2020.01.009 (2020).10.1016/j.cofs.2020.01.009
56. M’Barki A Bocquet L Stevenson A Linking rheology and printability for dense and strong ceramics by direct ink writing Sci. Rep. 2017 7 1 10 10.1038/s41598-017-06115-0 28127051
M’Barki, A., Bocquet, L. & Stevenson, A. Linking rheology and printability for dense and strong ceramics by direct ink writing. Sci. Rep. 7, 1–10. 10.1038/s41598-017-06115-0 (2017).28127051 10.1038/s41598-017-06115-0
57. Biswas N Cheow YL Tan CP Siow LF Physical, rheological and sensorial properties, and bloom formation of dark chocolate made with cocoa butter substitute (CBS) Lwt 2017 82 420 428 10.1016/j.lwt.2017.04.039
Biswas, N., Cheow, Y. L., Tan, C. P. & Siow, L. F. Physical, rheological and sensorial properties, and bloom formation of dark chocolate made with cocoa butter substitute (CBS). Lwt 82, 420–428. 10.1016/j.lwt.2017.04.039 (2017).10.1016/j.lwt.2017.04.039
58. Glicerina V Balestra F Rosa MD Romani S Rheological, textural and calorimetric modifications of dark chocolate during process J. Food Eng. 2013 119 173 179 10.1016/j.jfoodeng.2013.05.012
Glicerina, V., Balestra, F., Rosa, M. D. & Romani, S. Rheological, textural and calorimetric modifications of dark chocolate during process. J. Food Eng. 119, 173–179. 10.1016/j.jfoodeng.2013.05.012 (2013).10.1016/j.jfoodeng.2013.05.012
59. Mantihal S Prakash S Godoi FC Bhandari B Effect of additives on thermal, rheological and tribological properties of 3D printed dark chocolate Food Res. Int. 2019 119 161 169 10.1016/j.foodres.2019.01.056 30884645
Mantihal, S., Prakash, S., Godoi, F. C. & Bhandari, B. Effect of additives on thermal, rheological and tribological properties of 3D printed dark chocolate. Food Res. Int. 119, 161–169. 10.1016/j.foodres.2019.01.056 (2019).30884645 10.1016/j.foodres.2019.01.056
60. de Boer GB de Weerd C Thoenes D Goossens HW Laser diffraction spectrometry: Fraunhofer diffraction versus mie scattering Particle Particle Syst. Character. 1987 4 14 19 10.1002/ppsc.19870040104
de Boer, G. B., de Weerd, C., Thoenes, D. & Goossens, H. W. Laser diffraction spectrometry: Fraunhofer diffraction versus mie scattering. Particle Particle Syst. Character. 4, 14–19. 10.1002/ppsc.19870040104 (1987).10.1002/ppsc.19870040104
