Articulo de referencia

Agrupamiento cuántico

El agrupamiento cuántico (QC) es una clase de algoritmos de agrupamiento de datos que utilizan herramientas conceptuales y matemáticas de la mecánica cuántica . QC pertenece a l...

El agrupamiento cuántico (QC) es una clase de algoritmos de agrupamiento de datos que utilizan herramientas conceptuales y matemáticas de la mecánica cuántica . QC pertenece a la familia de algoritmos de agrupamiento basados ​​en densidad , donde los clústeres se definen por regiones de mayor densidad de puntos de datos.

QC fue desarrollado por primera vez por David Horn y Assaf Gottlieb en 2001. [ 1 ]

Algoritmo de agrupamiento cuántico original

Dado un conjunto de puntos en un espacio de datos n-dimensional, QC representa cada punto con una distribución gaussiana multidimensional , con un ancho (desviación estándar) sigma, centrada en la ubicación de cada punto en el espacio. Estas gaussianas se suman para crear una única distribución para todo el conjunto de datos. (Este paso es un ejemplo particular de estimación de densidad de kernel , a menudo denominada estimador de ventana de Parzen-Rosenblatt). Esta distribución se considera la función de onda cuántica para el conjunto de datos. En términos generales, la función de onda es una descripción generalizada de dónde es probable que haya puntos de datos en el espacio.

A continuación, la computación cuántica introduce el concepto de potencial cuántico ; utilizando la ecuación de Schrödinger independiente del tiempo , se construye una superficie de potencial cuya función de onda del conjunto de datos es una solución estable. Los detalles de la superficie de potencial son más robustos ante cambios en sigma (el ancho de las gaussianas) que los detalles correspondientes en la función de onda; esta ventaja es una de las motivaciones iniciales para el desarrollo de la computación cuántica.

La superficie potencial se considera el "paisaje" del conjunto de datos, donde los puntos "bajos" del paisaje corresponden a regiones de alta densidad de datos. Posteriormente, QC utiliza el descenso de gradiente para desplazar cada punto de datos "cuesta abajo" en el paisaje, lo que provoca que los puntos se agrupen en mínimos cercanos , revelando así clústeres dentro del conjunto de datos.

QC has a single main hyperparameter, which is the width sigma of the Gaussian distribution around each data point. For sufficiently small sigma, every data point will define its own depression in the landscape, and no points will move, thus creating no clusters. For sufficiently large sigma, the landscape becomes a single smooth bowl, and every data point will cluster together at the single global minimum in the landscape. Exploring the range of sigma values between these extremes yields information about the inherent structure of the data set, including hierarchy of structure; smaller sigma values reveal more fine-grained local structure, and larger sigma values reveal overall global structure. The QC algorithm does not specify a preferred or 'correct' value of sigma.

Dynamic Quantum Clustering

Developed by Marvin Weinstein and David Horn in 2009,[2]Dynamic Quantum Clustering (DQC) extends the basic QC algorithm in several ways.

Quantum evolution, non-local gradient descent, and tunneling

DQC uses the same potential landscape as QC, but it replaces classical gradient descent with quantum evolution. To do this, each data point is again represented by its individual wave function (a multidimensional Gaussian distribution with width sigma). The time-dependent Schrödinger equation is then used to compute each wave function's evolution over time in the given quantum potential. More precisely, a small time-step value is introduced, and the evolution of the wave function is calculated repeatedly at each time step, with a new expected location for the data point calculated after each step. This process builds a trajectory through the data space for each point; the evolution continues until all points have stopped moving.

Importantly, the Ehrenfest theorem from quantum mechanics states that this quantum evolution does, in fact, equate to the point moving downhill in the potential landscape, in expectation. The "in expectation" part is important because, unlike in classical physics, the point's motion is not influenced only by the gradient of the potential at the point's location; instead, the point's wave function extends over the entire landscape (with the Gaussian centered at the point's location), and a complex interaction between the wave function and the potential determines the point's motion. As a loose analogy: regions of the landscape that are below the point's current location 'attract' the point—the more so the lower the region is, but the less so the farther away from the point it is. In the same way, higher regions of the landscape 'repel' the point.

