Research ArticleOPTICS

Quantum coherence–driven self-organized criticality and nonequilibrium light localization

See allHide authors and affiliations

Science Advances  16 Mar 2018:
Vol. 4, no. 3, eaaq0465
DOI: 10.1126/sciadv.aaq0465


Self-organized criticality emerges in dynamical complex systems driven out of equilibrium and characterizes a wide range of classical phenomena in physics, geology, and biology. We report on a quantum coherence–controlled self-organized critical transition observed in the light localization behavior of a coherence-driven nanophotonic configuration. Our system is composed of a gain-enhanced plasmonic heterostructure controlled by a coherent drive, in which photons close to the stopped-light regime interact in the presence of the active nonlinearities, eventually synchronizing their dynamics. In this system, on the basis of analytical and corroborating full-wave Maxwell-Bloch computations, we observe quantum coherence–controlled self-organized criticality in the emergence of light localization arising from the synchronization of the photons. It is associated with two first-order phase transitions: one pertaining to the synchronization of the dynamics of the photons and the second pertaining to an inversionless lasing transition by the coherent drive. The so-attained light localization, which is robust to dissipation, fluctuations, and many-body interactions, exhibits scale-invariant power laws and absence of finely tuned control parameters. We also found that, in this nonequilibrium dynamical system, the effective critical “temperature” of the system drops to zero, whereupon one enters the quantum self-organized critical regime.


The self-organization of many nonequilibrium complex systems (1) toward an “ordered” state is a profound concept in basic science, ranging from biochemistry to physics (24). Examples include the group movement of flocks of birds (5), motions of human crowds (6), neutrino oscillations in the early universe (7), and the formation of shapes (“morphogenesis”) in biological organisms (8, 9). An intriguing trait of this nonequilibrium, driven-dissipative systems (2, 3) is that their self-organization can lead them to a phase transition and to critical behavior—a phenomenon known as self-organized criticality (SOC) (10). Unlike equilibrium phase-transition phenomena, such as superconductivity or ferromagnetism, where an exogenous control parameter (for example, temperature or pressure) needs to be precisely tuned for the phase transition to occur, no such fine-tuning is needed in SOC systems (1013): They can self-organize and reach their critical state even when driven far away from it. Similar to all critical phenomena (14), at exactly the phase-transition point, the SOC systems are described by scale-invariant power laws with well-defined critical exponents (15). Biological extinction events, dynamics of granular or mechanical media, and nonlinear processes in astrophysics, in magnetospheric, planetary, solar, and stellar physics, and in cosmology are examples of dynamical systems exhibiting SOC behavior (1013, 16).

As is evident from the above examples, most SOC systems studied until now have been classical, driven by a control parameter equivalent to temperature in conventional phase transitions (14) (“thermally” driven). It is thus intriguing to identify quantum-coherently controlled systems that exhibit a phase transition in a self-organized way, wherein the critical point behaves as an attractor of the dynamics. In such a nonequilibrium, self-organized critical system, it may also be possible to enter the regime of quantum SOC (1113, 1720)—similar to the regime of quantum criticality in equilibrium condensed matter phase transitions (2125). For a complex system to enter such a regime, several requirements should be fulfilled (1017, 2123). These include (i) the dropping of the control parameter (effective or equivalent “temperature”) to zero upon application of an external coherent field (2123), (ii) the absence of finely tuned control parameters (10, 11), (iii) automated attraction toward a nonequilibrium steady state (15, 16), (iv) emergence of scaling and power laws with well-defined critical exponents (1012, 15), (v) statistical independence of the events (“avalanches”; compare section S4) (1013), and (vi) a nonlinear coherent growth phase upon exceeding a critical threshold (1013).

To address the above issues, and profiting from the principle of universality in statistical physics, we introduce and study here a semiclassical many-body system in active nanophotonics that exhibits all characteristics of classical SOC but can be quantum-coherently driven. Specifically, we investigate whether light can exhibit a quantum self-organized phase transition to localization in fluctuating dissipative (2, 3) media where simple cavity effects are absent. Many-body interactions (26, 27) in the driven-dissipative nanostructure considered (Fig. 1) arise from the competition between many photonic modes in the presence of nonlinear gain (Fig. 2 and Materials and Methods) and, as shown later on, are vital for robust localization to emerge. Likewise, quantum effects in this system arise from the presence of a coherent drive Ωa mixing two quantum states in the active core layer of the system (Fig. 1 and section S3), leading to their quantum-coherent superposition—similar to quantum critical systems (for example, ferromagnets) in condensed matter where the application of a classical magnetic field realigns the spins in the direction of the field, leading to a quantum-coherent superposition of “up” and “down” spins (for example, |i=(|i+|i)/2), that is, to a “quantum paramagnet” (2125).

Fig. 1 The conceived coherence-driven photonic nanosystem.

