Abstract
We theoretically and numerically study the efficiency of Brillouin-based opto-acoustic data storage in a photonic waveguide in the presence of thermal noise and laser phase noise. We compare the physics of the noise processes and how they affect different storage techniques, examining both amplitude and phase storage schemes. We investigate the effects of storage time and pulse properties on the quality of the retrieved signal and find that phase storage is less sensitive to thermal noise than amplitude storage.
© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
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]. In backward SBS, a pump field at frequency $\omega _1$ is injected at one end of the waveguide, and a counter-propagating field at the Brillouin Stokes frequency $\omega _2 = \omega _1 - \Omega$, where $\Omega$ is the Brillouin frequency shift, is injected at the other; the interaction of these two fields creates a strong acoustic wave at frequency $\Omega$. When the Stokes power is comparable in size to the pump power, the pump energy can be completely depleted and some of it transferred to the acoustic field, along with its amplitude and phase information. A recent application of this effect is opto-acoustic memory storage (see Fig. 1): an optical data pulse at the pump frequency interacts with an optical “write” pulse (at the Stokes frequency), creating an acoustic hologram where the original data is temporarily stored [6–8]; an optical ”read” pulse at the Stokes frequency can then be used to regenerate the original data pulse. Brillouin-based opto-acoustic information storage has been demonstrated experimentally in fibres [6] as well as more recently in on-chip experiments [7–12], with storage times ranging up to 40 ns [12]. Opto-acoustic storage, however, will inevitably be limited by noise, which exists because of the presence of thermal phonons in the waveguide, as well as being an inherent feature of the lasers used in the data and read/write pulses. At room temperatures, it is known that thermal noise significantly degrades the quality of Brillouin processes [13–18], and it can be expected that it will also place limits on information retrieval efficiency in Brillouin storage experiments, measured as the amount of power that can be retrieved from the acoustic field after a certain time. This has been experimentally measured [6,11], however while noise has been observed as a feature of these experiments, it is not yet clear how noise impacts the accuracy of the information retrieval. Furthermore, most studies (except, notably, [10]) have focused on amplitude encoded storage, in which bits of information are stored in individual pulses. Alternatively, phase encoded storage may offer higher storage efficiency and be less sensitive to noise. A quantitative understanding of the impact of noise on the storage of both amplitude and phase information is needed for the further development of practical Brillouin-based storage devices.
In this paper, we apply our previous theoretical [18] and numerical [19] models to the simulation of Brillouin storage in the presence of thermal noise. These studies focused on the noise properties of optical pulses in the case of SBS amplification, with the pump containing a lot more power than the input Stokes seed. Here, we investigate the effect of noise on both amplitude and phase storage of information, by using a small pump power and a large Stokes seed power. We quantify the retrieval efficiency and accuracy in the form of the packet error rate (PER) of 8-bit sequences. We find that phase storage offers a significant improvement in the duration with which information can be stored without degradation due to thermal noise. We examine this effect in more detail by computing the effect of thermal noise on the amplitude and phase of a phase-encoded signal, and find that although the variance in amplitude and phase increases at the same rates, the phase information is more robust to noise in accordance with the additive-white-Gaussian-noise (AWGN) model of discrete communication theory.
2. Numerical simulation of noise
2.1 Brillouin coupled mode equations
We consider a waveguide of length $L$ oriented along the $z$ coordinate axis, in which an optical data stream propagates in the positive $z$ direction. The interaction between the two optical fields and the acoustic field can be described by the coupled system [18,19]
The thermal noise is represented by a zero-mean space-time white noise function, with auto-correlation function given by $\langle R(z,t)R(z',t')\rangle = \delta (z-z')\delta (t-t')$. The parameter $k_B T\alpha _{\textrm {ac}}$ is the strength of the thermal noise, which arises from the fluctuation-dissipation theorem [18]. The functions describing the data stream and the read-write pulses are given by $a_{\textrm {data}}(t)$, $a_{\textrm {read}}(t)$ and $a_{\textrm {write}}(t)$, which describe the values of the optical envelopes at $z=0$ (data) and $z=L$ (read/write). To incorporate phase noise from the lasers used to generate these fields, the envelope fields at the ends of the waveguide are specified by stochastic boundary conditions, namely $a_1(0,t)=a_{\textrm {data}}(t) e^{i\phi _1(t)}$ for the data pulse and $a_2(L,t)=\left [a_{\textrm {read}}(t) + a_{\textrm {write}}(t)\right ]e^{i\phi _2(t)}$. The random phase terms $\phi _{1,2}$ are statistically uncorrelated Brownian motions described by
where $\textrm {d}{W(t)}$ is a standard Wiener increment [22], and $\Delta \nu _L$ is the laser’s intrinsic linewidth [19].In non-storage SBS setups [18,19] the pump has higher input power $P_{\textrm {p0}}$ than the Stokes $P_{\textrm {s0}}$ and thus behaves as an amplifier. However, in the case of Brillouin storage, we require $P_{\textrm {p0}} \ll P_{\textrm {s0}}$ so that the read/write pulse can completely deplete the pump, which contains the data to be stored. Optimum storage requires two conditions on the read/write pulse. First, the read/write pulses must be at least as short in duration as the shortest data pulse. Second, the pulse area for the read/write pulse, defined as [6,23]
The storage efficiency tells us how effective the storage system is, and may be defined as the ratio of the total output data power $|a_{\textrm {out}}(t)|^{2}$ to the total input data power $|a_{\textrm {data}}(t)|^{2}$ [10]:
2.2 Solution to the coupled mode equations
We numerically solve Eq. (1)$-$(3) following the numerical method described in [19]. In brief, we note that the propagation distance of the acoustic wave over the time-scale of the SBS interaction is very small [13]. We therefore apply the limit $\partial _z b \rightarrow 0$ in Eq. (3), and solve the equation in time via the Green function method [18]. The thermal noise function is simulated starting at thermal equilibrium at $t=0$, by drawing random numbers from a steady-state probability distribution corresponding to an Ornstein-Uhlenbeck process [19,26]. We then substitute the formal solution $b(z,t)$ into Eq. (1) and (2), and integrate the two equations using a two-step iteractive method: first, the optical fields are translated in space by a distance $\Delta z = v\Delta t$ while the acoustic field remains stationary [19]. Then, the optical fields at each point in space $z$ are updated by integrating the equations in time with an Euler-Mayurama scheme [27]. This two-step process is iterated over multiple time steps to simulate the dynamics of the input pulses.
For our simulations we use the parameters summarized in Table 1. We assume a high gain chalcogenide waveguide, of the type used in previous SBS experiments [28]. We store individual 8-bit packets one at a time, consisting of all 256 possible unique 8-bit sequences. Each packet corresponds to 8 individual pulses in the amplitude storage case, and 4 phases in the phase storage case, as shown in Fig. 2.
It should be noted that the model used here (Eq. (1)$-$(3)) makes the assumption of the slowly-varying envelope approximation (SVEA) and rotating-wave approximation (RWA) [18,21]. Because the acoustic frequencies used in SBS experiments are typically in the range of a few GHz [1,29], the optical pulses simulated must be longer than a few hundred picoseconds in duration so that these approximations remain valid [30]. Therefore, we limit our simulations to pulses no shorter than 100 ps, for both for the data and read/write cases. We have chosen to focus on the storage of 8-bit sequences because they can fit inside a 30 cm chalcogenide waveguide without being too short to violate the SVEA and RWA for a 7.8 GHz Brillouin frequency shift, since the acoustic pulses generated from 300 ps optical bits are approximately 3.7 cm in length.
2.3 Data encoding
The two storage schemes used in this paper — amplitude and phase storage — are summarized in Fig. 2. Each bit is defined as having a duration $\tau _{\textrm {bit}}$. For amplitude storage, we encode ones into Gaussian pulses of the same peak power $P_{\textrm {p0}}$, while zeros are represented by gaps of duration $\tau _{\textrm {bit}}$. For phase storage, we use the same scheme as in quadrature phase-shift keying (QPSK) with gray coding [31], where bit pairs are assigned a unique phase, namely $11 = 45^{\circ }$, $01 = 135^{\circ }$, $00 = 225^{\circ }$ and $10 = 315^{\circ }$. For a given input information packet, we quantify the retrieval accuracy in both storage schemes via the packet error rate (PER). This is similar to the bit error rate (BER) of a binary stream, except that we count correct 8-bit packets as opposed to counting individual correct bits. Therefore, the PER is the ratio of correctly retrieved packets with respect to the input data, and has a value $0 \leq \textrm {PER}\leq 1$.
In the amplitude storage case, we encode bits=1 into Gaussian pulses of full-width at half-maximum (FWHM) $\tau _1$, while bits=0 are represented by gaps of duration $\tau _{\textrm {bit}}$ in the data sequence. By default, we choose all pulses to have a phase of 0$^{\circ }$. In the retrieval stage, the output data power $|a_{\textrm {out}}(t)|^{2}$ is separated into equal intervals of length $\tau _{\textrm {bit}}$. We use a dynamic threshold technique: initially, a threshold power $P_{\textrm {thresh}}$ is set. Then, we record the output power at the center of the bit period, and a single bit is read as 0 if $|a_{\textrm {out}}^{(n)}(t)|^{2}<P_{\textrm {thresh}}$ and as 1 if $|a_{\textrm {out}}^{(n)}(t)|^{2}\geq P_{\textrm {thresh}}$ [32], for all 256 data packets. This process is repeated for different values of $P_{\textrm {thresh}}$ until the total PER has been minimized.
In the phase storage case, we encode bit-pairs into 4 different phases, as shown in Fig. 3(a), as part of an analytic smooth rectangular pulse (ASR) [33] with FWHM $\tau _{\textrm {data}}$. We then apply a Gaussian filter to remove non-physical discontinuities; to emulate the effect of the modulators used in SBS experiments we have chosen a filter width of 5 GHz, resulting in a smooth phase transition between the bit-pairs (see Fig. 2(f)). In the retrieval stage, $|a_{\textrm {out}}(t)|^{2}$ is separated into equal intervals of length $2\tau _{\textrm {bit}}$. We extract the phase from the amplified output signal as $\phi _{\textrm {out}}(t) = \tan ^{-1}\left (\textrm {Im}\left [a_{\textrm {out}}(t)\right ]/\textrm {Re}\left [a_{\textrm {out}}(t)\right ]\right )$. We assume ideal detector conditions and record the phase at the center of the $2\tau _{\textrm {bit}}$ period, such that if $0^{\circ }<\phi _{\textrm {out}}(t) <90^{\circ }$, we read the output bit-pair as 11, if $90^{\circ }<\phi _{\textrm {out}}(t) <180^{\circ }$ we read the bit-pair as 01 and so on. Additionally, we use another measure of signal integrity based on the constellation diagram data: let $z_{\textrm {out}}^{(m)}$ represent a single point on the constellation diagram corresponding to the $m$th bit-pair in the binary sequence (such as 00), with magnitude $\left |z_{\textrm {out}}^{(m)}\right |$ and phase $\phi _{z}^{\textrm {m}}$. The variances in each can be found via
3. Results and discussion
3.1 Effect of thermal noise
We investigate the effect of varying the read/write pulse peak power between 0.5$-$10.5 W, while maintaining a read/write pulse width of 100 ps. We include thermal noise into the waveguide at 300 K, but neglect the input laser phase noise by setting the laser linewidth $\Delta \nu _L$ to zero. In the amplitude storage case, we use Gaussian data pulses with 10 mW peak power and pulse width 150 ps, while in the phase storage case we use a rectangular pulse of duration 2.44 ns, and phase intervals of duration 600 ps. Each bit of optical information is 300 ps in length for both amplitude and phase storage. The results of the simulations are shown in Fig. 4. First, we observe in Fig. 4(a) and (b) that the storage efficiency in both storage schemes is higher as we increase the peak write pulse power, as the pulse area is lower than the optimum value given by Eq. (5). Second, we see in Fig. 4(c) and (d) that lower peak write powers lead to increased PER, thus reducing the maximum storage time achievable in both encoding schemes. This occurs because at lower peak read/write powers — which also corresponds to the lower storage efficiency regime — the coherent output data field is less distinguishable from the amplified spontaneous noise arising from the interaction between the read/write pulse and the thermal background fluctuations in the waveguide [18,19]. Consequently, this increases the probability of retrieving the wrong data sequence at the waveguide output ($z=L$).
Similarly, we see an increase in PER in both AS and PS schemes with increasing $\tau _{\textrm {sto}}$, as shown in Fig. 4(c) and (d). This occurs because in the time between the write and read-process, the acoustic wave containing the stored information decays at a rate $1/\tau _a$. As $\tau _{\textrm {sto}}$ increases, the acoustic wave gets closer to the background thermal noise, increasing the fluctuations in the retrieved data field during the read-process. This effect is more clearly illustrated in Fig. 5, where the constellation points spread out in both radial and angular directions, indicating an increase in both amplitude and phase noise in the retrieved data field. In Fig. 6 we see that the rate of increase in the variance of the phase and amplitude of each encoding scheme is the same for the first 40 nanoseconds. However, in Fig. 4(f) we see that the PER in the phase storage case begins to increase at longer $\tau _{\textrm {sto}}$ compared to the amplitude storage case, suggesting that phase-encoded data is more robust to thermal noise and hence allows for longer storage times. This occurs because the phase variations are primarily constrained to a single quadrant on the constellation plots (as shown in Fig. 5), whereas the amplitude variations reach the noise floor more rapidly. Consequently, the probability of detecting the wrong bit of information in the amplitude is higher compared to the phase. This is the same observation that would be expected from the additive white Gaussian noise (AWGN) model in discrete communication theory, where the probability of bit error for phase-shift-keying (PSK) is lower than for amplitude-shift-keying (ASK) [34,35]. In addition, phase encoding can transfer more bits per symbol, which allows for the use of longer optical pulses without loss of data capacity.
3.2 Effect of laser phase noise
We now investigate the effect of adding input laser phase noise, at a linewidth of 100 kHz, with the same pulse and thermal noise parameters as before. We generate the phase noise in the data pulses and read/write pulses independently so they are statistically uncorrelated, but have the same mean and variance properties. Figure 7 shows the results for storage efficiency and PER for both encoding schemes. The plots in Fig. 7(a)$-$(d) look very similar to Fig. 4(a)$-$(d), indicating that laser phase noise does not have a significant impact on the storage efficiency or PER. This is consistent with previous work on SBS noise in short pulses [19], where it was found that laser phase noise only has a significant impact on the SBS process when the laser coherence time ($\tau _{\textrm {coh}} = 1/\pi \Delta \nu _L$) is comparable in magnitude to the SBS interaction time. Next, the constellation diagram results in Fig. 8 reveal that the input laser phase noise broadens the variance in the phase of the retrieved data field, making the distribution of constellation points slightly more elliptical compared to the thermal noise only case. This is also shown in Fig. 9, where the rates of increase for the amplitude and phase variances still remain approximately the same for the first 40 nanoseconds of storage time. This occurs because at this regime of powers for the data pulses, the SBS process acts as a linear amplifier, hence the total noise in the retrieved data is a linear combination of the waveguide thermal noise and the laser phase noise from the inputs.
4. Conclusion
We have numerically simulated the Brillouin storage of different data packets with thermal and laser noise, using amplitude storage and phase storage techniques in a photonic waveguide. Through these computer simulations, we have shown that phase encoded storage allows for longer storage times than amplitude encoded storage. This is because phase encoding is more robust to noise than amplitude encoding, in accordance with the additive-white-Gaussian-noise model of discrete communications theory [34]. It is therefore possible to increase Brillouin storage time by encoding information into the phase of the data field, without having to change the waveguide or laser properties. Furthermore, because phase storage techniques can encode more bits per symbol than amplitude storage [31,34], one can use longer optical pulses without loss of data capacity.
Funding
UTS Australian Government Research Training Program Scholarship (00099F); Macquarie University Research Fellowship Scheme (MQRF0001036); Australian Research Council (DP160101691, DP200101893).
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.
Data availability
Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.
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. A. Kobyakov, M. Sauer, and D. Chowdhury, “Stimulated Brillouin scattering in optical fibers,” Adv. Opt. Photonics 2(1), 1–59 (2010). [CrossRef]
3. 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]
4. R. W. Boyd, Nonlinear optics (Elsevier, 2003).
5. 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]
6. Z. Zhu, D. J. Gauthier, and R. W. Boyd, “Stored light in an optical fiber via stimulated Brillouin scattering,” Science 318(5857), 1748–1750 (2007). [CrossRef]
7. V. Kalosha, W. Li, F. Wang, L. Chen, and X. Bao, “Frequency-shifted light storage via stimulated Brillouin scattering in optical fibers,” Opt. Lett. 33(23), 2848–2850 (2008). [CrossRef]
8. M. Merklein, B. Stiller, and B. J. Eggleton, “Brillouin-based light storage and delay techniques,” J. Opt. 20(8), 083003 (2018). [CrossRef]
9. C.-H. Dong, Z. Shen, C.-L. Zou, Y.-L. Zhang, W. Fu, and G.-C. Guo, “Brillouin-scattering-induced transparency and non-reciprocal light storage,” Nat. Commun. 6(1), 6193–6196 (2015). [CrossRef]
10. 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]
11. M. Merklein, B. Stiller, K. Vu, P. Ma, S. J. Madden, and B. J. Eggleton, “On-chip broadband nonreciprocal light storage,” Nanophotonics 10(1), 75–82 (2020). [CrossRef]
12. B. Stiller, M. Merklein, C. Wolff, K. Vu, P. Ma, S. J. Madden, and B. J. Eggleton, “Coherently refreshing hypersonic phonons for light storage,” Optica 7(5), 492–497 (2020). [CrossRef]
13. R. W. Boyd, K. Rzazewski, and P. Narum, “Noise initiation of stimulated Brillouin scattering,” Phys. Rev. A 42(9), 5514–5521 (1990). [CrossRef]
14. 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]
15. 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]
16. P. Kharel, R. Behunin, W. Renninger, and P. Rakich, “Noise and dynamics in forward Brillouin interactions,” Phys. Rev. A 93(6), 063806 (2016). [CrossRef]
17. 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]
18. O. A. Nieves, M. D. Arnold, M. Steel, M. K. Schmidt, and C. G. Poulton, “Noise and pulse dynamics in backward stimulated Brillouin scattering,” Opt. Express 29(3), 3132–3146 (2021). [CrossRef]
19. O. A. Nieves, M. D. Arnold, M. J. Steel, M. K. Schmidt, and C. G. Poulton, “Numerical simulation of noise in pulsed Brillouin scattering,” J. Opt. Soc. Am. B 38(8), 2343–2352 (2021). [CrossRef]
20. 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]
21. 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]
22. W. Horsthemke, “Noise induced transitions,” in Non-equilibrium dynamics in chemical systems, (Springer, 1984), pp. 44–49.
23. M. Dong and H. G. Winful, “Area dependence of chirped-pulse stimulated Brillouin scattering: implications for stored light and dynamic gratings,” J. Opt. Soc. Am. B 32(12), 2514–2519 (2015). [CrossRef]
24. H. G. Winful, “Chirped Brillouin dynamic gratings for storing and compressing light,” Opt. Express 21(8), 10039–10047 (2013). [CrossRef]
25. L. Allen and J. H. Eberly, Optical resonance and two-level atoms, vol. 28 (John Wiley & Sons., 1975).
26. G. E. Uhlenbeck and L. S. Ornstein, “On the theory of the Brownian motion,” Phys. Rev. 36(5), 823–841 (1930). [CrossRef]
27. P. Kloeden and E. Platen, Numerical Solution of Stochastic Differential Equations (Springer-Verlag, Berlin Heidelberg, 1992), 1st ed.
28. 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]
29. C. Wolff, M. Smith, B. Stiller, and C. Poulton, “Brillouin scattering—theory and experiment: tutorial,” J. Opt. Soc. Am. B 38(4), 1243–1269 (2021). [CrossRef]
30. J. Piotrowski, M. K. Schmidt, B. Stiller, C. G. Poulton, and M. J. Steel, “Picosecond acoustic dynamics in stimulated Brillouin scattering,” Opt. Lett. 46(12), 2972–2975 (2021). [CrossRef]
31. S. Faruque, Radio frequency modulation made easy (Springer, 2017), pp. 69–83.
32. D. Walsh, D. Moodie, I. Mauchline, S. Conner, W. Johnstone, and B. Culshaw, “Practical bit error rate measurements on fibre optic communications links in student teaching laboratories,” in 9th International Conference on Education and Training in Optics and Photonics (ETOP), Marseille, France, Paper ETOP021, (2005).
33. E. Granot, “Propagation of chirped rectangular pulses in dispersive media: analytical analysis,” Opt. Lett. 44(19), 4745–4748 (2019). [CrossRef]
34. D. R. Smith, Digital transmission systems (Springer science & business media, 2012), pp. 370–381.
35. D. Bala, “Analysis the probability of bit error performance on different digital modulation techniques over AWGN channel using MATLAB,” J. Electr. Eng. Electron. Control. Comput. Sci. 7, 9–18 (2021).