Next Article in Journal
Microstructural Evolution and Mechanical Properties of a Ni-Based Alloy with High Boron Content for the Pre-Sintered Preform (PSP) Application
Previous Article in Journal
Recent Advances in Bromine Complexing Agents for Zinc–Bromine Redox Flow Batteries
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transition Temperature of Color Change in Thermochromic Systems and Its Description Using Sigmoidal Models

Department of Material Engineering, Faculty of Textile Engineering, Technical University of Liberec, 46117 Liberec, Czech Republic
*
Author to whom correspondence should be addressed.
Materials 2023, 16(23), 7478; https://doi.org/10.3390/ma16237478
Submission received: 4 November 2023 / Revised: 27 November 2023 / Accepted: 30 November 2023 / Published: 2 December 2023

Abstract

:
Background: Symmetric sigmoidal models with four parameters based on an idealized S/Z-shaped curve are commonly used to analyze the optical parameters of thermochromic materials. However, our experimental findings show that this approach leads to systematic errors involving the incorrect estimation of the transition temperature or the possibility of a virtual indication of the hysteresis nature of a reversible thermochromic change. For this reason, we sought to find a five-parameter model that would appropriately avoid this problem. Methods: Two commercial thermochromic pigments were used for the test and applied to a textile substrate at different concentrations. The optical properties were measured using reflectance spectrophotometry and then converted to Kubelka–Munk function values and colorimetric coordinates. The following statistics were used to assess the quality of the selected sigmoidal models: coefficient of determination, R2; adjusted coefficient of determination, AR2; root mean square error, RMSE; and Akaike Information Criterion, AIC. Results: The four-parameter models were compared with each other and with the five-parameter models using nested F-tests based on residual variance to obtain a statistical measure of superior performance. For all thermochromic color change data examined, the five-parameter models resulted in significantly better fitting. It could be shown that the five-parameter model showed significantly higher accuracy and precision in determining the transition temperature, like non-sigmoidal quantification methods. Conclusions: We concluded that the asymmetric five-parameter model is a valuable extension of the symmetric model in the investigation of thermochromic color changes, providing better parameter estimates and a new approach to investigating the mechanisms contributing to the asymmetry of the thermochromic curve.

1. Introduction

Chromic fabrics represent a specific area of indicators and sensory systems. Color change makes it easy to identify both the different stimuli and their intensity. In the case of thermochromic textiles, this involves temperature measurement or the identification of a specific pre-determined temperature. There are currently two systems that are based on a change of color depending on a specific temperature. The first of these systems is liquid crystals, which have the advantage of covering a temperature range of approximately 20 °C and are typically used as contact thermometers that can be placed on different surfaces [1,2,3,4,5,6,7]. The second system is usually a three-component composite system consisting of Leuco dye, a developer, and a suitable solvent [8,9,10,11,12]. This system can be reduced to a two-component system if the substance that acts as a solvent also acts as a developer [13]. Currently, a wide range of inks is available with a relatively wide temperature range where thermochromic change occurs, typically from −15 °C to 65 °C depending on the application [1,2,14,15,16,17,18]. In general, thermochromic materials can exhibit both irreversible and reversible color changes [18,19], and this paper focuses on the issue of reversible change. The thermochromic changes are schematically illustrated in Figure 1.
In the case of systems that exhibit irreversible color change, they are usually inorganic compounds [20,21], while in the case of reversible systems, they may be both inorganic and organic compounds or composite systems [22,23]. An accurate description of the color transition allows not only for the comparison of different thermochromic pigments but also for measuring possible changes in their behavior due to use. In the case of the use of composite thermochromic inks containing Leuco dye, a developer, and a suitable solvent, it is possible to find that the individual compositions differ from each other by examining the different average color change temperatures (the average temperature at which color change is observed) and the range of temperatures at which color change occurs [24,25].
If a thermochromic system is to be used as an indicator of reaching a certain temperature, it is preferable that the temperature range over which the color transition is achieved is as narrow as possible [26,27]. In addition, it is necessary to study the change in color response of the system during use. Most often, the number of cycles of thermochromic color change that the thermochromic ink can perform is studied. The influence of light, other factors, and effects of use on thermochromic systems are studied [28].
The color change in the case of composite thermochromic inks is caused by the phase change of the substance used as a solvent. This solvent is converted from a solid phase to a liquid phase, causing the developer and Leuco dye to separate and the resulting system to become colorless. In other words, it is possible to observe the change of state of the composite system from one state to another, while in published works usually only the record of the experimental points of the color change is evaluated, which has the shape of an S or a Z curve (this depends on the optical parameter measured), or the initial and final phases of the color change, where the corresponding temperatures are determined using a simple graphical procedure [14,29]. The principle of the graphic determination of the initial and final temperatures of the thermochromic color change consists of finding the intersections of a line through the linear part of the color transition data between two states, with the lines corresponding to the initial and final equilibrium of the measured color response of the thermochromic system. These changes can be analytically described using parametric functions such as third-order and higher-order polynomials, cotangential, sigmoidal, or logistical functions [30]. The advantage of parametric functions is the possibility of discovering the physical nature of the described phenomenon and comparing both the quality of the interleaving of experimental data and the study of possible errors associated with the appropriate determination of the parameters such as the initial, final, and mean temperature of the thermochromic transition [31]. As will be shown below, the most appropriate approach for this purpose is the use of sigmoidal models based on the original Sigmoidal Boltzmann Equation (SBE). This new approach to the analysis of thermochromic color change data also allows for reduced errors in the determination of the individual temperatures mentioned above when the analyzed data had an S-shape of color change versus a Z-shape and vice versa.

2. Materials and Methods

