## Abstract

In this tutorial, three examples of stochastic systems are considered: a strongly damped oscillator, a weakly damped oscillator and an undamped oscillator (integrator) driven by noise. The evolution of these systems is characterized by the temporal correlation functions and spectral densities of their displacements, which are determined and discussed. Damped oscillators reach steady stochastic states. Their correlations are decreasing functions of the difference between the sample times, and their spectra have peaks near their resonance frequencies. An undamped oscillator never reaches a steady state. Its energy increases with time, and its spectrum is sharply peaked at low frequencies. The required mathematical methods and physical concepts are explained on a just-in-time basis, and some theoretical pitfalls are mentioned. The insights one gains from studies of oscillators can be applied to a wide variety of physical systems, such as atom and semiconductor lasers, which will be discussed in a subsequent tutorial.

© 2021 Optical Society of America

## 1. INTRODUCTION

A wide variety of physical systems are driven by noise. Systems of interest to us include atom and semiconductor lasers [14]. To model lasers, one has to solve equations for the upper-level electon density (or number), the optical power (or photon flux), and the optical phase. These equations are nonlinear and are driven by three different noise sources. Although power and phase fluctuations have been modeled successfully, the complexity of the mathematics obscures the underlying physics, which is relatively simple. In this tutorial, three simple stochastic systems (a strongly damped oscillator, a weakly damped oscillator and an undamped oscillator driven by noise) are studied, with only a moderate amount of mathematics. One can apply the insights gained from this study of oscillators to the aforementioned physical systems.

This tutorial focuses on systems whose evolution is governed by stochastic differential equations of the form

$${d_t}A = L(A) + R(t),$$
where $A$ is a dependent variable (an oscillator displacement or a complex wave amplitude), $L$ is a deterministic linear function of its argument, and $R$ is a random function of time (noise source). In addition to linear systems, there are many nonlinear systems for which Eq. (1) governs the initial growth of the amplitude, or perturbations of the steady-state amplitude. Because the amplitude is driven by noise, it too is a random function of time (or has a random component). The properties of the amplitude are characterized by its first two moments in the time and frequency domains. The first temporal moment $M(t) = \langle A(t)\rangle$, where $\langle \;\rangle$ denotes an ensemble average, is called the mean, and the second moment $C(t,t^\prime) = \langle {A^*}(t)A(t^\prime)\rangle$ is called the two-time correlation. The variance $\langle |A(t{)|^2}\rangle - {\langle A(t)\rangle ^2}$ is specified by the equal-time correlation and the mean. Similar definitions of mean, variance, and correlation apply to amplitude vectors (in which case $L$ is a matrix).

Fourier transforms are usually defined for integrable functions, which tend to zero as time tends to infinity. However, in the examples considered, noise is always present, so the driven amplitudes are not described by integrable functions. To avoid singularities in the theory, we use the finite-time transforms

$$A(\omega) \def\LDeqtab{}= \int_0^T A(t){e^{i\omega t}}{\rm d}t,$$
$$A(t) \def\LDeqtab{}= \int_{- \infty}^\infty A(\omega){e^{- i\omega t}}{\rm d}\omega /2\pi .$$
In one interpretation of these equations, $A(t)$ is a Fourier series, in which the frequencies are separated by $2\pi /T$. However, if this frequency difference is much smaller than any frequency of interest, the series can be replaced by an integral. In another interpretation, $A(\omega)$ is the infinite-time transform of the truncated function ${A_T}(t) = A(t)H(t)H(T - t)$, where $H$ is the Heaviside step function. This tutorial is based on the latter interpretation. Although we introduced finite integration times to regularize the mathematics, they represent reality, in the sense that measurements of physical systems are taken over finite time intervals, so one can think of $T$ as the measurement time.

Some properties of Fourier transforms are reviewed in Appendix A. Of particular importance are the $\delta$-function identity

$$\int_{- \infty}^\infty {e^{\pm i\omega t}}{\rm d}\omega /2\pi = \delta (t)$$
and the Parseval equation
$$\int_0^T |A(t{)|^2}{\rm d}t = \int_{- \infty}^\infty |A(\omega {)|^2}{\rm d}\omega /2\pi ,$$
which follows from Eqs. (2) and (4).

In this paper, $A(t)$ is dimensionless, so $A(\omega)$ has units of time. However, in some applications (electronics and photonics), $|A(t{)|^2}$ has units of power, so its time integral has units of energy, and $|A(\omega {)|^2}$ has units of energy multiplied by time (divided by frequency). For this reason, $|A(\omega {)|^2}$ is called the spectral energy density, and $|A(\omega {)|^2}/T$ is called the spectral power density. We will use the descriptive terms power and energy whenever it is helpful to do so. The Parseval equation states that the time- and frequency-domain formulas for the energy of a system are equivalent. This equivalence provides a way to check the consistency of time- and frequency-domain analyses.

The two-frequency correlation $C(\omega ,\omega ^\prime) = \langle {A^*}(\omega)A(\omega ^\prime)\rangle$, where ${A^*}(\omega)$ is an abbreviation for ${[A(\omega)]^*}$, is the double Fourier transform of the two-time correlation $C(t,t^\prime)$. In turn, the temporal correlation is the double inverse transform of the frequency correlation. Some consequences of this relation are discussed in Appendix B. In the main text, it is sufficient to consider the spectral energy density $S(\omega) = C(\omega ,\omega)$, which is specifed by the equation

$$S(\omega) = \int_0^T \int_0^T C(t,t^\prime){e^{i\omega (t^\prime - t)}}{\rm d}t{\rm d}t^\prime .$$
Thus, if the temporal correlation is known, the spectrum can be calculated (at least in principle).

This tutorial is organized as follows: In Section 2, the properties of white and colored noise are defined and discussed. In Sections 35, the behaviors of strongly damped, weakly damped and undamped oscillators driven by noise are studied. Damped oscillators attain stochastic steady states, whereas undamped oscillators do not. Formulas are derived for the temporal correlations of the displacements and their associated spectra. These formulas are nonsingular, but only because the integration (measurement) time is nonsingular. In Section 6, the specific results obtained in the preceding sections are discussed in the broader context of stochastic processes [58]. Finally, in Section 7, the main results of the tutorial are summarized. For convenience, the tutorial also contains appendices on Fourier transforms and the Wiener–Khinchin theorem [9,10].

## 2. NOISE

The most common model of noise is white noise, which is defined by the properties

$$\langle R(t)\rangle \def\LDeqtab{}= 0,$$
$$\langle R(t)R(t^\prime)\rangle \def\LDeqtab{}= 0,$$
$$\langle {R^*}(t)R(t^\prime)\rangle \def\LDeqtab{}= \sigma \delta (t - t^\prime),$$
where $\sigma$ is the source strength. $R(t)dt$ is the impulse applied to a realization of the system during the time interval $dt$. The first equation states that the average impulse is zero (positive and negative, real and imaginary impulses are equally likely). The second and third equations state that different impulses are independent, and the real and imaginary parts of the same impulse have the same strength ($\sigma /2$), but are otherwise independent.

Although this noise model is common, the assumption of zero correlation time is problematic. First, the noise power $\langle |R(t{)|^2}\rangle = \sigma \delta (0)$ is undefined (infinite), as is the associated noise energy. Second, the spectral energy density

$$\begin{split}\langle |R(\omega {)|^2}\rangle& = \int_0^T \int_0^T \langle {R^*}(t)R(t^\prime)\rangle {e^{i\omega (t^\prime - t)}}{\rm d}t{\rm d}t^\prime \\ &= \int_0^T \sigma {\rm d}t = \sigma T.\end{split}$$
The spectral density is finite, because the integration time is finite, but the spectral energy, which is the frequency integral of the density, is still infinite.

One can resolve these problems by assuming that the correlation time is very short, but nonzero. Colored noise also has properties (7) and (8). The correlation function

$$\langle {R^*}(t)R(t^\prime)\rangle = (\sigma \nu /2){e^{- \nu |t - t^\prime |}}$$
is exponentially localized and has a correlation (coherence) time of $1/\nu$. (Other localized functions are also suitable.) The normalized correlation, which is correlation (11) divided by $(\sigma \nu /2)$, is plotted as a function of the time difference $\tau = t^\prime - t$ in Fig. 1(a). The noise power $\langle |R(t{)|^2}\rangle = \sigma \nu /2$ is finite, because $\nu$ is finite, and the noise energy
$$\int_0^T |R(t{)|^2}{\rm d}t = \int_0^T (\sigma \nu /2){\rm d}t = \sigma \nu T/2$$
is finite, because $\nu$ and $T$ are both finite.

Fig. 1. Temporal correlation function (top) and spectral energy density (bottom) of colored noise. The blue density curve represents the first term in Eq. (16), whereas the red curve denotes both terms. For $\nu T = 10$, the finite-time correction is noticeable.

The spectral energy density

$$\langle |R(\omega {)|^2}\rangle = \int_0^T \int_0^T (\sigma \nu /2){e^{- \nu |t - t^\prime |}}{e^{i\omega (t^\prime - t)}}{\rm d}t{\rm d}t^\prime .$$
On the right side of Eq. (13), the $t^\prime$-integral is
$$\begin{split}\int_0^t {e^{- (\nu + i\omega)(t - t^\prime)}}{\rm d}t^\prime + \int_t^T {e^{- (\nu - i\omega)(t^\prime - t)}}{\rm d}t^\prime \\= \frac{{1 - {e^{- (\nu + i\omega)t}}}}{{\nu + i\omega}} + \frac{{1 - {e^{- (\nu - i\omega)(T - t)}}}}{{\nu - i\omega}}.\end{split}$$
The constant terms in Eq. (14) sum to $2\nu /({\nu ^2} + {\omega ^2})$, the $t$-integral of which is $2\nu T/({\nu ^2} + {\omega ^2})$. The $t$-integral of the variable terms is
$${-}\frac{{1 - {e^{- (\nu + i\omega)T}}}}{{{{(\nu + i\omega)}^2}}} - \frac{{1 - {e^{- (\nu - i\omega)T}}}}{{{{(\nu - i\omega)}^2}}}.$$
The terms proportional to ${e^{- \nu T}}$ are exponentially small, and the other terms sum to ${-}2({\nu ^2} - {\omega ^2})/({\nu ^2} + {\omega ^2}{)^2}$. This contribution is smaller than the previous one by a factor of order $\nu T$, and its frequency integral is zero (as are the integrals of the exponential terms). By combining the preceding results, one obtains the spectrum
$$\!\!\!\langle |R(\omega {)|^2}\rangle = \sigma \nu \left[{\frac{{\nu T}}{{{\nu ^2} + {\omega ^2}}} - \frac{{{\nu ^2} - {\omega ^2}}}{{{{({\nu ^2} + {\omega ^2})}^2}}}} \!\right] \sim \frac{{\sigma {\nu ^2}T}}{{{\nu ^2} + {\omega ^2}}}.\!$$
The normalized spectral density, which is density (16) divided by $\sigma T$, is plotted as a function of frequency in Fig. 1(b). Notice that spectrum (10) is constant, whereas spectrum (16) varies with frequency, which is why formulas (9) and (11) are said to define white and colored noise, respectively. The frequency integral of the Lorentzian function $1/({\nu ^2} + {\omega ^2})$ is $\pi /\nu$, so the integrated spectrum
$$\int_{- \infty}^\infty \langle |R(\omega {)|^2}\rangle {\rm d}\omega /2\pi = \sigma \nu T/2.$$
Equation (17) is identical to Eq. (12). For colored noise with an exponential correlation function, the spectrum is (asymptotically) Lorentzian and the Parseval relation is satisfied.

Notice that for white noise, the spectral power density $\langle |R(\omega {)|^2}\rangle /T = \sigma$, which is the two-sided-in-time Fourier transform of correlation function (9) with respect to the time difference $t - t^\prime$. For colored noise, $\langle |R(\omega {)|^2}\rangle /T \sim \sigma {\nu ^2}/({\nu ^2} + {\omega ^2})$, which is the transform of correlation function (11) with respect to the time difference. These results are not coincidental. They are manifestations of the Wiener–Khinchin theorem, which states that if the correlation function depends on the time difference $t - t^\prime$ (rather than $t$ and $t^\prime$ separately), then the spectral power density is (asymptotically) equal to the two-sided Fourier transform of the correlation function with respect to the time difference. This important result is proved in Appendix C.

In summary of this section, the energy provided by a noise source is only finite if the measurement (reception) time and the noise bandwidth are both finite. If the correlation time of colored noise is much shorter than any other time scale relevant to the receiving system, its effects on that system should be similar to those of white noise. For this reason, one usually analyzes the effects of white noise first, because the analysis is simpler, and checks the results for singularities. If none is present, then the white-noise model is sufficient.

## 3. STRONGLY DAMPED OSCILLATOR

Consider the response of a driven oscillator (or wave), which is governed by the displacement (or amplitude) equation

$${d_t}A = - ({\nu _r} + i{\nu _i})A + {R_a}(t),$$
where ${\nu _r}$ and ${\nu _i}$ are damping and detuning (frequency-shift) parameters, respectively. The Langevin source term ${R_a}$ [11] is a random function of time with the properties (7)–(9). For example, Eq. (18) governs the evolution of a cavity-mode amplitude in a below-threshold laser. Let $A(t) = B(t)\exp (- i{\nu _i}t)$. Then the transformed amplitude satisfies the stochastic differential equation
$${d_t}B = - {\nu _r}B + {R_b}(t),$$
where ${R_b}(t) = {R_a}(t)\exp (i{\nu _i}t)$. It is easy to verify that ${R_b}$ has the same properties as ${R_a}$. In statistical physics, Eq. (19) defines an Ornstein–Uhlenbeck process [5,8,12,13].

Equation (19) is linear in the dependent variable, so one can use the superposition principle to write its solution in the Green-function form

$$B(t) = \int_0^t G(t - s)R(s){\rm d}s,$$
where $G(t)$ is the solution of the equation
$${d_t}G = - {\nu _r}G + \delta (t),$$
together with the initial condition $G(0) = 0$. The solution of Eq. (21) is $H(t){e^{- {\nu _r}t}}$, from which it follows that the solution of Eq. (20) is
$$B(t) = \int_0^t {e^{- \nu (t - s)}}R(s){\rm d}s,$$
where $R$ and $\nu$ are abbreviations of ${R_b}$ and ${\nu _r}$, respectively. Although the impulse $R(s)ds$ is a random function of $s$, its subsequent effect on the amplitude is determined by a homogenous linear equation, which one can solve using regular calculus. The stochastic differential equations that will be considered in Sections 4 and 5 are also linear, so regular calculus will be used throughout this tutorial. For nonlinear equations driven by noise (such as those that govern laser power and phase fluctuations), one has to choose between Ito and Stratonovich calculus [5]. It is fortunate that some interesting stochastic systems can be discussed without the extra complexity of stochastic calculus.

It follows from solution (20) and property (7) that $\langle B(t)\rangle = 0$. The amplitude is a random function of time, with a mean of zero. Its variance is $\langle |B(t{)|^2}\rangle - {\langle |B(t)|\rangle ^2}$. For systems such as this one, which have means of zero, the variance equals the power $P(t) = \langle |B(t{)|^2}\rangle$. It follows from solution (20) that

$$\begin{split}\langle |B(t{)|^2}\rangle& = \int_0^t \int_0^t \langle {R^*}(s)R(s^\prime)\rangle {e^{- \nu (2t - s - s^\prime)}}{\rm d}s{\rm d}s^\prime \\ &= \sigma \int_0^t {e^{- 2\nu (t - s)}}{\rm d}s\\& = \sigma (1 - {e^{- 2\nu t}})/2\nu .\end{split}$$
Initially, the power is zero, because the oscillator starts from rest. Subsequently, the power increases monotonically with time. For times longer than the damping time, it tends to its steady-state value $\sigma /2\nu$, which is proportional to the source strength and inversely proportional to the damping rate. The steady-state power is finite, because the damping rate is nonzero. The energy
$$\int_0^T \langle |B(t{)|^2}\rangle {\rm d}t = \sigma (T/2\nu - 1/4{\nu ^2}),$$
where a term proportional to ${e^{- 2\nu T}}$ was neglected. The energy is finite because the damping rate is nonzero and the measurement time is finite. The second term on the right side of Eq. (24), which is associated with the transient power deficit, is smaller than the first term by a factor of $2\nu T$. Although this energy deficit is small, it is noticeable (not exponentially small).

The two-time correlation function

$$\langle {B^*}(t)B(t^\prime)\rangle = \int_0^t \int_0^{{t^\prime}} \langle {R^*}(s)R(s^\prime)\rangle {e^{- \nu (t + t^\prime - s - s^\prime)}}{\rm d}s{\rm d}s^\prime .$$
It follows from the Cauchy–Schwartz inequality that the magnitude of the correlation is less than or equal to the power. Because the source function is $\delta$-correlated in time, the noise applied between $t$ and $t^\prime$ does not contribute to the correlation, so one can replace the upper limits of integration by $\min (t,t^\prime)$. Suppose that $t^\prime \gt t$. Then
$$\begin{split}\langle {B^*}(t)B(t^\prime)\rangle &= \sigma \int_0^t {e^{- \nu (t + t^\prime - 2s)}}{\rm d}s\\ &= \sigma [{e^{- \nu (t^\prime - t)}} - {e^{- \nu (t^\prime + t)}}]/2\nu .\end{split}$$
Because correlation (26) is real, it must be a symmetric function of $t$ and $t^\prime$. One can make it so by replacing $t^\prime - t$ with ${|t^\prime - t|}$. Correlation (26) is illustrated in Fig. 2. The ridge on the main diagonal represents strong correlations.

Fig. 2. Two-time correlation function of a strongly damped oscillator. For times that are longer than the damping time, the correlation depends on only the time difference $t^\prime - t$.

For times ($t$ and $t^\prime$) that are longer than the damping time, the second term on the right side of Eq. (26) is negligible, and the correlation is a function of only the time difference $\tau = t^\prime - t$. The correlation does not depend on the time origin, so $C(t,t + \tau) = C(0,\tau)$. In abbreviated notation,

$$C(\tau) \sim \sigma {e^{- \nu |\tau |}}/2\nu .$$
The normalized correlation, which is correlation (27) divided by $\sigma /2\nu$, is plotted as a function of the time difference in Fig. 3(a). It decreases exponentially with the difference (positive or negative). This behavior reflects the fact that each noise impulse only affects the cumulative displacement for a time of the order of the damping time [Eq. (21)].

Fig. 3. Temporal correlation function (top) and spectral energy density (bottom) of a strongly damped oscillator. The blue density curve represents the first term in Eq. (28), whereas the red curve denotes all three terms. For $\nu T = 10$, the finite-time and transient corrections are noticeable.

The correlation function (26) has two contributions: The first (longer-lasting) term is proportional to the correlation function of colored noise, which was specified in Eq. (11). Its contribution to the spectrum was stated in Eq. (16). The contribution to the spectrum of the second (transient) term is ${-}\sigma /2\nu ({\nu ^2} + {\omega ^2})$, plus exponentially small terms. Hence, the spectrum of a strongly damped oscillator is

$$\!\!\!\langle |B(\omega {)|^2}\rangle = \frac{{\sigma T}}{{{\nu ^2} + {\omega ^2}}} - \frac{{\sigma ({\nu ^2} - {\omega ^2})}}{{\nu {{({\nu ^2} + {\omega ^2})}^2}}} - \frac{\sigma}{{2\nu ({\nu ^2} + {\omega ^2})}}.\!$$
The normalized spectral density, which is density (28) divided by $\sigma T/{\nu ^2}$, is plotted as a function of frequency in Fig. 3(b). This spectrum is also asymptotically Lorentzian. Although the transient contribution (deficit) is small, it is noticeable. [Compare Figs. 1(b) and 3(b).] The integral of the second term on the right side of Eq. (28) is zero, and the integral of $1/({\nu ^2} + {\omega ^2})$ is $\pi /\nu$. Hence, the spectral energy
$$\int_{- \infty}^\infty \langle |B(\omega {)|^2}\rangle {\rm d}\omega /2\pi = \sigma (T/2\nu - 1/4{\nu ^2}).$$
Equation (29) is consistent with Eq. (24). Although the noise spectral energy is infinite (Section 2), the displacement spectral energy is finite. Damping acts like a frequency filter, accepting only noise components with frequencies between (roughly) ${-}2\nu$ and $2\nu$. In doing so, it converts white noise to colored noise.

