Fact-checked by Grok 2 weeks ago

Pink noise

Pink noise, also known as 1/f noise or flicker noise, is a random signal or process in which the power spectral density decreases inversely with frequency, such that the power per octave remains constant across the spectrum. This characteristic distinguishes it from white noise, which has equal power across all frequencies, and brown noise, which has power spectral density proportional to 1/f². The term "pink" arises from an analogy to visible light spectra, where pink light contains more lower-frequency (red) components compared to white light's uniform distribution. In terms of auditory perception, pink noise produces a balanced, natural-sounding rumble similar to steady rainfall, ocean waves, or wind through trees, with emphasis on lower frequencies that makes it less harsh than white noise. Its power spectral density follows the form S(f) \propto 1/f^\alpha where \alpha \approx 1, often exhibiting slight variations around this exponent in real-world measurements. Pink noise is prevalent in electronic systems as flicker noise, arising from material imperfections and surface effects in semiconductors, and it dominates at low frequencies in devices like resistors and transistors. Pink noise finds wide applications in engineering and audio testing, where it is used to evaluate the frequency response of speakers, amplifiers, and room acoustics due to its equal energy per octave, providing a more realistic simulation of musical content than white noise. Some studies in sleep research have suggested that timed pink noise stimulation can enhance slow-wave sleep, improve memory consolidation, and reduce brain wave complexity, potentially benefiting individuals with mild cognitive impairment, though results are mixed and evidence inconsistent. Additionally, it appears naturally in biological systems, such as heartbeat variability and neural firing patterns, reflecting underlying fractal or self-similar processes.

Fundamentals

Definition

Pink noise is a stochastic signal whose power spectral density (PSD) is proportional to $1/f, where f is the frequency, resulting in equal power distributed across each octave of frequency. This distinguishes it from other colored noises. In comparison, white noise exhibits a flat PSD, delivering equal power per hertz across the frequency spectrum. Brown noise, conversely, follows a PSD proportional to $1/f^2, concentrating greater energy at lower frequencies than pink noise. Formally, the PSD of pink noise is expressed as S(f) = C / |f|^\beta, where C is a constant and \beta \approx 1. It is also known as 1/f noise, a term commonly used in physics to describe this power-law spectral behavior.

Spectral Characteristics

Pink noise exhibits a power spectral density (PSD) that decreases inversely with frequency, following a form proportional to $1/f, where f is the frequency. This results in a roll-off of approximately -10 dB per decade in the frequency domain. In a log-log plot of the PSD, this characteristic manifests as a straight line with a slope of -1, providing a clear visual signature that distinguishes pink noise across frequency scales. The $1/f spectrum ensures an equal distribution of energy across octaves or decades, which contributes to a perceptually balanced sound, as it aligns more closely with the logarithmic sensitivity of human hearing compared to the flat spectrum of white noise. To measure the PSD of pink noise, the time-domain signal is transformed into the frequency domain using the Fourier transform, often via methods like the periodogram or Welch's averaged periodogram for improved estimation accuracy. In real-world approximations, the exponent \beta in the generalized $1/f^\beta form deviates slightly from the ideal value of 1, commonly ranging from 0.5 to 1.5 depending on the source and measurement conditions.

Properties

Amplitude and Power Spectra