The nonequilibrium system we consider is a quantum-coherently driven active nanophotonic heterostructure in its stopped-light regime. It consists of a 26-nm-thick quantum well (QW) bounded by aluminum (Al), with the upper cladding having a thickness of 50 nm. The QW is described by a three-level scheme (57), as shown in the lower inset. An incoherent incident pump (shown in red) sends electrons from level |1> to level |3> with a pump rate g, whereas a coherent drive (shown in white) couples the levels |2> and |3> with a Rabi frequency Ωa. Photons emitted in the |2>-|1> transition stimulate coherently the stopped-light (“heavy photon”) mode that the heterostructure supports (see Fig. 4). The pump beam deposits gain over the middle region highlighted translucently. Numerical Maxwell-Bloch-Langevin computations (61, 62) reveal the emergence of a sinc-shaped nonequilibrium-localized light field (shown in false color), in agreement with theory (see Fig. 2B, right).

Fig. 2 Synchronization in classical mechanics and in optics.

(A) The dynamics and emergence of synchronization in many-body classical systems can be modeled using coupled mechanical oscillators (metronomes) (16). When a collection of metronomes is placed on a fixed basis (top), in which case they do not couple, they oscillate randomly (“disordered phase”), never synchronizing. In contrast, when the metronomes are placed on a basis that can freely move (bottom), in which case the metronomes couple, they eventually synchronize without any external intervention (that is, “spontaneously”), reaching an “ordered” state where they all move in phase. (B) In the active nanophotonic heterostructure that we study (Fig. 1), a large number of oscillating lasing modes couple nonlinearly in the gain region (translucent in Fig. 1) just as the metronomes of (A) couple nonlinearly via the moving base. Each lasing mode can be considered a “particle” in a many-body nonequilibrium system (27) and, on its own, is delocalized, oscillating in space with a spatial frequency βn (n = 1, 2,…) (left). These {βn} modes initially oscillate out of phase, giving rise to a noise pattern resulting from their superposition and interference (left). However, owing to their nonlinear coupling, they eventually synchronize in space after a transient time interval, giving rise to the emergence of a localized sinc function–shaped supermode (right). The onset of spatial synchronization of the {βn} modes is associated with a nonequilibrium first-order phase transition (see main text and Materials and Methods).

Localization in optics and photonics is usually achieved by using a barrier (optical potential) to confine light, for example, in microstructured dielectric cavities or plasmonic nanoparticles and metamaterials. In all cases, this type of localization is thus known as “potential localization” (2830): The light wave is localized, owing to the sheer presence of a refractive index barrier, without any phase transition occurring. Most existing light localization active structures (such as, for example, optical microcavity or photonic crystal lasers, light-emitting diodes, and modulators) normally involve intricate optogeometric configurations, are difficult to integrate with waveguides and further multiplexing, and are not CMOS (complementary metal-oxide semiconductor)–compatible nor suitable for scalable, high-volume fabrication (for example, often requiring the creation of air bridge structures or fabrication using electron beam lithography, which is an accurate but time-consuming process) (29). Thus, it is desirable from a technological perspective to design light-localizing structures that can, ideally, be in the form of planar heterostructures, fully compatible with integrated CMOS platforms (for example, suitable for stepper photolithography), and minimalistic—not requiring the fabrication of well-defined longitudinal mirrors/boundaries to confine light. Furthermore, until now, the only known nonpotential light localization mechanism is the Anderson transition in strongly disordered media (2832) where in three dimensions an equilibrium phase transition occurs from diffusion to exponential localization. However, the so-attained nonpotential light localization is not robust: Dissipative losses (33, 34), nonlinear photon interactions, or time-varying refractive indices (32) can all destroy localization and restore diffusion (see also section S1) (28, 30, 32).

Here, we report on a coherence-driven photonic nanosystem that can dynamically localize light far from thermodynamic equilibrium (35) by interfering many nonlinearly coupled photonic modes in the presence of gain. In sharp contrast to conventional light localization schemes, this is achieved in the absence of refractive index barriers (non–potential-well localization); that is, there is a phase transition to localization—similarly to the Anderson scheme. However, although there is a phase transition to localization in both our case and the Anderson configuration, our scheme is inherently different: It does not involve disordered media but longitudinally uniform, active plasmonic heterostructures; thus, the localization mechanism here fundamentally differs from that in disordered media (where coherent backscattering by disorder leads to destructive interference of waves returning to points in the medium along time-reversed paths). Furthermore, unlike the Anderson scheme, the attained localization in our case is robust to the simultaneous presence of dissipation, many-body nonlinear interactions, and time-varying media parameters. In this quantum-coherently driven, nonequilibrium configuration, we observe all six of the previously outlined characteristics that a complex system must have to enter the regime of quantum SOC. Moreover, we identify a supercritical region in the phase diagram that is inaccessible to all classically driven SOC systems, emerging only if quantum effects drive the SOC nanosystem.