One can also work in the frequency domain, in which the source function has the properties

$$\langle R(\omega)\rangle \def\LDeqtab{}= \int_0^T \langle R(t)\rangle {e^{i\omega t}}{\rm d}t = 0,$$
$$\!\!\!\langle |R(\omega {)|^2}\rangle \def\LDeqtab{}= \int_0^T \!\int_0^T \!\langle {R^*}(t)R(t^\prime)\rangle {e^{i\omega (t^\prime - t)}}{\rm d}t{\rm d}t^\prime = \sigma T.\!$$
Other properties of the source function are described in Appendix B. It is a basic property of Fourier transforms that
$$\int_0^T [{d_t}B(t)]{e^{i\omega t}}{\rm d}t = B(T){e^{i\omega T}} - B(0) - i\omega \int_0^T \!B(t){e^{i\omega t}}{\rm d}t.$$
In this paper, the oscillators start from rest, so $B(0) = 0$. If they were not driven, then their responses would be transient, and $B(T)$ would tend to zero as $T \to \infty$. However, they are driven by noise, which is always present, so $B(T) \ne 0$, and one should apply the rule that the Fourier transform of ${d_t}B$ is the transform of $B$ multiplied by ${-}i\omega$ with caution.

By transforming Eq. (19), applying the initial condition $B(0) = 0$ and omitting $B(T)$, one finds that

$$B(\omega) = R(\omega)/(\nu - i\omega).$$
The convolution theorem (Appendix A) shows that the frequency-domain solution (33) is consistent with the time-domain solution (22). [In Appendix D, the inverse transform of $1/(\nu - i\omega)$ is shown to be $H(t){e^{- \nu t}}$.] By combining  Eqs. (31) and (33), one obtains the spectrum
$$\langle |B(\omega {)|^2}\rangle = \sigma T/({\nu ^2} + {\omega ^2}).$$
The spectrum in Eq. (34) is consistent with the first term in Eq. (28). It misses the second and third terms, which are of relative order $1/\nu T$. (The second term is the finite-time correction to the Lorentzian, and the third term is the transient deficit.) Thus, for a strongly damped oscillator, which attains a steady stochastic state, the standard frequency-domain analysis, which omits $B(T)$, is only asymptotically correct.

Results (27) and (28) apply to the transformed amplitude $B$. For the amplitude $A$, the corresponding results are

$$\langle {A^*}(0)A(\tau)\rangle \def\LDeqtab{}\sim (\sigma /2{\nu _r}){e^{- {\nu _r}|\tau | - i{\nu _i}\tau}},$$
$$\langle |A(\omega {)|^2}\rangle /T \def\LDeqtab{}\sim \sigma /[\nu _r^2 + {(\omega - {\nu _i})^2}],$$
respectively. Although the imaginary part of the rate coefficient shifts the phase of the correlation and the peak of the spectrum, the duration of the correlation and width of the spectrum depend on only the real part of the coefficient. Notice that formulas (35) and (36) satisfy the Wiener–Khinchin theorem.

## 4. WEAKLY DAMPED OSCILLATOR

In this section, the response of a weakly damped oscillator to a noise source is studied. The material is important for two reasons. First, oscillatory behavior is ubiquitous in science and engineering. Second, the response of a weakly damped oscillator is intermediate between the responses of a strongly damped oscillator (Section 3) and an undamped oscillator (Section 5). However, the physical results of this section are similar to those of Section 3 (after transients have died out, the oscillator displacement reaches a stochastic steady state in which the temporal correlation depends on only the time difference, and the spectrum is peaked near the natural frequency), but require much more mathematical effort to obtain. Consequently, we recommend that readers proceed directly to Section 5 and return to Section 4 later.

Let $X$ and $V$ be the oscillator displacement and velocity, respectively. Then the oscillator equations are

$${d_t}X \def\LDeqtab{}= V,$$
$${d_t}V \def\LDeqtab{}= - 2{\nu _0}V - \omega _0^2X + \epsilon R(t),$$
where ${\nu _0}$ and ${\omega _0}$ are damping and frequency parameters, respectively. In Eq. (38), $\epsilon$ is a book-keeping parameter with units of frequency, which allows the source function $R$ to have the same properties in this section as it did in the previous ones [Eqs. (7)–(9)]. The oscillator equations are similar to the ones that govern electron-number ($N$) and photon-number ($P$) fluctuations in an above-threshold laser [3,4]. $X$ corresponds to $P$ and $V$ corresponds to $N$.

Equations (37) and (38) are linear, so their solutions can be written in the forms

$$X(t) = \epsilon \int_0^t {G_x}(t - s)R(s){\rm d}s,$$
$$V(t) = \epsilon \int_0^t {G_v}(t - s)R(s){\rm d}s,$$
where the time-domain Green functions
$${G_x}(t) \def\LDeqtab{}= H(t)\sin ({\omega _r}t){e^{- {\nu _0}t}}/{\omega _r},$$
$${G_v}(t) \def\LDeqtab{}= H(t)[\cos ({\omega _r}t) - {\nu _0}\sin ({\omega _r}t)/{\omega _r}]{e^{- {\nu _0}t}},$$
and the resonance frequency ${\omega _r} = (\omega _0^2 - \nu _0^2{)^{1/2}}$. Both response functions decrease with time, but in an oscillatory manner (rather than a monotonic one). Notice that ${G_v} = {d_t}{G_x}$ for positive times (as it should do).

It follows from Eqs. (9), (39), and (41) that the squared displacement (power) is

$$\begin{split}\langle {X^2}(t)\rangle &= {\epsilon ^2}\sigma \int_0^t G_x^2(t - s){\rm d}s\\ &= ({\epsilon ^2}\sigma /2\omega _r^2)\int_0^t [1 - \cos (2{\omega _r}s)]{e^{- 2{\nu _0}s}}{\rm d}s.\end{split}$$
By integrating by parts twice, one can derive the identities
$$\!\!\!\!\!\int_a^b \cos (\alpha s){e^{- \beta s}}{\rm d}s \def\LDeqtab{}= \frac{{[\alpha \sin (\alpha s) - \beta \cos (\alpha s)]{e^{- \beta s}}|_a^b}}{{{\alpha ^2} + {\beta ^2}}},$$
$$\!\!\!\!\!\int_a^b \sin (\alpha s){e^{- \beta s}}{\rm d}s \def\LDeqtab{}= - \frac{{[\alpha \cos (\alpha s) + \beta \sin (\alpha s)]{e^{- \beta s}}|_a^b}}{{{\alpha ^2} + {\beta ^2}}}.\!$$
Notice that if $b$ is large (or infinite), then only the $a$-contributions are significant. Notice also that if $\alpha = 2{\omega _r}$ and $\beta = 2{\nu _0}$, then ${\alpha ^2} + {\beta ^2} = 4\omega _0^2$. By combining Eqs. (43) and (44), one finds that
$$\begin{split}\langle {X^2}(t)\rangle& = \frac{{{\epsilon ^2}\sigma}}{{2\omega _r^2}}\left\{{\frac{{1 - {e^{- 2{\nu _0}t}}}}{{2{\nu _0}}} - \frac{{{\omega _r}\sin (2{\omega _r}t){e^{- 2{\nu _0}t}}}}{{2\omega _0^2}}} \right.\\&\quad\left. {- \frac{{{\nu _0}[1 - \cos (2{\omega _r}t){e^{- 2{\nu _0}t}}]}}{{2\omega _0^2}}} \right\}.\end{split}$$
Initially, the power is zero. After a time of the order of the damping time $1/{\nu _0}$, the transient contributions die out, and the power attains its steady-state value
$$\langle {X^2}\rangle = \frac{{{\epsilon ^2}\sigma}}{{2\omega _r^2}}\left({\frac{1}{{2{\nu _0}}} - \frac{{{\nu _0}}}{{2\omega _0^2}}} \right) = \frac{{{\epsilon ^2}\sigma}}{{4{\nu _0}\omega _0^2}}.$$

The energy is the time integral of the power. The component integrals are

$$\def\LDeqtab{}(1/2{\nu _0})\int_0^T (1 - {e^{- 2{\nu _0}t}}){\rm d}t = T/2{\nu _0} - 1/4\nu _0^2,$$
$$\def\LDeqtab{}({\omega _r}/2\omega _0^2)\int_0^T \sin (2{\omega _r}t){e^{- 2{\nu _0}t}}{\rm d}t = ({\omega _r}/2\omega _0^2{)^2},$$
$$\begin{split}({\nu _0}/2\omega _0^2)\int_0^T [1 - \cos (2{\omega _r}t){e^{- 2{\nu _0}t}}]{\rm d}t\\\quad = ({\nu _0}/2\omega _0^2)(T - {\nu _0}/2\omega _0^2),\end{split}$$
where terms proportional to ${e^{- 2{\nu _0}T}}$ were omitted. In Eq. (48), the transient contribution is smaller than the steady-state contribution by the factor $2{\nu _0}T$, and in Eq. (50), it is smaller by the factor $2\omega _0^2T/{\nu _0} = 2{\nu _0}T{({\omega _0}/{\nu _0})^2}$. It follows from these results that the energy
$$\int_0^T \langle {X^2}(t)\rangle dt \sim \frac{{{\epsilon ^2}\sigma T}}{{2\omega _r^2}}\left({\frac{1}{{2{\nu _0}}} - \frac{{{\nu _0}}}{{2\omega _0^2}}} \right) = \frac{{{\epsilon ^2}\sigma T}}{{4{\nu _0}\omega _0^2}},$$
which is just the steady-state power multiplied by the measurement time.

It also follows from Eqs. (9), (39), and (41) that the temporal correlation

$$\begin{split}\langle X(t)X(t^\prime)\rangle &= {\epsilon ^2}\sigma \int_0^t {G_x}(t - s){G_x}(t^\prime - s){\rm d}s\\ &= ({\epsilon ^2}\sigma /2\omega _r^2)\int_0^t \{\cos [{\omega _r}(t^\prime - t)]\\&\quad - \cos [{\omega _r}(t^\prime + t - 2s)]\} {e^{- {\nu _0}(t^\prime + t - 2s)}}{\rm d}s.\end{split}$$
By making the substitution $u = t + t^\prime - 2s$, one can rewrite the first integral in Eq. (52) as
$$\begin{split}{I_1}(t,t^\prime)& = \cos [{\omega _r}(t^\prime - t)]\int_{t^\prime - t}^{t^\prime + t} {e^{- {\nu _0}u}}{\rm d}u/2\\& = \cos [{\omega _r}(t^\prime - t)][{e^{- {\nu _0}(t^\prime - t)}} - {e^{- {\nu _0}(t^\prime + t)}}]/2{\nu _0}.\end{split}$$
By repeating this substitution and using identity Eq. (44), one can rewrite the second integral as
$$\begin{split}{I_2}(t,t^\prime)& = \int_{t^\prime - t}^{t + t^\prime} \cos ({\omega _r}u){e^{- {\nu _0}u}}{\rm d}u/2\\ &= \{{\omega _r}\sin [{\omega _r}(t^\prime + t)] - {\nu _0}\cos [{\omega _r}(t^\prime + t)]\} \\&\quad \times {e^{- {\nu _0}(t^\prime + t)}}/2\omega _0^2 - \{{\omega _r}\sin [{\omega _r}(t^\prime - t)] \\&\quad- {\nu _0}\cos [{\omega _r}(t^\prime - t)]\} {e^{- {\nu _0}(t^\prime - t)}}/2\omega _0^2.\end{split}$$
It follows from Eqs. (52)–(54) that the correlation
$$C(t,t^\prime) = ({\epsilon ^2}\sigma /2\omega _r^2)[{I_1}(t,t^\prime) - {I_2}(t,t^\prime)].$$
Formula (55) was derived for positive time differences. Because the correlation is real, it must be a symmetic function of $t$ and $t^\prime$. One can make it so by replacing $t^\prime - t$ with $|t^\prime - t|$ in Eqs. (53) and (54). Correlation (55) is illustrated in Fig. 4. Notice that the diagonal ridge is flanked by two ravines.

Fig. 4. Two-time correlation function of a weakly damped oscillator for ${\nu _0}/{\omega _0} = 0.2$. For times that are longer than the damping time, the correlation depends on only the time difference $t^\prime - t$.

For times ($t$ and $t^\prime$) that are longer than the damping time $1/{\nu _0}$, the second contributions to integrals (53) and (54) can be neglected. By omitting these terms and combining the cosine terms that remain, one obtains the time-asymptotic correlation

$$C(\tau) \sim \frac{{{\epsilon ^2}\sigma}}{{4{\nu _0}\omega _0^2}}\left[{\cos ({\omega _r}\tau) + \frac{{{\nu _0}}}{{{\omega _r}}}\sin ({\omega _r}|\tau |)} \right]{e^{- {\nu _0}|\tau |}},$$
where $\tau = t^\prime - t$. The normalized correlation, which is correlation (56) divided by $({\epsilon ^2}\sigma /4\omega _0^2)$, is plotted as a function of the time difference in Fig. 5(a).

Fig. 5. Temporal correlation function (top) and spectral energy density (bottom) of a weakly damped oscillator for ${\nu _0}/{\omega _0} = 0.2$.

The spectral energy density is defined by Eqs. (6) and (56), the second of which can be rewritten in the form

$$C(\tau) \propto \frac{{{({\nu _0} + i{\omega _r}){e^{- ({\nu _0} - i{\omega _r})|\tau |}} - ({\nu _0} - i{\omega _r}){e^{- ({\nu _0} + i{\omega _r})|\tau |}}} }}{{2i{\omega _r}}}.$$
The terms to be transformed have the same form as the term in Eq. (11), with the modified damping rates $\nu = {\nu _0} - i{\omega _r}$ and ${\nu ^*} = {\nu _0} + i{\omega _r}$. Notice that $|\nu {|^2} = \omega _0^2$. Consider the first contribution to the double integral in Eq. (6). It follows from Eqs. (16) and (57) that this contribution is proportional to
$$\frac{1}{{2i{\omega _r}}}\left[{\frac{{{\nu ^*}2\nu T}}{{{\nu ^2} + {\omega ^2}}} - \frac{{\nu 2{\nu ^*}T}}{{{{({\nu ^*})}^2} + {\omega ^2}}}} \right] = \frac{{2\omega _0^2T}}{{{\omega _r}}}{\rm{Im}}\left[{\frac{1}{{{\nu ^2} + {\omega ^2}}}} \right].$$
By multiplying the right side of Eq. (58) by the factor ${\epsilon ^2}\sigma /4{\nu _0}\omega _0^2$, and using the identities ${\rm{Im}}[({\nu ^*}{)^2}] = 2{\nu _0}{\omega _r}$ and $|{\nu ^2} + {\omega ^2}{|^2} = (\omega _0^2 - {\omega ^2}{)^2} + 4\nu _0^2{\omega ^2}$, one obtains the contribution
$${S_1}(\omega) = \frac{{{\epsilon ^2}\sigma T}}{{{{(\omega _0^2 - {\omega ^2})}^2} + 4\nu _0^2{\omega ^2}}}.$$
The formula for the second contribution can be determined in the same way. However, it follows from Eq. (16) that the second contribution is smaller than the first by a factor of order $\nu T$, so it can be omitted. (In this discussion, ${\omega _0}$, ${\omega _r}$, and ${\nu _0}$ are all of the same order.)

Now consider the transient terms, which were omitted from correlation (56). When one Fourier transforms the transient term in Eq. (53), the double integral involves exponentials with the arguments

$$\begin{split}&{-}{\nu _0}(t^\prime + t) + i(\omega \pm {\omega _r})(t^\prime - t) \\&= - [{\nu _0} + i(\omega \pm {\omega _r})]t - [{\nu _0} - i(\omega \pm {\omega _r})]t^\prime .\end{split}$$
Hence, the double integrals (${\pm}$) factorize into single integrals and evaluate to $1/[\nu _0^2 + {(\omega \pm {\omega _r})^2}]$, plus exponentially small terms. When one transforms the transient term in Eq. (54), the double integral involves exponentials with the arguments
$$\begin{split}&{-}({\nu _0} \mp i{\omega _r})(t^\prime + t) + i\omega (t^\prime - t) \\&= - ({\nu _0} \mp i{\omega _r} + i\omega)t - ({\nu _0} \mp i{\omega _r} - i\omega)t^\prime .\end{split}$$
Hence, the double integrals evaluate to $1/[({\nu _0} \mp i{\omega _r}{)^2} + {\omega ^2}]$, plus exponentially small terms. Thus, the transient contributions to the spectrum are also smaller than contribution (58) by factors of order ${\omega _0}T$ and can be omitted. The normalized spectral density, which is density (59) divided by ${\epsilon ^2}\sigma T/\omega _0^4$, is plotted as a function of frequency in Fig. 5(b). Spectrum (59) has maxima of ${\epsilon ^2}\sigma T/4\nu _0^2\omega _r^2$ at $\omega = \pm {(\omega _0^2 - 2\nu _0^2)^{1/2}}$ and decreases like $1/{\omega ^4}$ as $\omega \to \infty$.

The frequency-domain versions of Eqs. (37) and (38), which are based on the standard derivative rule, can be written in the matrix form

$$\left[ {\begin{array}{*{20}{c}}{ - i\omega }&{ - 1}\\{\omega _0^2}&{ - i\omega + 2{\nu _0}}\end{array}} \right]\left[ {\begin{array}{*{20}{c}}{X(\omega )}\\{V(\omega )}\end{array}} \right] = \left[ {\begin{array}{*{20}{c}}0\\{\epsilon R(\omega )}\end{array}} \right].$$
Equation (62) has the solution
$$\left[{\begin{array}{*{20}{c}}{X(\omega)}\\{V(\omega)}\end{array}} \right] = \frac{1}{{\Delta (\omega)}}\left[{\begin{array}{*{20}{c}}{\epsilon R(\omega)}\\{- i\omega \epsilon R(\omega)}\end{array}} \right],$$
where the determinant
$$\Delta (\omega) = \omega _0^2 - 2i{\nu _0}\omega - {\omega ^2}.$$
The convolution theorem (Appendix A) shows that the frequency-domain solutions for $X$ and $V$ are consistent with the time-domain solutions (39) and (40). [In Appendix D, the inverse transform of $1/\Delta (\omega)$ is shown to be $H(t)\sin ({\omega _r}t){e^{- {\nu _0}t}}/{\omega _r}$.]

By combining Eqs. (31) and (63), one obtains the spectral energy density

$$\langle |X(\omega {)|^2}\rangle = \frac{{{\epsilon ^2}\sigma T}}{{{{(\omega _0^2 - {\omega ^2})}^2} + 4\nu _0^2{\omega ^2}}}.$$
One can calculate the spectral energy by doing a contour integral [14], the result of which is
$$\int_{- \infty}^\infty \langle |X(\omega {)|^2}\rangle {\rm d}\omega /2\pi = {\epsilon ^2}\sigma T/4{\nu _0}\omega _0^2.$$
Notice that the spectral density is nonsingular, because the measurement time is finite, and the spectral energy is finite, because the damping rate is nonzero. Even weak damping filters the noise source. By comparing Eqs. (65) and (66) to Eqs. (59) and (51), respectively, one finds that for a weakly damped oscillator, which attains a stochastic steady state, the standard frequency-domain analysis, which omits $X(T)$, is asymptotically correct. Equations (56) and (65) satisfy the Wiener–Khinchin theorem.

