Skip to main content
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian
  • Log in
  • My Cart

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses

New Research In

Physical Sciences

Featured Portals

  • Physics
  • Chemistry
  • Sustainability Science

Articles by Topic

  • Applied Mathematics
  • Applied Physical Sciences
  • Astronomy
  • Computer Sciences
  • Earth, Atmospheric, and Planetary Sciences
  • Engineering
  • Environmental Sciences
  • Mathematics
  • Statistics

Social Sciences

Featured Portals

  • Anthropology
  • Sustainability Science

Articles by Topic

  • Economic Sciences
  • Environmental Sciences
  • Political Sciences
  • Psychological and Cognitive Sciences
  • Social Sciences

Biological Sciences

Featured Portals

  • Sustainability Science

Articles by Topic

  • Agricultural Sciences
  • Anthropology
  • Applied Biological Sciences
  • Biochemistry
  • Biophysics and Computational Biology
  • Cell Biology
  • Developmental Biology
  • Ecology
  • Environmental Sciences
  • Evolution
  • Genetics
  • Immunology and Inflammation
  • Medical Sciences
  • Microbiology
  • Neuroscience
  • Pharmacology
  • Physiology
  • Plant Biology
  • Population Biology
  • Psychological and Cognitive Sciences
  • Sustainability Science
  • Systems Biology
Research Article

Musical rhythm spectra from Bach to Joplin obey a 1/f power law

Daniel J. Levitin, Parag Chordia, and Vinod Menon
PNAS March 6, 2012 109 (10) 3716-3720; https://doi.org/10.1073/pnas.1113828109
Daniel J. Levitin
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: daniel.levitin@mcgill.ca menon@stanford.edu
Parag Chordia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Vinod Menon
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: daniel.levitin@mcgill.ca menon@stanford.edu
  1. Edited* by Dale Purves, Duke University Medical Center, Durham, NC, and approved January 10, 2012 (received for review August 31, 2011)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Abstract

Much of our enjoyment of music comes from its balance of predictability and surprise. Musical pitch fluctuations follow a 1/f power law that precisely achieves this balance. Musical rhythms, especially those of Western classical music, are considered highly regular and predictable, and this predictability has been hypothesized to underlie rhythm's contribution to our enjoyment of music. Are musical rhythms indeed entirely predictable and how do they vary with genre and composer? To answer this question, we analyzed the rhythm spectra of 1,788 movements from 558 compositions of Western classical music. We found that an overwhelming majority of rhythms obeyed a 1/fβ power law across 16 subgenres and 40 composers, with β ranging from ∼0.5–1. Notably, classical composers, whose compositions are known to exhibit nearly identical 1/f pitch spectra, demonstrated distinctive 1/f rhythm spectra: Beethoven's rhythms were among the most predictable, and Mozart's among the least. Our finding of the ubiquity of 1/f rhythm spectra in compositions spanning nearly four centuries demonstrates that, as with musical pitch, musical rhythms also exhibit a balance of predictability and surprise that could contribute in a fundamental way to our aesthetic experience of music. Although music compositions are intended to be performed, the fact that the notated rhythms follow a 1/f spectrum indicates that such structure is no mere artifact of performance or perception, but rather, exists within the written composition before the music is performed. Furthermore, composers systematically manipulate (consciously or otherwise) the predictability in 1/f rhythms to give their compositions unique identities.

  • musical structure
  • 1/f distributions
  • fractal mathematics
  • temporal perception

Musical behaviors—singing, dancing, and playing instruments—date back to Neanderthals (1), and have been a part of every human culture as far back as we know (2, 3). People experience great enjoyment and pleasure from music (4⇓–6), and music theorists have argued that this enjoyment stems in part from the structural features of music, such as the generation and violation of expectations (5, 7⇓⇓⇓–11). Mathematics has often been used to characterize, model, and understand music, from Schenkerian analysis (12, 13) to neural topography (14) and geometric models of tonality (15, 16). One particular mathematical relation that has received attention in music is the 1/f distribution, which Mandelbrot (17) termed “fractal.”

1/f distributions have been found to be a key feature of a number of natural and sensory phenomena. In analyzing the frequency of several natural disasters, including earthquakes, landslides, floods, and terrestrial meteor impacts, Hsü (18) found an inverse log-log linear (fractal) relation between the frequency and the intensity of the events:

Embedded Image