The power spectral density (PSD) of pink noise follows a power-law relationship given by S(f) \propto \frac{1}{f}, where f represents the frequency and the proportionality holds for f > 0. This form arises in the idealized model of a stationary random process where the expected power per unit frequency decreases inversely with frequency. A key consequence is that the total power remains constant across logarithmic frequency intervals, such as octaves; the integrated power over the band from f to $2f is \int_f^{2f} S(g) \, dg \propto \int_f^{2f} \frac{dg}{g} = \ln 2, independent of the center frequency f. The amplitude spectrum of pink noise, denoting the magnitude of the Fourier transform components |A(f)|, scales as |A(f)| \propto \frac{1}{\sqrt{f}}. This relationship stems directly from the definition of PSD, as S(f) \propto |A(f)|^2, implying that the amplitude must decrease with the square root of the inverse frequency to yield the observed power-law decay. In the continuous-time domain, the PSD of pink noise is formally defined for a wide-sense stationary process as the Fourier transform of its autocorrelation function R(\tau):
S(f) = \int_{-\infty}^{\infty} R(\tau) e^{-j 2 \pi f \tau} \, d\tau = \frac{C}{|f|} ,
where C > 0 is a constant, assuming the process is band-limited to avoid divergence at f = 0. Deriving this from Fourier transform assumptions involves considering the expected value of the squared magnitude of the finite-time Fourier transform X_T(f) = \int_{-T/2}^{T/2} x(t) e^{-j 2 \pi f t} \, dt, normalized by the observation length T: as T \to \infty, \frac{1}{T} \mathbb{E} \left[ |X_T(f)|^2 \right] \to S(f) \propto \frac{1}{|f|}. This expectation captures the average energy distribution, but pure 1/f noise requires careful handling due to non-integrable low-frequency contributions leading to infinite total variance. For discrete-time signals, which are common in digital signal processing, the PSD is estimated from finite-length samples using methods like Welch's approach to approximate the continuous ideal. Welch's method segments the time series into overlapping windows, computes the periodogram (squared magnitude of the discrete Fourier transform) for each segment, and averages the results to reduce estimation variance:
\hat{S}(f_k) = \frac{1}{M} \sum_{m=1}^M \frac{1}{U_m} \left| \sum_{n \in \text{window } m} x w e^{-j 2 \pi f_k n} \right|^2 ,
where M is the number of segments, w is a window function, and U_m normalizes for window energy; this yields a smoothed estimate converging to S(f_k) \propto 1/f_k at discrete frequencies f_k, particularly effective for pink noise's uneven power distribution.

Statistical Distribution

Pink noise signals exhibit a probability distribution of instantaneous values that is typically Gaussian, characterized by a mean of zero and finite variance, as they are frequently modeled as linear transformations of Gaussian white noise. This Gaussian nature arises because linear filtering preserves the marginal distribution of the input noise, resulting in a stationary process where the amplitude probability density function (PDF) follows a normal distribution. However, the long-range correlations inherent to pink noise can lead to empirical distributions with heavier tails in finite-length realizations, particularly when low-frequency components dominate sample statistics. The moments of the distribution reflect this Gaussian framework: the first moment (mean) is zero by symmetry, the second moment (variance) is finite due to practical low-frequency cutoffs in the 1/f spectrum that prevent divergence, and higher-order moments—such as the third (skewness, zero) and fourth (kurtosis, approximately 3)—are influenced by the amplified low-frequency power, which can cause elevated variability in estimated moments for short signals compared to uncorrelated white noise. In natural occurrences of 1/f noise, such as in electronic conductors, these moments may show subtle deviations from ideal Gaussian values owing to underlying physical mechanisms like defect fluctuations. Empirical investigations often employ histogram analysis to assess the amplitude distribution, revealing that pink noise deviates slightly from the exact Gaussianity observed in white noise due to temporal correlations that cluster extreme values. For example, histograms of 1/f noise in thin metal films display near-Gaussian shapes but with marginally broader tails, attributed to intermittent large-amplitude excursions driven by correlated dynamics. Kolmogorov-Smirnov tests applied to synthetically generated pink noise signals confirm a strong fit to the Gaussian distribution, with test statistics typically yielding p-values above 0.05 for sufficiently long sequences, though shorter signals may fail the test due to undersampling of low-frequency contributions. These tests underscore the robustness of the Gaussian assumption while highlighting correlation-induced artifacts in practical measurements.

Autocorrelation Function

