## 1. Introduction

Amplified Spontaneous Emission (ASE) is a pervasive phenomenon ever present in any optical system which exhibits gain. Its origin lies in the existence of the induced emission, symmetric of absorption, postulated by Einstein [

1] in his seminal paper on radiation. In the field of lasers and optical amplifiers, which has grown a few decades later, it has progressively gained importance as a competitor [

2], precursor [

3] or substitute [

4] of lasing and broadband noise source [

5] or pulsed source in optical amplifiers (cf. [

6] and references therein for numerous applications).

From the laser point of view, fundamental aspects of ASE begin with the realization of the existence of a threshold for the spontaneous amplification [

7,

8]. Overcoming this threshold amounts to enabling the amplification of a signal at the end of a pumped sample to be spontaneously amplified through propagation in the medium—whence the term ASE. In such a case a “mirrorless laser” can be built [

9], and one of its early experimental realization is the (transversely excited, high pressure)

${N}_{2}$ laser [

10,

11]. As already mentioned, and to be expected, systems with high gain—which make for excellent oscillators—have been very well sought after for their ease in realizing a laser [

12,

13], as well as the more versatile dye-based counterparts, for which a thorough characterization of the combined ASE and coherent emission has been conducted [

14] (corrections to initially established relations, mostly applicable to solid state lasers, have been later determined [

15]). An alternative representation of ASE, seen as amplified noise, has been also proposed [

16].

At the macroscopic scale, ASE is responsible for astrophysical laser sources [

8,

17], on the one hand, and has been exploited to strongly improve the efficiency of the Free Electron Laser, using a self-amplification regime [

18,

19,

20], or for UV amplification [

21]. ASE, exploited as a biocompatible nonlinear, scattering induced amplifier [

22], is the natural precursor of random lasing [

23,

24] and is predicted to appear in concomitance with a reduction in the number of degrees of freedom in such systems [

25].

From the point of view of field monochromaticity, a large amount of ASE is detrimental in conjuction with coherent emission in bulk lasers [

2] or fiber amplifiers [

26]. However, it enables pulse shortening [

27] and very broadband lasers and amplifiers [

28], whose current application is the generation of ultrashort pulses [

29] approaching the attosecond regime. In this sense, ASE is strongly beneficial to the realization of sources with extreme performance. Current interest, at an entirely different lengthscale, covers the micro- and nanoscopic regime, where on the one hand solid-state microdisk lasers suffer from the incoherent emission component [

30], while ASE is used as a measure of the potential for obtaining lasing action in nanodevices [

31,

32] even before achieving actual lasing—its usefulness is therefore in validating the selection of candidates for nanolasing emission. Strong technological interest in ASE appears in the investigation of colloidal systems [

33] and their threshold engineering [

34] as well as in organic lasers [

35,

36] as reviewed in [

37].

A very interesting regime in superluminescent diodes has been identified, where the superpossonian statistics approaches the lasing regime as the current supplied to the device is increased [

38]. In addition to the potential technological interest proposed in the original paper, this observation may suggest the emission of frequent, ultrashort ASE pulses which amount to a quasi-continuous emission. Whether this is a generic kind of emission is a question which is now of great interest in relation with photon condensation. Without entering into a detailed discussion of this line of research, which escapes the scope of this contribution, we highlight this emerging topic, susceptible of leading to a more complete description of the physics of radiation in its evolution from the entirely spontaneous emitter relaxation into the purely coherent (lasing) emission.

The observation of the condensation of photons in interaction with a resonant medium in a cavity [

39], and of its different statistical and non-stationary properties [

40,

41], has prompted the investigation of its occurrence even in other systems [

42,

43], hinting to a phenomenon more general than what initially considered (even without considering polariton and excitons which lie outside the realm of this discussion). A generalized thermodynamical treatment of Kirchoff’s law [

44] enables the description of subsystems with different temperatures—thus out of equilibrium—bridging the gap between the traditional thermal point of view and more modern devices. Its applications are numerous and, in addition to enabling more efficient, tailored coupling between emitters and space, offers one more step in the description of the radiation field starting from the ideal spontaneous emission (one single emitter in vacuum) all the way to the non-equilibrium situation, eventually represented by the laser [

45]. Potential interest for futuristic high-speed communications rests on the development of these concepts [

46].

A complete overview of spontaneous emission is outside the scope of the present contribution; however, aspects which are important for our scopes include the filtering action of a cavity [

47] or the relation between spontaneous emission and gain in Quantum Well semiconductor devices [

48]. Fundamental features of the relationship between the so-called “Dicke superradiance” and ASE are discussed in [

49,

50], and provide the foundations of the transition between the spontaneous and the stimulated regimes. Finally, it is worth pointing out recent investigations on the coherence properties of ASE and the possibility for obtaining lasing action—through ASE—in a cavity-less system [

51,

52]. A good overview of the physics of small-sized lasers and of many of their threshold features—including discussions on noise properties and quantum-mechanical considerations—can be found in [

53] (and references therein).