where f is the temporal frequency, M is a parameter indexing the intensity of the events, c is a constant of proportionality, and D is the fractal dimension. The well-known Richardson Effect (19, 20), which began the modern study of fractals, states that measurements of natural phenomena are characterized by 1/f noise. 1/f fluctuations are also a prominent feature of human cognition. Neurons in primary visual cortex were found to exhibit higher gain, and the spike responses exhibit higher coding efficiency and information transmission rates for 1/f signals (21), and the fluctuations of voltage across the resting membrane of myelinated nerve fibers show a 1/f spectrum (22); these findings suggest that human sensory and neural systems evolved to encode certain regularities in the physical world (23, 24), in this case, those that manifest self-similarity.

The spectral power of such signals decays exponentially with frequency (f) as [1/f]β (where β is the spectral exponent). Although large values of β (>2) indicate greater long-range correlations, and hence, highly predictable signals, very small βs (<0.5) indicate highly unpredictable signals: An extreme example is white noise (β = 0), the structure of which is entirely unpredictable (25). On the other hand, musical pitch and loudness fluctuations are known to have an intermediate range of β values (1 < β < 2), and this range has been suggested to indicate a level of predictability that is optimal to our musical experience (26, 27).

In contrast, musical rhythms, note onsets, and durations in Western classical music are considered highly regular and predictable, so much so that the term “rhythmic” has come to be synonymous with “recurring,” “regular,” or “periodic.” Because music has a beat and is based on repetition, it has been said that “what” the next musical event will be is not always easy to guess, but “when” it is likely to happen can be easily predicted (11). In fact, this comforting predictability has been suggested to underlie rhythm's fundamental contribution to the aesthetic experience of music (7, 28, 29).

It has been previously suggested that music can be characterized by a fractal geometry (30⇓–32), obeying a 1/f power law. This suggestion was demonstrated for the amplitude (26, 27) and pitch structure (33, 34) of music. Here we show that the temporal/rhythmic properties of music across a wide range of compositions also obey the 1/f law, or fractal relation.

We analyzed the spectral structure of musical rhythms (indexed by note onsets, and using a multitaper spectral analysis) to investigate whether these are indeed as regular and predictable as commonly claimed. Whereas previous researchers have studied only a handful of compositions by Bach, Mozart, and European folk songs (26, 27, 33, 34), we sought to increase the generalizability and importance of such findings by leveraging the power of a large database of notated musical compositions, spanning three centuries and dozens of composers and styles. To that end, we analyzed the scores of 1,788 musical movements from 558 compositions (independent sections of larger compositions of Western classical music). Here, we demonstrate that just like musical pitch and loudness, musical rhythms also obey a 1/f power law with β ranging from ∼0.5–1. Furthermore, we show that composers, even from the same musical era, exhibit systematically different slopes (βs) in their 1/f rhythm spectra: the resulting differences in rhythmic predictability would permit them to uniquely identify their compositions, and to distinguish their works from those of their contemporaries.

Results

To create a representation for musical rhythms in time, we used a raster representation: note onsets were represented as pulses (“spikes”) of a point process, and durations as the intervals between successive pulses. The power spectrum S(f) of the rhythm rasters was then computed for each piece using the multitaper method for point processes (35), commonly used for spectral estimation of spike trains in neurobiology datasets (Materials and Methods).

Previous observations on the relationship between the correlation structure of signals and the slopes of their 1/f spectra have been limited to continuous-time signals (such as pitch and loudness fluctuations) (25⇓–27, 36). Hence, we first verified the validity of these observations for rasterized rhythms with 1/f spectra. We generated simulated 1/f rhythm sequences with different values of the spectral exponent, β (see Materials and Methods for details). Simulated rhythm sequences with small values of β (~0, characteristic of white noise) lacked both short and long range correlations, and hence, were essentially unpredictable (Fig. 1 A and B). Sequences with intermediate values of β (~1) showed an interesting balance of short and long range correlations, and a level of predictability that was intermediate between the two more extreme cases described here (Fig. 1 C and D). Sequences with large values of β (~2) varied gradually over time, showed correlations over long timescales, and hence, were highly predictable over these timescales (Fig. 1 E and F).

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

