Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2511.13622 · Physical Review A

Modeling Quantum Noise in Nanolasers using Markov Chains

Matias Bundgaard-Nielsen, Gian Luca Lippi, Jesper Mørk

quant-phphysics.opticsSubmitted Nov 17, 2025• Announced Nov 18, 2025• 0 citations
v1 · Nov 17, 2025v2 · Feb 2, 2026v3 · Mar 16, 2026
arXivHTMLPDF

Abstract

The random nature of spontaneous emission leads to unavoidable fluctuations in a laser's output. This is often included through random Langevin forces in laser rate equations, but this approach falls short for nanolasers. In this paper, we show that the laser quantum noise can be quantitatively computed for a very broad class of lasers by starting from simple and intuitive rate equations and merely assuming that the number of photons and excited electrons only takes discrete values. The success of the model is explained by showing that it constitutes a Markov chain, which can be derived from the full master equations. We show that in the many-photon limit, the model simplifies to Langevin equations. We perform an extensive comparison of different approaches for computing quantum noise in lasers, identifying the best approach for different system sizes, ranging from nanolasers to macroscopic lasers, and different levels of excitation, i.e., cavity photon number. In particular, we find that the numerical solution to the Langevin equations is inaccurate below the laser threshold, while the laser Markov chain model, on the other hand, is accurate for all pump values and laser sizes when collective emitter effects are excluded.

Journal reference
Phys. Rev. A 113, 043704, 2026
DOI
10.1103/f6f5-8dv8
Comments
17 Pages, 7 Figures, submitted for PRA

Related papers

Stochastic Approach to the Quantum Noise of a Single-Emitter Nanolaser

Matias Bundgaard-Nielsen, Emil V. Denning, Marco Saldutti +1 more

It is shown that the quantum intensity noise of a single-emitter nanolaser can be accurately computed by adopting a stochastic interpretation of the standard rate equation model. The only assumption made is that the emitter excitation and photon number are stochastic variables with integer values. This extends the validity of rate equations beyond the mean-field limit and avoids using the standard Langevin approach, which is shown to fail for few emitters. The model is validated by comparison to full quantum simulations of the relative intensity noise and second-order intensity correlation function, g^{(2)}(0). Surprisingly, even when the full quantum model displays vacuum Rabi oscillations, which are not accounted for by rate equations, the intensity quantum noise is correctly predicted by the stochastic approach. Adopting a simple discretization of the emitter and photon populations, thus, goes a long way in describing quantum noise in lasers. Besides providing a versatile and easy-to-use tool for modeling emerging nanolasers, these results provide insight into the fundamental nature of quantum noise in lasers.

arXiv:2301.11815•Physical Review Letters•Jun 22, 2023•11 citations

Multiplicative Noise in Nanolaser Amplification: From Discrete Emission Events to the Additive-Noise Limit

G. Sonnino, Ph. Peeters, P. Nardone +1 more

Laser amplification kicks off with just a few photons in the cavity, a regime where standard additive Langevin noise cannot capture the real fluctuations. To understand how this amplification begins, we map out a single-mode photon-carrier model that links macroscopic lasers directly into nanolasers. Stationary rate equations reveal that stimulated emission takes over across a tight pump window. In nanolasers, however, spontaneous processes coexist across a much broader range. Because steady-state averages miss the individual events that lead to the macroscopic photon number, we frame the stochastic dynamics around the actual discrete transitions, deriving a state-dependent diffusion approximation. The destruction of an excited carrier concurrent with the generation of a photon creates a negative photon-carrier cross-correlation. This unavoidably leads to a correlated, multiplicative Gaussian-noise framework. We track the resulting discrete dynamics with an established, event-based Laser Markov Chain (LMC) which stochastically samples pump, relaxation, emission, and cavity-loss events keeping all physical quantities as non-negative integers. Comparing its predictions to those of our multiplicative model proves that the latter reproduces excellently well the statistical behavior of photons and carrier, as long as a continuous diffusion description holds. Once the emission drops into the few-photon regime, however, the discrete model becomes indispensable. On the other hand, at large photon counts relative fluctuations shrink to a level where the state-dependence of the diffusion matrix can be ignored, to leading order, thus returning to the classic Langevin description of the system's stochastic dynamics.