In this contribution, we concentrate on pre-threshold photon bursts that have been experimentally observed in micro- and nanolasers, whose existence can be (mostly) predicted by a stochastic modelling approach based on Einstein’s semiclassical theory of radiation [

1]. At the present stage, no definite answers exist on their actual nature, aside from the fact that they appear unusual from a traditional laser physics point of view. We compare observed and predicted features, and comment on the technological limitations which constrain their experimental detection. Topological considerations, based on standard laser models, provide guidance in interpreting the observations and in reconstructing their physical origin. Second-order zero-delay autocorrelations are the most widespread means of characterizing very weak signals, such as those emitted by nanolasers, and have been widely employed also for thermal or quasi-thermal sources [

38]. A simple analysis of the autocorrelation of a pulsed source complements the guidance provided by the topological analysis and helps securing a grip onto some of the observed characteristics of pre-threshold photon bursts.

This paper is work in progress and aims at proposing more questions than it can provide answers. As such, it is possible that several Readers may have diverging views from the writer, but also that the latter may revise his position on some of the exposed points in the near or far future. Thus Readers should take this contribution not so much as a finished product, but, rather, as a stimulus to delve more deeply in the points that are raised and in the much more numerous ones which are only hinted upon as related issues.

The paper touches on different aspects of photon bursts, directly and indirectly observed in micro- and nanolaser experiments, predicted by fully stochastic models and (mostly) missed by differential models. The presentation is divided into 6 main sections, followed by conclusions (

Section 8). The question of the nature of a laser threshold in micro- and nanolaser is discussed in

Section 2 starting from the proof of its equivalence to a phase transition for macroscopic lasers and leading into the paradoxical case of the

thresholdless nanolaser. The experimental constraints, whose origin resides in the technological limitations of current instrumentation, are highlighted together with the potentially crucial role that models may play in directing the investigations on this topic. Since photon bursts appear in the smooth transition region from incoherent to coherent emission, this point is central to the discussion. The (scarce) collected evidence on photon bursts is summarized in

Section 3, with reference to the publications in which this is found, and a physical role in the incoherent→coherent transition is proposed (

Section 3.2). Since direct measurements of the temporal evolution of laser emission are possible only with microdevices, most of the information is gathered through the use of autocorrelation functions obtained through photon counting. A simple model,

Section 4, is used to show that the presence photon bursts can be statistically identified through the appearance of superthermal statistics (of course in the single-mode emission regime, as it is well-known that mode competition leads to strong bunching). The model, based on a sequence of pulses of variable amplitude sitting on an adjustable background, is used to show the limitations of this approach, where superthermal bunching disappears when the background grows too large or when the pulses become too frequent. A brief overview of the two main classes of models (

Section 5), differential and stochastic (

Section 5.2), is given, summarizing their main characteristics and their ability (

Section 5.1) to include some of the crucial aspects for a thorough description of the physics of laser threshold in small devices. The predictions of different implementations of stochastic models are summarized in

Section 5.3 in relationship to photon bursts. Up until here, the paper summarizes observations (both experimental and numerical) and examines the characteristics of models or analysis tools in relationship to their potential for describing photon bursts. Some clearer understanding, accompanied by predictive power and insight into the dynamics, is offered by a topological analysis of the phase space which underlies the description of laser physics is presented in

Section 6. Here, the eigenvalues emerging from the linear stability analysis are examined for nano-, micro-, and macrolasers, together with their accompanying eigenvectors. This analysis enables a clearer view of the stability of each kind of laser and of the possibility for adequate noise to give rise to photon bursts. One last section examines the features of a phenomenological modifications of the Rate Equations, very recently proposed, where photon bursts have been observed: a first in the ability for a differential model in predicting the emission of large photon bunches. Through the topological analysis (

Section 7), we see how the proposed modification may enable the more modest amount of noise which accompanies the differential description to give rise to photon bursts.

The schematic representation of

Figure 1 shows in graphical form the different topics which are going to be addressed in the discussion. On a typical nanolaser response curve (black line) which shows the photon number as a function of pump, we qualitatively indicate the

threshold region, in which photon bursts are observed (stochastically computed dark green trace). A physical interpretation on their origin will be offered on the basis of the phase space properties, showing that before continuous coherent emission any fluctuation (purple arrow) possesses a small likelihood of relaxing into the photon field, and thus easily grows extremely large, with the consequent emission of a photon burst (in a qualitatively similar way of what is achieved in laser Q-switching). Along the way, we examine the experimental aspects of the problem on the basis of existing results, but also considering the experimental limitations and offering some insight into the main analysis tool: the second order autocorrelation. The current models’ predictive abilities in regard to photon bursts are also reviewed, considering both differential and stochastic methods.

In the following, we will refer to spontaneous emission in the traditional sense used in laser physics, rather than in the thermal radiation sense partly discussed in this introduction. The term laser threshold is used to denote the true threshold point, as defined in a phase transition (but also in bifurcation theory), but also to identify the region in which the electromagnetic field evolves from the emission of (mainly) incoherent photons to the coherent regime. Threshold region is used when emphasizing the extent of the transition in pump range.