The 1/f power law in simulated rhythm spectra. (A) (Upper) Representative rhythm sequence the spectra of which obey a power law of the form [1/f]β for β = 0.0 generated from Eq. 1. (Scale bar, 10 samples.) (Lower) Rhythm rasters (n = 50 realizations, vertical axis) for β = 0.0 (horizontal axis is time). (Scale bar, 1 s). The complete lack of correlations at all timescales is characteristic of white noise. (B) Spectrum of the rhythm rasters from A. Logarithm of spectral power is plotted against logarithm of frequency. Light gray lines represent rhythm spectrum for the individual rasters (n = 100). Dark gray points represent the average spectrum. Black line represents the linear fit to the spectrum in the frequency range of 0.01 to 1 Hz (delineated by dotted vertical gray lines). Dashed black line represents extrapolation of the linear fit to other frequencies. (C) Rhythm sequence and rasters as in A, but for β = 1.0. Variations in the rhythm and range of correlations are intermediate between that of A and E. (D) Spectrum of rasters from C. Other conventions are as in B. (E) Rhythm sequence and rasters as in A, but for β = 2.0. Both sequence and rasters show slow variations (horizontal axis) and long-range correlations. (F) Spectrum of rasters from E. Other conventions are as in B. Scale bars from A also apply to panels C and E. Spectra in panels B, D, and F are displaced arbitrarily along the y axis for clarity of presentation.