We note that a variety of light-based many-body systems (36), as well as of open, driven-dissipative quantum systems (37), have been reported in the past, such as cold atomic ensembles and quantum gases in optical cavities (38), arrays of nonlinear cavities and superconducting circuits (enabling quantum simulators) (39, 40), and exciton-polaritons in semiconductor microcavities (41, 42), to name only a few. These systems exhibit a rich variety of dissipative, nonequilibrium dynamics [for example, described by a dissipative Gross-Pitaevskii equation (41)] or even quantum phase transitions such as, for example, open-system realizations of the Dickie quantum phase transition (38). However, these systems typically do not exhibit classical or quantum SOC because they do not fulfill one or more of the aforementioned six general characteristics (1013) of this type of phase transitions. In particular, those systems typically require precise tuning of a range of parameters, such as, for example, the pumping strength of a far-detuned laser beam, atomic detuning, cavity detuning, atom-field coupling strength, etc., without which dynamical instabilities are prevented. Even with the right choice of parameters, those systems typically exhibit self-organization (rather than SOC), similar to Bénard cell convection or to pattern formation far from equilibrium (2, 3, 8, 9), such as polarization gratings induced by intense femtosecond light irradiation of amorphous materials (42), involving system-wide processes that operate with long-range interactions, rather than statistically independent events (avalanches) arising from next-neighbor interactions (1113). For the above reasons, those systems are not classified as SOC in the pertinent literature. In contrast, as shown here in the following, our quantum-coherently driven light-localizing nanophotonic system inherently exhibits all generic characteristics of SOC systems, and for a judicious choice of parameters, it may also be possible for the system to enter the previously described regime of quantum SOC.


The quantum-coherently driven complex nanosystem

To see how robust nonpotential light localization and quantum coherence–driven SOC can arise, consider an active nanophotonic heterostructure, as shown in Fig. 1. It consists of a deep–ultraviolet (UV) QW [for example, AlGaN/AlN (43, 44)] bounded on both sides by a nanoplasmonic medium, Al (45, 46). The active nonlinearities in the QW allow for many-body interactions among the supported photonic states (Fig. 2) (27), whereas the nanoplasmonic layers (47) provide, in the stopped-light regime, the positive-feedback mechanism necessary for the system to reach criticality in a self-organized manner without the presence of an actual optical potential in the longitudinal direction (see also section S1).

In our analysis of this complex system, we first need to determine an accurate expression for the macroscopic permittivity of the Al layers, εAl(ω), and for this, we first focus on the electron-electron interactions (26) in the cover and substrate regions. In Materials and Methods, we present the details of the derived many-electron description of the polarization process in the above layers in the deep-UV regime. We find excellent agreement between experimentally extracted data for the permittivity of Al, εAl(ω), from Ehrenreich et al. (46) and a fit using our analytically derived expression, verifying the accuracy of this description of the deep-UV plasmonic medium. Note from Fig. 3 that the standard Drude theory fails to predict Re{εAl(ω)} < 0 for ω > ~13 × 1015 rad/s and, in particular, in the deep-UV region, Re{ω} ~ 16.8 × 1015 rad/s, which is the region of interest to us here.

Fig. 3 Many-electron characteristics of the nanoplasmonic layers of the structure of Fig. 1.

(A) Feynman diagrams illustrating schematically the polarization process (top) and renormalization of the electron-electron interaction strength, Ve, to Veff (bottom) in the considered plasmonic medium (Al). 1PI, one-particle irreducible. (B) Comparison between the real part of the macroscopic permittivity of Al (−Re{εAl}), as derived from the microscopic theory presented here (Eq. M1), the experimental data by Ehrenreich et al. (46), and the standard Drude formula (45).

Next, we calculate the complex frequency (48, 49) (complex ω) photonic states that the heterostructure supports (see Materials and Methods). Figure 4 presents an example of such a complex band structure calculation for the TM2 photonic mode, showing the real and imaginary parts of the angular frequency ω as a function of the (real) longitudinal propagation constant β. Note from Fig. 4 (bottom) that there are two points, one at β = 0 and a second at β ~ 30 × 106 m−1, where the group velocity vg → 0, making the band very flat for β < 40 × 106 m−1 and inducing a large effective mass of the photonic states in that region (“heavy photons”).

Fig. 4 Complex band diagram of the nanophotonic heterostructure.

Complex frequency dispersion band of the TM2 mode. Also shown schematically in the upper panel with a filled Lorentzian shape is the bandwidth of the QW gain (57), whereas highlighted are the excited {βn} states. The lower panel shows the group velocity vg of the TM2 mode, which becomes zero at two points in the region where β < 40 × 106 m−1. The inset in the lower panel presents an intuitive ray picture understanding of how positive feedback arises in the considered self-organized complex nanosystem.