## 2. Laser Threshold at the Micro- and Nanoscale

The question of laser threshold, a non-issue in the macroscopic devices used for decades, has come to the forefront with the considerable cavity volume reduction introduced by the Vertical Cavity Surface Emitting Laser (VCSEL) design [

54], which led, in less than two decades, to ultralow currents for coherent emission [

55]. The impressive reduction in cavity volume achieved by these devices brought a substantial increase in the fraction of spontaneous emission coupled into the lasing mode (

$\beta $ factor) with the consequent linear relationship between input and output for

$\beta =1$. As a result, based on the traditional macroscopic laser point of view, threshold was extrapolated to a zero pump value; the devices fitting this condition became therefore known as

thresholdless lasers [

56] or

zero-threshold lasers [

57,

58].

While in macroscopic lasers threshold is compatible with a phase transition definition [

59,

60,

61,

62], with nearly unmeasurable deviations from the thermodynamic limit, the strong finite size effects typical of nanocavities introduce a smooth change from the non-lasing to the lasing solution [

56,

63] and open questions about the nature of the transition [

64]. The extended pump range over which the transition takes place demonstrates an

unfolding of the thermodynamical discontinuity and the suggest the name of threshold region for the corresponding pump interval. The practical consequence of the uncertainty surrounding a clear definition of threshold in small devices has led to the call for more stringent criteria to identify lasing [

65], particularly in the technologically very active field of organic devices [

66]. From a quantum-statistical point of view, the definition of threshold is legitimate only in the thermodynamic limit and for devices which approach such a regime (thus macroscopic lasers, only), as clearly demonstrated in a seminal paper [

67]. However, the rigorous analysis does not provide any useful answers to the practical question of whether and when does the emission of a small source become (sufficiently) coherent, at least to be accepted as a laser. This question carries far reaching consequences, since many of the applications for which lasers are selected require a (sufficiently) good level of coherence, which has to be ensured for their usefulness. Luckily, statistical theories are not the only tools to define a transition and a very recent contribution has shown that through bifurcation analysis it is possible—in a suitably extended model (cf.

Section 5.1)—to predict the existence of a true laser threshold even in a thresholdless laser [

68].

In spite of recent advances, a clear description of the physics of the radiation emitted by a nanolaser in the transition from incoherent to coherent emission is still lacking. Two concomitant reasons contribute to the uncertainty. On the one hand, measurements are very difficult in the threshold region due to the extremely reduced photon flux which, coupled to extremely short evolution timescales, imposes a currently unsurmountable double constraint on instrumentation: high sensitivity and large bandwidth. Experimental schemes to circumvent the limitations vary from devising alternative (somewhat larger) sources with (hopefully) similar (enough) characteristics [

69], to better detectors [

70] or to more advanced correlation techniques [

38,

71]. However, none can entirely compensate for the required extreme performance, and the information that can be experimentally gathered remains partial. On the other hand, models have so far failed to provide a clear picture of the physics of radiation in the threshold region, leaving the more detailed predictions to fully numerical, stochastic approaches. The compounded experimental and theoretical difficulties have left the investigations in a state where progress is slow and difficult. A brief discussion of the current state of modelling is offered in

Section 5.

In this contribution we focus our attention onto the region which precedes coherent emission (even a very noisy one, with a broadband radiofrequency spectrum [

72]) where we observe pulses whose features resemble those of ASE. So far the number of reports concerning this kind of emission at the micro- and nanoscale is relatively small, nonetheless it has been seen—directly or indirectly—in different kinds of devices and at different laser scales. An interesting question, for future consideration, concerns the possible relationship between this regime of emission and what can be expected from thermal sources in a condensed photon phase. From a practical point of view, spontaneous photon bursts may bring along alternative ways of encoding information in small devices at lower power consumption, provided suitable schemes are employed [

73].

## 4. Relationship between Bursts and ${\mathit{g}}^{\left(\mathbf{2}\right)}\left(\mathbf{0}\right)$—Experimental Considerations

Some simple mathematical considerations help highlighting some of the issues which connect the (indirect) observations of photon bursts with the zero-delay second-order autocorrelation. Since the basic properties that we intend to point out are general, we can base our observations on a simple model where we compute the autocorrelation of a single rectangular pulse defined by

where

T represents the period,

A is the amplitude of the signal and

$\alpha $ represents the fraction of the maximum

A by which the bottom level (background) is raised. Its zero-delay second-order autocorrelation, defined by [

85]

where

$\langle \xb7\rangle $ stands for the temporal average, straightforwardly gives the following expression:

Two limits can be readily obtained:

$\alpha =0$ and

$\alpha =1$, which correspond to the disappearance of background and disappearance of the peak, respectively. The former gives

while from the latter it is immediate to see that

