The classification of higher-order photon emission becomes important with more methods being developed for deterministic multiphoton generation. The widely used second-order correlation g(2) is not sufficient to determine the quantum purity of higher photon Fock states. Traditional characterization methods require a large amount of photon detection events, which leads to increased measurement and computation time. Here, we demonstrate a machine learning model based on a 2D Convolutional Neural Network (CNN) for rapid classification of multiphoton Fock states up to |3⟩ with an overall accuracy of 94%. By fitting the g(3) correlation with simulated photon detection events, the model exhibits an efficient performance particularly with sparse correlation data, with 800 co-detection events to achieve an accuracy of 90%. Using the proposed experimental setup, this CNN classifier opens up the possibility for quasi-real-time classification of higher photon states, which holds broad applications in quantum technologies.

Recently, quantum light has been undergoing rapid advancements and playing a pivotal role in the development of quantum technologies.1–5 At the quantum level, the photon nature of light is characterized as discrete packets of energy,6,7 offering remarkable precision, sensitivity, and enhanced communication security beyond what classical optics can achieve8–10 and thereby facilitating the application of quantum systems in diverse fields such as metrology,11–14 computer science,15–17 and communication.18–21 

However, due to the typically low emission rate and detection inefficiency, experimental observation on quantum emitters often necessitates a lengthy experiment time and produces large datasets with a high level of noise,22–24 which renders the fitting process computationally expensive.25,26 The identification of multiphoton states is commonly addressed by assembling multiple single-photon detectors that are time-correlated.27–31 This exacerbates the challenge of data analysis, as traditional computational techniques such as the Levenberg–Marquardt (L-M) method32 require extensive co-detection events to achieve a satisfactory accuracy.26 

Over the past decade, novel data-driven formalisms such as Machine Learning (ML) have introduced new possibilities in quantum photonics experiments.33–37 Specialized in analyzing large and sparse datasets, ML models have provided speedup by orders of magnitude in certain quantum measurements38,39 and show potential in overcoming the inherent limitations of conventional fitting methods particularly in the low-photon flux regime.40 For example, a Convolutional Neural Network (CNN)-based algorithm was developed for rapid classification of single photon emitters in the NV center of nanodiamonds.41 Compared to the L-M method, the accuracy is improved with the CNN model by recognizing subtle features extracted from sparse correlation data. A single artificial neuron model was developed to reduce the required average number of photons down to less than one for distinguishing thermal light from coherent light in low-light measurements.42 A study by Cortes et al.43 demonstrated that employing statistical learning methods for the reconstruction of g(2) data can substantially accelerate the data acquisition process from few-shot measurements.

While considerable efforts were directed toward single-photon emitters,44–47 the emission of multiple, indistinguishable photons also becomes favorable for quantum systems,48,49 making them promising candidates to exert further influence on various quantum applications such as Boson sampling.50 As the commonly adopted g(2) correlation proves inadequate when detecting the photon “superbunching” in higher Fock states, it necessitates the introduction of higher-order correlation.49,51,52

In this study, we present a 2D CNN based ML model for rapid classification of multiphoton states, including photon Fock states up to |3⟩ and coherent states of laser emission. The time-dependent photon detection data are simulated, and by mixing each Fock state with the corresponding coherent state, the quantum purity of emitters can be manipulated. g(3) correlation is performed on the simulation data and fitted using a supervised machine learning model to return the photon state classification results. Through model training and optimization, the average accuracy of classification surpasses 90% for all the Fock states, with an overall accuracy of 94%. The model exhibits an effective performance with sparse datasets, with only 800 photon detection events to achieve a 90% average accuracy. Finally, we propose an experimental setup for quasi-real-time photon state classification accelerated by the ML model. For the first time, a 2D CNN algorithm is employed for identifying multiphoton states and shows enhanced accuracy and data efficiency.

To simulate photon correlation experiment within an extended Hanbury Brown and Twiss (HBT) scheme53 shown in Fig. 7, the Monte Carlo method is used to generate photon detection events with the arrival timestamps. The simulation model consists of three primary parts: photon stream emission, transport, and detection as the output. This model is built upon the TensorFlow Probability (TFP) extension and the TensorFlow Distribution (TFD) probabilistic model.54 

To emulate various quantum pureness of emitters and experimental factors that cause deviation on g(k)(0), we incorporate each Fock state |n⟩ with a corresponding coherent state |α⟩, with both states sharing the same average photon number. The expectation of photon number for a coherent state is given by α|n̂|α=|α|2, hence α=n̄.55 A pure coherent state with a temporally random sequence of photons produces the correlation g(k)(0) to always be 1, implying no correlation between the detection events.6 When mixing a Fock state and the corresponding coherent state, the variance of the photon number is increased, while the average photon number remains the same. Hence, the photon distribution of such a mixture lies between a delta function characterizing the Fock state and a Poissonian distribution for the coherent state. In addition, the g(k)(0) values of the superposition state are calculated as follows:
(1a)
(1b)
where n means the photon Fock state, p represents the probability of the Fock state in the mixture, and (1 − p) denotes the portion of the coherent state. By manipulating the contribution of each state, the lower or upper bound of g(k)(0) can be reached, representing an ideal quantum emitter at state |n⟩ or a coherent source at state |α⟩ with α=n, respectively.

Figure 1 illustrates the simulation of an imperfect |2⟩ state emitter as an example. To simulate the photon emission, the TFD categorical function is used to generate a list of light source labels, shown by the top row in Fig. 1(a). Each label represents either a quantum emission in the |n⟩ state (red) or a coherent laser emission |α⟩ (blue) with an average photon number of n. While for an ideal quantum emitter, no laser labels will be included, non-ideal quantum emitters contain a mixture of both quantum emission labels and laser labels that are randomly distributed within the list. The quality of quantum emitters can be assessed by the portion of quantum emission labels in the list, controlled by a simulation parameter called “quantum light probability” (QLP) that ranges from 0 to 1. In practical experiments, the significance of QLP reflects a cumulative effect from various factors such as the background signal of classical light, which impact the correlation result g(k)(0). As shown in Fig. 1, the |2⟩ state emitter with a QLP of 0.5 results in an even distribution of half quantum emission labels and half laser labels. Each light label is then replaced by an integer to represent the number of photons in each emission, illustrated by the second row in Fig. 1(a). For quantum emission labels, it is straightforward to fill in an integer 2 (in red) to represent the emission of two identical photons from the |2⟩ emitter. The laser labels are substituted by using the TFD Poisson function, which generates integers (in blue) in a Poissonian distribution with the average value of 2. This adjustment comes from the photon number for each emission from a coherent laser that follows a Poissonian distribution shown by the blue plot on the right of Fig. 1(a), instead of being a constant for quantum emission represented by the red plot.

FIG. 1.

Simulation of the HBT experiment, g(3) correlation algorithm, and rapid Fock state classification based on machine learning. (a) Schematic view of the Monte Carlo simulation. The categorical function produces a list of label markers distinguishing between quantum light and laser, where in this instance, “n” represents a photon at Fock state |2⟩, while “α” means a coherent state |α⟩ with the average number of photons set to 2. The portion of quantum light labels and laser labels in the list is determined by the quantum light probability, set here at 50%, indicating an equal mixture of both. Each marker is then replaced by the actual photon number on the second row, depending on its distribution: each quantum light “n” is replaced by 2, from its delta distribution shown by the red plot on the right, while each “α” is replaced with an integer generated by a Poissonian distribution function with an average of 2, as illustrated by the blue plot. (b) g(3) correlation algorithm. In the purple rectangle, the second row with photon numbers is equally split into three rows using a multinomial function, representing photon events detected by three virtual detectors labeled as “d1,” “d2,” and “d3.” The yellow lines linking detection events from different detectors illustrate the algorithm of g(3) correlation associated with a particular data point marked in purple on the g(3) heat map to the right. Here, τ12 and τ13, representing the time differences between detectors 1 and 2 and detectors 1 and 3 are −1 and 2, respectively. (c) CNN-based Fock state classification. The g(3) correlation results are fed into a CNN model that is pre-trained with similar simulated data, and the model yields the classification result.

FIG. 1.