We now pump a small region of the active nanoplasmonic heterostructure with a Gaussian incident beam (white beam in Fig. 1), depositing sufficient gain to exceed the dissipative losses (given by Im{ω} in Fig. 4) throughout the {β ∈ R: β < 40 × 106 m−1} TM2 band, indicated with a dashed colored circle in the dispersion diagram of Fig. 4. Each βn point in that part of the band may now become a lasing state (an optical oscillator) (50, 51), with all of them coupling nonlinearly in the gain region where they overlap spatially; that is, we have a many-body system made of a continuous ensemble of coupled optical nano-oscillators, similar to the prototypical SOC model of coupled mechanical oscillators (see Fig. 2) (16). The optical oscillators are not identical here because, although they have nearly the same oscillating frequency ω (because the band is flat in the above {β} region), they have slightly different losses (red dashed curve in Fig. 4). Nonetheless, we find below, in Fig. 5A, that all {β} states become unstable and spatially synchronize, similar to the synchronization in time of the metronomes of Fig. 2A (bottom).

Fig. 5 Spontaneous synchronization and emergence of attractors.

(A) Onset of synchronization: temporal dynamics of two coupled {βn} states showing repeated cycles of asynchronous (AS) motion, followed by synchronization (S), and then temporal coherence buildup (CB). In all cases, the onset of spontaneous synchronization precedes that of activity (CB). (B) Emergence of attractors: dynamic evolution in the phase space (control parameter/inversion versus order parameter/photon number) of the localized light field of Fig. 1. After a transient period (0 < τ < 3), during which the light field remains robustly localized and the permittivity of the pumped core region dynamically changes (the control parameter ρ22 − ρ11 ∝ Im{εcopumped}), the system is attracted to a stable point, robust to perturbations. The inset shows the same result but without quantum coherence driving the system.

To clearly see the emergence of positive feedback in our structure, without which the system cannot reach criticality in a self-organized way, we visualize in Fig. 4 (bottom) the trajectory of a guided light ray in that regime by calculating analytically (section S2) the total time-averaged power flow in the structure. This gives the following compact expression: Ptot = (1/4)ExmaxHymaxweff, where Exmax (Hymax) is the maximum value of the Ex field (Hy field) component in the core layer, and weff is the effective thickness of the heterostructure. In the region where the TM2 band in Fig. 4 is flat and vg approaches zero, Ptot must also approach zero (52, 53) so that weff → 0; that is, a light ray becomes trapped, forming a double light cone, as illustrated in the lower inset of Fig. 4. This feedback-forming mechanism is crucial for filtering out of out-of-phase fluctuations and the emergence of robust attractors, as we see below.

Semiclassical mean-field model

We ascertain the regime of quantum coherence–assisted SOC on the basis of a semiclassical mean-field model. We apply a coherent driving field Ωa (red beam in Fig. 1) to the transition |2> → |3> of the three-level active core medium (QW) (5458), as shown in the inset of Fig. 1. This creates an asymmetry between absorption and stimulated emission on the |2> → |1> transition and maintains the quantum coherence on the |2> → |3> transition. In contrast, the role of the incoherent pump, g, is simply to excite electrons to the upper state |3>, which then decay to states |1> and |2>. We design the nanostructure such that the |2> → |1> transition is resonant with the {Re{ω} ≈ 16.8 × 1015 rad/s; β ∈ R: β < 40 × 106 m−1} TM2 band. The Heisenberg equation of motion for the slowly varying amplitude α0βn of an individual βn state is given by (see section S3)α.0βn=Γnα0βn+ipρ21(p)Ωb(p)(1)where |α0βn|2 is the number of photons in the mode βn. The decay of excitation of the mode βn is quantified by the parameter Γn, whereas the polarization of the active (gain) medium is given by the nondiagonal matrix elements, such as (section S3)ρ.21=Γ21ρ21iΩb(ρ22ρ21)+iΩa*ρ31(2)

A key ingredient of Eq. 2 is the presence of the matrix element ρ31, which indicates the coherence between levels |3> and |1> that develops as a result of the coherent drive field Ωa. This coherence between the atomic states gives rise to quantum interference that is destructive for absorption and constructive for emission processes (55, 56, 58). The electric field associated with α0βnis given byE(z)=in=0N1{2πωn/V)1/2ei[(β1+nδβ)z+φn]αnt=0c.c.}(3)where we note that the phases φn are different for each mode. Eventually, because of the nonlinear coupling of the βn modes in the pumped region (Fig. 1), the modes synchronize as is shown in Fig. 5A, that is, they all obtain the same phase φn, which can be set to zero without loss of generality. The supermode arising from the synchronization of all of these modes has a sinc-shaped spatial profile given by (see Materials and Methods)I(2/Lloc)2sin c2(L/Lloc)(4)where Lloc = 2/Δβ is the localization length. Finally, as shown in Materials and Methods for the present system, but as it is also known for active (lasing) systems in general (50, 51, 55), the control parameter (“effective temperature”) in this configuration is the amount of gain in the pumped region, that is, the population inversion ρ33+ρ22ρ11.