Thus, in the total absence of a background (unlikely in an experiment, where there is at least some detection noise) the autocorrelation is inversely proportional to the width of the pulse and can diverge as this quantity goes to zero. Although the limit does not hold in a real system, it is easy to see how a single pulse can give rise to superthermal statistics, which occurs every time $\delta <0.5$. When, instead, the pulse disappears ($\alpha =1$) the statistics becomes that of Poisson, or coherent radiation.

From an experimental point of view, the intermediate cases are most interesting, especially because the number of pulses may be variable and background is unavoidable. First of all, the consideration of a single pulse in the period T (and consideration of a periodic function itself) is not a restriction. In practice, the observation window is of finite size, which can be assumed to be the duration T (periodicity is not a necessary hypothesis in the computation we have performed). The presence of more than one pulse (as long as we consider them to be all the same height) amounts to changing the fraction $\delta T$ of the time interval occupied by the pulses themselves. The computation is insensitive to the presence of N pulses of duration $\delta T$ or of a single pulse of duration $N\delta T$, as long as $N\delta \le 1$. The pulses can even be different in duration, what matters is only the fraction of the time interval T which they fill.

Figure 3 clearly shows the inverse relationship with the duty cycle in the absence of an offset (baseline

$\alpha =0$, black line) and the degradation of the expected value of

${g}^{\left(2\right)}\left(0\right)$ as a function of the duty cycle

$\delta $. It is interesting to notice that

${g}^{\left(2\right)}\left(0\right)$ decreases rapidly as the pulse width (or number of pulses) progressively fill the (measurement) time interval

T and fall below 2 (subthermal statistics) when

$\delta \ge 0.5$, as already stated. More interesting is the fact that the background (noise) rapidly reduces the autocorrelation, which becomes subthermal for any value of

$\delta $ when

$\alpha \gtrsim 0.2$ (magenta line in

Figure 3). While in itself there is no particular reason for preferring

${g}^{\left(2\right)}\left(0\right)>2$, its practical usefulness is clear, since it helps identifying a regime of photon bursts in experiments where

${g}^{\left(2\right)}\left(0\right)$ is the only easily accessible quantity. Whenever the statistics becomes subthermal it is impossible to distinguish between a regime with burst emission and a dynamically different one (photon bunched in a subthermal way, for instance, dynamical oscillations [

72,

86], or background effects [

84]). This simple example immediately shows the limitations of the autocorrelation measurement in the identification of bursts and highlights the difficulties inherent in the presence of too large a number of pulses (fractional filling of the measurement interval

T) and of the influence of background.

The picture does not substantially change, up to numerical factors (e.g., in the value of

$\delta $ and

$\alpha $ for which the statistics becomes subthermal), if the pulses are of unequal height. In practice, experiments will also suffer from bandwidth limitations which are going to broaden the pulses, an instrumental factor which contributes to reducing the actual value of

${g}^{\left(2\right)}\left(0\right)$. Amplitude reduction, due to detector sensitivity, also contributes to a degradation of

${g}^{\left(2\right)}\left(0\right)$, especially since the background may not be affected (if due to electronics). Overall, these factors compound and it is not uncommon to observe experimentally measured functional dependences in the autocorrelation where a growth from shot noise (

${g}^{\left(2\right)}\left(0\right)=1$) towards thermal statistics is observed as the pump is increased, followed by a decrease towards Poisson statistics (not necessarily reached, cf. e.g., Figure 2c in [

87]). The subthermal value of

${g}^{\left(2\right)}\left(0\right)$ at low pump is to be ascribed to the strong filtering of the instrumentation (and insufficient sampling, also observed in some stochastic numerical simulations,

Section 5.2). Of course, such functional dependences cannot be univocally interpreted: it is impossible to tell whether the failure to reach thermal statistics is simply due to instrumental limitations or whether a burst regime remains undetected, due to possible confusion with a noisy but coherent signal. Indeed, a subthermal, but superpossonian, regime on the decreasing slope of

${g}^{\left(2\right)}\left(0\right)$ may be due to statistical properties of the radiation, to dynamics, to frequent photon bursts which fill the observation window, or to instrumental background: the multiplicity of causes prevents a univocal interpretation of the observation. It is, however, interesting to notice that this simple picture may help interpret the observation of

${g}^{\left(2\right)}\left(0\right)<2$ in a regime of ASE with decreasing values as the power supplied to the superluminescent diode is increased [

38]: frequent, short ASE pulses may very well produce autocorrelation values compatible with the observations.

Table 1 summarizes the main results of this section, by identifying the main consequence of each feature and explaining the influence that it has on the measurement. These constraints are crucial in experiments, even though they control the computation of the autocorrelation even from signals obtained from a numerical integration. In the latter case, however, most of the constraints can be lifted (up to the difficulty of obtaining sufficient statistics below threshold), thus we focus on the experimental side.

As already explained, photon bursts can occur with a subthermal