Solution (63) implies that $V(\omega) = - i\omega X(\omega)$. By doing another contour integral [14], one finds that the spectral energy

$$\int_{- \infty}^\infty \langle |V(\omega {)|^2}\rangle {\rm d}\omega /2\pi = {\epsilon ^2}\sigma T/4{\nu _0}.$$
In steady state, the displacement and velocity variances are the associated spectral energies divided by the measurement time. Hence,
$$\langle {V^2}\rangle = {\epsilon ^2}\sigma /4{\nu _0} = \omega _0^2\langle {X^2}\rangle .$$
The average kinetic and potential (spring) energies are equal.

For a weakly damped oscillator, most of the spectral energy is contained in the peaks near ${\pm}{\omega _0}$. Let $\omega = \pm {\omega _0} + \delta$, where $\delta$, ${\nu _0} \ll {\omega _0}$. Then, near each peak, the spectral density

$$S(\delta) \approx {\epsilon ^2}\sigma T/4\omega _0^2(\nu _0^2 + {\delta ^2}).$$
The normalized spectral density, which is density (69) divided by ${\epsilon ^2}\sigma T/\omega _0^4$, is plotted as a function of frequency in Fig. 6. Spectrum (69) has a maximum of ${\epsilon ^2}\sigma T/4\nu _0^2\omega _0^2$ and an effective width of ${\pm}{\nu _0}$, and contains the spectral energy ${\epsilon ^2}\sigma T/8{\nu _0}\omega _0^2$. [The integral of $1/(\nu _0^2 + {\delta ^2})$ is $\pi /{\nu _0}$.] The sum of the spectral energies of the two Lorentzian peaks equals the total spectral energy Eq. (66). These results suggest that the strongly-damped-oscillator equation (19) models the behavior of a resonantly driven, weakly damped oscillator.

Fig. 6. Spectral energy density of a weakly damped oscillator for ${\nu _0}/{\omega _0} = 0.2$. The blue curve represents formula (65), whereas the red curve represents the Lorentzian formula (69).

It is easy to establish this connection. By combining Eqs. (37) and (38), one obtains the second-order equation

$$(d_{\textit{tt}}^2 + 2{\nu _0}{d_t} + \omega _0^2)X = \epsilon R(t).$$
Suppose that the displacement
$$X(t) = B(t){e^{- i{\omega _0}t}} + {B^*}(t){e^{i{\omega _0}t}},$$
where the amplitude $B$ varies slowly compared to the phase factor. Then, by substituting ansatz (71) in Eq. (70), and omitting the terms $d_{\textit{tt}}^2B$ and $2{\nu _0}{d_t}B$, one finds that
$$({d_t} + {\nu _0})B \approx {R_b}(t),$$
where ${R_b} = i\epsilon R(t){e^{i{\omega _0}t}}/2{\omega _0}$. This source function satisfies an equation similar to Eq. (9), in which the source strength ${\sigma _b} = {\epsilon ^2}\sigma /4\omega _0^2 = {\sigma _v}/4\omega _0^2$. Thus, a strongly damped oscillator can be considered as the near-resonance limit of a weakly damped oscillator.

## 5. UNDAMPED OSCILLATOR

Now consider the response of an undamped oscillator (integrator) that is driven by noise. The amplitude equation is

$${d_t}A = {R_a}(t),$$
where ${R_a}$ has the properties (7)–(9). The solution of Eq. (73) is
$$A(t) = \int_0^t R(s){\rm d}s.$$
This solution can be written in the form of Eq. (20), with $G(t) = H(t)$. Solution (74) describes a Wiener process, which is also called Brownian motion [5,8]. It follows from solution (74) that the correlation function
$$\begin{split}\langle {A^*}(t)A(t^\prime)\rangle &= \int_0^t \int_0^{{t^\prime}} \langle {R^*}(s)R(s^\prime)\rangle {\rm d}s^\prime {\rm d}s\\& = \sigma \min (t,t^\prime).\end{split}$$
Noise that is supplied to the oscillator between $t$ and $t^\prime$ affects the later displacement, but not the earlier one, so it cannot affect the correlation. Correlation (75) is illustrated in Fig. 7. Unlike a damped oscillator, at no time does the correlation become a function of the time difference $t^\prime - t$, so the Wiener–Khinchin theorem does not apply to an undamped oscillator. The power $\langle |A(t{)|^2}\rangle$ increases linearly with time, and the energy
$$\int_0^T \langle |A(t{)|^2}\rangle {\rm d}t = \sigma {T^2}/2$$
is finite, because the measurement time is finite. Notice that the energy increases more rapidly with $T$ for an undamped oscillator than for a damped one [Eq. (24)], because the undamped oscillator integrates the effects of the source.

Fig. 7. Two-time correlation function of an undamped oscillator. The correlation never becomes a function of the time difference $t^\prime - t$.

Despite the fact that the oscillator response is nonstationary, the spectral energy density can be determined. It follows from Eqs. (6) and (75) that

$$\langle |A(\omega {)|^2}\rangle = \int_0^T \int_0^T \sigma \min (t,t^\prime){e^{i\omega (t^\prime - t)}}{\rm d}t{\rm d}t^\prime .$$
The $t^\prime$-integral in Eq. (77) is
$$\begin{split}&\int_0^t t^\prime {e^{i\omega t^\prime}}{\rm d}t^\prime + \int_t^T t{e^{i\omega t^\prime}}{\rm d}t^\prime\\& = \frac{{t{e^{i\omega t}}}}{{i\omega}} + \frac{{{e^{i\omega t}} - 1}}{{{\omega ^2}}} + \frac{{t\left({{e^{i\omega T}} - {e^{i\omega t}}} \right)}}{{i\omega}}\\ &= \frac{{{e^{i\omega t}} - 1}}{{{\omega ^2}}} + \frac{{t{e^{i\omega T}}}}{{i\omega}},\end{split}$$
and the $t$-integral that follows is
$$\begin{split}&\int_0^T \left[{\frac{{1 - {e^{- i\omega t}}}}{{{\omega ^2}}} + \frac{{t{e^{i\omega (T - t)}}}}{{i\omega}}} \right]{\rm d}t \\&= \frac{T}{{{\omega ^2}}} + \frac{{{e^{- i\omega T}} - 1}}{{i{\omega ^3}}} + \frac{T}{{{\omega ^2}}} + \frac{{1 - {e^{i\omega T}}}}{{i{\omega ^3}}}\\& = \frac{{2T[1 - {\rm{sinc}}(\omega T)]}}{{{\omega ^2}}},\end{split}$$
where ${\rm{sinc}}(x) = \sin (x)/x$. Multiplying this result by $\sigma$ produces the spectrum
$$\langle |A(\omega {)|^2}\rangle = \frac{{2\sigma T[1 - {\rm{sinc}}(\omega T)]}}{{{\omega ^2}}}.$$
The derivation of Eq. (80) was based on Eqs. (2) and (74). One can check it by using the alternative starting equation
$$\begin{split}A(\omega) &= \int_0^T R(s)\int_s^T {e^{i\omega t}}{\rm d}t{\rm d}s\\ &= \int_0^T \frac{{{e^{i\omega T}} - {e^{i\omega s}}}}{{i\omega}}R(s){\rm d}s,\end{split}$$
which one obtains by changing the order of integration. It follows from Eq. (81) that
$$\begin{split}\langle |A(\omega {)|^2}\rangle& = \int_0^T \int_0^T \frac{{{e^{- i\omega T}} - {e^{- i\omega s}}}}{{- i\omega}}\frac{{{e^{i\omega T}} - {e^{i\omega s^\prime}}}}{{i\omega}}\\&\quad\times \langle{R^*}(s)R(s^\prime)\rangle {\rm d}s{\rm d}s^\prime \\& = \sigma \int_0^T \frac{{{e^{- i\omega T}} - {e^{- i\omega s}}}}{{- i\omega}}\frac{{{e^{i\omega T}} - {e^{i\omega s}}}}{{i\omega}}{\rm d}s\\& = \frac{{2\sigma}}{{{\omega ^2}}}\int_0^T \{1 - \cos [\omega (s - T)]\} {\rm d}s\\& = \frac{{2\sigma T[1 - {\rm{sinc}}(\omega T)]}}{{{\omega ^2}}},\end{split}$$
as required. The alternative derivation of the spectrum is shorter than the original one.

For low frequencies, $S(\omega) \approx \sigma {T^3}/3$, which is finite because $T$ is finite, whereas for high frequencies, $S(\omega) \approx 2\sigma T/{\omega ^2}$. The spectral energy

$$\int_{- \infty}^\infty |A(\omega {)|^2}\frac{{{\rm d}\omega}}{{2\pi}} = \frac{{\sigma {T^2}}}{\pi}\int_{- \infty}^\infty \frac{{1 - {\rm{sinc}}(x)}}{{{x^2}}}{\rm d}x.$$
By integrating by parts twice, one finds that the $x$-integral is $\pi /2$. Hence, the spectral energy
$$\int_{- \infty}^\infty |A(\omega {)|^2}{\rm d}\omega /2\pi = \sigma {T^2}/2.$$
Equation (84) is consistent with Eq. (76). Although the oscillator response is nonstationary, it satisfies the Parseval relation (as it should do). It is interesting to note that one can convert Eq. (84) to Eq. (29) by replacing one factor of $T$ by $1/\nu$. This result reflects the fact that a detector with measurement time $T$ acts like a frequency filter with bandwidth $\nu = 1/T$.

It is sometimes helpful to replace the complicated function $[1 - {\rm{sinc}}(x)]/{x^2}$ by the simple Lorentzian $1/({a^2} + {x^2})$, which is finite for $x = 0$ and tends to $1/{x^2}$ as $x \to \infty$. This Lorentzian has a maximum of $1/{a^2}$ and an integral of $\pi /a$. By choosing $a{= 6^{1/2}}$, one reproduces the maximum correctly, but underestimates the integral. In contrast, by choosing $a = 2$, one overestimates the maximum, but reproduces the integral correctly. We prefer the second choice, because the associated spectrum satisfies the Parseval relation. The exact and approximate spectral densities, divided by $2\sigma {T^3}$, are plotted as functions of frequency in Fig. 8.

Fig. 8. Spectral energy density of an undamped oscillator. The blue curve denotes the exact spectrum (80), whereas the green and red curves denote Lorentzian approximations with $a{= 6^{1/2}}$ and 2, respectively.

If we had solved the problem in the frequency domain, using the standard derivative rule, we would have found that

$$A(\omega) = iR(\omega)/\omega .$$
It follows from this solution that the spectrum
$$\langle |A(\omega {)|^2}\rangle = \sigma T/{\omega ^2}.$$
Although these results were easy to obtain, they are wrong in several ways. First, the inverse transform of $i/\omega$ is ${\rm{sign}}(t)/2$ (Appendix D), not $H(t)$, so solution (85) is not consistent with solution (74). Second, not only is spectrum (86) infinite at the origin, but its integral is infinite, and it underestimates spectrum (80) by a factor of two at high frequencies. These are all serious failures, which were caused by a naive application of the standard derivative rule. The derivative rule is applied carefully in Appendix E.

When singularities occur in the analysis of a physical system, it is common to add some damping ($\alpha$) to the system, redo the analysis, and take the limit as $\alpha \to 0$. By following this prescription, one would obtain Eq. (85), with $\omega$ replaced by $\omega + i\alpha$, and Eq. (86), with ${\omega ^2}$ replaced by ${\omega ^2} + {\alpha ^2}$. The inverse transform of $i/(\omega + i\alpha)$ is $H(t){e^{- \alpha t}}$ (Appendix D), which tends to $H(t)$ as $\alpha \to 0$. This result is consistent with Eq. (74). However, the integral of the spectrum is $\sigma T/2\alpha$, which tends to $\infty$ as $\alpha \to 0$. For an undamped oscillator, the standard prescription does not work.

In this section and the previous ones, only the spectra $S(\omega)$ were calculated [Eqs. (10), (16), (28), (59), and (80)]. The associated frequency correlations $C(\omega ,\omega ^\prime)$ are calculated in Appendix F, for completeness.

## 6. DISCUSSION

In Sections 35, the systems (processes) considered were characterized by their means and two-time correlations (from which their variances can be deduced). These quantities provided useful information about the processes, but did not describe them completely. In this section, the probability distributions of the displacements are discussed and the processes are characterized.

Let $X$ denote a random variable (process). Then the probability that $X$ has a value between $x$ and $x + dx$ is $P(x)dx$, where $P$ is the probability distribution function. The process is said to have Gaussian statistics if the probability distribution is the Gaussian function

$$P(x) = \exp [- {(x -\mu)^2}/2{\sigma ^2}]/(2\pi {\sigma ^2}{)^{1/2}},$$
where $\mu$ and $\sigma$ are parameters. In the rest of this section, we will use the same notation for $X$ and $x$. The $n$th moment $\langle {x^n}\rangle$ is defined as the integral
$$\langle {x^n}\rangle = \int_{- \infty}^\infty {x^n}P(x){\rm d}x.$$
For distribution function (87), the first moment (mean) $\langle x\rangle =\mu$, and the second moment $\langle {x^2}\rangle = {\sigma ^2} + {\mu ^2}$, from which it follows that the variance $\langle {x^2}\rangle - {\langle x\rangle ^2} = {\sigma ^2}$.

The characteristic function

$$C(k) = \langle {e^{\textit{ikx}}}\rangle = \int_{- \infty}^\infty P(x){e^{\textit{ikx}}}{\rm d}x$$
is the Fourier transform of the distribution function. It is useful because
$$\langle {x^n}\rangle = \mathop {\lim}\limits_{k \to 0} {(- i{\partial _k})^n}C(k)$$
and differentiation is usually easier than integration. The distribution function is the inverse transform of the characteristic function, so one can use either function to characterize the process. If the distribution function is Gaussian, then the characteristic function
$$C(k) = \exp (ik\mu - {\sigma ^2}{k^2}/2)$$
is also Gaussian. Notice that the coefficient of ${k^2}$ in Eq. (91) is proportional to the variance. If $x$ is measured relative to its mean, the $\mu$-terms in Eqs. (87) and (91) are absent.

Now let $X = [{x_1},{x_2} \ldots {x_n}{]^t}$ denote multiple random variables (displacement and velocity, for example) and $K = [{k_1},{k_2} \ldots {k_n}]$ denote multiple transform variables. Then this multiple-variable process is said to be Gaussian, with zero means, if its characteristic function

$$C(K) = \exp (- {K^t}MK/2),$$
where $M = [{m_{\textit{ij}}}]$ is a matrix. The moments
$$\langle x_1^{{n_1}}\rangle \def\LDeqtab{}= \mathop {\lim}\limits_{{k_1} \to 0} {(- i{\partial _{{k_1}}})^{{n_1}}}C(K),$$
$$\langle x_1^{{n_1}}x_2^{{n_2}}\rangle \def\LDeqtab{}= \mathop {\lim}\limits_{{k_1} \to 0} \mathop {\lim}\limits_{{k_2} \to 0} {(- i{\partial _{{k_1}}})^{{n_1}}}{(- i{\partial _{{k_2}}})^{{n_2}}}C(K),$$
and higher moments are defined by similar formulas. By applying formulas (93) and (94), one finds that the variances $\langle x_i^2\rangle = {m_{\textit{ii}}}$ and the correlations $\langle {x_i}{x_j}\rangle = {m_{\textit{ij}}}$. For this reason, $M$ is called the covariance matrix. By inverse transforming the characteristic function, one obtains the distribution function
$$P(X) = \frac{{\exp (- {X^t}{M^{- 1}}X/2)}}{{{{[(2\pi {)^n}\det (M)]}^{1/2}}}}.$$
Notice that the coefficient matrix in formula Eq. (95) is the inverse of the covariance matrix (and vice versa). Formula Eq. (95) is standard. A convenient derivation of it is provided in Appendix A of [15].

For many stochastic processes, it is reasonable to assume that the noise source has Gaussian statistics. To be precise, define the (Wiener) increment, or impulse, $W(\tau) = \int _0^\tau R(s){\rm d}s$, where $R(s)$ is the source function and $\tau$ is a short time interval. Then, for white noise, the impulses have a Gaussian distribution function with mean $\langle W(\tau)\rangle = 0$ and variance $\langle |W(\tau {)|^2}\rangle = \sigma \tau$, where $\sigma$ is the source strength. This assumption is common for two complementary reasons. First, suppose that $\hat a$ is the amplitude operator of a quantum mode, in which case $\langle {\hat a^\dagger}\hat a\rangle$ is the expected number of quanta. Then, if the mode is in a coherent state, the expected quadratures have Gaussian probability distributions with variances of 1/2 [15,16]. If the mode experiences linear amplification or attenuation, the quadrature means and variances change, but their distributions remain Gaussian. One can reproduce this behavior semi-classically by adding noise terms with Gaussian statistics to the classical amplitude equation. (Gain and loss both degrade the signal-to-noise ratio.) Second, if a classical diplacement is driven by many independent impulses, the central limit theorem [5,8] ensures that the displacement statistics are Gaussian, even if the impulse statistics are not.

A time-dependent random process is said to be strictly stationary if the $n$th moment $\langle x({t_1})x({t_2}) \ldots x({t_n})\rangle$ is independent of the origin of time. In other words, one can replace ${t_1}$, ${t_2}$, $\ldots ;{t_n}$ by ${t_1} - \tau$, ${t_2} - \tau$, $\ldots \;{t_n} - \tau$ without changing the moment. A process is said to be weakly (or wide-sense) stationary if the first two moments $\langle x({t_1})\rangle$ and $\langle x({t_1})x({t_2})\rangle$ are independent of the time origin. For such a process, the mean is constant, and the two-time correlation is a function of only the time difference ${t_2} - {t_1}$. Because the mean is constant, one can simplify the analysis of the process by measuring $x(t)$ relative to the mean.

A time-dependent process is said to have Gaussian statistics if its $n$-sample characteristic function has the form (92). If such a process is weakly stationary, it is also strictly stationary, because the correlation condition applies to arbitrary ${t_1}$ and ${t_2}$ (or ${t_i}$ and ${t_j}$). Hence, in the covariance matrix, each entry ${m_{\textit{ij}}}$ is a function of ${t_j} - {t_i}$, which is independent of the time origin. Consequently, so also is the $n$th moment, which is determined completely by the covariance matrix.

In Sections 35, the oscillator displacements (22), (39), and (74) are sums of Gaussian impulses multipled by deterministic response (Green) functions, so the displacements also have Gaussian statistics, which are specified completely by their means (of zero), and correlations (27), (56), and (75), from which the variances follow. In steady state, strongly and weakly damped oscillations are weakly stationary processes. Hence, as explained above, they are also stationary. In contrast, undamped oscillation is a nonstationary process, because the variance of the displacement increases linearly with time. Nonetheless, undamped oscillation has stationary increments: The displacement difference

$$\delta X(t,\tau) = X(t + \tau) - X(t) = \int_t^{t + \tau} R(s){\rm d}s.$$
Because the source term in Eq. (96) is stationary, the integral (increment) does not depend on the time origin $t$, so the displacement difference is a function of $\tau$ only: $\delta X(\tau) = \int _0^\tau R(s){\rm d}s$. Thus, the analyses of Sections 35 provide examples of the three main classes of stationarity, about which much is known [7].

## 7. SUMMARY

In this tutorial, three examples of stochastic systems were discussed: a strongly damped oscillator, a weakly damped oscillator and an undamped oscillator (integrator) driven by noise. The dynamics of these systems are characterized by the temporal correlation functions and spectral densities of their displacements [$\langle {A^*}(t)A(t^\prime)\rangle$ and $\langle |A(\omega {)|^2}\rangle$, respectively], which were determined and discussed.