From the above discussion, it is thus evident that the succession of elements used in the present model to attain quantum coherence–driven SOC and nonequilibrium light localization is as follows: (i) A pump field deposits gain (population inversion ρ33+ρ22ρ11) in a finite region in the core of the heterostructure (Fig. 1); (ii) the modes βn destabilize and start oscillating, each with an amplitude α0βn(t) (Eq. 1); (iii) the modes initially oscillate out of phase (Fig. 2B, left) but eventually synchronize (Figs. 2B, right, and 5A), owing to their nonlinear coupling in the gain region; (iv) the synchronization of the modes gives rise to a sinc-shaped localized hotspot (Fig. 1 and Eq. 4); (v) the light intensity in this hotspot is controlled by the population inversion ρ33+ρ22ρ11, which is a function of the coherent field Ωa (section S3 and Eq. 2); (vi) for a specific value of Ωa, the critical population inversion ρ33c+ρ22cρ11c [“critical temperature” (50, 51, 55) of the system] drops to zero (Fig. 6B), which is a requirement for entering the regime of quantum SOC.

Fig. 6 Emergence of scaling and toward the quantum self-organized critical regime.

(A) Order parameter (number of photons) versus the pump rate, g, with and without quantum coherence driving the nanosystem. The nonlinearly interacting {βn} lasing modes self-organize into a critical state, with a second-order (first-order) phase transition to a supercritical phase being observed without (with) the quantum-coherent drive. (B) Emergence of scale-invariant (fractal) power law in the order parameter versus the normalized coherent drive Ωa=Ωa/(ρ33c+ρ22cρ11c), where ρ33c+ρ22cρ11c is the critical inversion. The computed critical isotherm exponent is δ ~ 3.07. The upper inset shows the variation of the critical inversion with the applied coherent field Ωa. For another suitably selected choice of parameters (see section S3), the control parameter (that is, the critical inversion ρ33c+ρ22cρ11c, which plays the role of the effective temperature here) reduces to zero for a sufficiently strong drive Ωa (lower inset). That point is the quantum self-organized critical point (QSCP) of our nonequilibrium dynamical system, somewhat analogous to the quantum critical point (QCP) in equilibrium static systems in condensed matter (2125).

Onset of activity

There is a further crucial point discerned from Fig. 5A: The onset of synchronization is followed in time by the onset of “activity” (59), that is, here, by buildup of temporal coherence (CB), which then destroys synchronization, leading to AS motion (AS2 region in Fig. 5A), followed again by synchronization (S2) and then CB (CB2). We thus have repeated cycles of “AS” motion→leading to “S” motion→leading to “CB”→leading to “AS” motion, until the nonequilibrium steady state is reached. Note that, even at a steady state, this semiclassical many-body (27) nanosystem is still far from thermodynamic equilibrium because the inverted active medium is described by an occupation Boltzmann distribution function with a negative effective temperature (Tactive-medium < 0), whereas the nanostructure has a temperature much smaller than ω/kB (kB being Boltzmann’s constant); hence, we may effectively put Tnanostructure ≈ 0. Thus, the light field is coupled to two reservoirs with very different temperatures, and for this reason, the complex nanosystem is far from thermal equilibrium. Note also that, similar to the sandpile cellular automaton SOC model (10, 15), where although the current of incoming particles is known, the control parameter (slope increase of the sandpile) varies randomly and does not require external fine-tuning, also in our case, although the pump rate may be known, the control parameter [atomic inversion (v)] cannot be precisely known or tuned, owing to spontaneous emission; in both cases, the complex system is automatically attracted to the nonequilibrium self-organized critical point (Fig. 5B and Materials and Methods).

The evolution in the phase space of the robustly (and nonpotentially) localized hotspot, which, as we saw, arises from the superposition of all {β} states, is reported in Fig. 5B, showing the emergence of a steady-state nonequilibrium attractor, robust with respect to noise; that is, the self-organized critical point becomes an attractor for the dynamics. Note also from Fig. 5B the marked evolutionary advantage of the quantum coherence–driven complex nanosystem compared with the case where the system is only classically driven (inset in Fig. 5B): First, steady state is reached faster, and second, the number of photons generated, that is, the intensity of the nonequilibrium-localized hotspot, increases by more than an order of magnitude.