The autocorrelation function of pink noise is intimately linked to its power spectral density (PSD) through the Wiener-Khinchin theorem, which establishes that the autocorrelation is the inverse Fourier transform of the PSD. This theorem applies to wide-sense stationary processes and highlights how the 1/f PSD of pink noise results in non-trivial temporal correlations. In the one-dimensional case, the autocorrelation R(τ) for pink noise (β=1) exhibits a slow power-law decay for large lags, approximately R(τ) ∼ c / |τ| with oscillatory terms, derived from the inverse Fourier transform of the 1/f PSD with appropriate cutoffs to ensure stationarity. This slow decay reflects long-range correlations, leading to signal persistence where fluctuations tend to maintain their direction over extended periods, distinguishing pink noise from white noise's delta-like autocorrelation. For two-dimensional pink noise, as encountered in spatial signals like images, the autocorrelation under radial symmetry takes the form R(τ_x, τ_y) ∝ 1 / \sqrt{τ_x^2 + τ_y^2}, arising from the two-dimensional inverse Fourier transform of the 1/f PSD. This 1/r decay implies isotropic long-range spatial correlations, contributing to the fractal-like structure observed in natural textures modeled by pink noise.

Generation

Synthetic Methods

Synthetic methods for generating pink noise aim to produce signals with a power spectral density (PSD) that decreases as 1/f, where f is frequency, typically by processing white noise or directly constructing the spectrum in the frequency domain. One widely used approach is the Voss-McCartney algorithm, which approximates the 1/f spectrum by summing contributions from multiple white noise sources, each band-limited to successive octaves. The method maintains an array of random values (often 30 for audio-range coverage spanning about 20 kHz to 1/20 Hz) and a running total; at each step, a new uniform random number is generated and added to the total at a specific "level" corresponding to an octave, with levels updated at rates that halve with each higher octave to ensure equal power per octave. This creates the desired spectral slope without explicit filtering, making it efficient for real-time generation. The algorithm originates from work on 1/f noise in music and has been implemented in various audio software and hardware. Another common technique is the FFT-based method, which generates pink noise by creating a frequency-domain representation with amplitudes scaled by 1/√f (to yield the 1/f power spectrum) and random phases, then applying the inverse fast Fourier transform (IFFT) to obtain the time-domain signal. This method provides high accuracy over a finite bandwidth and is suitable for offline generation of long sequences, though it requires padding for periodic artifacts and is computationally intensive for real-time use. In software, these algorithms are readily implemented in libraries like Python's NumPy and SciPy or MATLAB. For the FFT-based approach in Python, the process involves generating white noise, computing its real FFT (rfft), scaling the positive frequency components by 1/√f (excluding DC), and applying the inverse real FFT (irfft) to obtain the pink noise signal, which can then be normalized. Pseudocode for this is: 
import numpy as np
from numpy.fft import rfft, irfft

def generate_pink_noise_fft(n_samples, fs=44100):
    white_noise = np.random.randn(n_samples)
    freqs = np.fft.rfftfreq(n_samples, 1/fs)
    scale = np.zeros_like(freqs)
    scale[1:] = 1 / np.sqrt(freqs[1:])
    pink_fft = rfft(white_noise) * scale
    signal = irfft(pink_fft)[:n_samples]
    return signal / np.std(signal)  # Normalize
This produces a zero-mean signal with the target spectrum. For the Voss-McCartney method, implementations often use a fixed array size and bit-shifting for efficient level updates. MATLAB equivalents employ similar randn for white noise input and filter for IIR-based approximations or direct spectral methods. Hardware generation typically involves analog circuits or digital signal processors (DSPs). In analog setups, white noise from a zener diode or transistor base is passed through a multi-stage RC low-pass filter designed for a -3 dB per octave roll-off, often using op-amps to approximate the 1/f response over the audio band (20 Hz to 20 kHz). For example, a three-pole filter with appropriately spaced time constants can achieve close to the ideal slope. Digital hardware, such as DSP chips in audio interfaces or synthesizers, implements the Voss-McCartney or IIR filtering algorithms in real-time, offering low latency and precise control.

Natural Sources

