Abstract
Continuous variable-based quantum cryptography (CV-QKD) is an emerging field in quantum information science, offering unprecedented security for communication protocols by harnessing the principles of quantum mechanics. However, ocean environments pose unique challenges to quantum communication due to their distinct properties and characteristics. This work investigates the impact of turbulence on the transmission of Gaussian light beams used in a continuous variable-based quantum key distribution system for underwater quantum communication. The objective is to quantitatively analyze the induced losses and propose methodologies to mitigate their effects. To achieve this, we adopt the widely accepted ABCD matrix formalism, which provides a comprehensive framework for characterizing the propagation of optical beams through different media. Moreover, a numerical simulation framework is developed to assess the resulting losses and evaluate the performance of the proposed system. The implications of these numerical simulation frameworks for the design and optimization of quantum communication systems for oceanic environments are thoroughly discussed.









Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.Data Availability
Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.
References
Zeng, Z., Fu, S., Zhang, H., Dong, Y., Cheng, J.: A survey of underwater optical wireless communications. IEEE Commun. Surv. Tutor. 19(1), 204–238 (2016)
Hoeher, P.A., Sticklus, J., Harlakin, A.: Underwater optical wireless communications in swarm robotics: a tutorial. IEEE Commun. Surv. Tutor. 23(4), 2630–2659 (2021)
Cossu, G., Sturniolo, A., Messa, A., Scaradozzi, D., Ciaramella, E.: Full-fledged 10base-t ethernet underwater optical wireless communication system. IEEE J. Sel. Areas Commun. 36(1), 194–202 (2017)
Bennett, C., Brassard, G. i.: “Proc. ieee int. conf. on computers, systems and signal processing, bangalore, india,” in Proc. of the IEEE Int. Conf. on Computers Systems and Signal Processing Bangalore India, 175, (1984)
Shor, P.W., Preskill, J.: Simple proof of security of the bb84 quantum key distribution protocol. Phys. Rev. Lett. 85(2), 441 (2000)
Mayers, D.: Unconditional security in quantum cryptography. J. ACM (JACM) 48(3), 351–406 (2001)
Loepp, S., Wootters, W.K.: Protecting information: from classical error correction to quantum cryptography. Cambridge University Press, Cambridge (2006)
Grosshans, F., Grangier, P.: Continuous variable quantum cryptography using coherent states. Phys. Rev. Lett. 88(5), 057902 (2002)
Srikara, S., Thapliyal, K., Pathak, A.: Continuous variable b92 quantum key distribution protocol using single photon added and subtracted coherent states. Quantum Inf. Process. 19, 1–16 (2020)
Pirandola, S., Andersen, U.L., Banchi, L., Berta, M., Bunandar, D., Colbeck, R., Englund, D., Gehring, T., Lupo, C., Ottaviani, C., et al.: Advances in quantum cryptography. Adv. Opt. Photonics 12(4), 1012–1236 (2020)
Ralph, T.C.: Continuous variable quantum cryptography. Phys. Rev. A 61(1), 010303 (1999)
Diamanti, E., Leverrier, A.: Distributing secret keys with quantum continuous variables: principle, security and implementations. Entropy 17(9), 6072–6092 (2015)
Laudenbach, F., Pacher, C., Fung, C.-H.F., Poppe, A., Peev, M., Schrenk, B., Hentschel, M., Walther, P., Hübel, H.: Continuous-variable quantum key distribution with gaussian modulation-the theory of practical implementations. Adv. Quantum Technol. 1(1), 1800011 (2018)
Ghalaii, M., Pirandola, S.: Quantum communications in a moderate-to-strong turbulent space. Commun. Phys. 5(1), 38 (2022)
Heim, B., Peuntinger, C., Killoran, N., Khan, I., Wittmann, C., Marquardt, C., Leuchs, G.: Atmospheric continuous-variable quantum communication. New J. Phys. 16(11), 113018 (2014)
Sharma, V., Banerjee, S.: Analysis of atmospheric effects on satellite-based quantum communication: a comparative study. Quantum Inf. Process. 18, 1–24 (2019)
Sharma, V., Banerjee, S.: “Analysis of quantum key distribution based satellite communication,” in 2018 9th International Conference on Computing, Communication and Networking Technologies (ICCCNT), 1–5, IEEE, (2018)
Meena, Muskan, R., Banerjee, S.: “Analysing qber and secure key rate under various losses for satellite based free space qkd,” (2023)
Huang, P., He, G., Fang, J., Zeng, G.: Performance improvement of continuous-variable quantum key distribution via photon subtraction. Phys. Rev. A 87(1), 012317 (2013)
Ye, W., Zhong, H., Wu, X., Hu, L., Guo, Y.: Continuous-variable measurement-device-independent quantum key distribution via quantum catalysis. Quantum Inf. Process. 19, 1–22 (2020)
Li, Z., Zhang, Y., Wang, X., Xu, B., Peng, X., Guo, H.: Non-gaussian postselection and virtual photon subtraction in continuous-variable quantum key distribution. Phys. Rev. A 93(1), 012310 (2016)
Giuliano, G.: Underwater optical communication systems. PhD thesis, University of Glasgow, (2019)
Gussen, C.M., Diniz, P.S., Campos, M.L., Martins, W.A., Costa, F.M., Gois, J.N.: A survey of underwater wireless communication technologies. J. Commun. Inf. Syst. 31(1), 242–255 (2016)
Tian, P., Liu, X., Yi, S., Huang, Y., Zhang, S., Zhou, X., Hu, L., Zheng, L., Liu, R.: High-speed underwater optical wireless communication using a blue gan-based micro-led. Opt. Express 25(2), 1193–1201 (2017)
Wang, J., Lu, C., Li, S., Xu, Z.: 100 m/500 mbps underwater optical wireless communication using an nrz-ook modulated 520 nm laser diode. Opt. Express 27(9), 12171–12181 (2019)
Zhang, L., Wang, Z., Wei, Z., Chen, C., Wei, G., Fu, H., Dong, Y.: Towards a 20 gbps multi-user bubble turbulent noma uowc system with green and blue polarization multiplexing. Opt. Express 28(21), 31796–31807 (2020)
Dai, Y., Chen, X., Yang, X., Tong, Z., Du, Z., Lyu, W., Zhang, C., Zhang, H., Zou, H., Cheng, Y., et al.: 200-m/500-mbps underwater wireless optical communication system utilizing a sparse nonlinear equalizer with a variable step size generalized orthogonal matching pursuit. Opt. Express 29(20), 32228–32243 (2021)
Fei, C., Wang, Y., Du, J., Chen, R., Lv, N., Zhang, G., Tian, J., Hong, X., He, S.: 100-m/3-gbps underwater wireless optical transmission using a wideband photomultiplier tube (pmt). Opt. Express 30(2), 2326–2337 (2022)
Duntley, S.Q.: Light in the sea. JOSA 53(2), 214–233 (1963)
Gilbert, G., Stoner, T., Jernigan, J.: “Underwater experiments on the polarization, coherence, and scattering properties of a pulsed blue-green laser,” in Underwater Photo Optics I, 7, 8–14, SPIE, (1966)
Borah, D.K., Voelz, D.G.: Estimation of laser beam pointing parameters in the presence of atmospheric turbulence. Appl. Opt. 46(23), 6010–6018 (2007)
Muskan, A. Dutta., Banerjee, S., Pathak, A.: “Analysis for satellite-based high-dimensional extended B92 and high-dimensional BB84 quantum key distribution,” arXiv preprint[SPACE]arXiv:2311.00309, (2023)
Coronel, J., Elayoubi, K., Al Ahmadi, A., Al Hosani, S., Al Blooshi, A., Al Ameri, R., Cordette, S., Bouchalkha, A., Alameri, J., Matras, G. et al., “Characterization of infrared laser beam through atmospheric optical turbulence in laboratory environment,” in Free-Space Laser Communications XXXV, 12413, pp. 455–467, SPIE, (2023)
Wu, H., van Iersel, M.: “Experimental study on the effects of oam beams propagating through atmospheric turbulence,” in Laser Communication and Propagation through the Atmosphere and Oceans XII, 12691, pp. 132–142, SPIE, (2023)
Tang, S., Zhang, X., Dong, Y.: “Temporal statistics of irradiance in moving turbulent ocean,” in 2013 MTS/IEEE OCEANS-Bergen, 1–4, IEEE, (2013)
Andrews, L.C., Phillips, R.L.: Laser beam Propag. Rand. Med. Second Edition, Laser Beam Propagation Through Random Media (2005)
Zhang, J., Kou, L., Yang, Y., He, F., Duan, Z.: Monte-carlo-based optical wireless underwater channel modeling with oceanic turbulence. Opt. Commun. 475, 126214 (2020)
Enghiyad, N., Sabbagh, A.G.: Impulse response of underwater optical wireless channel in the presence of turbulence, absorption, and scattering employing monte carlo simulation. JOSA A 39(1), 115–126 (2022)
Geldard, C. T., Thompson, J., Popoola, W. O.: “Effects of turbulence induced scattering on underwater optical wireless communications,” arXiv preprint[SPACE]arXiv:2008.01152, (2020)
Xu, D., Yue, P., Yi, X., Liu, J.: Improvement of a monte-carlo-simulation-based turbulence-induced attenuation model for an underwater wireless optical communications channel. JOSA A 39(8), 1330–1342 (2022)
Ji, X., Yin, H., Liang, Y., Wang, J., et al.: Analysis of aperture averaging effect and communication system performance of wireless optical channels with from weak to strong turbulence in natural turbid water. Opt. Commun. 528, 129018 (2023)
Hou, W., Woods, S., Jarosz, E., Goode, W., Weidemann, A.: Optical turbulence on underwater image degradation in natural environments. Appl. Opt. 51, 2678–2686 (2012)
Nootz, G., Jarosz, E., Dalgleish, F.R., Hou, W.: Quantification of optical turbulence in the ocean and its effects on beam propagation. Appl. Opt. 55, 8813–8820 (2016)
Quan, X., Fry, E.S.: Empirical equation for the index of refraction of seawater. Appl. Opt. 34(18), 3477–3480 (1995)
Kogelnik, H.: Imaging of optical modes-resonators with internal lenses. Bell Syst. Tech. J. 44(3), 455–494 (1965)
Kogelnik, H.: On the propagation of gaussian beams of light through lenslike media including those with a loss or gain variation. Appl. Opt. 4, 1562–1569 (1965)
Tu, J., Yeoh, G. H., Liu, C., Tao, Y.: Computational fluid dynamics: a practical approach. Elsevier, (2023)
Jaruwatanadilok, S.: Underwater wireless optical communication channel modeling and performance evaluation using vector radiative transfer theory. IEEE J. Sel. Areas Commun. 26(9), 1620–1627 (2008)
Malpani, P., Alam, N., Thapliyal, K., Pathak, A., Narayanan, V., Banerjee, S.: Lower-and higher-order nonclassical properties of photon added and subtracted displaced fock states. Ann. Phys. 531(2), 1800318 (2019)
Meena, R., Banerjee, S.: Characterization of quantumness of non-gaussian states under the influence of gaussian channel. Quantum Inf. Process. 22(8), 298 (2023)
Malpani, P., Thapliyal, K., Alam, N., Pathak, A., Narayanan, V., Banerjee, S.: Impact of photon addition and subtraction on nonclassical and phase properties of a displaced fock state. Opt. Commun. 459, 124964 (2020)
Zhong, H., Guo, Y., Mao, Y., Ye, W., Huang, D.: Virtual zero-photon catalysis for improving continuous-variable quantum key distribution via gaussian post-selection. Sci. Rep. 10(1), 17526 (2020)
Xiang, Y., Wang, Y., Ruan, X., Zuo, Z., Guo, Y.: Improving the discretely modulated underwater continuous-variable quantum key distribution with heralded hybrid linear amplifier. Phys. Scr. 96(6), 065103 (2021)
McNeil, G.T.: Metrical fundamentals of underwater lens system. Opt. Eng. 16(2), 128–139 (1977)
Griffiths, D. J.: “Introduction to electrodynamics,” (2005)
Leverrier, A., Grangier, P.: Continuous-variable quantum-key-distribution protocols with a non-gaussian modulation. Phys. Rev. A 83(4), 042312 (2011)
Peng, Q., Chen, G., Li, X., Liao, Q., Guo, Y.: Performance improvement of underwater continuous-variable quantum key distribution via photon subtraction. Entropy 21(10), 1011 (2019)
Zhao, Y., Zhang, Y., Xu, B., Yu, S., Guo, H.: Continuous-variable measurement-device-independent quantum key distribution with virtual photon subtraction. Phys. Rev. A 97(4), 042328 (2018)
Navascués, M., Grosshans, F., Acin, A.: Optimality of gaussian attacks in continuous-variable quantum cryptography. Phys. Rev. Lett. 97(19), 190502 (2006)
García-Patrón, R., Cerf, N.J.: Unconditional optimality of gaussian attacks against continuous-variable quantum key distribution. Phys. Rev. Lett. 97(19), 190503 (2006)
Wolf, M.M., Giedke, G., Cirac, J.I.: Extremality of gaussian quantum states. Phys. Rev. Lett. 96(8), 080502 (2006)
Farwell, N. H., Korotkova, O.: “Multiple phase-screen simulation of oceanic beam propagation,” in Laser communication and propagation through the atmosphere and oceans III, 9224, pp. 374–380, SPIE, (2014)
Johnson, L.J., Green, R.J., Leeson, M.S.: Underwater optical wireless communications: depth dependent variations in attenuation. Appl. Opt. 52, 7867–7873 (2013)
Raymont, J. E.: Plankton & productivity in the oceans: Volume 1: Phytoplankton. Elsevier, (2014)
Uitz, J., Claustre, H., Morel, A., Hooker, S. B.: “Vertical distribution of phytoplankton communities in open ocean: An assessment based on surface chlorophyll,” Journal of Geophysical Research: Oceans, 111(C8), (2006)
Hooker, Stanford, B.: “Seawifs technical report series: an overview of seawifs and ocean color,” NASA Technical Memorandum, I, (1992)
Kameda, T., Matsumura, S.: Chlorophyll biomass off sanriku, northwestern pacific, estimated by ocean color and temperature scanner (octs) and a vertical distribution model. J. Oceanogr. 54, 509–516 (1998)
Pontbriand, C., Farr, N., Ware, J., Preisig, J., Popenoe, H.: “Diffuse high-bandwidth optical communications,” in OCEANS 2008, 1–4, IEEE, (2008)
Bricaud, A., Babin, M., Morel, A., Claustre, H.: Variability in the chlorophyll-specific absorption coefficients of natural phytoplankton: analysis and parameterization. J. Geophys. Res. Oceans 100(C7), 13321–13332 (1995)
Acknowledgements
The author would like to thank CSIR for the fellowship support. SB acknowledges support from Interdisciplinary Cyber Physical Systems (ICPS) program of the Department of Science and Technology (DST), India, Grant No.:DST/ICPS/QuST/Theme-1/2019/6. SB also acknowledges support from the Interdisciplinary Research Platform (IDRP) on Quantum Information and Computation (QIC) at IIT Jodhpur.
Author information
Authors and Affiliations
Contributions
1. Ramniwas Meena - Conceptualized and designed the research model. - Contributed to the analyse design and methodology. - Drafted the initial manuscript and revised it critically for intellectual content. - Approved the final version of the manuscript. 2. Subhashish Banerjee - Provided guidance throughout the research process. - Assisted in data interpretation, critically revised the manuscript, and provided valuable intellectual input. - Approved the final version of the manuscript. Both authors have read and approved the final manuscript for submission.
Corresponding author
Ethics declarations
Conflict of interest
The authors declare no conflict of interest related to this research.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendices
Appendix A
1.1 Absorption
Different substances in ocean have specific absorption spectra, meaning they absorb light at different wavelengths. For example, chlorophyll, a pigment found in marine plants and phytoplankton, absorbs light primarily in the blue and red parts of the spectrum while reflecting green light, giving the ocean its characteristic blue color. The absorption properties of ocean have significant implications for marine life, as they determine the available light for photosynthesis and affect the distribution of marine organisms in the water column. They are also essential for understanding the optical properties of the ocean, which are critical for applications such as remote sensing, oceanography, and underwater communication. The full absorption spectrum multiplied by their respective concentrations, such that [63]
where \(a_w\) is the pure water absorption coefficient in \(m^{-1}\), \(a_f^{0}\) is the specific absorption coefficient of fulvic acid in \(m^{-1}\), \(a_h^{0}\) is the specific absorption coefficient of humic acid in \(m^{-1}\), \(a_c^{0}\) is the specific absorption coefficient of chlorophyll in \(m^{-1}\), \(C_f\) is the concentration of fulvic acid in \(mg/m^3\), \(C_h\) is the concentration of humic acid in \(mg/m^3\), \(C_c\) is the concentration of chlorophyll-a in \(mg/m^3\), \(k_f\) is the fulvic acid exponential coefficient (\(k_f = 0.0189 m^{-1}\)), and \(k_h\) is the humic acid exponential coefficient (\(k_h = 0.01105 m^{-1}\)).
1.1.1 Scattering
The scattering spectra in ocean are impacted by two main biological factors: scattering caused by pure water and scattering resulting from particulate substances. The latter can be further divided into two categories: small and large particles, each exhibiting distinct statistical distributions and scattering properties. These factors collectively contribute to the overall scattering behavior in the ocean, influencing how light interacts with the water and its constituents. The full form of the equation is given as [63]:
where the pure water scattering coefficient is represented as \(b_w\) in \(m^{-1}\), \(b_{s}^{0}\) is the scattering coefficient for small particulate matter in \({m^{2}/{g}}\), the scattering coefficient for large particulate matter is denoted as \({b_{l}^{0}}\) and is also measured in \({m^2}/g\), \({C_s}\) represents the concentration of small particles in \(g/{m^{3}}\), and \({C_l}\) represents the concentration of large particles in \(g/{m^{3}}\).
Chlorophyll-a is a predominant substance found in phytoplankton, a group of microscopic organisms. These photosynthesizing organisms thrive in the photic or euphotic zone of the ocean, where sunlight can penetrate. For efficient photosynthesis, phytoplankton requires an adequate supply of nutrients, which are usually more abundant in coastal areas due to land runoff and upwelling of subsurface water into the photic zone [64]. The spectral dependencies for the scattering coefficients of small and large particulate matter are given by the following equations:
Scattering contributes much less to the overall attenuation coefficient than absorption, through it is much greater in particulate-rich areas.
Monitoring the near-surface chlorophyll concentration in the ocean is essential, and the NASA SeaWiFS (Sea-viewing Wide Field-of-View Sensor) project employs ocean color observation and quantification to determine chlorophyll concentration [66]. Higher chlorophyll concentrations are typically observed along the equator, east-facing coastlines, and in high latitude regions. In the open ocean, typical chlorophyll concentrations range from 0.01 to 4.0 \({\mathrm{mg/m}}^{3}\), while near-shore levels can be as high as 60 \({\mathrm{mg/m}}^{3}\). The chlorophyll profile over a depth d(m) from the surface \(C_c(d)\) can be modeled as a Gaussian curve that includes five numerically determined parameters and has the generic form [67]:
where \(B_0\) is the background chlorophyll concentration on the surface, \(\textit{S}\) is the vertical gradient of concentration, which is always negative due to the slow decrease in chlorophyll concentration with depth, h is the total chlorophyll above the background level and standard deviation of chlorophyll concentration (\(\sigma \)):
Haltrin proposed a simplified one-parameter model for attenuation, which establishes the relationship between the concentrations of different particulates [68]. These particulate concentrations were determined numerically in relation to the chlorophyll concentration and are expressed as follows:
Recall that the specific attenuation of chlorophyll also contained a chlorophyll concentration term; so now too has depth dependency:
where \(A(\lambda )\) and \(B(\lambda )\) represent empirical constants, a full list of which is given in [69].
Appendix B: covariance matrix of m-Photon-subtracted TMSV state
Suppose \(\gamma ^{(m)}_{AB_2}\) represents the covariance matrix of \(\rho ^{(m)}_{AB_2}\), and it has the following formula:
where \(V_A\) and \(V_B\) represent the variances in modes \(A\) and \(B_2\), respectively. The covariance between the quadratures of modes \(A\) and \(B_2\) is denoted as \(\phi _{AB}\). The matrices \(\mathbb {I}\) and \(\sigma _z\) are defined as diagonal matrices with elements \((1,1)\) and \((1,-1)\), respectively.
Suppose \(x'\), \(p'\) are the heterodyne measurement results of mode A, and x is the homodyne measurement result of mode B. Then,
where
the weighting function W and probability \(P(x',p')\) are described in Sect. 4.
After simplifying Eqs. (B2), the resulting expressions are as follows:
where \(\bar{V}\) is defined as \(\bar{V}=\int {x'^{2} W\;P({x',p'}) \,dx' \,dp'}\) and further calculations yield:
Here, \(m\) represents the subtracted photon, \(\zeta \) is the squeezing parameter of the two-mode squeezed vacuum (TMSV) source, and \(T\) is the transmittance of the beam splitter (\(BS_1\)).
Appendix C
From Eq.(4.6), we have the covariance matrix for the final output state,
where \(\mathbb {I}\) is diag(1, 1), and \(\sigma _z\) is \(diag(1,-1)\), and Alice always uses heterodyne detection. The mutual information between Alice and Bob can be defined as follows:
where \(V_a=(V_1+1)/2\), \(V_b=V_2\), and
If we assume that Eve can purify the entire system, the Holevo quantity for homodyne measurements can be expressed as:
where S(AB) is a function of the symplectic eigenvalues \(\lambda _{1,2}\) of output matrix, which is
where
and
Here we have used the notations \(\Delta =V^{2}_{1}+V^{2}_{2}-2\phi ^{2}\), and \(D=V_1 V_2-\phi ^2\). Also, \(S(A|B)=G\left( \frac{{\lambda _3 - 1}}{2}\right) \) is a function of the symplectic eigenvalue \(\lambda _3\) of the covariance matrix \(\gamma ^{b}_A\) of the A mode after Bob’s homodyne detection, where \(\lambda _3 = \sqrt{V_1 \left( V_1 - \frac{\phi ^2}{V_2}\right) } \). The covariance matrix \(\gamma ^{b}_A\) can be described as follow [13]:
where \(\gamma _A=V_1 \mathbb {I}\), \(\Sigma _c=V_2 \sigma _z\), \(\Pi _{x}=diag(1,0)\) (in case x is measured) and \(\Pi _{p}=diag(0,1)\) (in case p is measured).
Appendix D: Channel transmission including excess noise
Suppose Alice and Bob share an input covariance matrix, represented as:
Here, initially \(V_A=V_B\). The channel loss can be modeled using a beam splitter characterized by its transmittance \( T_c \). To incorporate the noise introduced by the beam splitter, an additional noise mode is added to the covariance matrix. This noise mode is represented by a \( 2 \times 2 \) unit matrix and is integrated into the total state through direct summation, resulting in the following total covariance matrix [13]:
where \( n = N \mathbb {I} \).
The symplectic representation of the beam splitter matrix is given by:
with \( t = \sqrt{T_c}\; \mathbb {I} \) and \( r = \sqrt{1-T_c}\; \mathbb {I} \). To match the dimensions of the covariance matrix \( V_{out} \), a \( 2 \times 2 \) unitary matrix representing the beam splitter’s action on Alice’s mode is appended, yielding:
In this configuration, the beam splitter acts on Bob’s mode and the noise mode, while leaving Alice’s mode unchanged. The transformation of a Gaussian state, represented by a covariance matrix \( V \), under the action of a symplectic operator \( S \), follows the rule:
Substituting the values, this transformation expands to:
This yields the transformed covariance matrix:
From this result, the variance of the quadrature in Bob’s mode is:
Defining the noise input as:
the covariance matrix takes its final form:
This is Eq. (4.6).
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Meena, R., Banerjee, S. Continuous variable-based quantum communication in the ocean. Quantum Inf Process 24, 46 (2025). https://doi.org/10.1007/s11128-025-04664-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s11128-025-04664-2