To better observe the effect of the coherent control on the considered nonequilibrium nanosystem, we have derived and show in Fig. 6A its phase diagram, both with and without the external coherent drive. Note how the continuous, second-order phase transition in the absence of quantum coherence is transformed to a discontinuous first-order phase transition when the coherent field drives the SOC nanosystem. We find from Fig. 6A that the presence of quantum coherence not only preserves but also actually enhances SOC: The magnitude of the order parameter markedly increases under the quantum-coherent drive, and furthermore, there is a region (highlighted with pale blue in Fig. 6A) that was subcritical in the absence of quantum coherence but becomes supercritical in the presence of quantum coherence.

Finally, Fig. 6B presents two of the key results of the present study: First, similar to classical SOC systems (15), we find (in mean-field approximation) that the order parameter obeys a scale-invariant power law as a function of the normalized quantum-coherent drive Ωa with a “critical isotherm” (14, 50, 51, 55) exponent δ ~ 3.07. Our quasi-1D system belongs to the Ising universality class, but unlike equilibrium static systems in condensed matter, here, the critical point is reached by starting far from thermodynamic equilibrium in a self-organized manner. The presence of an attractor in Fig. 5B is insensitive to the parameters of the model. This robustness to dissipation, fluctuations, and many-body interactions is the fundamental difference of the present type of nonpotential light localization compared to light localization in disordered media. Second, similar to quantum criticality in equilibrium static systems, which is reached when the critical temperature Tc drops to zero upon application of an external magnetic field (2123), the lower inset of Fig. 6B shows that, for a suitable choice of parameters (section S3), the effective critical temperature of our nonequilibrium dynamical system falls to zero for sufficiently strong coherent drive (Ωa ~ 0.6). Exactly that point, where the total critical inversion Δρc (critical effective temperature; see Materials and Methods) falls to zero, is a quantum self-organized critical point, analogous to the quantum critical point in equilibrium phase-transition phenomena (compare Fig. 1) (2124) but without the need for fine-tuning.

In conclusion, the main findings of the present study can be summarized as follows: First, we have identified a quantum coherence–driven self-organized phase transition in an active nanophotonic system (Figs. 1 and 4B). Second, we have shown that quantum coherence bestows a marked advantage in the evolution of self-organized systems toward a phase transition; namely, it accelerates the reaching of the critical point (Fig. 5B) and enhances the order parameter (Fig. 6A). Third, we have identified a scheme where unconventional light localization is attained: Similar to the Anderson scheme in disordered media, there is phase transition to localization (that is, the localization is not trivial because it is not aided by cavity-like effects); however, unlike the Anderson scheme, we do not deploy disordered media (but active, longitudinally uniform, plasmonic heterostructures), and the attained dynamical localization remains unharmed by and robust to dissipative effects, nonlinear many-body interactions, and time-varying media parameters (Figs. 2B and 5). The fact that there is no need to design optical potentials in this scheme for robustly localizing light [that is, the fact that the conceived nanostructure is minimalistic (planar and longitudinally uniform)] strongly suggests that it is also appealing from an applied perspective for a host of applications in condensed matter photonics and nanophotonics.


Many-electron polarization process in the cover and substrate layers

The response of the electron gas to a spatiotemporally varying weak potential Ve(r, t) = Ve(q, ω)ei(qr−ωt)eγt (γ → 0) is governed by the Hamiltonian (26) H=p,sεpa^ps+a^ps+sΦ^s+(r)Ve(r,t)Φ^s(r)d3r, where εp = 2k2/2m is the energy of a single-electron state, a^ps(+) is the annihilation (creation) operator for a state with wave vector p and spin s, and Φ^s(+)(r)=1V1/2pa^ps(+)e()ip.r is the corresponding electron field operators. In the Heisenberg representation, the equation of motion of the relevant density operator is, thus, given by id(ρ^p,q,s)/dt=(εp+qεq)ρ^p,q,s+(a^p,s+a^p,sa^p+q,s+a^p+q,s)Ve(q,ω) from where upon taking the thermal average (<Q^>= Tr(Q^eβH)/Z, β being the thermodynamic beta, and Z being the canonical partition function), we found the dynamical response function χ0(q,ω)=(1/V)p,s(np+q,s0np,s0)/(εp+qεqωiγ), where np,s0=<a^p,s+a^p,s>. For the polarization process in the plasmonic medium, we needed to sum all one-particle irreducible Feynman diagrams that can be inserted into an interaction wiggle [an example of which (for iχ0) is shown in Fig. 3A]. In the absence of a photon propagator at this stage, the aforementioned sum will simply renormalize the potential, changing the interaction strength from −iVe to the effective potential −iVeff. Carrying out the sum to infinity in the usual way (see Fig. 3A), we got −iVeff = −iVe + (−iVe)(iχ0)(−iVe) + (−iVe)(iχ0)(−iVe) (iχ0)(−iVe) + … = (−iVe)/(1 − Veχ0), from where the equivalent permittivity may be defined as ε(q, ω) = 1 − Veχ0 = 1 − 4πe2χ0(q, ω)/q2. Here, we sought to identify the smallest possible nanovolume over which nonequilibrium light localization of oscillatory (that is, not surface) states may be achieved in the considered Al nanoplasmonic heterostructure; hence, we needed to consider the highest possible incident light frequencies (in this case, in the deep-UV regime), for which not only conduction band (that is, sp band) but also d band to sp band electronic transitions are important. In particular, Al (and other non-noble metals) contains multiple valence electrons that, for the high-energy deep-UV incident photons that we shall consider here, can transition between parallel or near-parallel one-electron conduction bands (sp band to sp band). Thus, assuming that |q| ≪ kF (kF being the radius of the Fermi sphere in k space around k = 0), expanding χ0(q, ω) in q, taking into account that εp+q = εp + q·pεp+…, np+q0=np0+(n0/ε)q·pεp+, and n0/εp=δ(εpεF) (at T = 0), withpεp=υFk/k being the Fermi velocity, and using standard second-order perturbation theory, we arrived at the following generic expression for the nonnoble metal (Al) dynamic permittivityεAl(ω)=lim|q|0ε(q,ω)4πie2ωm2i,jf(Ei)f(Ej)EiEj<i|pa|j><j|pb|i>iω+1/τ+(i/)(EiEj)(M1)where the sum is over the valence and conduction band states i and j, f(Ei) and f(Ej) are the Fermi functions (appearing owing to the Pauli principle), “a” is the direction upon which the current density flows upon application of an electric field in the direction “b,” τ is the relaxation time, and < i | pa | j > is the magnitude of the momentum matrix element coupling the valence band state i and the conduction band state j. Figure 3B reports an excellent agreement between experimentally extracted data for εAl(ω) by Ehrenreich et al. (46) and a fit using Eq. M1 above, verifying the accuracy of the above many-electron theory.