After transients associated with the initial conditions have died out, strongly and weakly damped oscillators attain stochastic steady states, in which the temporal correlations (27) and (56) depend on the time difference $|t - t^\prime |$, rather than $t$ and $t^\prime$ separately. The correlations are decreasing functions of the time difference, because the effects of noise impulses only influence the displacements for times of the order of the damping time. (These systems have finite memory.) Because these systems (processes) are weakly stationary, the Wiener–Khinchin theorem applies, and their spectral energy densities (28) and (59) are the Fourier transforms of their temporal correlations with respect to the time difference, multiplied by the measurement time $T$. Although the idealized noise model on which these results are based (white noise) has infinite energy, the spectral energy densities and integrated spectral energies are finite, because damping filters the applied noise.

In contrast, an undamped oscillator never reaches a steady state. The oscillator power increases monotonically with time, because the oscillator integrates the effects of the applied noise impulses. (This system has infinite memory.) The correlation (75) is a function of $\min (t,t^\prime)$, because the noise applied between $t$ and $t^\prime$ has no effect on the correlation of the earlier and later displacements. The spectrum (80) has a sharp peak at zero frequency, but the spectral energy density and integrated spectral energy are finite if the measurement time is finite. Because this system is nonstationary, the Wiener–Khinchin theorem does not apply.

For each system, the dynamics were analyzed in the time domain (by solving differential equations) and in the frequency domain (by solving algebraic equations). The frequency-domain analyses were based on the standard derivative rule that the Fourier transform of ${d_t}A$ is the transform of $A$ multiplied by ${-}i\omega$. For damped oscillators, the frequency-domain analyses miss the transient displacements, but predict the time-asymptotic steady states accurately. (The relative errors in the spectral densities and energies are of order $1/\nu T$, where $\nu$ is the damping rate.) However, for an undamped oscillator, the frequency-domain analysis produces results that are incorrect and singular. These failures are caused by the omission of the final value $A(T)$ in the derivative rule. When one models driven systems, in which the displacements do not tend to zero as time tends to infinity, one should use the derivative rule with caution.

The mathematical methods required to study randomly driven oscillators and the physical insights one gains from such a study can be applied to a variety of problems, in fields such as biology, chemistry, finance, and optimization [1720]. In a subsequent tutorial [21], we will use the methods and insights discussed herein to study laser linewidths (phase fluctuations). This tutorial was restricted to linear oscillators subject to additive noise. A review of linear and nonlinear oscillators, subject to additive and multiplicative noise, can be found in [22].

## APPENDIX A: FOURIER TRANSFORMS

Let ${A_2}(t)$ be an integrable two-sided function of time (defined for positive and negative times). Then the forward and backward Fourier transforms are defined by the equations

$${A_2}(\omega) \def\LDeqtab{}= \int_{- \infty}^\infty {A_2}(t){e^{i\omega t}}{\rm d}t,$$
$${A_2}(t) \def\LDeqtab{}= \int_{- \infty}^\infty {A_2}(\omega){e^{- i\omega t}}{\rm d}\omega /2\pi .$$
In both integrals, the independent variables are real and range from ${-}\infty$ to $\infty$. Now let ${B_2}(t)$ be another integrable two-sided function. Then the mathematical cross-correlation of ${A_2}$ and ${B_2}$ is defined by the integral
$${C_2}(\tau) = \int_{- \infty}^\infty A_2^*(t){B_2}(t + \tau){\rm d}t,$$
where the time difference $\tau$ is arbitrary. (The physical correlation involves integration from zero to $T$ and division by $T$, so that the correlation has the same units as the product of the amplitudes.) Let ${C_2}(\omega)$ be the Fourier transform of the auto-correlation with respect to the time difference. Then the correlation theorem states that
$${C_2}(\omega) = A_2^*(\omega){B_2}(\omega),$$
where $A_2^*(\omega)$ is an abbreviation for ${[{A_2}(\omega)]^*}$. For the special case in which the functions are identical, ${C_2}(\omega) = |{A_2}(\omega {)|^2}$. The Fourier transform of the temporal auto-correlation is the spectral energy density. Likewise, the inverse transform of the spectrum is the auto-correlation. For the special case in which $\tau = 0$, the correlation theorem reduces to the Parseval equation
$$\int_{- \infty}^\infty |{A_2}(t{)|^2}{\rm d}t = \int_{- \infty}^\infty |{A_2}(\omega {)|^2}{\rm d}\omega /2\pi ,$$
which states that the time- and frequency-domain definitions of energy are consistent.

The convolution is defined by the integral

$${D_2}(\tau) = \int_{- \infty}^\infty {A_2}(t){B_2}(\tau - t){\rm d}t,$$
where $\tau$ is arbitrary. In contrast to the correlation integral, ${A_2}$ is not conjugated, and ${B_2}$ is a function of $\tau - t$, not $\tau + t$. The convolution theorem states that
$${D_2}(\omega) = {A_2}(\omega){B_2}(\omega).$$
It is used in analyses of linear systems, in which context ${A_2}$ is the driving term, and ${B_2}$ is the response (Green) function.

This tutorial involves the solutions of initial-value problems, which are defined for nonnegative times. If these solutions were integrable, one could derive one-sided versions of the correlation and convolution theorems from the two-sided versions by replacing ${A_2}(t)$ with ${A_1}(t) = {A_2}(t)H(t)$, where $H$ is the Heaviside step function. However, in the problems considered, the dependent variables are driven by noise, which is always present, so the solutions are not integrable. One can prevent singularities in the analyses of these problems by limiting the integration (measurement) time to the large, but finite, value $T$, and one can derive finite-time versions of the correlation and convolution theorems from the one-sided versions by replacing ${A_1}(t)$ with ${A_T}(t) = {A_1}(t)H(T - t) = {A_2}(t)H(t)H(T - t)$. In this appendix, finite-time results are derived directly from the two-sided results stated above. One can deduce the corresponding one-sided results by letting $T \to \infty$ and assuming integrability.

The finite-time amplitudes ${A_T}(t)$ and ${B_T}(t)$ are only nonzero for $0 \le t \le T$, so the Fourier transform

$${A_T}(\omega) = \int_0^T {A_T}(t){e^{i\omega t}}{\rm d}t$$
and the correlation
$${C_T}(\tau) = \int_0^T A_T^*(t){B_T}(t + \tau){\rm d}t.$$
On the right side of Eq. (A9), the second function is nonzero for ${-}t \le \tau \le T - t$. Because the extreme values of $t$ are zero and $T$, the extreme values of $\tau$ are ${-}T$ and $T$. Hence, the Fourier transform
$${C_T}(\omega) = \int_{- T}^T {C_T}(\tau){e^{i\omega \tau}}{\rm d}\tau .$$
The integration region for the correlation theorem is illustrated in Fig. 9(a).

Fig. 9. Integration regions in the $t - \tau$ plane for the correlation theorem (left) and the convolution theorem (right). The horizontal axis represents $t$ and the vertical axis represents $\tau$.

By combining Eqs. (A9) and (A10), and reversing the order of integration, one finds that

$$\begin{split}{C_T}(\omega) &= \int_0^T \int_{- T}^T A_T^*(t){e^{- i\omega t}}{B_T}(t + \tau){e^{i\omega (t + \tau)}}{\rm d}\tau {\rm d}t\\ &= \int_0^T A_T^*(t){e^{- i\omega t}}\int_{- t}^{T - t} {B_T}(t + \tau){e^{i\omega (t + \tau)}}{\rm d}\tau {\rm d}t\\& = \int_0^T A_T^*(t){e^{- i\omega t}}dt\int_0^T {B_T}(t^\prime){e^{i\omega t^\prime}}{\rm d}t^\prime \\ &= A_T^*(\omega){B_T}(\omega).\end{split}$$
The second line of Eq. (A11) was obtained from the first by using the step functions that are implicit in ${B_T}$. It is also instructive to consider the inverse transform
$$\begin{split}{C_T}(\tau) &= \int_{- \infty}^\infty \int_0^T \int_0^T A_T^*(t){B_T}(t^\prime){e^{i\omega (t^\prime - t - \tau)}}{\rm d}t^\prime {\rm d}t{\rm d}\omega /2\pi \\ &= \int_0^T \int_0^T A_T^*(t){B_T}(t^\prime)\delta (t^\prime - t - \tau){\rm d}t^\prime {\rm d}t\\& = \int_0^T A_T^*(t){B_T}(t + \tau){\rm d}t.\end{split}$$
The second line of Eq. (A12) shows that ${C_T}(\tau)$ is defined for ${-}T \le \tau \le T$, as stated above.

The finite-time convolution

$${D_T}(\tau) = \int_0^T {A_T}(t){B_T}(\tau - t){\rm d}t.$$
On the right side of Eq. (A13), the second function is nonzero for $t \le \tau \le T + t$, so the formal upper limit $T$ can be replaced by the effective limit $\tau$. Because the extreme values of $t$ are zero and $T$, the extreme values of $\tau$ are zero and $2T$. Hence, the Fourier transform
$${D_T}(\omega) = \int_0^{2T} {D_T}(\tau){e^{i\omega \tau}}{\rm d}\tau .$$
The integration region for the convolution theorem is illustrated in Fig. 9(b).

By combining Eqs. (A13) and (A14), and reversing the order of integration, one finds that

$$\begin{split}{D_T}(\omega)& = \int_0^T \int_0^{2T} {A_T}(t){e^{i\omega t}}{B_T}(\tau - t){e^{i\omega (\tau - t)}}{\rm d}\tau {\rm d}t,\\ &= \int_0^T {A_T}(t){e^{i\omega t}}\int_t^{T + t} {B_T}(\tau - t){e^{i\omega (\tau - t)}}{\rm d}\tau {\rm d}t\\ &= \int_0^T {A_T}(t){e^{i\omega t}}{\rm d}t\int_0^T {B_T}(t^\prime){e^{i\omega t^\prime}}{\rm d}t^\prime \\& = {A_T}(\omega){B_T}(\omega).\end{split}$$
The second line of Eq. (A15) was obtained from the first by using the step functions that are implicit in ${B_T}$. One also finds that the inverse transform
$$\begin{split}{D_T}(\tau) &= \int_{- \infty}^\infty \int_0^T \int_0^T {A_T}(t){B_T}(t^\prime){e^{i\omega (t^\prime + t - \tau)}}{\rm d}t^\prime {\rm d}t{\rm d}\omega /2\pi \\[-3pt]& = \int_0^T \int_0^T {A_T}(t){B_T}(t^\prime)\delta (t^\prime + t - \tau){\rm d}t^\prime {\rm d}t\\[-3pt]& = \int_0^T {A_T}(t){B_T}(\tau - t){\rm d}t.\end{split}$$
The second line of Eq. (A16) shows that ${D_T}(\tau)$ is defined for $0 \le \tau \le 2T$, as stated above. Thus, the correlation and convolution theorems for finite integration times are similar to those for infinite times. The only differences are the limits of integration in definitions (A10) and (A14).

The preceding discussion of the convolution theorem was based on the assumption that the functions $A$ and $B$ are both defined for arguments between zero and $T$. Solutions (20), (39), and (74) were written in the Green-function form

$$A(t) = \int_0^t G(t - s)R(s){\rm d}s,$$
where $0 \le s \le t$ (causality), and $t \le T$ (finite time). The transformed amplitude
$$\begin{split}{A_T}(\omega) &= \int_0^T \int_0^t G(t - s)R(s){e^{i\omega t}}{\rm d}s{\rm d}t\\[-3pt]& = \int_0^T \int_s^T G(t - s)R(s){e^{i\omega (t - s + s)}}{\rm d}t{\rm d}s\\[-3pt]& = \int_0^T R(s){e^{i\omega s}}\int_0^{T - s} G(\tau){e^{i\omega \tau}}{\rm d}\tau {\rm d}s,\end{split}$$
where $\tau = t - s$. The right side of Eq. (A18) is not a product of two independent finite-time transforms.

For a strongly damped oscillator, the transformed Green function ${G_T}(\omega) = 1/(\nu - i\omega)$, and for a weakly damped oscillator, ${G_T}(\omega) = \epsilon /(\omega _0^2 - 2i{\nu _0}\omega - {\omega ^2})$, where exponentially small terms were omitted. By squaring these functions and multiplying by the noise spectral density $\sigma T$, one obtains the largest contribution to spectrum (28) and spectrum (59), respectively. Thus, for damped oscillators, formula (A15) produces spectra that are asymptotically correct.

For an undamped oscillator, $G(t) = H(t)$, so its transform

$${G_T}(\omega) = ({e^{i\omega T}} - 1)/i\omega ,$$
and the associated spectrum
$$\langle |{A_T}(\omega {)|^2}\rangle = 2\sigma T[1 - \cos (\omega T)]/{\omega ^2},$$
where nothing was omitted. Formula (A20) is much better than the naive formula $\sigma T/{\omega ^2}$. It has a maximum of $\sigma {T^3}$ at $\omega = 0$, an integral of $2\pi \sigma {T^2}$, and is proportional to $2\sigma T/{\omega ^2}$ (with the correct factor of two). However, it differs from the exact formula (80), which has a maximum of $\sigma {T^3}/3$ and an integral of $\pi \sigma {T^2}$. Thus, for the systems of interest to us, the convolution theorem does not apply as written, and for an undamped oscillator, the relative error is of order one, which does not tend to zero as $T$ tends to infinity. It is for this reason that we did not use the standard formula $A(\omega) = G(\omega)R(\omega)$ to calculate spectra.

## APPENDIX B: DOUBLE FOURIER TRANSFORMS

Equation (6) shows that one can always determine the spectral energy density of a stochastic process from the two-time correlation. However, one cannot always determine the temporal correlation from the spectrum, because there is only one frequency variable ($\omega)$, but two time variables ($t$ and $t^\prime$) are required. (If one already knows that the process is stationary, then one can apply the Wiener–Khinchin theorem, which is stated in Appendix C.) In this appendix, the relations between the time- and frequency-domain characterizations of stochastic processes are discussed.

Define the two-time correlation $C(t,t^\prime) = \langle {A^*}(t)A(t^\prime)\rangle$ and the two-frequency correlation $C(\omega ,\omega ^\prime) = \langle {A^*}(\omega)A(\omega ^\prime)\rangle$, and suppose that both functions are integrable. Then they are linked by the double Fourier transforms

$$C(\omega ,\omega ^\prime) \def\LDeqtab{}= \int_{- \infty}^\infty \int_{- \infty}^\infty C(t,t^\prime){e^{i\omega ^\prime t^\prime - i\omega t}}{\rm d}t{\rm d}t^\prime ,$$
$$C(t,t^\prime) \def\LDeqtab{}= \int_{- \infty}^\infty \int_{- \infty}^\infty C(\omega ,\omega ^\prime){e^{i\omega t - i\omega ^\prime t^\prime}}{\rm d}\omega {\rm d}\omega ^\prime /(2\pi {)^2}.$$
For reference, recall that
$$\int_{- \infty}^\infty {e^{\pm i\omega t}}{\rm d}\omega /2\pi = \delta (t),$$
from which it follows that
$$\int_{- \infty}^\infty {e^{\pm i\omega t}}{\rm d}t = 2\pi \delta (\omega).$$

Now consider some examples. Suppose that the correlation,

$$C(t,t^\prime) = S(t)\delta (t - t^\prime),$$
is $\delta$-correlated in time. Then its transform,
$$C(\omega ,\omega ^\prime) = S(\omega ^\prime - \omega),$$
is a function of the frequency difference $\omega ^\prime - \omega$. For the special case in which $S(t) = \sigma$, where $\sigma$ is constant, $C(\omega ,\omega ^\prime) = 2\pi \sigma \delta (\omega - \omega ^\prime)$. Conversely, if the correlation,
$$C(\omega ,\omega ^\prime) = 2\pi S(\omega)\delta (\omega - \omega ^\prime),$$
is $\delta$-correlated in frequency, then its inverse transform,
$$C(t,t^\prime) = S(t^\prime - t),$$
is a function of the time difference $t^\prime - t$. For the special case in which $S(\omega) = \sigma$, then $C(t,t^\prime) = \sigma \delta (t - t^\prime)$, as stated above. If we had used frequency $f$ rather than angular frequency $\omega$, the factors of $2\pi$ in the preceding results would have been absent, because $a\delta (ax) = \delta (x)$. The bidirectional relations described by Eqs. (B5)–(B8) are well known. However, they are also problematic. In the first relation, the power $\langle |A(t{)|^2}\rangle = S(t)\delta (0)$ is undefined, whereas in the second relation, the spectrum $\langle |A(\omega {)|^2}\rangle = 2\pi S(\omega)\delta (0)$ is undefined.

Because the displacements and correlation functions associated with driven processes are not integrable, one replaces infinite-time integrals by finite ones. In the one-sided formalism, the amplitude ${A_T}(t)$ is nonzero for $0 \le t \le T$. The frequency correlation

$${C_T}(\omega ,\omega ^\prime) = \int_0^T \int_0^T {C_T}(t,t^\prime){e^{i\omega ^\prime t^\prime - i\omega t}}{\rm d}t{\rm d}t^\prime ,$$
and the time-integral identity
$$\int_0^T {e^{\pm i\omega t}}{\rm d}t = ({e^{i\omega T}} - 1)/i\omega = 2\pi {\delta _T}(\omega),$$
where ${\delta _T}$ is the complex $\delta$-function that is discussed in Appendix D (see Fig. 11). Notice that $2\pi {\delta _T}(0) = T$.

Now reconsider the previous examples with finite time integrals. If the temporal correlation ${C_T}(t,t^\prime) = {S_T}(t)\delta (t - t^\prime)$, then its transform

$${C_T}(\omega ,\omega ^\prime) = \int_0^T {S_T}(t){e^{i(\omega ^\prime - \omega)t}}{\rm d}t = {S_T}(\omega ^\prime - \omega).$$
For the special case in which ${S_T}(t) = \sigma$, ${C_T}(\omega ,\omega ^\prime) = 2\pi \sigma {\delta _T}(\omega - \omega ^\prime)$, where the effective $\delta$-function was defined in Eq. (B10). The associated spectrum $\langle |{A_T}(\omega {)|^2}\rangle = \sigma T$. If the correlation ${C_T}(t,t^\prime) = {S_T}(t^\prime - t)$, where ${-}T \le t^\prime - t \le T$, then its transform
$$\begin{split}{C_T}(\omega ,\omega ^\prime) &= \int_0^T \int_0^T {S_T}(t^\prime - t){e^{i\omega ^\prime (t^\prime - t)}}{e^{i(\omega ^\prime - \omega)t}}{\rm d}t{\rm d}t^\prime \\[-3pt] &= \int_0^T {e^{i(\omega ^\prime - \omega)t}}\int_{- t}^{T - t} {S_T}(\tau){e^{i\omega ^\prime \tau}}{\rm d}\tau {\rm d}t\\[-3pt] &= \int_{- T}^T {S_T}(\tau){e^{i\omega ^\prime \tau}}\int_{\max (0, - \tau)}^{\min (T - \tau ,T)} {e^{i(\omega ^\prime - \omega)t}}{\rm d}t{\rm d}\tau \\[-3pt]& \sim 2\pi {S_T}(\omega ^\prime){\delta _T}(\omega ^\prime - \omega).\end{split}$$
The integration region in the $t - \tau$ plane is illustrated in Fig. 9(a). A related integral is considered carefully in Appendix C (see Fig. 10). For now, it is sufficient to know that if $T$ is large and ${S_T}(\tau)$ is a decreasing function of $|\tau |$, then Eq. (B12) is asymptotically correct. Equations (B10)–(B12) are equivalent to Eqs. (B4)–(B8), including the factors of $2\pi$. The only difference between the two sets of equations is that the standard $\delta$-function, which is singular, is replaced by the effective $\delta$-function, which is nonsingular. Using finite time intervals removes the singularities in the spectra.

## APPENDIX C: WIENER–KHINCHIN THEOREM