Simulation of the HBT experiment, g(3) correlation algorithm, and rapid Fock state classification based on machine learning. (a) Schematic view of the Monte Carlo simulation. The categorical function produces a list of label markers distinguishing between quantum light and laser, where in this instance, “n” represents a photon at Fock state |2⟩, while “α” means a coherent state |α⟩ with the average number of photons set to 2. The portion of quantum light labels and laser labels in the list is determined by the quantum light probability, set here at 50%, indicating an equal mixture of both. Each marker is then replaced by the actual photon number on the second row, depending on its distribution: each quantum light “n” is replaced by 2, from its delta distribution shown by the red plot on the right, while each “α” is replaced with an integer generated by a Poissonian distribution function with an average of 2, as illustrated by the blue plot. (b) g(3) correlation algorithm. In the purple rectangle, the second row with photon numbers is equally split into three rows using a multinomial function, representing photon events detected by three virtual detectors labeled as “d1,” “d2,” and “d3.” The yellow lines linking detection events from different detectors illustrate the algorithm of g(3) correlation associated with a particular data point marked in purple on the g(3) heat map to the right. Here, τ12 and τ13, representing the time differences between detectors 1 and 2 and detectors 1 and 3 are −1 and 2, respectively. (c) CNN-based Fock state classification. The g(3) correlation results are fed into a CNN model that is pre-trained with similar simulated data, and the model yields the classification result.

Close modal
With the substitution of photon numbers, the emulated photon stream is fed into a virtual multi-port beam splitter (BS) using the TFD multinomial function. Photons in each emission packet are split into three paths with a multinomial distribution, as shown in the purple rectangle in Fig. 1(b). The split photon events are recorded by virtual photon detectors [noted by “d1,” “d2,” and “d3” in Fig. 1(b)] and output as simulation results, including the number of photons detected each time, along with the corresponding timestamp. For convenience and consistency, we assume that all light emitters are periodically excited by a pulsed source, with the pulse period serving as the unit of time throughout this simulation. For simplicity, the photon number of each detection can be resolved by the virtual detectors. Experimental imperfections such as optical loss, dark counts, and detection inefficiency are not directly simulated, as these imperfections can be addressed either through post-selections with respect to the triggering signal or through the use of low-loss photon detectors,56,57 without significantly altering the overall photon distribution. The total number of detection events is another parameter that represents the experiment time and determines the data sparsity of g(3) correlation,
(2)

The detected photon events are analyzed using the g(3) function, by Eq. (2). The time differences τ12 and τ13 between detectors 1 and 2 and detectors 1 and 3, respectively, are variables to compute the normalized g(3)(τ12, τ13) correlation, shown by the heat map on the right of Fig. 1(b). The data point at (τ12 = −1, τ13 = 2) marked in purple on the heat map is chosen to demonstrate the g(3) algorithm depicted with the yellow stripes in Fig. 1(b). With these specific time differences, the d2 vector is shifted forward by 1 unit relative to d1, while the d3 vector is shifted backward by 2 units. The dot product of the three vectors is divided by the vector length and normalized by its mean, shown by the purple-marked pixel in Fig. 1(b). The g(3) function implements the above algorithm for all τ12 and τ13 variables, ranging from −16 to 16 with an interval of 1, and returns a two-dimensional matrix of g(3) correlation results. A 2D CNN model is developed for photon state classification based on the correlation results, shown in Fig. 1(c). In the given example, the |2⟩ state quantum emitter is successfully identified among four options: laser within a coherent state and quantum emitters in |1⟩, |2⟩, and |3⟩ states, respectively.

According to Eq. (1), to characterize a multiphoton Fock state |n⟩, a k-th order correlation with k > n is required to ensure that the central point g(n)(0) remains zero. Although the presented algorithm is extensible for high-order correlations, due to the rapid scaling of computational complexity with the order of correlation,26  g(3) is chosen for categorizing and characterizing multiphoton states within a manageable computation time. Taking the time difference τ as the variable for g(2), the g(3) function has two variables τ12 and τ13 represented by x and y axes, being the time differences between detectors 1 and 2 and detectors 1 and 3, respectively. The third time difference τ23, which can be recursively derived from the other two as τ23 = τ13τ12, is not considered as a variable in this context.

Table I provides theoretical g(2)(0) and g(3)(0) values for Fock states |1⟩, |2⟩, and |3⟩ when mixed with corresponding coherent states |α⟩. The quantum light probability (QLP) represents the portion of Fock states in the mixture, e.g., QLP = 1 denotes an ideal quantum state, while QLP = 0 denotes a coherent state. During the simulation, coherent states with an average photon number of 3 or less are modeled, which is often achieved by strongly attenuating a classical source, such as a laser, when measuring quantum emission.58 The g(3) results in Fig. 2 align well with the g|2(3) values in Table I: First, the zero-delay value g|2(3)(0) is 0, as the order of correlation 3 is greater than the value of the Fock state 2. Second, the data within the cross-pattern and the anti-diagonal that are in green approach 0.5 due to the special values of time differences. For the vertical (or horizontal) stroke of the cross, the time difference τ12 (or τ13) is 0. Meanwhile, for the anti-diagonal (bottom left to top right), the two time differences are equal, signifying that τ13 is 0. With one time difference being zero, the g(3) correlation can be treated to be equivalent to a g(2), which correlates only two detectors with one variable of time difference. Hence, the value of 0.5 on the heatmap can be explained with g|2(2)(0)=0.5. Finally, the remaining data points are normalized with respect to the average photon number, typically resulting in a normalization factor of 1. To better visualize the 2D g(3) matrix, three cross sections marked with colored arrowheads on the heatmap are plotted bar charts on the right side of Fig. 2. Each represents the evolution of g(3) as a function of only τ12, with τ13 being set at special values: τ13 = −τ12 for black bars, τ13 = 0 for purple, and τ13 = τ12 for yellow, corresponding to the main diagonal, horizontal stroke in the cross-pattern, and the anti-diagonal, on the heatmap.

TABLE I.

Theoretical g(2) and g(3) values of Fock states when mixing with the corresponding coherent states. The quantum light probability indicates the portion of quantum light in the simulation when it is mixed with a corresponding coherent state. A probability of 0.5 (in green) signifies a light source composed of equal parts of quantum light and laser, which is set to be the critical point distinguishing between a quantum emitter and a coherent laser. Values below 0.5 (in blue) or above 0.5 (in red) are labeled as a laser or quantum light, respectively.

 
 
FIG. 2.

An example of g(3) correlation results from simulated data for a |2⟩ Fock state without mixing with a coherent state. Left: Normalized g(3) correlation in a heat map, as a function of τ12 and τ13, which represent the time differences between detectors 1 and 2 and detectors 1 and 3, respectively. Right: 2D bar plots of g(3) correlation, as a function of only τ12, while τ13 is set to be either negative τ12, 0, or equal to τ12, corresponding to black, magenta, and yellow bars. The trace of each bar plot is indicated by the arrows within the same color on the heat map.

FIG. 2.

An example of g(3) correlation results from simulated data for a |2⟩ Fock state without mixing with a coherent state. Left: Normalized g(3) correlation in a heat map, as a function of τ12 and τ13, which represent the time differences between detectors 1 and 2 and detectors 1 and 3, respectively. Right: 2D bar plots of g(3) correlation, as a function of only τ12, while τ13 is set to be either negative τ12, 0, or equal to τ12, corresponding to black, magenta, and yellow bars. The trace of each bar plot is indicated by the arrows within the same color on the heat map.

Close modal

In contrast to the second-order correlation, which has only one critical point g(2) (0) for categorization, the third-order correlation g(3) features multiple critical points: one at g(3)(0) and multiple at g(3)(0, τ ≠ 0), effectively serving as g(2)(0). In accordance with Eq. (1), the complete g(3) cross correlation enables the identification of unknown Fock state emission and the assessment of its quantum purity. Even for Fock states higher than |2⟩, although g(3) (0) is no longer 0, it can still be determined from the g(3) (0) and g(2)(0) values. Meanwhile, the dynamics of excitons can be further explored through other correlation algorithms, such as two-time second-order auto-correlation.59 In addition, the two-dimensional nature of the g(3) data proves advantageous for CNN models that are adept at extracting spatial features from images.60,61