quant-pharXiv:2610.25520•Oct 20, 2026

Quantum theory of a single-emitter nanolaser

Ларионов Николай Владимирович, Mikhail I. Kolobov

Using two equivalent methods, the Heisenberg-Langevin and the master equation, we investigate the properties of a single-emitter nanolaser: an incoherently pumped single two-level system interacting with a single cavity mode of finite finesse. We are interested in a stationary generation regime when the mean number of photons inside the laser cavity is much larger than unity. In this case we linearize the Heisenberg-Langevin equations around the stationary classical solution and obtain analytical results for the linewidth, the quadrature fluctuation spectra, and the intracavity Mandel $Q$ parameter. As the second method we use the master equation for the density operator of the laser field. Using the coherent-state representation of the density matrix we obtain the stationary solution for the Glauber $P$ function of the laser mode and investigate the quantum statistics of the laser field and its threshold behavior. We compare the results of both methods with numerical simulations and find very good agreement. In the end we use our theoretical results for analyzing experiment with single-ion nanolasers.

arXiv:1301.5652•Physical Review A•Jul 26, 2013•16 citations

Properties of nanolasers based on few discrete emitters

Anders Lund, P. Nielsen, Michael Lorke +2 more

The laser has evolved from table size apparatuses to truly nano sized devices, in much the same way that computer chips have been continuously minimized. The few-emitter nanolaser represents an extreme in terms of size. The emitters can be either atoms or quantum dots, and they are coupled to a high-Q optical cavity. This system has previously been studied for the one and two emitters case, see e.g. [1, 2]. The emitters and cavity are modelled as two-level systems and a harmonic oscillator, respectively. The coupled system is modelled using the Jaynes-Cummings Hamiltonian and the two-level systems are pumped incoherently by a rate P. Solutions are found using the corresponding master equation. However, with cavity populations exceeding 100 and several emitters, the dimension of the Hilbert space of the system becomes too large to handle efficiently on a conventional computer. E.g. for four emitters and 100 photon states the density matrix has more than 2.5 × 106 elements. We have been able to simplify the problem significantly by adiabatically eliminating the photon-assisted polarizations and the correlations between emitters and cavity [3]. This results in a set of rate equations for the population of the cavity, na, and the occupation of the emitters. Fig. 1a) shows na for up to 4 emitters coupled to a cavity, as a function of the pumping rate P. The figure also shows results from a full master equation solution, and the correspondence is very good for large values of P. In Fig. 1b) the second order correlation functions g(2)(0) are shown, as obtained using the full model. The second order correlation functions become 1 for sufficiently large values of P signifying the onset of lasing. However, for larger values of P the laser quenches and the emission becomes chaotic (g(2)(0) = 2). This quenching effect is well-known in two-level models [3], but might not be realistic for semiconductor quantum dots. A proposed application for these nanolasers is in on-chip optical interconnects. In this application the modulation bandwidth is important. We have calculated the modulation bandwidth using a small-signal analysis of the rate equations, which has the advantage compared the full model of yielding a semianalytical result. The results are shown in Fig. 1c). The nanolaser has a high bandwidth, which is mostly limited by the cavity lifetime (/ 1/), but shows a prominent drop during the lasing regime. The modulation bandwidth again becomes limited by as the output becomes chaotic.

•Jan 1, 2012•0 citations

Stochastic and probabilistic equations for three- and four-level lasers: tutorial

Colin J. McKinstrie

In this tutorial, Langevin stochastic equations and Markov probability equations are used to model electron- and photon-number fluctuations in three- and four-level lasers. Equations are derived for the moments of the electron and photon numbers (means, variances and correlations). Both approaches produce the same moment equations. In the Langevin approach, the moments of the noise terms must be specified by other means, whereas in the Markov approach, they are determined self-consistently and satisfy the shot-noise rule: For each process that is modeled by the rate equations, the driving terms in the variance equations equal the moduli of the associated terms in the mean equations. The driving terms in the correlation equations have the same magnitudes as the variance terms, but can be positive or negative, depending on whether the changes in the electron and photon numbers are correlated or anti-correlated, respectively. Formulas are derived for the relative intensity noise and its spectrum. The consequences of these results for three- and four-level lasers are discussed.