In Section 2, the spectral power densities of white and colored noise were shown to be the Fourier transforms of their temporal correlation functions with respect to the time difference. In Sections 3 and 4, this relation, which is called the Wiener–Khinchin theorem, was demonstrated for strongly and weakly damped oscillators that have attained stochastic steady states. This theorem is similar to the auto-correlation theorem of Appendix A. The difference is that the correlation theorem involves the time integral of an integrable function, whereas the Wiener–Khinchin theorem involves the ensemble average of a nonintegrable function.

The inverse transform of the spectral energy density is

$$\begin{split}&\int_{- \infty}^\infty \langle |{A_T}(\omega {)|^2}\rangle {e^{- i\omega \tau}}{\rm d}\omega /2\pi\\&\quad = \int_{- \infty}^\infty \int_0^T \int_0^T \langle A_T^*(t){A_T}(t^\prime)\rangle {e^{i\omega (t^\prime - t - \tau)}}{\rm d}t{\rm d}t^\prime {\rm d}\omega /2\pi \\&\quad = \int_0^T \int_0^T \langle A_T^*(t){A_T}(t^\prime)\rangle \delta (t^\prime - t - \tau){\rm d}t{\rm d}t^\prime \\ &\quad= \int_{\max (0, - \tau)}^{\min (T - \tau ,T)} \langle A_T^*(t){A_T}(t + \tau)\rangle {\rm d}t,\end{split}$$
where ${-}T \le \tau \le T$. For a damped system, which attains a weakly stationary state, the correlation depends on only the time difference [$C(t,t + \tau) = C(0,\tau) = C(\tau)$]. For such a system, the correlation is constant along horizontal lines in Fig. 10. It follows from the Cauchy–Schwartz inequality that $|C(\tau)| \lt C(0)$, so the largest contribution to the integral in Eq. (C1) comes from the line $\tau = 0$. For a damped system, the correlation is a decreasing function of $|\tau |$, with decay rate $\nu$. Only the lines with $|\tau | \lt 1/\nu$ contribute significantly to the integral. The contributing region is almost a rectangle of width $T$ and height $2/\nu$, so
$$S(\omega) \sim T\int_{- 1/\nu}^{1/\nu} C(\tau){e^{i\omega \tau}}{\rm d}\tau .$$
The integration error associated with the missing triangles is negative, and is of relative order $1/\nu T$. Because regions with $|\tau | \gt 1/\nu$ do not contribute significantly to the integral, one can extend the limits of integration and rewrite Eq. (C2) in the standard form
$$S(\omega)/T \sim \int_{- T}^T C(\tau){e^{i\omega \tau}}{\rm d}\tau .$$
Thus, the spectral power density is asymptotic to the Fourier transform of the temporal correlation function with respect to the time difference.

Fig. 10. Integration region in the $t - \tau$ plane for the Wiener–Khinchin theorem. The horizontal axis represents $t$ and the vertical axis represents $\tau$. The shaded region contributes to the integral.

For white noise, only the line $\tau = 0$ contributes to the integral in Eq. (C1), so Eq. (C3) is exact. For colored noise, Eq. (C3) is asymptotic, as written. For a damped system that starts from rest, the values of the correlation $C(t,t + \tau)$ are smaller than the steady-state value $C(\tau)$ for $0 \lt t \lt 1/\nu$ and comparable to it for $1/\nu \lt t \lt T$. The lengths of the left line segments (one line for each value of $\tau$) are much shorter than those of the right segments, so Eq. (C3) also applies when the transient response is included.

## APPENDIX D: SPECIFIC FOURIER TRANSFORMS

In this appendix, the forward and backward Fourier transforms of some functions that appear in the main text are calculated. For a strongly damped oscillator (Section 3), the time-domain Green function