Pink noise emerges naturally in various physical systems through mechanisms that produce power spectral densities scaling as 1/f, where f is frequency. In fluid dynamics, turbulence exhibits fractal properties that result in pink noise characteristics, as the spatial and temporal fluctuations follow self-affine scaling behaviors akin to fractional Brownian motion. Similarly, fractal processes in natural phenomena, such as wind speed variations, generate pink noise due to the persistence in time series data, where high or low values tend to cluster, leading to long-range correlations. In electronic systems, flicker noise—a form of pink noise—arises from random trapping and release of charge carriers at material interfaces, causing fluctuations in carrier mobility. Specific examples illustrate these mechanisms in environmental and material contexts. River flow fluctuations often display pink noise, with spectral colors shifting from red to pink across networks due to hydrological variability and watershed scale effects, where larger basins exhibit more persistent, 1/f-like dynamics. In semiconductors and resistors, flicker noise manifests from charge trapping in defects, producing low-frequency variations that dominate over white thermal noise at sub-Hz frequencies. Self-similar structures in nature underpin many instances of pink noise through power-law spectra. These scaling laws arise from self-organized criticality, where systems evolve to critical states producing avalanches with 1/f distributions, as seen in geophysical processes like earthquakes or fluid instabilities. Such self-similarity ensures that the noise spectrum remains invariant across scales, reflecting underlying fractal geometries in the system's dynamics. Recent studies have highlighted quantum effects contributing to pink noise in advanced materials like graphene. For example, in hBN-encapsulated graphene transistors, high-field conditions amplify 1/f noise through bremsstrahlung-like quantum processes, linking it to electron-phonon interactions in two-dimensional quantum materials.

Applications

Audio Engineering and Testing

In audio engineering, pink noise serves as a fundamental test signal for evaluating and optimizing sound systems due to its equal energy distribution across octaves, which mimics the spectral balance of many natural and musical sounds. This property allows engineers to assess frequency response in a way that aligns with perceptual loudness, providing a more representative analysis than white noise, which has uniform power per hertz. For room response testing, pink noise is employed to measure the overall frequency response of acoustic environments and sound reinforcement systems. It is played through speakers while a microphone captures the output, enabling the identification of room modes, resonances, and decay characteristics via tools like real-time analyzers that display octave-band levels. Compared to swept sine methods, which generate a logarithmic frequency sweep to derive impulse responses and better isolate nonlinear distortions, pink noise offers a simpler, broadband excitation that averages the response over time but can be less precise in reverberant spaces due to overlapping frequency interactions. Swept sines are often preferred in modern acoustic measurements for their higher signal-to-noise ratio and ability to deconvolve the room transfer function, yet pink noise remains valuable for quick, on-site calibrations where computational resources are limited. In equalization processes, pink noise is used to calibrate speakers and headphones for a flat response per octave, ensuring balanced playback across the audible spectrum. By routing pink noise through the system and adjusting EQ filters until the output measures equal loudness in each octave band—typically using a spectrum analyzer—engineers achieve a response that compensates for inherent speaker irregularities and room influences. This technique is standard in professional setups, as the noise's 3 dB per octave roll-off in power spectral density naturally aligns with logarithmic frequency scaling in human hearing. The application of pink noise in recording studios dates back to the 1970s, when it became a common tool for mixer calibration and system alignment, particularly with the rise of multitrack recording and Dolby noise reduction systems. Engineers would generate pink noise bursts or continuous signals to set gain staging and balance channels, ensuring consistent headroom and frequency neutrality during mixing sessions. This practice helped standardize studio monitoring, reducing variability in playback across different environments. Professional standards from organizations like the Audio Engineering Society (AES) and the International Electrotechnical Commission (IEC) endorse pink noise as a reference signal for audio testing. For instance, AES2-2012 specifies bandpass-filtered pink noise with a 24 dB/octave Butterworth response and controlled crest factor for measuring loudspeaker performance, while IEC 60268-5 guidelines incorporate it for evaluating sound system maximum linear sound levels and frequency response in controlled conditions. These standards ensure reproducibility in tests for power handling, distortion, and overall system fidelity.

Electronics and Timekeeping

