Abstract
We theoretically and numerically study the effects of thermal noise on pulses in backwards stimulated Brillouin scattering (SBS). Using a combination of stochastic calculus and numerical methods, we derive a theoretical model that can be used to quantitatively predict noise measurements. We study how the optical pulse configuration, including the input powers of the pump and Stokes fields, pulse durations and interaction time, affects the noise in the output Stokes field. We investigate the effects on the noise of the optical loss and waveguide length, and we find that the signal-to-noise ratio can be significantly improved, or reduced, for specific combinations of waveguide properties and pulse parameters.
© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
Backward stimulated Brillouin scattering (SBS) is an opto-acoustic process that results from the coherent interaction between two counter-propagating optical fields and an acoustic wave inside an optical waveguide [1–5]. The coupling between these fields occurs either through electrostriction or radiation pressure [6,7], in which the electromagnetic fields induce very small strains in the material, which in turn stimulate the generation of coherent phonons via photoelasticity or via physical motion of the waveguide boundaries [6–8]. This interaction has been used in various applications, including narrow-band RF and optical signal filtering [9,10], phase conjugation and precision spectroscopy [1], novel laser sources [11,12], and in recent experiments in opto-acoustic pulse storage and memory [13]. An intrinsic property of SBS is noise: the interaction between the pump and the thermal phonon background of the waveguide generates spontaneous Stokes photons, which combine with the input signal to produce a noisy output signal. SBS noise was first studied by Boyd and Gaeta [14,15] and Ferreira et al. [16]. These studies assumed that the optical and acoustic fields are all weakly guided, and furthermore examined the effects of noise in the quasi-continuous wave regime. Later work [17] extended the distributed fluctuating source model from Boyd and Gaeta to include the effect of optical loss, with the addition of an acousto-optic effective area using scalar modal fields, and was again restricted to quasi-CW pulses and weak guidance. More recent studies have analyzed the effects of thermal noise in the context of SBS lasers and micro-ring resonators [18–21]. A dynamic, fully-vectorial treatment that can be applied to predict noise in modern integrated SBS waveguide experiments [2,22] is important for the further development of SBS signal-processing applications.
In this paper, we derive a theoretical model for SBS noise for pulses in integrated optical waveguides with high-contrast materials. We focus on amplitude noise arising from the thermal phonon field, because in backwards SBS experiments this process has been observed to be significantly stronger than phase noise [23,24]. We explore the dependence of the noise on the waveguide and pulse properties. We find that in the regime where the pump is undepleted by the SBS process, it is possible to derive analytic results for the absolute level of noise, as well as for the signal-to-noise ratio, even for time-varying pulses of arbitrary shape. Since numerical computation of ensemble averages is expensive, an analytic model that computes the SBS noise properties of pulses is extremely useful for the design of SBS-active devices. We use this analytic model to explore the impact of pulse and waveguide properties on the noise and optical signal-to-noise ratio (OSNR), and discuss the implications of these results in SBS experiments, where it is often advantageous to minimize the Stokes noise.
2. Theory and formalism
We consider backward SBS interactions in a waveguide of finite length $L$ oriented along the $z$-axis (see Fig. 1): a pump pulse with angular frequency $\omega _1$ is injected into the waveguide at $z=0$ and propagates in the positive $z$-direction, while a signal pulse is injected at $z=L$ and propagates in the negative $z$-direction. The signal pulse has frequency centred around the Brillouin Stokes frequency $\omega _2=\omega _1-\Omega$, which is down-shifted from the pump by the Brillouin shift $\Omega$, and has a spectral width that lies entirely within the Brillouin linewidth. When these two pulses meet, the pump transfers energy to the signal field, resulting in coherent amplification of the signal around the Brillouin frequency. At the same time, as the pump moves through the waveguide, it interacts with the thermal phonon field and generates spontaneous Stokes photons which also propagate in the negative $z$-direction. This spontaneous Stokes field combines with the coherent signal to form a noisy amplified output field centered around the Stokes frequency.
To compute the noise transferred to the signal from the thermal field, we first solve the coupled equations of the pump $\mathbf {E}_1$, Stokes $\mathbf {E}_2$ and acoustic $\mathbf {U}$ fields. We write these fields using the ansatz
The waveguide has optical loss $\alpha$ and acoustic attenuation $\alpha _{\textrm {ac}}$ (in units of m$^{-1}$), with optical group velocity $v$ and acoustic group velocity $v_a$. The coefficients $Q_{1,2,a}$ represent the coupling strength of the SBS interaction, and can be calculated from the optical and acoustic modes of the waveguide [25]. From conservation of energy, we have $Q_2 = Q_1^{*}$ and $Q_a=Q_1$ [7]. These coupling constants contain all the information about the vectorial modal fields, including polarization.
The function $R(z,t)$ represents a complex-valued space-time white noise with zero mean $\left \langle R(z,t)\right \rangle = 0$ and auto-correlation function $\left \langle R(z,t) R^{*}(z',t') \right \rangle = \delta (z-z')\delta (t-t')$, where $\delta$ is the Dirac-delta function. This type of delta correlation arises in the standard treatment of noise in optical fibers [26], and arises from the effective continuum of phonon modes supported by the waveguide and coupled through dissipation to the coherent excitation of the elastic field. The strength of the noise $\sigma$ gives the conversion of coherent elastic energy in the medium into thermal phonons, in the absence of any nonlinear optical effects [14], and is given by (see Appendix A)
where $k_B T$ is the thermal energy per phonon mode. This quantity is directly proportional to the acoustic loss in the medium $\alpha _{\textrm {ac}}$ [27] in accordance with the fluctuation-dissipation theorem [28].2.1 Noise under the undepleted pump approximation
An analytic solution to Eqs. (4)–(6) exists under the undepleted pump approximation, which is valid for the case in which the signal power $P_2$ remains much smaller than the pump power $P_1$. We also assume that the depletion of the pump by the thermal field is negligible compared to optical loss [16]. The equation for the pump envelope is therefore
We compute the total, noisy, Stokes power by considering the Stokes envelope as a sum of two contributions: a coherent signal field, which is the Stokes in the absence of any thermal noise, and a spontaneous Stokes field which contains all the contributions from the background thermal noise. This is expressed by separating the two terms in Eq. (17) as $a_2(z,t) = a_2^{\textrm {sig}}(z,t) + a_2^{\textrm {spo}}(z,t)$ where
This product of integrals can be cast as a double-integral by applying the stochastic Fubini theorem [32]. A sufficient condition for this to hold is that the integrand in (20) be square-integrable and finite [33]. This condition holds since $D(z,t)$ reaches a steady-state after a short time $t$ [34] and both $a_1$ and $G_N$ are bounded functions [30]. The spontaneous power becomes
The relation $\left \langle a_1(z,t) D(z,t)\right \rangle = \langle a_1(z,t)\rangle \langle D(z,t)\rangle = 0$ holds since the laser noise in $a_1$ and the thermal noise in $D$ are not coupled to each other, so they are statistically independent. Therefore, the ensemble-averaged spontaneous Stokes is
Once this spontaneous power has been calculated, a useful measure to quantify the integrity of the Stokes signal detected at the waveguide output is the OSNR. This is computed by taking the ratio of the average coherent signal power to spontaneous noise power over a specified bandwidth [35]. In time-domain, the OSNR may be calculated over an arbitrary but suitably long time period $T_0$ using the time averaged powers
3. Results for a chalcogenide SBS chip
Figure 2 shows the ensemble-averaged signal power and spontaneous Stokes power computed using Eqs. (22) and (25). These calculations assume a dielectric waveguide at temperature 300 K of length 23.7 cm and SBS gain coefficient $g_0 = 423$ m$^{-1}$W$^{-1}$, which corresponds to the SBS chip structure used in [36] at zero detuning ($\Delta = 0$). We consider a Gaussian pump pulse with peak power 110 mW and FWHM of 200 ps, and a Gaussian signal pulse with peak power 10 nW and FWHM of 200 ps. The optical pump frequency is fixed at 200 THz, while the acoustic frequency is set to 7.6 GHz, with a Brillouin linewidth of $\Delta \nu _B=30$ MHz. The linear optical loss is 0.05 dB/cm. Equations (4)–(6) are solved numerically with a split-step integration scheme that calculates the drift of the pulses first, and then uses an Euler-Mayurama step for computing the noisy nonlinear interaction between the pump and Stokes.
For these simulations, the numerical grid consists of 600 points in $z$ and 900 points in $t$. These values were chosen, using convergence tests in the absence of noise, to ensure an error in the gain of under $1\%$. We see in Fig. 2(a) that as the pump propagates through the waveguide, a noisy Stokes field is generated that trails the leading edge of the pump pulse. This field combines with the input signal, with the resulting noisy Stokes field further amplified by the pump. Figure 2(b) shows the total Stokes power as a function of $z$ at $t=1.8$ ns, where the dashed black line represents the analytic solution to Eq. (25) added to the coherent Stokes power in Eq. (22), and the blue line is a numerical ensemble average over 20 independent runs of the numerical method. Equation (25) contributes approximately 3 nW of Stokes power across the length of the waveguide, while Eq. (22) contributes 10 nW at $z=23.7$ cm which is amplified as it moves towards $z=0$.
We now use Eq. (25) to study the effect on the noise of the tunable properties of the pump such as peak power and pulse duration, and how this may interact with fixed parameters such as waveguide length and optical loss. We consider a Gaussian signal pulse with FWHM of 200 ps, and a rectangular pump pulse of width $\tau _p$ that varies between 10 ps$-$10 ns. The total pulse energy of the input pump field is kept constant, while the peak power $P_{\textrm {p0}}$ varies with the pulse width, ranging from $P_{\textrm {p0}}= 5$ mW ($\tau _p=10$ ns) to $P_{\textrm {p0}}= 5$ W ($\tau _p=10$ ps). In all the calculations, the output Stokes field at $z=0$ is averaged over a time period of 800 ps, which is four times the signal FWHM and contains more than 99.7% of the pulse’s energy. The results for the spontaneous Stokes power, signal power and OSNR are shown in Fig. 3. We observe in Fig. 3(a) and (b) that the Stokes noise and signal powers remain constant for pump widths $\tau _p \lesssim 3$ ns before decreasing. In this low optical loss case, the near-constant values of the Stokes noise and signal power can be attributed to the constant energy in the pump pulse, as shown in Fig. 4: for short pulse-lengths, the signal and the pump interact entirely within the waveguide. As the pump pulse becomes longer than the waveguide and its peak power decreases, the Stokes field is able to interact with less of its energy inside the medium, resulting in less amplification of the Stokes. The signal and spontaneous fields fall off at different rates with respect to $\tau _p$ because the signal is proportional to the net gain squared $|G_N|^{2}$ via Eq. (22), while the spontaneous Stokes power is proportional to the pump power $P_1$ multiplied by $|G_N|^{2}$ via Eq. (25). The gain $|G_N|^{2}$ decreases with larger $\tau _p$ because more of the pump pulse lies outside of the interaction region, and at the same time $P_1$ decreases while its peak power is lowered for longer $\tau _p$, causing the spontaneous Stokes to decrease more rapidly than the signal. This results in an OSNR figure which increases with $\tau _p$, as shown in Fig. 3(c).
For the high optical loss case shown in Fig. 3(d)$-$(f), we observe a maximum in both the Stokes noise and signal powers for a critical value of $\tau _p$, before decreasing with longer pump lengths as in the lower loss case. This results in a minimum in the OSNR; this minimum occurs because the optical loss is dominant in the regions inside the waveguide that are not covered by the pump (see Fig. 5). As the pump becomes wider, the pump coverage inside the waveguide increases, compensating for the losses incurred at the uncovered regions. This results in higher amplification of the Stokes despite the peak pump power decreasing with larger $\tau _p$. However, this increased amplification vanishes once the pump pulse becomes longer than the waveguide, in which case the Stokes power decreases in a similar manner to the low-loss case.
The emergence of the OSNR minimum in Fig. 3(f) is a function of the optical loss $\alpha$, the pump pulse width $\tau _p$ and the waveguide length $L$. This is shown as a contour plot in Fig. 6; we see that for smaller losses there is no maximum peak in the spontaneous noise power, and therefore no minimum in the OSNR. Below this value the Stokes power and OSNR follow the behavior of Fig. 3(a)$-$(c), staying approximately constant for short pump pulses, then decreasing and increasing respectively for higher values of $\tau _p$. However, once the loss reaches a critical value, the noise power develops a maximum at approximately the point in which the interaction time $\tau _i$ between pump and signal is equal to the transit time $L/v$. This corresponds to a minimum in the OSNR, because making $\tau _p$ longer than the transit time results in a reduced level of noise, since the spontaneous Stokes decreases more rapidly than the amplified signal. The value $\tau _{\textrm {max}}$ that maximizes the noise corresponds to the pump-length when the interaction time $\tau _i$ satisfies $\tau _i = L/v$ and can be computed as follows: first, we compute an effective pulse-width $\tau _{\textrm {eff}}$ using the equation
This noise behavior is not restricted to rectangular pump pulses. In Fig. 6(c) and (d), we observe the same behavior when using a Gaussian pump pulse, and again the same maximum in the Stokes noise can be observed. Here, some of the pump’s energy is more localized near the central peak, while the long tails of the pulse contain a very small portion of the energy. This limits the region inside the waveguide in which the amplification from the pump is able to compensate for the optical losses, and results in a maximum Stokes noise peak that occurs at a smaller $\tau _p$ compared to the rectangular pump case.
Finally, we compute the effects of the noise and OSNR in the case where the pump pulse has a fixed peak power. Both the spontaneous Stokes noise and signal power increase with $\tau _p$, reaching a steady state for $\tau _p$ longer than the transit time of the waveguide. This results in the OSNR decreasing until reaching a constant value for larger $\tau _p$, as shown in Fig. 7(b). This occurs due to the pump pulse covering more of the waveguide as $\tau _p$ gets larger: initially there is more energy injected into the SBS medium, but as soon as the pump pulse becomes comparable in size to the waveguide; the amount of pump energy that interacts with the Stokes field reaches a constant value. This then leads to a steady-state for both the signal and spontaneous Stokes powers.
4. Validity of the analytic model
We now estimate the region of validity of our analytic model, by considering an SBS chalcogenide chip with length 23.7 cm, with a Gaussian pump pulse of $P_{\textrm {p0}} = 300$ mW and fixed optical loss of $\alpha = 0.3$ dB/cm. The energy in both the pump and signal input fields is given by the expressions $E_{1} = \int _{-\infty }^{\infty } P_1(0,t) \textrm {d}{t}$ and $\quad E_{2} = \int _{-\infty }^{\infty } P_2(L,t) \textrm {d}{t}$ respectively. For each ratio of energies $E_2/E_1$ where $E_2 \leq E_1$ we compute the relative error in the OSNR using both the analytic model in Eq. (25) and the numerical solver for Eqs. (4)–(6).
In Fig. 8(a) we see that the ONSR error remains below 2% when the signal energy relative to the pump is small, but increases suddenly for larger signal powers. This occurs because of pump depletion, which was neglected in Eq. (25) in deriving the analytic model: as the pump becomes depleted, some of its energy gets transferred to the spontaneous Stokes as well as to the signal, enhancing the coupling between the fields. This results in some of the spontaneous Stokes energy feeding back into the pump, making the multiplicative noise more dominant than the additive noise in Eq. (25). This means that as the pump becomes depleted, there is less spontaneous Stokes generated from the thermal background than coherent Stokes amplified from the input signal, even though both fields increase as the pump power increases. This is in line with previous studies [37,38] which considered variance in the reflected Stokes beam in the absence of an input signal (i.e. for purely spontaneous SBS): these studies found that the relative strength of the spontaneous Stokes field was reduced compared to the coherent field once the pump power began to be depleted. For the analytic model, the departure from the undepleted regime results in an error in the OSNR, as shown in Fig. 8(a). The analytic model is therefore expected to yield accurate results provided the pump remains undepleted. To extend the study of the noise dynamics of SBS interactions beyond the undepleted pump regime, and to study other features such as phase noise, it is necessary to solve Eqs. (4)–(6) numerically. This can be done using a split-step integration scheme, in which the propagation of the optical fields is performed first along $z$, followed by an integration in time at each point $z$ of the nonlinear interaction between $a_1$ and $a_2$ and their coupling to the thermal noise $D(z,t)$. A full study of this numerical method is beyond the scope of this paper.
5. Conclusion
In this paper, we have developed a general mathematical model for solving the dynamical SBS equations in the presence of thermal noise. We present a simple analytic model in the undepleted pump regime whereby the amplitude noise in the Stokes field can be computed for pulses of arbitrary shape and size, which is useful in Brillouin signal processing. The analytic results show that in the case of a constant energy pump field and lossy media, the OSNR has a minimum peak near the region in which the interaction time of the pulses matches the transit time of the signal in the waveguide. This occurs as the spontaneous Stokes noise increases towards a maximum value, which results from longer pump pulses compensating for the linear optical losses in the medium. Once the pump pulse becomes longer than the waveguide, the loss dominates again as less pump energy fits inside the SBS medium. This behaviour of the OSNR is mediated by the pump energy, pump shape, waveguide length and optical loss, making it important to choose the right parameter combination to maximize the OSNR in a specific device.
Appendix A. Derivation of the noise strength
In this Appendix, we find an expression for the noise strength parameter $\sigma$ in Eq. (7) by considering the acoustic envelope field $b_j$ at a set of discrete positions along the waveguide $z_j$ separated by a small width $\Delta z$, such that $b_j=b(z_j)$, and assuming $\partial _z b$ is negligible due to the small propagation length of the phonons [4]. The thermal noise field in each slice is $\tilde {R}_j$, and we assume it to be statistically independent between any two distinct slices. Additionally, we assume a weak opto-acoustic coupling so that the $b_j$ field is mainly dominated by the thermal noise. This leads to the following equation for the acoustic envelope field in the $j$th slice in the waveguide
Appendix B. Full derivation of the solution to the Stokes field in the undepleted pump approximation
In this Appendix we show the complete derivation of the analytic solution for the Stokes envelope field in Eq. (17). Starting with Eq. (15) and by using the transformation of variables $\xi = z$ and $\tau = t + z/v$, one moves into a frame of reference which moves along with the Stokes field. The partial derivatives in the new coordinate system $(\xi ,\tau )$ are thus
Funding
Australian Research Council (DP160101691, DP200101893); Macquarie University (MQRF0001036); University of Technology Sydney (00099F).
Acknowledgments
The authors acknowledge funding from the Australian Research Council (ARC) (Discovery Projects DP160101691, DP200101893), the Macquarie University Research Fellowship Scheme (MQRF0001036) and the UTS Australian Government Research Training Program Scholarship (00099F). Part of the numerical calculations were performed on the UTS Interactive High Performance Computing (iHPC) facility.
Disclosures
The authors declare no conflicts of interest.
References
1. B. J. Eggleton, C. G. Poulton, P. T. Rakich, M. J. Steel, and G. Bahl, “Brillouin integrated photonics,” Nat. Photonics 13(10), 664–677 (2019). [CrossRef]
2. R. Pant, D. Marpaung, I. V. Kabakova, B. Morrison, C. G. Poulton, and B. J. Eggleton, “On-chip stimulated Brillouin scattering for microwave signal processing and generation,” Laser Photonics Rev. 8(5), 653–666 (2014). [CrossRef]
3. L. Brillouin, “Diffusion de la lumière et des rayons X par un corps transparent homogène,” Ann. Phys. 9(17), 88–122 (1922). [CrossRef]
4. R. W. Boyd, Nonlinear optics (Elsevier, 2003).
5. A. Kobyakov, M. Sauer, and D. Chowdhury, “Stimulated Brillouin scattering in optical fibers,” Adv. Opt. Photonics 2(1), 1–59 (2010). [CrossRef]
6. P. T. Rakich, C. Reinke, R. Camacho, P. Davids, and Z. Wang, “Giant enhancement of stimulated Brillouin scattering in the subwavelength limit,” Phys. Rev. X 2(1), 011008 (2012). [CrossRef]
7. C. Wolff, M. J. Steel, B. J. Eggleton, and C. G. Poulton, “Stimulated Brillouin scattering in integrated photonic waveguides: Forces, scattering mechanisms, and coupled-mode analysis,” Phys. Rev. A 92(1), 013836 (2015). [CrossRef]
8. J. E. Sipe and M. J. Steel, “A Hamiltonian treatment of stimulated Brillouin scattering in nanoscale integrated waveguides,” New J. Phys. 18(4), 045004 (2016). [CrossRef]
9. A. Choudhary, Y. Liu, D. Marpaung, and B. J. Eggleton, “On-chip Brillouin filtering of RF and optical signals,” IEEE J. Sel. Top. Quantum Electron. 24(6), 1–11 (2018). [CrossRef]
10. H. Jiang, L. Yan, W. Pan, B. Luo, and X. Zou, “Ultra-high speed RF filtering switch based on stimulated Brillouin scattering,” Opt. Lett. 43(2), 279–282 (2018). [CrossRef]
11. B. Steinhausser, A. Brignon, E. Lallier, J.-P. Huignard, and P. Georges, “High energy, single-mode, narrow-linewidth fiber laser source using stimulated Brillouin scattering beam cleanup,” Opt. Express 15(10), 6464–6469 (2007). [CrossRef]
12. W. Loh, A. A. Green, F. N. Baynes, D. C. Cole, F. J. Quinlan, H. Lee, K. J. Vahala, S. B. Papp, and S. A. Diddams, “Dual-microcavity narrow-linewidth Brillouin laser,” Optica 2(3), 225–232 (2015). [CrossRef]
13. M. Merklein, B. Stiller, K. Vu, S. J. Madden, and B. J. Eggleton, “A chip-integrated coherent photonic-phononic memory,” Nat. Commun. 8(1), 574–577 (2017). [CrossRef]
14. R. W. Boyd, K. Rzazewski, and P. Narum, “Noise initiation of stimulated Brillouin scattering,” Phys. Rev. A 42(9), 5514–5521 (1990). [CrossRef]
15. A. L. Gaeta and R. W. Boyd, “Stochastic dynamics of stimulated Brillouin scattering in an optical fiber,” Phys. Rev. A 44(5), 3205–3209 (1991). [CrossRef]
16. M. Ferreira, J. Rocha, and J. Pinto, “Analysis of the gain and noise characteristics of fiber Brillouin amplifiers,” Opt. Quantum Electron. 26(1), 35–44 (1994). [CrossRef]
17. A. Kobyakov, S. Darmanyan, and D. Chowdhury, “Exact analytical treatment of noise initiation of SBS in the presence of loss,” Opt. Commun. 260(1), 46–49 (2006). [CrossRef]
18. W. Loh, S. B. Papp, and S. A. Diddams, “Noise and dynamics of stimulated Brillouin scattering microresonator lasers,” Phys. Rev. A 91(5), 053843 (2015). [CrossRef]
19. P. Kharel, R. Behunin, W. Renninger, and P. Rakich, “Noise and dynamics in forward Brillouin interactions,” Phys. Rev. A 93(6), 063806 (2016). [CrossRef]
20. R. O. Behunin, N. T. Otterstrom, P. T. Rakich, S. Gundavarapu, and D. J. Blumenthal, “Fundamental noise dynamics in cascaded-order Brillouin lasers,” Phys. Rev. A 98(2), 023832 (2018). [CrossRef]
21. N. T. Otterstrom, R. O. Behunin, E. A. Kittlaus, Z. Wang, and P. T. Rakich, “A silicon Brillouin laser,” Science 360(6393), 1113–1116 (2018). [CrossRef]
22. W. Zhang and R. A. Minasian, “Widely tunable single-passband microwave photonic filter based on stimulated Brillouin scattering,” IEEE Photonics Technol. Lett. 23(23), 1775–1777 (2011). [CrossRef]
23. M. Pelusi, T. Inoue, and S. Namiki, “Noise characterization of Brillouin amplified narrowband carriers for coherent communications,” in 2018 European Conference on Optical Communication (ECOC), (IEEE, 2018), pp. 1–3.
24. A. Zarifi, M. Merklein, Y. Liu, A. Choudhary, B. J. Eggleton, and B. Corcoran, “Broadband sbs filter for optical carrier recovery applications in telecommunication systems,” in 2020 Conference on Lasers and Electro-Optics (CLEO), (IEEE, 2020), pp. 1–2.
25. B. C. Sturmberg, K. B. Dossou, M. J. Smith, B. Morrison, C. G. Poulton, and M. J. Steel, “Finite element analysis of stimulated Brillouin scattering in integrated photonic waveguides,” J. Lightwave Technol. 37(15), 3791–3804 (2019). [CrossRef]
26. P. Drummond and J. F. Corney, “Quantum noise in optical fibers. I. Stochastic equations,” J. Opt. Soc. Am. B 18(2), 139–152 (2001). [CrossRef]
27. H. L. Pécseli, Fluctuations in physical systems (Cambridge University, 2000).
28. J. Lima and A. Plastino, “On the classical energy equipartition theorem,” Braz. J. Phys. 30(1), 176–180 (2000). [CrossRef]
29. R. Kubo, “Statistical-mechanical theory of irreversible processes. i. general theory and simple applications to magnetic and conduction problems,” J. Phys. Soc. Jpn. 12(6), 570–586 (1957). [CrossRef]
30. K. Itô, “On a formula concerning stochastic differentials,” Nagoya Math. J. 3, 55–65 (1951). [CrossRef]
31. P. Drábek and G. Holubová, Elements of partial differential equations (Walter de Gruyter GmbH & Co KG, 2014).
32. J. Van Neerven and M. C. Veraar, “On the stochastic Fubini theorem in infinite dimensions,” in Stochastic Partial Differential Equations and Applications VII, Levico Terme, 2004, in: Lecture Notes Pure Appl. Math., (CRC, 2005), p. 323–336.
33. M. Veraar, “The stochastic Fubini theorem revisited,” Stochastics 84(4), 543–551 (2012). [CrossRef]
34. R. Kubo, “The fluctuation-dissipation theorem,” Rep. Prog. Phys. 29(1), 306255 (1966). [CrossRef]
35. G. P. Agrawal, Fiber-optic communication systems (John Wiley & Sons, 2012), 4th ed.
36. Y. Xie, A. Choudhary, Y. Liu, D. Marpaung, K. Vu, P. Ma, D.-Y. Choi, S. Madden, and B. J. Eggleton, “System-level performance of chip-based Brillouin microwave photonic bandpass filters,” J. Lightwave Technol. 37(20), 5246–5258 (2019). [CrossRef]
37. A. Fotiadi, R. Kiyan, O. Deparis, P. Mégret, and M. Blondel, “Statistical properties of stimulated Brillouin scattering in single-mode optical fibers above threshold,” Opt. Lett. 27(2), 83–85 (2002). [CrossRef]
38. M. Bowers and N. Luzod, “Stimulated Brillouin scattering in optical fibers with end reflections excited by broadband pump waves,” Opt. Eng. 58(10), 1 (2019). [CrossRef]