$$G(t) = H(t){e^{- \nu t}},$$
where the Heaviside step function
$$H(t) = \left\{{\begin{array}{*{20}{c}}0&{{\rm{if}}}&{t \lt 0,}\\{1/2}&{{\rm{if}}}&{t = 0,}\\1&{{\rm{if}}}&{t \gt 0.}\end{array}} \right.$$
Notice that this Green function is integrable. It follows from Eqs. (A1) and (D1) that the frequency-domain Green function
$$G(\omega) = 1/(\nu - i\omega) = (\nu + i\omega)/({\nu ^2} + {\omega ^2}).$$

If one uses contour integration to evaluate the inversion integral (A2), one needs to complete the contour in the upper (lower) half-plane for $t \lt 0$ ($t \gt 0$). The Green function in Eq. (D2) is analytic in the upper half-plane, but has a pole in the lower half-plane at $\omega = - i\nu$. Hence, for $t \lt 0$, $G(t) = 0$. For $t = 0$, the imaginary part of the inversion integral is zero, whereas the real part is 1/2 [because the integral of $1/({\nu ^2} + {\omega ^2})$ is $\pi /\nu$]. For $t \gt 0$, the inversion integral is ${-}2\pi i/2\pi$ multiplied by the residue $i{e^{- \nu t}}$, which is ${e^{- \nu t}}$. Hence, the inverse transform of the frequency-domain Green function equals the time-domain Green function (as it should do).

For an undamped oscillator (Section 5), the time-domain Green function

$$G(t) = H(t).$$
Unfortunately, this Green function is not integrable. The related function $H(t){e^{- \alpha t}}$ is integrable when $\alpha$ is positive, so one can try to fix the non-integrability problem by using the limit function
$${G_\alpha}(t) = \mathop {\lim}\limits_{\alpha \to 0} H(t){e^{- \alpha t}}.$$
It follows from Eqs. (D3) and (D5) that
$${G_\alpha}(\omega) = \mathop {\lim}\limits_{\alpha \to 0} (\alpha + i\omega)/({\alpha ^2} + {\omega ^2}).$$
How should one interpret this limit function? The real part has a peak at $\omega = 0$, of height $1/\alpha$ and width $\alpha$, decreases like $1/{\omega ^2}$ as $\omega \to \infty$ and its integral is $\pi$. The imaginary part is an odd function of $\omega$, decreases like $1/\omega$ as $\omega \to \infty$ and its integral is zero. Hence, as $\alpha \to 0$,
$${G_\alpha}(\omega) \to \pi \delta (\omega) + i/\omega .$$
Formula (D7) is standard [23]. It is illustrated in Fig. 11.

The backward transform of the first term on the right side of Eq. (D7) is 1/2, for all times. The inversion integral for the second term is

$$\int_{- \infty}^\infty [i\cos (\omega t) + \sin (\omega t)]{\rm d}\omega /2\pi \omega .$$
The real part of this integral is ${\rm{sign}}(t)/2$ [because the integral of ${\rm{sinc}}(x)$ is $\pi$ and sine is an odd function of $t$], whereas the imaginary part has a principal value of zero. Thus, the inverse transform of the frequency-domain limit function equals the time-domain limit function $H(t)$. The limit procedure used to obtain Eq. (D7) is sensible.

Another important function (which was just mentioned) is the sign function

$$\Sigma (t) = \left\{{\begin{array}{*{20}{c}}{- 1}&{{\rm{if}}}&{t \lt 0,}\\0&{{\rm{if}}}&{t = 0,}\\1&{{\rm{if}}}&{t \gt 0.}\end{array}} \right.$$
Notice that $\Sigma (t) = H(t) - H(- t)$ is the odd extension of the step function (and 1 is the even extension). This function is also not integrable, so one defines the limit function
$${\Sigma _\alpha}(t) = \mathop {\lim}\limits_{\alpha \to 0} \Sigma (t){e^{- \alpha |t|}}.$$
It follows from Eqs. (D3) and (D10) that
$$\begin{split}{\Sigma _\alpha}(\omega) &= \mathop {\lim}\limits_{\alpha \to 0} [1/(\alpha - i\omega) - 1/(\alpha + i\omega)]\\& = \mathop {\lim}\limits_{\alpha \to 0} 2i\omega /({\alpha ^2} + {\omega ^2}).\end{split}$$
The function on the right side of Eq. (D11) is zero at $\omega = 0$, has large extremes at $\omega = \pm \alpha$, and decreases like $2i/\omega$ as $\alpha \to 0$. Hence, as $\alpha \to 0$,
$${\Sigma _\alpha}(\omega) \to 2i/\omega .$$
Formula (D12) is also standard [23]. It is illustrated (divided by two) in Fig. 11(b), and its inverse transform was discussed after Eq. (D8).

Fig. 11. Real (top) and imaginary (bottom) parts of the Fourier transform of the step function. The dashed curves represent formula (D7), whereas the solid curves represent formula (D13).

One can also try to fix the non-integrability problems by using finite-time transforms, in which context the truncated step function ${H_T}(t) = H(t)H(T - t)$. Its transform

$$\begin{split}{H_T}(\omega) &= ({e^{i\omega T}} - 1)/i\omega \\ &= \sin (\omega T)/\omega + i[1 - \cos (\omega T)]/\omega .\end{split}$$
This transform (divided by $T$) is illustrated in Fig. 11. The real part of the transform has a peak at $\omega = 0$, of height $T$ and width $\pi /T$, and its integral is $\pi$. The imaginary part is an odd function of $\omega$, its local average decreases like $1/\omega$ as $\omega \to \infty$, and its integral is zero. Hence, for large values of $T$,
$${H_T}(\omega) \approx \pi \delta (\omega) + i/\omega .$$
Formula Eq. (D14) is consistent with formula Eq. (D7).

Likewise, the truncated sign function ${\Sigma _T}(t) = \Sigma (t)H(t + T)H(T - t)$, and its transform

$$\begin{split}{\Sigma _T}(\omega) &= ({e^{i\omega T}} - 1)/i\omega - ({e^{- i\omega T}} - 1)/(- i\omega)\\ &= 2i[1 - \cos (\omega T)]/\omega .\end{split}$$
This transform (divided by $2T$) is illustrated in Fig. 11(b). The function on the right side of Eq. (D15) was discussed after Eq. (D13). For large values of $T$,
$${\Sigma _T}(\omega) \approx 2i/\omega .$$
Formula (D16) is consistent with formula (D12).

Now consider the inverse transforms of formulas (D13) and (D15). If the time-integration domain is extended slightly to the left of ${-}T$ and slightly to the right of $T$, the end values of the step function are zero, in which case one can apply the rule that the transform of ${d_t}A(t)$ is the transform of $A(t)$ multiplied by ${-}i\omega$. Conversely, the inverse transform of $iA(\omega)/\omega$ is the integral $\int _{- T}^tA(s)ds$. Hence, the inverse transform of ${H_T}(\omega)$ is

$$\begin{split}&\int_{- T}^t \int_{- \infty}^\infty (1 - {e^{i\omega T}}){e^{- i\omega s}}{\rm d}\omega {\rm d}s/2\pi \\ &\quad= \int_{- T}^t [\delta (s) - \delta (s - T)]ds = {H_T}(t).\end{split}$$
Likewise, the inverse transform of ${\Sigma _T}(\omega)$ is
$$\begin{split}&\int_{- T}^t \int_{- \infty}^\infty (2 - {e^{i\omega T}} - {e^{- i\omega T}}){e^{- i\omega s}}{\rm d}\omega {\rm d}s/2\pi \\&\quad= \int_{- T}^t [- \delta (s + T) + 2\delta (s) - \delta (s - T)]{\rm d}s = {\Sigma _T}(t).\end{split}$$
Thus, in the finite-time approach, the inverse transforms of the frequency-domain functions equal the original time-domain functions (as they should do).

Suppose that $F(t)$ describes a causal response, which equals zero for $t \lt 0$. Then its Fourier transform $F(\omega)$ is analytic in the upper half-plane and tends to zero as $|\omega | \to \infty$. By integrating along a contour that consists of the real $\omega ^\prime$ axis, a small semicircle around the point $\omega ^\prime = \omega$, and a large semicircle in the upper half-plane, one obtains the Kramers–Kronig equation

$$F(\omega) = \frac{i}{\pi}P\int_{- \infty}^\infty \frac{{F(\omega ^\prime)}}{{\omega - \omega ^\prime}}{\rm d}\omega ^\prime ,$$
where $P$ denotes a principal value [24,25]. The integral on the right side of Eq. (D19), without the factor of $i$, is called a Hilbert transform. By splitting Eq. (D19) into real and imaginary parts, one finds that the real and imaginary parts of $F(\omega)$ are Hilbert transforms of each other. They contain the same information, so one cannot exist without the other. If $F(t)$ is real, then the real part of $F(\omega)$ is an even function of frequency (the cosine transform), whereas the imaginary part is an odd function of frequency (the sine transform). Transforms (D3), (D7), (D12), (D13), and (D15) all have these properties.

Let $u$ denote the upper limit $\infty$ or $T$. Then it follows from Eqs. (2) and (3) that

$$\begin{split}F(\omega) = \int_0^u \int_{- \infty}^\infty F(\omega ^\prime){e^{- i\omega ^\prime t}}({\rm d}\omega ^\prime /2\pi){e^{i\omega t}}{\rm d}t\\ = \int_{- \infty}^\infty F(\omega ^\prime)\int_0^u {e^{i(\omega - \omega ^\prime)t}}(dt/2\pi){\rm d}\omega ^\prime .\end{split}$$
The left and right sides of Eq. (D20) are consistent if and only if
$$\int_0^u {e^{i(\omega - \omega ^\prime)t}}{\rm d}t/2\pi = {\delta _u}(\omega - \omega ^\prime).$$
The left side of Eq. (D21) is the Fourier transform of a step function divided by $2\pi$. Unlike common representations of the $\delta$-function, which are real, the natural $\delta$-function for one-sided-in-time transforms is complex. It follows from Eqs. (D7) and (D13) that
$${\delta _\infty}(\omega) \def\LDeqtab{}= \delta (\omega)/2 + i/2\pi \omega ,$$
$${\delta _T}(\omega) \def\LDeqtab{}= \sin (\omega T)/2\pi \omega + i[1 - \cos (\omega T)]/2\pi \omega .$$
By substituting formula (D22) in Eq. (D20) and using Eq. (D19), one finds that
$$\begin{split}F(\omega)& = \int_{- \infty}^\infty [\delta (\omega - \omega ^\prime) + i/\pi (\omega - \omega ^\prime)]F(\omega ^\prime){\rm d}\omega ^\prime /2\\ &= [F(\omega) + F(\omega)]/2.\\[-1.3pc]\end{split}$$
Thus, the $\delta$-function defined by Eq. (D22) sifts any function of frequency that is the infinite-time transform of a causal function. A similar remark applies to the $\delta$-function defined by Eq. (D23).

## APPENDIX E: DERIVATIVE RULE

One can reconcile the frequency- and time-domain analyses of Sections 35 if one applies the derivative rule carefully. We will prove this statement for an undamped oscillator, because the frequency-domain analysis of this system was the most problematic.

First, recall that the time-domain solution of Eq. (73) is

$$A(t) = \int_0^t R(s){\rm d}s,$$
from which it follows that
$$\langle |A(t{)|^2}\rangle = \int_0^t \int_0^t \langle {R^*}(s)R(s^\prime)\rangle {\rm d}s{\rm d}s^\prime = \sigma t.$$
Second, observe that the Fourier transform of Eq. (73) is
$${-}i\omega A(\omega) = R(\omega) - A(T){e^{i\omega T}},$$
from which it follows that
$$\begin{split}\!\!{\omega ^2}\langle |A(\omega {)|^2}\rangle &= \langle |R(\omega {)|^2}\rangle - 2{\rm{Re}}[\langle {R^*}(\omega)A(T){e^{i\omega T}}\rangle] \!\!\\&\quad +\langle |A(T{)|^2}\rangle .\end{split}$$
According to Eqs. (10) and (E2), the first and third terms on the right side of Eq. (E4), $\langle |R(\omega {)|^2}\rangle$ and $\langle |A(T{)|^2}\rangle$, both equal $\sigma T$. The second term
$$\begin{split}\langle {R^*}(\omega)A(T){e^{i\omega T}}\rangle &= \int_0^T \int_0^T \langle {R^*}(s)R(s^\prime)\rangle {e^{i\omega (T - s)}}{\rm d}s{\rm d}s^\prime \\ &= \sigma ({e^{i\omega T}} - 1)/i\omega .\\[-1.3pc]\end{split}$$
By assembling these results, one finds that
$$\begin{split}{\omega ^2}\langle |A(\omega {)|^2}\rangle& = \sigma T - 2\sigma \sin (\omega T)/\omega + \sigma T\\ &= 2\sigma T[1 - {\rm{sinc}}(\omega T)].\end{split}$$
Equation (E6) is equivalent to Eq. (80). It shows that the terms required to correct Eq. (86) come from the final value $A(T)$, which is often omitted.

This proof of consistency depends on prior knowledge of the time-domain solution. It is also possible to work entirely in the frequency domain. By setting $\omega = 0$ in Eq. (E3), one finds that $A(T) = R(0) = \sigma T$. Hence, one can rewrite Eq. (E3) in the alternative form

$${-}i\omega A(\omega) = R(\omega) - R(0){e^{i\omega T}},$$
in which the terms on the right side are known. Equation (E6) follows directly from Eq. (E7).

## APPENDIX F: SPECIFIC FREQUENCY CORRELATIONS

Let $R(t)$ be the source function for noise, which has properties (7) and (8). For white noise, the two-time correlation

$$\langle {R^*}(t)R(t^\prime)\rangle = \sigma \delta (t - t^\prime).$$
White noise is $\delta$-correlated in time. The associated two-frequency correlation
$$\begin{split}\langle {R^*}(\omega)R(\omega ^\prime)\rangle &= \int_0^T \int_0^T \langle {R^*}(t)R(t^\prime)\rangle {e^{i(\omega ^\prime t^\prime - \omega t)}}{\rm d}t{\rm d}t^\prime \\[-2pt]& = \sigma \int_0^T {e^{i(\omega ^\prime - \omega)t}}{\rm d}t\\[-2pt] &= \sigma [{e^{i(\omega ^\prime - \omega)T}} - 1]/i(\omega ^\prime - \omega)\\[-2pt] &= 2\pi \sigma {\delta _T}(\omega ^\prime - \omega),\end{split}$$
where ${\delta _T}$ is the effective $\delta$-function that was introduced in Appendix D. (See Fig. 11.) Thus, white noise is also $\delta$-correlated in frequency, as stated after Eq. (B11). Notice that interchanging $\omega$ and $\omega ^\prime$ conjugates the correlation. Notice also that $\langle {R^*}(\omega)R(\omega ^\prime)\rangle \to \sigma T$ as $\omega ^\prime \to \omega$. This result is consistent with Eq. (10).

Now let $A(t)$ be the displacement of a damped oscillator. In Sections 3 and 4, and Appendix B, it was shown that $A(\omega) \sim G(\omega)R(\omega)$, where $A(\omega)$ and $G(\omega)$ are the transformed displacement and Green function, respectively. For such a system, the frequency correlation

$$C(\omega ,\omega ^\prime) \sim 2\pi \sigma {G^*}(\omega)G(\omega ^\prime){\delta _T}(\omega ^\prime - \omega).$$
Thus, damped oscillators are also (asymptotically) $\delta$-correlated in frequency. Notice that $C(\omega ^\prime ,\omega) = {C^*}(\omega ,\omega ^\prime)$.

For colored noise, the temporal correlation

$$\langle {R^*}(t)R(t^\prime)\rangle = (\sigma \nu /2){e^{- \nu |t - t^\prime |}}.$$
The associated frequency correlation
$$\langle {R^*}(\omega)R(\omega ^\prime)\rangle = (\sigma \nu /2)\int_0^T \int_0^T {e^{- \nu |t - t^\prime |}}{e^{i(\omega ^\prime t^\prime - \omega t)}}{\rm d}t{\rm d}t^\prime .$$
On the right side of Eq. (F5), the $t^\prime$-integral is
$$\begin{split}&{e^{- (\nu + i\omega)t}}\int_0^t {e^{(\nu + i\omega ^\prime)t^\prime}}{\rm d}t^\prime + {e^{(\nu - i\omega)t}}\int_t^T {e^{- (\nu - i\omega ^\prime)t^\prime}}{\rm d}t^\prime \\[-2pt]&\quad= \frac{{{e^{i(\omega ^\prime - \omega)t}} - {e^{- (\nu + i\omega)t}}}}{{(\nu + i\omega ^\prime)}} + \frac{{{e^{i(\omega ^\prime - \omega)t}} - {e^{- (\nu - i\omega ^\prime)T}}{e^{(\nu - i\omega)t}}}}{{(\nu - i\omega ^\prime)}}.\end{split}$$
The $t$-integral of the first term on the right side of Eq. (F6) is
$$\frac{{{e^{i(\omega ^\prime - \omega)T}} - 1}}{{i(\omega ^\prime - \omega)(\nu + i\omega ^\prime)}} - \frac{{1 - {e^{- (\nu + i\omega)T}}}}{{(\nu + i\omega)(\nu + i\omega ^\prime)}}$$
and the $t$-integral of the second term is
$$\frac{{{e^{i(\omega ^\prime - \omega)T}} - 1}}{{i(\omega ^\prime - \omega)(\nu - i\omega ^\prime)}} - \frac{{{e^{i(\omega ^\prime - \omega)T}} - {e^{- (\nu - i\omega ^\prime)T}}}}{{(\nu - i\omega)(\nu - i\omega ^\prime)}}.$$
By adding 1 to, and subtracting 1 from, the numerator of the second term in integral (F8), and rewriting $i(\omega ^\prime - \omega)$ as $(\nu - i\omega) - (\nu - i\omega ^\prime)$, one obtains the contribution
$$\begin{split}&\frac{{{e^{i(\omega ^\prime - \omega)T}} - 1}}{{i(\omega ^\prime - \omega)}}\left[{\frac{1}{{\nu + i\omega ^\prime}} + \frac{1}{{\nu - i\omega ^\prime}} - \frac{{i(\omega ^\prime - \omega)}}{{(\nu - i\omega)(\nu - i\omega ^\prime)}}} \right]\\&\quad = \frac{{{e^{i(\omega ^\prime - \omega)T}} - 1}}{{i(\omega ^\prime - \omega)}}\frac{{2\nu + i(\omega ^\prime - \omega)}}{{(\nu + i\omega ^\prime)(\nu - i\omega)}}.\\[-1.3pc]\end{split}$$
The remaining contribution is
$$\frac{{{e^{- (\nu + i\omega)T}} - 1}}{{(\nu + i\omega)(\nu + i\omega ^\prime)}} + \frac{{{e^{- (\nu - i\omega ^\prime)T}} - 1}}{{(\nu - i\omega)(\nu - i\omega ^\prime)}}.$$
Notice that interchanging $\omega$ and $\omega ^\prime$ conjugates contributions (F9) and (F10). If one neglects the exponentially small terms in the latter contribution, it reduces to
$${-}\frac{{2({\nu ^2} - \omega \omega ^\prime)}}{{({\nu ^2} + {\omega ^2})[{\nu ^2} + {{(\omega ^\prime)}^2}]}}.$$
By combining these contributions and multiplying the result by $\sigma \nu /2$, one obtains the frequency correlation
$$\begin{split}C(\omega ,\omega ^\prime) &= \frac{{\sigma \nu}}{2}\left\{\frac{{{e^{i(\omega ^\prime - \omega)T}} - 1}}{{i(\omega ^\prime - \omega)}}\frac{{2\nu + i(\omega ^\prime - \omega)}}{{(\nu + i\omega ^\prime)(\nu - i\omega)}} \right.\\&\quad\left.- \frac{{2({\nu ^2} - \omega \omega ^\prime)}}{{({\nu ^2} + {\omega ^2})[{\nu ^2} + {{(\omega ^\prime)}^2}]}} \right\}.\end{split}$$
Notice that Eq. (F12) is consistent with Eq. (16).

The normalized correlation, which is correlation (F12) divided by $\sigma T$, is illustrated in Fig. 12. Its real part is a symmetric function of $\omega$ and $\omega ^\prime$, whereas its imaginary part is an asymmetric function. The first term on the right side of Eq. (F12), which is proportional to the effective $\delta$-function, is only significant for a narrow range of frequencies ($|\omega ^\prime - \omega | \sim 1/T$). Although the second term is significant for a wide range of frequencies ($\omega$, $\omega ^\prime \sim \nu$), it is smaller than the first term by a factor of $\nu T$. In this sense, colored noise is $\delta$-correlated in frequency.

Fig. 12. Real part (top) and imaginary part (bottom) of the frequency correlation of colored noise for $\nu T = 10$. For $|\omega ^\prime - \omega | \gg 1/T$, the correlation is weak. For $\omega ^\prime = \omega$, the correlation is real and equals the spectrum shown in Fig. 1(b).

Fig. 13. Real part (top) and imaginary part (bottom) of the frequency correlation of an undamped oscillator. For $|\omega ^\prime - \omega | \gg 1/T$, the correlation is weak. For $\omega ^\prime = \omega$, the correlation is real and equals the spectrum shown in Fig. 8.

In Section 5, it was shown that the temporal correlation of an undamped oscillator is $\sigma \min (t,t^\prime)$. One can calculate the associated frequency correlation by doing the $t$- and $t^\prime$-integrals in Eq. (B1) directly. Alternatively, by combining Eqs. (81) and (F1), one finds that

$$\begin{split}C(\omega ,\omega ^\prime) &= \sigma \int_0^T \frac{{{e^{- i\omega T}} - {e^{- i\omega s}}}}{{- i\omega}}\frac{{{e^{i\omega ^\prime T}} - {e^{i\omega ^\prime s}}}}{{i\omega ^\prime}}{\rm d}s\\ &= (\sigma /\omega \omega ^\prime)\int_0^T \{{e^{i(\omega ^\prime - \omega)T}} - {e^{i\omega ^\prime s - i\omega T}}\\&\quad - {e^{i\omega ^\prime T - i\omega s}} + {e^{i(\omega ^\prime - \omega)s}}\} {\rm d}s.\end{split}$$
The combined integral of the first and fourth terms on the right side of Eq. (F13) is
$$T{e^{i(\omega ^\prime - \omega)T}} + \frac{{{e^{i(\omega ^\prime - \omega)T}} - 1}}{{i(\omega ^\prime - \omega)}},$$
whereas the combined integral of the second and third terms is
$$\frac{{{e^{- i\omega T}}(1 - {e^{i\omega ^\prime T}})}}{{i\omega ^\prime}} + \frac{{{e^{i\omega ^\prime T}}({e^{- i\omega T}} - 1)}}{{i\omega}}.$$
Notice that interchanging $\omega$ and $\omega ^\prime$ conjugates contributions (F14) and (F15). By combining them and multiplying the result by $\sigma /\omega \omega ^\prime$, one obtains the frequency correlation
$$\begin{split}C(\omega ,\omega ^\prime)& = \frac{\sigma}{{\omega \omega ^\prime}}\left[{T{e^{i(\omega ^\prime - \omega)T}} + \frac{{{e^{i(\omega ^\prime - \omega)T}} - 1}}{{i(\omega ^\prime - \omega)}}} \right.\\ &\quad\left. + \frac{{{e^{- i\omega T}}(1 - {e^{i\omega ^\prime T}})}}{{i\omega ^\prime}} + \frac{{{e^{i\omega ^\prime T}}({e^{- i\omega T}} - 1)}}{{i\omega}} \right].\end{split}$$
For very low frequencies ($\omega$, $\omega ^\prime \ll 1/T$), the zeroth- and first-order terms in the square brackets cancel, so the correlation is not singular, and the second-order term is ${T^3}\omega \omega ^\prime /3$. Hence, $C(\omega ,\omega ^\prime) \to \sigma {T^3}/3$ as $\omega$, $\omega ^\prime \to 0$. Notice that Eq. (F16) is consistent with Eq. (80).

The normalized correlation, which is correlation (F16) divided by $\sigma {T^3}/3$, is illustrated in Fig. 13. Its real part is a symmetric function of $\omega$ and $\omega ^\prime$, whereas its imaginary part is an asymmetric function. If one sections the real and imaginary correlation surfaces along the lines $\omega ^\prime = - \omega$, the curves that result are similar to those in Fig. 11. Although the heights and depths of the peaks and troughs are independent of $T$, the range of frequencies for which the correlation is significant decreases as $T$ increases. In this sense, the fluctuations of an undamped oscillator are $\delta$-correlated in frequency.

## Funding

Natural Sciences and Engineering Research Council of Canada.

## Acknowledgment

We thank the reviewers for bringing to our attention [17,18,20,22]. TJS and ASH were supported by NSERC.

## Disclosures

The authors declare no conflicts of interest.

## Data Availability

No data were generated or analyzed in this research.

## REFERENCES

1. A. E. Siegman, Lasers (University Science Books, 1986).

2. P. W. Milonni and J. H. Eberly, Lasers (Wiley, 1988).

3. G. P. Agrawal and N. K. Dutta, Semiconductor Lasers, 2nd ed. (Van Nostrand Reinhold, 1993).

4. L. A. Coldren, S. W. Corzine, and M. L. Masanovic, Diode Lasers and Photonic Integrated Circuits, 2nd ed. (Wiley, 2012).

5. C. W. Gardiner, Handbook of Stochastic Methods, 2nd ed. (Springer, 1985).

6. J. W. Goodman, Statistical Optics (Wiley, 1985).

7. A. Papoulis, Probability, Random Variables and Stochastic Processes, 3rd ed. (McGraw-Hill, 1991).

8. D. S. Lemons, An Introduction to Stochastic Processes in Physics (Johns Hopkins University, 2002).

9. N. Wiener, “Generalized harmonic analysis,” Acta Math. 55, 117–258 (1930). [CrossRef]

10. A. Khintchine, “Korrelationstheorie der stationären stochastischen Prozesse,” Math. Ann. 109, 604–615 (1934). [CrossRef]

11. D. S. Lemons and A. Gythiel, “Paul Langevin’s 1908 paper ‘On the Theory of Brownian Motion’ [‘Sur la théorie du mouvement brownien,’ C. R. Acad. Sci. (Paris) 146, 530–533 (1908)],” Am. J. Phys. 65, 1079–1081 (1997). [CrossRef]

12. L. S. Ornstein, “On the Brownian motion,” Proc. Acad. Amst. 21,96–108 (1919).

13. G. E. Uhlenbeck and L. S. Ornstein, “On the theory of Brownian motion,” Phys. Rev. 36, 823–841 (1930). [CrossRef]

14. C. J. McKinstrie, “Stochastic and probabilistic equations for three- and four-level lasers,” J. Opt. Soc. Am. B 37, 1333–1358 (2020). [CrossRef]  The referenced integrals are in Appendix E.

15. C. J. McKinstrie and J. P. Gordon, “Field fluctuations produced by parametric processes in fibers,” IEEE J. Sel. Top. Quantum Electron. 18, 958–969 (2012). [CrossRef]

16. W. H. Louisell, Quantum Statistical Properties of Radiation (Wiley, 1973), Sec. 2.5.

17. H. C. Berg, Random Walks in Biology (Princeton University, 1983).

18. N. G. van Kampen, Stochastic Processes in Physics and Chemistry, 3rd ed. (Elsevier, 2007).

19. M. Baxter and A. Rennie, Financial Calculus (Cambridge, 1996).

20. K. Marti, Stochastic Optimization Methods, 3rd ed. (Springer, 2015).

21. C. J. McKinstrie, T. J. Stirling, and A. S. Helmy, “Laser linewidths: tutorial,” J. Opt. Soc. Am. B 38, 3837–3848 (2021). [CrossRef]

22. M. Gitterman, The Noisy Oscillator (World Scientific, 2005).

23. R. N. Bracewell, The Fourier Transform and its Applications, 3rd ed. (McGraw-Hill, 2000).

24. B. Y. K. Hu, “Kramers–Kronig in two lines,” Am. J. Phys. 57, 821 (1989). [CrossRef]

25. J. D. Jackson, Classical Electrodynamics, 3rd ed. (Wiley, 1999).

### References

• View by:

1. A. E. Siegman, Lasers (University Science Books, 1986).
2. P. W. Milonni and J. H. Eberly, Lasers (Wiley, 1988).
3. G. P. Agrawal and N. K. Dutta, Semiconductor Lasers, 2nd ed. (Van Nostrand Reinhold, 1993).
4. L. A. Coldren, S. W. Corzine, and M. L. Masanovic, Diode Lasers and Photonic Integrated Circuits, 2nd ed. (Wiley, 2012).
5. C. W. Gardiner, Handbook of Stochastic Methods, 2nd ed. (Springer, 1985).
6. J. W. Goodman, Statistical Optics (Wiley, 1985).
7. A. Papoulis, Probability, Random Variables and Stochastic Processes, 3rd ed. (McGraw-Hill, 1991).
8. D. S. Lemons, An Introduction to Stochastic Processes in Physics (Johns Hopkins University, 2002).
9. N. Wiener, “Generalized harmonic analysis,” Acta Math. 55, 117–258 (1930).
[Crossref]
10. A. Khintchine, “Korrelationstheorie der stationären stochastischen Prozesse,” Math. Ann. 109, 604–615 (1934).
[Crossref]
11. D. S. Lemons and A. Gythiel, “Paul Langevin’s 1908 paper ‘On the Theory of Brownian Motion’ [‘Sur la théorie du mouvement brownien,’ C. R. Acad. Sci. (Paris) 146, 530–533 (1908)],” Am. J. Phys. 65, 1079–1081 (1997).
[Crossref]
12. L. S. Ornstein, “On the Brownian motion,” Proc. Acad. Amst. 21,96–108 (1919).
13. G. E. Uhlenbeck and L. S. Ornstein, “On the theory of Brownian motion,” Phys. Rev. 36, 823–841 (1930).
[Crossref]
14. C. J. McKinstrie, “Stochastic and probabilistic equations for three- and four-level lasers,” J. Opt. Soc. Am. B 37, 1333–1358 (2020).The referenced integrals are in Appendix E.
[Crossref]
15. C. J. McKinstrie and J. P. Gordon, “Field fluctuations produced by parametric processes in fibers,” IEEE J. Sel. Top. Quantum Electron. 18, 958–969 (2012).
[Crossref]
16. W. H. Louisell, Quantum Statistical Properties of Radiation (Wiley, 1973), Sec. 2.5.
17. H. C. Berg, Random Walks in Biology (Princeton University, 1983).
18. N. G. van Kampen, Stochastic Processes in Physics and Chemistry, 3rd ed. (Elsevier, 2007).
19. M. Baxter and A. Rennie, Financial Calculus (Cambridge, 1996).
20. K. Marti, Stochastic Optimization Methods, 3rd ed. (Springer, 2015).
21. C. J. McKinstrie, T. J. Stirling, and A. S. Helmy, “Laser linewidths: tutorial,” J. Opt. Soc. Am. B 38, 3837–3848 (2021).
[Crossref]
22. M. Gitterman, The Noisy Oscillator (World Scientific, 2005).
23. R. N. Bracewell, The Fourier Transform and its Applications, 3rd ed. (McGraw-Hill, 2000).
24. B. Y. K. Hu, “Kramers–Kronig in two lines,” Am. J. Phys. 57, 821 (1989).
[Crossref]
25. J. D. Jackson, Classical Electrodynamics, 3rd ed. (Wiley, 1999).

#### 2012 (1)

C. J. McKinstrie and J. P. Gordon, “Field fluctuations produced by parametric processes in fibers,” IEEE J. Sel. Top. Quantum Electron. 18, 958–969 (2012).
[Crossref]

#### 1997 (1)

D. S. Lemons and A. Gythiel, “Paul Langevin’s 1908 paper ‘On the Theory of Brownian Motion’ [‘Sur la théorie du mouvement brownien,’ C. R. Acad. Sci. (Paris) 146, 530–533 (1908)],” Am. J. Phys. 65, 1079–1081 (1997).
[Crossref]

#### 1989 (1)

B. Y. K. Hu, “Kramers–Kronig in two lines,” Am. J. Phys. 57, 821 (1989).
[Crossref]

#### 1934 (1)

A. Khintchine, “Korrelationstheorie der stationären stochastischen Prozesse,” Math. Ann. 109, 604–615 (1934).
[Crossref]

#### 1930 (2)

N. Wiener, “Generalized harmonic analysis,” Acta Math. 55, 117–258 (1930).
[Crossref]

G. E. Uhlenbeck and L. S. Ornstein, “On the theory of Brownian motion,” Phys. Rev. 36, 823–841 (1930).
[Crossref]

#### 1919 (1)

L. S. Ornstein, “On the Brownian motion,” Proc. Acad. Amst. 21,96–108 (1919).

#### Agrawal, G. P.

G. P. Agrawal and N. K. Dutta, Semiconductor Lasers, 2nd ed. (Van Nostrand Reinhold, 1993).

#### Baxter, M.

M. Baxter and A. Rennie, Financial Calculus (Cambridge, 1996).

#### Berg, H. C.

H. C. Berg, Random Walks in Biology (Princeton University, 1983).

#### Bracewell, R. N.

R. N. Bracewell, The Fourier Transform and its Applications, 3rd ed. (McGraw-Hill, 2000).

#### Coldren, L. A.

L. A. Coldren, S. W. Corzine, and M. L. Masanovic, Diode Lasers and Photonic Integrated Circuits, 2nd ed. (Wiley, 2012).

#### Corzine, S. W.

L. A. Coldren, S. W. Corzine, and M. L. Masanovic, Diode Lasers and Photonic Integrated Circuits, 2nd ed. (Wiley, 2012).

#### Dutta, N. K.

G. P. Agrawal and N. K. Dutta, Semiconductor Lasers, 2nd ed. (Van Nostrand Reinhold, 1993).

#### Eberly, J. H.

P. W. Milonni and J. H. Eberly, Lasers (Wiley, 1988).

#### Gardiner, C. W.

C. W. Gardiner, Handbook of Stochastic Methods, 2nd ed. (Springer, 1985).

#### Gitterman, M.

M. Gitterman, The Noisy Oscillator (World Scientific, 2005).

#### Goodman, J. W.

J. W. Goodman, Statistical Optics (Wiley, 1985).

#### Gordon, J. P.

C. J. McKinstrie and J. P. Gordon, “Field fluctuations produced by parametric processes in fibers,” IEEE J. Sel. Top. Quantum Electron. 18, 958–969 (2012).
[Crossref]

#### Gythiel, A.

D. S. Lemons and A. Gythiel, “Paul Langevin’s 1908 paper ‘On the Theory of Brownian Motion’ [‘Sur la théorie du mouvement brownien,’ C. R. Acad. Sci. (Paris) 146, 530–533 (1908)],” Am. J. Phys. 65, 1079–1081 (1997).
[Crossref]

#### Hu, B. Y. K.

B. Y. K. Hu, “Kramers–Kronig in two lines,” Am. J. Phys. 57, 821 (1989).
[Crossref]

#### Jackson, J. D.

J. D. Jackson, Classical Electrodynamics, 3rd ed. (Wiley, 1999).

#### Khintchine, A.

A. Khintchine, “Korrelationstheorie der stationären stochastischen Prozesse,” Math. Ann. 109, 604–615 (1934).
[Crossref]

#### Lemons, D. S.

D. S. Lemons and A. Gythiel, “Paul Langevin’s 1908 paper ‘On the Theory of Brownian Motion’ [‘Sur la théorie du mouvement brownien,’ C. R. Acad. Sci. (Paris) 146, 530–533 (1908)],” Am. J. Phys. 65, 1079–1081 (1997).
[Crossref]

D. S. Lemons, An Introduction to Stochastic Processes in Physics (Johns Hopkins University, 2002).

#### Louisell, W. H.

W. H. Louisell, Quantum Statistical Properties of Radiation (Wiley, 1973), Sec. 2.5.

#### Marti, K.

K. Marti, Stochastic Optimization Methods, 3rd ed. (Springer, 2015).

#### Masanovic, M. L.

L. A. Coldren, S. W. Corzine, and M. L. Masanovic, Diode Lasers and Photonic Integrated Circuits, 2nd ed. (Wiley, 2012).

#### McKinstrie, C. J.

C. J. McKinstrie and J. P. Gordon, “Field fluctuations produced by parametric processes in fibers,” IEEE J. Sel. Top. Quantum Electron. 18, 958–969 (2012).
[Crossref]

#### Milonni, P. W.

P. W. Milonni and J. H. Eberly, Lasers (Wiley, 1988).

#### Ornstein, L. S.

G. E. Uhlenbeck and L. S. Ornstein, “On the theory of Brownian motion,” Phys. Rev. 36, 823–841 (1930).
[Crossref]

L. S. Ornstein, “On the Brownian motion,” Proc. Acad. Amst. 21,96–108 (1919).

#### Papoulis, A.

A. Papoulis, Probability, Random Variables and Stochastic Processes, 3rd ed. (McGraw-Hill, 1991).

#### Rennie, A.

M. Baxter and A. Rennie, Financial Calculus (Cambridge, 1996).

#### Siegman, A. E.

A. E. Siegman, Lasers (University Science Books, 1986).

#### Uhlenbeck, G. E.

G. E. Uhlenbeck and L. S. Ornstein, “On the theory of Brownian motion,” Phys. Rev. 36, 823–841 (1930).
[Crossref]

#### van Kampen, N. G.

N. G. van Kampen, Stochastic Processes in Physics and Chemistry, 3rd ed. (Elsevier, 2007).

#### Wiener, N.

N. Wiener, “Generalized harmonic analysis,” Acta Math. 55, 117–258 (1930).
[Crossref]

#### Acta Math. (1)

N. Wiener, “Generalized harmonic analysis,” Acta Math. 55, 117–258 (1930).
[Crossref]

#### Am. J. Phys. (2)

D. S. Lemons and A. Gythiel, “Paul Langevin’s 1908 paper ‘On the Theory of Brownian Motion’ [‘Sur la théorie du mouvement brownien,’ C. R. Acad. Sci. (Paris) 146, 530–533 (1908)],” Am. J. Phys. 65, 1079–1081 (1997).
[Crossref]

B. Y. K. Hu, “Kramers–Kronig in two lines,” Am. J. Phys. 57, 821 (1989).
[Crossref]

#### IEEE J. Sel. Top. Quantum Electron. (1)

C. J. McKinstrie and J. P. Gordon, “Field fluctuations produced by parametric processes in fibers,” IEEE J. Sel. Top. Quantum Electron. 18, 958–969 (2012).
[Crossref]

#### Math. Ann. (1)

A. Khintchine, “Korrelationstheorie der stationären stochastischen Prozesse,” Math. Ann. 109, 604–615 (1934).
[Crossref]

#### Phys. Rev. (1)

G. E. Uhlenbeck and L. S. Ornstein, “On the theory of Brownian motion,” Phys. Rev. 36, 823–841 (1930).
[Crossref]

L. S. Ornstein, “On the Brownian motion,” Proc. Acad. Amst. 21,96–108 (1919).

#### Other (16)

W. H. Louisell, Quantum Statistical Properties of Radiation (Wiley, 1973), Sec. 2.5.

H. C. Berg, Random Walks in Biology (Princeton University, 1983).

N. G. van Kampen, Stochastic Processes in Physics and Chemistry, 3rd ed. (Elsevier, 2007).

M. Baxter and A. Rennie, Financial Calculus (Cambridge, 1996).

K. Marti, Stochastic Optimization Methods, 3rd ed. (Springer, 2015).

A. E. Siegman, Lasers (University Science Books, 1986).

P. W. Milonni and J. H. Eberly, Lasers (Wiley, 1988).

G. P. Agrawal and N. K. Dutta, Semiconductor Lasers, 2nd ed. (Van Nostrand Reinhold, 1993).

L. A. Coldren, S. W. Corzine, and M. L. Masanovic, Diode Lasers and Photonic Integrated Circuits, 2nd ed. (Wiley, 2012).

C. W. Gardiner, Handbook of Stochastic Methods, 2nd ed. (Springer, 1985).

J. W. Goodman, Statistical Optics (Wiley, 1985).

A. Papoulis, Probability, Random Variables and Stochastic Processes, 3rd ed. (McGraw-Hill, 1991).

D. S. Lemons, An Introduction to Stochastic Processes in Physics (Johns Hopkins University, 2002).

M. Gitterman, The Noisy Oscillator (World Scientific, 2005).

R. N. Bracewell, The Fourier Transform and its Applications, 3rd ed. (McGraw-Hill, 2000).

J. D. Jackson, Classical Electrodynamics, 3rd ed. (Wiley, 1999).

### Data Availability

No data were generated or analyzed in this research.

### Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

### Figures (13)

Fig. 1. Temporal correlation function (top) and spectral energy density (bottom) of colored noise. The blue density curve represents the first term in Eq. (16), whereas the red curve denotes both terms. For $\nu T = 10$, the finite-time correction is noticeable.
Fig. 2. Two-time correlation function of a strongly damped oscillator. For times that are longer than the damping time, the correlation depends on only the time difference $t^\prime - t$.
Fig. 3. Temporal correlation function (top) and spectral energy density (bottom) of a strongly damped oscillator. The blue density curve represents the first term in Eq. (28), whereas the red curve denotes all three terms. For $\nu T = 10$, the finite-time and transient corrections are noticeable.
Fig. 4. Two-time correlation function of a weakly damped oscillator for ${\nu _0}/{\omega _0} = 0.2$. For times that are longer than the damping time, the correlation depends on only the time difference $t^\prime - t$.
Fig. 5. Temporal correlation function (top) and spectral energy density (bottom) of a weakly damped oscillator for ${\nu _0}/{\omega _0} = 0.2$.
Fig. 6. Spectral energy density of a weakly damped oscillator for ${\nu _0}/{\omega _0} = 0.2$. The blue curve represents formula (65), whereas the red curve represents the Lorentzian formula (69).
Fig. 7. Two-time correlation function of an undamped oscillator. The correlation never becomes a function of the time difference $t^\prime - t$.
Fig. 8. Spectral energy density of an undamped oscillator. The blue curve denotes the exact spectrum (80), whereas the green and red curves denote Lorentzian approximations with $a{= 6^{1/2}}$ and 2, respectively.
Fig. 9. Integration regions in the $t - \tau$ plane for the correlation theorem (left) and the convolution theorem (right). The horizontal axis represents $t$ and the vertical axis represents $\tau$.
Fig. 10. Integration region in the $t - \tau$ plane for the Wiener–Khinchin theorem. The horizontal axis represents $t$ and the vertical axis represents $\tau$. The shaded region contributes to the integral.
Fig. 11. Real (top) and imaginary (bottom) parts of the Fourier transform of the step function. The dashed curves represent formula (D7), whereas the solid curves represent formula (D13).
Fig. 12. Real part (top) and imaginary part (bottom) of the frequency correlation of colored noise for $\nu T = 10$. For $|\omega ^\prime - \omega | \gg 1/T$, the correlation is weak. For $\omega ^\prime = \omega$, the correlation is real and equals the spectrum shown in Fig. 1(b).
Fig. 13. Real part (top) and imaginary part (bottom) of the frequency correlation of an undamped oscillator. For $|\omega ^\prime - \omega | \gg 1/T$, the correlation is weak. For $\omega ^\prime = \omega$, the correlation is real and equals the spectrum shown in Fig. 8.

### Equations (178)

$d t A = L ( A ) + R ( t ) ,$
$A ( ω ) = ∫ 0 T A ( t ) e i ω t d t ,$
$A ( t ) = ∫ − ∞ ∞ A ( ω ) e − i ω t d ω / 2 π .$
$∫ − ∞ ∞ e ± i ω t d ω / 2 π = δ ( t )$
$∫ 0 T | A ( t ) | 2 d t = ∫ − ∞ ∞ | A ( ω ) | 2 d ω / 2 π ,$
$S ( ω ) = ∫ 0 T ∫ 0 T C ( t , t ′ ) e i ω ( t ′ − t ) d t d t ′ .$
$⟨ R ( t ) ⟩ = 0 ,$
$⟨ R ( t ) R ( t ′ ) ⟩ = 0 ,$
$⟨ R ∗ ( t ) R ( t ′ ) ⟩ = σ δ ( t − t ′ ) ,$
$⟨ | R ( ω ) | 2 ⟩ = ∫ 0 T ∫ 0 T ⟨ R ∗ ( t ) R ( t ′ ) ⟩ e i ω ( t ′ − t ) d t d t ′ = ∫ 0 T σ d t = σ T .$
$⟨ R ∗ ( t ) R ( t ′ ) ⟩ = ( σ ν / 2 ) e − ν | t − t ′ |$
$∫ 0 T | R ( t ) | 2 d t = ∫ 0 T ( σ ν / 2 ) d t = σ ν T / 2$
$⟨ | R ( ω ) | 2 ⟩ = ∫ 0 T ∫ 0 T ( σ ν / 2 ) e − ν | t − t ′ | e i ω ( t ′ − t ) d t d t ′ .$
$∫ 0 t e − ( ν + i ω ) ( t − t ′ ) d t ′ + ∫ t T e − ( ν − i ω ) ( t ′ − t ) d t ′ = 1 − e − ( ν + i ω ) t ν + i ω + 1 − e − ( ν − i ω ) ( T − t ) ν − i ω .$
$− 1 − e − ( ν + i ω ) T ( ν + i ω ) 2 − 1 − e − ( ν − i ω ) T ( ν − i ω ) 2 .$
$⟨ | R ( ω ) | 2 ⟩ = σ ν [ ν T ν 2 + ω 2 − ν 2 − ω 2 ( ν 2 + ω 2 ) 2 ] ∼ σ ν 2 T ν 2 + ω 2 .$
$∫ − ∞ ∞ ⟨ | R ( ω ) | 2 ⟩ d ω / 2 π = σ ν T / 2.$
$d t A = − ( ν r + i ν i ) A + R a ( t ) ,$
$d t B = − ν r B + R b ( t ) ,$
$B ( t ) = ∫ 0 t G ( t − s ) R ( s ) d s ,$
$d t G = − ν r G + δ ( t ) ,$
$B ( t ) = ∫ 0 t e − ν ( t − s ) R ( s ) d s ,$
$⟨ | B ( t ) | 2 ⟩ = ∫ 0 t ∫ 0 t ⟨ R ∗ ( s ) R ( s ′ ) ⟩ e − ν ( 2 t − s − s ′ ) d s d s ′ = σ ∫ 0 t e − 2 ν ( t − s ) d s = σ ( 1 − e − 2 ν t ) / 2 ν .$
$∫ 0 T ⟨ | B ( t ) | 2 ⟩ d t = σ ( T / 2 ν − 1 / 4 ν 2 ) ,$
$⟨ B ∗ ( t ) B ( t ′ ) ⟩ = ∫ 0 t ∫ 0 t ′ ⟨ R ∗ ( s ) R ( s ′ ) ⟩ e − ν ( t + t ′ − s − s ′ ) d s d s ′ .$
$⟨ B ∗ ( t ) B ( t ′ ) ⟩ = σ ∫ 0 t e − ν ( t + t ′ − 2 s ) d s = σ [ e − ν ( t ′ − t ) − e − ν ( t ′ + t ) ] / 2 ν .$
$C ( τ ) ∼ σ e − ν | τ | / 2 ν .$
$⟨ | B ( ω ) | 2 ⟩ = σ T ν 2 + ω 2 − σ ( ν 2 − ω 2 ) ν ( ν 2 + ω 2 ) 2 − σ 2 ν ( ν 2 + ω 2 ) .$
$∫ − ∞ ∞ ⟨ | B ( ω ) | 2 ⟩ d ω / 2 π = σ ( T / 2 ν − 1 / 4 ν 2 ) .$
$⟨ R ( ω ) ⟩ = ∫ 0 T ⟨ R ( t ) ⟩ e i ω t d t = 0 ,$
$⟨ | R ( ω ) | 2 ⟩ = ∫ 0 T ∫ 0 T ⟨ R ∗ ( t ) R ( t ′ ) ⟩ e i ω ( t ′ − t ) d t d t ′ = σ T .$
$∫ 0 T [ d t B ( t ) ] e i ω t d t = B ( T ) e i ω T − B ( 0 ) − i ω ∫ 0 T B ( t ) e i ω t d t .$
$B ( ω ) = R ( ω ) / ( ν − i ω ) .$
$⟨ | B ( ω ) | 2 ⟩ = σ T / ( ν 2 + ω 2 ) .$
$⟨ A ∗ ( 0 ) A ( τ ) ⟩ ∼ ( σ / 2 ν r ) e − ν r | τ | − i ν i τ ,$
$⟨ | A ( ω ) | 2 ⟩ / T ∼ σ / [ ν r 2 + ( ω − ν i ) 2 ] ,$
$d t X = V ,$
$d t V = − 2 ν 0 V − ω 0 2 X + ϵ R ( t ) ,$
$X ( t ) = ϵ ∫ 0 t G x ( t − s ) R ( s ) d s ,$
$V ( t ) = ϵ ∫ 0 t G v ( t − s ) R ( s ) d s ,$
$G x ( t ) = H ( t ) sin ⁡ ( ω r t ) e − ν 0 t / ω r ,$
$G v ( t ) = H ( t ) [ cos ⁡ ( ω r t ) − ν 0 sin ⁡ ( ω r t ) / ω r ] e − ν 0 t ,$
$⟨ X 2 ( t ) ⟩ = ϵ 2 σ ∫ 0 t G x 2 ( t − s ) d s = ( ϵ 2 σ / 2 ω r 2 ) ∫ 0 t [ 1 − cos ⁡ ( 2 ω r s ) ] e − 2 ν 0 s d s .$
$∫ a b cos ⁡ ( α s ) e − β s d s = [ α sin ⁡ ( α s ) − β cos ⁡ ( α s ) ] e − β s | a b α 2 + β 2 ,$
$∫ a b sin ⁡ ( α s ) e − β s d s = − [ α cos ⁡ ( α s ) + β sin ⁡ ( α s ) ] e − β s | a b α 2 + β 2 .$
$⟨ X 2 ( t ) ⟩ = ϵ 2 σ 2 ω r 2 { 1 − e − 2 ν 0 t 2 ν 0 − ω r sin ⁡ ( 2 ω r t ) e − 2 ν 0 t 2 ω 0 2 − ν 0 [ 1 − cos ⁡ ( 2 ω r t ) e − 2 ν 0 t ] 2 ω 0 2 } .$
$⟨ X 2 ⟩ = ϵ 2 σ 2 ω r 2 ( 1 2 ν 0 − ν 0 2 ω 0 2 ) = ϵ 2 σ 4 ν 0 ω 0 2 .$
$( 1 / 2 ν 0 ) ∫ 0 T ( 1 − e − 2 ν 0 t ) d t = T / 2 ν 0 − 1 / 4 ν 0 2 ,$
$( ω r / 2 ω 0 2 ) ∫ 0 T sin ⁡ ( 2 ω r t ) e − 2 ν 0 t d t = ( ω r / 2 ω 0 2 ) 2 ,$
$( ν 0 / 2 ω 0 2 ) ∫ 0 T [ 1 − cos ⁡ ( 2 ω r t ) e − 2 ν 0 t ] d t = ( ν 0 / 2 ω 0 2 ) ( T − ν 0 / 2 ω 0 2 ) ,$
$∫ 0 T ⟨ X 2 ( t ) ⟩ d t ∼ ϵ 2 σ T 2 ω r 2 ( 1 2 ν 0 − ν 0 2 ω 0 2 ) = ϵ 2 σ T 4 ν 0 ω 0 2 ,$
$⟨ X ( t ) X ( t ′ ) ⟩ = ϵ 2 σ ∫ 0 t G x ( t − s ) G x ( t ′ − s ) d s = ( ϵ 2 σ / 2 ω r 2 ) ∫ 0 t { cos ⁡ [ ω r ( t ′ − t ) ] − cos ⁡ [ ω r ( t ′ + t − 2 s ) ] } e − ν 0 ( t ′ + t − 2 s ) d s .$
$I 1 ( t , t ′ ) = cos ⁡ [ ω r ( t ′ − t ) ] ∫ t ′ − t t ′ + t e − ν 0 u d u / 2 = cos ⁡ [ ω r ( t ′ − t ) ] [ e − ν 0 ( t ′ − t ) − e − ν 0 ( t ′ + t ) ] / 2 ν 0 .$
$I 2 ( t , t ′ ) = ∫ t ′ − t t + t ′ cos ⁡ ( ω r u ) e − ν 0 u d u / 2 = { ω r sin ⁡ [ ω r ( t ′ + t ) ] − ν 0 cos ⁡ [ ω r ( t ′ + t ) ] } × e − ν 0 ( t ′ + t ) / 2 ω 0 2 − { ω r sin ⁡ [ ω r ( t ′ − t ) ] − ν 0 cos ⁡ [ ω r ( t ′ − t ) ] } e − ν 0 ( t ′ − t ) / 2 ω 0 2 .$
$C ( t , t ′ ) = ( ϵ 2 σ / 2 ω r 2 ) [ I 1 ( t , t ′ ) − I 2 ( t , t ′ ) ] .$
$C ( τ ) ∼ ϵ 2 σ 4 ν 0 ω 0 2 [ cos ⁡ ( ω r τ ) + ν 0 ω r sin ⁡ ( ω r | τ | ) ] e − ν 0 | τ | ,$
$C ( τ ) ∝ ( ν 0 + i ω r ) e − ( ν 0 − i ω r ) | τ | − ( ν 0 − i ω r ) e − ( ν 0 + i ω r ) | τ | 2 i ω r .$
$1 2 i ω r [ ν ∗ 2 ν T ν 2 + ω 2 − ν 2 ν ∗ T ( ν ∗ ) 2 + ω 2 ] = 2 ω 0 2 T ω r I m [ 1 ν 2 + ω 2 ] .$
$S 1 ( ω ) = ϵ 2 σ T ( ω 0 2 − ω 2 ) 2 + 4 ν 0 2 ω 2 .$
$− ν 0 ( t ′ + t ) + i ( ω ± ω r ) ( t ′ − t ) = − [ ν 0 + i ( ω ± ω r ) ] t − [ ν 0 − i ( ω ± ω r ) ] t ′ .$
$− ( ν 0 ∓ i ω r ) ( t ′ + t ) + i ω ( t ′ − t ) = − ( ν 0 ∓ i ω r + i ω ) t − ( ν 0 ∓ i ω r − i ω ) t ′ .$
$[ − i ω − 1 ω 0 2 − i ω + 2 ν 0 ] [ X ( ω ) V ( ω ) ] = [ 0 ϵ R ( ω ) ] .$
$[ X ( ω ) V ( ω ) ] = 1 Δ ( ω ) [ ϵ R ( ω ) − i ω ϵ R ( ω ) ] ,$
$Δ ( ω ) = ω 0 2 − 2 i ν 0 ω − ω 2 .$
$⟨ | X ( ω ) | 2 ⟩ = ϵ 2 σ T ( ω 0 2 − ω 2 ) 2 + 4 ν 0 2 ω 2 .$
$∫ − ∞ ∞ ⟨ | X ( ω ) | 2 ⟩ d ω / 2 π = ϵ 2 σ T / 4 ν 0 ω 0 2 .$
$∫ − ∞ ∞ ⟨ | V ( ω ) | 2 ⟩ d ω / 2 π = ϵ 2 σ T / 4 ν 0 .$
$⟨ V 2 ⟩ = ϵ 2 σ / 4 ν 0 = ω 0 2 ⟨ X 2 ⟩ .$
$S ( δ ) ≈ ϵ 2 σ T / 4 ω 0 2 ( ν 0 2 + δ 2 ) .$
$( d tt 2 + 2 ν 0 d t + ω 0 2 ) X = ϵ R ( t ) .$
$X ( t ) = B ( t ) e − i ω 0 t + B ∗ ( t ) e i ω 0 t ,$
$( d t + ν 0 ) B ≈ R b ( t ) ,$
$d t A = R a ( t ) ,$
$A ( t ) = ∫ 0 t R ( s ) d s .$
$⟨ A ∗ ( t ) A ( t ′ ) ⟩ = ∫ 0 t ∫ 0 t ′ ⟨ R ∗ ( s ) R ( s ′ ) ⟩ d s ′ d s = σ min ( t , t ′ ) .$
$∫ 0 T ⟨ | A ( t ) | 2 ⟩ d t = σ T 2 / 2$
$⟨ | A ( ω ) | 2 ⟩ = ∫ 0 T ∫ 0 T σ min ( t , t ′ ) e i ω ( t ′ − t ) d t d t ′ .$
$∫ 0 t t ′ e i ω t ′ d t ′ + ∫ t T t e i ω t ′ d t ′ = t e i ω t i ω + e i ω t − 1 ω 2 + t ( e i ω T − e i ω t ) i ω = e i ω t − 1 ω 2 + t e i ω T i ω ,$
$∫ 0 T [ 1 − e − i ω t ω 2 + t e i ω ( T − t ) i ω ] d t = T ω 2 + e − i ω T − 1 i ω 3 + T ω 2 + 1 − e i ω T i ω 3 = 2 T [ 1 − s i n c ( ω T ) ] ω 2 ,$
$⟨ | A ( ω ) | 2 ⟩ = 2 σ T [ 1 − s i n c ( ω T ) ] ω 2 .$
$A ( ω ) = ∫ 0 T R ( s ) ∫ s T e i ω t d t d s = ∫ 0 T e i ω T − e i ω s i ω R ( s ) d s ,$
$⟨ | A ( ω ) | 2 ⟩ = ∫ 0 T ∫ 0 T e − i ω T − e − i ω s − i ω e i ω T − e i ω s ′ i ω × ⟨ R ∗ ( s ) R ( s ′ ) ⟩ d s d s ′ = σ ∫ 0 T e − i ω T − e − i ω s − i ω e i ω T − e i ω s i ω d s = 2 σ ω 2 ∫ 0 T { 1 − cos ⁡ [ ω ( s − T ) ] } d s = 2 σ T [ 1 − s i n c ( ω T ) ] ω 2 ,$
$∫ − ∞ ∞ | A ( ω ) | 2 d ω 2 π = σ T 2 π ∫ − ∞ ∞ 1 − s i n c ( x ) x 2 d x .$
$∫ − ∞ ∞ | A ( ω ) | 2 d ω / 2 π = σ T 2 / 2.$
$A ( ω ) = i R ( ω ) / ω .$
$⟨ | A ( ω ) | 2 ⟩ = σ T / ω 2 .$
$P ( x ) = exp ⁡ [ − ( x − μ ) 2 / 2 σ 2 ] / ( 2 π σ 2 ) 1 / 2 ,$
$⟨ x n ⟩ = ∫ − ∞ ∞ x n P ( x ) d x .$
$C ( k ) = ⟨ e ikx ⟩ = ∫ − ∞ ∞ P ( x ) e ikx d x$
$⟨ x n ⟩ = lim k → 0 ⁡ ( − i ∂ k ) n C ( k )$
$C ( k ) = exp ⁡ ( i k μ − σ 2 k 2 / 2 )$
$C ( K ) = exp ⁡ ( − K t M K / 2 ) ,$
$⟨ x 1 n 1 ⟩ = lim k 1 → 0 ⁡ ( − i ∂ k 1 ) n 1 C ( K ) ,$
$⟨ x 1 n 1 x 2 n 2 ⟩ = lim k 1 → 0 ⁡ lim k 2 → 0 ⁡ ( − i ∂ k 1 ) n 1 ( − i ∂ k 2 ) n 2 C ( K ) ,$
$P ( X ) = exp ⁡ ( − X t M − 1 X / 2 ) [ ( 2 π ) n det ( M ) ] 1 / 2 .$
$δ X ( t , τ ) = X ( t + τ ) − X ( t ) = ∫ t t + τ R ( s ) d s .$
$A 2 ( ω ) = ∫ − ∞ ∞ A 2 ( t ) e i ω t d t ,$
$A 2 ( t ) = ∫ − ∞ ∞ A 2 ( ω ) e − i ω t d ω / 2 π .$
$C 2 ( τ ) = ∫ − ∞ ∞ A 2 ∗ ( t ) B 2 ( t + τ ) d t ,$
$C 2 ( ω ) = A 2 ∗ ( ω ) B 2 ( ω ) ,$
$∫ − ∞ ∞ | A 2 ( t ) | 2 d t = ∫ − ∞ ∞ | A 2 ( ω ) | 2 d ω / 2 π ,$
$D 2 ( τ ) = ∫ − ∞ ∞ A 2 ( t ) B 2 ( τ − t ) d t ,$
$D 2 ( ω ) = A 2 ( ω ) B 2 ( ω ) .$
$A T ( ω ) = ∫ 0 T A T ( t ) e i ω t d t$
$C T ( τ ) = ∫ 0 T A T ∗ ( t ) B T ( t + τ ) d t .$
$C T ( ω ) = ∫ − T T C T ( τ ) e i ω τ d τ .$
$C T ( ω ) = ∫ 0 T ∫ − T T A T ∗ ( t ) e − i ω t B T ( t + τ ) e i ω ( t + τ ) d τ d t = ∫ 0 T A T ∗ ( t ) e − i ω t ∫ − t T − t B T ( t + τ ) e i ω ( t + τ ) d τ d t = ∫ 0 T A T ∗ ( t ) e − i ω t d t ∫ 0 T B T ( t ′ ) e i ω t ′ d t ′ = A T ∗ ( ω ) B T ( ω ) .$
$C T ( τ ) = ∫ − ∞ ∞ ∫ 0 T ∫ 0 T A T ∗ ( t ) B T ( t ′ ) e i ω ( t ′ − t − τ ) d t ′ d t d ω / 2 π = ∫ 0 T ∫ 0 T A T ∗ ( t ) B T ( t ′ ) δ ( t ′ − t − τ ) d t ′ d t = ∫ 0 T A T ∗ ( t ) B T ( t + τ ) d t .$
$D T ( τ ) = ∫ 0 T A T ( t ) B T ( τ − t ) d t .$
$D T ( ω ) = ∫ 0 2 T D T ( τ ) e i ω τ d τ .$
$D T ( ω ) = ∫ 0 T ∫ 0 2 T A T ( t ) e i ω t B T ( τ − t ) e i ω ( τ − t ) d τ d t , = ∫ 0 T A T ( t ) e i ω t ∫ t T + t B T ( τ − t ) e i ω ( τ − t ) d τ d t = ∫ 0 T A T ( t ) e i ω t d t ∫ 0 T B T ( t ′ ) e i ω t ′ d t ′ = A T ( ω ) B T ( ω ) .$
$D T ( τ ) = ∫ − ∞ ∞ ∫ 0 T ∫ 0 T A T ( t ) B T ( t ′ ) e i ω ( t ′ + t − τ ) d t ′ d t d ω / 2 π = ∫ 0 T ∫ 0 T A T ( t ) B T ( t ′ ) δ ( t ′ + t − τ ) d t ′ d t = ∫ 0 T A T ( t ) B T ( τ − t ) d t .$
$A ( t ) = ∫ 0 t G ( t − s ) R ( s ) d s ,$
$A T ( ω ) = ∫ 0 T ∫ 0 t G ( t − s ) R ( s ) e i ω t d s d t = ∫ 0 T ∫ s T G ( t − s ) R ( s ) e i ω ( t − s + s ) d t d s = ∫ 0 T R ( s ) e i ω s ∫ 0 T − s G ( τ ) e i ω τ d τ d s ,$
$G T ( ω ) = ( e i ω T − 1 ) / i ω ,$
$⟨ | A T ( ω ) | 2 ⟩ = 2 σ T [ 1 − cos ⁡ ( ω T ) ] / ω 2 ,$
$C ( ω , ω ′ ) = ∫ − ∞ ∞ ∫ − ∞ ∞ C ( t , t ′ ) e i ω ′ t ′ − i ω t d t d t ′ ,$
$C ( t , t ′ ) = ∫ − ∞ ∞ ∫ − ∞ ∞ C ( ω , ω ′ ) e i ω t − i ω ′ t ′ d ω d ω ′ / ( 2 π ) 2 .$
$∫ − ∞ ∞ e ± i ω t d ω / 2 π = δ ( t ) ,$
$∫ − ∞ ∞ e ± i ω t d t = 2 π δ ( ω ) .$
$C ( t , t ′ ) = S ( t ) δ ( t − t ′ ) ,$
$C ( ω , ω ′ ) = S ( ω ′ − ω ) ,$
$C ( ω , ω ′ ) = 2 π S ( ω ) δ ( ω − ω ′ ) ,$
$C ( t , t ′ ) = S ( t ′ − t ) ,$
$C T ( ω , ω ′ ) = ∫ 0 T ∫ 0 T C T ( t , t ′ ) e i ω ′ t ′ − i ω t d t d t ′ ,$
$∫ 0 T e ± i ω t d t = ( e i ω T − 1 ) / i ω = 2 π δ T ( ω ) ,$
$C T ( ω , ω ′ ) = ∫ 0 T S T ( t ) e i ( ω ′ − ω ) t d t = S T ( ω ′ − ω ) .$
$C T ( ω , ω ′ ) = ∫ 0 T ∫ 0 T S T ( t ′ − t ) e i ω ′ ( t ′ − t ) e i ( ω ′ − ω ) t d t d t ′ = ∫ 0 T e i ( ω ′ − ω ) t ∫ − t T − t S T ( τ ) e i ω ′ τ d τ d t = ∫ − T T S T ( τ ) e i ω ′ τ ∫ max ( 0 , − τ ) min ( T − τ , T ) e i ( ω ′ − ω ) t d t d τ ∼ 2 π S T ( ω ′ ) δ T ( ω ′ − ω ) .$
$∫ − ∞ ∞ ⟨ | A T ( ω ) | 2 ⟩ e − i ω τ d ω / 2 π = ∫ − ∞ ∞ ∫ 0 T ∫ 0 T ⟨ A T ∗ ( t ) A T ( t ′ ) ⟩ e i ω ( t ′ − t − τ ) d t d t ′ d ω / 2 π = ∫ 0 T ∫ 0 T ⟨ A T ∗ ( t ) A T ( t ′ ) ⟩ δ ( t ′ − t − τ ) d t d t ′ = ∫ max ( 0 , − τ ) min ( T − τ , T ) ⟨ A T ∗ ( t ) A T ( t + τ ) ⟩ d t ,$
$S ( ω ) ∼ T ∫ − 1 / ν 1 / ν C ( τ ) e i ω τ d τ .$
$S ( ω ) / T ∼ ∫ − T T C ( τ ) e i ω τ d τ .$
$G ( t ) = H ( t ) e − ν t ,$
$H ( t ) = { 0 i f t < 0 , 1 / 2 i f t = 0 , 1 i f t > 0.$
$G ( ω ) = 1 / ( ν − i ω ) = ( ν + i ω ) / ( ν 2 + ω 2 ) .$
$G ( t ) = H ( t ) .$
$G α ( t ) = lim α → 0 ⁡ H ( t ) e − α t .$
$G α ( ω ) = lim α → 0 ⁡ ( α + i ω ) / ( α 2 + ω 2 ) .$
$G α ( ω ) → π δ ( ω ) + i / ω .$
$∫ − ∞ ∞ [ i cos ⁡ ( ω t ) + sin ⁡ ( ω t ) ] d ω / 2 π ω .$
$Σ ( t ) = { − 1 i f t < 0 , 0 i f t = 0 , 1 i f t > 0.$
$Σ α ( t ) = lim α → 0 ⁡ Σ ( t ) e − α | t | .$
$Σ α ( ω ) = lim α → 0 ⁡ [ 1 / ( α − i ω ) − 1 / ( α + i ω ) ] = lim α → 0 ⁡ 2 i ω / ( α 2 + ω 2 ) .$
$Σ α ( ω ) → 2 i / ω .$
$H T ( ω ) = ( e i ω T − 1 ) / i ω = sin ⁡ ( ω T ) / ω + i [ 1 − cos ⁡ ( ω T ) ] / ω .$
$H T ( ω ) ≈ π δ ( ω ) + i / ω .$
$Σ T ( ω ) = ( e i ω T − 1 ) / i ω − ( e − i ω T − 1 ) / ( − i ω ) = 2 i [ 1 − cos ⁡ ( ω T ) ] / ω .$
$Σ T ( ω ) ≈ 2 i / ω .$
$∫ − T t ∫ − ∞ ∞ ( 1 − e i ω T ) e − i ω s d ω d s / 2 π = ∫ − T t [ δ ( s ) − δ ( s − T ) ] d s = H T ( t ) .$
$∫ − T t ∫ − ∞ ∞ ( 2 − e i ω T − e − i ω T ) e − i ω s d ω d s / 2 π = ∫ − T t [ − δ ( s + T ) + 2 δ ( s ) − δ ( s − T ) ] d s = Σ T ( t ) .$
$F ( ω ) = i π P ∫ − ∞ ∞ F ( ω ′ ) ω − ω ′ d ω ′ ,$
$F ( ω ) = ∫ 0 u ∫ − ∞ ∞ F ( ω ′ ) e − i ω ′ t ( d ω ′ / 2 π ) e i ω t d t = ∫ − ∞ ∞ F ( ω ′ ) ∫ 0 u e i ( ω − ω ′ ) t ( d t / 2 π ) d ω ′ .$
$∫ 0 u e i ( ω − ω ′ ) t d t / 2 π = δ u ( ω − ω ′ ) .$
$δ ∞ ( ω ) = δ ( ω ) / 2 + i / 2 π ω ,$
$δ T ( ω ) = sin ⁡ ( ω T ) / 2 π ω + i [ 1 − cos ⁡ ( ω T ) ] / 2 π ω .$
$F ( ω ) = ∫ − ∞ ∞ [ δ ( ω − ω ′ ) + i / π ( ω − ω ′ ) ] F ( ω ′ ) d ω ′ / 2 = [ F ( ω ) + F ( ω ) ] / 2.$
$A ( t ) = ∫ 0 t R ( s ) d s ,$
$⟨ | A ( t ) | 2 ⟩ = ∫ 0 t ∫ 0 t ⟨ R ∗ ( s ) R ( s ′ ) ⟩ d s d s ′ = σ t .$
$− i ω A ( ω ) = R ( ω ) − A ( T ) e i ω T ,$
$ω 2 ⟨ | A ( ω ) | 2 ⟩ = ⟨ | R ( ω ) | 2 ⟩ − 2 R e [ ⟨ R ∗ ( ω ) A ( T ) e i ω T ⟩ ] + ⟨ | A ( T ) | 2 ⟩ .$
$⟨ R ∗ ( ω ) A ( T ) e i ω T ⟩ = ∫ 0 T ∫ 0 T ⟨ R ∗ ( s ) R ( s ′ ) ⟩ e i ω ( T − s ) d s d s ′ = σ ( e i ω T − 1 ) / i ω .$
$ω 2 ⟨ | A ( ω ) | 2 ⟩ = σ T − 2 σ sin ⁡ ( ω T ) / ω + σ T = 2 σ T [ 1 − s i n c ( ω T ) ] .$
$− i ω A ( ω ) = R ( ω ) − R ( 0 ) e i ω T ,$
$⟨ R ∗ ( t ) R ( t ′ ) ⟩ = σ δ ( t − t ′ ) .$
$⟨ R ∗ ( ω ) R ( ω ′ ) ⟩ = ∫ 0 T ∫ 0 T ⟨ R ∗ ( t ) R ( t ′ ) ⟩ e i ( ω ′ t ′ − ω t ) d t d t ′ = σ ∫ 0 T e i ( ω ′ − ω ) t d t = σ [ e i ( ω ′ − ω ) T − 1 ] / i ( ω ′ − ω ) = 2 π σ δ T ( ω ′ − ω ) ,$
$C ( ω , ω ′ ) ∼ 2 π σ G ∗ ( ω ) G ( ω ′ ) δ T ( ω ′ − ω ) .$
$⟨ R ∗ ( t ) R ( t ′ ) ⟩ = ( σ ν / 2 ) e − ν | t − t ′ | .$
$⟨ R ∗ ( ω ) R ( ω ′ ) ⟩ = ( σ ν / 2 ) ∫ 0 T ∫ 0 T e − ν | t − t ′ | e i ( ω ′ t ′ − ω t ) d t d t ′ .$
$e − ( ν + i ω ) t ∫ 0 t e ( ν + i ω ′ ) t ′ d t ′ + e ( ν − i ω ) t ∫ t T e − ( ν − i ω ′ ) t ′ d t ′ = e i ( ω ′ − ω ) t − e − ( ν + i ω ) t ( ν + i ω ′ ) + e i ( ω ′ − ω ) t − e − ( ν − i ω ′ ) T e ( ν − i ω ) t ( ν − i ω ′ ) .$
$e i ( ω ′ − ω ) T − 1 i ( ω ′ − ω ) ( ν + i ω ′ ) − 1 − e − ( ν + i ω ) T ( ν + i ω ) ( ν + i ω ′ )$
$e i ( ω ′ − ω ) T − 1 i ( ω ′ − ω ) ( ν − i ω ′ ) − e i ( ω ′ − ω ) T − e − ( ν − i ω ′ ) T ( ν − i ω ) ( ν − i ω ′ ) .$
$e i ( ω ′ − ω ) T − 1 i ( ω ′ − ω ) [ 1 ν + i ω ′ + 1 ν − i ω ′ − i ( ω ′ − ω ) ( ν − i ω ) ( ν − i ω ′ ) ] = e i ( ω ′ − ω ) T − 1 i ( ω ′ − ω ) 2 ν + i ( ω ′ − ω ) ( ν + i ω ′ ) ( ν − i ω ) .$
$e − ( ν + i ω ) T − 1 ( ν + i ω ) ( ν + i ω ′ ) + e − ( ν − i ω ′ ) T − 1 ( ν − i ω ) ( ν − i ω ′ ) .$
$− 2 ( ν 2 − ω ω ′ ) ( ν 2 + ω 2 ) [ ν 2 + ( ω ′ ) 2 ] .$
$C ( ω , ω ′ ) = σ ν 2 { e i ( ω ′ − ω ) T − 1 i ( ω ′ − ω ) 2 ν + i ( ω ′ − ω ) ( ν + i ω ′ ) ( ν − i ω ) − 2 ( ν 2 − ω ω ′ ) ( ν 2 + ω 2 ) [ ν 2 + ( ω ′ ) 2 ] } .$
$C ( ω , ω ′ ) = σ ∫ 0 T e − i ω T − e − i ω s − i ω e i ω ′ T − e i ω ′ s i ω ′ d s = ( σ / ω ω ′ ) ∫ 0 T { e i ( ω ′ − ω ) T − e i ω ′ s − i ω T − e i ω ′ T − i ω s + e i ( ω ′ − ω ) s } d s .$
$T e i ( ω ′ − ω ) T + e i ( ω ′ − ω ) T − 1 i ( ω ′ − ω ) ,$
$e − i ω T ( 1 − e i ω ′ T ) i ω ′ + e i ω ′ T ( e − i ω T − 1 ) i ω .$
$C ( ω , ω ′ ) = σ ω ω ′ [ T e i ( ω ′ − ω ) T + e i ( ω ′ − ω ) T − 1 i ( ω ′ − ω ) + e − i ω T ( 1 − e i ω ′ T ) i ω ′ + e i ω ′ T ( e − i ω T − 1 ) i ω ] .$