${g}^{\left(2\right)}\left(0\right)$ (cf. e.g., Figure 4b in [

76]), but similar autocorrelation values can also be obtained with strong oscillations (Figure 4c in [

76]). The lack of bijectivity prevents all conclusions. Frequent bursts increase

$\delta $ (

Figure 3), while a degradation in the signal-to-noise ratio (or signal-to-background ratio), represented by increasing values of

$\alpha $ (e.g., along a vertical cut in

Figure 3), reduces again

${g}^{\left(2\right)}\left(0\right)$. Not surprisingly, in the detection system high sensitivity and low noise are therefore crucial, a difficult set of requirements to match for any kind of micro- or nanolaser. For simplicity, the table has been prepared with a detection scheme in mind where a linear detector samples the full temporal signal with sufficient resolution. In photon counting the requirements are the same; there, the signal-to-noise ratio is represented by the dark count rate and the integration window (number of samples required to obtain a good signal on top of noise). The latter can often be large (overcoming very easily hundreds of thousands of measurements) and require long measurement times, with the consequent stability issues for the experiment. The longer the time windows, the stronger the perturbations, thus the difficulties are equivalent to those outlined in the table, even though they may appear in a different form. The same holds for measurements carried out far below threshold, where the occurrence of a burst is extremely unlikely: the accumulated background noise takes over and the autocorrelation is governed by the shot noise (cf. e.g., Figure 2c in [

87]). Finally, in photon counting the detection bandwidth is replaced by the finite detector averaging time.

## 6. Phase Space Information

While only a proper stochastic simulation is currently capable of reproducing the observed dynamics [

84], a topological investigation of the phase space properties of the rate equations, with the added contribution of the average spontaneous emission, provides useful information. This is not surprising, given that the underlying topological structure on which the random walk takes place (e.g., [

67,

114]) is the same as the one of the stochastic laser simulators [

83] (both derive from the semiclassical description of radiation [

1]).

The traditional rate equations with the spontaneous emission contribution to the lasing mode read for a Quantum Well system [

96] (equivalent to a “four-level” laser model):

where

n and

N represent the photon and carrier number (or population inversion), respectively,

${\mathsf{\Gamma}}_{c}$ and

$\gamma $ are the relaxation rates for the intra-cavity photons and for the population inversion, respectively,

P is the pump rate and

$\beta $ is the previously defined fraction of spontaneous emission coupled into the lasing mode. This set of REs includes the average contribution of the spontaneous emission to the number of coherent photons in the cavity mode through the term

$\beta \gamma N$ (Equation (

7)). The steady state (equilibrium) values for the photon number and the population inversion, as directly obtained from Equations (

7) and (8), are:

where the over-strike represents equilibrium.

The linear stability analysis (lsa) [

97] of Equations (

7) and (8), performed around the steady-state values for the two variables, is readily obtained from:

where we have defined the perturbations

The full analysis is performed here for

${\mathsf{\Gamma}}_{c}={10}^{11}{s}^{-1}$ and

$\gamma ={10}^{9}{s}^{-1}$. The “threshold”, Equation (11), is the one defined in [

67].

#### 6.1. Eigenvalue Analysis

The

lsa provides information on the laser’s return (or lack thereof) to its stationary emission state when a perturbation is applied to its two variables. The eigenvalues give a measure of the state’s robustness and are shown in

Figure 4 for different values of

$\beta $.

For ease of comparison, the horizontal scale plots the pump rate relative to the threshold value (Equation (11)) for the matching laser (i.e., normalized to

$\beta $). The first observation is that below a critical pump rate, there are two distinct real eigenvalues; however, they collapse into a pair of complex conjugate eigenvalues for instance at

$\frac{P}{{P}_{thr}}\approx 2.1$ for

$\beta =1$ (black curve). From this point on, the real parts of the two eigenvalues coincide, since they distinguish themselves only for their imaginary parts (not shown). There are two important remarks which arise from this figure: (1) as the fraction of spontaneous emission coupled into the lasing mode decreases, the eigenvalues become complex conjugate at an earlier (normalized) pump rate; (2) as

$\beta $ decreases, the eigenvalue closer to the real axis becomes progressively less negative, thus indicating that the corresponding laser is less stable. Notice that the traditional rate equations model devoid of the spontaneous emission contribution predicts one zero eigenvalue at threshold, matching the indifferent stability responsible for the divergence of fluctuations and for the presence of a phase transition [

59,

60,

61,

62]. This corresponds to a device where

$\beta \to 0$. Here, the finiteness of

$\beta $ stabilizes the solution, giving rise to a so-called

imperfect bifurcation [

99] whose features have a certain resemblance to those of an

avoided level crossing in quantum mechanics [

120].

Figure 4 clearly shows that the stabilization effect is much stronger at the nanoscopic than at the microscopic scale, since a more negative eigenvalue represents a faster convergence to 0 of the corresponding perturbation (Equations (

14) and (15)).

As an aside, it is interesting to notice the displacement of the threshold, compared to the standard predicted value [

67], which corresponds to 1 on the (normalized) pump axis (

Figure 4, from Equation (11)). For