In electronic circuits, pink noise, also known as flicker noise, arises predominantly from material imperfections and carrier trapping in semiconductor devices such as transistors, manifesting as a power spectral density proportional to $1/f. This noise type significantly impacts the performance of oscillators, where it undergoes upconversion to contribute to close-in phase noise. In metal-oxide-semiconductor (MOS) transistors, the intrinsic $1/f noise has worsened with scaling, particularly in FinFET technologies, leading to enhanced phase noise in the $1/f^3 region of the spectrum due to nonlinear mixing and waveform asymmetries during oscillation. The $1/f noise from switching transistors in oscillators primarily determines the phase noise floor in the low-offset frequency regime, often dominating over thermal noise contributions. Techniques to model this include extending Kurokawa's theory by treating $1/f perturbations as small-signal dc offsets that cyclostationarily modulate the carrier, resulting in a phase noise spectrum shaped by the oscillator's impulse sensitivity function. For instance, harmonic distortion in van der Pol oscillators can further amplify flicker noise upconversion through modulation of the output voltage waveform. In precision timekeeping applications, such as atomic clocks, pink noise establishes a fundamental limit on frequency stability, appearing as a flicker frequency modulation that sets a floor in Allan variance measurements. The Allan variance, defined as \sigma_y^2(\tau) = \frac{1}{2} \langle (y_{k+1} - y_k)^2 \rangle where y_k is the fractional frequency average over interval \tau, converges to a constant value for flicker frequency noise, reflecting the $1/f power-law spectrum in frequency fluctuations. This noise floor is evident in long-term stability assessments of cesium and optical atomic clocks, where it bounds the achievable accuracy beyond which quantum projection noise or other white noise types may dominate at shorter averaging times. Allan deviation plots, the square root of the Allan variance, characteristically exhibit a flat region for pink noise-dominated regimes, corresponding to the $1/f^3 phase noise spectral density in the time domain. In such plots, the deviation \sigma_y(\tau) remains independent of averaging time \tau for flicker frequency noise, distinguishing it from white frequency noise (which decreases as $1/\sqrt{\tau}) and random walk frequency noise (which increases as \sqrt{\tau}). This signature is critical for diagnosing stability limits in ultra-stable quartz oscillators and atomic standards, as seen in NIST evaluations of primary frequency references. To mitigate $1/f noise effects in analog-to-digital converters (ADCs) used in electronic systems, techniques such as chopping and feedback stabilization are employed to shift or suppress the noise spectrum. Chopping involves modulating the input signal to higher frequencies where $1/f noise is negligible, then demodulating post-amplification, effectively removing offset, drift, and flicker components while introducing manageable ripple that can be filtered. Feedback methods, including autozeroing in precision amplifiers, periodically sample and subtract low-frequency noise, reducing $1/f contributions in dc signal chains without excessive bandwidth limitations. These approaches are particularly vital in high-resolution ADCs for instrumentation, where unmitigated pink noise can degrade effective number of bits at low input frequencies.

Biological and Physiological Contexts

In healthy human subjects, heart rate variability (HRV) exhibits a power spectral density that approximates 1/f scaling, characteristic of pink noise, reflecting adaptive cardiovascular dynamics and long-range correlations in interbeat intervals. This 1/f-like spectrum, with a slope near -1 in the low-frequency range (0.0001–0.1 Hz), indicates robust autonomic regulation and resilience to perturbations in young, healthy individuals. Electroencephalography (EEG) recordings during resting states in healthy brains often display pink noise signatures, with 1/f power spectra linked to self-organized criticality, a state optimizing information processing and network flexibility. This criticality manifests as scale-free dynamics across frequency bands, enabling efficient neural avalanches and adaptability, as observed in both scalp and intracranial EEG data. Human gait stride intervals and postural sway fluctuations in healthy adults follow 1/f noise patterns, demonstrating fractal scaling that supports stable yet flexible locomotion and balance. These long-range correlations, quantified by detrended fluctuation analysis exponents around 0.8–1.0, arise from multifractal interactions in neuromuscular control, promoting efficient adaptation to environmental demands. Deviations from pink noise patterns in biological signals serve as biomarkers for neurological disorders, particularly Parkinson's disease (PD). In PD patients, HRV shows reduced 1/f scaling and lower complexity, indicating autonomic dysfunction and impaired adaptability. EEG in PD exhibits altered criticality with steeper or shallower 1/f slopes, correlating with cognitive and motor decline up to advanced stages. Gait and postural sway in PD lose fractal 1/f dynamics, shifting toward more periodic or random fluctuations, which predicts fall risk and disease progression.