Transfer matrix and argument principle methods

Assuming that a light wave in the i-th layer oscillates as φi1cos[κi(xxi−1)] + φi2(uii)sin[κi(xxi−1)], with u = εii) for transverse magnetic (transverse electric) modes and that the light wave in the terminating Nth layer is of the form αNeγN(xxN1), the photonic states obey a general dispersion relation of the form (48) F(ω, β2) = φN2 − φN1N/uN) = 0, where γN = (β2k02nN2)1/2, β being the longitudinal propagation constant, k0 being the vacuum wave number, and φN1 and φN2 being the oscillatory wave constants. The M roots, w, of the (analytic) function F in a bound domain K with smooth boundary ∂K can be found using the argument principle method (APM), according to which 12πiKwF(w)F(w)dw=i=1Mwi, with being a positive integer. The number of roots in the domain K is identified by setting = 0. Upon isolation of a root, we may set = 1 to locate its exact position on K. For the applicability of the APM, the complex function F(ω, β2) must be analytic and without poles in the domain K; however, in its above form, F(ω, β2) has two branch points at β2 = k02n12 and β2 = k02nN2, resulting in a four-sheeted Riemann surface [where, for a given value of β2, F(ω, β2) takes on four possible values]. To unfold the four-sheeted Riemann surface (β2 plane) and map it onto a new complex variable z (H(z) = F[ξ(z)], with ξ representing the original complex solution plane that is mapped onto the new complex variable z), we introduced two functions, g(z) and h(z), demanding that g2(z) = β2 − ξ(z)n12 and h2(z) = β2 − ξ(z)nN2, which can also be written as g(z) = n1(ez + Rez) and h(z) = nN(ezRez), where R = β2(nN2n12)/(4nN2n12); hence, ξ(z) = Se2zR2e−2z, with S = β2(nN2 + n12)/(2nN2n12), from where the angular frequency may be calculated as ω(z) = c[ξ(z)]1/2.

First-order phase transition to synchronization of the {βn} states