While the g(3) critical values can typically be determined using normal fitting methods, the experimental data often exhibit significant sparsity due to the limited number of photon detection events, entailing the introduction of advanced data fitting techniques. To efficiently categorize photon states by fitting the g(3) data, a machine learning model is developed based on the open-source API Keras.62 The model architecture and layer operations are illustrated in Fig. 3. The model mainly comprises 2D convolution layers (Conv2D), 2D max pooling layers (MaxPool), 2D global average pooling layers (AvePool), and dense layers. The output shape of each layer is noted in parentheses.

FIG. 3.

Schematics of the CNN model. Conv2D is for the 2D convolution layer, MaxPool is for the 2D max pooling layer, AvePool is for the 2D global average pooling layer, and ×2 represents a layer or sequence of operations (in brackets) that is repeated twice. The shape of output data from each layer is written in parentheses. The dashed lines illustrate the connection between input data points and respective output data during layer operations. The initial input data consist of g3 correlation results in a 2D matrix, shown by the white cells, while red cells represent zero paddings. The Conv2D layer produces 3D matrices of data with a depth of 64, equal to the number of convolutional kernels (blue dashed lines represent one of the kernels). The data shape is condensed to one dimension at the AvePool layer, and the final dense layer outputs a score vector, indicating the confidence of predictions for each light source category.

FIG. 3.

Schematics of the CNN model. Conv2D is for the 2D convolution layer, MaxPool is for the 2D max pooling layer, AvePool is for the 2D global average pooling layer, and ×2 represents a layer or sequence of operations (in brackets) that is repeated twice. The shape of output data from each layer is written in parentheses. The dashed lines illustrate the connection between input data points and respective output data during layer operations. The initial input data consist of g3 correlation results in a 2D matrix, shown by the white cells, while red cells represent zero paddings. The Conv2D layer produces 3D matrices of data with a depth of 64, equal to the number of convolutional kernels (blue dashed lines represent one of the kernels). The data shape is condensed to one dimension at the AvePool layer, and the final dense layer outputs a score vector, indicating the confidence of predictions for each light source category.

Close modal

Before passing to any layers, the correlation results are truncated and rescaled to improve the model’s performance. The first 32 elements in each dimension of the original 33 × 33 g(3) matrix are retained so that the truncated matrix, with a shape of 32 × 32, aligns with the preference of CNN layers for input data dimensions that are multiples of 2. In cases of simulation with few-shot data, where normalized g(3) results may contain extremely high values, the rescaling ensures that all the g(3) elements fall within the range of 0–1. The preprocessed datasets, depicted by the white square cells in the first diagram of Fig. 3, have a spatial dimension of 32 × 32 and a depth of 1.

As the core of the CNN model, a 2D convolutional layer connecting to the input layer extracts the spatial features. The input data are zero-padded on the outer border, as indicated by the red square cells in Fig. 3, to prevent rapid degradation of the information at original borders and maintain the output dimension of the Conv2D layer at 32 × 32. Each Conv2D layer contains 64 convolutional kernels, and each kernel is a 3 × 3 matrix that spatially slides across the data array with a step of 1. The blue dashed boxes in Fig. 3 illustrate the sum product calculation of a convolutional kernel, with the dashed lines connecting to the cells representing the results. The above operation is performed by all the 64 kernels, with each kernel capturing a unique feature from the input, producing the Conv2D output with a depth of 64. The zero-padding is also applied to the Conv2D output data to maintain the layer dimension. Two of Conv2D layers are sequentially stacked and followed by a 2D max pooling layer, which calculates the maximum values of each pooling patch to spatially downsample the data dimension. Three patches are depicted by the red dashed boxes in Fig. 3 for visualization, and this process is repeated for all the 64 slices in depth. The layer dimension is halved with a step size of 2 for max pooling, resulting in an output shape of 16 × 16 × 64. The number of trainable parameters is significantly reduced by the downsampling so that a relatively low model capacity is maintained without losing essential information.63 The architecture enclosed in the square bracket in Fig. 3, composed of two consecutive convolutional layers and one MaxPooling layer, is executed twice for learning hierarchical features.

The following global average pooling layer computes the average for each slice in depth, shown by the gray cells in the figure. Each average reflects the importance of an extracted feature from the neural network. This array of feature scores is analyzed by two dense layers, with each neuron fully connected to all the neurons from the previous layer. The output layer computes scores for four potential photon states: the coherent state, as well as Fock states |1⟩, |2⟩, and |3⟩. Each score indicates the model’s confidence in predicting a specific photon state, and the state with the highest score is determined as the outcome. As the QLP represents the deviation of correlation results from the theoretical values of an ideal |n⟩ quantum emitter due to experimental factors, the exact value is not included in the machine learning outcome. Instead, by discretizing QLPs into bins, the regression task is reformulated into a classification problem, which aligns better with the CNN’s specialty.64 

Among the wide variety of machine learning methods, algorithms with straightforward architectures, such as k-Nearest Neighbors65 (k-NNs), Logistic Regression66 (LG), and Gradient Boosting67 (GB), have been implemented for classification problems68–70 and have also been applied to quantum and photonics research.71,72 While these simple classifiers offer better interpretability73 and lower computational requirements74 compared to the proposed 2D-CNN model, they are less suitable and may have limited performance for photon correlation problems. As shown in the supporting information of Ref. 41, for a binary decision problem of single-photon emission based on 1D g(2) correlation, the CNN model demonstrates the highest overall fitting accuracy, outperforming other algorithms including GB, and a voting classifier that combines LG and k-NN. In addition, the accuracy of the CNN improves with more training data, while neither LR nor k-NN shows further noticeable improvements.

In order to enhance the experimental determination of photon states, it is essential to achieve fast and accurate predictions from sparse and noisy experimental data, while the cost of the model training process is not the main concern. One drawback of simple ML algorithms such as k-NN is that prediction with large training datasets becomes computationally expensive due to its iterative comparison process.75 Algorithms such as LG and GB require upfront feature engineering to adapt to the non-linear relationship between g(3) results and photon Fock states. This process is particularly challenging for real-time fitting with noisy correlation data due to signal obscuration and reduced feature quality. In contrast, the complexity of neural networks for decision-making becomes a trade-off for efficiently capturing spatial hierarchies and patterns, especially with sufficient training data.74,76,77 By adopting architectures of well-developed models and utilizing novel hyperparameter optimization frameworks such as KerasTuner,78 CNN becomes a feasible and efficient ML technique for physicists.

This CNN model consists of 11 layers (excluding the input layer) and contains over 150 000 trainable parameters. Except for the dense layers, batch normalization is applied to the output of each layer (not shown in the diagram). The Adam method is chosen as the model optimizer, and all CNN layers use ReLU (Rectified Linear Unit) as the activation function. The g(3) data are simulated for |1⟩, |2⟩, and |3⟩ Fock states, with two variable parameters: QLP and the number of detection events. The QLP has 21 possible values ranging from 0 to 1 with intervals of 0.05, to represent the quantum pureness of the emitter. Qualified quantum emitters require a QLP greater than or equal to 0.5, while a QLP less than 0.5 will be categorized as a coherent state. The number of detection events that simulates the experiment time ranges from 100 to 100 000. Values from 100 to 10 000 are taken at intervals of 100, while values from 10 000 to 100 000 are taken at intervals of 1000, totaling 190 values. A total of 11 970 cases are simulated, with 100 measurements for each case, resulting in 11 970 000 g(3) correlation datasets. 70% of the datasets are used for model training, with a batch size of 32 and a single epoch, given the substantial data that do not necessitate iterations. 20% of the data are validation sets for hyperparameter tuning, while the remaining 10% are used for testing the accuracy. Considering this study’s emphasis on improving fitting accuracy for sparse data, the model’s evaluation primarily relies on the required number of photon detection events to achieve satisfactory predictions, instead of computational performance metrics such as computation time.

