3.1. Evolution of Neuroaesthetic Variables Throughout the Renaissance
In this study, we were interested in the drift of values of neuroaesthetic variables in relation to the passing of time and the evolution of art. Because specialized brain mechanisms constrain these variables, one may expect that they remain relatively constant over time. In contrast, a recent study demonstrated that a certain degree of flexibility appears to exist with respect to neuroaesthetic variables [
4]. Therefore, they could potentially evolve across different periods of art. We thus asked whether changes across art periods can occur in the absence of evolution of neuroaesthetic variables. To answer this question, we first performed computational measurements of asymmetry, imbalance, and complexity (normalized entropy—[
3]). In particular, we probed the changes that happened to these variables in Italy and in the rest of Europe. Our study involved a time span of close to 200 years of Renaissance. In
Figure 1, we see the results divided to the periods of Early Renaissance, High Renaissance, and Mannerism.
Our analysis revealed no significant changes in symmetry across the Early Renaissance, High Renaissance, and Mannerism (
Figure 1a). In contrast, imbalance rose significantly between the Early and High Renaissance (
Figure 1b, two-way ANOVA and post-hoc two-sided
t-test, 298 d.f., t = 3.19,
p < 0.002). The degree of imbalance grew by almost 30% in the span of 80 years. Complexity (i.e., normalized entropy) also evolved over time. We found falls in Complexities of Order 1 and Order 2 between the Early and High Renaissance (
Figure 1c,d). These falls were significant for both Order 1 (432 d.f., t = 6.90,
p < 2 × 10
−11) and Order 2 (449 d.f., t = 7.26,
p < 2 × 10
−12). The falls reduced complexities of both orders by about 10%. Interestingly, however, no changes occurred in imbalance or complexity from High Renaissance to Mannerism. Hence, all changes in neuroaesthetic variables took place during the Early Renaissance. Finally, although we detected temporal changes in these variables throughout the Renaissance, we found no significant differences between Italy and the rest of Europe.
In conclusion, although symmetry was constant throughout the Renaissance, balance and complexities fell during the Early Renaissance.
3.2. Abrupt Transitions at the End of the 15th Century
To quantify these results further and to compare top artists with the other painters in our dataset, we produced scatter plots of the data. An example of the analysis of these plots appears in
Figure 2 for Complexity of Order 2 (normalized spatial entropy).
The basic scatter plot for Complexity of Order 2 appears in
Figure 2a, in which each point represents an individual painting. The abscissas correspond to the years of painting completion and the ordinates are the measured complexities. The data show great variability in each moment of portrait evolution. To quantify the variability, we calculated the ratio between the standard deviation and the mean of the Complexity of Order 2. Overall, this variability ratio was 15%. Despite the variability, the results in
Figure 2a confirmed and extended the temporal trends in
Figure 1. One observes that the Complexity of Order 2 falls during the Renaissance. The Kendall’s
correlation coefficient was statistically significantly negative for Complexity of Order 2 (Kendall’s
= −0.211,
p < 3 × 10
−11).
We attempted to characterize the fall of Complexity of Order 2 throughout the Renaissance by fitting four models (Equations (1)–(4);
Figure 2b).
Figure 1d had suggested that this fall was nonlinear and thus, we attempted exponential and error-function fits (Equations (3)–(5)). The latter seemed especially relevant, because no fall was apparent from the High Renaissance to the Mannerism. For completeness, we also attempted constant and linear fits. All the fits were statistically robust, by first extracting median complexities in small sections of the data (black dashed line in
Figure 2b—obtained from the scatter plot with medians from 25 paintings).
The median Complexity-of-order-2 curve appeared to exhibit an abrupt fall around 1490. Not surprisingly, therefore, the Error-function model (Equation (4)) provided the best fit to the data. For example, the sums of squared errors for the optimal fits were 0.036, 0.014, 0.013, and 0.0068 for the Constant, Linear, Exponential, and Error-function models respectively. However, that the fit was better for the Error-function model was perhaps not surprising, because it had more parameters and could subsume some of the other models. Hence, we tested the quality of the fits with a regression test for arbitrary fits (
Section 2.4; [
21,
22]). This test considered the number of parameters of the models. The test first plotted model predictions against the data and then analyzed the statistics of the resulting linear regression. The predictions were of positive correlation, with an intercept of 0 and a slope of 1. For all models, except the Constant one, we could not reject the null hypothesis that the correlation was positive. But the correlation was highest for the Error-function model (
, and
for the Constant, Linear, Exponential, and Error-function models respectively). Furthermore, we could reject that the intercept was zero for the Constant, Linear, and Exponential models (
p < 0.0001,
p < 0.005, and
p < 0.008 respectively). In contrast, we could not reject this null hypothesis for the Error-function fit. Similarly, although we could not reject that the slope was 1 for the Error-function model, we could reject this null hypothesis for the Constant, Linear, and Exponential models (
p < 0.0001,
p < 0.005, and
p < 0.008 respectively). Consequently, the Error-function model provided a superior fit than did the others. This superiority was true for all the fits in this article for data exhibiting trends. Moreover, we could not reject the Error-function model for any of these data.
The excellent error-function fit reinforced the conclusion of an abrupt fall of Complexity of Order 2 around 1490. The optimal transient year parameter (
in Equation (4)) was 1493. In addition, the optimal transition was indeed abrupt as shown by the red line
Figure 2b. However, the transition was not as abrupt as suggested by the red line. This line was obtained by fitting the curve of medians from 25 paintings, corresponding to a span of 12 years around 1493, namely [1488–1500]. Therefore, all that we could say was that 12 years was the upper bound for the duration of the transition. We call this time window ([1488–1500] for this example) the upper-bound transition interval.
Such an abrupt transition was surprising, because it was not immediately apparent in the scatter plot (
Figure 2a). We thus wished to obtain model-independent evidence for such a transition. This evidence is what
Figure 2c,d show. In
Figure 2c, we show the Kendall’s
correlation coefficients for three non-overlapping section of the data. The middle section has 50 paintings around 1493. The other sections have all the paintings before and after the middle section. As the figure shows, although the middle section is far smaller than are the others, only it has a statistically significantly negative Kendall’s
. The Kendall’s
for the middle section is −0.275, with
p < 0.004 that the correlation coefficient is 0. Comparison of this Kendall’s
with that obtained for the entire data (−0.211) suggests that most of the fall of Complexity of Order 2 happens during the period encompassed by the middle section [1478–1502]. This result is compatible with the upper-bound transition interval estimated above.
Further confirmation of the conclusion of abrupt transition appears in
Figure 2d. This figure plots the Kendall’s
’s and their standard errors for non-overlapping, consecutive sections of the data comprising 55 paintings. Only one point in the plot is statistically significantly negative, namely, the one centered on 1492.
In conclusion, our data indicate an abrupt fall in Complexity of Order 2 in the last decade of the 15th century. Similar analyses have shown abrupt transitions in most other variables studied in this paper, except as indicated.
3.3. Dynamics of Neuroaesthetic Variables and the Top Painters
We extended the scatter-plot analysis of
Section 3.2 to Asymmetry, Imbalance, and Complexity of Order 1. In particular, we superimposed on the scatter plots temporal trend lines to help quantify the time courses of the drift of values of neuroaesthetic variables (Equations (4) and (5)). Finally, we colored the points of the top artists (
Section 2.2) to compare them with peers from their periods. The results appear in
Figure 3.
As for Complexity of Order 2 (
Figure 2),
Figure 3 shows great variability in each moment of portrait evolution. This variability holds for the results of both the whole group of painters and the work of the leading masters. The ratio between the standard deviation and the mean is 28%, 73%, and 13% for asymmetry, imbalance, and Complexity of Order 1 respectively. Consequently, the variability is lowest for complexities (see also
Section 3.2) and highest for balance.
Also, like
Figure 2, despite the variability, the results in
Figure 3 confirmed and extended the temporal trends in
Figure 1. The Kendall’s
correlation coefficient was not significantly different from zero for asymmetry (see its robust constant regression in
Figure 3a). However, the coefficient was statistically significantly positive for Imbalance (Kendall’s
= 0.0738,
p < 0.03). In contrast, the coefficient was statistically significantly negative for Complexity of Order 1 (Kendall’s
= −0.204,
p < 2 × 10
−10), as it was for Order 2 (
Section 3.2). The Error-function regression lines (Equation (4)) attempted to capture these positive and negative tendencies in the evolution of the neuroaesthetic variables. The line rose for imbalance (
Figure 3b; obtained from the scatter plot with medians from 30 paintings). But the line fell for complexities (
Figure 3c,d; obtained from the scatter plot with medians from 25 paintings). The best-fit transition times for Imbalance and Complexities of Order 1 were 1467 and 1502 respectively. The corresponding upper-bound transition intervals (see discussion of
Figure 2b) were [1451–1473] for imbalance and [1499–1505] for Complexity of Order 1. Therefore, the transition was about 25 years earlier for imbalance than Complexity of Order 2 (see also
Section 3.2). The transition was also slower for imbalance. Moreover, the transition may have been about 10 years later for Complexity of Order 1 than of Order 2.
The top painters did not generally appear to behave differently from their peers in terms of neuroaesthetic variables. Statistical comparisons of the values of their neuroaesthetic variables show relatively equal distribution above and below the trend lines. Titian was the only exception for Complexity of Order 1 (t = 3.48, 17 d.f., and p < 0.003). In turn, two painters were exceptions for Complexity of Order 2. They were van Eyck (t = 3.42, 5 d.f., and p < 0.02) and Titian (t = 3.11, 17 d.f., and p < 0.007). These painters produced portraits that were less complex than were those of the peers, as evaluated by the Error-function model.
In sum, balance and complexity declined abruptly towards the end of the 15th century, but this fall was not generally due to the top painters of those times. Together with the results of
Section 3.1 and
Section 3.2, we thus established that neuroaesthetic variables were not constant throughout the Renaissance.
3.4. Evolution of the Thickness of Balance Lines
To understand the decline of balance over time, we must start from the definition of imbalance. We defined it relative to the midline of the canvas. Elsewhere, we also considered the position of balance, i.e., the place for which the integrals of intensities to the right and left of it were equal [
4]. Thus, the decline of balance in the Renaissance meant that the distance between the midline and the position of balance tended to increase over time (see examples in
Figure 4a,b). However, this decline did not imply a rise in the sloppiness of balancing different parts of the portrait. Painters could continue to balance portraits delicately but simply do it in a position of balance away from the midline. In an earlier publication, we reported that painters in the Early Renaissance not only balanced their portraits, but also did so at every row of the canvas [
4]. We thus decided to test if this delicate form of balance was also diminished as the Renaissance progressed. To do so, we measured the positions of balance at every height of the painting. All these points together formed the Balance Line (
Figure 4c,d). We then measured the thickness of this line as a fraction of the horizontal size of the canvas. The results appear in
Figure 4e,f.
The results show that the thickness of the balance line also rises in the Renaissance. We can appreciate an example by comparing a portrait by Domenico Veneziano in the Early Renaissance with one by Giovanni Battista Moroni during Mannerism (
Figure 4c,d, respectively). In Veneziano’s portrait, the balance line shows that the distributions of intensities on the two sides of the midline of the canvas are similar. The balance line is close to midline at every height analyzed. In contrast, in Moroni’s portrait, the balance line has more variation across vertical positions. Hence, the balance line in Moroni’s portrait has more thickness (0.165) than in Veneziano’s (0.018). Thus, Moroni was “sloppier” in balancing different parts of his portrait than was Veneziano.
This difference held when we analyzed the thicknesses of balance lines throughout the Renaissance. In the Early Renaissance, the thickness of the balance line was significantly lower in Italy than in the rest of Europe (two-way ANOVA and post-hoc two-sided post-hoc
t-test, 149 d.f.,
t = 4.05,
p < 9 × 10
−5—
Figure 4e).
Afterwards the thickness of the balance line grew in Italy from the Early to High Renaissance (211 d.f.,
t = 7.00,
p < 4 × 10
−11—
Figure 4e), catching up with the values in the rest of Europe. In contrast, the thickness of the balance line was statistically constant in the rest of Europe throughout the Renaissance. Therefore, portraits in the rest of Europe were more prescient of future trends of balance than were Italian ones.
The scatter plots showed that the thicknesses of the balance lines (
Figure 4f) followed a trend like that of imbalance (
Figure 3d). The data in
Figure 4f show great variability of thicknesses in each moment of portrait evolution. The ratio between the standard deviation and the mean of the thicknesses of balance lines was 50%. Despite the variability of thicknesses, the results in
Figure 4f confirmed and extended the temporal trends in
Figure 4e. Portraits tended to be carefully balanced in the Early Renaissance, but exhibit sloppier balances in the High Renaissance and Mannerism. Accordingly, the Kendall’s
correlation coefficient was statistically significantly positive for the thicknesses of balance lines (Kendall’s
= 0.190,
p < 1.21 × 10
−9). The best-fit transition time was 1467, confirming that most change happened in the Early Renaissance. However, the change was much slower for the thickness of balance line than for other aesthetic variables. Consequently, its change was not abrupt. Finally, as for imbalance (
Figure 3b), top painters did not generally produce portraits with thicker balance lines than those of peers (
Figure 4f). In conclusion, as the Early Renaissance progressed, painters tended to become “sloppier” in balancing different parts of the portrait.
3.5. Evolution of Spatial Complexity
How are we to understand the decline of complexity over time? The Complexities of Order 1 and 2 in
Figure 1 and
Figure 3 have different types of interpretation [
4]. Complexity of Order 1 measures the normalized entropy in the distribution of intensities in the image. In turn, Complexity of Order 2 begins from Complexity of Order 1 and then discounts the reduction of entropy due to spatial organization. Consequently, if we want to isolate the loss of complexity due to spatial organization alone, we must calculate Complexity of Order 1 minus Complexity of Order 2. We call this quantity the Spatial Simplicity [
4]. The temporal drift of the Spatial Simplicity throughout the Renaissance appears in
Figure 5.
Portrait paintings tended to become spatially simpler as the Renaissance progressed. Thus, in High Renaissance and Mannerist periods, spatial complexity was lower than in the Early Renaissance (
Figure 5a;
t = 4.38, 447 d.f.,
p < 0.00002). However, as for
Figure 1, although we detected temporal changes in Spatial Simplicity throughout the Renaissance, we found no significant differences between Italy and the rest of Europe. The scatter plots confirmed the rise of spatial simplicity (
Figure 5b). Accordingly, the Kendall’s
correlation coefficient was statistically significantly positive for spatial simplicity (
τ = 0.0994,
p < 0.002). The Error-function regression lines (Equation (4)) rose abruptly for Spatial Simplicity (
Figure 5b; obtained from the scatter plot with medians from 35 paintings). The best-fit transition time was 1486, with the upper-bound transition interval lasting 17 years, namely, [1477–1494]. Finally, four of the seven top painters produced portraits with different spatial-simplicity distributions than those of their peers (
Figure 5b). Botticelli (
t = 2.55, 8 d.f.,
p < 0.04), van Eyck (
t = 3.19, 4 d.f.,
p < 0.04), and Raphael (
t = 3.76, 18 d.f.,
p < 0.002) exhibited more spatial simplicity than did their peers. In contrast, Caravaggio exhibited less (
t = 3.60, 7 d.f.,
p < 0.009).
Hence, spatial complexity (the component of entropy due to spatial organization) fell abruptly towards the end of the Early Renaissance. This fall mimicked the decline of the complexity due to the distribution of intensities (Complexity of Order 1—
Figure 3c).
3.6. Chiaroscuro and the Fall of Complexity
The decline of complexity over time (
Figure 1c,d,
Figure 2,
Figure 3c,d and
Figure 5) was surprising to us. We had expected complexity to increase as paintings became more realistic in the Renaissance. Ideas that evolved throughout the Renaissance, such as naturalism, anatomical studies, linear perspective, and aerial perspective should perhaps have made portraits more complex. Therefore, we wondered why complexity fell. We hypothesized that portrait paintings got simpler with the invention of
chiaroscuro. It introduced large dominant regions with fewer colors and homogeneous intensities. We can appreciate an example of such regions by comparing a portrait by Andrea Mantegna in the Early Renaissance with one by Caravaggio during Mannerism (
Figure 6a,b, respectively). Mantegna’s portrait has no dominant regions in terms of blacks and white, and thus has no or very little
chiaroscuro. In contrast, Caravaggio’s portrait is a good example of
chiaroscuro, with some bright regions contrasting against a large, dark background. Hence, the indices of
chiaroscuro extremes (Equation (8)) and
chiaroscuro intensities (Equation (11)) in Mantegna’s portrait (0.33 and 0.094 respectively) were lower than those in Caravaggio’s (0.75 and 0.83 respectively). Thus, our measurements confirm the common knowledge that Caravaggio used
chiaroscuro more than did Mantegna (and most other painters—[
8]). A quantitative study of these
chiaroscuro indices across the Renaissance appears in
Figure 6c–f.
Figure 6c revealed that the index of
chiaroscuro extremes (Equation (8)) rose between the Early Renaissance and the High Renaissance periods (two-way ANOVA and post-hoc two-sided
t-test, 303 d.f.,
t = 8.28,
p < 4 × 10
−15). However, no such rise occurred from the High Renaissance to Mannerism.
Figure 6d also showed that this index grew abruptly towards the end of the Early Renaissance (Kendall’s
= 0.231,
p < 3 × 10
−13; best-fit transition time = 1487, upper-bound transition interval lasting 15 years, namely, [1479, 1494]). These findings were replicated for the index of
chiaroscuro intensities (Equation (11)) in
Figure 6e,f (313 d.f.,
t = 5.92,
p < 9 × 10
−9; Kendall’s
= 0.194,
p < 8 × 10
−10; best-fit transition time = 1481, transition interval lasting 14 years, namely, [1476, 1490]). Another finding in Early Renaissance was the statistical similarity of Italy with the rest of Europe in terms of
chiaroscuro tendencies. Finally, we again detected little difference from top painters and their peers (
Figure 6d,f). The only exceptions were van Eyck for the index of
chiaroscuro extremes (
t = 3.75, 5 d.f.,
p < 0.02), and Caravaggio for the index of
chiaroscuro intensities (
t = 6.67, 7 d.f.,
p < 3 × 10
−4). Both van Eyck and Caravaggio exhibited more chiaroscuro than did their peers.
That the degree of
chiaroscuro usage went up in the Renaissance was consistent with our hypothesis for the decline of complexity. However, we still had to demonstrate that more
chiaroscuro in a portrait tended to lead to less complexity. In
Figure 7, we classified portrait paintings in three compositional concepts that could affect complexity: linear perspective, aerial perspective, and
chiaroscuro (
Section 2.3). We also included a class for those portraits that do not belong to any of these three categories. Finally, we quantified Complexity of Order 1, spatial simplicity, and the index of
chiaroscuro extremes in these four categories.
In
Figure 7a, we observe that the prevalence of
chiaroscuro portraits increases as time progresses in the Renaissance. This increase occurs specially from the Early to the High Renaissance (Fisher’s exact test, odds ratio = 0.30,
p < 5 × 10
−4). Although we categorized these portraits by hand (
Section 2.3),
Figure 7b supported the idea that our
chiaroscuro category was correct. The index of
chiaroscuro extremes was higher for this category than was for the others (one-way ANOVA followed by post-hoc one-sided
t-tests; linear perspective, 114 d.f.,
t = 6.37,
p < 3 × 10
−9; aerial perspective, 151 d.f.,
t = 6.52,
p < 6 × 10
−10; None, 381 d.f.,
t = 8.28,
p < 2 × 10
−15). Consequently, the use of
chiaroscuro techniques increased over time. In contrast, the prevalence of portraits with the other tested compositional categories, namely, linear and aerial perspective, was statistically constant throughout the Renaissance. Hence, of the compositional elements studied,
chiaroscuro is the only candidate available to explain the fall of complexity over time. Is
chiaroscuro contributing to the simplification of portraits? The answer to this question appears in
Figure 7c,d. The former shows that Complexity of Order 1 is significantly lower in portrait paintings with
chiaroscuro than in portraits in the other categories (linear perspective, 106 d.f.,
t = 3.83,
p < 2 × 10
−4; aerial perspective, 143 d.f.,
t = 6.09,
p < 5 × 10
−9; None, 378 d.f.,
t = 5.40,
p < 6 × 10
−8). Consequently,
chiaroscuro reduces the complexity of the distribution of intensities. Furthermore, in
Figure 7d, we see that spatial simplicity is significantly lower in portraits with
chiaroscuro than is in portraits of the other categories (linear perspective, 106 d.f.,
t = 3.30,
p < 2 × 10
−3; aerial perspective, 142 d.f.,
t = 2.57,
p < 6 × 10
−3; None, 372 d.f.,
t = 3.37,
p < 5 × 10
−4). Therefore,
chiaroscuro tends to reduce the spatial complexity of portraits. We conclude that the reduction in complexity in the Renaissance may be due to the rise of
chiaroscuro.