Informe de grado modalidad semillero investigación: métricas de distinguibilidad de estados cuánticos e introducción a los algoritmos cuánticos variacionales
En este informe se presenta el trabajo desarrollado como estudiante miembro del semillero de investigación en Información Cuántica y Computación Cuántica, y que fue presentado en eventos científicos como ponencias orales y posters. El trabajo está dividido en dos partes principales. En primer lugar, se llevó a cabo un análisis de distinguibilidad de estados cuánticos por medio de la comparación de la distancia traza y la fidelidad, como métricas en el campo de la información cuántica. Se discutió la dinámica de estados cuánticos de dos qubits mediante la evolución temporal de dichas métricas al considerar la acción de ruido de amortiguación de amplitud sobre un estado inicial. En segundo lugar, se estudió la técnica de los algoritmos cuánticos variacionales (VQAs) como un enfoque para mejorar la precisión de las computadoras cuánticas. Los VQAs aprovechan procesadores clásicos para entrenar circuitos cuánticos mediante redes neuronales y luego se miden utilizando un procesador cuántico. Se presentaron varios sistemas cuánticos simulados utilizando estos algoritmos, comparando la efectividad de la técnica al simular sistemas sin optimización y sistemas optimizados. Se examinó el comportamiento de la función de costo de cada circuito cuántico para diferentes dinámicas de los sistemas simulados. Para realizar las optimizaciones de los circuitos, se utilizó la biblioteca Pennylane en Python, la cual proporciona diversos optimizadores variacionales cuánticos.
En este informe se presenta el trabajo desarrollado como estudiante miembro del semillero de investigación en Información Cuántica y Computación Cuántica, y que fue presentado en eventos científicos como ponencias orales y posters. El trabajo está dividido en dos partes principales. En primer lugar, se llevó a cabo un análisis de distinguibilidad de estados cuánticos por medio de la comparación de la distancia traza y la fidelidad, como métricas en el campo de la información cuántica. Se discutió la dinámica de estados cuánticos de dos qubits mediante la evolución temporal de dichas métricas al considerar la acción de ruido de amortiguación de amplitud sobre un estado inicial. En segundo lugar, se estudió la técnica de los algoritmos cuánticos variacionales (VQAs) como un enfoque para mejorar la precisión de las computadoras cuánticas. Los VQAs aprovechan procesadores clásicos para entrenar circuitos cuánticos mediante redes neuronales y luego se miden utilizando un procesador cuántico. Se presentaron varios sistemas cuánticos simulados utilizando estos algoritmos, comparando la efectividad de la técnica al simular sistemas sin optimización y sistemas optimizados. Se examinó el comportamiento de la función de costo de cada circuito cuántico para diferentes dinámicas de los sistemas simulados. Para realizar las optimizaciones de los circuitos, se utilizó la biblioteca Pennylane en Python, la cual proporciona diversos optimizadores variacionales cuánticos.
En un principio, la Teoría de la Información fue concebida en forma abstracta, en base a una formulación matemática. Sin embargo, la habilidad y eficiencia para procesar o transmitir información dependerá en última instancia del sistema físico elegido para codificar esa información. Es así que, al codificar la información en las propiedades físicas de sistemas cuyo comportamiento se describe por las leyes de la mecánica cuántica, se abren nuevas posibilidades que son irrealizables clásicamente. Entre ellas, la posibilidad de generar algoritmos cuánticos para ciertos cálculos que requieren un número de pasos significativamente menor que cualquier algoritmo clásico, reduciendo su complejidad. Entre las características fundamentales de la mecánica cuántica que posibilitan estos desarrollos, sobresalen el principio de superposición y el entrelazamiento cuántico. Una de las implementaciones más accesibles para la comunicación cuántica y el procesamiento cuántico de la información emplea fotones individuales como portadores de la información. Estos pueden viajar entre dos estaciones distantes siendo poco afectados por el ruido que introduce el entorno. Además, poseen varios grados de libertad en los cuales es posible codificar la información y pueden manipularse con un alto grado de control utilizando tecnología estándar. En este trabajo se describen algunas aplicaciones de la implementación fotónica, realizadas en colaboración con el Laboratorio de Procesado de Imágenes (LPI) del Depto. de Física de la UBA. En primer lugar, se describe la determinación experimental de medidas de correlaciones cuánticas en estados no puros de un par de fotones correlacionados en polarización, mediante técnicas de tomografía de estados cuánticos. Luego se describe la utilización de moduladores espaciales de luz programables para entrelazar la polarización con grados de libertad espaciales, y así poder simular los denominados estados historia de evoluciones temporales cuánticas. Se muestra que este último proveeun esquema eficientepara el cálculo de promedios temporales de observables del sistema.
Se pretende que el alumno desarrolle un proyecto de investigación asociado a la simulación numérica de la dinámica y control de sistemas estocásticos con potenciales aplicaciones en el novedoso campo de las tecnologías cuánticas. Objetivo 1) Utilizar software adaptado al análisis y simulación numérica de ecuaciones diferenciales estocásticas. Objetivo 2) Modelización y simulación de sistemas lineales estocásticos análogos a sistemas cuánticos lineales. Objetivo 3) Planteamiento de protocolos de control óptimo sobre estos sistemas. Objetivo 4) Simulación de la dinámica y control óptimo de un sistema cuántico lineal. \nEl control es una rama interdisciplinar de la ingeniería y las matemáticas, aplicada al comportamiento de los sistemas dinámicos. Un controlador manipula la entrada o estado inicial del sistema para lograr un efecto deseado en un estado posterior o salida. Hoy en día, prácticamente todo sistema mecánico y electrónico posee un sistema de control asociado, por sencillo que sea. Un sistema de control simple consiste en una acción sobre el sistema y la obtención del resultado esperado sobre el mismo. Sin embargo, para evitar un final abierto, los sistemas realizan, por lo regular, mediciones intermedias sobre sus salidas y con esa información el sistema es retroalimentado para modificar la acción de control reiniciando el ciclo de entrada-salida. Por otra parte, el desarrollo de nuevas tecnologías basadas en dispositivos cuyas leyes vienen dictadas por la mecánica cuántica se ha convertido en una prioridad para el desarrollo tecnológico actual y diversas aplicaciones tales como las comunicaciones, sensores y computadores cuánticos esperan a su desarrollo para ser implementadas. Esto exige que la rama del control cuántico se desarrolle. Sin embargo, algunas cualidades de la mecánica cuántica no permiten pasar fácilmente el esquema del control clásico al ámbito cuántico. Clásicamente, es posible adquirir toda la información necesaria y emplearla de manera eficiente para retroalimentar el sistema y especializar la acción de control al estado actual. Por el contrario, en un sistema cuántico no toda la información es asequible mediante una medición y, además, todo intento de medir el sistema lo perturba de manera notable. Técnicamente, las ecuaciones dinámicas en un proceso tal, son ecuaciones diferenciales estocásticas no lineales en lugar de las ecuaciones diferenciales deterministas de los sistemas clásicos. En este proyecto, el alumno implementará un software de simulación de sistemas estocásticos adaptado a las peculiaridades de sistemas cuánticos lineales para a continuación, ensayar la validez de diversos protocolos de control óptimo propuestos para estos sistemas
Desde su origen, la teoria de la informacion ha tenido fuertes conexiones con la mecanica estadistica. El mismo termino entropia de la informacion fue elegido por Shannon a partir del termino usado en termodinamica, bajo sugerencia de Von Neumann. La misma relacion se ha desarrollado recientemente entre las correspondientes generalizaciones cuanticas de ambas teorias, que son la informacion cuantica, que estudia el almacenamiento y la manipulacion de la informacion en sistemas cuanticos, y la fisica de muchos cuerpos. Las conexiones entre los dos campos aumentan cada dia y van en las dos direcciones. Las herramientas e ideas de la informacion cuantica ayudan a solucionar problemas abiertos en teoria de la materia condensada, y nuevos modelos de muchos cuerpos se desarrollan para aplicaciones de la informacion cuantica. Al mismo tiempo la implementacion y el control experimental de pequenos sistemas cuanticos ha mejorado de forma espectacular, aumentando la posibilidad de que estos experimentos se puedan llevar a cabo a escala mas grande, abriendo la posibilidad de aplicaciones practicas.
Es reciente el interes de la comunidad de investigadores en informacion cuantica por estos modelos, y se ha empezado a ver este tipo de evoluciones no solamente como un problema que resolver, sino como un recurso que explotar. Este fundamental cambio de perspectiva ha abierto una nueva linea de investigacion en el campo, incorporando ideas sobre como construir evoluciones artificiales en las cuales la disipacion trabaja en nuestro favor en vez de en nuestra contra, incluso protegiendo el sistema de otros ruidos incontrolados. De esta manera se han planteado nuevos problemas interesantes, como cual es el poder computacional de estos modelos y que condiciones pueden garantizar la resistencia contra el ruido externo.
Esta tesis es parte del esfuerzo para mejorar nuestro conocimiento de estos modelos.
Nos hemos centrado en el estudio de las propiedades de evoluciones disipativas de sistemas cuanticos de espines en un reticulo, bajo una hipotesis que llamamos equilibracion rapida. Consideramos evoluciones disipativas con un unico punto fijo, y que comprimen todo el espacios de estados iniciales en un entorno decreciente del punto fijo. La escala temporal a la cual se tiene esta compresion, es decir el tiempo que tenemos que esperar para que cualquier estado inicial se encuentre arbitrariamente cercano al punto fijo, se llama tiempo de equilibracion del proceso. La condicion de equilibracion rapida pide que este tiempo sea logaritmico en el tamano del sistema.
La contribucion principal de esta tesis es demostrar que la equilibracion rapida tiene importantes consecuencias en las propiedades del sistema. Este sera estable bajo perturbaciones externas y sus puntos fijos satisfaran una ley de area para la informacion mutua.
La parte mas tecnica de este trabajo se basa en las llamadas cotas de Lieb Robinson. Muy conocidas en el caso de sistemas cerrados, pero igualmente validas para sistemas abiertos, las cotas de Lieb Robinson formalizan la idea que la informacion en un sistema de muchos cuerpos solo puede moverse con una velocidad finita, dado que su propagacion esta mediada por las interacciones locales. Por lo tanto, controlar la escala temporal con la cual la evolucion converge permite a su vez controlar la escala espacial a la cual se pueden crear correlaciones.
La computación cuántica ha impulsado múltiples líneas de investigación en los últimos años, debido al interés conjunto de la comunidad científica y de la industria tecnológica en diseñar ordenadores cuánticos capaces de ejecutar circuitos con alta precisión. Como resultado, nuevos enfoques como el MBQC (computación cuántica basada en mediciones) han surgido como alternativas viables al modelo estándar de computación cuántica. De igual manera, empresas tecnológicas como Google o IBM continúan impulsando de manera sostenida el desarrollo del hardware para computación cuántica, logrando avances significativos como el procesador Heron de 156 cúbits de IBM o el chip Willow de Google para la corrección de errores cuánticos. En particular, la introducción de Qiskit como software libre, por parte de IBM, brinda a los usuarios la posibilidad de realizar experimentos de computación cuántica utilizando ordenadores cuánticos reales, lo que sugiere la importancia de Qiskit dentro de este campo. Por tal motivo, el objetivo de este trabajo es diseñar una plantilla en Python con Qiskit para aplicar el protocolo MQBC en diferentes circuitos y ejecutarlos en un ordenador cuántico de IBM. Esta plantilla se basa en el modelo One-Way, un enfoque particular del MQBC, aplicado a un clúster lineal de 4 cúbits, y utiliza técnicas como la Postselección, Reetiquetado y Tomografía cuántica para caracterizar los estados de salida generados. Con ello, se busca implementar 5 circuitos cuánticos correspondientes a diferentes configuraciones clústers lineales y bidimensionales, a partir del estado clúster lineal de 4 cúbits. Asimismo, se mostrarán los pasos para implementar un circuito que detalla una búsqueda cuántica de dos cúbits. Para verificar la validez de los resultados, se calculará la fidelidad para cada estado de salida, considerando los circuitos teóricos descritos en el paper Experimental One-Way Quantum Computing. Finalmente, para las configuraciones de clúster bidimensionales, se analizará el grado de entrelazamiento para confirmar la generación de estados entrelazados a partir de cúbits codificados inicialmente separables. En adición, esta tesis tiene como finalidad servir como documento de referencia para facilitar el aprendizaje de Qiskit en futuros trabajos de computación cuántica.
Este trabajo proporciona una visión generalizada del estado actual de los algoritmos genéticos cuánticos (QGAs), mostrando los avances realizados en esta área de investigación los últimos 24 años. Los QGAs combinan conceptos de la computación cuántica y los algoritmos genéticos clásicos (CGAs), lo que les permite abordar problemas complejos de optimización y búsqueda de manera eficiente. Se presentan los principales hallazgos y contribuciones de estos algoritmos cuánticos destacando las tendencias y los enfoques más prometedores, así como los desafíos y limitaciones que deben superarse. Se presentan nuevos enfoques y técnicas de implementación de QGAs, incluyendo operadores genéticos cuánticos y esquemas de codificación eficientes que contribuyen a mejorar el rendimiento y la convergencia de los algoritmos. Se comparan los QGAs y otros enfoques similares, como los CGAs y los algoritmos cuánticos puros, destacando las ventajas y desventajas relativas de los QGAs en comparación a sus versiones clásicas. Se muestra también una implementación de QGA utilizando la biblioteca Qiskit. Se presentan la selección de los métodos usados para la generación de la población inicial, el cruzamiento y la mutación de las distintas poblaciones de los circuitos cuánticos simulados en los experimentos realizados, ejemplificando las ventajas significativas que estos pueden traer en comparación con los enfoques clásicos.
La motivacion de esta tesis es estudiar el entrelazamiento cuantico en general y bajo que circunstancias es maximo en particular. En primer lugar, se analizan desigualdades de Bell desde un punto de vista operacional, se presentan nuevas desigualdades para particulas de dimension mayor que dos y se estudian los estados que las violan maximamente. Dichos estados son del tipo GHZ para las desigualdades que involucran qubits. En el caso de experimentos con qutrits o particulas dimension local mayor, los estados que violan maximamente estas desigualdades son deformaciones de estados GHZ generalizados. Este resultado demuestra que aunque maximo entrelazamiento y la no-localidad son conceptos muy proximos, no son equivalentes. En segundo lugar, se estudia el entrelazamiento multipartito en cadenas de espines. En particular, se usa como figura de merito el hiperdeterminante y los invariantes polinomicos que lo forman, S y T. Estas figuras cuantifican un tipo concreto de entrelazamiento cuatripartito, tal y como se demuestra de un analisis de estados cuanticos bien conocidos como el GHZ o el W. Cuando aplicamos esta figura a una cadena de espines con interaccion tipo Ising, observamos un pico pronunciado cerca de la transicion de fase cuantica. Resultados similares se observan para los modelos de XXZ y de Haldane-Shastry para el caso de los invariantes S y T. Con todo, concluimos que estas figuras de merito son sensibles a las transiciones de fase cuanticas. La segunda parte de la tesis se focaliza en el campo de la computacion cuantica. Dicho campo ha experimentado un crecimiento relevante en los ultimos anos, debido en parte a la mejora tecnologica sobre el control de experimentos cuanticos. Con motivo de este crecimiento, diversas empresas han empezado a desarrollar prototipos de ordenadores cuanticos. Por esta razon, una de las lineas de investigacion en computacion cuantica consiste en proponer metodos para poner a prueba estas maquinas. Uno de los metodos presentados en esta tesis es la simulacion exacta de una cadena de espines con interaccion de tipo XY. Se explicita como debe de ser el circuito cuantico que lleva a cabo dicha simulacion y se propone la simulacion tanto de estados termicos como de la evolucion en el tiempo. Debido a que el modelo XY es resoluble analiticamente, se pueden comparar los resultados obtenidos de la simulacion con el ordenador cuantico con la solucion teorica. Al correr un circuito para cuatro qubits en dos ordenadores de la empresa IBM y uno de la empresa Rigetti computing, obtenemos resultados peores a los esperados. Con ello, concluimos que hay fuentes de error que no se estan teniendo en cuenta en general y que cobran relevancia incluso en circuitos tan pequenos. Ademas de este circuito, proponemos la simulacion de estados altamente entrelazados. La implementacion de estos circuitos constituye un test tan duro como necesario para un ordenador cuantico, debido a que la ventaja de los algoritmos cuanticos radica en la generacion de entrelazamiento. Para finalizar, estudiamos la generacion de entrelazamiento al nivel mas fundamental: la fisica de particulas. Obtenemos que la interaccion de QED viene dictada por el maximo entrelazamiento en aquellos procesos a “tree-level”. Aplicamos la misma filosofia en procesos que involucren corrientes debiles neutras (tambien a “tree-level”) y con ello obtenemos que el angulo de Weinberg debe ser de pi/6 si imponemos que las particulas finales deben estar maximamente entrelazadas, valor muy cercano al experimental.