The average accuracy of classification results is shown in Fig. 4, as a function of variables, including the photon state, QLP, and the number of detection events. The left bar chart displays the average accuracy for all cases in each photon state, with the coherent state plotted in blue, and the Fock states |1⟩, |2⟩, and |3⟩ are plotted in red. An accuracy of 90% is achieved for all photon states, with classification of the coherent state and Fock state |3⟩ being higher compared to |1⟩ and |2⟩. The middle and right plots describe the accuracy over all the photon states as functions of QLP and the number of events. A significant drop in accuracy is observed when the QLP approaches 0.5. This is because QLP = 0.5 is defined as the threshold distinguishing between quantum light (Fock state) labels and non-quantum light (coherent state) labels. Near this decision boundary, the similarity in the correlation results posts a challenge for accurate classification, which is further indicated in Fig. 5. As shown on the right, the average accuracy across all light categories improves from 72% to 95% as the number of simulation events increases from 100 to 10 000. The accuracy boundary of 90%, achieved with only 800 events, is depicted by the green vertical line.

FIG. 4.

Accuracy of the CNN-based photon state classifier. Left: Overall accuracy for each light source category. Blue: Laser with a coherent state. Red: Quantum light Fock states. Middle: Averaged accuracy over all light source categories as a function of quantum light probability, which indicates the portion of quantum light in the simulation when it is mixed with a corresponding coherent state. A probability of 0.5 signifies a light source composed of equal parts of quantum light and laser, while 0.0 or 1.0 indicates a pure laser or pure quantum light, respectively. Right: Averaged accuracy over all light source categories as a function of number of detection events. The number of detections is in logarithmic scale for better visualization. The green vertical line indicates the 90% accuracy boundary.

FIG. 4.

Accuracy of the CNN-based photon state classifier. Left: Overall accuracy for each light source category. Blue: Laser with a coherent state. Red: Quantum light Fock states. Middle: Averaged accuracy over all light source categories as a function of quantum light probability, which indicates the portion of quantum light in the simulation when it is mixed with a corresponding coherent state. A probability of 0.5 signifies a light source composed of equal parts of quantum light and laser, while 0.0 or 1.0 indicates a pure laser or pure quantum light, respectively. Right: Averaged accuracy over all light source categories as a function of number of detection events. The number of detections is in logarithmic scale for better visualization. The green vertical line indicates the 90% accuracy boundary.

Close modal
FIG. 5.

Accuracy distribution of photon Fock state classification, as functions of quantum light probability and the number of detections for each Fock state case. From top to bottom: Accuracy distributions for light sources, comprising |3⟩, |2⟩, or |1⟩ Fock states of quantum light, respectively, each combined with its corresponding coherent state |α⟩, where the average photon number is 3, 2, or 1, respectively. For each case, a quantum light probability of 0.5 signifies a light source composed of equal parts of quantum light and laser, while 0.0 or 1.0 indicates a pure laser or pure quantum light, respectively. Green vertical line: boundary of 90% averaged accuracy over all light source categories, at 800 detection events.

FIG. 5.

Accuracy distribution of photon Fock state classification, as functions of quantum light probability and the number of detections for each Fock state case. From top to bottom: Accuracy distributions for light sources, comprising |3⟩, |2⟩, or |1⟩ Fock states of quantum light, respectively, each combined with its corresponding coherent state |α⟩, where the average photon number is 3, 2, or 1, respectively. For each case, a quantum light probability of 0.5 signifies a light source composed of equal parts of quantum light and laser, while 0.0 or 1.0 indicates a pure laser or pure quantum light, respectively. Green vertical line: boundary of 90% averaged accuracy over all light source categories, at 800 detection events.

Close modal

To explicitly illustrate the model’s performance on a case-by-case basis, Fig. 5 presents accuracy heatmap for the three Fock states, where each pixel is colored according to the average accuracy of all cases with a specific QLP and the number of events. The accuracy of coherent state classification is included in each Fock state, where QLP is less than 0.5, instead of a separate plot. In general, the accuracy is significantly enhanced with an increasing number of events across most QLPs. A consistent accuracy drop occurs at the decision boundary with QLP approaching 0.5, which aligns with the trend in the middle figure of Fig. 4. For all three cases, the performance decreases at the left of the green line (the 90% overall accuracy boundary from Fig. 4, right), where datasets contain fewer detection events and have QLPs higher than 0.5. While a relatively low accuracy seems reasonable with few-shot data, the “vertical” asymmetry of accuracy comes from the large data sparsity and uncertainty with higher QLPs. In photon statistics, coherent states with a larger variance are more likely to emit higher photon states,6 which produce recognizable features in the correlation. For higher QLPs, the occurrence of coherent states drops down, consequently leading to less distinctive features in the correlation results for recognition. As the number of simulation events increases, the g(3) results tend to stabilize for both coherent states and Fock states, which leads to a quick disappearance of the triangular region with low accuracy.

For the |1⟩ case, the low-accuracy triangle vanishes earlier than |2⟩ and |3⟩ with growing detection events. According to Table I, the |1⟩ state exhibits larger variations in g(2)(0) and g(3)(0) values with changing QLP, which contributes greatly to the CNN classification. Conversely, the smaller differences for |2⟩ and |3⟩ make them less conducive to recognition with fewer events, resulting in lower rate of increase in accuracy. At the decision boundary with a QLP of 0.5, the accuracy of the |3⟩ state exhibits a higher rate of increase compared to |1⟩. This is because the relatively higher number of photons from the |3⟩ state facilitates pattern recognition by the model and consequently requires fewer events for accurate classification.

For the |2⟩ state, the accuracy shows a decrease with more detection events, particularly at intermediate QLPs. Some mispredictions of coherent states can be attributed to the confusion near the decision boundary; as mentioned earlier, the majority of mispredictions involve the |3⟩ state, potentially caused by overfitting.79,80 Given that |2⟩ and |3⟩ exhibit similar values of g(2)(0), the decision boundary g(2)(0) = 0.83 [determined by Eq. (1)] for distinguishing |3⟩ and its coherent state is misused for identifying |2⟩. Consequently, for the |2⟩ cases with g(2)(0) greater than 0.83 (or the QLP exceeds 0.3), misclassification as |3⟩ results in a noticeable drop in accuracy, even with increased numbers of detection events. While they possess different g(3)(0) values, the contribution of this single element to the classifier is overshadowed by the abundance of g(2)(0) elements in the correlation matrix. Given that typically low photon co-detection events are involved in quantum correlation experiments,58 this paper primarily focuses on datasets containing fewer than 10 000 detection events.

To visualize the prediction distribution over photon states, a 3D bar chart of the confusion matrix is shown in Fig. 6. The diagonal elements in green indicate the average accuracy for each state, while the off-diagonal elements in yellow describe the distribution of misclassification. The highest accuracy achieved is 98.7% for the coherent state, while the accuracy for Fock states remains above 90%. This is attributed to more datasets that are simulated for coherent states, which boosts the ML algorithm during the training process. The model showcases a balanced performance in recognizing the Fock states, with the most common incorrect prediction being the coherent state, owing to their similarity near the decision boundary. Another notable misclassification is observed between |2⟩ and |3⟩, which arises from the minor difference in the g(3)(0) values as discussed earlier. Other elements in the error matrix are below 0.1% and can be considered negligible.

FIG. 6.

Confusion matrix of the classifier in a 3D bar plot. The correct predictions are represented by the diagonal elements in green, while the off-diagonal elements, indicating incorrect predictions, are depicted in yellow. Blue labels: Laser with a coherent state. Red labels: Quantum light Fock states.

FIG. 6.

Confusion matrix of the classifier in a 3D bar plot. The correct predictions are represented by the diagonal elements in green, while the off-diagonal elements, indicating incorrect predictions, are depicted in yellow. Blue labels: Laser with a coherent state. Red labels: Quantum light Fock states.

Close modal

In contrast to the fitting methods that can be implemented directly, the machine learning algorithm requires an additional training phase prior to data analysis. However, this trade-off can yield substantial enhancement of fitting accuracy and efficiency.39 By achieving an average accuracy of 90% with only 800 simulated events, the ML model demonstrates the capability for reliable classification without requiring lengthy measurement time. Another advantage of ML algorithms is their ability to be optimized for recognizing specific data pattern.81 As 2D CNNs have been extensively developed for recognizing spatial features,61 the integration of g(3) with 2D CNN holds great potential for a novel quantum light classifier, especially in the multiphoton regime.