Scientific and Environmental Phenomena

In gravitational wave detection, pink noise, also known as 1/f noise or flicker noise, arises primarily from instrumental sources such as electronics, photodetectors, and laser systems in detectors like the Laser Interferometer Gravitational-Wave Observatory (LIGO). This noise manifests as a power spectral density proportional to 1/f, contributing to the low-frequency noise floor that limits sensitivity in the detection band, particularly below 100 Hz, where it can mask weak signals from astrophysical events. Seismic disturbances couple into the interferometer arms, introducing additional low-frequency components that exhibit 1/f-like characteristics through environmental vibrations, further complicating signal extraction in LIGO's strain measurements. In climate dynamics, global surface temperature records display pink noise signatures over multidecadal timescales, characterized by a power spectrum S(f) ∝ f^{-1}, reflecting intrinsic variability independent of anthropogenic forcing. Analysis of monthly-averaged temperatures from 1901 to 2012 reveals this pink noise mode as the dominant empirical orthogonal function, explaining fluctuations such as the early 21st-century global warming hiatus through natural internal dynamics. Such 1/f spectra in temperature time series, spanning decades to centuries, arise from coupled ocean-atmosphere interactions and persist across spatial scales, distinguishing them from white noise in short-term weather patterns. Anomalous diffusion in porous media often involves waiting time distributions following a power-law ψ(t) ∝ t^{-(1+α)} with 0 < α < 1, leading to subdiffusive behavior where the mean squared displacement grows slower than linearly with time, ⟨x²(t)⟩ ∝ t^α. This non-Fickian transport, observed in solute migration through heterogeneous soils or rocks, generates 1/f noise in the temporal power spectrum of particle displacements due to long trapping events in pore structures. In continuous time random walk models applied to such media, these heavy-tailed waiting times capture the scale-free persistence of immobilization, resulting in pink noise-like fluctuations that deviate from classical Brownian motion. Recent analyses of exoplanet light curves from missions like TESS and JWST have identified pink noise as a key component in variability studies, particularly from 2023 onward, where correlated red noise (γ ≈ 1.2) models stellar activity and instrumental effects to improve transit detection accuracy. In 2023 wavelet-based methods decomposed light curves into pink noise-dominated components, enabling better separation of planetary signals from host star granulation and detector systematics in high-cadence data. In JWST observations, 1/f noise primarily appears as an instrumental artifact in near-infrared detectors, requiring correction techniques such as stripe noise removal to enhance transmission spectroscopy accuracy, as demonstrated in 2024 benchmark analyses of hot Jupiters.

History and Naming

Etymology

