Authors: Felix Schoeller, Leonardo Christov-Moore, Caitlin Lynch, Thomas Diot, Nicco Reggente
Categories: Social and Political Sciences, absorption, aesthetic chills, individual differences, machine learning, positive affect exposure, AcademicSubjects/MED00010, AcademicSubjects/SCI00010, AcademicSubjects/SOC00010, PNAS_Nexus/psych-soc
Source: PNAS Nexus
Why does the same experience elicit strong emotional responses in some individuals while leaving others largely indifferent? Is the variance influenced by who people are (personality traits), how they feel (emotional state), where they come from (demographics), or a unique combination of these? In this 2,900+ participants study, we disentangle the factors that underlie individual variations in the universal experience of aesthetic chills, the feeling of cold and shivers down the spine during peak experiences. Here, we unravel the interplay of psychological and sociocultural dynamics influencing self-reported chills reactions. A novel technique harnessing mass data mining of social media platforms curates the first large database of ecologically sourced chills-evoking stimuli. A combination of machine learning techniques (LASSO and SVM) and multilevel modeling analysis elucidates the interacting roles of demographics, traits, and states factors in the experience of aesthetic chills. These findings highlight a tractable set of features predicting the occurrence and intensity of chills—age, sex, pre-exposure arousal, predisposition to Kama Muta (KAMF), and absorption (modified tellegen absorption scale [MODTAS]), with 73.5% accuracy in predicting the occurrence of chills and accounting for 48% of the variance in chills intensity. While traditional methods typically suffer from a lack of control over the stimuli and their effects, this approach allows for the assignment of stimuli tailored to individual biopsychosocial profiles, thereby, increasing experimental control and decreasing unexplained variability. Further, they elucidate how hidden sociocultural factors, psychological traits, and contextual states shape seemingly “subjective” phenomena.
Keywords: aesthetic chills, machine learning, absorption, positive affect exposure, individual differences
We have all experienced it—sharing a joke, a movie scene, or an inspirational speech with a friend, expecting it to move them as profoundly as it moved us, only to be met with a lukewarm reaction. “Beauty is in the eye of the beholder,” they say. But what accounts for such vastly different reactions to the same stimulus? Are some innately more prone to being moved than others? How do factors like personality and upbringing shape our susceptibility? Though methodologically daunting, disentangling the complex factors underlying these individual differences is key to harnessing the therapeutic promise of peak emotional experiences (1–5). One such universal peak experience is aesthetic chills (henceforth “chills”)—an intense psychophysiological response to specific rewarding or aversive stimuli associated with shivers, goosebumps, and a feeling of cold down the neck and/or spine (6–12). Here chills serve as a somatic marker for peak emotional experiences (13–19)—i.e. bodily sensations that act as emotional, conditioned responses representing complex reinforcer values related to prior experiences (7, 20). Chills have been a research topic of growing interest in the past decade (7, 21), notably for their effects on positive affect (22, 23), reward learning (6, 13, 22, 24), memory and attentional processes (25, 26), prosocial tendencies (27–29), and as a nonpharmacological substitute for dopaminergic-related illnesses, and mood disorders (4, 7, 22). Substantial individual differences exist in both propensity for and intensity of chills (24, 30–40), which hinders progress in understanding and utilizing chills as therapeutic tools for mental health conditions and enhancing positive emotional experiences broadly.
Here, we leverage a large cohort (N = 2,937) and a robust, multimethod approach to disentangle the influences of demographics, traits, and stimulus on chills likelihood and intensity. Participants completed validated measures assessing demographics, personality, and emotional state before exposure to musical and audiovisual stimuli (see Fig. 1), which were empirically confirmed across large samples to elicit chills beforehand in the population of interest (41, 42). Advances in data science have revolutionized the scientific study of subjective phenomena by permitting the analysis of massive datasets of publicly available naturalistic reactions to stimuli provided by social media (43, 44). Here, we used a novel technique leveraging online social platforms and crowdsourcing to curate and empirically validate a stimulus set for eliciting aesthetic chills. Through developing natural language processing algorithms that data mine YouTube and Reddit comments for somatic markers of chills and goosebumps, we extracted 40 emotion-evoking audiovisual clips collectively identified by users to produce shivers (41, 42).
Fig. 1. The interplay of demographics, traits, and states predict the occurrence of aesthetic chills in response to chills stimuli. Each scale represents distinct layers of influence on the individual's chills response. Demographics encompass age, education, gender, and political orientation. Traits include personality aspects like extraversion, conscientiousness, and absorption. States denote emotional factors like valence and arousal.
To disentangle the multifactorial influences shaping self-reported chills, we followed a multilevel statistical approach. First, we used univariate analyses to establish baseline demographic and contextual differences (a.k.a., states) between those who report chills vs. those who do not. Subsequently, we used multivariate regression to model trait predictors of chills reactions. We then performed latent class analysis (LCA) to reveal hidden “cultures” reflecting sociocultural response patterns after accounting for stimuli, states, and trait effects (Fig. 1 and Table 1). This approach is effective for characterizing the dataset and identifying clusters of subjects who are more prone to experience chills or more intense chills. However, it does not have predictive value on its own. For this purpose, we leveraged machine learning techniques with the goal to move toward providing researchers interested in chills with an instrument that would allow them to administer the right stimuli at the right time to the right person, and thereby allow operationalized control of the emotion in experimental settings.
Accordingly, in order to assess the ecological validity of using individual characteristics to predict the occurrence of chills in response to audiovisual stimuli, we leveraged machine learning within a cross-validation framework to predict both the chills occurrence (binary; support vector machine [SVM] approach) and chills intensity (continuous; least absolute shrinkage, and selection operator [LASSO] approach). These approaches not only yield an assessment of predictive power but also identify the most parsimonious combination of demographic, psychological, and contextual factors that constitute an informative feature set for predicting response in unseen individuals.
The true value of models related to human emotions lies not only in their ability to describe but also in their potential to reduce uncertainty about future emotional responses. If a person's characteristics, current state, and demographic information indeed influence their emotional reactions to stimuli, then having this information before exposure to the stimuli should allow for accurate predictions about their responses. Successfully implementing this approach could enable personalized content curation aimed at maximizing the likelihood of eliciting chills, which holds significant implications for both research and clinical applications. We hypothesized that individuals who score higher on measures of openness to experience, absorption, and dispositional positive emotions would report more frequent and intense experiences of chills in response to validated audiovisual stimuli compared to individuals scoring lower on these measures. Specifically, we predicted that participants’ scores on the NEO-FFI Openness domain, the modified tellegen absorption scale (MODTAS), and the dispositional positive emotions scale (DPES) would show significant positive correlations with both the likelihood of experiencing chills and the self-reported intensity of chills experiences across stimuli.
The study employs a cross-sectional correlational research design to comprehensively investigate the underlying factors influencing individual variability in the experience of aesthetic chills. Its primary objective is to unveil the intricate interplay of psychological and sociocultural dynamics in shaping chills reactions in a diverse cohort of N ≍ 3,000 participants from Southern California. Through a multimethod approach, validated scales assessing demographics, personality traits, emotional tendencies, and psychological characteristics are administered to participants. Subsequently, participants are exposed to a validated set of 40 emotion-evoking audiovisual stimuli. The study employs a two-step combination of categorical and continuous data analysis techniques, including (i) univariate and multivariate regression, LCA, and (ii) LASSO regression, and SVM (see the Statistical methods section) within a cross-validation framework. These analyses collectively enable the exploration of relationships between demographics, personality, and context and the likelihood of experiencing chills. The study's cross-sectional nature enables the identification of associations without manipulating variables, contributing to a comprehensive understanding of the complex interplay between psychological traits and sociocultural influences in shaping the subjective experience of aesthetic chills.
Participants were recruited for this study through an online platform (Qualtrics.com) with a focus on individuals residing in Southern California. Before proceeding with the study, participants underwent an initial screening to confirm their geographical location and provide their informed consent. Participants were then asked to provide demographic information including gender, education level, and age. Additionally, participants were queried about their political orientation and whether they were affiliated with any political party. In order to assess participants’ emotional state, they were prompted to indicate their levels of valence and arousal. Subsequently, participants completed several questionnaires including the disposition positive affect (DPES), NEO five-factor inventory (NEO-FFI), MODTAS, and Kama Muta questionnaire (KAMF). Participants were then pseudorandomly assigned to one of 40 stimuli. After exposure to the randomly assigned stimulus, participants were asked to report their emotional state in terms of valence and arousal once again. They were also asked to indicate their level of liking for the video, whether they had seen the video previously, whether they experienced chills while watching the video, and if so, to rate the frequency and intensity of their chills. Participants were also asked whether the video reminded them of a personal experience, and if they experienced goosebumps or tears, they were asked to indicate what elicited those responses. Upon completion of the study, participants were thanked for their participation and provided with appropriate remuneration for their time and effort. The average duration of each experiment was ∼37 min.
A total of 3,259 participants initially took part in the experiment. Following data cleanup (see details in Supplementary material), the experiment involved a diverse group of 2,937 participants, all of whom hailed from Southern California. The gender distribution was fairly balanced, with 54.24% identifying as female and 41.44% as male. In terms of political affiliation, the largest group identified as Democrats (50.66%), followed by Republicans (21.59%), and Independents (14.81%). Notably, 11.64% of participants did not specify a political affiliation, while a small proportion (1.19%) identified with other political affiliations. In regards to racial identity, the majority of participants (68.44%) identified as White or Caucasian. The second largest racial group was those identifying as other (11.37%), followed by American Indian/Native American or Alaska Native (4.97%), and Black or African American (1.46%). A small percentage of participants (0.31%) identified as Asian. Given this broad demographic range, this dataset provides a rich, representative sample for examining the phenomena under investigation in the context of Southern California.
The DPES (45) measures one's dispositional tendencies to feel positive emotions towards others in their daily lives. The DPES consists of seven subscales, including joy, contentment, pride, love, compassion, amusement, and awe. The scale comprises a total of 38 items, with each subscale containing five or six items. Participants rate their agreement with statements such as “On a typical day, many events make me happy” on a seven-point Likert scale ranging from “1—strongly disagree” to “7—strongly agree.” The total score is obtained by averaging the item responses, yielding a range from 1 to 7, with higher scores indicating greater levels of positive emotion.
The NEO five-factor inventory (NEO-FFI-3) (46) is a widely used personality assessment tool based on the five-factor model (FFM) of personality. It is designed to measure five broad dimensions of neuroticism, extraversion, openness to experience, agreeableness, and conscientiousness. The NEO-FFI-3 consists of 60 items, with 12 items for each of the five personality factors. Participants respond to statements such as “I worry a lot” or “I am talkative” on a Likert-type scale, indicating the extent to which each statement reflects their own personality traits. The NEO-FFI-3 provides researchers and practitioners with a concise and reliable measure of the key dimensions of personality, allowing for a comprehensive understanding of individuals’ typical patterns of behavior, emotions, and cognition. The constituent items of NEO-FFI include negative affect, self-reproach, positive affect, sociability, activity, aesthetic interests, intellectual interests, unconventionality, nonantagonistic orientation, prosocial orientation, orderliness, goal-striving, and dependability.
The MODTAS (47) is a 34-item multidimensional measure that assesses imaginative involvement and the tendency to become mentally absorbed in everyday activities (48). We employed a modified version, the MODTAS (47), which has a Likert-scaled response format. Subjects are asked to rate the frequency of their experiences on a five-point scale, ranging from 0 (never) to 4 (very often). MODTAS has been found to have a clearer covariance structure than the original tellegen absorption scale (TAS). It consists of five correlated primary synesthesia, altered states of consciousness, aesthetic involvement in nature, imaginative involvement, and apparent experiences of extra-sensory perception.
The KAMF (49–51) is a seven items scale measures predisposition for Kama Muta (Sanskrit, “moved by love”) an emotion described as “being moved,” “heart-warming,” “stirring,” or “being emotionally touched.” It is a primarily positive emotion that is experienced as a feeling of buoyancy, security, and warmth in the chest, and may be accompanied by goosebumps or tears.
A crowdsourced database of chills-eliciting stimuli was constructed as a specialized open-access database tailored to the goals of this study (41). Leveraging online social media platforms, including YouTube and Reddit, a Python-based tool was developed to curate a collection of stimuli with the potential to elicit aesthetic chills. A breadth-first search algorithm (52) navigated the platforms, using a predefined dictionary of somatic markers such as “frisson,” “chill,” “goosebump,” and “gooseflesh” to identify these markers within user comments. Only videos with over 100 comments and at least 10 occurrences of these markers were deemed suitable for inclusion. This process resulted in the acquisition of 204 potential stimuli, of which the top 50, based on the highest somatic marker frequency, were selected for further analysis. Online validation was conducted with over 3,500+ participants, involving recruitment from Prolific, Qualtrics, and an online interface with embedded attention checks to ensure data quality. Ultimately, this systematic approach yielded a curated set of 40 stimuli, equally divided between audio and audiovisual formats (see details in Ref. (41)).
We employed a two-step comprehensive, stepwise approach to identify the factors influencing aesthetic chills (see Table 1 for overview). The progression ensures we accurately isolate specific influences like psychological traits or stimulus types, thereby avoiding confounding effects such as mistakenly attributing chills to general emotional responses instead of specific stimuli characteristics.
Fig. 2. Clustering of beta coefficients across pre-stimulus logistic regression analyses of most predictive traits against chills likelihood (a-d). Absorption (MODTAS) and prestimulus arousal seem to differentially drive chills likelihood than extraversion and KAMF, i.e. while all are strong predictors when considering all stimuli together, either absorption/arousal, or extraversion/KAMF dominantly drive chills in specific sets of videos, indicating possible thematic clusters appealing to different axes of individual variation.
Note that in order to account for multiple testing in those three first sections, we employed the Bonferroni correction method. Eighty-seven tests were executed for the main analysis (61 for both parts 1 and 2; 26 for part 3), and tests were considered significant for P-values < 0.0005. Furthermore, we included additional analyses using the standard 0.05 significance threshold in the Supplementary material. These comprised 172 tests (145 for regression analysis and 27 for LCA). For these supplementary analyses, we did not apply the Bonferroni correction, as they served as sensitivity analyses to support the main findings. These analyses were primarily regression analyses where we varied the fitting models, rather than examining statistical associations between variables.
This procedure was executed 10,000 times, with the RMSE averaged across each iteration. The process was replicated over a range of 100 lambda values spaced logarithmically between 0.0001 and 100. To mitigate overfitting, the lambda value (λ) yielding the lowest RMSE was then utilized in a leave-25%-out, 10,000-fold cross-validation, which generated the results reported herein. Mirroring the SVM efforts above, we also repeated these efforts with the three levels of trait summary metrics (individual responses, subscale, and highest summary level) and subjected resultant P-values to a Bonferroni corrected threshold (P < 0.017) to warrant subsequent AIC and BIC tests for parsimony. Beta coefficients allocated to each feature during each fold were averaged across the iterations and ordered by their absolute values to gauge their relative significance. The top five beta coefficients for the most parsimonious model are reported.
The study was conducted in compliance with the Helsinki Declaration and was granted an exemption following review by Advarra IRB (Pro00068209). Participants provided informed consent and were informed about the purpose of the research, their right to decline or withdraw, and confidentiality limits. Participants were also given the opportunity to ask questions and receive answers about the phenomenon under study.
Using a classical logistic regression model, the stimulus “Hallelujah Choir (Audio)” had the highest odds ratio (OR) for inducing chills at 6.09 (95% CI: 3–12.84; P < 0.01; Fig. 3). Close contenders were “think too much feel too little” with an OR of 3.29 (95% CI: 1.68–6.57) and “Great Dictator (Audio)” with an OR of 3.22 (95% CI: 1.65–6.43), both also significant at P < 0.01. When comparing audio-only and audiovisual stimuli had an OR of 1.14 (95% CI: 0.99–1.32), which was not statistically significant (P = 0.07), indicating no significant difference in the likelihood of experiencing chills between audio and video stimuli.
Fig. 3. Stimuli probability of eliciting chills a) and correlation heatmaps of psychological scales b). The bar chart illustrates the log odds of each stimulus eliciting chills when compared to “Motorcycle Diaries (Audio)” (the lowest chills ratio stimulus), which is used as the reference stimulus. The gradient color represents a visual ranking from top to bottom, with stimuli at the top having the highest probability of eliciting chills.
As described in the Materials and methods section, in order to account for the previous results and identify factors associated with chills independently of stimulus presentation, we employed a mixed-effects regression. Higher arousal levels showed increased odds of experiencing chills, with a univariate OR of 1.19 (P < 0.01). Similarly, higher valence levels were associated with an increased likelihood, displaying a univariate OR of 1.1 (P < 0.01). Elevated mood also significantly predicted chills, with a univariate OR of 1.49 (P < 0.01). Prior exposure to the stimulus drastically increased the odds, as shown by a univariate OR of 3.75 (P < 0.01) and a multivariate OR of 2.2 (P < 0.01). Overall, higher arousal, valence, and mood, as well as prior exposure to the stimulus, were strongly associated with a greater likelihood of experiencing aesthetic chills (Fig. 4).
Fig. 4. Influence of state and prior exposure on aesthetic chills.
Age groups 35–44 showed an increased OR of 2.1 in the univariate model (P < 0.01) and an OR of 1.36 in the multivariate model (P = 0.07). Individuals with a graduate or professional degree also had higher odds (univariate OR = 2.81, P < 0.01; multivariate OR = 1.15, P = 0.42; see Table 2). Males were more likely than females to experience chills (univariate OR = 1.85, P < 0.01; multivariate OR = 1.44, P < 0.01). Finally, participants who identified as Latino had increased odds in the univariate model (OR = 1.39, P < 0.01), although this was not significant in the multivariate model (P = 0.06). Overall, being in the age group of 35–44, having a higher educational level, and being male were identified as strong predictors for experiencing aesthetic chills.
In summary, being a middle-aged adult, having a graduate or professional degree, being male, identifying as Latino, and having democratic political preferences are independent factors that increase the likelihood of experiencing aesthetic chills, although not all of these remain significant in multivariate analysis (see Table 2).
Among the NEO-FFI personality traits, higher levels of extraversion were associated with increased odds of experiencing chills, with a univariate OR of 1.11 (P < 0.01) and a multivariate OR of 1.56 (P < 0.01). Similarly, higher levels of conscientiousness were linked to an increased likelihood, displaying a univariate OR of 1.07 (P < 0.01). Agreeableness showed a negative association with a univariate OR of 0.93 (P < 0.01). Other psychological characteristics like synesthesia (OR = 2.53, P < 0.01), Asc (OR = 2.62, P < 0.01), aesthetic involvement (OR = 2.49, P < 0.01), imaginative involvement (OR = 3, P < 0.01), and DPES total score (OR = 1.04, P < 0.01) also significantly increased the likelihood of experiencing chills. In summary, higher levels of Extraversion, Conscientiousness, and certain psychological characteristics were strongly associated with a greater likelihood of experiencing aesthetic chills (Fig. 5).
Fig. 5. Contribution of personality to chills. The plot illustrates a cube constructed in a three-dimensional space with dimensions corresponding to MODTAS, extraversion, and DPES values, where data points colored in blue (representing “yes”) and red (representing “no”) indicate the presence of two distinct clusters.
Demographically, this class was predominantly aged between 35 and 44 years. A significant 77.2% of individuals in this class held a graduate or professional degree, distinguishing them educationally (P < 0.01). In terms of personality traits, class D scored notably higher in extraversion (42.8), and conscientiousness (28.1), compared to the overall sample averages (P < 0.01). Moreover, this class led in specific psychological scales, most notably achieving the highest total MODTAS score of 128.6, significantly higher than the overall average (P < 0.01). These results are summarized in Table 3.
This class consisted of 311 participants and scored the highest on the KAMF scale with a mean of 4.9 (SD = 1.3). On the MODTAS scale, class 7 also led in synesthesia (4.2), aesthetic involvement (4.2), and imaginative involvement (4.3). In terms of the DPES scores, this class had the highest mean scores in the categories of pride (31.3, SD = 12.3) and awe (49.0, SD = 20.8). Demographically, the age group of 35–44 made up 40.5% of this class. Individuals with a graduate or professional degree constituted 28.0% of the class, and 18.6% identified as Republican. When comparing sociodemographic and psychological traits across the seven groups, significant differences were observed in all these characteristics, each showing a P-value of less than 0.01 (see Table S10).
Using the feature set described above (highest-level summary trait metrics; 15 total features), the classification scheme yielded a highly significant prediction accuracy (68.2%; P < 0.001). Changing the feature set to include each trait's subscales (31 features), where available, did appear to quantitatively improve classification accuracy (70.5%; P < 0.001). Leveraging each trait's individual questionnaire responses (141 features) resulted in the highest classification accuracy (74%; P < 0.001). All three tests survived Bonferroni correction (P < 0.017). We employed AIC and BIC to compare two classification models, incorporating both accuracy and the number of features used. The criteria were calculated using the models’ classification accuracy as a surrogate for likelihood, the number of observations, and the number of features. The original feature set was deemed to be preferred over the more extensive feature sets using both AIC and BIC comparisons. Accordingly, subsequent investigations into the importance of features utilized classification results with the original feature set composed of the highest-level summary trait metrics. The RFE revealed that no single feature significantly impacted the accuracy of the model according to the AIC and BIC (Table 4), indicating a robust, multivariate dependency within the model.
The exhaustive, exploratory feature-set combination analysis resulted in a total of 2,054 analyses, representing the number of unique combinations across all demographic, trait, and state measures available. The highest classification accuracy observed was 73.5% (P < 0.001) when age, sex, arousal, KAMF, and MODTAS were included in the feature set. Both AIC and BIC comparisons revealed that the model with these five features yielded significantly higher classification values than the model above that used all the highest-level trait metrics (15 total features). The five highest classification accuracy from this feature-set search are reported alongside the features included are displayed in Table 5. All reported values passed the Bonferroni threshold identified for this exhaustive analysis (P < 2 × 10^−5^). Upon a preliminary review, it's evident that while the feature sets have distinct compositions and yield comparable classification accuracies, only KAMF, MODTAS, and Sex consistently appear across the highest-performing sets, suggesting their potential cruciality for successful classification.
In line with the predefined exclusion criteria, each fold of the leave-25%-out-cross-validation comprised 2,097 participants in the training set and 700 participants in the test set. The lambda value (λ) that minimized the RMSE was 0.248. Utilizing the feature set defined above accounted for a significant portion of the variance (R = 0.482; R2 = 0.232, P < 0.00001) with a meaningfully low RMSE of 22.01, suggesting that, on-average, the prediction was off by ∼22% in predicting an individual's chills intensity. Mirroring the SVM efforts and retraining the model with trait subscales, where available, did not improve the model's predictive power (R = 0.4579; RMSE = 22.354). Neither did employ individual questionnaire responses (R = 0.4293; RMSE = 23.1480). All these models produced P-values that passed Bonferroni correction, but the original model (with summary-level metrics) was the most parsimonious. The highest-ranking features from the utilized feature set (age, sex, education, political preference, arousal, valence, mood, DPES, KAMF, MODTAS, and NEO-FFI) in the LASSO model by absolute beta values were MODTAS (β = 6.2878), sex (β = 3.205), KAMF (β = 2.635), political preference (β = −1.081), and prestimulus-exposure mood (β = 1.08).
This study utilized a multilevel statistical method to provide a comprehensive understanding of the factors contributing to the experience of aesthetic chills. Among the stimuli tested, the “Hallelujah Choir (Audio)” was notably effective in eliciting chills. Moreover, age, educational level, and gender were identified as strong predictors of experiencing this phenomenon. Crucially, arousal at baseline seems to be a key predictor confirming the findings of Mori and Iwanaga (39). The results were consistent whether applying the classical or mixed method, even when introducing a fixed effect on prior exposure. The LCA revealed hidden classes of individuals who are significantly more likely to experience aesthetic chills.
Certain latent classes—specifically, class D in the demographic LCA and class 7 in the psychological LCA—exhibited a notably high prevalence of experiencing aesthetic chills. Class D, which consisted predominantly of individuals aged between 35 and 44 years and with a higher educational level, had an 88.2% rate of experiencing aesthetic chills. This class also scored significantly higher in personality traits such as extraversion and conscientiousness and led in the MODTAS psychological scale. Similarly, class 7, which consisted of 311 participants, led to experiencing chills at a rate of 86.5%. This class also scored highest in multiple psychological scales, including KAMF, and in the MODTAS subscales of synesthesia, aesthetic involvement, and imaginative involvement. Interestingly, the age group of 35–44 also constituted a significant portion of class 7. These observations may point to a specific profile of individuals who are more predisposed to experiencing aesthetic chills, suggesting the role of certain demographic and psychological factors in this unique emotional experience. The high prevalence of aesthetic chills in these hidden classes raises important questions for future research, especially in understanding the underlying mechanisms that contribute to such heightened emotional experiences. These findings could also be instrumental in clustering patients who are most likely to benefit from chills-based interventions and providing each patient with personalized exposure.
Interestingly, class D and class 7 not only exhibited a high prevalence of aesthetic chills but also scored significantly higher on psychological scales related to emotional and aesthetic involvement, such as MODTAS and KAMF. It is plausible that individuals in these classes have a unique way of processing emotional and sensory information, making them more receptive to the aspects of stimuli that elicit chills (24, 30, 31, 34, 35, 37, 38).
Preliminary evidence suggests that musical chills are mediated by interoceptive awareness (53–55), and subjects experiencing musical chills display improved emotional awareness compared to those who did not, speaking to the role of interoceptive inference in the chills process (7, 55). To further elucidate this phenomenon, future research could explore the neural mechanisms underlying these differences in emotional processing. Investigating whether these individuals have distinct neural pathways or activations in response to aesthetic stimuli could offer valuable insights into why they experience chills more frequently and intensely. Sachs et al. (32) found that individual differences in chills and musical reward sensitivity are tied to white matter connectivity between sensory processing areas in the superior temporal gyrus and emotional and social processing areas in the insula and medial prefrontal cortex. Another study by Williams et al. (34, 35, 38) found that participants reporting greater proneness to aesthetic chills responses exhibited significantly higher connectivity between the default network and sensory and motor cortices, higher connectivity between the ventral default and salience networks, and decreased connectivity between the cerebellum and somatomotor cortex. The susceptibility to aesthetic chills in these latent classes could be indicative of underlying neural mechanisms that facilitate greater sensory access to emotional and social reward systems (22, 24).
Building on McCrae's (30) suggestion that aesthetic chills may serve as a universal marker of “openness to experience,” our LCA instead revealed different personality traits of extraversion and conscientiousness. These traits, along with elevated scores on psychological scales like MODTAS and KAMF, could serve as more nuanced indicators of “openness.” Indeed, the finding that MODTAS and KAMF were the two trait metrics contained in the feature set that yielded the highest predictive accuracy for the occurrence of chills suggests that, in this context, these measures are more sensitive than the NEO-FFI. Silvia and Nusbaum (31) also found that openness to experience and expertise in the arts are significant contributors to experiencing aesthetic chills. Given that both aesthetic chills and Extraversion are mediated by dopaminergic reward pathways (9, 12, 24, 25), it stands to reason that class D's heightened propensity for chills is underpinned by these dopaminergic mechanisms (56).
The latent classes can be likened to distinctive microcultures within the larger population of Southern California. Remarkably, the classes that exhibit heightened chills sensitivity are characterized by specific demographic and psychological traits that are commonly observed in this region. The concept of “cultural resonance” might explain how the demographic and psychological attributes of individuals in these latent classes align with their surrounding culture (57–59). This would explain why these individuals may tend to manifest greater extroversion and heightened tendencies for immersion in daily experiences. Cultural resonance also potentially facilitates the emergence of positive emotions and strong reactions to cultural artifacts, as their individual attributes and values correspond with the prevalent cultural norms of their social context (60–64).
While the application of LCA provided robust insights into the subpopulations more predisposed to aesthetic chills, the study's limitations warrant cautious interpretation. First, there could be diluted effects from averaging together heterogeneous forms of chills. Indeed, the use of a mixed-effects model typically tends to attenuate the estimates of each regression coefficient since these models are more restrictive compared to conventional models, and this strengthens our confidence in the independent association of each parameter with chills. However, the geographical confinement to Southern California and the cross-sectional nature of the data collection constrain the extrapolation of these findings to broader demographics and over time scales. Given these caveats, future research should aim for more geographically diverse and longitudinal designs to validate and extend our observations.
Considering gender differences in personality across the ten aspects of the big five and that much larger effect sizes have been observed at the level of the aspects for both gender and political orientation differences, future work should investigate these trends further (65). It is noteworthy that multivariate methods, as opposed to univariate methods, are known to unveil larger effects when applied to gender analyses, and our findings align with this broader trend (66). Recognizing this, we incorporated this consideration into our study, employing a comprehensive statistical methodology to account for this effect.
Considering the strong genetic basis of the absorption trait (67) and Bignardi et al.'s finding (40) that 36% of the variance in feeling aesthetic chills can be explained by additive genetic factors and the remaining 64% by environmental sources of variation, the incorporation of genetic and neuroimaging methodologies could enrich our understanding of these complex emotional phenomena. Thus, while LCA serves as a powerful tool for dissecting heterogeneity in aesthetic experiences, these findings should invite more comprehensive studies including physiological measures, especially given the diluted effects from averaging together heterogenous forms of chills.
This study sheds light on the age-old question of why the same experience can evoke profoundly different reactions in different individuals. The obtained results mark a significant step forward in the development of personalized stimuli for experiments in neuroimaging and affective neuroscience. Our findings suggest that the experience of aesthetic chills is far from random; it is shaped by a complex interplay of predictable individual and cultural factors. Variable collected before stimulus exposure such as age, sex, arousal, KAMF, and MODTAS consistently prognosticated chills occurrence, irrespective of the specific stimulus presented. This suggests that knowing these factors beforehand can help predict if a stimulus, known for inducing chills, will actually have that effect on participants, thereby increasing experimental control and decreasing unexplained variability. Future research could focus on creating tailored models for individual stimuli and participant populations, aiming to best match participants with stimuli most likely to elicit chills in each individual participant. Further, the notion of “cultural resonance” to account for the chills latent class identified provides some answers to the question posed at the beginning of this article. Yes, beauty may be in the eye of the beholder, but it is also deeply embedded in the cultural and psychological fabric of the beholder. As we move forward in this line of research, the insights gained hold promise not only for understanding the complexities of human emotion but also for the clinical applications of these powerful emotional experiences.
Felix Schoeller, Institute for Advanced Consciousness Studies, Santa Monica, CA 90403, USA; Media Lab, Massachusetts Institute of Technology, Cambridge, MA 02139, USA.
Leonardo Christov-Moore, Institute for Advanced Consciousness Studies, Santa Monica, CA 90403, USA.
Caitlin Lynch, Institute for Advanced Consciousness Studies, Santa Monica, CA 90403, USA.
Thomas Diot, Department of Psychiatry, GHU Paris Psychiatrie et Neurosciences, Paris 75010, France.
Nicco Reggente, Institute for Advanced Consciousness Studies, Santa Monica, CA 90403, USA.
Supplementary material is available at PNAS Nexus online.
This work was funded by Tiny Blue Dot Foundation and the Mordechai Walder Charitable Fund.
Conceptualization: F.S., L.C.-M., N.R., C.L.; F.S., L.C.-M., N.R., C.L.; F.S., L.C.-M., N.R., C.L.; F.S.; T.D., F.S., L.C.-M., N.R.; all.
All data are available in the main text or the Supplementary material. The ChillsDB 2.0 dataset is released under a Creative Commons Attribution 4.0 International (CC BY 4.0) license on FigShare. This allows others to freely share and adapt the dataset as long as appropriate credit is given to the original creators by citing the published paper (7). The dataset is divided into two .csv files available under a CC BY 4.0 license on the associated FigShare (ChillsDB2). For a comprehensive understanding of each column, researchers are advised to refer to the Header Explanation File. All code supporting these analytical efforts is included in the following repository. Note: Requires the LibSVM toolbox. https://github.com/Institute-for-Advanced-Consciousness/E4-F01.
All data are available in the main text or the Supplementary material. The ChillsDB 2.0 dataset is released under a Creative Commons Attribution 4.0 International (CC BY 4.0) license on FigShare. This allows others to freely share and adapt the dataset as long as appropriate credit is given to the original creators by citing the published paper (7). The dataset is divided into two .csv files available under a CC BY 4.0 license on the associated FigShare (ChillsDB2). For a comprehensive understanding of each column, researchers are advised to refer to the Header Explanation File. All code supporting these analytical efforts is included in the following repository. Note: Requires the LibSVM toolbox. https://github.com/Institute-for-Advanced-Consciousness/E4-F01.