$\beta \ge {10}^{-2}$ the eigenvalues become complex conjugate at

$P>{P}_{th}$, i.e., in the upper region of the characteristic S-curve of laser emission (cf.

Figure 5). The new quantum-mechanical model [

68] also predicts the appearance of the coherent field bifurcation close to the upper branch, signalling a displacement of the threshold towards larger pump values. This strengthens the case for a strong role of photon bursts in the transition between spontaneous and coherent emission. Extensions based on the statistical theory of phase transitions, instead, fail, as they predict a threshold at the mid-point of the steeper region of growth (i.e.,

$\frac{P}{{P}_{th}}=1$) [

67]).

Figure 6 shows the equivalent eigenvalues for a borderline microlaser (

$\beta ={10}^{-4}$) and for a macroscopic laser (

$\beta ={10}^{-6}$), even though not a very large one (a large semiconductor, for instance, or a solid state microdevice). The difference is striking, as in the macrosopic laser one eigenvalue comes much closer to zero than in all other devices, proving the progressive convergence towards the thermodynamic limit. As in

Figure 4, the microlaser displays a (small) dip at pump values close to “threshold”, while the macroscopic one reduces its stability in this region, inverting the trend observed at the micro- and nanoscale and approaching the thermodynamic limit.

While the lessening of the stability with decreasing

$\beta $ is not surprising when seen from the perspective of the thermodynamic limit, intuitive considerations could have also justified the opposite conclusion. Based on the fact that the relative amount of spontaneous emission coupled into the lasing mode increases with

$\beta $ (i.e., going towards the nanoscale) one could have argued for a nanolaser’s accrued sensitivity to noise and, with that, for a decrease in stability. The lsa belies this conclusion and advances potentially interesting consequences in applications, for instance in telecommunications [

121]. This consideration has the merit of clarifying an important difference: noise and intrinsic stability are not equivalent. We can therefore conclude that

nanolasers are intrinsically noisier devices, due to their low photon number, but they are structurally more dynamically robust than their micro- and macroscopic counterparts.#### 6.2. Eigenvector Analysis

The origin of spontaneous spiking (photon bursts) in the transition from a purely spontaneous emission regime to the coherent one can be better understood from the analysis of the eigenvector [

84].

Figure 7 shows the normalized component of the eigenvector which matches the photon number

n (i.e., the

$\nu $ component in Equation (

12)) for the least stable of the two eigenvalues (

Figure 4) for different values of

$\beta $ ranging from nano- to microlasers. As can be easily evinced from the curves, the other component (

$\mu $) has an amplitude close to 1 (except for

$\beta =1$ and

$\frac{P}{{P}_{thr}}>1$), given that the eigenvector is normalized. The second eigenvector (not shown), matching the more negative eigenvalue, is (also) almost entirely aligned with the population variable (

N) [

121], thus the two eigenvectors are nearly parallel below threshold. This observation has far-reaching implications, since it proves that, in spite of their origin (

$\nu $ or

$\mu $) fluctuations are strongly decoupled from the photon number variable: noise affects the population inversion and can make it deviate strongly from its equilibrium position, with only a weak transfer onto the photon field. This implies that the population inversion can grow quite far from its stationary value with a practically negligible consequence on the photon number. In other words, the transfer of a fluctuation onto the photon number is an unlikely process which takes place only occasionally, thus hindering the stimulated emission process.

The (near) absence of coupling to the photon field allows for the growth of population fluctuations of considerable amplitude. Thus, when a “transfer” towards the photon field takes place, the accumulated excess population enables the violent growth of the photon number through stimulated emission, and a photon burst ensues, in a way similar to a laser Q-switch [

50]. This is the origin of the observed pulses.

Physical considerations on the spontaneous emission can be used to give an alternative interpretation of this same picture. It is well-known that the spontaneous coupling of emitting dipoles to the radiation field is very weak, due to their extremely small size compared to the wavelength (

$O\left({10}^{-3}\right))$, which renders them extremely inefficient antennas. The well-known consequence is the weak population’s relaxation rate (also known as fluorescence rate) which, for the systems we consider here, is about five to six orders of magnitude smaller than the absorption (and stimulated emission [

1]) rate. This is the reason why below threshold—i.e., in a regime where stimulated emission is disfavoured—the population can grow to excessively large amounts. It is this built-up population excess which enables the transient ignition of stimulated emission (cf. bilinear term in Equations (

7) and (8)) responsible for the photon bursts.

Comparison of the photon component of the least stable eigenvector for different laser sizes is also illuminating. For the extreme nanolaser (

$\beta =1$), where all spontaneous emission is channeled into the lasing mode (black curve in

Figure 7) the relative weight of the photon eigendirection is stronger. This can be understood in terms of the transfer efficiency between population inversion and photon number, given that

**all** spontaneous emission ends up into the lasing mode. From a more mathematical point of view, this remark matches the fact that the coupling coefficient in Equations (

7) and (8) is larger thanks to