Thus, quantum evolution for each point acts as a form of non-local gradient descent in the potential. This non-locality creates the possibility of tunneling, where a point will seem to ignore or pass through a potential barrier on its way toward some lower minimum. The biggest problem in non-convex gradient descent is often the existence of many small and uninteresting local minima where points can get stuck as they descend. (This problem tends to get worse as the number of dimensions increases, which is part of the curse of dimensionality.) DQC's use of non-local gradient descent and tunneling presents a solution to this problem.

DQC introduces two new hyperparameters: the time step, and the mass of each data point (which controls the degree of tunneling behavior). Whereas tuning of sigma is integral to understanding any new data set, both time step and mass can usually be left at reasonable default values and still produce useful results.

An important downside of the quantum-evolution approach is that the time complexity of the evolution is now O(n3){\displaystyle O(n^{3})} in the number of data points, since interacting with the entire potential landscape is O(n2){\displaystyle O(n^{2})} for each point. For large data sets, the computation time would quickly become intractable. When needed, DQC addresses this problem by selecting a limited number of points from the data set to act as a basis (see next section).

Use of a limited basis

For a data set of n points, DQC creates a set of n quantum eigenstates for use in its underlying calculations; the eigenstates are orthonormal, and each one is a linear combination of the Gaussian distributions representing each data point.

For large n, the use of n eigenstates becomes computationally intractable, since the creation of the potential and the evolution of individual points are both O(n3){\displaystyle O(n^{3})}. To solve this problem, DQC allows for the selection of a more limited basis, as follows. To act as the basis, b data points (b < n) are chosen such that the basis points span the space occupied by the data set. (This can be done in multiple ways; the essential goal is to choose the basis points to be as far away from each other as possible.) DQC then constructs b eigenstates using these b basis points. These eigenstates represent the basis points perfectly; they are also used to represent all non-basis points, but those representations will be imperfect. The loss of information is relative to sigma; for a given basis, sigma must be chosen large enough that the basis can be used to represent the non-basis points to some reasonable degree of accuracy. In a sense, the size of the chosen basis can be thought of as the 'resolution' that is used to 'view' the structure of the data.

The largest reasonable basis size depends on available computing resources, and on how long one is willing to wait for results. As of 2020, without access to enterprise-level computing resources, the largest tractable basis size is typically in the range of 1,500–2,000 points.

Use of dynamic visualization

DQC calculates a trajectory for each data point, which allows for the creation of an animated ('dynamic') visualization in which all data points move along their trajectories simultaneously. These animations present information not just in the final destination of each point, but in each entire trajectory along the way. Notably, animations can reveal the existence of channels leading to a given cluster. (In the landscape metaphor, these structures can be thought of as riverbeds and lakes.) Although a given visualization is limited to at most 3 spatial dimensions, the appearance of channels and other structures creates the ability to 'see' what is happening in more than 3 dimensions.

Use of a PCA coordinate system is helpful for these visualizations; viewing the trajectories in the first 3 PCA dimensions packs as much information as possible into a single visualization.

A given 3D visualization of these trajectories is not an embedding of the trajectories in 3 dimensions. The trajectories have the same dimensionality as the data space, which is often much larger than 3; the visualization is simply a 3D view into a higher-dimensional motion.

The channels ('riverbeds') can have meaning in two different ways. First, they can be treated as subclusters, where different subclusters join the main cluster from different directions. Second, they can be treated as regressions: position along the channel at a given time (or, equivalently, order of arrival at the cluster center) may be correlated with some metadata of interest.

Applications

Variants of QC have been applied to real-world data in many fields, including biology,[1][2][3][4][5][6] geology,[3][7] physics,[3][4][8] finance,[3] engineering,[4] and economics.[9] With these applications, a comprehensive mathematical analysis to find all the roots of the quantum potential has also been worked out.[10]

References

  1. 12Horn, D.; Gottlieb, A. (2001). "Algorithm for Data Clustering in Pattern Recognition Problems Based on Quantum Mechanics". Physical Review Letters. 88 (1) 018702. Bibcode:2001PhRvL..88a8702H. doi:10.1103/PhysRevLett.88.018702. PMID 11800996.
  2. 12Weinstein, M.; Horn, D. (2009). "Dynamic quantum clustering: a method for visual exploration of structures in data". Physical Review E. 80 (6) 066117. arXiv:0908.2644. Bibcode:2009PhRvE..80f6117W. doi:10.1103/PhysRevE.80.066117. PMID 20365241. S2CID 10550999.
  3. 1234Weinstein, M.; Meirer, F.; Hume, A.; Sciau, Ph.; Shaked, G.; Hofstetter, R.; Persi, E.; Mehta, A.; Horn, D. (2013). "Analyzing Big Data with Dynamic Quantum Clustering". arXiv:1310.2700 [physics.data-an].
  4. 1 2 3 Scott, TC; Therani, M.; Wang, XM (2017). "Agrupamiento de datos con mecánica cuántica" . Matemáticas . 5 (1): 5. doi : 10.3390/math5010005 .
  5. Roche, K.; Weinstein, M.; Dunwoodie, LJ; Poehlman, WL; Feltus, FA (2018). "La clasificación de cinco tipos de tumores humanos revela biomarcadores específicos y genes de clasificación de fondo" . Scientific Reports . 8 (1): 8180. Bibcode : 2018NatSR...8.8180R . doi : 10.1038/s41598-018-26310- x . PMC 5970138. PMID 29802335 .  
  6. Casaña-Eslava, RV; Lisboa, PJG; Ortega-Martorell, S.; Jarman, HI; Martín-Guerrero, JD (2020). "Un marco probabilístico para la agrupación cuántica". Sistemas basados ​​en el conocimiento . 194 . arXiv : 1902.05578 . doi : 10.1016/j.knosys.2020.105567 . S2CID 213468799 . 
  7. Shaked, G. (2013). Agrupamiento cuántico de grandes conjuntos de datos (PDF) (M.Sc.).
  8. Weinstein, M.; Heifetz, A.; Klann, R. (2014). "Detección de fuentes nucleares en un estudio de búsqueda utilizando agrupamiento cuántico dinámico de datos espectrales de rayos gamma". The European Physical Journal Plus . 129 (11): 239. arXiv : 1406.0746 . Bibcode : 2014EPJP..129..239W . doi : 10.1140/epjp/i2014-14239-3 . S2CID 119217077 . 
  9. Decheng, F.; Jon, S.; Pang, C.; Dong, W.; Won, C. (2018). "Análisis de agrupamiento cuántico mejorado basado en la distancia ponderada y su aplicación" . Heliyon . 4 ( 11) e00984. Bibcode : 2018Heliy...400984D . doi : 10.1016/j.heliyon.2018.e00984 . PMC 6275214. PMID 30761372 .  
  10. Maignan, A.; Scott, TC (2021). "Un análisis exhaustivo del agrupamiento cuántico: encontrar todos los mínimos potenciales" (PDF) . Revista internacional de minería de datos y gestión del conocimiento . 11 (1): 33-54. doi : 10.5121/ijdkp.2021.11103 .
  • Vídeo (YouTube.com): "Explorando datos complejos de alta dimensión para descubrir estructuras ocultas", Marvin Weinstein (Charlas SETI)
  • Vídeo (YouTube.com): "Descubriendo estructuras asombrosas ocultas en grandes conjuntos de datos complejos y densos", Marvin Weinstein (Charlas SETI)
  • Vídeo (YouTube.com): "Agrupamiento cuántico: algoritmo de agrupamiento inspirado en la física", Sigalit Bechler
  • Vídeo (YouTube.com): "Información cuántica a partir de conjuntos de datos complejos", Marvin Weinstein (Charlas en Google)