With Fig. 7, we propose the implementation of the 2D CNN model in quasi-real-time photon state distribution measurements. Examples of different light emitters are included in Fig. 7(a). The model on the top shows single-photon emission from the exciton recombination in a two-level system, such as a semiconductor quantum dot. The red dots represent the single photons emitted periodically in the temporal regime, with the photon Fock state labeled as |1⟩. The middle section illustrates a two-photon emitter that emits two indistinguishable photons represented by the red dots in pairs upon each excitation, with a Fock state of |2⟩. The photon stream from a continuous-wave laser is shown at the bottom, where each emission occurs randomly in time, exhibiting a coherent photon state |α⟩. These listed light sources are coupled into an optical fiber system with an HBT configuration, indicated by the dark blue wavy lines in Fig. 7(b). The incident beam path is evenly split into four output paths by the 50:50 beam splitters (BS). Each output is attached to an avalanche photodiode (APD) that records the photon arrival time and transmits the data to the correlation electronics (CB). g(3) correlation data are computed online and fed into a Tensor Processing Unit (TPU), serving as a machine learning hardware accelerator. The ML model is pre-trained and optimized to provide real-time photon state classification results, as shown in Fig. 7(c). With the specialty of rapid fitting on sparse datasets, the measurement time on individual samples can be substantially reduced, and immediate feedback from the ML analysis enables real-time parameter optimization.

FIG. 7.

Experiment layout for measuring quasi-real-time Fock state distribution enabled by machine learning. (a) Light sources including quantum light emitters, specifically two-level systems in a semiconductor quantum dot, with potential photon Fock states such as |1⟩ or |2⟩ (red dots), and an attenuated laser in a coherent state |α⟩ (blue dots) with a Poissonian photon distribution. (b) Optical setup. The photon emission is coupled into an optical fiber and is then equally split into four paths by the 50:50 fiber Beam Splitter (BS). Each path leads to an Avalanche Photodiode (APD) that is connected to the Correlation Board (CB). (c) Correlation data will be analyzed by an optimized CNN model in the Tensor Processing Unit (TPU), serving as a specialized hardware accelerator, which returns the predicted light source category.

FIG. 7.

Experiment layout for measuring quasi-real-time Fock state distribution enabled by machine learning. (a) Light sources including quantum light emitters, specifically two-level systems in a semiconductor quantum dot, with potential photon Fock states such as |1⟩ or |2⟩ (red dots), and an attenuated laser in a coherent state |α⟩ (blue dots) with a Poissonian photon distribution. (b) Optical setup. The photon emission is coupled into an optical fiber and is then equally split into four paths by the 50:50 fiber Beam Splitter (BS). Each path leads to an Avalanche Photodiode (APD) that is connected to the Correlation Board (CB). (c) Correlation data will be analyzed by an optimized CNN model in the Tensor Processing Unit (TPU), serving as a specialized hardware accelerator, which returns the predicted light source category.

Close modal

With the development of ML software,54,62,82 further enhancement on the presented prototype becomes possible. For example, 2D locally connected layers can be a promising substitute for the CNN layers for capturing the pattern (shown on Fig. 2 left) that are spatially fixed in correlation data. In these layers, the parameter-sharing scheme between convolutional kernels is relaxed, which facilitates the recognition of specific structures.83 However, the model’s capacity will be increased by incorporating multiple such layers. Moreover, an ensemble of similar models can be developed, with each model trained for specific cases, to enhance the overall performance.

The presented simulation model can be tailored to simulate broader scenarios involving random coherent lights with low average photon numbers mixing with quantum emission. This adaptation is particularly relevant as quantum light detection often involves laser emission, as the excitation source of quantum emitters. In addition, this model can be further modified to compute dynamic correlations, such as the conditional auto-correlation function (CACF) for studying heralded photons.84 As the landscape of quantum emitters has been enriched by the emerging 2D materials with functional heterostructures and satisfactory fabrication efficiency such as graphene,85 TMDs (Transition Metal Dichalcogenides),47,86,87 and moiré superlattices,88 the scalability for photonic circuits is often impeded by the rare and random occurrence of quantum emitters,47 which are considered from defects or strains in these atomically thin layers.89,90 This study offers insights into Fock state mapping on 2D materials, rapidly determining the spatial location and assessing the quantum purity of potential quantum emitters.

Recently, novel methods for multiphoton preparation have been demonstrated. Reference 91 proposed a theoretical model of deterministic generation of large Fock states up to |100⟩ based on the resonant interaction of a coherent state with two-level systems. The experimental generation of eight indistinguishable photons using temporal-to-spatial demultiplexing with single photon sources was presented.92 The proposed g(3) correlation–2D CNN combination introduces new possibilities for categorizing and characterizing multiphoton states.

In addition, this approach opens up the possibility for optimizing photosynthesis through real-time quantum feedback. Reference 93 has shown that the implantation of biocompatible quantum dots (BQDs) enables quantum measurements for monitoring intracellular environments, effectively overcoming the persistent challenge of strong autofluorescence background in traditional spectroscopy of plant cells.94 

In this article, we present a machine learning model for quasi-real-time categorization of photon states with g(3) correlation and propose the implementation in experiments. This methodology introduces new feasible solutions for identifying quantum emitters and holds broad applications in quantum metrology within the field of nanophotonics.

The authors have no conflicts to disclose.