$\beta =1$. As

$\beta $ decreases there is a rapid reduction in coupling between the two variables, represented by the ever smaller component of the eigenvector projection along the photon axis. This holds true even close to the transition—identified by the end of the line traced in each figure—which corresponds to the transformation from a real into a complex eigenvector, or to the appearance of an imaginary part in the eigenvalue. In microlasers (e.g.,

$\beta ={10}^{-3}$) the amplitude of the photon number component of the eigenvector is only 1% even at threshold (

Figure 7).

The weak stability, which leads to the true phase transition in the thermodynamic limit, can be visualized by the transformation of a very large amount of stored energy into the field, giving rise to the (nearly) discontinuous transition. This is the picture that we normally use to describe threshold crossing by macroscopic lasers.

Figure 8 compares the photon component of the least stable eigenvector, as in

Figure 7, for a microlaser (

$\beta ={10}^{-4}$) and a macrolaser (

$\beta ={10}^{-6}$). The emerging feature is that the photon component remains zero until very close to threshold for the macroscopic laser (it is always very small for the microlaser) and that its value remains quite small (well below 1%) even at the transition. In addition, a rapid rotation of the eigenvector in its plane appears through a sign switch in the components for the microlaser (

Figure 9, green line—the diagonal lines join consecutive points). The macroscopic device shows a double rotation (dashed orange line), where the first signals the threshold crossing, and the second the end of the region where the eigenvalue is real (even above threshold). The fact that class B lasers [

117] possess a narrow region with real eigenvalues (i.e., no relaxation oscillations) above threshold has been known for a long time [

122], but can be observed only in devices where the relative time constants (

${\mathsf{\Gamma}}_{c}$ and

$\gamma $) are sufficiently close to permit it (i.e., not the CO

_{2} lasers considered in [

122]).

#### 6.3. Topological Conclusions

The combined information coming from the eigenvalues and the eigenvector analysis prove that the origin of the sudden photon spikes is dynamical and resides in the overall weak coupling between photons and population inversion below threshold—something which is to be expected, since the laser does not operate in that regime. The transition, however, is more or less sharp and sudden depending on the laser parameters: a nanolaser, more subject to intrinsic noise, is nonetheless less unstable and is expected therefore to display smaller and more frequent photon spikes in the so-called ASE region.

An interesting consequence of this analysis is that relaxation oscillations, believed not to exist in the thresholdless laser [

56] actually do occur even in these limiting devices, albeit only at large relative pump values

$\left(\frac{P}{{P}_{thr}}\gtrsim 2.1\right)$, as shown by the eigenvalues (

Figure 4). The physical argument advanced in favour of the lack of oscillations [

56]—absence of energy storing and immediate transformation of all the population inversion into photons—does not hold and is only a misinterpretation of the insufficient amount of coupling to initiate steady lasing at

$\frac{P}{{P}_{thr}}=1$ (

Figure 4). We also see, through the very small coupling at low pump rate (

$P\ll {P}_{thr}$) that the

zero-threshold connotation for

$\beta =1$ devices [

64] does not make any sense. This name originally came from the belief that the straight line relationship between input and output originated from the position of threshold at virtually zero pump. Although this is no longer considered to be true, the lsa gives clear evidence for a presence of threshold at

$\frac{P}{{P}_{thr}}>2.1$.

## 7. Symmetry Break between Spontaneous and Stimulated Processes

As already mentioned in

Section 5, differential models with Langevin noise do not predict the existence of photon bursts. Yet, the introduction of a phenomenological asymmetry in the relaxation constants, intended to model the aforementioned correlations in Quantum Dot devices, appears to show (small) bursts (Figure 5 in [

104]). Given the unusual nature of this result, it is worth trying to understand its origin, at least from a topological point of view. For this purpose, we transpose it onto the Quantum Well REs, such as Equations (

7) and (8) studied in

Section 6, so as to compare the changes introduced by the phenomenological assumption to the standard case we have examined in detail:

where

$\gamma $ now corresponds to the stimulated process and

${\gamma}_{s}$ to the spontaneous one. The symmetry is now (artificially) broken in the equation for the photon balance, since the stimulated and spontaneous rates are phenomenologically chosen to be different (the total rate in the equation for the population stays the same). We outline the changes, relative to

Section 6, by giving the equivalent expressions, which can be easily deduced from Equations (

16) and (17):

where

$\xi =\frac{{\gamma}_{s}}{\gamma}$ [

104] introduces the symmetry break (

$\xi >1$: superradiance;

$\xi <1$: subradiance). In the presence of the rate asymmetry the stability matrix becomes:

Concentrating on the case

$\beta ={10}^{-1}$, for concreteness, we look at the influence of

$\xi $ on the eigenvalues and on the “photon number component” of the least stable eigenvector, as in

Section 6.

Figure 10 shows the real part of the eigenvalue pair, which possess the same characteristic features as the standard REs, for different values of asymmetry