Next, we analyzed the spectral structure of musical rhythms of 1,788 movements from 558 Western classical music compositions. Selections spanned 16 genres and 44 composers (Table S1). Discarding pitch information, we extracted from each voice in each selection their rhythmic content (onsets and durations of notes and rests) that was encoded with the raster representation described above (Fig. 2A). The rhythm spectrum for an example piece (Beethoven's Quartet No. 1) (Fig. 2A) illustrates our central finding: Rhythm spectral power decreased exponentially with frequency as 1/f. In the log-domain, this meant that power decayed linearly with frequency (β = 0.85) (Fig. 2B). One might argue that any sequence of musical rhythms displays 1/f structure and that our experiment lacks a control condition. To address this argument, we conducted additional analyses in which each stimulus served as its own control, disrupting the temporal structure of the rhythm at all frequencies (Materials and Methods) by shuffling the note onsets globally, across the piece, but keeping note durations intact (37, 38). The spectrum of the resulting “shuffled” piece was flat and resembled white noise (β = 0.02) (Fig. 2C), indicating that the 1/f rhythm spectrum reflects the global structure of correlations across the entire piece, and that this structure is a consequence of the specific ordering of the rhythms, not their mere presence in the piece at random locations. The population of all 1,788 selections analyzed showed similar trends: The pooled distribution of βs was significantly different from its null (shuffled) distribution (Fig. 2C), and the median β (0.78, SD = 0.35) was not significantly different from 1. This finding demonstrates that, across compositions analyzed, 1/f spectra ubiquitously characterized the structure of musical rhythms. Converging results were obtained by two independent analysis approaches, Hurst analysis, and detrended fluctuation analysis (SI Materials and Methods).

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Musical rhythm spectra obey a 1/f power law. (A) Rasterized rhythm representation (Lower) showing note onsets extracted from Beethoven's Quartet Op. 18. No. 1 (score, Upper). The representation shown is schematic: actual durations were extracted from the Humdrum kern format (Materials and Methods). (B) (Left) The spectrum of the rhythm raster from A has power that decays linearly (in a log-scale) with frequency as 1/f (gray dots). The slope of the spectrum (spectral exponent or β) is 0.8. Colored segments show the sequence of durations (internote intervals). Black line represents the linear fit to the spectrum in the frequency range of 0.01 to 1 Hz (delineated by dotted vertical gray lines). Dashed line represents extrapolation of the linear fit to other frequencies. (Right) The spectrum of a sequence with the note onsets shuffled randomly, keeping durations intact. The shuffled spectrum is flat (β = 0.0). Other conventions are as shown (Left). (C) Distribution of rhythm spectral exponents pooled across genres (black) obtained by linear fits to individual pieces across the population of 1,788 pieces analyzed. Gray: spectral exponent distribution for the corresponding shuffled rhythms. Inverted triangles represent the distribution median. Dashed vertical line: β = 0.

Classical instrumental genres, such as quartets and sonatas (Fig. 3 A and B), and a variety of other folk, dance, and sung music spanning nearly four centuries of compositions all demonstrated 1/f spectra (Figs. S1 and S2 and Table S1). β-Distributions for these genres were significantly different from their respective null (shuffled) counterparts (P < 0.01, Wilcoxon signed rank test, Holm-Bonferroni correction for multiple comparisons) (Fig. 3 C and D and Table S1). We then analyzed how β-distributions varied across different genres of music (Materials and Methods). Symphonies and quartets had the largest β-values, indicating the most predictable rhythms (Fig. 3E and Fig. S1). On the other hand, folk and modern styles, such as mazurka and ragtime, had among the smallest β-values, indicating the least-predictable rhythms (Fig. 3E and Fig. S1).

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

The 1/f rhythm spectra are ubiquitous across genres. (A) Rhythm spectra for quartets. Average spectra (dark blue points) and linear fit (dark blue) to average spectrum in the frequency range of 0.01 to 1 Hz. Faded blue lines represent spectra of individual pieces. Gray data represent spectra of shuffled rhythms. Other conventions are as in Fig. 2B. (B) Distribution of rhythm spectral exponents obtained by linear fits to individual pieces (blue), and for the corresponding shuffled rhythms (gray). Inverted triangle represents median exponents. Dashed vertical line: β = 0. (C) Rhythm spectra for sonatas (red) and corresponding shuffled rhythms (gray). Other conventions are as in A. (D) Distribution of rhythm spectral exponents for sonatas (red) and corresponding shuffled rhythms (gray). Other conventions are as in B. (E) Distribution of rhythm spectral exponents for musical genres ordered from largest mean exponent to smallest. Larger exponents indicate correlations over longer timescales, and hence more predictable rhythms (vertical gray arrow). Circles are mean exponents, and error bars are 95% CI. Disjoint intervals indicate significantly different mean exponents (Tukey–Kramer HSD).

The β-values for compositions grouped by composer revealed a surprising trend: β-distributions varied widely across composers regardless of period or style (Fig. 4 and Table S2). Beethoven, Haydn, and Mozart demonstrated different 1/f rhythm spectra. Beethoven's rhythms exhibited the largest βs (and thus the most predictability), Mozart's βs were among the smallest (with the least predictability), whereas those of Haydn were intermediate between those of Mozart and Beethoven, and significantly different from either (Fig. 4 A, B, and E) [Tukey–Kramer honestly significantly different (HSD), P < 0.05; see Materials and Methods]. Monteverdi and Joplin (Table S2) exhibited nearly identical, overlapping rhythm β-distributions (Tukey–Kramer HSD, P > 0.05) (Fig. 4 C–E). Figs. 3E and 4E display 95% confidence intervals (CI); nonoverlapping CIs are significant at P < 0.05.

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Composers exhibit distinct 1/f rhythm spectra. (A) Average rhythm spectra for Beethoven (dark green), Haydn (violet), and Mozart (olive green): contemporary composers belonging to the Classical era (1750–1820). Other conventions are as in Fig. 3A. (B) Distribution of rhythm spectral exponents for compositions of Beethoven, Haydn, and Mozart. Color conventions are as in A. Other conventions are as in Fig. 3B. (C) Average rhythm spectra for Monteverdi (blue) and Joplin (green): composers separated by nearly three centuries of compositions. Other conventions are as in Fig. 3A. (D) Distribution of spectral exponents for compositions of Monteverdi and Joplin. Color conventions are as in C. Other conventions are as in Fig. 3B. (E) Distribution of spectral exponents for composers ordered from largest mean exponent to smallest. Spectral exponents of Haydn, for example (dotted vertical lines, 95% CI), are significantly different from those of Beethoven and Mozart (P < 0.05, Tukey–Kramer HSD). Other conventions are as in Fig. 3E.

Discussion

Most listeners find music pleasing when it creates an optimal balance of predictability and surprise (7⇓–9, 28, 29). Whereas such a balance was previously attributed exclusively to 1/f structure in musical pitch and amplitude (26, 27, 33, 34), the present results demonstrate that musical rhythms also exhibit 1/f spectral structure. The 1/f structure allows us to quantify the range of predictability, self-similarity, or fractal-like structure within which listeners find aesthetic pleasure. Critically, compositions across four centuries and several subgenres of Western classical music all demonstrated 1/f structure in their rhythm spectra.

A previous study analyzing classical compositions from the 18th to 20th centuries reported nearly identical 1/f pitch structure among composers, with a very narrow range of spectral exponents (1.79 ≤ β ≤ 1.97) (8). The 1/f pitch structure could not systematically delineate one composer's work from another. In contrast, our findings demonstrate that 1/f rhythm spectral exponents varied widely and systematically among composers (0.48 ≤ β ≤ 1.05). Even composers belonging to the same musical era (1750–1820, the Classical era), such as Beethoven, Haydn, and Mozart, demonstrated distinctive 1/f rhythm spectra. Conversely, Monteverdi and Joplin, composers of entirely different musical eras and composing nearly three centuries apart, exhibited similar rhythm spectra. These results suggest a heretofore underappreciated importance of rhythm and hint at its even greater role than pitch in conveying the distinctive style of composers.

Human perception is known to be sensitive to 1/f structure in the environment (21, 22, 25, 36, 39, 40), and electrophysiological studies have suggested a preference of sensory neural coding for 1/f signals (21, 39, 41). How does this relate to music? The balance of expectations both realized and violated in pleasurable music requires a certain amount of instability in the temporal structure, and such dynamical instability may best be modeled by a power law for temporal fluctuation (42).

Cognitive psychologists (43) have noted that human performance in a variety of tasks, including the production of rhythmic sequences by finger tapping (44), fluctuates over time according to 1/f. Here, we reveal that this same structure characterizes the stage of human cognition before action: the written composition of temporal intervals that will become action plans only at a later date when performed by musicians. Thus, 1/f is not merely an artifact of performance, but exists in the written scores themselves. Perhaps composers can't help but produce 1/f rhythm spectra, perhaps musical conventions require this. From a psychological standpoint, the finding suggests that composers have internalized some of the regularities of the physical world as it interacts with biological systems (including the mind) to recreate self-similarity in works of musical art (cf. 24).

Our finding of 1/f rhythm structure in nearly four centuries of musical compositions may therefore be rooted in a fundamental propensity of our sensory and motor systems to both perceive and produce 1/f structure (44, 45).

Finally, an important methodological contribution of our study is the development of raster representations borrowed from the neuroscience literature, combined with a multitaper spectral analysis (SI Materials and Methods) to rigorously examine rhythm structure in music. Although pitch and loudness are musical attributes that can be represented as a continuous signal having a specific value at each instant of time (36), musical rhythms are essentially sequences of note durations (intervals of time) that, by definition, cannot be represented as a continuous signal with a specific value at each point in time. Representation of rhythmic structure as point processes allows us to examine spectral and temporal dynamics with unprecedented mathematical precision and rigor. The methods developed in our study are likely to useful for addressing many other important questions involving rhythmicity in music and speech across cultures.

Materials and Methods

Analysis of Musical Scores.

All musical scores were drawn from the Humdrum Kern database (46). Voices (instruments) and durations were extracted with custom scripts using the Humdrum Toolkit (47). Musical movements, which are self-contained sections of larger compositions, were treated as independent units (pieces). Onsets for each voice of each piece were converted into a time-aligned raster representation, providing a marker of interonset intervals (interevent durations). Onsets from different voices were merged (as shown in Fig. 2A). Rhythm spectra were then computed with the multitaper approach using the Chronux toolbox (35, 48). The multitaper approach (35) minimizes the variance (and spectral leakage) of the spectral estimate by premultiplying the data with spectrally concentrated orthogonal (Slepian) tapers (35, 49, 50): here, for estimating rhythm spectra, we used three tapers and a sampling rate of 10 Hz (which is well above the Nyquist rate for our temporal frequencies of interest; see next section) (51, 52). Only movements with more than 200 notes were included in the analysis to permit robust spectral estimation (36). We fit straight lines to model the variation in rhythm spectral power over two orders of magnitude of frequency (between 0.01 and 1 Hz) (36) using robust regression. (There are obvious limits of the 1/f linear fit: at one end it is limited by note durations, and at the other end it is limited by the length of the pieces/scores we analyzed.) To generate an empirical null distribution for statistical comparison, and as a control to verify and validate our spectral estimation on rhythm rasters, note onsets were shuffled randomly and rasterized: as expected, shuffled onset rasters showed flat spectra (characteristic of white noise) (Fig. 2 B and C). (If the mere inclusion of notes of specific durations in a musical piece, but not their order, were responsible for 1/f, the random arrangements of notes would still yield this power distribution; thus the shuffled version of the pieces forms a suitable null hypothesis.)

Spectral exponents for pieces were compared with their shuffled counterparts using the Wilcoxon signed rank test (last column, Table S1). To compare exponent distributions across genres and composers, we performed a one-way ANOVA and tested for significant differences between exponent means based on the Tukey–Kramer HSD criterion. Statistical analyses were performed in Matlab. We repeated the spectral analysis for each voice extracted separately, normalizing all pieces to the same total duration, as well as including pieces of all lengths (< 200 notes), and observed results similar to those reported here. Plots of the power spectra for individual composers and genres, as well as additional details on the compositions used in this study, can be found in SI Appendix 1 and SI Appendix 2.

Generation of Simulated 1/f Rhythms.

To generate duration sequences that obey a power law of the form [1/f]β we modeled the internote intervals, or durations, (τ) as a multiplicative stochastic point process (53):

Embedded Image

In this formula, γ is a relaxation factor, ε represents a normally distributed, zero-mean, unit-variance, white-noise process, and σ represents a scaling factor for the SD of the white noise. The relaxation factor (γ) was varied to generate duration sequences with different spectral exponents (β). Other parameters were set to fixed values (μ = 0.0, σ = 0.05). Larger values of the relaxation factor led to sequences with longer (persistent) history effects: such sequences showed correlations over longer time windows (Fig. 1A), and exhibited 1/f spectra with larger βs (Fig. 1B). Duration sequences with β = 2.0 were generated with γ = 0.025, whereas sequences with β = 1.0 were generated with γ = 0.0. Sequences with β = 0.0 (white noise) were generated with γ = 0.0, but without the first term on the right hand side of Eq. 1, thus eliminating all history effects. Rhythm rasters were then derived from the duration sequences by placing “ticks” on the time axis such that successive ticks in time were separated by successive values in the duration sequence. All simulations were performed in Matlab.

Acknowledgments

We thank Evan Balaban for assistance in developing the original idea for this study and Devarajan Sridharan for his help with generation and analysis of simulated rhythm sequences, rasterized representation of the musical rhythms from the Center for Computer Assisted Research in the Humanities database, multitaper spectral analysis, statistical analysis of rhythm spectra, and with figures and tables. This work was supported by National Science Foundation Grants BCS-0449927 (to V.M. and D.J.L.) and IIS-0855758 (to P.C.), by Natural Sciences and Engineering Research Council Grant 228175-10 and a Google Faculty Research Award (to D.J.L.).

Footnotes

  • ↵1To whom correspondence may be addressed. E-mail: daniel.levitin{at}mcgill.ca or menon{at}stanford.edu.
  • Author contributions: D.J.L. and V.M. designed research; D.J.L., P.C., and V.M. performed research; and D.J.L., P.C., and V.M. wrote the paper.

  • The authors declare no conflict of interest.

  • ↵*This Direct Submission article had a prearranged editor.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1113828109/-/DCSupplemental.

View Abstract

References

  1. ↵
    1. Mithen S
    (2006) The Singing Neanderthals: The Origins of Music, Language, Mind, and Body (Harvard Univ Press, Cambridge, MA).
  2. ↵
    1. Yi SW
    1. Cross I
    (1999) in Music, Mind and Science, ed Yi SW (Western Music Research Institute, Seoul National University, Seoul), pp 35–40.
  3. ↵
    1. Cross I
    (2001) Music, cognition, culture, and evolution. Ann N Y Acad Sci 930:28–42.
    OpenUrlPubMed
  4. ↵
    1. Plato
    (1987) The Republic (Penguin, London).
  5. ↵
    1. Huron D
    (2006) Sweet Anticipation: Music and the Psychology of Expectation (MIT Press, Cambridge, MA).
  6. ↵
    1. Seashore CE
    (1967) Psychology of Music (Dover, Mineola, NY).
  7. ↵
    1. Meyer L
    (1956) Emotion and Meaning in Music (Univ of Chicago Press, Chicago).
  8. ↵
    1. Lerdahl F,
    2. Jackendoff R
    (1983) A Generative Theory of Tonal Music (MIT, Cambridge).
  9. ↵
    1. Levitin DJ
    (2006) This Is Your Brain on Music (Dutton, New York).
  10. ↵
    1. Narmour E
    (1992) The Analysis and Cognition of Melodic Complexity: The Implication-Realization Model (Univ of Chicago Press, Chicago).
  11. ↵
    1. Levitin DJ
    (2010) Why music moves us. Nature 464:834–835.
    OpenUrlCrossRef
  12. ↵
    1. Forte A,
    2. Gilbert SE
    (1982) Introduction to Schenkerian Analysis: Form and Content in Tonal Music (Norton, New York).
  13. ↵
    1. Schenker H
    (1954) Harmony, trans Borgese E (Univ of Chicago Press, Chicago).
  14. ↵
    1. Janata P,
    2. et al.
    (2002) The cortical topography of tonal structures underlying Western music. Science 298:2167–2170.
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Callender C,
    2. Quinn I,
    3. Tymoczko D
    (2008) Generalized voice-leading spaces. Science 320:346–348.
    OpenUrlAbstract/FREE Full Text
  16. ↵
    1. Tymoczko D
    (2006) The geometry of musical chords. Science 313:72–74.
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Mandelbrot BB
    (1977) The Fractal Geometry of Nature (Freeman, New York).
  18. ↵
    1. Hsü KJ
    (1983) Actualistic catastrophism: Address of the retiring president of the international association of sedimentologists. Sedimentology 30(1):3–9.
    OpenUrlCrossRef
  19. ↵
    1. Mandelbrot B
    (1967) How long is the coast of britain? Statistical self-similarity and fractional dimension. Science 156:636–638.
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Richardson LF
    (1961) The problem of contiguity: An appendix to statistic of deadly quarrels. GenSys: Yearbk Soc Adv Gen Sys Thry 6(139):129–187.
    OpenUrl
  21. ↵
    1. Yu Y,
    2. Romero R,
    3. Lee TS
    (2005) Preference of sensory neural coding for 1/f signals. Phys Rev Lett 94:108103.
    OpenUrlCrossRefPubMed
  22. ↵
    1. Verveen AA,
    2. Derksen HE
    (1968) Fluctuation phenomena in nerve membrane. Proc IEEE 56:906–916.
    OpenUrlCrossRef
  23. ↵
    1. Shepard RN
    (1994) Perceptual-cognitive universals as reflections of the world. Psychon Bull Rev 1(1):2–28.
    OpenUrlCrossRef
  24. ↵
    1. Shepard RN
    (2001) Perceptual-cognitive universals as reflections of the world. Behav Brain Sci 24:581–601, discussion 652–671.
    OpenUrlCrossRefPubMed
  25. ↵
    1. Schmuckler MA,
    2. Gilden DL
    (1993) Auditory perception of fractal contours. J Exp Psychol Hum Percept Perform 19:641–660.
    OpenUrlCrossRefPubMed
  26. ↵
    1. Voss RF,
    2. Clarke J
    (1978) “1/f noise” in music: Music from 1/f noise. J Acoust Soc Am 63:258–263.
    OpenUrlCrossRef
  27. ↵
    1. Voss RF,
    2. Clark J
    (1975) 1/f noise in music and speech. Nature 258:317–318.
    OpenUrlCrossRef
  28. ↵
    1. Bernstein L
    (1959) The Joy of Music (Simon and Schuster, New York).
  29. ↵
    1. Bernstein L
    (1976) The Unanswered Question: Six Talks at Harvard (Charles Eliot Norton Lectures) (Harvard Univ Press, Cambridge, MA).
  30. ↵
    1. Campbell P
    (1986) The music of digital computers. Nature 324:523–528.
    OpenUrlCrossRef
  31. ↵
    1. Campbell P
    (1987) Is there such a thing as fractal music? Nature 325:766.
    OpenUrlPubMed
  32. ↵
    1. Schroeder MR
    (1987) Is there such a thing as fractal music? Nature 325:765–766.
    OpenUrlPubMed
  33. ↵
    1. Hsü KJ,
    2. Hsü AJ
    (1990) Fractal geometry of music. Proc Natl Acad Sci USA 87:938–941.
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Hsü KJ,
    2. Hsü AJ
    (1991) Self-similarity of the “1/f noise” called music. Proc Natl Acad Sci USA 88:3507–3509.
    OpenUrlAbstract/FREE Full Text
  35. ↵
    1. Mitra P,
    2. Bokil H
    (2008) Observed Brain Dynamics (Oxford Univ Press, New York).
  36. ↵
    1. Boon JP,
    2. Decroly O
    (1995) Dynamical systems theory for music dynamics. Chaos 5:501–508.
    OpenUrlCrossRefPubMed
  37. ↵
    1. Abrams DA,
    2. et al.
    (2011) Decoding temporal structure in music and speech relies on shared brain resources but elicits different fine-scale spatial patterns. Cereb Cortex 21:1507–1518.
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Levitin DJ,
    2. Menon V
    (2003) Musical structure is processed in “language” areas of the brain: A possible role for Brodmann Area 47 in temporal coherence. Neuroimage 20:2142–2152.
    OpenUrlCrossRefPubMed
  39. ↵
    1. Patel AD,
    2. Balaban E
    (2000) Temporal patterns of human cortical activity reflect tone sequence structure. Nature 404:80–84.
    OpenUrlCrossRefPubMed
  40. ↵
    1. Garcia-Lazaro JA,
    2. Ahemd B,
    3. Schnupp JW
    (2006) Dynamical substructure of coordinated rhythmic movements. Curr Biol 16:264–271.
    OpenUrlCrossRefPubMed
  41. ↵
    1. Wu D,
    2. Li C-Y,
    3. Yao D-Z
    (2009) Scale-free music of the brain. PLoS ONE 4:e5915.
    OpenUrlCrossRefPubMed
  42. ↵
    1. Bak P,
    2. Tang C,
    3. Wiesenfeld K
    (1987) Self-organized criticality: An explanation of the 1/f noise. Phys Rev Lett 59:381–384.
    OpenUrlCrossRefPubMed
  43. ↵
    1. Wagenmakers E-J,
    2. Farrell S,
    3. Ratcliff R
    (2004) Estimation and interpretation of 1/falpha noise in human cognition. Psychon Bull Rev 11:579–615.
    OpenUrlCrossRefPubMed
  44. ↵
    1. Gilden DL,
    2. Thornton T,
    3. Mallon MW
    (1995) 1/f noise in human cognition. Science 267:1837–1839.
    OpenUrlAbstract/FREE Full Text
  45. ↵
    1. Schmidt RC,
    2. Beek PJ,
    3. Treffner PJ,
    4. Turvey MT
    (1991) Dynamical substructure of coordinated rhythmic movements. J Exp Psychol Hum Percept Perform 17:635–651.
    OpenUrlCrossRefPubMed
  46. ↵
    1. Huron D
    (2010) Humdrum Kern database, http://kern.humdrum.net. Accessed August 15, 2010.
  47. ↵
    1. Huron D
    (1998) Humdrum Toolkit, http://humdrum.org/Humdrum/. Accessed August 15, 2010.
  48. ↵
    (2008) Mitra-Laboratory (MatLab, Cold Spring Harbor, NY) Available at http://chronux.org/chronux/. Accessed August 15, 2010.
  49. ↵
    1. Slepian D
    (1978) Prolate spheroidal wave functions, Fourier analysis, and uncertainty – V: The discrete case. Bell Syst Tech J 57:1371–1430.
    OpenUrl
  50. ↵
    1. Thompson DJ
    (1982) Spectrum estimation and harmonic analysis. Proc IEEE 70:1055–1096.
    OpenUrlCrossRef
  51. ↵
    1. Nyquist H
    (1928) Certain topics in telegraph transmission theory. Trans AIEE 47:617–644.
    OpenUrl
  52. ↵
    1. Shannon CE
    (1937) Communication in the presence of noise. Proc Inst Radio Eng 37:10–21.
    OpenUrl
  53. ↵
    1. Kaulakys B,
    2. Gontis V,
    3. Alaburda M
    (2005) Point process model of 1/f noise vs a sum of lorentzians. Phys Rev Lett E71(5):1–11.
    OpenUrl
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Musical rhythm spectra from Bach to Joplin obey a 1/f power law
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Musical rhythm spectra from Bach to Joplin obey a 1/f power law
Daniel J. Levitin, Parag Chordia, Vinod Menon
Proceedings of the National Academy of Sciences Mar 2012, 109 (10) 3716-3720; DOI: 10.1073/pnas.1113828109

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Musical rhythm spectra from Bach to Joplin obey a 1/f power law
Daniel J. Levitin, Parag Chordia, Vinod Menon
Proceedings of the National Academy of Sciences Mar 2012, 109 (10) 3716-3720; DOI: 10.1073/pnas.1113828109
Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley
Proceedings of the National Academy of Sciences: 109 (10)
Table of Contents

Submit

Sign up for Article Alerts

Article Classifications

  • Social Sciences
  • Psychological and Cognitive Sciences
  • Physical Sciences
  • Physics

Jump to section

  • Article
    • Abstract
    • Results
    • Discussion
    • Materials and Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Abstract depiction of a guitar and musical note
Science & Culture: At the nexus of music and medicine, some see disease treatments
Although the evidence is still limited, a growing body of research suggests music may have beneficial effects for diseases such as Parkinson’s.
Image credit: Shutterstock/agsandrew.
Large piece of gold
News Feature: Tracing gold's cosmic origins
Astronomers thought they’d finally figured out where gold and other heavy elements in the universe came from. In light of recent results, they’re not so sure.
Image credit: Science Source/Tom McHugh.
Dancers in red dresses
Journal Club: Friends appear to share patterns of brain activity
Researchers are still trying to understand what causes this strong correlation between neural and social networks.
Image credit: Shutterstock/Yeongsik Im.
Yellow emoticons
Learning the language of facial expressions
Aleix Martinez explains why facial expressions often are not accurate indicators of emotion.
Listen
Past PodcastsSubscribe
Goats standing in a pin
Transplantation of sperm-producing stem cells
CRISPR-Cas9 gene editing can improve the effectiveness of spermatogonial stem cell transplantation in mice and livestock, a study finds.
Image credit: Jon M. Oatley.

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490