The nomenclature for types of noise draws from an analogy to the colors of light in the visible spectrum. White noise features a flat power spectral density across all frequencies, akin to white light that contains equal intensities of all visible wavelengths. Pink noise, with its power spectral density proportional to 1/f (decreasing by 3 dB per octave), evokes pink light, which appears rose-tinted due to a relative emphasis on lower-frequency (redder) components compared to higher ones. Brown noise, exhibiting a steeper 1/f² decay, corresponds to brown light dominated by even lower frequencies, giving it a darker, more bass-heavy character. The term "pink noise" first appeared in print in 1961, introduced by audio engineer G. A. Briggs in his book A to Z in Audio, where he defined it as random noise components spanning the full audio band with equal energy per octave, distinguishing it from white noise's equal energy per hertz. This usage emerged within audio engineering contexts to describe signals useful for testing and calibration, reflecting the field's growing interest in spectrally balanced test tones during the mid-20th century. In physics and electronics, pink noise is synonymous with "1/f noise," a term highlighting its inverse frequency dependence, and "flicker noise," which originated from early observations of irregular fluctuations resembling the flickering of light in vacuum tube devices. The "flicker" label was coined in 1926 by physicist Walter Schottky, building on measurements by J. B. Johnson in 1925 of low-frequency excess noise in electronic components beyond thermal (white) noise levels. Early literature also referred to it as "excess noise" to denote its deviation from expected thermal or shot noise in systems like resistors and transistors. In oceanography and geophysics, terminology diverges slightly; a 1/f spectrum is sometimes called "pink noise," but a 1/f² spectrum—equivalent to brown noise—is more commonly termed "red noise" to describe persistent, low-frequency variability in phenomena like sea surface temperatures or atmospheric fluctuations.

Historical Development

The study of pink noise, also known as 1/f noise due to its power spectral density scaling inversely with frequency, traces its roots to early empirical observations in electronic materials during the mid-20th century. In the 1960s, researchers began documenting this type of noise in resistors and semiconductors, distinguishing it from white noise and thermal noise. A pivotal contribution came from F. N. Hooge, who in 1969 proposed an empirical law describing 1/f noise as a fundamental bulk property of homogeneous materials, where the noise voltage spectral density is proportional to 1/(f N), with N representing the total number of charge carriers. This law, known as Hooge's empirical relation, shifted understanding from surface effects to intrinsic material fluctuations and laid the groundwork for noise analysis in solid-state physics. Building on these foundations, Benoit Mandelbrot extended the theoretical framework in the 1960s and 1970s by linking 1/f noise to fractal processes and long-range dependence in stochastic time series. His seminal 1968 paper with John Van Ness introduced fractional Brownian motion models, which generate spectra approximating 1/f^α with α ≈ 1, bridging direct current signals and Brownian motion while explaining the "wild self-affinity" observed in natural fluctuations. Through the 1970s and into the 1980s, Mandelbrot's collected works further connected these models to multifractal structures, influencing fields from geophysics to economics by demonstrating how 1/f noise arises in self-similar, scale-invariant systems. This period marked a shift from empirical descriptions to mathematical formalisms, emphasizing pink noise's ubiquity in complex systems. A major theoretical milestone occurred in 1987 with Per Bak, Chao Tang, and Kurt Wiesenfeld's introduction of self-organized criticality (SOC), which provided a dynamical explanation for 1/f noise in spatially extended systems. Their sandpile model illustrated how systems naturally evolve to a critical state, producing avalanches with power-law distributions that yield flicker noise signatures without external tuning. This framework unified observations across disciplines, suggesting pink noise emerges from inherent instabilities near criticality, and inspired subsequent models in seismology, biology, and neural networks. In the 1990s, pink noise gained prominence in biological research, particularly through analyses of heart rate variability (HRV). Studies revealed that healthy HRV exhibits 1/f scaling, reflecting efficient physiological adaptability, as opposed to the more rigid white noise-like patterns in pathological states. Seminal work by the European Society of Cardiology's Task Force highlighted long-term HRV components dominated by 1/f fluctuations, linking them to autonomic nervous system dynamics and fractal regulation in cardiovascular health. This application expanded pink noise research into physiology, underscoring its role as an indicator of system resilience. Entering the 2010s, pink noise research intersected with quantum technologies, where it poses challenges and probes for decoherence in quantum information systems. Reviews emphasized 1/f noise's origins in defect traps and charge fluctuations within solid-state qubits, influencing error rates in quantum computing and necessitating noise mitigation strategies. Concurrently, post-2010 interdisciplinary advances incorporated machine learning for modeling and suppressing pink noise, with deep neural networks trained to detect 1/f signatures in time series data, enhancing signal processing in noisy environments like random telegraph signals. These developments highlight ongoing evolution toward predictive tools in quantum and data-driven contexts, addressing gaps in earlier theoretical models.