$\xi $. The red curve matches the corresponding one in

Figure 4 since

$\xi =1$ corresponds to the symmetric rate configuration studied in

Section 6. Superradiance (

$\xi >1$) stabilizes somewhat the nanolaser by rendering more negative the eigenvalue closer to 0 (green curve), while subradiance (

$\xi <1$) renders it a little more susceptible to noise by bringing it closer to 0 (black curve).

The eigenvector analysis confirms the destabilizing quality of subradiance (

Figure 11) as opposed to the superradiance-induced stabilization. The amplitude of the normalized photon component of the least stable eigenvector (matching the upper eigenvalue in

Figure 10) becomes progressively smaller as the symmetry break proceeds from super- to subradiant interactions (colour coding as in

Figure 10). As already discussed in

Section 6 this amounts to saying that the component of this eigenvector aligned with the photon number variable

n is small (extremely small in the subradiant case) over a broad range of pump values. Hence, this eigenvector is nearly aligned with the population variable

N (the other eigenvector—not shown—is also aligned with

N). Thus, the fluctuations which occur in this eigendirection—due to the closeness to 0 of the eigenvalue (

Figure 10)—influence almost exclusively the population

N, which may undergo large excursions, since the only relaxation channel is the spontaneous relaxation (partially inhibited by subradiance). As discussed at length in

Section 6, if the weak coupling to the photon variable takes place in the presence of a large excess in

N, then a photon burst takes place; more readily in the case of subradiant interaction.

The analysis shows that the symmetry break quantitatively affects the phase space structure and the photon burst appearance; no qualitative—e.g., structural—changes are introduced into the modified phase space. For the standard REs (for Quantum Wells) this would amount to larger photon bursts (more frequent depending on the parameter details). However, since photon bursts and superthermal statistics were already predicted under these conditions [

72,

84,

114,

115,

116,

119], the symmetry break does not seem to bring anything new. On the other hand, it is possible that for a sufficiently large asymmetry the integration of the REs model may also produce bursts, thanks to the laser’s closeness to the instability threshold which may be reached even with perturbations of Langevin type [

123].

For completeness, it is important to remark that in Equations (

16) and (17) we have attributed the asymmetry only to the decay into the lasing mode (no specific indication is given in [

104] as to the choice made). The extension of the asymmetry (in particular the subradiant case,

$\xi <1$ for our purposes) to the whole relaxation from the excited state (i.e., also to the last term of Equation (17)), would have lowered even more strongly the stability of the device, both in the eigenvalue and in the eigenvector components (results not shown). This way, the small amount of Langevin noise—proportional to the photon number (very small below threshold), could even more easily produce photon bursts.

## 8. Conclusions

The objective of this contribution has been the presentation of a summary of the current observations of photon bursts and their comparison with the properties of ASE. No conclusion is drawn about the true nature of these bursts, whence the quotes around Amplified Spontaneous Emission in the title, which undoubtedly require further evidence and investigation.

Photon bursts are seen to precede the establishment of coherence in micro- and nanolaser and therefore take a significance which is not necessarily present in their macroscopic counterparts. Even though their role is not fully characterised, they appear as an instance of a nascent instability, as shown through the analysis of the phase space properties of the most basic physical model. An interesting question, not yet fully answered, is the inability of differential models of predicting the appearance of photon bursts—up to the recent exception [

104] which probably benefits from the destabilization found in

Section 7—while the intrinsic granularity and randomness arising from numerical stochastic approaches give an excellent match to the observations, independently of the model details. The different aspects of the stochastic predictions are summarized and compared to the main conclusions which can be drawn from the current experimental knowledge of photon spikes. The experimentally relevant question of the identification of photon bursts through the second-order autocorrelation—the most widespread detection tool—is addressed with the help of a simple model, which explains the role of the various parameters (burst duration, relative height of background to peak, filling factor in the observation time interval) on the actual identification of a superthermal statistics.

Devices with very few QDs are now consistently built, with good control of their numbers and properties [

124] and limitations exist on the minimum number of emitters necessary for lasing [

111,

112]. Interesting questions arise in the dynamical interaction among photons and emitters in the below-threshold region, given the very restricted number of QDs actually interacting (as

matter) with the thermal or pseudo-thermal radiation inside the cavity. In other words, it would be interesting to understand whether the thermodynamic limits taken in macroscopic out-of-equilibrium systems hold [

44]. This question is certainly one that will need to be addressed in a near future in order to shed further light on the physics of radiation in interaction with matter.

On the list of forthcoming work on this topic the following items ought to appear: a better characterization of photon bursts, including their coherence time and their statistical properties (such as repetition rate, amplitude statistics to probe for possible systematic or deterministic aspects); the evolution of the photon bursts into the cw (noisy) coherent emission; the possible relationship between burst features and cavity parameters, currently hinted to by stochastic modelling; the possible role that these bursts may play as a connection between the radiation statistics in thermodynamic equilibrium below and above the transparency threshold and the coherent, lasing regime.