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:2610.25520

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

G. Sonnino, Ph. Peeters, P. Nardone, G. L. Lippi

quant-phphysics.opticsSubmitted Oct 18, 2026• Announced Oct 20, 2026
arXivHTMLPDF

Abstract

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.

Comments
15 pages, 8 figures

Related papers

Modeling Quantum Noise in Nanolasers using Markov Chains

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

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. While the approach has seen previous success, we here derive it rigorously from an open quantum system master equation, whereas it was previously introduced only on phenomenological grounds. We further 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 show that below the laser threshold, stochastic fluctuations in the numerical solution to the Langevin equations can drive populations to unphysical negative values, requiring the introduction of population bounds, which in turn skew the noise statistics, leading to inaccuracies. The Laser Markov Chain model, on the other hand, is accurate for all pump values and laser sizes when collective emitter effects are excluded.

quant-pharXiv:2511.13622•Physical Review A•Nov 18, 2025•0 citations

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

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

Coherence buildup and laser thresholds from nanolasers to macroscopic lasers

Mark Anthony Carroll, Giampaolo D’Alessandro, Gian Luca Lippi +2 more

We detail the derivation of nanolaser models that include coherent and incoherent variables and predict the existence of a laser threshold, irrespective of cavity size and emitter number, for both single- and multielectron systems. The growth in photon number in the lasing mode is driven by an increase in correlation between absorption and emission processes, leading to the onset of self-sustained stimulated emission (laser threshold), followed in turn by a correlation decrease and ending with the dominance of coherent emission. The first-order coherence ${g}^{(1)}$ steadily increases, as the pump grows towards the laser threshold value, and reaches unity at or beyond threshold. The transition toward coherent emission becomes increasingly sharp as the numbers of emitters and coupled electromagnetic cavity modes increase, continuously connecting, in the thermodynamic limit, the physics of nano- and macroscopic lasers at threshold. Our predictions are in remarkable agreement with experiments whose first-order coherence measurements have so far been explained only phenomenologically. A consistent evaluation of different threshold indicators provides a tool for a correct interpretation of experimental measurements at the onset of laser action.

arXiv:2302.08824•Physical Review A•Jun 26, 2023•10 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

Efficient stochastic simulation of rate equations and photon statistics of nanolasers

Emil C. André, Jesper Mørk, Martijn Wubs

Based on a rate equation model for single-mode two-level lasers, two algorithms for stochastically simulating the dynamics and steady-state behaviour of micro- and nanolasers are described in detail. Both methods lead to steady-state photon numbers and statistics characteristic of lasers, but one of the algorithms is shown to be significantly more efficient. This algorithm, known as Gillespie's First Reaction Method (FRM), gives up to a thousandfold reduction in computation time compared to earlier algorithms, while also circumventing numerical issues regarding time-increment size and ordering of events. The FRM is used to examine intra-cavity photon distributions, and it is found that the numerical results follow the analytics exactly. Finally, the FRM is applied to a set of slightly altered rate equations, and it is shown that both the analytical and numerical results exhibit features that are typically associated with the presence of strong inter-emitter correlations in nanolasers.

•Oct 26, 2020•17 citations
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.
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. While the approach has seen previous success, we here derive it rigorously from an open quantum system master equation, whereas it was previously introduced only on phenomenological grounds. We further 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 show that below the laser threshold, stochastic fluctuations in the numerical solution to the Langevin equations can drive populations to unphysical negative values, requiring the introduction of population bounds, which in turn skew the noise statistics, leading to inaccuracies. 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.
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.
We detail the derivation of nanolaser models that include coherent and incoherent variables and predict the existence of a laser threshold, irrespective of cavity size and emitter number, for both single- and multielectron systems. The growth in photon number in the lasing mode is driven by an increase in correlation between absorption and emission processes, leading to the onset of self-sustained stimulated emission (laser threshold), followed in turn by a correlation decrease and ending with the dominance of coherent emission. The first-order coherence g(1){g}^{(1)}g(1) steadily increases, as the pump grows towards the laser threshold value, and reaches unity at or beyond threshold. The transition toward coherent emission becomes increasingly sharp as the numbers of emitters and coupled electromagnetic cavity modes increase, continuously connecting, in the thermodynamic limit, the physics of nano- and macroscopic lasers at threshold. Our predictions are in remarkable agreement with experiments whose first-order coherence measurements have so far been explained only phenomenologically. A consistent evaluation of different threshold indicators provides a tool for a correct interpretation of experimental measurements at the onset of laser action.
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.
Based on a rate equation model for single-mode two-level lasers, two algorithms for stochastically simulating the dynamics and steady-state behaviour of micro- and nanolasers are described in detail. Both methods lead to steady-state photon numbers and statistics characteristic of lasers, but one of the algorithms is shown to be significantly more efficient. This algorithm, known as Gillespie's First Reaction Method (FRM), gives up to a thousandfold reduction in computation time compared to earlier algorithms, while also circumventing numerical issues regarding time-increment size and ordering of events. The FRM is used to examine intra-cavity photon distributions, and it is found that the numerical results follow the analytics exactly. Finally, the FRM is applied to a set of slightly altered rate equations, and it is shown that both the analytical and numerical results exhibit features that are typically associated with the presence of strong inter-emitter correlations in nanolasers.