Journal of the Optical Society of America B•Feb 3, 2020•7 citations

Collective Effects in Nanolasers Explained by Generalized Rate Equations

Игорь Евгеньевич Проценко, Emil C. André, Alexander V. Uskov +2 more

We study the stationary photon output and statistics of small lasers. Our closed-form expressions clarify the contribution of collective effects due to the interaction between quantum emitters. We generalize laser rate equations and explain photon trapping: a decrease of the photon number output below the lasing threshold, derive an expression for the stationary cavity mode autocorrelation function $g_2$, which implies that collective effects may strongly influence the photon statistics. We identify conditions for coherent, thermal and superthermal radiation, the latter being a unique fingerprint for collective emission in lasers. These generic analytical results agree with recent experiments, complement numerical results, and provide insight into and design rules for nanolasers.

arXiv:1709.08200•Sep 24, 2017•7 citations
The random nature of spontaneous emission leads to unavoidable fluctuations in a laser's output. This is often included through random Langevin forces in laser rate equations, but this approach falls short for nanolasers. In this paper, we show that the laser quantum noise can be quantitatively computed for a very broad class of lasers by starting from simple and intuitive rate equations and merely assuming that the number of photons and excited electrons only takes discrete values. The success of the model is explained by showing that it constitutes a Markov chain, which can be derived from the full master equations. We show that in the many-photon limit, the model simplifies to Langevin equations. We perform an extensive comparison of different approaches for computing quantum noise in lasers, identifying the best approach for different system sizes, ranging from nanolasers to macroscopic lasers, and different levels of excitation, i.e., cavity photon number. In particular, we find that the numerical solution to the Langevin equations is inaccurate below the laser threshold, while the laser Markov chain model, on the other hand, is accurate for all pump values and laser sizes when collective emitter effects are excluded.
It is shown that the quantum intensity noise of a single-emitter nanolaser can be accurately computed by adopting a stochastic interpretation of the standard rate equation model. The only assumption made is that the emitter excitation and photon number are stochastic variables with integer values. This extends the validity of rate equations beyond the mean-field limit and avoids using the standard Langevin approach, which is shown to fail for few emitters. The model is validated by comparison to full quantum simulations of the relative intensity noise and second-order intensity correlation function, g^{(2)}(0). Surprisingly, even when the full quantum model displays vacuum Rabi oscillations, which are not accounted for by rate equations, the intensity quantum noise is correctly predicted by the stochastic approach. Adopting a simple discretization of the emitter and photon populations, thus, goes a long way in describing quantum noise in lasers. Besides providing a versatile and easy-to-use tool for modeling emerging nanolasers, these results provide insight into the fundamental nature of quantum noise in lasers.
Laser amplification kicks off with just a few photons in the cavity, a regime where standard additive Langevin noise cannot capture the real fluctuations. To understand how this amplification begins, we map out a single-mode photon-carrier model that links macroscopic lasers directly into nanolasers. Stationary rate equations reveal that stimulated emission takes over across a tight pump window. In nanolasers, however, spontaneous processes coexist across a much broader range. Because steady-state averages miss the individual events that lead to the macroscopic photon number, we frame the stochastic dynamics around the actual discrete transitions, deriving a state-dependent diffusion approximation. The destruction of an excited carrier concurrent with the generation of a photon creates a negative photon-carrier cross-correlation. This unavoidably leads to a correlated, multiplicative Gaussian-noise framework. We track the resulting discrete dynamics with an established, event-based Laser Markov Chain (LMC) which stochastically samples pump, relaxation, emission, and cavity-loss events keeping all physical quantities as non-negative integers. Comparing its predictions to those of our multiplicative model proves that the latter reproduces excellently well the statistical behavior of photons and carrier, as long as a continuous diffusion description holds. Once the emission drops into the few-photon regime, however, the discrete model becomes indispensable. On the other hand, at large photon counts relative fluctuations shrink to a level where the state-dependence of the diffusion matrix can be ignored, to leading order, thus returning to the classic Langevin description of the system's stochastic dynamics.
Using two equivalent methods, the Heisenberg-Langevin and the master equation, we investigate the properties of a single-emitter nanolaser: an incoherently pumped single two-level system interacting with a single cavity mode of finite finesse. We are interested in a stationary generation regime when the mean number of photons inside the laser cavity is much larger than unity. In this case we linearize the Heisenberg-Langevin equations around the stationary classical solution and obtain analytical results for the linewidth, the quadrature fluctuation spectra, and the intracavity Mandel QQQ parameter. As the second method we use the master equation for the density operator of the laser field. Using the coherent-state representation of the density matrix we obtain the stationary solution for the Glauber PPP function of the laser mode and investigate the quantum statistics of the laser field and its threshold behavior. We compare the results of both methods with numerical simulations and find very good agreement. In the end we use our theoretical results for analyzing experiment with single-ion nanolasers.
The laser has evolved from table size apparatuses to truly nano sized devices, in much the same way that computer chips have been continuously minimized. The few-emitter nanolaser represents an extreme in terms of size. The emitters can be either atoms or quantum dots, and they are coupled to a high-Q optical cavity. This system has previously been studied for the one and two emitters case, see e.g. [1, 2]. The emitters and cavity are modelled as two-level systems and a harmonic oscillator, respectively. The coupled system is modelled using the Jaynes-Cummings Hamiltonian and the two-level systems are pumped incoherently by a rate P. Solutions are found using the corresponding master equation. However, with cavity populations exceeding 100 and several emitters, the dimension of the Hilbert space of the system becomes too large to handle efficiently on a conventional computer. E.g. for four emitters and 100 photon states the density matrix has more than 2.5 × 106 elements. We have been able to simplify the problem significantly by adiabatically eliminating the photon-assisted polarizations and the correlations between emitters and cavity [3]. This results in a set of rate equations for the population of the cavity, na, and the occupation of the emitters. Fig. 1a) shows na for up to 4 emitters coupled to a cavity, as a function of the pumping rate P. The figure also shows results from a full master equation solution, and the correspondence is very good for large values of P. In Fig. 1b) the second order correlation functions g(2)(0) are shown, as obtained using the full model. The second order correlation functions become 1 for sufficiently large values of P signifying the onset of lasing. However, for larger values of P the laser quenches and the emission becomes chaotic (g(2)(0) = 2). This quenching effect is well-known in two-level models [3], but might not be realistic for semiconductor quantum dots. A proposed application for these nanolasers is in on-chip optical interconnects. In this application the modulation bandwidth is important. We have calculated the modulation bandwidth using a small-signal analysis of the rate equations, which has the advantage compared the full model of yielding a semianalytical result. The results are shown in Fig. 1c). The nanolaser has a high bandwidth, which is mostly limited by the cavity lifetime (/ 1/), but shows a prominent drop during the lasing regime. The modulation bandwidth again becomes limited by as the output becomes chaotic.
In this tutorial, Langevin stochastic equations and Markov probability equations are used to model electron- and photon-number fluctuations in three- and four-level lasers. Equations are derived for the moments of the electron and photon numbers (means, variances and correlations). Both approaches produce the same moment equations. In the Langevin approach, the moments of the noise terms must be specified by other means, whereas in the Markov approach, they are determined self-consistently and satisfy the shot-noise rule: For each process that is modeled by the rate equations, the driving terms in the variance equations equal the moduli of the associated terms in the mean equations. The driving terms in the correlation equations have the same magnitudes as the variance terms, but can be positive or negative, depending on whether the changes in the electron and photon numbers are correlated or anti-correlated, respectively. Formulas are derived for the relative intensity noise and its spectrum. The consequences of these results for three- and four-level lasers are discussed.
We study the stationary photon output and statistics of small lasers. Our closed-form expressions clarify the contribution of collective effects due to the interaction between quantum emitters. We generalize laser rate equations and explain photon trapping: a decrease of the photon number output below the lasing threshold, derive an expression for the stationary cavity mode autocorrelation function g2g_2g2​, which implies that collective effects may strongly influence the photon statistics. We identify conditions for coherent, thermal and superthermal radiation, the latter being a unique fingerprint for collective emission in lasers. These generic analytical results agree with recent experiments, complement numerical results, and provide insight into and design rules for nanolasers.