The synchronization in space (rather than only in time) of the {βn} states shown in Figs. 2B and 5A, leading to robust nonpotential localization, is a first-order phase-transition phenomenon arising from their nonlinear coupling. To better observe this, we started by considering the general equation governing the coupling between the field amplitudes αn in the presence of dispersion and nonlinearity (27, 60) α.n=(χnξniδn)αn+(υsat+iυnonl)pp+q=nαpαp*αq+SFn, where χn is the gain at the frequency of the βn state, δn is the phase accumulated over a cycle (associated with group-velocity dispersion), υsat is associated with the saturable amplitude, and υnonl is associated with the third-order nonlinearity. The frequency-dependent gain saturation may, as usual, be modeled by χn=χ0,n/(1+nχ0,nαnαn*/nχ0,nsatαnsatαnsat*), where the constant nχ0,nsatαnsatαnsat* term represents the saturation power. It is now a simple exercise to show that the field amplitudes αn obey, equivalently, the equation α.n=iH/αn*H/αn*+SFn, where H=υnonlpp+q=nαpαp*αqαn*nδnαnαn* and H=υsatpp+q=nαpαp*αqαn*+nξnαnαn* nχ0,nsatαnsatαnsat*ln(nχ0,nsatαnsatαnsat*+nχ0,nαnαn*). From Fig. 4, we can see that the {β} band is very flat (that is, the group-velocity dispersion is negligible), and hence, we may put δn = 0. Furthermore, if we assume that the effect of self-phase modulation is inconsequential (υnonl = 0) here, we may immediately recognize from the above expression for αn that −H″ can be associated with the Gibbs free energy, G. Upon expanding the logarithmic term in H″ to second order, we may therefore arrive at the following expression for G: G=n(ξnχ0,nψsatnχ0,nsatαnsatαnsat*)αnαn*+(12nχ0,nsatαnsatαnsat*υsat)pp+q=nαpαp*αqαn*,where ψsat= ln(nχ0,nsatαnsatαnsat*). This expression for G is very similar to the Landau expression for the free energy of a system spontaneously undergoing a first-order phase transition: The coefficient of the second-order term (αnαn*) can become negative for high values of the gain (χ0,n), signifying the onset of spontaneous symmetry-breaking, and for high values of the nonlinear saturable amplitude, υsat, the coefficient of the fourth-order term of G is also negative, thereby describing a first-order phase transition.

Emergence of sinc-shaped localized hotspot

As seen from Fig. 5A, initially, two quantum-coherently driven βn states are completely AS, but they eventually synchronize, owing to their nonlinear coupling in the pumped active region—similar to the eventual synchronization of the motion of the pendulum metronomes in Fig. 2A. As has just been demonstrated analytically above, the onset of synchronization is associated here with a nonequilibrium first-order phase transition, with the role of temperature played by the atomic inversion (50, 51). The instantaneous electric field along z at an arbitrary time instant t = 0 is given by the superposition of all βn states E(z)=in=0N1{2πωn/V)1/2ei[(β1+nδβ)z+φn]αnt=0c.c.} and is initially (before synchronization) a noise signal because φn is different for each βn state (Fig. 2B, left). To see how a nonrandom localized field pattern emerges upon synchronization (“optical morphogenesis”), we calculated E(z) at a later time instant t = t′ at which the fields at a point z = L oscillate coherently, that is, φn = φ = 0, E(z=L)=An=0N1einδβLeiβ1Lβ1=0,δβ0A0βmax=ΔβeiβLdβ=A(eiΔβL1)/(iL). Therefore, in the pumped region, the field intensity of the localized hotspot upon synchronization will be IE2Δβ2[(eiΔβL1)/(ΔβL)].[(eiΔβL1)/(ΔβL)]=Δβ2sinc2(ΔβL/2), where Δβ = βmax ≈ 40 × 106 m−1, which can also be written asI(2/Lloc)2sinc2(L/Lloc)(M2)

with Lloc = 2/Δβ being the localization length.

Note that these expressions for I and Lloc are completely independent of the shape and width of the pump and coherent-drive beams because they should be for true localization. The characteristics of the nonequilibrium-localized hotspot depend only on the intrinsic characteristics of the active nanostructure and arise from the interaction and spontaneous synchronization of the βn states. The flatter the TM2 band is in the reciprocal space [that is, the more βn states participate in the cooperative process (that is, the larger the Δβ is)], the “tighter” in space and more intense will the hotspot become. To further corroborate the validity of Eq. M2, Fig. 1 shows a steady-state snapshot from full-wave Maxwell-Bloch-Langevin time domain simulations (61, 62) of the spatial shape of the emergent localized light field, which has a sinc function dependence in the longitudinal direction. Note that this is not a potential localization: It does not arise from jumps in either the real or the imaginary part of the permittivity of the active region (see section S1).


Supplementary material for this article is available at

table S1. Comparison between the two nonpotential (nontrivial) light localization schemes.

fig. S1. Electron probability density for various potential configurations (63).

fig. S2. Schematic illustration of the considered nanosystem along with its macroscopic characteristics.

fig. S3. Negative Goos-Hänchen phase shift of a light ray upon reflection at a dielectric/plasmonic.

section S1. Supplementary text on wave localization

section S2. Effective thickness of and feedback formation in the plasmonic heterostructure

section S3. Quantum coherence–driven three-level active medium

section S4. Statistical independence of the nonequilibrium localization events

References (6367)

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Funding: This work was supported by the Office of Naval Research Multidisciplinary University Research Initiative program under Grant No. N00014-13-1-0649 and Samsung Electronics. Author contributions: K.L.T., P.K.J., and X.Z. initiated the project. K.L.T. and P.K.J. developed the theory and performed the simulations. K.L.T., P.K.J., Y.W., and X.Z. contributed to the discussions and to the writing of the paper. X.Z. supervised the work at all stages. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

Stay Connected to Science Advances

Navigate This Article