Guangpeng Xu: Data curation (lead); Formal analysis (lead); Software (lead); Validation (lead); Visualization (lead); Writing – original draft (lead). Jeffrey Carvalho: Writing – original draft (supporting). Chiran Wijesundara: Writing – original draft (supporting). Tim Thomay: Conceptualization (lead); Funding acquisition (lead); Project administration (lead); Software (supporting); Supervision (lead); Visualization (supporting); Writing – original draft (supporting).

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1.
W.
Zhang
,
D.-S.
Ding
,
Y.-B.
Sheng
,
L.
Zhou
,
B.-S.
Shi
, and
G.-C.
Guo
, “
Quantum secure direct communication with quantum memory
,”
Phys. Rev. Lett.
118
,
220501
(
2017
).
2.
S.-J.
Yang
,
X.-J.
Wang
,
X.-H.
Bao
, and
J.-W.
Pan
, “
An efficient quantum light–matter interface with sub-second lifetime
,”
Nat. Photonics
10
,
381
384
(
2016
).
3.
S.
Mittal
,
E. A.
Goldschmidt
, and
M.
Hafezi
, “
A topological source of quantum light
,”
Nature
561
,
502
506
(
2018
).
4.
X.
He
,
H.
Htoon
,
S. K.
Doorn
,
W. H. P.
Pernice
,
F.
Pyatkov
,
R.
Krupke
,
A.
Jeantet
,
Y.
Chassagneux
, and
C.
Voisin
, “
Carbon nanotubes as emerging quantum-light sources
,”
Nat. Mater.
17
,
663
670
(
2018
).
5.
J.-Y.
Hu
,
B.
Yu
,
M.-Y.
Jing
,
L.-T.
Xiao
,
S.-T.
Jia
,
G.-Q.
Qin
, and
G.-L.
Long
, “
Experimental quantum secure direct communication with single photons
,”
Light: Sci. Appl.
5
,
e16144
(
2016
).
6.
A.
Fox
,
Quantum Optics: An Introduction
,
Oxford Master Series in Physics
(
Oxford University Press
,
Oxford
,
2006
).
7.
A.
Einstein
, “
Über einen die Erzeugung und Verwandlung des Lichtes betreffenden heuristischen Gesichtspunkt
,”
Ann. Phys.
322
,
132
148
(
1905
).
8.
J. L.
O’Brien
,
A.
Furusawa
, and
J.
Vučković
, “
Photonic quantum technologies
,”
Nat. Photonics
3
,
687
695
(
2009
).
9.
J.
Wang
,
F.
Sciarrino
,
A.
Laing
, and
M. G.
Thompson
, “
Integrated photonic quantum technologies
,”
Nat. Photonics
14
,
273
284
(
2020
).
10.
P.
Lodahl
,
S.
Mahmoodian
,
S.
Stobbe
,
A.
Rauschenbeutel
,
P.
Schneeweiss
,
J.
Volz
,
H.
Pichler
, and
P.
Zoller
, “
Chiral quantum optics
,”
Nature
541
,
473
480
(
2017
).
11.
J.
Ye
,
H. J.
Kimble
, and
H.
Katori
, “
Quantum state engineering and precision metrology using state-insensitive light traps
,”
Science
320
,
1734
1738
(
2008
).
12.
O.
Pinel
,
J.
Fade
,
D.
Braun
,
P.
Jian
,
N.
Treps
, and
C.
Fabre
, “
Ultimate sensitivity of precision measurements with intense Gaussian quantum light: A multimodal approach
,”
Phys. Rev. A
85
,
010101
(
2012
).
13.
J. P.
Dowling
and
K. P.
Seshadreesan
, “
Quantum optical technologies for metrology, sensing, and imaging
,”
J. Lightwave Technol.
33
,
2359
2370
(
2015
).
14.
E.
Polino
,
M.
Valeri
,
N.
Spagnolo
, and
F.
Sciarrino
, “
Photonic quantum metrology
,”
AVS Quantum Sci.
2
,
024703
(
2020
).
15.
I.
Walmsley
, “
Light in quantum computing and simulation: Perspective
,”
Opt. Quantum
1
,
35
40
(
2023
).
16.
T. D.
Ladd
,
F.
Jelezko
,
R.
Laflamme
,
Y.
Nakamura
,
C.
Monroe
, and
J. L.
O’Brien
, “
Quantum computers
,”
Nature
464
,
45
53
(
2010
).
17.
P.
Walther
,
K. J.
Resch
,
T.
Rudolph
,
E.
Schenck
,
H.
Weinfurter
,
V.
Vedral
,
M.
Aspelmeyer
, and
A.
Zeilinger
, “
Experimental one-way quantum computing
,”
Nature
434
,
169
176
(
2005
).
18.
R.
Ursin
,
F.
Tiefenbacher
,
T.
Schmitt-Manderbach
,
H.
Weier
,
T.
Scheidl
,
M.
Lindenthal
,
B.
Blauensteiner
,
T.
Jennewein
,
J.
Perdigues
,
P.
Trojek
,
B.
Ömer
,
M.
Fürst
,
M.
Meyenburg
,
J.
Rarity
,
Z.
Sodnik
,
C.
Barbieri
,
H.
Weinfurter
, and
A.
Zeilinger
, “
Entanglement-based quantum communication over 144 km
,”
Nat. Phys.
3
,
481
486
(
2007
).
19.
L.-M.
Duan
,
M. D.
Lukin
,
J. I.
Cirac
, and
P.
Zoller
, “
Long-distance quantum communication with atomic ensembles and linear optics
,”
Nature
414
,
413
418
(
2001
).
20.
Z.-S.
Yuan
,
X.-H.
Bao
,
C.-Y.
Lu
,
J.
Zhang
,
C.-Z.
Peng
, and
J.-W.
Pan
, “
Entangled photons and quantum communication
,”
Phys. Rep.
497
,
1
40
(
2010
).
21.
J.-W.
Pan
,
C.
Simon
,
C.
Brukner
, and
A.
Zeilinger
, “
Entanglement purification for quantum communication
,”
Nature
410
,
1067
1070
(
2001
).
22.
A.
Kuhn
,
M.
Hennrich
, and
G.
Rempe
, “
Deterministic single-photon source for distributed quantum networking
,”
Phys. Rev. Lett.
89
,
067901
(
2002
).
23.
M.
Hennrich
,
T.
Legero
,
A.
Kuhn
, and
G.
Rempe
, “
Photon statistics of a non-stationary periodically driven single-photon source
,”
New J. Phys.
6
,
86
(
2004
).
24.
E.
Moreau
,
I.
Robert
,
J. M.
Gérard
,
I.
Abram
,
L.
Manin
, and
V.
Thierry-Mieg
, “
Single-mode solid-state single photon source based on isolated quantum dots in pillar microcavities
,”
Appl. Phys. Lett.
79
,
2865
2867
(
2001
).
25.
Z.
Lin
,
L.
Schweickert
,
S.
Gyger
,
K. D.
Jöns
, and
V.
Zwiller
, “
Efficient and versatile toolbox for analysis of time-tagged measurements
,”
J. Instrum.
16
,
T08016
(
2021
).
26.
Z. A.
Kudyshev
,
V. M.
Shalaev
, and
A.
Boltasseva
, “
Machine learning for integrated quantum photonics
,”
ACS Photonics
8
,
34
46
(
2021
).
27.
D.
Höckel
,
L.
Koch
, and
O.
Benson
, “
Direct measurement of heralded single-photon statistics from a parametric down-conversion source
,”
Phys. Rev. A
83
,
013802
(
2011
).
28.
D.
Elvira
,
X.
Hachair
,
V.
Verma
,
R.
Braive
,
G.
Beaudoin
,
I.
Robert-Philip
,
I.
Sagnes
,
B.
Baek
,
S. W.
Nam
,
E. A.
Dauler
,
I.
Abram
,
M. J.
Stevens
, and
A.
Beveratos
, “
Higher-order photon correlations in pulsed photonic crystal nanolasers
,”
Phys. Rev. A
84
,
061802
(
2011
).
29.
M. J.
Stevens
,
S.
Glancy
,
S. W.
Nam
, and
R. P.
Mirin
, “
Third-order antibunching from an imperfect single-photon source
,”
Opt. Express
22
,
3244
3260
(
2014
).
30.
G. A.
Steudle
,
S.
Schietinger
,
D.
Höckel
,
S. N.
Dorenbos
,
I. E.
Zadeh
,
V.
Zwiller
, and
O.
Benson
, “
Measuring the quantum nature of light with a single source and a single detector
,”
Phys. Rev. A
86
,
053814
(
2012
).
31.
M.
Jönsson
and
G.
Björk
, “
Evaluating the performance of photon-number-resolving detectors
,”
Phys. Rev. A
99
,
043822
(
2019
).
32.
D. W.
Marquardt
, “
An algorithm for least-squares estimation of nonlinear parameters
,”
J. Soc. Ind. Appl. Math.
11
,
431
441
(
1963
).
33.
K.
Yao
,
R.
Unni
, and
Y.
Zheng
, “
Intelligent nanophotonics: Merging photonics and artificial intelligence at the nanoscale
,”
Nanophotonics
8
,
339
366
(
2019
).
34.
J.
Zhou
,
B.
Huang
,
Z.
Yan
, and
J.-C. G.
Bünzli
, “
Emerging role of machine learning in light-matter interaction
,”
Light: Sci. Appl.
8
,
84
(
2019
).
35.
Z. A.
Kudyshev
,
A. V.
Kildishev
,
V. M.
Shalaev
, and
A.
Boltasseva
, “
Machine learning–assisted global optimization of photonic devices
,”
Nanophotonics
10
,
371
383
(
2020
).
36.
S.
Yu
,
F.
Albarrán-Arriagada
,
J. C.
Retamal
,
Y.-T.
Wang
,
W.
Liu
,
Z.-J.
Ke
,
Y.
Meng
,
Z.-P.
Li
,
J.-S.
Tang
,
E.
Solano
,
L.
Lamata
,
C.-F.
Li
, and
G.-C.
Guo
, “
Reconstruction of a photonic qubit state with reinforcement learning
,”
Adv. Quantum Technol.
2
,
1800074
(
2019
).
37.
J.
Wang
,
S.
Paesani
,
R.
Santagati
,
S.
Knauer
,
A. A.
Gentile
,
N.
Wiebe
,
M.
Petruzzella
,
J. L.
O’Brien
,
J. G.
Rarity
,
A.
Laing
, and
M. G.
Thompson
, “
Experimental quantum Hamiltonian learning
,”
Nat. Phys.
13
,
551
555
(
2017
).
38.
R. S.
Hegde
, “
Deep learning: A new tool for photonic nanostructure design
,”
Nanoscale Adv.
2
,
1007
1023
(
2020
).
39.
P.
Freire
,
E.
Manuylovich
,
J. E.
Prilepsky
, and
S. K.
Turitsyn
, “
Artificial neural networks for photonic applications—From algorithms to implementation: Tutorial
,”
Adv. Opt. Photonics
15
,
739
834
(
2023
).
40.
S.-B.
Wu
,
Z.-M.
Li
,
J.
Gao
,
H.
Zhou
,
C.-S.
Wang
, and
X.-M.
Jin
, “
Classification of quantum correlation using deep learning
,”
Opt. Express
31
,
3479
3489
(
2023
).
41.
Z. A.
Kudyshev
,
S. I.
Bogdanov
,
T.
Isacsson
,
A. V.
Kildishev
,
A.
Boltasseva
, and
V. M.
Shalaev
, “
Rapid classification of quantum sources enabled by machine learning
,”
Adv. Quantum Technol.
3
,
2000067
(
2020
).
42.
C.
You
,
M. A.
Quiroz-Juárez
,
A.
Lambert
,
N.
Bhusal
,
C.
Dong
,
A.
Perez-Leija
,
A.
Javaid
,
R. d. J.
León-Montiel
, and
O. S.
Magaña-Loaiza
, “
Identification of light sources using machine learning
,”
Appl. Phys. Rev.
7
,
021404
(
2020
).
43.
C. L.
Cortes
,
S.
Adhikari
,
X.
Ma
, and
S. K.
Gray
, “
Accelerating quantum optics experiments with statistical learning
,”
Appl. Phys. Lett.
116
,
184003
(
2020
).
44.
S.
Khalid
and
F. P.
Laussy
, “
Perfect single-photon sources
,”
Sci. Rep.
14
,
2684
(
2024
).
45.
M. D.
Eisaman
,
J.
Fan
,
A.
Migdall
, and
S. V.
Polyakov
, “
Invited review article: Single-photon sources and detectors
,”
Rev. Sci. Instrum.
82
,
071101
(
2011
).
46.
I.
Aharonovich
,
D.
Englund
, and
M.
Toth
, “
Solid-state single-photon emitters
,”
Nat. Photonics
10
,
631
641
(
2016
).
47.
C.
Palacios-Berraquero
,
D. M.
Kara
,
A. R.-P.
Montblanch
,
M.
Barbone
,
P.
Latawiec
,
D.
Yoon
,
A. K.
Ott
,
M.
Loncar
,
A. C.
Ferrari
, and
M.
Atatüre
, “
Large-scale quantum-emitter arrays in atomically thin semiconductors
,”
Nat. Commun.
8
,
15093
(
2017
).
48.
F.
Dell’Anno
,
S.
De Siena
, and
F.
Illuminati
, “
Multiphoton quantum optics and quantum state engineering
,”
Phys. Rep.
428
,
53
168
(
2006
).
49.
K.
Laiho
,
T.
Dirmeier
,
M.
Schmidt
,
S.
Reitzenstein
, and
C.
Marquardt
, “
Measuring higher-order photon correlations of faint quantum light: A short review
,”
Phys. Lett. A
435
,
128059
(
2022
).
50.
H.
Wang
,
J.
Qin
,
X.
Ding
,
M.-C.
Chen
,
S.
Chen
,
X.
You
,
Y.-M.
He
,
X.
Jiang
,
L.
You
,
Z.
Wang
,
C.
Schneider
,
J. J.
Renema
,
S.
Höfling
,
C.-Y.
Lu
, and
J.-W.
Pan
, “
Boson sampling with 20 input photons and a 60-mode interferometer in a 1014-dimensional Hilbert space
,”
Phys. Rev. Lett.
123
,
250503
(
2019
).
51.
P.
Grünwald
, “
Nonquantum information gain from higher-order correlation functions
,”
Phys. Rev. Res.
2
,
023147
(
2020
).
52.
E.
Zubizarreta Casalengua
,
J. C.
López Carreño
,
E.
del Valle
, and
F. P.
Laussy
, “
Structure of the harmonic oscillator in the space of n-particle Glauber correlators
,”
J. Math. Phys.
58
,
062109
(
2017
).
53.
R. H.
Brown
and
R. Q.
Twiss
, “
Correlation between photons in two coherent beams of light
,”
Nature
177
,
27
29
(
1956
).
54.
M.
Abadi
,
A.
Agarwal
,
P.
Barham
,
E.
Brevdo
,
Z.
Chen
,
C.
Citro
,
G. S.
Corrado
,
A.
Davis
,
J.
Dean
,
M.
Devin
,
S.
Ghemawat
,
I.
Goodfellow
,
A.
Harp
,
G.
Irving
,
M.
Isard
,
Y.
Jia
,
R.
Jozefowicz
,
L.
Kaiser
,
M.
Kudlur
,
J.
Levenberg
,
D.
Mane
,
R.
Monga
,
S.
Moore
,
D.
Murray
,
C.
Olah
,
M.
Schuster
,
J.
Shlens
,
B.
Steiner
,
I.
Sutskever
,
K.
Talwar
,
P.
Tucker
,
V.
Vanhoucke
,
V.
Vasudevan
,
F.
Viegas
,
O.
Vinyals
,
P.
Warden
,
M.
Wattenberg
,
M.
Wicke
,
Y.
Yu
, and
X.
Zheng
, “
TensorFlow: Large-scale machine learning on heterogeneous distributed systems
,” arXiv:1603.04467.
55.
R.
Loudon
,
The Quantum Theory of Light
(
Oxford University Press
,
Oxford
,
2000
).
56.
T.
Schapeler
,
N.
Lamberty
,
T.
Hummel
,
F.
Schlue
,
M.
Stefszky
,
B.
Brecht
,
C.
Silberhorn
, and
T. J.
Bartley
, “
Electrical trace analysis of superconducting nanowire photon-number-resolving detectors
,”
Phys. Rev. Appl.
22
,
014024
(
2024
).
57.
A. E.
Lita
,
A. J.
Miller
, and
S. W.
Nam
, “
Counting near-infrared single-photons with 95% efficiency
,”
Opt. Express
16
,
3032
3040
(
2008
).
58.
T.
Thomay
,
S. V.
Polyakov
,
O.
Gazzano
,
E.
Goldschmidt
,
Z. D.
Eldredge
,
T.
Huber
,
V.
Loo
, and
G. S.
Solomon
, “
Simultaneous, full characterization of a single-photon state
,”
Phys. Rev. X
7
,
041036
(
2017
).
59.
E. B.
Flagg
,
S. V.
Polyakov
,
T.
Thomay
, and
G. S.
Solomon
, “
Dynamics of nonclassical light from a single solid-state quantum emitter
,”
Phys. Rev. Lett.
109
,
163601
(
2012
).
60.
M.
Egmont-Petersen
,
D.
de Ridder
, and
H.
Handels
, “
Image processing with neural networks—A review
,”
Pattern Recognit.
35
,
2279
2301
(
2002
).
61.
Y.
LeCun
,
K.
Kavukcuoglu
, and
C.
Farabet
, “
Convolutional networks and applications in vision
,” in
Proceedings of 2010 IEEE International Symposium on Circuits and Systems
(
IEEE
,
2010
), pp.
253
256
.
62.
F.
Chollet
, .,
2024
, https://keras.io.
63.
J.
Nagi
,
F.
Ducatelle
,
G. A.
Di Caro
,
D.
Cireşan
,
U.
Meier
,
A.
Giusti
,
F.
Nagi
,
J.
Schmidhuber
, and
L. M.
Gambardella
, “
Max-pooling convolutional neural networks for vision-based hand gesture recognition
,” in
2011 IEEE International Conference on Signal and Image Processing Applications (ICSIPA)
(
IEEE
,
2011
), pp.
342
347
.
64.
L.
Stewart
,
F.
Bach
,
Q.
Berthet
, and
J.-P.
Vert
, “
Regression as classification: Influence of task formulation on neural network features
,” arXiv:2211 05641.
65.
N. S.
Altman
, “
An introduction to kernel and nearest-neighbor nonparametric regression
,”
Am. Stat.
46
,
175
185
(
1992
).
66.
S.
Dreiseitl
and
L.
Ohno-Machado
, “
Logistic regression and artificial neural network classification models: A methodology review
,”
J. Biomed. Inf.
35
,
352
359
(
2002
).
67.
A.
Natekin
and
A.
Knoll
, “
Gradient boosting machines, a tutorial
,”
Front. Neurorobotics
7
,
21
(
2013
).
68.
J.
Laaksonen
and
E.
Oja
, “
Classification with learning k-nearest neighbors
,” in
Proceedings of International Conference on Neural Networks (ICNN’96)
(
IEEE
,
1996
), Vol.
3
, pp.
1480
1483
.
69.
A.
Mucherino
,
P. J.
Papajorgji
, and
P. M.
Pardalos
, “
K-nearest neighbor classification
,” in
Data Mining in Agriculture
, edited by
A.
Mucherino
,
P. J.
Papajorgji
, and
P. M.
Pardalos
(
Springer
,
2009
), pp.
83
106
.
70.
W.
Cheng
and
E.
Hüllermeier
, “
Combining instance-based learning and logistic regression for multilabel classification
,”
Mach. Learn.
76
,
211
225
(
2009
).
71.
V. T.
Hai
,
P. H.
Chuong
, and
P. T.
Bao
, “
New approach of KNN algorithm in quantum computing based on new design of quantum circuits
,”
Informatica
46
,
95
(
2022
).
72.
S.
Gulia
,
M. T.
Beig
,
R.
Vatsa
, and
Y.
Sharma
, “
The dynamic integration of computational approaches and machine learning for cutting-edge solutions in photonics
,”
Interactions
245
,
81
(
2024
).
73.
P.
Linardatos
,
V.
Papastefanopoulos
, and
S.
Kotsiantis
, “
Explainable AI: A review of machine learning interpretability methods
,”
Entropy
23
,
18
(
2020
).
74.
M. J.
Kearns
,
The Computational Complexity of Machine Learning
(
MIT Press
,
1990
).
75.
P.
Cunningham
and
S. J.
Delany
, “
K-nearest neighbour classifiers—A tutorial
,”
ACM Comput. Surv.
54
,
128:1
128:25
(
2022
).
76.
D. C.
Ciresan
,
U.
Meier
,
J.
Masci
,
L. M.
Gambardella
, and
J.
Schmidhuber
, “
Flexible, high performance convolutional neural networks for image classification
,” in
Proceedings of the Twenty-Second International Joint Conference on Artificial Intelligence
(
ACM
,
2011
); Vol. 2, pp.
1237
1242
77.
Q.
Li
,
W.
Cai
,
X.
Wang
,
Y.
Zhou
,
D. D.
Feng
, and
M.
Chen
, “
Medical image classification with convolutional neural network
,” in
2014 13th International Conference on Control Automation Robotics and Vision (ICARCV)
(
IEEE
,
2014
), pp.
844
848
.
78.
T.
O’Malley
,
E.
Bursztein
,
J.
Long
,
F.
Chollet
,
H.
Jin
,
L.
Invernizzi
et al, Kerastuner,
2019
, https://github.com/keras-team/keras-tuner.
79.
R.
Chauhan
,
K. K.
Ghanshala
, and
R.
Joshi
, “
Convolutional neural network (CNN) for image detection and recognition
,” in
2018 First International Conference on Secure Cyber Computing and Communication (ICSCCC)
(
IEEE
,
2018
), pp.
278
282
.
80.
L.
Alzubaidi
,
J.
Zhang
,
A. J.
Humaidi
,
A.
Al-Dujaili
,
Y.
Duan
,
O.
Al-Shamma
,
J.
Santamaría
,
M. A.
Fadhel
,
M.
Al-Amidie
, and
L.
Farhan
, “
Review of deep learning: Concepts, CNN architectures, challenges, applications, future directions
,”
J. Big Data
8
,
53
(
2021
).
81.
M. I.
Jordan
and
T. M.
Mitchell
, “
Machine learning: Trends, perspectives, and prospects
,”
Science
349
,
255
260
(
2015
).
82.
A.
Paszke
,
S.
Gross
,
F.
Massa
,
A.
Lerer
,
J.
Bradbury
,
G.
Chanan
,
T.
Killeen
,
Z.
Lin
,
N.
Gimelshein
,
L.
Antiga
,
A.
Desmaison
,
A.
Kopf
,
E.
Yang
,
Z.
DeVito
,
M.
Raison
,
A.
Tejani
,
S.
Chilamkurthy
,
B.
Steiner
,
L.
Fang
,
J.
Bai
, and
S.
Chintala
, “
PyTorch: An imperative style, high-performance deep learning library
,” in
Advances in Neural Information Processing Systems
(
Curran Associates, Inc.
,
2019
), Vol.
32
.
83.
J.
Bruna
,
W.
Zaremba
,
A.
Szlam
, and
Y.
LeCun
, “
Spectral networks and locally connected networks on graphs
,” arXiv:1312.6203.
84.
T.-J.
Shih
,
W.-K.
Huang
,
Y.-M.
Lin
,
K.-B.
Li
,
C.-Y.
Hsu
,
J.-M.
Chen
,
P.-Y.
Tu
,
T.
Peters
,
Y.-F.
Chen
, and
I. A.
Yu
, “
Universal relation between the conditional auto-correlation function and the cross-correlation function of biphotons
,”
Opt. Express
32
,
13657
13671
(
2024
).
85.
S.
Zhao
,
J.
Lavie
,
L.
Rondin
,
L.
Orcin-Chaix
,
C.
Diederichs
,
P.
Roussignol
,
Y.
Chassagneux
,
C.
Voisin
,
K.
Müllen
,
A.
Narita
,
S.
Campidelli
, and
J.-S.
Lauret
, “
Single photon emission from graphene quantum dots at room temperature
,”
Nat. Commun.
9
,
3470
(
2018
).
86.
T. T.
Tran
,
K.
Bray
,
M. J.
Ford
,
M.
Toth
, and
I.
Aharonovich
, “
Quantum emission from hexagonal boron nitride monolayers
,”
Nat. Nanotechnol.
11
,
37
41
(
2016
).
87.
A.
Srivastava
,
M.
Sidler
,
A. V.
Allain
,
D. S.
Lembke
,
A.
Kis
, and
A.
Imamoğlu
, “
Optically active quantum dots in monolayer WSe2
,”
Nat. Nanotechnol.
10
,
491
496
(
2015
).
88.
H.
Yu
,
G.-B.
Liu
,
J.
Tang
,
X.
Xu
, and
W.
Yao
, “
Moiré excitons: From programmable quantum emitter arrays to spin-orbit–coupled artificial lattices
,”
Sci. Adv.
3
,
e1701696
(
2017
).
89.
Y.-M.
He
,
G.
Clark
,
J. R.
Schaibley
,
Y.
He
,
M.-C.
Chen
,
Y.-J.
Wei
,
X.
Ding
,
Q.
Zhang
,
W.
Yao
,
X.
Xu
,
C.-Y.
Lu
, and
J.-W.
Pan
, “
Single quantum emitters in monolayer semiconductors
,”
Nat. Nanotechnol.
10
,
497
502
(
2015
).
90.
P.
Tonndorf
,
R.
Schmidt
,
R.
Schneider
,
J.
Kern
,
M.
Buscema
,
G. A.
Steele
,
A.
Castellanos-Gomez
,
H. S. J.
van der Zant
,
S.
Michaelis de Vasconcellos
, and
R.
Bratschitsch
, “
Single-photon emission from localized excitons in an atomically thin semiconductor
,”
Optica
2
,
347
352
(
2015
).
91.
M.
Uria
,
P.
Solano
, and
C.
Hermann-Avigliano
, “
Deterministic generation of large Fock states
,”
Phys. Rev. Lett.
125
,
093603
(
2020
).
92.
L. M.
Hansen
,
L.
Carosini
,
L.
Jehle
,
F.
Giorgino
,
R.
Houvenaghel
,
M.
Vyvlecka
,
J. C.
Loredo
, and
P.
Walther
, “
Single-active-element demultiplexed multi-photon source
,”
Opt. Quantum
1
,
1
5
(
2023
).
93.
S. J.
Rosenthal
,
J. C.
Chang
,
O.
Kovtun
,
J. R.
McBride
, and
I. D.
Tomlinson
, “
Biocompatible quantum dots for biological applications
,”
Chem. Biol.
18
,
10
24
(
2011
).
94.
M.
Monici
, “
Cell and tissue autofluorescence research and diagnostic applications
,”
Biotechnol. Annu. Rev.
11
,
227
256
(
2005
).