Two commercial thermochromic inks from Matsui International, Co., Kyoto, Japan, which are suitable for screen printing, were used to verify the usefulness of the considered sigmoidal function for estimating the thermochromic color curve. Gold Orange and Fast Blue, Type 37, for which Matsui states that the ink color is visible below 32 °C and the color disappears above 41 °C. Concentrations of thermochromic ink: 10, 50, 100, 150, and 300 g·kg−1 were used. A proven complex thickener system containing glycerin, Lukosan S (defoamer from Lučební závody a.s., Kolín, Czech Republic), Sokrat 492 (anionic binder from Chemické závody, Sokolov, Czech Republic), Acramin BA (binder from Bayer AG, Leverkusen, Germany, Lambicol L 90 (thickener from Lamberti S.p.A., Gallarate (VA) Italy), ammonia and water was used as printing paste.

2.1. Printing

The sample was printed on a Johannes Zimmer laboratory printing machine MINI-MD-R/F (Zimmer Maschinenbau GmbH, Klagenfurt, Austria) designed for screen printing. The printing was carried out using a flat screen and a roller squeegee. The squeegee was moved over the stencil using electromagnets placed under the printing blanket. The laboratory printing machine was set up to move the squeegee at a speed of 3 m·min−1 with a pressure of 3 N·cm−1. To minimize overprinting, only one squeegee stroke was carried out. Printing was followed by gradual drying of the samples on filter papers until they were dry to the touch. This was followed by fixation of the samples in an HS 122 hot-air oven at 80 °C for 5 min.

2.2. Measurement of Spectral and Color Parameters

The measurements were designed as isothermal measurements using two heated plates, one of which had a constant temperature of 20 °C using a Minichiller thermostat from Peter Huber Kältemaschinenbau SE, Offenburg, Germany and the second heated plate was gradually adjusted to temperatures of 30 to 45 °C in 1 °C steps using a Julabo F25-HE thermostat (JULABO GmbH, Seelbach, Germany). The actual reflectance measurement was carried out on a Datacolor D100 spectrophotometer (Lawrenceville, NJ, USA, aperture diameter 20 mm, mode: di:8° [32]), which was set to the vertical position, the sample was placed face down on the aperture measurement port and the corresponding heated plate was placed on top of the sample. First, the sample response was measured at 20° and then the thermal plate at a higher temperature was attached, and time measurements were made at intervals of 5 s to determine the kinetics of the color change and to determine the equilibrium of the color change corresponding to a condition where, on ten consecutive measurements, the deviation between measurements was no greater than 0.1 DE *. The result is the color response data for the two thermochromic inks tested and their different concentrations.

2.3. Spectral and Color Data Evaluation

The observed reflectance data were used to determine the Kubelka–Munk function values and tristimulus values XYZ, which represent CIE 1931 color space [33]:
( K / S ) λ = [ 1 β R ( λ ) ] 2 2 β R ( λ )
where βR(λ) is the reflected radiance factor [34], and K and S are Kubelka–Munk absorption and scattering coefficients.
X = k 360 830 ε ( λ ) β R ( λ ) x ¯ ( λ ) d λ
Y = k 360 830 ε ( λ ) β R ( λ ) y ¯ ( λ ) d λ
Z = k 360 830 ε ( λ ) β R ( λ ) z ¯ ( λ ) d λ
where k is the constant chosen such that Y = 100 for object, for which βR(λ) equals 1 for all wavelengths of selected wavelength interval; ε(λ) is spectral power of used light source, typically CIE (CIE—International Commission on Illumination (http://cie.co.at)) D65 illuminant and x ¯ ( λ ) ,   y ¯ ( λ ) , z ¯ ( λ ) are color matching functions of CIE standard observer, in this case, version 1931—2°. These functions are the numerical description of the chromatic response of the CIE standard observer—an average human observer in the range of 360–830 nm.
From these tristimulus values, color coordinates are computed in approximate uniform color space CIELAB [33] with use of following relations:
L * = 116 f ( Y / Y 0 ) 16
a * = 500 [ f ( X / X 0 ) f ( Y / Y 0 ) ]
b * = 200 [ f ( Y / Y 0 ) f ( Z / Z 0 ) ]
where function f equals to (τ/τ0)1/3 if (τ/τ0) > (24/116)3 or (841/108) (τ/τ0) + 16/116 if (τ/τ0) (24/116)3; τ are mentioned tristimulus values X, Y, Z of measured specimen; and τ0 is related to tristimulus values of ideal white diffuser-object for which βR(λ) equals 1 for all wavelengths of selected wavelength interval. The function f ensures that the transformation used from tristimulus values X, Y, Z corresponds as closely as possible to the nature of human color perception, with individual colors approximately evenly scaled.
The coordinates of the CIELAB color space form a rectangular coordinate system; however, a derived cylindrical coordinate system, which we refer to as CIELCH, is preferable in terms of visual color evaluation [33,35]. In this system, in addition to the lightness L*, two other coordinates are used: the chroma C* and the hue angle h°.
C * = ( a * ) 2 + ( b * ) 2
h ° = a r c t g ( b * a * )
Using these coordinates, it is possible to consider the perception of individual differences using the human eye, so that an observer can usually distinguish the difference in brightness (10), chromaticity (11), and hue (12). The sample at the beginning of the measurement is usually used as a standard, and the change in color of the sample as the temperature increases or decreases is considered to be a corresponding batch related to the temperature used.
L * = L b * L s *
C * = C b * C s *
H * = 2 C b * . C s * sin ( h ° 2 )
where subscript b refers to the batch and subscript s to the value relative to the standard.
The advantage of color difference evaluation is that it can be used in cases where the temperature indicator based on thermochromic inks will consist of two parts, one being color stable and the other consisting of a color changing thermochromic system. In this case, the two parts of the temperature indicator appear to be color identical below the temperature of thermochromic change and color different above the temperature of the color change. The human eye is more sensitive to these changes than to the color change of one uniform colored area [36].
In this work, the standard used was the measurement of the test sample of thermochromic ink at 20 °C before the start of the warm-up phase of the test. The measurements that corresponded to the equilibrium color of the sample at the selected test temperature were then used as the batch, as discussed in the previous section.

2.4. Mathematical Modeling

The thermochromic color change in thermochromic inks containing dye and developer is characterized using a sigmoid transition along with an inflection point, as is the case with many other phase transitions [37,38,39,40]. A frequently used function in this case is the sigmoidal function proposed by Boltzmann (1879) [37], which has the following univariate binary form: y = 1/(1 + ex). In our case, we used a modified version that allows us to determine the inflection point of the color transition and temperature distances (the difference between the temperature at which the color change of the thermochromic system begins and the temperature at which it ends). This form is called Sigmoidal Boltzmann Equation (SBE):
y = y b + ( y f y b ) A + e c x α
where yb and yf are equilibrium values of the dependent variable (K/S, ΔL*, ΔC*, ΔH* and ΔE*) before and after transition; c is critical value of stimuli x; in our case transition temperature TT, which is also abscissa of inflection point.
A is a parameter that adjusts the position of the critical value relative to the geometric location of the transition between two states. If A = 1, the critical value is in the middle of the transition, and it is true that yc = (yb + yf)/2. In cases where A is not equal to 1 but not very different from 1 the critical value of the transition is located near the mid [point]. Because in the case of calculating the fit of the experimental data with the model (13) with an alternative value of the parameter A, the calculations proved to be unstable in some parameters, a fixed value of A = 1 was chosen, when the calculations were stable, and the statistical criteria confirmed the validity of the model.
Parameter α describes the properties of the color transition slope and identifies a possible discontinuity in the monitored process. If α is close to zero, then the first derivative of model (13) approaches infinity and discontinuous transition is observed. Parameter α is also a curvature parameter that indicates over what temperature range the color of the thermochromic material changes. In the case of sensing materials, it is important that this range is as narrow as possible. In the case that the transition temperature is not significantly different from the temperature determined as the center of the color transition between the two states of the thermochromic sensor, Equation (13) is applicable for such determinations. The final version of SBE is the following for the Kubelka–Munk function, with a similar solution for the colorimetric parameters ΔL*, ΔC*, ΔH* and ΔE*.
( K / S ) = ( K / S ) b + ( K / S ) f ( K / S ) b 1 + e c T α
The transition point generally assumes that 50% of the thermochromic ink has undergone color transformation. However, in addition to the transit temperature, the temperatures of interest are those showing 10% and 90% transformation of the thermochromic ink, which are usually found at the beginning and end of the approximately linear part of SBE, where the curved parts of the function leading asymptotically to the limiting values end or begin. To address these cases, we can make the following modification of Equation (14), where the temperature of 10% or 90% thermochromic change is denoted as TR.
e c T R α = ( K / S ) f ( K / S ) b ( K / S ) R ( K / S ) b 1
Using simple modifications of Equation (17) we obtain:
c T R α = l n [ ( K / S ) f ( K / S ) b ( K / S ) R ( K / S ) b 1 ]
T R = c α { l n [ ( K / S ) f ( K / S ) b ( K / S ) R ( K / S ) b 1 ] }
Equation (19) has a solution only under the condition that [ ( K S ) f ( K S ) b ( K S ) R ( K S ) b 1 ] > 0 .

3. Results and Discussion

The reflectance values of the measured samples were converted to the Kubelka–Munk function (1) and colorimetric parameters (2–12) according to the above equations. The colorimetric parameters and the CIE Yxy and CIELAB coordinates are shown graphically in the graphs in Figure 2. These graphs show that during the heating of the thermochromic fabrics, the observed samples are gradually decolorized, which is manifested as a decrease in the excitation purity for the CIE xy diagram and in the chromaticity in plane a* b* for the CIELAB color space. At the same time, there is an increase in brightness Y and lightness L*. In the case of the a* b* plane of the CIELAB color space, it is also seen that the decoloring process follows a curve and not a straight line as one might assume. The reason for this behavior is that the concentration dependence of the sample position in the color space is affected in real colorants both by the substrate itself, which is not ideally achromatic, and by the change in the spectral reflectance waveform, where, after a certain minimum reflectance is reached, with a further increase in concentration there is a decrease in reflectance around the dominant wavelength of the colorant, which ultimately results in a decrease in chroma and a certain hue shift. If we express the values of the Kubelka–Munk function (K/S) as a function of temperature, we obtain a typical sigmoidal curve for thermochromic color change with well resolved individual concentrations of the tested thermochromic samples, as can be seen in the plots in Figure 3.
At the same time, one can see the fit of the experimental points using a symmetric SBE. Table 1 and Table 2 show the results for SBE (14) for both tested inks, individual concentrations, and estimates of the temperatures at which 50%, 10% and 90% decolorization of the tested inks occurs. The estimates of temperatures TR10 and TR90 are then used to determine the basic temperature interval in which a significant thermochromic color change occurs. Its width, denoted as temperature distance ΔT, then indicates the suitability of the thermochromic system in terms of its sensory application. As already mentioned, the aim is to keep the temperature distance in which the color change occurs as narrow as possible. Both inks approximately fall into the manufacturer’s declared category referred to as 37, where the temperature range can be 32 °C up to 41 °C. (https://matsui-color.com/ourproducts/specialty-products/chromicolor, accessed 25 November 2023). However, their temperature distances are differently sized and differently dependent on concentration, i.e., ΔT depends both on the ink used and, above all, on its concentration.
From the graph in Figure 4A, the thermochromic ink ‘Gold Orange’ has a significantly higher dependence of ΔT on temperature (essentially an exponential decrease approaching 2.5 °C) compared to the ink ‘Fast Blue’. From the plot in Figure 4B, it is also evident that the temperature distance ΔT is the inverse of the slope coefficient α. Expressing the thermal distance ΔT from Equation (17) for the case that ΔT = TR90TR10, we obtain Equation (18):
T = 2 α | l n [ ( K / S ) f ( K / S ) b ( K / S ) R ( K / S ) b 1 ] |
Equation (18) has the following solution, since the difference (K/S)f (K/S)b represents 100% of the color change of the thermochromic ink. If we consider the difference of the measured values at 10% conversion of the thermochromic ink as the situation where the difference (K/S)R (K/S)b represents 90% and at 90% conversion as the difference when (K/S)R (K/S)b represents 10%, then we obtain the following solution of Equation (18):
Conversion   10 % :   T = 2 α | l n [ 100 90 1 ] | = 2 α × | 2.1860178 | = 2 α × 2.1860178
Conversion   90 % :   T = 2 α | l n [ 100 10 1 ] | = 2 α × | 2.1860178 | = 2 α × 2.1860178
And Equation (18) can be rewritten into a simple form:
T = 2 α . ln ( 10 9 1 )
It was mentioned that the solution of Equations (14)(20) given for the Kubelka–Munk function (K/S) can also be applied on colorimetric color parameters (ΔL*, ΔC*, ΔH*, and ΔE*). Their advantage is that they can be measured both with simple colorimeters and, for example, with digital cameras. Moreover, this evaluation is closer to the perception of the human eye, which is important in terms of the application of thermochromic systems as simple visual indicators of a certain temperature. In addition to the asymmetric progression of the individual difference curves, we can also note the difference in the contribution of the partial differences (ΔL*, ΔC*, and ΔH*) to the overall color difference between the compared inks, Figure 5 and Figure 6. In the case of the Fast Blue ink, we see that the contribution of the difference in lightness ΔL* has the highest influence on the total color difference ΔE*, then for the Gold Orange ink we see that the main component of the total color difference ΔE* is the difference in chroma ΔC*.
An interesting pattern in the measured data can also be seen in the hue difference ΔH*, which reaches a maximum in the transition temperature area and then returns to the minimum difference. Also, in the case of the Fast Blue ink, we can see that the proportion of the hue difference has a relatively small effect on the total color difference ΔE*. From this, we can conclude that the evaluation of the transit temperature from the progression of the hue difference dependence will be the least effective and, therefore, in the next section, we will concentrate on the analysis of the progression of the remaining differences, i.e., ΔE*, ΔL*, and ΔC*.

3.1. Asymmetric Models of Sigmoidal Fit

Over the years, several asymmetric sigmoidal functions have been developed. The first tested equation was published by Gompertz [41] and is an example of a mathematical model for time series; however, it can also be used for the case of general sigmoidal trends including the color change of the thermochromic ink developer tested in this work. All the following equations can also be written for ΔL* and ΔC*.
E * = E b * + ( E f * E b * ) e e ( c T α )
The Gompertz equation is a special case of a generalized logistical function that is also known as the Richards [42] function, which was included as a second asymmetric model:
E * = E b * + ( E f * E b * ) [ 1 + ν e ( c T α ) ] 1 ν
where ν is constant of asymmetry.
The next pair of equations are the Double Sigmoidal Boltzmann Equations DSBE [43]. The first of these, denoted as DSBE1, has the following form:
E * = E b * + E f * [ p 1 + e ( c A T α ) + 1 p 1 + e ( c B T β ) ]
where β is the second curvature parameter and CA/CB are the critical stimulus values (the temperatures related to the first and second part of the curvature of the transition function).
The second Double Sigmoidal Boltzmann Equation, denoted as DSBE2 [44], is written as follows:
E * = E b * 1 + e ( c A T α ) + E f * E b * 1 + e ( c B T β )
The last model tested was BARO5, which is a five-parameter model extended from the four-parameter SBE model with the addition of a second curvature parameter. This model is used to assess arterial baroreflex using vasoactive agents allowing investigators to collect pairs of data over a wide range of blood pressures and reflex responses [45].
E * = E b * + E f * E b * 1 + f e ( c T α ) + ( 1 f ) e ( c T β )
where f defines a logistic weighting function varying smoothly between 0 and 1 centered around c (V50), and if α and β are of the same sign, the mean curvature f is given by:
f = 1 1 + e ( c T γ )
γ = 2 α β | α + β |

3.2. Assessment of the Best Model Using Statistical Tests

The most common way of comparing the quality of equations modelling a particular dependence is to use the coefficient of determination R2 [46], which is calculated according to the following equation:
R 2 = 1 R S S T S S
with RSS being the residual sum of squares R S S = i = 1 n ( y i y ^ i ) 2 and TSS the total sum of squares T S S = i = 1 n ( y i y ¯ i ) 2 , where n represents number of observations.
Ideally, RSS = 0 and the coefficient of determination reaches 1. Except for Gompertz’s model, the other tested models of the asymmetric sigmoidal function are nested models, i.e., they are extensions of the original four parameters of SBE (13) using an additional parameter or a doubling. As a result, the F-test, which quantifies the relative decrease in the sum of squares when moving from a simpler model to a more general model, can be used to select the best fitting model.
F = ( R S S 4 R S S 5 ) ( d f 4 d f 5 ) R S S 5 d f 5
where RSS4 is the residual sum of squares of a four-parameter fit; RSS5 is the residual sum of squares of a five-parameter fit; and df4 and df5 are the degrees of freedom for a four- and five-parameter fit.
Another way to test the relative quality of the individual models is the Akaike Information Criterion (AIC). This criterion tests the relative quality of a model, i.e., in relation to other models. The lower the AIC, the better the model fits the data and the higher its ranking [47]. For the case of small data sets, which is also the case for thermochromic color change measurements with a 1 °C step, the recommended AIC form is as follows:
A I C = n ln ( R S S i n ) + 2 k i + 2 k i ( k i + 1 ) n k i 1 ,   n k i + 1
where ki is the number of model parameters plus one.
The AIC primarily rewards improvement in the goodness of fit to the experimental data, while to account for increased uncertainty, the AIC also includes a penalty that is an increasing function of the number of parameters.
Like the AIC, the Adjusted Coefficient of Determination AR2 also handicaps the increase in the number of parameters in the model when assessing the quality of the fit to the experimental data, and its formula is as follows:
A R 2 = 1 ( 1 R 2 )   n 1 n p 1
where p is the number of parameters.
The last criterion used in this paper is the root-mean-squared-error RMSE, which allows comparing different models using absolute RMSE values, where the smaller the RMSE, the better the tested model fit.
A R 2 = 1 ( 1 R 2 )   n 1 n p 1
All asymmetric models were tested against the SBE (except for the Gompertz model), i.e., in a nested F-test, the asymmetric model match was considered to be an alternative hypothesis (H1) to the null hypothesis (H0), i.e., the SBE match. Thus, the best performing model is indicated in bold in the table. The tables below also allow for a comparison of all the models using all the statistical criteria mentioned.
For readability, only the results for the total color difference ΔE* and one concentration of thermochromic inks are shown here, the other results are available from the authors upon request. All tests and parameter estimates were performed using GraphPad Prism 10 statistical software for macOS (version 10.1.1) and were also cross-checked using OriginPro 2018 on the Microsoft Windows platform (version b.9.5.1.195).
One of the reasons for testing asymmetric models of thermochromic color change in this work is the difference in transition temperature estimates using spectrophotometric data and colorimetric evaluation. Note in Table 3 and Table 4 that the transition temperature estimate TT by the symmetric SBE model is 38.0 °C for the total color difference ΔE* and the Fast Blue thermochromic ink, which is 1.3 °C higher than the TT estimate using the Kubelka–Munk function (Table 1). In the case of the Gold Orange (Table 2) ink, the difference between the two estimates of TT is lower, but here again the TT estimate using the total color difference ΔE* is 0.8 °C higher. This is a problem, of course, because in the case of one original data set of the reflected radiance factor, we expect the estimate of the transition temperature to be the same for both the Kubelka–Munk and the colorimetric data, since both are based on this data set.
Table 3 and Table 4 show that the DSBE fit using the Double Sigmoidal Boltzmann Equations was the most effective (the first in order is marked in bold). In the case of the Fast Blue thermochromic ink, this was the DSBE2 model, whereas in the case of the Gold Orange thermochromic ink it was the DSBE1 model. In the case of the DSBE2 and the Fast Blue ink, we see that this model was the best performing in terms of AIC and AR2. In terms of comparison using R2 and RSME, this model was second in order. For this reason, a comparison using the F-test was made between this model, DSBE2, and BARO5.
The difference between the TT estimates decreases to 1 °C for the Fast Blue thermochromic ink and to 0.3 °C for the Gold Orange thermochromic ink. Thus, the question remains why the colorimetric data lead to a higher temperature estimate than the Kubelka–Munk function. In the plots in Figure 7, we see that the fit to the experimental data is tighter in the case of the asymmetric sigmoidal models than in the case of the SBE. Thus, the same comparison needs to be made in the case of the fit for the Kubelka–Munk function (K/S) data. The results of the evaluation of each model for both pigments are summarized in Table 5 and Table 6. It is evident that the use of asymmetric models results in a shift in the TT estimation compared to the SBE. For the Fast Blue thermochromic ink, model DSBE1 was evaluated as the most effective for (K/S), in contrast to the total color difference data where it was model DSBE2. In the case of both inks, we can see that the data for model DSBE2 is not filled because the calculation of the estimation of some parameters was unstable.
From this it can be concluded that this model is suitable for the case of S-curves, not for Z-curves such as the curve (K/S). For the Gold Orange thermochromic ink, BARO5 was judged to be the best in case (K/S), but the difference in performance between this model and DSBE1 is very small. Similarly, we can evaluate the ranking of the individual models for the Fast Blue thermochromic ink and the overall color difference data, where DSBE2 was evaluated as the best performing model, but the results of DSBE1 are similar. The fit of the experimental data to the Kubelka–Munk function and the color difference ΔΕ* using the DSBE1 model is shown in the plots in Figure 8. The fit of the experimental data with the DSBE1 model is tighter than that of the four-parameter SBE model. This makes up for the difference between the adaptation of the spectral data (K/S) and the colorimetric data ΔΕ*. On the other hand, it should be considered that the double sigmoid models provide estimates of two temperatures, which may differ significantly from each other, as demonstrated in the case of the Fast Blue ink, where the difference is 2.4 °C. This is due to the nature of the data, which for the curve (K/S) leads to two displaced areas for the calculation of the inflection point, as documented in the graph in Figure 8.

4. Conclusions

The results discussed showed that the symmetric Sigmoidal Boltzmann Equation SBE may be a source of underestimating the position of the inflection point in case of descending Z-curves. Conversely, in the case of ascending S-type sigmoidal curves, the position of the inflection point is overestimated. As a result, there is a difference in the estimates of the transit temperature TT thermochromic color change determined using spectrophotometric and colorimetric data. Among the set of asymmetric sigmoidal models, the Double Sigmoidal Boltzmann Equation DSBE1 was the best performing for the evaluation of the Kubelka–Munk function data and the overall color difference. The resulting difference between the transition temperatures determined using the two mentioned data sets decreased from 1.3 °C to 0.4 °C for the Fast Blue thermochromic ink and from 0.8 °C to 0.2 °C for the Gold Orange thermochromic ink. A certain drawback of the DSBE1 method is the estimation of two critical temperatures, which can be influenced by local fluctuations in the measured data, resulting in a secondary inflection region of the curve that corresponds to the data. In such cases, the BARO5 model, which unifies the calculation of the critical temperature of the thermochromic color change for both asymmetric parts of the sigmoidal transition, seems preferable.

Author Contributions

M.V. (Martina Viková) was responsible for conceptualization and writing the original draft, M.V. (Michal Vik) performed the investigation and data analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Ministry of Education, Youth and Sports of the Czech Republic and the European Union—European Structural and Investment Funds in the framework of the Operational Program Research, Development and Education—project Hybrid Materials for Hierarchical Structures (HyHi, Reg. No. CZ.02.1.01/0.0/0.0/16_019/0000843).

Data Availability Statement

All experimental data are available from the authors.

Acknowledgments

We would like to thank Marcela Pechová, Tereza Pocová, and Branislav Budoš for their help in the preparation of the concentration lines of the thermochromic inks for the printed textiles and measurement.

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationship that could have influenced the work reported in this paper.

References

  1. Bamfield, P.; Hutchings, M. Chromic Phenomena, 3rd ed.; Royal Society of Chemistry: London, UK, 2018; pp. 75–92. [Google Scholar]
  2. Vik, M.; Periyasamy, A.P.; Viková, M. Chromic Materials, Fundamentals, Measurements and Applications; Apple Academic Press Inc.: Waretown, NJ, USA, 2018; pp. 35–155. [Google Scholar]
  3. van der Werff, L.; Kyratzis, I.L.; Robinson, A.; Cranston, R.; Peeters, G.; O’Shea, M.; Nichols, L. Thermochromic composite fibres containing liquid crystals formed via melt extrusion. J. Mater. Sci. 2013, 48, 5005–5011. [Google Scholar] [CrossRef]
  4. He, Y.; Sun, S.; Han, N.; Zhang, X.; Li, W. Thermal energy regulated and thermochromic composite film with temperature-sensitive “breathable” stomata. J. Mater. Sci. 2020, 55, 12921–12939. [Google Scholar] [CrossRef]
  5. MacLaren, D.C.; White, M.A. Design rules for reversible thermochromic mixtures. J. Mater. Sci. 2005, 40, 669–676. [Google Scholar] [CrossRef]
  6. Tatíčková, Z.; Kudláček, J.; Zoubek, M.; Kuchař, J. Behaviour of Thermochromic Coatings under Thermal Exposure. Coatings 2023, 13, 642. [Google Scholar] [CrossRef]
  7. Guo, X.; Daka, S.; Fan, M.; Lin, X.; Sun, W. Reversibly thermochromic wood. J. Mater. Sci. 2023, 58, 2188–2197. [Google Scholar] [CrossRef]
  8. Cheng, Y.; Zhang, X.; Fang, C.; Chen, J.; Wang, Z. Discoloration mechanism, structures and recent applications of thermochromic materials via different methods: A review. J. Mater. Sci. Technol. 2018, 34, 2225–2234. [Google Scholar] [CrossRef]
  9. Ribeiro, L.S.; Pinto, T.; Monteiro, A.; Soares, O.S.G.P.; Pereira, C.; Freire, C.; Pereira, M.F.R. Silica nanoparticles functionalized with a thermochromic dye for textile applications. J. Mater. Sci. 2013, 48, 5085–5092. [Google Scholar] [CrossRef]
  10. Hajzeri, M.; Bašnec, K.; Bele, M.; Gunde, M.K. Influence of developer on structural, optical and thermal properties of a benzofluoran-based thermochromic composite. Dye. Pigment. 2015, 113, 754–762. [Google Scholar] [CrossRef]
  11. Jin, Y.; Bai, Y.; Zhu, Y.; Li, X.; Ge, M. Thermosensitive luminous fiber based on reversible thermochromic crystal violet lactone pigment. Dye. Pigment. 2017, 146, 567–575. [Google Scholar] [CrossRef]
  12. Bašnec, K.; Perše, L.S.; Šumiga, B.; Huskic, M.; Meden, A.; Hladnik, A.; Podgornik, B.B.; Gunde, M.K. Relation between colour- and phase changes of a leuco dye-based thermochromic composite. Sci. Rep. 2018, 8, 5511. [Google Scholar] [CrossRef]
  13. Kim, I.J.; Ramalingam, M.; Son, Y.-A. Investigation of reversible self-thermochromism in microencapsulated fluoran-based materials. Dye. Pigment. 2018, 151, 64–74. [Google Scholar] [CrossRef]
  14. Kulčar, R.; Friškovec, M.; Hauptman, N.; Vesel, A.; Gunde, M.K. Colorimetric properties of reversible thermochromic printing inks. Dye. Pigment. 2010, 86, 271–277. [Google Scholar] [CrossRef]
  15. Aklujkar, P.S.; Kandasubramanian, B. A review of microencapsulated thermochromic coatings for sustainable building applications. J. Coat. Technol. Res. 2021, 18, 19–37. [Google Scholar] [CrossRef]
  16. Chowdhury, M.A.; Butola, B.S.; Joshi, M. Application of thermochromic colorants on textiles: Temperature dependence of colorimetric properties. Color. Technol. 2013, 129, 232–237. [Google Scholar] [CrossRef]
  17. Geng, X.; Li, W.; Wang, Y.; Lu, J.; Wang, J.; Wang, N.; Li, J.; Zhang, X. Reversible thermochromic microencapsulated phase change materials for thermal energy storage application in thermal protective clothing. Appl. Energy 2018, 217, 281–294. [Google Scholar] [CrossRef]
  18. Liu, B.; Rasines Mazo, A.; Gurr, P.A.; Qiao, G.G. Reversible Nontoxic Thermochromic Microcapsules. ACS Appl. Mater. Interfaces 2020, 12, 9782–9789. [Google Scholar] [CrossRef] [PubMed]
  19. Seeboth, A.; Ruhmann, R.; Mühling, O. Thermotropic and Thermochromic Polymer Based Materials for Adaptive Solar Control. Materials 2010, 3, 5143–5168. [Google Scholar] [CrossRef]
  20. Zhang, W.; Ji, X.; Al-Hashimi, M.; Wang, C.; Fang, L. Feasible fabrication and textile application of polymer composites featuring dual optical thermoresponses. Chem. Eng. J. 2021, 419, 129553. [Google Scholar] [CrossRef]
  21. Nguyen, D.K.; Bach, Q.V.; Lee, J.H.; Kim, I.T. Synthesis and Irreversible Thermochromic Sensor Applications of Manganese Violet. Materials 2018, 11, 1693. [Google Scholar] [CrossRef]
  22. Bigger, S.W.; Gribben, M.; Rogers, A. Colour changes in pigmented Nylon-6 systems containing copper iodide. Polym. Degrad. Stab. 2008, 93, 1711–1714. [Google Scholar] [CrossRef]
  23. Pareja-Rivera, C.; Solis-Ibarra, D. Reversible and Irreversible Thermochromism in Copper-Based Halide Perovskites. Adv. Opt. Mater. 2021, 9, 2100633. [Google Scholar] [CrossRef]
  24. Li, Y.; Hui, B.; Li, G.; Li, J. Fabrication of smart wood with reversible thermoresponsive performance. J. Mater. Sci. 2017, 52, 7688–7697. [Google Scholar] [CrossRef]
  25. Kumar, S.; Qadir, A.; Maury, F.; Bahlawane, N. Visible Thermochromism in Vanadium Pentoxide Coatings. ACS Appl. Mater. Interfaces 2017, 9, 21447–21456. [Google Scholar] [CrossRef] [PubMed]
  26. Rougeau, L.; Picq, D.; Rastello, M.; Frantz, Y. New irreversible thermochromic polydiacetylenes. Tetrahedron 2008, 64, 9430–9436. [Google Scholar] [CrossRef]
  27. Garcia-Souto, M.D.P.; Dabnichki, P. Non-invasive and wearable early fever detection system for young children. Measurement 2018, 116, 216–229. [Google Scholar] [CrossRef]
  28. Cabral, I.; Santiago, D.; Steffens, F. Chromic textiles: Colour fastness properties and irreversible colour change behaviour of textiles screen printed with thermochromic, photochromic and hydrochromic colourants. Color. Technol. 2023, 139, 200–208. [Google Scholar] [CrossRef]
  29. Rožić, M.; Šegota, N.; Vukoje, M.; Kulčar, R.; Šegota, S. Description of Thermochromic Offset Prints Morphologies Depending on Printing Substrate. Appl. Sci. 2020, 10, 8095. [Google Scholar] [CrossRef]
  30. Ahn, S.J. Geometric fitting of parametric curves and surfaces. J. Inf. Process. Syst. 2008, 4, 153–158. [Google Scholar] [CrossRef]
  31. Milanov, P.; Koroleova, G.; Mavrevski, R.; Pencheva, N. Curve fitting problem: Torque–velocity relationship with polynomials and Boltzmann sigmoid functions. Acta Bioeng. Biomech. 2018, 20, 169–184. [Google Scholar]
  32. Carter, E.C.; Schanda, J.; Hirschler, R.; Jost, S.; Luo, M.R.; Melgosa, M.; Ohno, Y.; Pointer, M.R.; Rich, D.C.; Viénot, F.; et al. Colorimetry, 4th ed.; CIE 015:2018; CIE: Vienna, Austria, 2018. [Google Scholar] [CrossRef]
  33. Shamey, R. Encyclopedia of Color Science and Technology; Springer: Berlin/Heidelberg, Germany, 2023. [Google Scholar] [CrossRef]
  34. IEC 60050-845:2020; International Electrotechnical Vocabulary (IEV)—Part 845: Lighting. IEC: Geneva, Switzerland, 2020.
  35. Hunt, R.W.G.; Pointer, M.R. Measuring Colour; Wiley: Hoboken, NJ, USA, 2011; pp. 53–57. [Google Scholar]
  36. Abbasi, A.; Vik, M.; Viková, M. Color difference formulae evaluation by method of adjustment. World J. Eng. 2014, 11, 89–94. [Google Scholar] [CrossRef]
  37. Flores, J.C.; Palma-Chilla, L. Entropy behavior for isolated systems containing bounded and unbounded states: Latent heat at the inflection point. J. Phys. Commun. 2020, 4, 035002. [Google Scholar] [CrossRef]
  38. Xu, M.; Pan, X.; Wu, Y. Inflection point as a characteristic of the QCD critical point. Nucl. Phys. A 2014, 927, 69–77. [Google Scholar] [CrossRef]
  39. Dylewsky, D.; Lenton, T.M.; Scheffer, M.; Bury, T.M.; Fletcher, C.G.; Anand, M.; Bauch, C.T. Universal early warning signals of phase transitions in climate systems. J. R. Soc. Interface 2023, 20, 20220562. [Google Scholar] [CrossRef] [PubMed]
  40. Raditoiu, A.; Raditoiu, V.; Nicolae, C.A.; Raduly, M.F.; Amariutei, V.; Wagner, L.E. Optical and structural dynamical behavior of Crystal Violet Lactone—Phenolphthalein binary thermochromic systems. Dye. Pigment. 2016, 134, 69–76. [Google Scholar] [CrossRef]
  41. Gompertz, B. On the nature of the function expressive of the law of human mortality, and on a new mode of determining the value of life contingencies. Philos. Trans. R. Soc. Lond. 1825, 115, 513–583. [Google Scholar]
  42. Richards, F.J. A flexible growth function for empirical use. J. Exp. Bot. 1959, 10, 290–301. [Google Scholar] [CrossRef]
  43. El Aferni, A.; Guettari, M.; Tajouri, T. Mathematical model of Boltzmann’s sigmoidal equation applicable to the spreading of the coronavirus (Covid-19) waves. Environ. Sci. Pollut. Res. Int. 2021, 28, 40400–40408. [Google Scholar] [CrossRef]
  44. Fernandes, T.J.; Pereira, A.A.; Muniz, J.A. Double sigmoidal models describing the growth of coffee berries. Ciência Rural. 2017, 47, e20160646. [Google Scholar] [CrossRef]
  45. Ricketts, J.H.; Head, G.A. A five-parameter logistic equation for investigating asymmetry of curvature in baroreflex studies. Am. J. Physiol. -Regul. Integr. Comp. Physiol. 1999, 277, R441–R454. [Google Scholar] [CrossRef]
  46. Little, T.D. (Ed.) The Oxford Handbook of Quantitative Methods in Psychology; Statistical Analysis; Oxford Library of Psychology: Oxford, UK, 2013; Volume 2. [Google Scholar]
  47. Moore, D.S.; McCabe, G.P.; Graig, B.A. Introduction to the Practice of Statistics; Macmillan Learning: New York, NY, USA, 2021. [Google Scholar]
Figure 1. Schematic representation of thermochromic change. Part (A) of the figure is an illustration of irreversible thermochromic change and part (B) demonstrates reversible thermochromic change. The reversible thermochromic change in this case is an illustration of composite thermochromic inks, where d stands for the developer and the circle stands for the Leuco dye itself.
Figure 1. Schematic representation of thermochromic change. Part (A) of the figure is an illustration of irreversible thermochromic change and part (B) demonstrates reversible thermochromic change. The reversible thermochromic change in this case is an illustration of composite thermochromic inks, where d stands for the developer and the circle stands for the Leuco dye itself.
Materials 16 07478 g001
Figure 2. Representation of the positions of individual thermochromic samples (blue and orange) at different temperatures in the CIE Y, x, y color space (A) and in the CIELAB color space (B). The concentration of both samples was 150 g·kg−1. The vertical axes represent a change in brightness (Y) and lightness (L*), respectively, of both colors of thermochromic samples on CIE diagram x, y (A) and in CIELAB color space (B). Concentration of both samples was 150 g·kg−1.
Figure 2. Representation of the positions of individual thermochromic samples (blue and orange) at different temperatures in the CIE Y, x, y color space (A) and in the CIELAB color space (B). The concentration of both samples was 150 g·kg−1. The vertical axes represent a change in brightness (Y) and lightness (L*), respectively, of both colors of thermochromic samples on CIE diagram x, y (A) and in CIELAB color space (B). Concentration of both samples was 150 g·kg−1.
Materials 16 07478 g002
Figure 3. (K/S) transition as a function of temperature for inks (A) Fast Blue; (B) Gold Orange.
Figure 3. (K/S) transition as a function of temperature for inks (A) Fast Blue; (B) Gold Orange.
Materials 16 07478 g003
Figure 4. Dependence of temperature distance ΔT (A) and slope coefficient α (B) on the concentration of Fast Blue and Gold Orange thermochromic inks.
Figure 4. Dependence of temperature distance ΔT (A) and slope coefficient α (B) on the concentration of Fast Blue and Gold Orange thermochromic inks.
Materials 16 07478 g004
Figure 5. Measured values of total color difference ΔE* and lightness difference ΔL* in relation to sample temperature for 300 g·kg−1 Fast Blue (A) and 300 g·kg−1 Gold Orange (B) inks.
Figure 5. Measured values of total color difference ΔE* and lightness difference ΔL* in relation to sample temperature for 300 g·kg−1 Fast Blue (A) and 300 g·kg−1 Gold Orange (B) inks.
Materials 16 07478 g005
Figure 6. Measured values of chroma difference ΔC* and hue difference ΔH* in relation to sample temperature for 300 g·kg−1 Fast Blue (A) and 300 g·kg−1 Gold Orange (B) inks.
Figure 6. Measured values of chroma difference ΔC* and hue difference ΔH* in relation to sample temperature for 300 g·kg−1 Fast Blue (A) and 300 g·kg−1 Gold Orange (B) inks.
Materials 16 07478 g006
Figure 7. Fit of temperature dependence of total color difference ΔE* for 300 g·kg−1 Fast Blue (A) and 300 g·kg−1 Gold Orange (B) inks.
Figure 7. Fit of temperature dependence of total color difference ΔE* for 300 g·kg−1 Fast Blue (A) and 300 g·kg−1 Gold Orange (B) inks.
Materials 16 07478 g007
Figure 8. Fit of temperature dependence of (K/S) and total color difference ΔE* for 300 g·kg−1 Fast Blue (A) and 300 g·kg−1 Gold Orange (B) inks.
Figure 8. Fit of temperature dependence of (K/S) and total color difference ΔE* for 300 g·kg−1 Fast Blue (A) and 300 g·kg−1 Gold Orange (B) inks.
Materials 16 07478 g008
Table 1. Estimated parameters of Equations (14) and (17) for Fast Blue thermochromic ink.
Table 1. Estimated parameters of Equations (14) and (17) for Fast Blue thermochromic ink.
Concentration
[g·kg−1]
(K/S)b
[—]
(K/S)f
[—]
TT
[°C]
α
[—]
TR10
[°C]
TR90
[°C]
ΔT
[°C]
100.0400.48436.5−1.0534.238.84.6
500.0242.32936.6−0.9934.438.74.4
1000.0194.05836.6−0.9934.538.84.3
1500.0456.96736.8−0.9034.838.73.9
3000.07210.33036.7−0.8534.838.63.8
Table 2. Estimated parameters of Equations (14) and (17) for Gold Orange thermochromic ink.
Table 2. Estimated parameters of Equations (14) and (17) for Gold Orange thermochromic ink.
Concentration
[g·kg−1]
(K/S)b
[—]
(K/S)f
[—]
TT
[°C]
α
[—]
TR10
[°C]
TR90
[°C]
ΔT
[°C]
100.0810.63638.8−1.2636.041.65.5
500.1253.10838.4−0.9336.440.54.1
1000.1755.72138.3−0.7736.640.03.4
1500.0936.27439.0−0.6537.640.42.8
3000.2548.74639.2−0.5837.940.42.6
Table 3. Estimated TT and evaluation criteria of tested models for 300 g·kg−1 Fast Blue thermochromic ink and total color difference ΔE*.
Table 3. Estimated TT and evaluation criteria of tested models for 300 g·kg−1 Fast Blue thermochromic ink and total color difference ΔE*.
Tested EquationTT
[°C]
F
(DFn, DFd)
AIC
[—]
R2
[—]
AR2
[—]
RMSE
[—]
SBE (14)38.029.250.99560.99451.562
Gompertz (22)37.839.860.99150.98942.177
Richards (23)37.61.82 (1, 11) H032.130.99620.99491.447
DSBE1 * (24)37.7 (40.6)9.09 (3, 9) H127.520.99890.99820.778
DSBE2 * (25)37.7 (40.5)13.86 (2, 10) H120.010.99880.99830.805
BARO5 (26)37.97.19 (1,11) H126.540.99730.99641.215
* Two critical temperatures cA and cB are calculated in the DSBE equations, therefore, the second temperature is given in brackets.
Table 4. Estimated TT and evaluation criteria of tested models for 300 g·kg−1 Gold Orange thermochromic ink and total color difference ΔE*.
Table 4. Estimated TT and evaluation criteria of tested models for 300 g·kg−1 Gold Orange thermochromic ink and total color difference ΔE*.
Tested EquationTT
[°C]
F
(DFn, DFd)
AIC
[—]
R2
[—]
AR2
[—]
RMSE
[—]
SBE (14)40.0029.120.99660.99581.556
Gompertz (22)39.619.080.99820.99771.137
Richards (23)39.69.58 (1,11) H124.430.99820.99751.138
DSBE1 * (24)39.5 (40.4)27.62 (3, 9) H112.520.99970.99940.487
DSBE2 * (25)38.4 (38.9)4.60 (2, 10) H130.670.99820.99741.123
BARO5 (26)39.941.42 (1, 11) H19.4660.99930.99900.713
* Two critical temperatures cA and cB are calculated in the DSBE equations, therefore, second for the Gold Orange thermochromic ink, the DSBE1 was evaluated as the best performing model in terms of R2, AR2 and RSME, including the F-test against the DSBE2 and BARO5 models. At the same time, we can see that in both cases, the transition temperature estimates for thermochromic inks are lower than when using SBE.
Table 5. Estimated TT and evaluation criteria of tested models for 300 g·kg−1 Fast Blue thermochromic ink and Kubelka–Munk function (K/S).
Table 5. Estimated TT and evaluation criteria of tested models for 300 g·kg−1 Fast Blue thermochromic ink and Kubelka–Munk function (K/S).
Tested EquationTT
[°C]
F
(DFn, DFd)
AIC
[—]
R2
[—]
AR2
[—]
RMSE
[—]
SBE (14)36.7−23.090.99570.99470.305
Gompertz (22)36.8−11.160.99100.98880.442
Richards (23)37.27.07 (1, 11) H1−25.690.99740.99650.238
DSBE1 * (24)37.3 (34.9)7.25 (3, 9) H1−22.170.99880.99790.165
DSBE2 * (25)------
BARO5 (26)36.86.41 (1, 11) H1−25.100.99730.99630.242
* Two critical temperatures cA and cB are calculated in the DSBE equations, therefore the second temperature is given in brackets.
Table 6. Estimated TT and evaluation criteria of tested models for 300 g·kg−1 Gold Orange thermochromic ink and Kubelka–Munk function (K/S).
Table 6. Estimated TT and evaluation criteria of tested models for 300 g·kg−1 Gold Orange thermochromic ink and Kubelka–Munk function (K/S).
Tested EquationTT
[°C]
F
(DFn, DFd)
AIC
[—]
R2
[—]
AR2
[—]
RMSE
[—]
SBE (14)39.2−28.710.99580.99480.255
Gompertz (22)39.2−17.850.99180.98970.359
Richards (23)39.54.273 (1, 11) H0−28.630.99700.99590.217
DSBE1 * (24)39.3 (38.35)8.117 (3, 9) H1−29.100.99890.99810.133
DSBE2 * (25)------
BARO5 (26)39.210.60 (1, 11) H1−34.180.99790.99710.182
* The parameters of the DSBE1 model are not shown because their calculation was unstable.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Viková, M.; Vik, M. Transition Temperature of Color Change in Thermochromic Systems and Its Description Using Sigmoidal Models. Materials 2023, 16, 7478. https://doi.org/10.3390/ma16237478

AMA Style

Viková M, Vik M. Transition Temperature of Color Change in Thermochromic Systems and Its Description Using Sigmoidal Models. Materials. 2023; 16(23):7478. https://doi.org/10.3390/ma16237478

Chicago/Turabian Style

Viková, Martina, and Michal Vik. 2023. "Transition Temperature of Color Change in Thermochromic Systems and Its Description Using Sigmoidal Models" Materials 16, no. 23: 7478. https://doi.org/10.3390/ma16237478

APA Style

Viková, M., & Vik, M. (2023). Transition Temperature of Color Change in Thermochromic Systems and Its Description Using Sigmoidal Models. Materials, 16(23), 7478. https://doi.org/10.3390/ma16237478

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop