Una partícula clásica tiene una posición y un momento definidos, y por lo tanto se representa mediante un punto en el espacio de fases. Dada una colección ( conjunto ) de partículas, la probabilidad de encontrar una partícula en una posición determinada en el espacio de fases se especifica mediante una distribución de probabilidad, la densidad de Liouville. Esta interpretación estricta falla para una partícula cuántica, debido al principio de incertidumbre . En cambio, la distribución de Wigner de cuasiprobabilidad anterior desempeña un papel análogo, pero no satisface todas las propiedades de una distribución de probabilidad convencional; y, a la inversa, satisface propiedades de acotación que no están disponibles para las distribuciones clásicas.
Por ejemplo, la distribución de Wigner puede y normalmente toma valores negativos para estados que no tienen un modelo clásico, y es un indicador conveniente de interferencia mecánico-cuántica. (Véase más abajo una caracterización de estados puros cuyas funciones de Wigner no son negativas). Suavizar la distribución de Wigner a través de un filtro de tamaño mayor que ħ (por ejemplo, convolucionando con una gaussiana de espacio de fase, una transformada de Weierstrass , para producir la representación de Husimi, más abajo), da como resultado una función semidefinida positiva, es decir, se puede pensar que se ha convertido en una función semiclásica. [a]
Se puede demostrar (mediante la convolución con una gaussiana pequeña) que las regiones de ese valor negativo son "pequeñas": no pueden extenderse a regiones compactas mayores que unos pocos ħ y, por lo tanto, desaparecen en el límite clásico . Están protegidas por el principio de incertidumbre , que no permite una ubicación precisa dentro de regiones del espacio de fases menores que ħ y, por lo tanto, hace que esas " probabilidades negativas " sean menos paradójicas.
Definición y significado
La distribución de Wigner W ( x , p ) de un estado puro se define como
donde ψ es la función de onda, y x y p son la posición y el momento, pero podría ser cualquier par de variables conjugadas (por ejemplo, partes reales e imaginarias del campo eléctrico o frecuencia y tiempo de una señal). Nótese que puede tener soporte en x incluso en regiones donde ψ no tiene soporte en x ("pulsaciones").
Es simétrico en x y p :
donde φ es la función de onda del espacio-momento normalizada, proporcional a la transformada de Fourier de ψ .
Este límite desaparece en el límite clásico, ħ → 0. En este límite, W ( x , p ) se reduce a la densidad de probabilidad en el espacio de coordenadas x , normalmente muy localizada, multiplicada por las funciones δ en el momento: el límite clásico es "puntiagudo". Por lo tanto, este límite mecánico-cuántico excluye una función de Wigner que sea una función δ perfectamente localizada en el espacio de fases, como reflejo del principio de incertidumbre. [6]
10. La transformación de Wigner es simplemente la transformada de Fourier de las antidiagonales de la matriz de densidad, cuando esa matriz se expresa en base a una posición. [7]
Esto puede derivarse de la expresión para las funciones de onda de estados propios estáticos,
donde es el polinomio de Hermite -ésimo . De la definición anterior de la función de Wigner, al cambiar las variables de integración,
La expresión se deduce entonces de la relación integral entre los polinomios de Hermite y Laguerre. [8]
Ecuación de evolución de la función de Wigner
La transformación de Wigner es una transformación invertible general de un operador Ĝ en un espacio de Hilbert a una función g ( x , p ) en el espacio de fases y está dada por
Los operadores hermíticos se traducen en funciones reales. La inversa de esta transformación, del espacio de fases al espacio de Hilbert, se denomina transformación de Weyl :
La función de Wigner W ( x , p ) analizada aquí se ve así como la transformada de Wigner del operador de matriz de densidad ρ̂ . Por lo tanto, la traza de un operador con la transformada de Wigner de la matriz de densidad es la superposición integral del espacio de fase equivalente de g ( x , p ) con la función de Wigner.
donde H ( x , p ) es el hamiltoniano y {{⋅, ⋅}} es el corchete de Moyal . En el límite clásico, ħ → 0 , el corchete de Moyal se reduce al corchete de Poisson , mientras que esta ecuación de evolución se reduce a la ecuación de Liouville de la mecánica estadística clásica.
Formalmente, la ecuación clásica de Liouville se puede resolver en términos de las trayectorias de las partículas en el espacio de fases, que son soluciones de las ecuaciones clásicas de Hamilton. Esta técnica de resolución de ecuaciones diferenciales parciales se conoce como el método de las características . Este método se traslada a los sistemas cuánticos, donde las "trayectorias" de las características determinan ahora la evolución de las funciones de Wigner. La solución de la ecuación de evolución de Moyal para la función de Wigner se representa formalmente como
donde y son las trayectorias características sujetas a las ecuaciones cuánticas de Hamilton con condiciones iniciales y , y donde la composición del producto se entiende para todas las funciones de argumento.
Dado que la composición de funciones es completamente no local (el "fluido de probabilidad cuántica" se difunde, como observó Moyal), los vestigios de trayectorias locales en los sistemas cuánticos son apenas discernibles en la evolución de la función de distribución de Wigner. [b] En la representación integral de los productos, las operaciones sucesivas de ellos se han adaptado a una integral de trayectorias en el espacio de fases, para resolver la ecuación de evolución de la función de Wigner [9] (véase también [10] [11] [12] ). Esta característica no local de la evolución temporal de Moyal [13] se ilustra en la galería siguiente, para hamiltonianos más complejos que el oscilador armónico. En el límite clásico, la naturaleza de trayectoria de la evolución temporal de las funciones de Wigner se vuelve cada vez más clara. En ħ = 0, las trayectorias de las características se reducen a las trayectorias clásicas de partículas en el espacio de fases.
Ejemplos de evoluciones temporales de la función de Wigner
Estado puro en un potencial cuártico. Las líneas continuas representan el conjunto de niveles del hamiltoniano.
Tunelización de un paquete de ondas a través de una barrera de potencial. Las líneas continuas representan el conjunto de niveles del hamiltoniano.
Evolución a largo plazo de un estado mixto ρ en un pozo de potencial anarmónico. Los valores marginales se representan a la derecha ( p ) y en la parte superior ( x ).
Un estado mixto de equilibrio ρ (evoluciona hacia sí mismo), en el mismo potencial anarmónico.
Evolución temporal del oscilador armónico
Sin embargo, en el caso especial del oscilador armónico cuántico , la evolución es simple y parece idéntica al movimiento clásico: una rotación rígida en el espacio de fases con una frecuencia dada por la frecuencia del oscilador. Esto se ilustra en la galería siguiente. Esta misma evolución temporal ocurre con los estados cuánticos de los modos de luz , que son osciladores armónicos.
Ejemplos de evoluciones temporales de la función Wigner en un oscilador armónico cuántico
Distribución de cuasiprobabilidad de Wigner de un estado coherente.
Estados de gatos con 2, 3, 4 gatos. La separación entre gatos varía de 0,5, 1, 2, 4, mostrando una inferencia cada vez más nítida.
Un estado con muchos felinos, con 10 gatos separados en .
Fock afirma
Estados de Fock superpuestos, por ejemplo .
.
.
Estados comprimidos
Distribución de cuasiprobabilidad de Wigner de estados comprimidos, para cantidades variables de desplazamiento y cambio de fase. Cuando el cambio de fase es cero (columna más a la izquierda), la amplitud de oscilación de pico a pico es menos incierta, pero la fase (el momento en el que cruza el punto medio) es más incierta. Lo opuesto es cierto cuando el cambio de fase es de medio ciclo (columna más a la derecha).
Pureza variable, desde totalmente puro a totalmente mezclado.
Límite clásico
La función de Wigner permite estudiar el límite clásico , ofreciendo una comparación de la dinámica clásica y cuántica en el espacio de fases. [14] [15]
Se ha sugerido que el enfoque de la función de Wigner puede verse como una analogía cuántica de la formulación operacional de la mecánica clásica introducida en 1932 por Bernard Koopman y John von Neumann : la evolución temporal de la función de Wigner se acerca, en el límite ħ → 0, a la evolución temporal de la función de onda de Koopman-von Neumann de una partícula clásica. [16]
Los momentos de la función de Wigner generan promedios de operadores simetrizados, en contraste con el orden normal y el orden antinormal generados por la representación P de Glauber-Sudarshan y la representación Q de Husimi respectivamente. Por lo tanto, la representación de Wigner es muy adecuada para realizar aproximaciones semiclásicas en óptica cuántica [17] y teoría de campos de condensados de Bose-Einstein donde la ocupación de modos altos se acerca a un límite semiclásico. [18]
Positividad de la función de Wigner
Como ya se ha señalado, la función de Wigner del estado cuántico suele adoptar algunos valores negativos. De hecho, para un estado puro en una variable, si para todos y , entonces la función de onda debe tener la forma
para algunos números complejos con (teorema de Hudson [19] ). Nótese que se permite que sea complejo. En otras palabras, es un paquete de ondas gaussianas unidimensionales . Por lo tanto, los estados puros con funciones de Wigner no negativas no son necesariamente estados de incertidumbre mínima en el sentido de la fórmula de incertidumbre de Heisenberg ; más bien, dan igualdad en la fórmula de incertidumbre de Schrödinger , que incluye un término anticonmutador además del término conmutador. (Con una definición cuidadosa de las varianzas respectivas, todas las funciones de Wigner de estado puro conducen a la desigualdad de Heisenberg de todos modos).
En dimensiones superiores, la caracterización de estados puros con funciones de Wigner no negativas es similar; la función de onda debe tener la forma
donde es una matriz compleja simétrica cuya parte real es definida positiva, es un vector complejo y c es un número complejo. [20] La función de Wigner de cualquier estado de este tipo es una distribución gaussiana en el espacio de fases.
Soto y Claverie [20] dan una prueba elegante de esta caracterización, utilizando la transformada de Segal-Bargmann . El razonamiento es el siguiente. La función Q de Husimi de puede calcularse como la magnitud al cuadrado de la transformada de Segal-Bargmann de , multiplicada por una gaussiana. Mientras tanto, la función Q de Husimi es la convolución de la función de Wigner con una gaussiana. Si la función de Wigner de es no negativa en todas partes del espacio de fases, entonces la función Q de Husimi será estrictamente positiva en todas partes del espacio de fases. Por lo tanto, la transformada de Segal-Bargmann de no será cero en ninguna parte. Por lo tanto, por un resultado estándar del análisis complejo, tenemos
para alguna función holomorfa . Pero para pertenecer al espacio de Segal-Bargmann —es decir, para ser integrable al cuadrado con respecto a una medida gaussiana— debe tener como máximo un crecimiento cuadrático en el infinito. A partir de esto, se puede utilizar el análisis complejo elemental para demostrar que debe ser en realidad un polinomio cuadrático. Así, obtenemos una forma explícita de la transformada de Segal-Bargmann de cualquier estado puro cuya función de Wigner no sea negativa. A continuación, podemos invertir la transformada de Segal-Bargmann para obtener la forma reclamada de la función de onda de posición.
No parece haber ninguna caracterización simple de estados mixtos con funciones de Wigner no negativas.
La función de Wigner en relación con otras interpretaciones de la mecánica cuántica
Se ha demostrado que la función de distribución de cuasiprobabilidad de Wigner puede considerarse como una deformación ħ de otra función de distribución del espacio de fases que describe un conjunto de trayectorias causales de De Broglie-Bohm . [21] Basil Hiley ha demostrado que la distribución de cuasiprobabilidad puede entenderse como la matriz de densidad reexpresada en términos de una posición media y un momento de una "celda" en el espacio de fases, y la interpretación de De Broglie-Bohm permite describir la dinámica de los centros de dichas "celdas". [22] [23]
Existe una estrecha conexión entre la descripción de estados cuánticos en términos de la función de Wigner y un método de reconstrucción de estados cuánticos en términos de bases mutuamente imparciales . [24]
Usos de la función de Wigner fuera de la mecánica cuántica
En el modelado de sistemas ópticos como telescopios o dispositivos de telecomunicaciones de fibra, la función de Wigner se utiliza para salvar la brecha entre el simple trazado de rayos y el análisis de onda completa del sistema. Aquí p / ħ se reemplaza con k = | k | sen θ ≈ | k | θ en la aproximación de ángulo pequeño (paraxial). En este contexto, la función de Wigner es lo más cercano que se puede llegar a la descripción del sistema en términos de rayos en la posición x y el ángulo θ mientras que todavía se incluyen los efectos de la interferencia. [25] Si se vuelve negativa en cualquier punto, entonces el simple trazado de rayos no será suficiente para modelar el sistema. Es decir, los valores negativos de esta función son un síntoma del límite de Gabor de la señal de luz clásica y no de las características cuánticas de la luz asociadas con ħ .
En el análisis de señales , una señal eléctrica, vibración mecánica u onda sonora que varía con el tiempo se representa mediante una función de Wigner . En este caso, x se reemplaza por el tiempo y p / ħ se reemplaza por la frecuencia angular ω = 2π f , donde f es la frecuencia regular.
En la óptica ultrarrápida , los pulsos láser cortos se caracterizan con la función de Wigner utilizando las mismas sustituciones f y t que se mencionaron anteriormente. Los defectos de pulso, como el chirrido (el cambio de frecuencia con el tiempo), se pueden visualizar con la función de Wigner. Véase la figura adyacente.
En óptica cuántica, x y p / ħ se reemplazan por las cuadraturas X y P , los componentes real e imaginario del campo eléctrico (ver estado coherente ).
Otras distribuciones de cuasiprobabilidad relacionadas
La distribución de Wigner fue la primera distribución de cuasiprobabilidad que se formuló, pero le siguieron muchas más, formalmente equivalentes y transformables a partir de ella (véase Transformación entre distribuciones en análisis de tiempo-frecuencia ). Como en el caso de los sistemas de coordenadas, debido a la variación de sus propiedades, varias de ellas presentan diversas ventajas para aplicaciones específicas:
Sin embargo, en cierto sentido, la distribución de Wigner ocupa una posición privilegiada entre todas estas distribuciones, ya que es la única cuyo producto estrella requerido desaparece (se integra por partes hasta la unidad efectiva) en la evaluación de los valores esperados, como se ilustra arriba, y por eso puede visualizarse como una medida de cuasiprobabilidad análoga a las clásicas.
Nota histórica
Como se ha indicado, la fórmula de la función de Wigner se derivó de forma independiente varias veces en diferentes contextos. De hecho, aparentemente, Wigner no sabía que, incluso en el contexto de la teoría cuántica, ya la habían introducido previamente Heisenberg y Dirac [26] [27] , aunque de manera puramente formal: estos dos no comprendieron su importancia, ni la de sus valores negativos, ya que simplemente la consideraron como una aproximación a la descripción cuántica completa de un sistema como el átomo. (Por cierto, Dirac se convertiría más tarde en cuñado de Wigner, casándose con su hermana Manci ). Simétricamente, en la mayor parte de su legendaria correspondencia de 18 meses con Moyal a mediados de la década de 1940, Dirac no sabía que la función generadora de momentos cuánticos de Moyal era efectivamente la función de Wigner, y fue Moyal quien finalmente le llamó la atención al respecto. [28]
^ En concreto, dado que esta convolución es invertible, de hecho, no se ha sacrificado ninguna información y la entropía cuántica completa aún no ha aumentado. Sin embargo, si esta distribución de Husimi resultante se utiliza como una medida simple en una evaluación integral del espacio de fases de los valores esperados sin el producto estrella requerido de la representación de Husimi , entonces, en esa etapa, se ha perdido la información cuántica y la distribución es, efectivamente, semiclásica . Es decir, dependiendo de su uso en la evaluación de los valores esperados, la misma distribución puede servir como una función de distribución cuántica o clásica .
^ Las características cuánticas no deben confundirse con las trayectorias de la integral de trayectoria de Feynman o las trayectorias de la teoría de De Broglie-Bohm . Esta triple ambigüedad permite comprender mejor la postura de Niels Bohr , quien se opuso enérgicamente, pero de manera contraproducente, a la noción de trayectoria en la física atómica. En la Conferencia de Pocono de 1948, por ejemplo, le dijo a Richard Feynman : "... no se podía hablar de la trayectoria de un electrón en el átomo, porque era algo no observable". ("The Beat of a Different Drum: The Life and Science of Richard Feynman", de Jagdish Mehra (Oxford, 1994, pp. 245-248)). Argumentos de este tipo fueron ampliamente utilizados en el pasado por Ernst Mach en su crítica de una teoría atómica de la física y más tarde, en la década de 1960, por Geoffrey Chew , Tullio Regge y otros para motivar la sustitución de la teoría cuántica de campos locales por la teoría de la matriz S. Hoy en día, la física estadística basada completamente en conceptos atomísticos se incluye en los cursos estándar, la teoría de la matriz S pasó de moda, mientras que el método de la integral de trayectorias de Feynman ha sido reconocido como el método más eficiente en las teorías de calibre .
Referencias
^ EP Wigner (1932). "Sobre la corrección cuántica para el equilibrio termodinámico". Physical Review . 40 (5): 749–759. Bibcode :1932PhRv...40..749W. doi :10.1103/PhysRev.40.749. hdl : 10338.dmlcz/141466 .
^ HJ Groenewold (1946). "Sobre los principios de la mecánica cuántica elemental". Physica . 12 (7): 405–460. Bibcode :1946Phy....12..405G. doi :10.1016/S0031-8914(46)80059-4.
^ H. Weyl (1927). "Quantenmechanik und gruppentheorie". Zeitschrift für Physik . 46 (1–2): 1. Bibcode : 1927ZPhy...46....1W. doi :10.1007/BF02055756. S2CID 121036548.; H. Weyl, Gruppentheorie und Quantenmechanik (Leipzig: Hirzel) (1928); H. Weyl, La teoría de los grupos y la mecánica cuántica (Dover, Nueva York, 1931).
^ J. Ville, "Théorie et Applications de la Notion de Signal Analytique", Câbles et Transmission , 2 , 61–74 (1948).
^ Moyal, JE (1949). "La mecánica cuántica como teoría estadística". Actas matemáticas de la Cambridge Philosophical Society . 45 (1). Cambridge University Press (CUP): 99–124. Bibcode :1949PCPS...45...99M. doi :10.1017/s0305004100000487. ISSN 0305-0041. S2CID 124183640.
^ Curtright, TL; Zachos, CK (2012). "Mecánica cuántica en el espacio de fases". Boletín de Física de Asia Pacífico . 1 : 37. arXiv : 1104.5269 . doi :10.1142/S2251158X12000069. S2CID 119230734.; C. Zachos , D. Fairlie y T. Curtright , Mecánica cuántica en el espacio de fases (World Scientific, Singapur, 2005). ISBN 978-981-238-384-6 .
^ Hawkes, Peter W. (2018). Avances en imágenes y física electrónica . Academic Press. pág. 47. ISBN9780128155424.
^ Schleich, Wolfgang P. (9 de febrero de 2001). Óptica cuántica en el espacio de fases (1.ª ed.). Wiley. pág. 105. doi :10.1002/3527602976. ISBN978-3-527-29435-0.
^ B. Leaf (1968). "Transformada de Weyl en dinámica cuántica no relativista". Revista de Física Matemática . 9 (5): 769–781. Código Bibliográfico :1968JMP.....9..769L. doi :10.1063/1.1664640.
^ P. Sharan (1979). "Representación de integrales de trayectorias mediante productos estelares". Physical Review D . 20 (2): 414–418. Bibcode :1979PhRvD..20..414S. doi :10.1103/PhysRevD.20.414.
^ MS Marinov (1991). "Un nuevo tipo de integral de trayectorias en el espacio de fases". Physics Letters A . 153 (1): 5–11. Bibcode :1991PhLA..153....5M. doi :10.1016/0375-9601(91)90352-9.
^ B. Segev: Núcleos evolutivos para distribuciones de espacio de fases . En: MA Olshanetsky; Arkady Vainshtein (2002). Múltiples facetas de cuantificación y supersimetría: volumen conmemorativo de Michael Marinov. World Scientific. págs. 68–90. ISBN.978-981-238-072-2. Recuperado el 26 de octubre de 2012 .Véase especialmente la sección 5. "Integral de trayectoria para el propagador" en las páginas 86 a 89. También disponible en línea.
^ M. Oliva, D. Kakofengitis y O. Steuernagel (2018). "Los sistemas mecánicos cuánticos anarmónicos no presentan trayectorias en el espacio de fases". Physica A . 502 : 201–210. arXiv : 1611.03303 . Bibcode :2018PhyA..502..201O. doi :10.1016/j.physa.2017.10.047. S2CID 53691877.{{cite journal}}: CS1 maint: multiple names: authors list (link)
^ Véase, por ejemplo: Wojciech H. Zurek , Decoherence and the transition from quantum to classic – revisited , Los Alamos Science, 27, 2002, arXiv:quant-ph/0306072, pp. 15 y siguientes.
^ Véase, por ejemplo: C. Zachos, D. Fairlie, T. Curtright, Mecánica cuántica en el espacio de fases: una visión general con artículos seleccionados , World Scientific, 2005. ISBN 978-981-4520-43-0 .
^ Bondar, Denys I.; Cabrera, Renan; Zhdanov, Dmitry V.; Rabitz, Herschel A. (2013). "Distribución del espacio de fases de Wigner como función de onda". Physical Review A . 88 (5): 052108. arXiv : 1202.3628 . doi :10.1103/PhysRevA.88.052108. ISSN 1050-2947. S2CID 119155284.
^ Graham, R (1973). "Teoría estadística de inestabilidades en sistemas estacionarios fuera de equilibrio con aplicaciones a láseres y óptica no lineal". Springer Tracts in Modern Physics . Springer. pág. 1.
^ Blakie†, PB; Bradley†, AS; Davis, MJ; Ballagh, RJ; Gardiner, CW (1 de septiembre de 2008). "Dinámica y mecánica estadística de gases de Bose ultrafríos utilizando técnicas de campo c". Avances en Física . 57 (5): 363–455. arXiv : 0809.1487 . doi :10.1080/00018730802564254. ISSN 0001-8732 . Consultado el 1 de marzo de 2023 .
^ Hudson, Robin L. (1974). "¿Cuándo la densidad de cuasi-probabilidad de Wigner es no negativa?". Informes de física matemática . 6 (2): 249–252. Bibcode :1974RpMP....6..249H. doi :10.1016/0034-4877(74)90007-X.
^ ab F. Soto y P. Claverie, "¿Cuándo es no negativa la función de Wigner de sistemas multidimensionales?", Journal of Mathematical Physics 24 (1983) 97–100.
^ Dias, Nuno Costa; Prata, João Nuno (2002). "Trayectorias bohmianas y distribuciones cuánticas del espacio de fases". Physics Letters A . 302 (5–6): 261–272. arXiv : quant-ph/0208156v1 . Código Bibliográfico :2002PhLA..302..261D. doi :10.1016/s0375-9601(02)01175-1. ISSN 0375-9601. S2CID 39936409.
^ BJ Hiley: Descripciones del espacio de fases de los fenómenos cuánticos , en: A. Khrennikov (ed.): Teoría cuántica: reconsideración de los fundamentos–2 , págs. 267–286, Växjö University Press, Suecia, 2003 (PDF).
^ B. Hiley: Función característica de Moyal, la matriz de densidad y el idempotente de von Neumann (preimpresión).
^ FC Khanna, PA Mello, M. Revzen, Reconstrucción de estados mecánicos clásicos y cuánticos, arXiv:1112.3164v1 [quant-ph] (enviado el 14 de diciembre de 2011).
^ Bazarov, Ivan V. (3 de mayo de 2012). "Representación de la radiación de sincrotrón en el espacio de fases". Physical Review Special Topics - Accelerators and Beams . 15 (5). American Physical Society (APS): 050703. arXiv : 1112.4047 . Bibcode :2012PhRvS..15e0703B. doi : 10.1103/physrevstab.15.050703 . ISSN 1098-4402. S2CID 53489256.
^ W. Heisenberg (1931). "Über die inkohärente Streuung von Röntgenstrahlen". Physikalische Zeitschrift . 32 : 737–740.
^ Dirac, PAM (1930). "Nota sobre fenómenos de intercambio en el átomo de Thomas". Actas matemáticas de la Sociedad filosófica de Cambridge . 26 (3). Cambridge University Press (CUP): 376–385. Bibcode :1930PCPS...26..376D. doi : 10.1017/s0305004100016108 . ISSN 0305-0041. S2CID 97185766.
^ Ann Moyal, (2006), "Matemático inconformista: La vida y la ciencia de JE Moyal", ANU E-press, 2006, ISBN 1-920942-59-9 .
Lectura adicional
M. Levanda y V. Fleurov, "Función de cuasi-distribución de Wigner para partículas cargadas en campos electromagnéticos clásicos", Annals of Physics , 292 , 199–231 (2001). arXiv :cond-mat/0105137.
Enlaces externos
Implementación de la función Wigner en QuTiP.
Galería de Óptica Cuántica.
Software gratuito con licencia GPL para la distribución de archivos de señales de cuasiprobabilidad de Wigner.