Iñigo Liberal1 and Richard W. Ziolkowski2
1Department of Electrical, Electronic and Communications Engineering, Institute of Smart Cities (ISC), Public University of Navarre (UPNA), Pamplona, Spain
2Global Big Data Technologies Centre, School of Electrical and Data Engineering, University of Technology Sydney, Ultimo, New South Wales, Australia
Quantum technologies are defined by their ability to harness quantum mechanical effects in order to achieve performances that cannot be attained within the boundaries of classical physics. In doing so, quantum technologies hold the promise to revolutionize communication, sensing, simulation and computing technologies, thus bringing transformative advances to science, industry and society [1]. Photonic quantum technologies, for which light acts a fundamental resource, are arguably the major driving force in promoting basic quantum technologies to commercial systems. For example, quantum key distribution (QKD) is the most promising example of a commercial quantum communication system [2]. QKD enables a secure channel with the guarantee that interceptions in the communication can be detected, an invaluable property to distribute cryptographic keys. By using optical links, several companies currently offer commercial QKD systems [3–5], and satellite QKD communications have been demonstrated by state‐driven enterprises [6]. Similarly, quantum sensing takes advantage of nonclassical field correlations to enable metrology with an unprecedented precision. The use of nonclassical light at the Laser Interferometer Gravitational‐wave Observatory (LIGO) is perhaps the most iconic example of a practical quantum metrology system [7]. However, sensing via the optical readout of solid state spins have enabled a myriad of applications, fascinating examples including temperature mapping within a neuron [8], intra‐bacterial magnetic field mapping [9], nuclear magnetic resonance of single proteins [10], and nanoscale magnetic resonance imaging (MRI) [11], only to name a few. Quantum computing, while still a futuristic goal [12], is the most appealing quantum technology application for the general public. Despite the fact that many fundamental challenges must be addressed, several companies are investing in the development of a practical optical quantum computer [13, 14]. In fact, the first demonstration of quantum computation advantage was based on a photonic device [15]. These are only a few examples of how quantum photonic technologies are a growing reality that will only continue to increase its impact in the following years.
Nonclassical light states are the fundamental resource driving most of the aforementioned photonic quantum technologies. For example, single‐photon states are the basis of QKD communications systems, squeezed light states and N00N states empower sensing below shot noise, and single‐photon and cluster states give power to linear‐optics quantum computers. Therefore, the design of light sources capable of generating prescribed nonclassical light states is one of the most important aspects in developing photonic quantum technologies.
In an entirely different field, the vast body of research nowadays known simply as antenna theory has led to a high degree of sophistication in the ability to tailor the radiation characteristics of electromagnetic waves emitted by innumerable types of sources. Previous chapters of this book have focused on classical antennas and antenna arrays; they exemplify how this is indeed the case. Without any doubts, much of this work can be applied to catalyze the design of nonclassical light sources. Based on this perspective, and following the general themes of this book, we will discuss in this chapter how photon statistics can be engineered by designing the geometry of quantum antenna arrays. The main goal is to summarize the basic theory behind quantum antenna arrays, providing antenna designers with the tools to explore how their know‐how can be applied to the exciting field of quantum technologies. In doing so, we will explore intriguing nonclassical light sources with no classical counterparts, including directionally entangled photon bunches and perfectly isotropic single‐photon sources.
To achieve our goal, we will begin with an overview of quantum antenna array theory. While arrays of electrically small classical antennas were described in Chapter 10, the radiating elements herein are physically extremely small being quantum emitters, i.e. atoms. We then follow with a relevant introduction to the concepts of photon statistics, emphasizing directional correlations. These two fields are then combined to reveal how the design of linear arrays of quantum emitters can enable beamforming of their photon statistics. We emphasize the effects that can and cannot be observed in classical systems. Examples of first, second and Nth order correlation outcomes demonstrate how their directivities can be tailored by a quantum array's configuration. We apply these results to understand how one could in fact achieve an isotropic single‐photon source. We also address the difficulties associated with the very contrasting concept of quantum superdirectivity. A description of the technologies available today that are the most promising to realize these exciting quantum array results is given. Finally, we explore the many possible future directions of the field of quantum antenna arrays and the concurrent opportunities for electromagnetic scientists and engineers.
We start by introducing the basic quantum theory of radiation for an array of quantum emitters. We try throughout to establish the connections with classical antenna array theory that was the focus of most of the previous chapters of this book. In particular, the final expressions are written in forms that are readily connected to their classical counterparts.
As a first approach to quantum antenna arrays, we analyze the general scenario schematically depicted in Figure 10.1. An array of quantum emitters is depicted that is coupled to a macroscopic photonic environment through the fields radiated by its individual elements. This complex environment is characterized by a relative permittivity distribution, . Quantum excitation in such macroscopic environment will be modeled as a bath of polaritonic modes (i.e. a polariton is a quasi‐particle formed when a photon couples strongly with an exciton, e.g. a phonon, and, hence, it represents a joint state of light and matter) following the macroscopic quantum electrodynamics (QED) formalism (see, e.g. [16–18]). This approach allows for the analysis of lossy and dispersive media. We restrict the analysis in this chapter to the case of isotropic dielectric media, which is typical of most configurations at optical frequencies. However, generalizations of the same quantization procedure to account for magnetically polarizable, bi‐anisotropic, and nonlocal media are available [19]. While most quantum optics textbooks treat the decay of quantum emitters from the Schrödinger and/or interaction pictures, here we will adopt the Heisenberg picture because fully time‐evolving field operators bring a closer analogy to classical antenna theory.
The Hamiltonian of the system depicted in Figure 12.1 is given by the addition of the Hamiltonians for the array of quantum emitter's and the electromagnetic field plus polarizable matter (i.e. the polaritonic) system, as well as the interaction Hamiltonian between both systems:
The emitters are modeled as two‐level systems that are located at positions and have the transition frequency and dipole moment . This choice leads to the Hamiltonian
where is the Pauli operator with .
Excitations in the macroscopic photonic environment are modeled as a continuum of polariton modes for each position , frequency , and polarization with annihilation operators . They can also be written in vector form as and obey equal‐time bosonic commutation relations:
In this manner, the Hamiltonian for the macroscopic photonic environment can be compactly written as
Finally, we consider an interaction Hamiltonian within the electric dipole approximation
where , with and being the annihilation and creation operators associated with the nth emitter and with . Note that because we are assuming point‐like emitters with an electric dipole transition, the electric field operator is evaluated at the position of the emitters . This approximation is similar to assuming that the elements of a classical antenna array are electric Hertzian dipoles. However, similar to antenna theory, it is possible to further control the quantum radiated field by engineering the current distributions of the individual radiating emitters [20].
The electric field operator is given by
where we have introduced the response function
which is proportional to the classical dyadic Green's function [21]. Equation (12.6) intuitively conveys the idea that the electric field operator can be understood as the result of radiation by the polaritonic excitations in the macroscopic environment. Therefore, it is not surprising that the classical dyadic Green's functions play the role of a propagator in such a process. However, the antenna designer might be surprised that the propagation term also includes the square root of the imaginary part of the permittivity as a coupling term between fields and environmental excitations. Intuitively, it makes sense that the effective exchange of energy between fields and matter must be mediated by dissipation. From a more rigorous perspective, this term is exactly what is needed to reconstruct the macroscopic Maxwell's equations from the coupling of the free‐space fields and matter when the quantum mechanical canonical quantization of the fields is required and employed in the presence of matter [17]. The presence of the square root of the imaginary part of the permittivity also ensures that the fluctuation dissipation theorem, as it is often invoked in the semiclassical treatment of thermal emission [22], directly emerges from this framework.
The Hamiltonian formulated in Section 12.2.1 is the standard formulation of macroscopic QED [16–18]. Next, we derive an expression for the radiated field in the far zone that has very clear analogies with antenna theory. To this end, we compute the Heisenberg equation of motion, , for the polaritonic operator , and integrate it, leading to its decomposition into the free‐evolving and source parts
Physically, the free‐evolving part corresponds to the electromagnetic field and matter evolution that would be taking place in the absence of the interactions with the quantum emitters. It can be written as:
On the other hand, the source part corresponds to the excitations of the electromagnetic field and matter system induced by the interaction with the quantum emitter system. It is given by
with
By introducing (12.8) into (12.6), we find that the electric field operator can be similarly decomposed into free‐evolving and source driven parts
with
Introducing (12.11) into (12.6), using the completeness relation of the dyadic Green's function to compute the volume integral, and noting that is an odd function in , the contribution of each emitter to the source field can be written as a memory kernel acting on the electric dipole operator
with
Recalling that the electric dipole operator is , the decomposition of the source electric field operator into positive and negative frequency components is a useful approach:
with
and
Physically, the positive frequency component is associated with the radiation of one photon by the emitter (as a transmitter), while the negative frequency component represents the absorption of one photon by the emitter (as a receiver). A very compact expression for the electric field operator can now be obtained by applying the Laplace transform to these operators:
and then evaluating the integral over all frequencies via complex analysis techniques. The result is that Eq. (12.17) can be written as
The formulation up to this point has been developed for an arbitrary photonic environment. Particular cases can be addressed by introducing a specific form of the dyadic Green's function. Thus, we assume that the emitters are immersed in free‐space and the observation point for the electric field is in the far‐zone of this system to establish a closer connection with antenna theory.
The classical dyadic Green's function for this case is
with , , , and the unit vector . The first‐order approximation to is given by the well‐known expression: . We keep this first‐order approximation in the exponential and only the zero‐order approximation in the other terms. The dyadic Greens function then becomes
with . In this manner, we can rewrite the electric field operator as
Finally, we apply the inverse Laplace transform and get the general form for the source electric field in the far zone
with being the delayed time with respect to the origin of the coordinates.
Equation (12.24) describes the electric field operator in the far‐zone for an arbitrary time evolution of the quantum emitter operator. This formalism can be applied to a variety of configurations. We particularize it here to the case of the decay dynamics of noninteracting quantum emitters in the weak coupling regime. Thus, we will show that it exhibits a simple exponential decay. This case corresponds to the physical situation in which each emitter operates as an independent single‐photon source.
First, we analyze the Heisenberg equations of motion, , for the quantum emitter annihilation operator . They lead to the relation for the th emitter:
Next, we clarify the approximations involved in assuming a simple exponential decay rate.
We begin by assuming that the macroscopic photonic environment is initially in the vacuum state. We can then substitute for all relevant interactions. From an antenna theory perspective, this condition is similar to assuming that there are no external sources to the problem. In addition, we assume that the emitters are not interacting amongst themselves. This allows us to approximate . This is the typical assumption in antenna theory where each element of the array is driven independently, and the array geometry is designed so that the coupling between each of the antenna elements in the array is weak enough.
We also introduce the low‐excitation approximation (also known as the one‐photon correlation approximation). It consists of replacing . This approximation is similar to having started with a Hamiltonian within the rotating wave approximation. For the problem of initially‐excited noninteracting quantum emitters, this means that each emitter effectively sees only a single excitation. It is a very accurate approximation.
Even with all of these approximations, the model still is capable of describing both the weak and strong coupling regimes. These are the ones in which most optical systems operate. However, the approximation will not be valid in the ultra‐strong coupling regime where the number of excitations is not preserved.
With these approximations and the application of the Laplace transform, the equation of motion for reduces to
Using (12.20) and rearranging its terms we find the following expression for the Laplace‐transformed emitter operator:
where we have introduced the emitter's self‐energy :
Finally, the standard Born–Markov approximation consists of neglecting the dispersion of the self‐energy, i.e. . This approximation is accurate when the coupling between the emitter and the photonic environment is small. This situation is usually referred to as the weak‐coupling regime. The approximation is also equivalent to assuming that the interaction between the emitter and the photonic environment has no memory effects.
Within this approximation, the imaginary part of the self‐energy corresponds to the decay rate associated with the exponential decay. On the other hand, the real part of the self‐energy results in a frequency shift, .
Note that the real part of the dyadic Green's function diverges. As a consequence, the self‐energy must be properly renormalized to properly compute this frequency shift which is well‐known as the Lamb shift. However, since this small shift does not affect the overall decay dynamics, it is typically assumed that it was included in the transition frequency of the emitter, . Finally, within all of these approximations, we recover the usual exponential decay expression:
After detailing all of the approximations leading to the simple exponential decay, it might seem that it is a crude approximation that can only be applied in very restricted circumstances. However, the truth is that for most weekly interacting single‐photon sources, exponential decay is a very accurate approximation. On the other hand, observing non‐exponential decay dynamics is a challenging experimental task in most configurations.
Finally, by introducing (12.29) into (12.24) we find the expression for the positive frequency electric field operator that we will use in the rest of the chapter:
The term
where
The theory detailed in Section 12.2 shows that the temporal and spatial variations of the electric field operator are the same as those of a classical field. This outcome suggests that many aspects of the theory learned for the design of classical antenna arrays could be transferred to quantum antenna arrays. At the same time, it is expected that quantum antenna arrays will give access to effects that cannot be observed in classical systems. This is indeed the case.
Quantum effects fundamentally arise in the measurement of photon statistics, i.e. in coincidence photon measurements. Early experiments in the 1970s and 1980s confirmed that the measurements of photon statistics cannot be explained by a classical theory. They revealed numerous novel phenomena including antibunching [23], two‐photon cascade emission [24] and violation of Bell's inequalities [25]. These effects provide the most direct evidence that a quantum theory of light is needed. Moreover, correlations that cannot be achieved by classical systems are the major driving force behind the currently evolving quantum technologies. Being the most direct experimental proof for quantum mechanics, early experiments have been very critically analyzed, motivating researchers to identify loopholes invalidating the conclusions extracted from them [26–28]. Nevertheless, the improved and loophole‐free experimental setups continue to confirm the quantum theory [29–31].
Following the theme of antenna arrays, we outline in this section the basic theory needed to describe how the geometry of a quantum antenna array influences directional correlations in photon statistics. To this end, we assume that the quantum antenna array is surrounded by a number of photon counters located at positions , . We also assume that these devices operate via direct photon detection, though alternatives such as homodyne detection might be considered [32]. Therefore, the detectors can be phenomenologically modeled by using Glauber's theory of optical coherence [33, 34], which connects the probability of photon coincidence detections with field correlations.
Specifically, a measurement event can be modeled as a change in the quantum state of the detector from to , triggered by the absorption of a photon. Since the annihilation of environmental excitations is described by the positive frequency electric field operator, the probability amplitude for a photon with polarization to be detected at time and location is . Consequently, the total probability density of detecting a photon is proportional to the sum over all of the probabilities for each final state:
where we have made use of the identity operator to obtain the intermediate expression.
For a set of detectors, the procedure outlined earlier can be generalized to describe the coincidence measurement of photons with polarizations, at different times , and for different positions :
If we assume that the detectors are in the far‐zone of the array at positions it is clear from Eq. (12.31) that the photon statistics are not affected by the separation distance . Therefore, it is possible to define the probability density per unit and as:
The probability density described by Eq. (12.35) enables the study of interesting spatio‐temporal correlations. Importantly, it can be shown that any other photon statistic can be constructed from such a probability density [35, 36].
As classical antenna array theory and its practical applications emphasize, directivity is a major performance characteristic associated with any array. Consequently, we are most interested in the directional correlations that can be induced by the geometry of a quantum array. In order to emphasize the directional behavior of its emissions, independent of the times at which the measurements take place, we define the time integral of the probability density:
The integrated correlations represent the average number of ‐photon coincidence measurements for any time delay. They are a good figure of merit of how “bunched” the photons are in different directions after their emission.
The functions (12.36) include time integrals over rapidly oscillatory functions that are related to the transition frequency of the dipoles, . Unless the transition frequencies are very close to each other, these integrals will average to zero over time. As a consequence, interference phenomena might only be observable for short time intervals. In general, this property indicates that no quantum interference of any significance will take place unless the emitters have similar transition frequencies, i.e. . In such a case, the time‐integrated correlation (12.36) reduces to:
with . With this general form, time‐integrated correlation functions of arbitrary order can be evaluated for any array geometry and initial state of the emitters in the quantum array.
We illustrate the applicability of the theory outlined earlier by studying the photon statistics of a uniform vertical linear array of quantum emitters, an archetypical geometry in textbooks of antenna theory [37–40]. Specifically, we consider the geometry schematically depicted in Figure 12.2, where all of the emitters are located at the positions along the ‐axis with equal separation distance , and each has its transition dipole moment also oriented along the ‐axis, i.e. .
For this particular configuration, the far‐zone electric field operator has a single polarization component . Then, Eq. (12.24) reduces to
where . In the same manner the time‐integrated correlation function (12.36) can be compactly written as
where we have defined a generalized ‐order quantum array factor:
Similar to the classical array factor, the generalized array factor describes how the geometry of the array affects the photon statics of different order. For example, the case corresponds to the average number of photons measured in a given direction. The case provides a good figure of merit of how directionally bunched the emitted photons are. When each emitter acts as a single‐photon source, for .
The first‐order correlation, , physically predicts the average number of photons measured by a detector placed along the direction . The generalized array factor for reduces to
Equation (12.41) reveals that the first‐order generalized quantum array factor is very similar to the magnitude squared of the classical array factor [37–40]. This is not surprising since the classical power emission pattern should correspond to the average number of photons. At the same time, however, there is an important difference that results in nonclassical effects even in the average number of photons. Specifically, the initial‐time correlation plays the role of the signals feeding the th and th radiating elements. On the other hand, the importance difference is that the initial‐time correlation is not always necessarily factored as the product of two ‐numbers, e.g. as . As a consequence, access is granted to configurations forbidden in the classical case.
In the same manner that classical antenna arrays control the emission pattern by changing the magnitude and phases of the signals driving their antenna elements, quantum antenna arrays control the emission pattern by controlling the initial state of the system, i.e. how the array is initially excited to start the emission process. There are several experimental techniques that enable the preparation of a given initial state of a quantum system. For example, initialization laser pulses and incoherent pumping can be employed for initial states corresponding to product states (see, e.g. [41]), where each emitter is individually excited. More complex initialization schemes are required for entangled states. Specific phase profiles for them could be implemented, for instance, if they are based on laser‐assisted interactions [42, 43]).
We consider the single‐excitation state:
as an example of an initial entangled state. This is actually a maximally entangled state in which a single excitation is shared by all the elements of the array. By placing the same phase on each emitter, the result is similar to that of a classical array with a uniform excitation. In fact, we recover a factorization analogous to the classical case . The first‐order generalized array factor then reduces to
Despite the fact that we are using a nonclassical, maximally entangled state, it is found that the radiated field can be beamformed exactly as in the classical case. The probability amplitudes play exactly the same role as the currents driving the classical antenna elements. This is an exciting outcome because it means that all the machinery developed for the synthesis of classical antenna arrays can be directly applied to single‐excitation state initiated quantum arrays.
Next, we consider the initial state in which all emitters are initially excited, i.e.
Despite the fact that it is a factorized (not entangled) state, this configuration is very relevant from a practical standpoint. It represents the case in which all emitters are independently excited, which could be implemented in practice with electrically driven elements. Moreover, it leads to nonclassical photon statistics. Specifically, the correlator for this initial state reduces to . Thus, the generalized array factor takes the very simple form:
It can be concluded from this simple expression that the emitters do not interfere and, hence, do not modify the directional properties of the average number of photons. In other words, the emission pattern in terms of the average number of photons is identical to the individual emitter pattern, i.e. a function, and is independent of the number and geometrical arrangement of the emitters. The fact that the emission pattern is independent of the number of emitters and the geometry of the array is a purely quantum effect with no counterpart in classical antenna theory.
In order to illustrate this point, Figure 12.3 depicts the first‐order time‐integrated correlation for a uniform linear vertical array composed by identical emitters separated by one wavelength, i.e. for the emitter locations . The patterns are only depicted in the ‐plane due to the fact that the geometry is rotationally symmetric around the ‐axis. Classical linear vertical arrays are known for their enhancement of the directivity, focusing the emitter radiation into the plane perpendicular to the array axis (i.e. the broadside direction ) [37]. Following the theory earlier, Figure 12.3 confirms that the emitted photons will be constrained to an narrower set of directions than in the single emitter case (shown as a dashed black line) when the array is initially driven into a single‐excitation state. On the other hand, it is found that the emission pattern is identical to that of the individual elements when all three emitters are initially excited. This outcome confirms the peculiar quantum effect that the directionality on the average number of photons in the case in which all of the identical emitters are initially excited is not affected by their array configuration.
It can be concluded from the two examples earlier that it is possible to beamform photon statistics by selecting how the array is initially excited. While we have shown here the two main examples of single‐excitation and ‐excitation states, many other options are available and could be explored. Based on similar considerations, the possibility of designing a perfectly isotropic antenna is discussed in Section 12.5. Similarly, the difficulties in harnessing these extra degrees of freedom to obtain quantum superdirectivity beyond classical limits are summarized in Section 12.6.
The second‐order () time‐integrated correlation for the same linear array takes the form:
By definition, it represents the average number of two‐photon coincidences for any time delay. In practice, it is a good figure of merit of how directionally bunched the emitted photons are. Specifically, since Glauber's probability densities are nonexclusive quantities [35], implies that no photon coincidences of any order will be recorded along the directions and for the same decay process. Consequently, if has a highly directive pattern, it indicates that all of the emitted photons are bunched into a narrow set of directions.
For this reason, engineering a quantum antenna array with a highly‐directive second‐order correlation would lead to a nonclassical light source whose photons are emitted in directionally entangled bunches. Similar to the first‐order array factor, the second‐order array factor, , can be controlled by the initial state of the system, which affects the correlation function, , as well as by the geometry of the array.
Importantly, the angular patterns of the second‐order correlations do not have to match those obtained for the first‐order correlation. They can have an entirely different behavior. In order to illustrate this point, Figure 12.3b depicts the time‐integrated second order correlation for the same initial configurations associated with Figure 12.3a. To obtain these results, we fixed the first evaluation direction to be and then evaluated the correlation function as a function of .
First, for the single‐excitation state the initial‐time correlator trivially reduces to . Consequently, the second order correlation is identically zero for all evaluation angles, i.e. . This nonclassical effect is the directional counterpart to temporal antibunching. It describes the impossibility of measuring two excitations from a single‐excitation state. As was noted, for within our model. From a practical standpoint this condition implies that two different detectors would never detect a photon on the same decay run if the system was initialized with a single‐excitation state.
In contrast, we obtain for the ‐excitation state a nontrivial initial‐time correlation
Then, the second‐order generalized array factor can be written as:
In this manner, we find that there are directional effects in the second order correlation factor for an array in which all emitters are independently excited. They occur despite the absence of any directionality in the first‐order correlation. Therefore, it can be then concluded that the photons are probabilistically emitted with no directional preference other than the pattern of the single emitter, i.e. with no directionality for . On the other hand, the photons will all be measured bunched around a set of directions for a given decay process, i.e. with the directionality of .
This effect is somewhat analogous to two‐photon interference in beamsplitters. The two photons are equally probable to exit this system into any of the two output ports. Nevertheless, both are certain to exit into the same output port [32]. Here, a number of photons are equally probable (up to the individual emitter pattern) to exit the system on a continuum of output channels (i.e. the directions), but they are certain to be measured within a narrow interval of directions. This effect can be more clearly appreciated in Figure 12.4. The second‐order correlation is presented for three different given by (), (), and (). The angular patterns consist of lobes centered around each , plus additional grating lobes.
These photon statistic outcomes can then be tailored by changing the geometry of the array. To this end, different design techniques could be developed. They are very likely to exhibit both similarities and differences with respect to the synthesis techniques developed for classical antenna arrays. In order to illustrate this point, Figure 12.5 depicts as a function of the number of emitters and the distance of separation between them. In analogy with classical antenna arrays, we find that increasing the separation between the emitters results in the appearance of grating lobes, i.e. the emitted photons will be bunched around not one, but a set of multiple lobes. However, in contrast with classical antenna arrays, it is found that the directivity does not increase monotonically along with the number of emitters. Conversely, the results in the superposition of the single emitter pattern with that of a highly directive lobe. These results suggest that nonclassical techniques will have to be developed for the synthesis of quantum antenna arrays and each of their different photon correlations.
Examining more closely the generalized array factor Eq. (12.40), it is clear that computing the photon statistics becomes a very challenging task as the order increases. Having addressed the first‐ and second‐order correlations, here we address the ‐order photon correlations for an antenna array in which all emitters are initially excited, i.e. for which
It is the highest‐order nonzero correlation for an array of single photon sources. Interestingly, a tractable expression can be derived for this extreme case.
The main difficulty in computing the photon statistics stems from evaluating the ‐order correlator:
In order to evaluate the correlator, let with be a basis for all possible states available to the array of emitters. This basis can be arbitrary, but for convenience we define the first state of this basis as all emitter's being in the ground state, i.e. . Next, we find the following identity
We then make use of the identity operator to evaluate Eq. (12.52) as follows
where we define for and in the rest of the cases.
With this result, the generalized quantum array factor can be written as
Consequently, the definition of implies that the sum in Eq. (12.53) runs over all possible permutations of the position vectors of the following product of exponentials:
For the example of three emitters used in this chapter, the generalized array factor is explicitly:
This result can be understood intuitively in terms of the “which‐path” information that is available in the state of the quantum emitter subsystem. After measuring a ‐photon coincidence event, the emitter's subsystem must necessarily be in the ground state, . Therefore, the state of the emitters contain no information in terms of what emitter radiated the photon measured at a given detector. In other words, there is no which‐path information and, hence, quantum interference can occur. In this manner, each of the permutations in (12.40) can be understood as a different path that connects the detection in each direction to the emitter in the location. For these reasons, computing the photon statistics become “simpler.” This result also shows that the directivity of the photon statistics is expected to increase along with their order since less and less which‐path information is available as more photons are measured (see Figure 12.3).
We provide in this section yet another example of a quantum antenna device with functionalities that cannot be obtained with a classical antenna system. In particular, we demonstrate that it is possible to have a perfectly isotropic and temporally coherent single‐photon source with full‐polarization coverage [44]. Isotropic antennas are used as a reference in the very definition of directivity. However, it is a well‐known result from classical antenna theory that perfectly isotropic and temporally coherent radiators are not possible [45–48]. This conclusion is best visualized by the hairy ball theorem [49, 50], an algebraic topology result stating that there is no nonvanishing continuous tangent vector field over the entire surface of a sphere. In other words, “you cannot comb a hairy ball flat without creating a cowlick.” Therefore, perfect isotropy with a well‐defined polarization is impossible.
Despite this impossibility, antenna designers continue to seek efficient quasi‐isotropic antennas that approximate the properties of an isotropic radiator [51–54]. The motivation behind quasi‐isotropic antennas is that they provide full‐coverage, an invaluable property for any wireless technology in which transmitting and receiving devices cannot be aligned [55–57]. As a separate note, isotropic sources do exist in the form of thermal radiators [22], and multifrequency systems [58], at the cost of being temporally incoherent. Beyond electromagnetic systems, isotropic sources do exist in other physical systems such as in acoustics [59].
Nevertheless, an isotropic single‐photon source can be designed by using the emitter schematically depicted in Figure 12.6. It consists of a multi‐level system with a single excited state , and different ground states . The transition from the excited state to each ground place is characterized by the transition frequency and the dipole moment . The electric field operator takes the same form as that derived for an array of quantum emitters, albeit with all the emitters being located at the same position. Consequently, Eq. (12.38) reduces to
The probability density of measuring a photon with polarization at the position and time is given by Eq. (12.33). For the multi‐level emitter under consideration it can be explicitly written as
Again, the nonclassical properties of the emission pattern in the proposed configuration arise from the correlators in Eq. (12.57). In contrast with the classical case, the correlation does not need to be decomposed into the product of two complex numbers, opening up additional degrees of freedom. For a multi‐level emitter with a single excited state and multiple ground states, we have that . Interestingly, the result of such a correlator is that there is no interference between the field operators associated with each transition. The absence of interference can be intuitively understood as the result of having distinguishable ground states and thus access to perfect which‐path information. In this sense, the distinguishable ground states play the same role as passive observers similar to popular effects such as two‐slit experiments [60] and quantum erasers [61, 62].
The lack of interference is a nonclassical property arising from the superposition of distinguishable decay paths. It can be harnessed to design an isotropic single‐photon source. To this end, we particularize the multi‐level emitter to the case with three degenerate ground‐states with the same transition frequency () and with transition dipole moments of the same magnitude (), but oriented along each of the Cartesian axis, i.e. , and . For this configuration, the field correlation reduces to
In order to show that the emission will be isotropic for all arbitrary polarizations, we recall that an arbitrary polarization basis can be obtained from the usual basis via rotation by an angle and the introduction a phase shift (see Figure 12.6):
By using this decomposition, the probability density of measuring a photon with arbitrary polarization at time and location is found to be:
Equation (12.61) reveals that the probability density of detecting a photon is independent of the measurement direction, which is described by the angles and . Therefore, it can be concluded that the emission is isotropic, i.e. the same probability is obtained for all directions. Moreover, since the probability density is independent of the choice of the polarization vector , it can also be concluded that the same result holds for all polarizations. In this manner, it is demonstrated that a perfectly isotropic and with full polarization coverage single‐photon source does exist.
The emission from an isotropic single‐photon source resembles black‐body radiation from a heated object, i.e. the emitted radiation is both isotropic and unpolarized. The main difference is that isotropic single‐photon sources are coherent radiators. Thus, the emitted light corresponds to a single photon with a well‐defined wavepacket. For this reason, it is a coherent source for which interference phenomena can be observed. For example, interference fringes will appear if the emitted photon is sent through a Mach–Zehnder interferometer. In conclusion, isotropic single‐photon sources are a quantum light source with properties that cannot be obtained with classical antennas.
It was shown in earlier sections of this chapter that quantum antenna systems provide access to radiation patterns as the average number of photons per direction that cannot be obtained with classical antenna systems. This poses the question of whether or not quantum sources can be used to increase the directivity of antenna systems beyond their classical limits, leading to quantum superdirectivity. A simple proof seems to indicate that this is not possible [63]. However, researchers are known for relentlessly searching for ways to circumvent fundamental limits, and the situation might change in the future. Since the proof provides an additional insight into the mode of operation of quantum antenna arrays, we quickly reproduce it here.
For the sake of simplicity let us consider here an array of all identical emitters, even though the proof can be generalized to more complex scenarios. The starting point is then the first‐order time‐integrated correlation function – the average number of photons measured along the direction – discussed previously:
where is the first‐order generalized array factor:
The classical case is recovered for . On the other hand, does not need to be factorized as the product of two complex numbers for a quantum antenna array. Again, as a consequence, it provides additional degrees of freedom.
It was shown in Section 12.5 that such freedom can be used to design an isotropic single‐photon source. At the same time, the initial‐time correlation has some restrictions that limit the achievable radiation patterns. To this end, we define a matrix with elements and a vector with elements . We use them to express in matrix form as:
In contrast, the radiation pattern of a classical antenna array would be restricted to the more constrained form
Nevertheless, quantum antenna arrays also have some restrictions of their own. First, we note that and, hence, is a Hermitian matrix. In addition, since is a positive real number, we then have that is a positive semidefinite matrix. Importantly, this implies that all its eigenvalues are positive. Next, let and , , be the eigenvalues and associated eigenvectors of . Then the spectral decomposition of is given by
Subsequently, can be rewritten as:
Equation (12.67) reveals that the quantum generalized array factor describing the average number of photons measured as a function of the direction can be written as a linear combination of classical power patterns for the same array, but with different driving coefficients. From this perspective, the power behind quantum antenna arrays is that they enable the superposition of classical power radiation patterns. However, since is a positive definite matrix, its eigenvalues must be positive. Therefore, there cannot be destructive interference between the classical radiation patterns. This fact suggests that the directivity cannot be increased by the superposition. In many senses, the situation is very similar to polarization considerations in the maximization of the directivity of classical antennas [64]. The polarization degrees of freedom provide orthogonal channels in which the power can be distributed. However, they cannot be used to increase the directivity.
To demonstrate this explicitly, we first define the quantum directivity as follows:
Next, if we optimize the coefficients, we would obtain the same optimal coefficients for any , and therefore the sums in the numerator and denominator would compensate in obtaining the maximal directivity. Consequently, the directivity of a quantum antenna array in terms of the average number of photons cannot exceed the upper limit for a classical antenna array.
Having described the basic theory and its particularly interesting outcomes, we discuss in this section those material platforms that are currently the most promising candidates for implementing quantum antenna arrays. First, the basis of a quantum antenna array is the single‐photon sources that act as their individual elements. The history of single‐photon sources dates back to the 1970s with the first demonstration of a single‐photon source based on a beam of sodium atoms [23]. This effort was followed a decade later by experiments with trapped ions [65]. The first solid‐state single‐photon source was demonstrated in the 1990s using a single molecule dye trapped in a solid [66]. Since then, research on single‐photon sources has exploded, leading to their rapid development. Today, there is a variety of solid‐state single‐photon sources in the form of quantum dots (QDs) [67–70]; color centers in semiconductors, such as diamond, silicon carbide, and silicon [71–76]; and defects in two‐dimensional (2D) materials [77–80]. In general, solid‐state sources suffer from their coupling to a complex environment, e.g. phonons in a crystal lattice, charge accumulations, and interactions with nearby spins. This coupling limits their coherence times. Moreover, it poses difficulties in the fabrication of identical emitters. However, the quality of solid‐state single‐photon sources has been steadily improving and is rapidly approaching the required standards in terms of purity, indistinguishability and on‐demand operation.
Free from a disturbing environment, ultra‐cold atoms [81, 82] and trapped ions [83] are excellent single‐photon sources. Atomic‐based sources provide better single‐photon performance, longer coherence times and a finer control over the photon's wavepacket than their solid‐state counterparts do. Nevertheless, the integration of atomic‐based devices into commercial systems is likely to be restricted to very high‐end applications. Superconducting circuits [84, 85] exhibit the best performance in terms of the coupling strength to the environment and their coherence times. However, the low energy of microwave photons forces these circuits to operate at very low temperatures. The entire system must be contained within a cryogenic setup. Furthermore, nonlinear processes such as spontaneous parametric down conversion were identified early as very high quality nonclassical light sources [86–88], albeit at the cost of being bulky and nondeterministic in nature.
It must be noted that single‐photon sources seldom consist of a quantum emitter alone. The pioneering work of Purcell [89] showed that emission processes are affected by their photonic environments. Subsequently, the use of nanostructures to tailor light‐matter interactions and to enhance the performance of single‐photon sources has been aggressively pursued [69, 90]. Different strategies include the use of resonant cavities [91], a field usually referred to as cavity QED [90]. The resonant interaction with a cavity mode increases the brightness and efficiency of the source, locks the emitted photons into a well‐defined mode, and empowers the observation of reversible decay dynamics in the strong coupling regime. Commonly used resonators include nanopillars [92], nanodisks [93], photonic crystal [94, 95], superconducting [96], and plasmonic [97, 98] cavities. More recently, attention has been partially shifted to the integration of quantum emitters into waveguides [99–101]. Waveguides exhibit advantages similar to cavities in terms of the coupling efficiency to a specific mode and the Purcell enhancement associated with slow light, while at the same time providing intrinsic long‐range interactions that enable the exploration of collective effects [102, 103]. Light emission can not only be enhanced, but also suppressed by the environment. Examples include photonic crystals that exhibit a band gap [104, 105], closed cavities far from resonance [106, 107] and near‐zero‐index (NZI) media [108–111] which lead to the excitation of bound states [112–114] and the associated many‐body physics [115].
Beyond the design of individual emitters, the deterministic fabrication of arrays of quantum emitters with control over their position, orientation, and emission frequency still poses important technological challenges. Again, the most advanced technology is that of atomic‐based emitters. Intensive research on the trapping and cooling of atoms in atomic physics has enabled the deterministic generation of arrays of quantum emitters with hundreds of elements, including one‐dimensional (1D) [116], two‐dimensional (2D) [117], and three‐dimensional (3D) [118, 119] topologies with quite arbitrary geometries. Additional advantages of atomic‐based systems are that all emitters are ensured to be identical, and the initialization and coherent control of the quantum state of the array is much more advanced. We believe that all of the theory presented in earlier sections of this chapter is within the reach of atomic‐based technologies.
Unfortunately, atomic physics setups are not a viable option for most practical applications, and research is needed in the fabrication of quantum antenna arrays in solid‐state platforms. However, the fabrication of arrays of solid‐state emitters is much more challenging. As solid‐state emitters mostly consist of defects in a solid, the interaction with the complex solid‐state environment makes it very difficult to obtain identical emitters with the same orientation, emission wavelength, linewidth, etc. The accurate positioning of the emitters is also a challenging task. Nevertheless, there have been some promising demonstrations of arrays of color centers based on selectively growing on arrays of diamond nanopillars and nanopyramids [120, 121]. Moreover, emitters based on atomically‐thin semiconductors, e.g. tungsten diselenide () and hexagonal boron nitride (hBN), are currently showing the most promising advances toward the deterministic fabrication of arrays of solid‐state quantum emitters. In an atomically‐thin semiconductor, an emitter at a specific location can be induced by using nanostructures producing local strains [78, 122, 123], as well as by laser [124] and electron beam [125] irradiation. With these fabrication techniques, we believe that planar quantum antenna arrays based on atomically‐thin semiconductors should be feasible in the short term.
This chapter has outlined the basic theory of quantum antenna arrays, establishing connections with the antenna theory field whenever possible. While multiple opportunities can be identified at the fundamental and technological levels, quantum antenna arrays are still a nascent field, and it is difficult to predict how it will evolve in the following years. At the fundamental level, it has been demonstrated that quantum antenna arrays enable directional correlations that cannot be achieved with classical systems. However, the complete set of available correlations in the quantum regime is not known, and their impact on physical observables is yet to be discovered. Note that we have used a simple two‐level system model for the emitters. We have not discussed how information would be encoded in the photons emitted by them. Research on more complex emitters, e.g. emitters with a non‐dipolar current distributions and/or multi‐particle transitions, is likely to facilitate even better performance characteristics and other nontrivial phenomena. Information on the emitted photons might be encoded by using polarization and/or temporal control over them. At the same time, the directional correlations induced by the initial state of the systems are a form of embedded information that could be investigated for some practical information exchanges, i.e. communication‐enabled applications. In general, we anticipate that in the next years there will be a large number of theoretical efforts both at the very basic level, as well as in the development of design and synthesis procedures for quantum antenna arrays.
From a practical standpoint, technological applications of quantum antenna arrays are yet to be explored. Classical antenna arrays radiate into free‐space, being the basis of wireless technologies. However, current photonic quantum technologies are mostly based on guided platforms such as optical fibers and photonic integrated circuits. The main reason for the predominant use of guided setups is that these quantum optical systems utilize single or a very reduced number of photons. Therefore, propagation photon losses can be prohibitively high for most applications. Nevertheless, aggressive investments on quantum experiments at the space scale are being funded by state agencies, including the recent demonstrations of satellite QKD distribution, seem to indicate that there will be a future for wireless quantum communications.
We expect that the benefits of wireless connections will also be exploited at small scales and within integrated systems. For example, quantum sensing technologies based on solid‐state spins rely on their optical readout. Currently, optical readout tasks are carried out with bulky far field optics, an approach that is incompatible with a practical implementation outside the laboratory. They need to be addressed in a wireless fashion. Optical antenna arrays might solve this problem by enabling compact and efficient interfaces between solid‐state spins and integrated photonics. Moreover, with the size of solid‐state spins being at the atomic scale, dense arrays of spins cannot be individually resolved with far‐field optics. Having the need to address a sample with many of these emitters, quantum antenna array theory is likely to be needed for the design of advanced optical interfaces, where quantum correlations might be harnessed to improve the performance of their readout.
These are only a few examples of how future directions might look like for quantum antenna arrays. We re‐emphasize that quantum antenna arrays are an emerging field of science and engineering. Many aspects and applications are yet to be explored both from the fundamental and technology perspectives. There are numerous opportunities for electromagnetic researchers with quantum analysis, simulation, and measurement experiences to contribute their skills to this nascent field.